Gutzwiller density functional studies of FeAs-based superconductors: structure optimization and evidence for a three-dimensional Fermi surface

The Gutzwiller conjugate gradient minimization method for correlated electron systems

We review our recent work on the Gutzwiller conjugate gradient minimization method, anab initioapproach developed for correlated electron systems. The complete formalism has been outlined that allows for a systematic understanding of the method, followed by a discussion of benchmark studies of dimers, one- and two-dimensional single-band Hubbard models. In the end, we present some preliminary results of multi-band Hubbard models and large-basis calculations of F₂to illustrate our efforts to further reduce the computational complexity.

Correlation of Tunable CoSi4 Tetrahedron with the Superconducting Properties of LaCoSi

It is known that as the FeAs4 tetrahedron in the Fe-based superconductor is close to the regular tetrahedron, critical temperature (Tc) can be greatly increased. Recently, a Co-based superconductor of LaCoSi (4 K) with "111" structure was found. In this work, we improve the Tc of LaCoSi through structural regulation. Tc can be increased by the chemical substitution of Co by Fe, while the superconductivity is suppressed by the Ni substitution. The combined analysis of neutron and synchrotron X-ray powder diffractions reveals that the change of the Si-Co-Si bond angles of the CoSi4 tetrahedron is possibly responsible for the determination of superconducting properties. The Fe chemical substitution is favorable for the formation of the regular tetrahedron of CoSi4. The present new Co-based superconductor of LaCoSi provides a possible method to enhance the superconductivity performance of the Co-based superconductors via controlling Co-based tetrahedra similar to those well established in the Fe-based superconductors.

First-principles calculation of correlated electron materials based on Gutzwiller wave function beyond Gutzwiller approximation

We propose an approach that is under the framework of Gutzwiller wave function but goes beyond the commonly adopted Gutzwiller approximation to improve the accuracy and flexibility in treating the correlation effects. Detailed formalism is described for a dimer which is straightforwardly generalized later to more complicated periodic bulk systems. The accuracy of the approach is demonstrated by evaluating the potential energy curves of spin-singlet N₂ dimer, spin-triplet O₂ dimer, and 1D hydrogen chain. The computational workload of the approach can be easily handled by efficient parallel computing.

Benchmark of correlation matrix renormalization method in molecule calculations

We report benchmark calculations of the correlation matrix renormalization (CMR) approach for 23 molecules in the well-established G2 molecule set. This subset represents molecules with spin-singlet ground state in a variety of chemical bonding and coordination environments. The QUAsi-atomic minimal basis-set orbitals (QUAMBOs) are used as local orbitals in both CMR and full configuration interaction (FCI) calculations for comparison. The results obtained from the calculations are also compared with available experimental data. It is shown that the CMR method produces binding and dissociation energy curves in good agreement with the QUAMBO-FCI calculations as well as experimental results. The CMR benchmark calculations yield a standard deviation of 0.09 Å for the equilibrium bond length and 0.018 Hartree/atom for the formation energy, with a gain of great computational efficiency which scales like Hartree-Fock method.

Quantum fluctuations beyond the Gutzwiller approximation

We present a simple scheme to evaluate linear response functions including quantum fluctuation corrections on top of the Gutzwiller approximation. The method is derived for a generic multiband lattice Hamiltonian without any assumption about the dynamics of the variational correlation parameters that define the Gutzwiller wavefunction, and which thus behave as genuine dynamical degrees of freedom that add on those of the variational uncorrelated Slater determinant. We apply the method to the standard half-filled single-band Hubbard model. We are able to recover known results, but, as by-product, we also obtain few novel ones. In particular, we show that quantum fluctuations can reproduce almost quantitatively the behaviour of the uniform magnetic susceptibility uncovered by dynamical mean field theory, which, though enhanced by correlations, is found to be smooth across the paramagnetic Mott transition. By contrast, the simple Gutzwiller approximation predicts that susceptibility to diverge at the transition.

Efficient and accurate treatment of electron correlations with Correlation Matrix Renormalization theory

