Solving Dirac equations on a 3D lattice with inverse Hamiltonian and spectral methods

A new method to solve the Dirac equation on a 3D lattice is proposed, in which the variational collapse problem is avoided by the inverse Hamiltonian method and the fermion doubling problem is avoided by performing spatial derivatives in momentum space with the help of the discrete Fourier transform, i.e., the spectral method. This method is demonstrated in solving the Dirac equation for a given spherical potential in 3D lattice space. In comparison with the results obtained by the shooting method, the differences in single particle energy are smaller than $10^{-4}$~MeV, and the densities are almost identical, which demonstrates the high accuracy of the present method. The results obtained by applying this method without any modification to solve the Dirac equations for an axial deformed, non-axial deformed, and octupole deformed potential are provided and discussed.Comment: 18 pages, 6 figure

Possibility of Unconventional Pairing Due to Coulomb Interaction in Fe-Based Pnictide Superconductors: Perturbative Analysis of Multi-Band Hubbard Models

Possibility of unconventional pairing due to Coulomb interaction in iron-pnictide superconductors is studied by applying a perturbative approach to realistic 2- and 5-band Hubbard models. The linearized Eliashberg equation is solved by expanding the effective pairing interaction perturbatively up to third order in the on-site Coulomb integrals. The numerical results for the 5-band model suggest that the eigenvalues of the Eliashberg equation are sufficiently large to explain the actual high Tc for realistic values of Coulomb interaction and the most probable pairing state is spin-singlet s-wave without any nodes just on the Fermi surfaces, although the superconducting order parameter changes its sign between the small Fermi pockets. On the other hand the 2-band model is quite insufficient to explain the actual high Tc.Comment: 2 pages, 3 figures. Proceedings of the Intl. Symposium on Fe-Oxypnictide Superconductors (Tokyo, 28-29th June 2008

Promotion of cooperation induced by nonlinear attractive effect in spatial Prisoner's Dilemma game

We introduce nonlinear attractive effects into a spatial Prisoner's Dilemma game where the players located on a square lattice can either cooperate with their nearest neighbors or defect. In every generation, each player updates its strategy by firstly choosing one of the neighbors with a probability proportional to $\mathcal{A}^\alpha$ denoting the attractiveness of the neighbor, where $\mathcal{A}$ is the payoff collected by it and $\alpha$ ($\geq$0) is a free parameter characterizing the extent of the nonlinear effect; and then adopting its strategy with a probability dependent on their payoff difference. Using Monte Carlo simulations, we investigate the density $\rho_C$ of cooperators in the stationary state for different values of $\alpha$. It is shown that the introduction of such attractive effect remarkably promotes the emergence and persistence of cooperation over a wide range of the temptation to defect. In particular, for large values of $\alpha$, i.e., strong nonlinear attractive effects, the system exhibits two absorbing states (all cooperators or all defectors) separated by an active state (coexistence of cooperators and defectors) when varying the temptation to defect. In the critical region where $\rho_C$ goes to zero, the extinction behavior is power law-like $\rho_C$ $\sim$ $(b_c-b)^{\beta}$, where the exponent $\beta$ accords approximatively with the critical exponent ($\beta\approx0.584$) of the two-dimensional directed percolation and depends weakly on the value of $\alpha$.Comment: 7 pages, 4 figure