Self-consistent calculation of particle-hole diagrams on the Matsubara frequency: FLEX approximation

We implement the numerical method of summing Green function diagrams on the Matsubara frequency axis for the fluctuation exchange (FLEX) approximation. Our method has previously been applied to the attractive Hubbard model for low density. Here we apply our numerical algorithm to the Hubbard model close to half filling ($\rho = 0.40$), and for $T/t = 0.03$, in order to study the dynamics of one- and two-particle Green functions. For the values of the chosen parameters we see the formation of three branches which we associate with the a two-peak structure in the imaginary part of the self-energy. From the imaginary part of the self-energy we conclude that our system is a Fermi liquid (for the temperature investigated here), since Im$\Sigma(\vec{k},\omega) \approx w^2$ around the chemical potential. We have compared our fully self-consistent FLEX solutions with a lower order approximation where the internal Green functions are approximated by free Green functions. These two approches, i.e., the fully selfconsistent and the non-selfconsistent ones give different results for the parameters considered here. However, they have similar global results for small densities.Comment: seven pages, nine figures as ps files. Accepted in Int. J. Modern Phys. C (1997

Temperature and doping dependence of normal state spectral properties in a two-orbital model for ferropnictides

Using a second-order perturbative Green's functions approach we determined the normal state single-particle spectral function $A(\vec{k},\omega)$ employing a minimal effective model for iron-based superconductors. The microscopic model, used before to study magnetic fluctuations and superconducting properties, includes the two effective tight-binding bands proposed by S.Raghu et al. [Phys. Rev. B 77, 220503 (R) (2008)], and intra- and inter-orbital local electronic correlations, related to the Fe-3d orbitals. Here, we focus on the study of normal state electronic properties, in particular the temperature and doping dependence of the total density of states, $A(\omega)$, and of $A(\vec{k},\omega)$ in different Brillouin zone regions, and compare them to the existing angle resolved photoemission spectroscopy (ARPES) and previous theoretical results in ferropnictides. We obtain an asymmetric effect of electron and hole doping, quantitative agreement with the experimental chemical potential shifts as a function of doping, as well as spectral weight redistributions near the Fermi level as a function of temperature consistent with the available experimental data. In addition, we predict a non-trivial dependence of the total density of states with the temperature, exhibiting clear renormalization effects by correlations. Interestingly, investigating the origin of this predicted behaviour by analyzing the evolution with temperature of the k-dependent self-energy obtained in our approach, we could identify a number of specific Brillouin zone points, none of them probed by ARPES experiments yet, where the largest non-trivial effects of temperature on the renormalization are present.Comment: Manuscript accepted in Physics Letters A on Feb. 25, 201

Normal state electronic properties of LaO1−x_{1-x}Fx_{x}BiS2_{2} superconductors

A good description of the electronic structure of BiS$_{2}$-based superconductors is essential to understand their phase diagram, normal state and superconducting properties. To describe the first reports of normal state electronic structure features from angle resolved photoemission spectroscopy (ARPES) in LaO$_{1-x}$F$_{x}$BiS$_{2}$, we used a minimal microscopic model to study their low energy properties. It includes the two effective tight-binding bands proposed by Usui et al [Phys.Rev.B 86, 220501(R)(2012)], and we added moderate intra- and inter-orbital electron correlations related to Bi-($p_{Y}$, $p_{X}$) and S-($p_{Y}$, $p_{X}$) orbitals. We calculated the electron Green's functions using their equations of motion, which we decoupled in second-order of perturbations on the correlations. We determined the normal state spectral density function and total density of states for LaO$_{1-x}$F$_{x}$BiS$_{2}$, focusing on the description of the k-dependence, effect of doping, and the prediction of the temperature dependence of spectral properties. Including moderate electron correlations, improves the description of the few experimental ARPES and soft X-ray photoemission data available for LaO$_{1-x}$F$_{x}$BiS$_{2}$. Our analytical approximation enabled us to calculate the spectral density around the conduction band minimum at $\vec{k}_{0}=(0.45\pi,0.45\pi)$, and to predict the temperature dependence of the spectral properties at different BZ points, which might be verified by temperature dependent ARPES.Comment: 9 figures. Manuscript accepted in Physica B: Condensed Matter on Jan. 25, 201

Metal-Insulator Transition in the Two-Dimensional Hubbard Model at Half-Filling with Lifetime Effects within the Moment Approach

We explore the effect of the imaginary part of the self-energy, $Im\Sigma(\vec{k},\omega)$, having a single pole, $\Omega(\vec{k},\omega)$, with spectral weight, $\alpha(\vec{k})$, and quasi-particle lifetime, $\Gamma(\vec{k})$, on the density of states. We solve the set of parameters, $\Omega(\vec{k},\omega$), $\alpha(\vec{k})$, and $\Gamma(\vec{k})$ by means of the moment approach (exact sum rules) of Nolting. Our choice for $\Sigma(k,\omega)$, satisfies the Kramers - Kronig relationship automatically. Due to our choice of the self - energy, the system is not a Fermi liquid for any value of the interaction, a result which is also true in the moment approach of Nolting without lifetime effects. By increasing the value of the local interaction, $U/W$, at half-filling ($\rho = 1/2$), we go from a paramagnetic metal to a paramagnetic insulator, (Mott metal - insulator transition ($MMIT$)) for values of $U/W$ of the order of $U/W \geq 1$ ($W$ is the band width) which is in agreement with numerical results for finite lattices and for infinity dimensions ($D = \infty$). These results settle down the main weakness of the spherical approximation of Nolting: a finite gap for any finite value of the interaction, i.e., an insulator for any finite value of $U/W$. Lifetime effects are absolutely indispensable. Our scheme works better than the one of improving the narrowing band factor, $B(\vec{k})$, beyond the spherical approximation of Nolting.Comment: 5 pages and 5 ps figures (included

Effective t-J Hamiltonian for the Copper Oxides

Starting from the Emery model, which is assumed to describe the copper oxygen planes, and including direct oxygen hopping matrix elements, we have been able to derive the effective t-J Hamiltonian for the copper orbitals using the Linked Cluster Expansion Method up to fourth order in the hybridization matrix element.Comment: (ps version of the dvi file, resubmitted because previous uucompressed version was corrupted), 9 page