We present an efficient method for calculating the electronic structure and total energy of strongly correlated electron systems. The method extends the traditional Gutzwiller approximation for one-particle operators to the evaluation of the expectation values of two particle operators in the many-electron Hamiltonian. The method is free of adjustable Coulomb parameters, and has no double counting issues in the calculation of total energy, and has the correct atomic limit. We demonstrate that the method describes well the bonding and dissociation behaviors of the hydrogen and nitrogen clusters, as well as the ammonia composed of hydrogen and nitrogen atoms. We also show that the method can satisfactorily tackle great challenging problems faced by the density functional theory recently discussed in the literature. The computational workload of our method is similar to the Hartree-Fock approach while the results are comparable to high-level quantum chemistry calculations.

Calculation of the specific heat of optimally K-doped BaFe₂As₂

The calculated specific heat of optimally K-doped BaFe2As2 in density functional theory is about five times smaller than that found in the experiment. We report that by adjusting the potential on the iron atom to be slightly more repulsive for electrons improves the calculated heat capacity as well as the electronic band structure of Ba0.6K0.4Fe2As2. In addition, structural and magnetic properties are moved in the direction of experimental values. Applying the same correction to the antiferromagnetic state, we find that the electron-phonon coupling is strongly enhanced.

Dynamical correlations and screened exchange on the experimental bench: spectral properties of the cobalt pnictide BaCo2As2

Understanding the Fermi surface and low-energy excitations of iron or cobalt pnictides is crucial for assessing electronic instabilities such as magnetic or superconducting states. Here, we propose and implement a new approach to compute the low-energy properties of correlated electron materials, taking into account both screened exchange beyond the local density approximation and local dynamical correlations. The scheme allows us to resolve the puzzle of BaCo2As2, for which standard electronic structure techniques predict a ferromagnetic instability not observed in nature.

Correlation-induced self-doping in the iron-pnictide superconductor Ba2Ti2Fe2As4O

The electronic structure of the iron-based superconductor Ba2Ti2Fe2As4O (Tc(onset)=23.5  K) has been investigated by using angle-resolved photoemission spectroscopy and combined local density approximation and dynamical mean field theory calculations. The electronic states near the Fermi level are dominated by both the Fe 3d and Ti 3d orbitals, indicating that the spacer layers separating different FeAs layers are also metallic. By counting the enclosed volumes of the Fermi surface sheets, we observe a large self-doping effect; i.e., 0.25 electrons per unit cell are transferred from the FeAs layer to the Ti2As2O layer, leaving the FeAs layer in a hole-doped state. This exotic behavior is successfully reproduced by our dynamical mean field calculations, in which the self-doping effect is attributed to the electronic correlations in the 3d shells. Our work provides an alternative route of effective doping without element substitution for iron-based superconductors.

Many-body effects in iron pnictides and chalcogenides: nonlocal versus dynamic origin of effective masses

We apply the quasiparticle self-consistent GW approximation (QSGW) to some of the iron pnictide and chalcogenide superconductors. We compute Fermi surfaces and density of states, and find excellent agreement with experiment, substantially improving over standard band-structure methods. Analyzing the QSGW self-energy we discuss nonlocal and dynamic contributions to effective masses. We present evidence that the two contributions are mostly separable, since the quasiparticle weight is found to be essentially independent of momentum. The main effect of nonlocality is captured by the static but nonlocal QSGW effective potential. Moreover, these nonlocal self-energy corrections, absent in, e.g., dynamical mean field theory, can be relatively large. We show, on the other hand, that QSGW only partially accounts for dynamic renormalizations at low energies. These findings suggest that QSGW combined with dynamical mean field theory will capture most of the many-body physics in the iron pnictides and chalcogenides.

High energy pseudogap and its evolution with doping in Fe-based superconductors as revealed by optical spectroscopy

