• Physical properties of multi-orbital materials depend not only on the strength of the effective interactions among the valence electrons but also on their type. Strong correlations are caused by either Mott physics that captures the Coulomb repulsion among charges, or Hund physics that aligns the spins in different orbitals. We identify four energy scales marking the onset and the completion of screening in orbital and spin channels. The differences in these scales, which are manifest in the temperature dependence of the local spectrum and of the charge, spin and orbital susceptibilities, provide clear signatures distinguishing Mott and Hund physics. We illustrate these concepts with realistic studies of two archetypal strongly correlated materials, and corroborate the generality of our conclusions with a model Hamiltonian study.
  • We study the fundamental question of the lattice dynamics of a metallic ferromagnet in the regime where the static long range magnetic order is replaced by the fluctuating local moments embedded in a metallic host. We use the \textit{ab initio} Density Functional Theory(DFT)+embedded Dynamical Mean-Field Theory(eDMFT) functional approach to address the dynamic stability of iron polymorphs and the phonon softening with increased temperature. We show that the non-harmonic and inhomogeneous phonon softening measured in iron is a result of the melting of the long range ferromagnetic order, and is unrelated to the first order structural transition from the BCC to the FCC phase, as is usually assumed. We predict that the BCC structure is dynamically stable at all temperatures at normal pressure, and is only thermodynamically unstable between the BCC-$\alpha$ and the BCC-$\delta$ phase of iron.
  • Two-particle Green's functions and the vertex functions play a critical role in theoretical frameworks for describing strongly correlated electron systems. However, numerical calculations at two-particle level often suffer from large computation time and massive memory consumption. We derive a general expansion formula for the two-particle Green's functions in terms of an overcomplete representation based on the recently proposed "intermediate representation" basis. The expansion formula is obtained by decomposing the spectral representation of the two-particle Green's function. We demonstrate that the expansion coefficients decay exponentially, while all high-frequency and long-tail structures in the Matsubara-frequency domain are retained. This representation therefore enables efficient treatment of two-particle quantities and opens a route to the application of modern many-body theories to realistic strongly correlated electron systems.
  • Iridates abound with interesting magnetic behaviours because of their strong spin-orbit coupling. Sr$_3$NiIrO$_6$ brings together the spin-orbital entanglement of the Ir$^{4+}$ ion with a 3$d$ Ni cation and a one-dimensional crystal structure. It has a ferrimagnetic ground state with a 55 T coercive field. We perform a theoretical study of the magnetic interactions in this compound, and elucidate the role of anisotropic symmetric exchange as the source of its strong magnetic anisotropy. Our first-principles calculations reproduce the magnon spectra of this compound and predict a signature in the cross sections that can differentiate the anisotropic exchange from single-ion anisotopy.
  • We report infrared studies of AFe$_{2}$As$_{2}$ (A = Ba, Sr), two representative parent compounds of iron-arsenide superconductors, at magnetic fields (B) up to 17.5 T. Optical transitions between Landau levels (LLs) were observed in the antiferromagnetic states of these two parent compounds. Our observation of a $\sqrt{B}$ dependence of the LL transition energies, the zero-energy intercepts at B = 0 T under the linear extrapolations of the transition energies and the energy ratio ($\sim$ 2.4) between the observed LL transitions, combined with the linear band dispersions in two-dimensional (2D) momentum space obtained by theoretical calculations, demonstrates the existence of massless Dirac fermions in antiferromagnetic BaFe$_{2}$As$_{2}$. More importantly, the observed dominance of the zeroth-LL-related absorption features and the calculated bands with extremely weak dispersions along the momentum direction $k_{z}$ indicate that massless Dirac fermions in BaFe$_{2}$As$_{2}$ are 2D. Furthermore, we find that the total substitution of the barium atoms in BaFe$_{2}$As$_{2}$ by strontium atoms not only maintains 2D massless Dirac fermions in this system, but also enhances their Fermi velocity, which supports that the Dirac points in iron-arsenide parent compounds are topologically protected.
  • The metal-insulator transition (MIT) remains among the most thoroughly studied phenomena in solid state physics, but the complexity of the phenomena, which usually involves cooperation of many degrees of freedom including orbitals, fluctuating local moments, magnetism, and the crystal structure, have resisted predictive ab-initio treatment. Here we develop ab-initio theoretical method for correlated electron materials, based on Dynamical Mean Field Theory, which can predict the change of the crystal structure across the MIT at finite temperature. This allows us to study the coupling between electronic, magnetic and orbital degrees of freedom with the crystal structure across the MIT in rare-earth nickelates. We predict the free energy profile of the competing states, and the theoretical magnetic ground state configuration, which is in agreement with neutron scattering data, but is different from the magnetic models proposed before. The resonant elastic X-ray response at the K-edge, which was argued to be a direct probe of the charge order, is theoretically modelled within the Dynamical Mean Field Theory, including the core-hole interaction. We show that the line-shape of the measured resonant elastic X-ray response can be explained with the "site-selective" Mott scenario without real charge order on Ni sites.
  • We investigate the evolution of magnetic properties as a function of hydrogen doping in iron based superconductor LaFeAsO$_{1-x}$H$_x$ using the dynamical mean-field theory combined with the density-functional theory. We find that two independent consequences of the doping, the increase of the electron occupation and the structural modification, have the opposite effects on the strength of electron correlation and magnetism, resulting in the minimum of the calculated magnetic moment around the intermediate doping level as a function of $x$. Our result provides a natural explanation for the puzzling recent experimental discovery of the two separated antiferromagnetic phases at low and high doping limits. Furthermore, the increase of orbital occupation and correlation strength with the doping results in reduced orbital polarization of $d_{xz/yz}$ orbitals and the enhanced role of $d_{xy}$ orbital in the magnetism at high doping levels, and their possible implications to the superconductivity are discussed in line with the essential role of the magnetism.
  • We implemented combination of DMFT and $GW$ in its fully self-consistent way, one shot $GW$ approximation, and quasiparticle self-consistent scheme, and studied how well these combined methods perform on H$_2$ molecule as compared to more established methods such as LDA+DMFT. We found that most flavors of $GW$+DMFT break down in strongly correlated regime due to causality violation. Among $GW$+DMFT methods, only the self-consistent quasiparticle $GW$+DMFT with static double-counting, and a new method with causal double-counting, correctly recover the atomic limit at large H-atom separation. While some flavors of $GW$+DMFT improve the single-electron spectra as compared to LDA+DMFT, the total energy is best predicted by LDA+DMFT, for which the exact double-counting is known, and is static.
  • We use inelastic neutron scattering to study temperature and energy dependence of spin excitations in optimally P-doped BaFe2(As0.7P0.3)2 superconductor (Tc = 30 K) throughout the Brillouin zone. In the undoped state, spin waves and paramagnetic spin excitations of BaFe2As2 stem from antiferromagnetic (AF) ordering wave vector QAF= (1/-1,0) and peaks near zone boundary at (1/-1,1/-1) around 180 meV. Replacing 30% As by smaller P to induce superconductivity, low-energy spin excitations of BaFe2(As0.7P0.3)2form a resonance in the superconducting state and high-energy spin excitations now peaks around 220 meV near (1/-1,1/-1). These results are consistent with calculations from a combined density functional theory and dynamical mean field theory, and suggest that the decreased average pnictogen height in BaFe2(As0.7P0.3)2 reduces the strength of electron correlations and increases the effective bandwidth of magnetic excitations.
  • The theoretical understanding of emergent phenomena in quantum materials is one of the greatest challenges in condensed matter physics. In contrast to simple materials such as noble metals and semiconductors, macroscopic properties of quantum materials cannot be predicted by the properties of individual electrons. One of the examples of scientific importance is strongly correlated electron system. Neither localized nor itinerant behaviors of electrons in partially filled 3d, 4f, and 5f orbitals give rise to rich physics such as Mott insulators, high-temperature superconductors, and superior thermoelectricity, but hinder quantitative understanding of low-lying excitation spectrum. Here we present a new first-principles approach to strongly correlated solids. It is Q4 based on a combination of the quasiparticle self-consistent GW approximation and the dynamical mean-field theory. The sole input in this method is the projector to the set of correlated orbitals for which all local Feynman graphs are being evaluated. For that purpose, we choose very localized quasiatomic orbitals spanning large energy window, which contains most strongly hybridized bands, as well as upper and lower Hubbard bands. The self-consistency is carried out on the Matsubara axis. This method enables the first-principles study of Mott insulators in both their paramagnetic and antiferromagnetic phases. We illustrate the method on the archetypical charge transfer correlated insulators La2CuO4 and NiO, and obtain spectral properties and magnetic moments in good agreement with experiments.
  • The spin state transition in LaCoO$_3$ has eluded description for decades despite concerted theoretical and experimental effort. In this study, we approach this problem using charge consistent Density Functional Theory + Embedded Dynamical Mean Field Theory (DFT+DMFT). We show, from first principles, that LaCoO$_3$ cannot be described by a single, pure spin state at any temperature, but instead shows a gradual change in the population of higher spin multiples as temperature is increased. We explicitly elucidate the critical role of the lattice expansion and oxygen octahedral rotations in the spin state transition. We also show that the spin state transition and the metal-insulator transition in LaCoO$_3$ occur at different temperature scales. In addition, our results shed light on the importance of electronic entropy, which has so far been ignored in all first principles studies of this material.
  • We introduce a cluster DMFT (Dynamical Mean Field Theory) approach to study the normal state of the iron pnictides and chalcogenides. In the regime of moderate mass renormalizations, the self-energy is very local, justifying the success of single site DMFT for these materials and for other Hunds metals. We solve the corresponding impurity model with CTQMC (Continuous Time Quantum Monte-Carlo) and find that the minus sign problem is not severe in regimes of moderate mass renormalization.
  • We use neutron scattering to study spin excitations in single crystals of LiFe$_{0.88}$Co$_{0.12}$As, which is located near the boundary of the superconducting phase of LiFe$_{1-x}$Co$_{x}$As and exhibits non-Fermi-liquid behavior indicative of a quantum critical point. By comparing spin excitations of LiFe$_{0.88}$Co$_{0.12}$As with a combined density functional theory (DFT) and dynamical mean field theory (DMFT) calculation, we conclude that wave-vector correlated low energy spin excitations are mostly from the $d_{xy}$ orbitals, while high-energy spin excitations arise from the $d_{yz}$ and $d_{xz}$ orbitals. Unlike most iron pnictides, the strong orbital selective spin excitations in LiFeAs family cannot be described by anisotropic Heisenberg Hamiltonian. While the evolution of low-energy spin excitations of LiFe$_{1-x}$Co$_x$As are consistent with electron-hole Fermi surface nesting condition for the $d_{xy}$ orbital, the reduced superconductivity in LiFe$_{0.88}$Co$_{0.12}$As suggests that Fermi surface nesting conditions for the $d_{yz}$ and $d_{xz}$ orbitals are also important for superconductivity in iron pnictides.
  • We implemented the derivative of the free energy functional with respect to the atom displacements, so called force, within the combination of Density Functional Theory and the Embedded Dynamical Mean Field Theory. We show that in combination with the numerically exact quantum Monte Carlo (MC) impurity solver, the MC noise cancels to a great extend, so that the method can be used very efficiently for structural optimization of correlated electron materials. As an application of the method, we show how strengthening of the fluctuating moment in FeSe superconductor leads to a substantial increase of the anion height, and consequently to a very large effective mass, and also strong orbital differentiation.
  • Magnetoelectric (ME) effect is recognized for its utility for low-power electronic devices.Largest ME coefficients are often associated with phase transitions in which ferroelectricity is induced by magnetic order. Unfortunately, in these systems, large ME response is revealed only upon elaborate poling procedures. These procedures may become unnecessary in single-polar-domain crystals of polar magnets. Here we report giant ME effects in a polar magnet Fe2Mo3O8 at temperatures as high as 60 K.Polarization jumps of 0.3 {\mu}C/cm2, and repeated mutual control of ferroelectric and magnetic moments with differential ME coefficients on the order of 10$^4$ ps/m are achieved.Importantly, no electric or magnetic poling is needed, as necessary for applications.The sign of the ME coefficients can be switched by changing the applied "bias" magnetic field. The observed effects are associated with a hidden ferrimagnetic order unveiled by application of a magnetic field.
  • Dynamical mean field theory (DMFT) combined with the local density approximation (LDA) is widely used in solids to predict properties of correlated systems. In this paper, its application to one of the simplest strongly correlated systems, the hydrogen molecule H$_2$, is demonstrated to develop a parameter-free LDA+DMFT framework. We propose a method to calculate the exact intersection of LDA and DMFT that leads to highly accurate subtraction of the doubly counted correlation in both methods. When the exact double-counting treatment and a good projector to the correlated subspace are used, LDA+DMFT yields very accurate total energy and excitation spectrum of the H$_2$ molecule. We also discuss how this double-counting scheme can be extended to solid state calculations.
  • Combining density functional theory (DFT) and embedded dynamical mean-field theory (DMFT) methods, we study the metal-insulator transition in $R_2$Ir$_2$O$_7$ ($R$=Y, Eu, Sm, Nd, Pr, and Bi) and the topological nature of the insulating compounds. Accurate free energies evaluated using the charge self-consistent DFT+DMFT method reveal that the metal-insulator transition occurs for an A-cation radius between that of Nd and Pr, in agreement with experiments. The all-in-all-out magnetic phase, which is stable in the Nd compound but not the Pr one, gives rise to a small Ir$^{4+}$ magnetic moment of $\approx 0.5\mu_B$ and opens a sizable correlated gap. We demonstrate that within this state-of-the-art theoretical method, the insulating bulk pyrochlore iridates are topologically trivial.
  • We present a systematical theoretical study on the transport properties of an archetypal family of Hund's metals, Sr$_2$RuO$_4$,Sr$_3$Ru$_2$O$_7$, SrRuO$_3$ and CaRuO$_3$, within the combination of first principles density functional theory and dynamical mean field theory. The agreement between theory and experiments for optical conductivity and resistivity is good, which indicates that electron-electron scattering dominates the transport of ruthenates. We demonstrate that in the single-site dynamical mean field approach the transport properties of Hund's metals fall into the scenario of "resilient quasiparticles". We explains why the single layered compound Sr$_2$RuO$_4$ has a relative weak correlation with respect to its siblings, which corroborates its good metallicity.
  • Discovery of new transition metal compounds with large spin orbit coupling (SOC) coexisting with strong electron-electron correlation among the $d$ electrons is essential for understanding the physics that emerges from the interplay of these two effects. In this study, we predict a novel class of $J_{eff}=1/2$ Mott insulators in a family of fluoride compounds that are previously synthesized, but not characterized extensively. First principles calculations in the level of all electron Density Functional Theory + Dynamical Mean Field Theory (DFT+DMFT) indicate that these compounds have large Mott gaps and some of them exhibit unprecedented proximity to the ideal, $SU(2)$ symmetric $J_{eff}=1/2$ limit.
  • The stationary functional of the all-electron density functional plus dynamical mean field theory (DFT+DMFT) formalism to perform free energy calculations and structural relaxations is implemented for the first time. Here, the first order error in the density leads to a much smaller, second order error in the free energy. The method is applied to several well known correlated materials; metallic SrVO$_3$, Mott insulating FeO, and elemental Cerium, to show that it predicts the lattice constants with very high accuracy. In Cerium, we show that our method predicts the iso-structural transition between the $\alpha$ and $\gamma$ phases, and resolve the long standing controversy in the driving mechanism of this transition.
  • We propose a continuum representation of the Dynamical Mean Field Theory, in which we were able to derive an exact overlap between the Dynamical Mean Field Theory and band structure methods, such as the Density Functional Theory. The implementation of this exact double-counting shows improved agreement between theory and experiment in several correlated solids, such as the transition metal oxides and lanthanides. Previously introduced nominal double-counting is in much better agreement with the exact double-counting than most widely used fully localized limit formula.
  • Layered 5d transition metal dichalcogenide (TMD) IrTe2 is distinguished from the traditional TMDs (such as NbSe2) by the existence of multiple CDW-like stripe phases and superconductivity at low temperatures. Despite of intensive studies, there is still no consensus on the physical origin of the stripe phases or even the ground state modulation for this 5d material. Here, we present atomic-scale evidence from scanning tunneling microscopy and spectroscopy (STM/STS), that the ground state of IrTe2 is a q=1/6 stripe phase, identical to that of the Se-doped compound. Furthermore, our data suggest that the multiple transitions and stripe phases are driven by the intralayer Ir-Ir dimerization that competes against the interlayer Te-Te bonding. The competition results in a unified phase diagram with a series of hierarchical modulated stripe phases, strikingly similar to the renowned "devil's staircase" phenomena.
  • The solution of a generalized impurity model lies at the heart of electronic structure calculations with dynamical mean-field theory (DMFT). In the strongly-correlated regime, the method of choice for solving the impurity model is the hybridization expansion continuous time quantum Monte Carlo(CT-HYB). Enhancements to the CT-HYB algorithm are critical for bringing new physical regimes within reach of current computational power. Taking advantage of the fact that the bottleneck in the algorithm is a product of hundreds of matrices, we present optimizations based on the introduction and combination of two concepts of more general applicability: a) skip lists and b) fast rejection of proposed configurations based on matrix bounds. Considering two very different test cases with $d$ electrons, we find speedups of $\sim 25$ up to $\sim 500$ compared to the direct evaluation of the matrix product. Even larger speedups are likely with $f$ electron systems and with clusters of correlated atoms.
  • A combination of dynamical mean field theory and density functional theory, as implemented in Phys. Rev. B 81, 195107 (2010), is applied to both the early and late transition metal oxides. For fixed value of the local Coulomb repulsion, without fine tuning, we obtain the main features of these series, such as the metallic character of SrVO$_3$ and the the insulating gaps of LaVO$_3$, LaTiO$_3$ and La$_2$CO$_4$ which are in good agreement with experiment. The study highlights the importance of local physics and high energy hybridization in the screening of the Hubbard interaction and how different low energy behaviors can emerge from the a unified treatment of the transition metal series.
  • Using tight-binding models and first-principles calculations, we demonstrate the possibility to achieve a quantum anomalous Hall (QAH) phase on a two-dimensional square lattice, which can be realized in monolayers of double perovskites. We show that effective intersite spin-orbit coupling between eg orbitals can be induced perturbatively, giving rise to a QAH state. Moreover, the effective spin-orbit coupling can be enhanced by octahedral rotations. Based on first-principles calculations, we propose that this type of QAH state could be realized in La2MnIrO6 monolayers, with the size of the gap as large as 26 meV in the ideal case. We observe that the electronic structure is sensitive to structural distortions, and that an enhanced Hubbard U tends to stabilize the nontrivial gap.