We report optical spectroscopic measurements on electron- and hole-doped BaFe2As2. We show that the compounds in the normal state are not simple metals. The optical conductivity spectra contain, in addition to the free carrier response at low frequency, a temperature-dependent gap-like suppression at fairly high energy scale near 0.6 eV. This suppression evolves with the As–Fe–As bond angle induced by electron or hole doping. Furthermore, the feature becomes much weaker in the Fe-chalcogenide compounds. We elaborate that the feature is mainly caused by the strong Hund's rule coupling effect between the itinerant electrons and localized electron moment arising from the multiple Fe 3d orbitals. The coupling strength changes with the environment of the Fe atom. Our experiments demonstrate the coexistence of itinerant and localized electrons in iron-based compounds, which would then lead to a more comprehensive picture of the metallic magnetism in the materials.

Time-dependent Gutzwiller theory for multiband Hubbard models

Based on the variational Gutzwiller theory, we present a method for the computation of response functions for multiband Hubbard models with general local Coulomb interactions. The improvement over the conventional random-phase approximation is exemplified for an infinite-dimensional two-band Hubbard model where the incorporation of the local multiplet structure leads to a much larger sensitivity of ferromagnetism on the Hund coupling. Our method can be implemented into local-density approximation and Gutzwiller schemes and will therefore be an important tool for the computation of response functions for strongly correlated materials.

Incoherent c-axis interplane response of the iron chalcogenide FeTe(0.55)Se(0.45) superconductor from infrared spectroscopy

We report on the interplane c-axis electronic response of FeTe(0.55)Se(0.45) investigated by infrared spectroscopy. We find that the normal-state c-axis electronic response of FeTe(0.55)Se(0.45) is incoherent and bears significant similarities to those of mildly underdoped cuprates. The c-axis optical conductivity σ(c)(ω) of FeTe(0.55)Se(0.45) does not display well-defined Drude response at all temperatures. As temperature decreases, σ(c)(ω) is continuously suppressed. The incoherent c-axis response is found to be related to the strong dissipation in the ab-plane transport: a pattern that holds true for various correlated materials as well as FeTe(0.55)Se(0.45).

Antiferromagnetic order in multiband Hubbard models for iron pnictides

We investigate multiband Hubbard models for the three iron 3d t(2g) bands and the two iron 3d e(g) bands in LaOFeAs by means of the Gutzwiller variational theory. Our analysis of the paramagnetic ground state shows that neither Hartree-Fock mean-field theories nor effective spin models describe these systems adequately. In contrast to Hartree-Fock-type approaches, the Gutzwiller theory predicts that antiferromagnetic order requires substantial values of the local Hund's-rule exchange interaction. For the three-band model, the antiferromagnetic moment fits experimental data for a broad range of interaction parameters. However, for the more appropriate five-band model, the iron e(g) electrons polarize the t(2g) electrons and they substantially contribute to the ordered moment.

Measurement of the c-axis optical reflectance of AFe2As2 (A=Ba, Sr) single crystals: evidence of different mechanisms for the formation of two energy gaps

We present the c-axis optical reflectance measurement on single crystals of BaFe2As2 and SrFe2As2, the parent compounds of FeAs based superconductors. Different from the ab-plane optical response where two distinct energy gaps were observed in the spin-density-wave (SDW) state, only the smaller energy gap could be seen clearly for E∥c axis. The very pronounced energy gap structure seen at a higher energy scale for E∥ab plane is almost invisible. We propose a novel picture for the band structure evolution across the SDW transition and suggest different driving mechanisms for the formation of the two energy gaps.

Electron correlation and spin density wave order in iron pnictides

We study the correlation effects on the electronic structure, orbital and spin density wave (SDW) order in Fe pnictides. Using the multiorbital Hubbard model and Gutzwiller projection, we show that correlation effects are essential to stabilize the metallic SDW phase for the intermediate correlation strengths appropriate for pnictides. We find that the ordered moments depend sensitively on Hund's rule coupling J but weakly on the intraorbital Coulomb repulsion U, varying from 0.3μB to 1.5μB in the range J=0.3-0.8  eV for U=3-4  eV. We study the phase diagram, the evolution of the Fermi surface with the ordered moment, the effects of electron doping, and compare to recent experiments.