First-principles approach to the electronic structure of strongly correlated systems: combining the GW approximation and dynamical mean-field theory

Internal consistency of multi-tier G W+EDMFT

The multi-tier G W+EDMFT scheme is an ab-initio method for calculating the electronic structure of correlated materials. While the approach is free from ad-hoc parameters, it requires a selection of appropriate energy windows for describing low-energy and strongly correlated physics. In this study, we test the consistency of the multi-tier description by considering different low-energy windows for a series of cubic SrXO3 (X = V, Cr, Mn) perovskites. Specifically, we compare the 3-orbital t 2g model, the 5-orbital t 2g  + e g model, the 12-orbital t 2g + O p model, and (in the case of SrVO3) the 14-orbital t 2g  + e g  + O p model and compare the results to available photoemission and X-ray absorption measurements. The multi-tier method yields consistent results for the t 2g and t 2g  + e g low-energy windows, while the models with O p states produce stronger correlation effects and mostly agree well with experiment, especially in the unoccupied part of the spectrum. We also discuss the consistency between the fermionic and bosonic spectral functions and the physical origin of satellite features, and present momentum-resolved charge susceptibilities.

Challenges with relativistic GW calculations in solids and molecules

For molecules and solids containing heavy elements, accurate electronic-structure calculations require accounting not only for electronic correlations but also for relativistic effects. In molecules, relativity can lead to severe changes in the ground-state description. In solids, the interplay between both correlation and relativity can change the stability of phases or it can lead to an emergence of completely new phases. Traditionally, the simplest illustration of relativistic effects can be done either by including pseudopotentials in non-relativistic calculations or alternatively by employing large all-electron basis sets in relativistic methods. By analyzing different electronic properties (band structure, equilibrium lattice constant and bulk modulus) in semiconductors and insulators, we show that capturing the interplay of relativity and electron correlation can be rather challenging in Green's function methods. For molecular problems with heavy elements, we also observe that similar problems persist. We trace these challenges to three major problems: deficiencies in pseudopotential treatment as applied to Green's function methods, the scarcity of accurate and compact all-electron basis sets that can be converged with respect to the basis-set size, and linear dependencies arising in all-electron basis sets, particularly when employing Gaussian orbitals. Our analysis provides detailed insight into these problems and opens a discussion about potential approaches to mitigate them.

Restoring translational symmetry in periodic all-orbital dynamical mean-field theory simulations

Dynamical mean-field theory (DMFT) and its cluster extensions provide an efficient Green's function formalism to simulate spectral properties of periodic systems at the quantum many-body level. However, traditional cluster DMFT breaks translational invariance in solid-state materials, and the best strategy to capture non-local correlation effects within cluster DMFT remains elusive. In this work, we investigate the use of overlapping atom-centered impurity fragments in recently-developed ab initio all-orbital DMFT, where all local orbitals within the impurity are treated with high-level quantum chemistry impurity solvers. We demonstrate how the translational symmetry of the lattice self-energy can be restored by designing symmetry-adapted embedding problems, which results in an improved description of spectral functions in two-dimensional boron nitride monolayers and graphene at the levels of many-body perturbation theory (GW) and coupled-cluster theory. Furthermore, we study the convergence of self-energy and density of states as the embedding size is systematically expanded in one-shot and self-consistent DMFT calculations.

Mott-Insulator State of FeSe as a Van der Waals 2D Material Is Unveiled

We undertook a comprehensive investigation of the electronic structure of FeSe, known as a Hund metal, and found that it is not uniquely defined. Through accounting for all two-particle irreducible diagrams constructed from electron Green's function G and screened Coulomb interaction W in a self-consistent manner, a Mott-insulator phase of 2D-FeSe is unveiled. The metal-insulator transition is driven by the strong on-site Coulomb interaction in its paramagnetic phase, accompanied by the weakening of both local and nonlocal screening effects on the Fe-3d orbitals. Our results suggest that Mott physics may play a pivotal role in shaping the electronic, optical, and superconducting properties of monolayer or nanostructured FeSe.

Local and Nonlocal Electronic Correlations at the Metal-Insulator Transition in the Two-Dimensional Hubbard Model

Elucidating the physics of the single-orbital Hubbard model in its intermediate-coupling regime is a key missing ingredient to our understanding of metal-insulator transitions in real materials. Using recent nonperturbative many-body techniques that are able to interpolate between the spin-fluctuation-dominated Slater regime at weak coupling and the Mott insulator at strong coupling, we obtain the momentum-resolved spectral function in the intermediate regime and disentangle the effects of antiferromagnetic fluctuations and local electronic correlations in the formation of an insulating state. This allows us to identify the Slater and Heisenberg regimes in the phase diagram, which are separated by a crossover region of competing spatial and local electronic correlations. We identify the crossover regime by investigating the behavior of the local magnetic moment, shedding light on the formation of the insulating state at intermediate couplings.

Efficient Full-Frequency GW Calculations Using a Lanczos Method

The GW approximation is widely used for reliable and accurate modeling of single-particle excitations. It also serves as a starting point for many theoretical methods, such as its use in the Bethe-Salpeter equation (BSE) and dynamical mean-field theory. However, full-frequency GW calculations for large systems with hundreds of atoms remain computationally challenging, even after years of efforts to reduce the prefactor and improve scaling. We propose a method that reformulates the correlation part of the GW self-energy as a resolvent of a Hermitian matrix, which can be efficiently and accurately computed using the standard Lanczos method. This method enables full-frequency GW calculations of material systems with a few hundred atoms on a single computing workstation. We further demonstrate the efficiency of the method by calculating the defect-state energies of silicon quantum dots with diameters up to 4 nm and nearly 2,000 silicon atoms using only 20 computational nodes.

Correlated Electronic Structure of La_{3}Ni_{2}O_{7} under Pressure

Recently, superconductivity with a T_{c} up to 78 K has been reported in bulk samples of the bilayer nickelate La_{3}Ni_{2}O_{7} at pressures above 14 GPa. Important theoretical tasks are the formulation of relevant low-energy models and the clarification of the normal state properties. Here, we study the correlated electronic structure of the high-pressure phase in a four-orbital low-energy subspace using different many-body approaches: GW, dynamical mean field theory (DMFT), extended DMFT (EDMFT) and GW+EDMFT, with realistic frequency-dependent interaction parameters. The nonlocal correlation and screening effects captured by GW+EDMFT result in an instability toward the formation of charge stripes, with the 3d_{z^{2}} as the main active orbital. We also comment on the potential relevance of the rare-earth self-doping pocket, since hole doping suppresses the ordering tendency.

Coulomb engineering of two-dimensional Mott materials

Two-dimensional materials can be strongly influenced by their surroundings. A dielectric environment screens and reduces the Coulomb interaction between electrons in the two-dimensional material. Since in Mott materials the Coulomb interaction is responsible for the insulating state, manipulating the dielectric screening provides direct control over Mottness. Our many-body calculations reveal the spectroscopic fingerprints of such Coulomb engineering: we demonstrate eV-scale changes to the position of the Hubbard bands and show a Coulomb engineered insulator-to-metal transition. Based on our proof-of-principle calculations, we discuss the (feasible) conditions under which our scenario of Coulomb engineering of Mott materials can be realized experimentally.

First principles investigation of screened Coulomb interaction and electronic structure of low-temperature phase TaS2

By means of ab initio computation schemes, we examine the electronic screening, Coulomb interaction strength, and the electronic structure of a quantum spin liquid candidate monolayer TaS2 in its low-temperature commensurate charge-density-wave phase. Not only local (U) but non-local (V) correlations are estimated within random phase approximation based on two different screening models. Using GW + EDMFT (GW plus extended dynamical mean-field theory) method, we investigate the detailed electronic structure by increasing the level of non-local approximation from DMFT (V=0) to EDMFT and GW + EDMFT.

Bootstrap Embedding on a Quantum Computer

We extend molecular bootstrap embedding to make it appropriate for implementation on a quantum computer. This enables solution of the electronic structure problem of a large molecule as an optimization problem for a composite Lagrangian governing fragments of the total system, in such a way that fragment solutions can harness the capabilities of quantum computers. By employing state-of-art quantum subroutines including the quantum SWAP test and quantum amplitude amplification, we show how a quadratic speedup can be obtained over the classical algorithm, in principle. Utilization of quantum computation also allows the algorithm to match─at little additional computational cost─full density matrices at fragment boundaries, instead of being limited to 1-RDMs. Current quantum computers are small, but quantum bootstrap embedding provides a potentially generalizable strategy for harnessing such small machines through quantum fragment matching.

Breakdown of the scaling relation of anomalous Hall effect in Kondo lattice ferromagnet USbTe

The interaction between strong correlation and Berry curvature is an open territory of in the field of quantum materials. Here we report large anomalous Hall conductivity in a Kondo lattice ferromagnet USbTe which is dominated by intrinsic Berry curvature at low temperatures. However, the Berry curvature induced anomalous Hall effect does not follow the scaling relation derived from Fermi liquid theory. The onset of the Berry curvature contribution coincides with the Kondo coherent temperature. Combined with ARPES measurement and DMFT calculations, this strongly indicates that Berry curvature is hosted by the flat bands induced by Kondo hybridization at the Fermi level. Our results demonstrate that the Kondo coherence of the flat bands has a dramatic influence on the low temperature physical properties associated with the Berry curvature, calling for new theories of scaling relations of anomalous Hall effect to account for the interaction between strong correlation and Berry curvature.

Quantum embedding theories to simulate condensed systems on quantum computers

Quantum computers hold promise to improve the efficiency of quantum simulations of materials and to enable the investigation of systems and properties that are more complex than tractable at present on classical architectures. Here, we discuss computational frameworks to carry out electronic structure calculations of solids on noisy intermediate-scale quantum computers using embedding theories, and we give examples for a specific class of materials, that is, solid materials hosting spin defects. These are promising systems to build future quantum technologies, such as quantum computers, quantum sensors and quantum communication devices. Although quantum simulations on quantum architectures are in their infancy, promising results for realistic systems appear to be within reach.

Mott versus Hybridization Gap in the Low-Temperature Phase of 1T-TaS_{2}

We address the long-standing problem of the ground state of 1T-TaS_{2} by computing the correlated electronic structure of stacked bilayers using the GW+EDMFT method. Depending on the surface termination, the semi-infinite uncorrelated system is either band insulating or exhibits a metallic surface state. For realistic values of the on-site and inter-site interactions, a Mott gap opens in the surface state, but it is smaller than the gap originating from the bilayer structure. Our results are consistent with recent scanning tunneling spectroscopy measurements for different terminating layers, and with our own photoemission measurements, which indicate the coexistence of spatial regions with different gaps in the electronic spectrum. By comparison to exact diagonalization data, we clarify the interplay between Mott insulating and band insulating behavior in this archetypal layered system.

Simplified GW/BSE Approach for Charged and Neutral Excitation Energies of Large Molecules and Nanomaterials

Inspired by Grimme's simplified Tamm-Dancoff density functional theory approach [Grimme, S. J. Chem. Phys. 2013, 138, 244104], we describe a simplified approach to excited-state calculations within the GW approximation to the self-energy and the Bethe-Salpeter equation (BSE), which we call sGW/sBSE. The primary simplification to the electron repulsion integrals yields the same structure as with tensor hypercontraction, such that our method has a storage requirement that grows quadratically with system size and computational timing that grows cubically with system size. The performance of sGW is tested on the ionization potential of the molecules in the GW100 test set, for which it differs from ab initio GW calculations by only 0.2 eV. The performance of sBSE (based on the sGW input) is tested on the excitation energies of molecules in Thiel's set, for which it differs from ab initio GW/BSE calculations by about 0.5 eV. As examples of the systems that can be routinely studied with sGW/sBSE, we calculate the band gap and excitation energy of hydrogen-passivated silicon nanocrystals with up to 2650 electrons in 4678 spatial orbitals and the absorption spectra of two large organic dye molecules with hundreds of atoms.

The plain and simple parquet approximation: single-and multi-boson exchange in the two-dimensional Hubbard model

Abstract

The parquet approach to vertex corrections is unbiased but computationally demanding. Most applications are therefore restricted to small cluster sizes or rely on various simplifying approximations. We have recently shown that the bosonization of the parquet diagrams provides interpretative and algorithmic advantages over the original purely fermionic formulation. Here, we present first results of the numerical implementation of this method by applying it to the half-filled Hubbard model on the square lattice at weak coupling. The improved algorithmic performance allows us to evaluate the parquet approximation for a 16 × 16 lattice, retaining the full momentum and frequency structure of the various vertex functions. We discuss their symmetries and consider parametrizations of their momentum dependence using the truncated-unity approximation.

Graphical abstract

Exploring Coupled Cluster Green's Function as a Method for Treating System and Environment in Green's Function Embedding Methods

Within the self-energy embedding theory (SEET) framework, we study the coupled cluster Green's function (GFCC) method in two different contexts: as a method to treat either the system or the environment present in the embedding construction. Our study reveals that when GFCC is used to treat the environment we do not see improvement in total energies in comparison to the coupled cluster method itself. To rationalize this puzzling result, we analyze the performance of GFCC as an impurity solver with a series of transition metal oxides. These studies shed light on the strength and weaknesses of such a solver and demonstrate that such a solver gives very accurate results when the size of the impurity is small. We investigate if it is possible to achieve a systematic accuracy of the embedding solution when we increase the size of the impurity problem. We found that in such a case, the performance of the solver worsens, both in terms of finding the ground state solution of the impurity problem and the self-energies produced. We concluded that increasing the rank of GFCC solver is necessary to be able to enlarge impurity problems and achieve a reliable accuracy. We also have shown that natural orbitals from weakly correlated perturbative methods are better suited than symmetrized atomic orbitals (SAO) when the total energy of the system is the target quantity.

Designing and controlling the properties of transition metal oxide quantum materials

This Perspective addresses the design, creation, characterization and control of synthetic quantum materials with strong electronic correlations. We show how emerging synergies between theoretical/computational approaches and materials design/experimental probes are driving recent advances in the discovery, understanding and control of new electronic behaviour in materials systems with interesting and potentially technologically important properties. The focus here is on transition metal oxides, where electronic correlations lead to a myriad of functional properties including superconductivity, magnetism, Mott transitions, multiferroicity and emergent behaviour at picoscale-designed interfaces. Current opportunities and challenges are also addressed, including possible new discoveries of non-equilibrium phenomena and optical control of correlated quantum phases of transition metal oxides.

Spatial non-locality of electronic correlations beyond GW approximation

The question of spatial locality of electronic correlations beyond GW approximation is one of the central issues of the famous combination of GW and dynamical mean field theory, GW + DMFT. In this work, the above question is addressed directly (for the first time) by performing calculations with and without assumption of locality of the corresponding diagrams. For this purpose we use sc(GW + G3W2) approach where the higher order part (G3W2) is evaluated with fully momentum dependent Green's function G and screened interaction W and with 'local' variant, where the single site approximation is assumed for both G and W. For all three materials studied in this work (NiO,α-Ce, LiFeAs), we have found the spatial non-locality effects to be strong. For NiO and LiFeAs they, in fact, are decisive for the proper evaluation of vertex corrections. The results of this study have direct impact on our understanding of approximations made in practical implementations of GW + DMFT method, where all diagrams beyond GW (DMFT part) are assumed to be local. Taking into account the fact that the first diagrams beyond GW represent the most important contribution also in GW + DMFT calculations, we conclude that the basic assumption of GW + DMFT, namely the locality of diagrams evaluated in the DMFT part, is not as good as it is believed to be.

Fock-Space Embedding Theory: Application to Strongly Correlated Topological Phases

A many-body wave function can be factorized in Fock space into a marginal amplitude describing a set of strongly correlated orbitals and a conditional amplitude for the remaining weakly correlated part. The marginal amplitude is the solution of a Schrödinger equation with an effective Hamiltonian that can be viewed as embedding the marginal wave function in the environment of weakly correlated electrons. Here, the complementary equation for the conditional amplitude is replaced by a generalized Kohn-Sham equation, for which an orbital-dependent functional approximation is shown to reproduce the topological phase diagram of a multiband Hubbard model as a function of crystal field and Hubbard parameters. The roles of band filling and interband fluctuations are elucidated.

Fast Green's Function Method for Ultrafast Electron-Boson Dynamics

The interaction of electrons with quantized phonons and photons underlies the ultrafast dynamics of systems ranging from molecules to solids, and it gives rise to a plethora of physical phenomena experimentally accessible using time-resolved techniques. Green's function methods offer an invaluable interpretation tool since scattering mechanisms of growing complexity can be selectively incorporated in the theory. Currently, however, real-time Green's function simulations are either prohibitively expensive due to the cubic scaling with the propagation time or do neglect the feedback of electrons on the bosons, thus violating energy conservation. We put forward a computationally efficient Green's function scheme which overcomes both limitations. The numerical effort scales linearly with the propagation time while the simultaneous dressing of electrons and bosons guarantees the fulfillment of all fundamental conservation laws. We present a real-time study of the phonon-driven relaxation dynamics in an optically excited narrow band-gap insulator, highlighting the nonthermal behavior of the phononic degrees of freedom. Our formulation paves the way to first-principles simulations of electron-boson systems with unprecedented long propagation times.

Electronic-structure methods for materials design

The accuracy and efficiency of electronic-structure methods to understand, predict and design the properties of materials has driven a new paradigm in research. Simulations can greatly accelerate the identification, characterization and optimization of materials, with this acceleration driven by continuous progress in theory, algorithms and hardware, and by adaptation of concepts and tools from computer science. Nevertheless, the capability to identify and characterize materials relies on the predictive accuracy of the underlying physical descriptions, and on the ability to capture the complexity of realistic systems. We provide here an overview of electronic-structure methods, of their application to the prediction of materials properties, and of the different strategies employed towards the broader goals of materials design and discovery.

Advanced First-Principle Modeling of Relativistic Ruddlesden—Popper Strontium Iridates

In this review, we provide a survey of the application of advanced first-principle methods on the theoretical modeling and understanding of novel electronic, optical, and magnetic properties of the spin-orbit coupled Ruddlesden–Popper series of iridates Srn+1IrnO3n+1 (n = 1, 2, and ∞). After a brief description of the basic aspects of the adopted methods (noncollinear local spin density approximation plus an on-site Coulomb interaction (LSDA+U), constrained random phase approximation (cRPA), GW, and Bethe–Salpeter equation (BSE)), we present and discuss select results. We show that a detailed phase diagrams of the metal–insulator transition and magnetic phase transition can be constructed by inspecting the evolution of electronic and magnetic properties as a function of Hubbard U, spin–orbit coupling (SOC) strength, and dimensionality n, which provide clear evidence for the crucial role played by SOC and U in establishing a relativistic (Dirac) Mott–Hubbard insulating state in Sr2IrO4 and Sr3Ir2O7. To characterize the ground-state phases, we quantify the most relevant energy scales fully ab initio—crystal field energy, Hubbard U, and SOC constant of three compounds—and discuss the quasiparticle band structures in detail by comparing GW and LSDA+U data. We examine the different magnetic ground states of structurally similar n = 1 and n = 2 compounds and clarify that the origin of the in-plane canted antiferromagnetic (AFM) state of Sr2IrO4 arises from competition between isotropic exchange and Dzyaloshinskii–Moriya (DM) interactions whereas the collinear AFM state of Sr3Ir2O7 is due to strong interlayer magnetic coupling. Finally, we report the dimensionality controlled metal–insulator transition across the series by computing their optical transitions and conductivity spectra at the GW+BSE level from the the quasi two-dimensional insulating n = 1 and 2 phases to the three-dimensional metallic n=∞ phase.

Mott Metal-Insulator Transition from Steady-State Density Functional Theory

We present a computationally efficient method to obtain the spectral function of bulk systems in the framework of steady-state density functional theory (i-DFT) using an idealized scanning tunneling microscope (STM) setup. We calculate the current through the STM tip and then extract the spectral function from the finite-bias differential conductance. The fictitious noninteracting system of i-DFT features an exchange-correlation (XC) contribution to the bias which guarantees the same current as in the true interacting system. Exact properties of the XC bias are established using Fermi-liquid theory and subsequently implemented to construct approximations for the Hubbard model. We show for two different lattice structures that the Mott metal-insulator transition is captured by i-DFT.

Multiorbital Quantum Impurity Solver for General Interactions and Hybridizations

We present a numerically exact inchworm Monte Carlo method for equilibrium multiorbital quantum impurity problems with general interactions and hybridizations. We show that the method, originally developed to overcome the dynamical sign problem in certain real-time propagation problems, can also overcome the sign problem as a function of temperature for equilibrium quantum impurity models. This is shown in several cases where the current method of choice, the continuous-time hybridization expansion, fails due to the sign problem. Our method therefore enables simulations of impurity problems as they appear in embedding theories without further approximations, such as the truncation of the hybridization or interaction structure or a discretization of the impurity bath with a set of discrete energy levels, and eliminates a crucial bottleneck in the simulation of ab initio embedding problems.

Effect of Nonlocal Correlations on the Electronic Structure of LiFeAs

We investigate the role of nonlocal correlations in LiFeAs by exploring an ab initio-derived multiorbital Hubbard model for LiFeAs via the two-particle self-consistent (TPSC) approach. The multiorbital formulation of TPSC approximates the irreducible interaction vertex to be an orbital-dependent constant, which is self-consistently determined from local spin and charge sum rules. Within this approach, we disentangle the contribution of local and nonlocal correlations in LiFeAs and show that in the local approximation one recovers the dynamical mean field theory result. The comparison of our theoretical results to most recent angular-resolved photoemission spectroscopy and de Haas-van Alphen data shows that nonlocal correlations in LiFeAs are decisive to describe the measured spectral function A(k[over →],ω), Fermi surface, and scattering rates. These findings underline the importance of nonlocal correlations and benchmark different theoretical approaches for iron-based superconductors.

Bandwidth renormalization due to the intersite Coulomb interaction

The theory of correlated electrons is currently moving beyond the paradigmatic Hubbard U, towards the investigation of intersite Coulomb interactions. Recent investigations have revealed that these interactions are relevant for the quantitative description of realistic materials. Physically, intersite interactions are responsible for two rather different effects: screening and bandwidth renormalization. We use a variational principle to disentangle the roles of these two processes and study how appropriate the recently proposed Fock treatment of intersite interactions is in correlated systems. The magnitude of this effect in graphene is calculated based on cRPA values of the intersite interaction. We also apply the variational principle to benzene and find effective parameters comparable to those obtained by ab initio density matrix downfolding.

Magnetic force theory combined with quasi-particle self-consistent GW method

We report a successful combination of magnetic force linear response theory with quasiparticle self-consistent GW method. The self-consistently determined wavefunctions and eigenvalues can just be used for the conventional magnetic force calculations. While its formulation is straightforward, this combination provides a way to investigate the effect of GW self-energy on the magnetic interactions which can hardly be quantified due to the limitation of current GW methodology in calculating the total energy difference in between different magnetic phases. In ferromagnetic 3d elements, GW self-energy slightly reduces the d bandwidth and enhances the interactions while the same long-range feature is maintained. In antiferromagnetic transition-metal monoxides, QSGW significantly reduces the interaction strengths by enlarging the gap. Orbital-dependent magnetic force calculations show that the coupling between e _g and the nominally-empty 4s orbital is noticeably large in MnO which is reminiscent of the discussion for cuprates regarding the role of Cu-4s state. This combination of magnetic force theory with quasiparticle self-consistent GW can be a useful tool to study various magnetic materials.

Electronic correlations in dense iron: from moderate pressure to Earth's core conditions

We discuss the role of dynamical many-electron effects in the physics of iron and iron-rich solid alloys under applied pressure on the basis of recent ab initio studies employing the dynamical mean-field theory (DMFT). We review in detail two particularly interesting regimes: first, a moderate pressure range up to 60 GPa and, second, the ultra-high pressure of about 360 GPa expected inside the solid inner core of Earth. Electronic correlations in iron under the moderate pressure of several tens GPa are discussed in the first section. DMFT-based methods predict an enhancement of electronic correlations at the pressure-induced body-centered cubic α to hexagonal close-packed [Formula: see text] phase transition. In particular, the electronic effective mass, scattering rate and electron-electron contribution to the electrical resistivity undergo a step-wise increase at the transition point. One also finds a significant many-body correction to the [Formula: see text]-Fe equation of state, thus clarifying the origin of discrepancies between previous DFT studies and experiment. An electronic topological transition is predicted to be induced in [Formula: see text]-Fe by many-electron effects; its experimental signatures are analyzed. The next section focuses on the geophysically relevant pressure-temperature regime of the Earth's inner core (EIC) corresponding to the extreme pressure of 360 GPa combined with temperatures up to 6000 K. The three iron allotropes ([Formula: see text], [Formula: see text] and face-centered-cubic [Formula: see text]) previously proposed as possible stable phases at such conditions are found to exhibit qualitatively different many-electron effects as evidenced by a strongly non-Fermi-liquid metallic state of [Formula: see text]-Fe and an almost perfect Fermi liquid in the case of [Formula: see text]-Fe. A recent active discussion on the electronic state and transport properties of [Formula: see text]-Fe at the EIC conditions is reviewed in details. Estimations for the dynamical many-electron contribution to the relative phase stability are presented. We also discuss the impact of a Ni admixture, which is expected to be present in the core matter. We conclude by outlining some limitation of the present DMFT-based framework relevant for studies of iron-base systems as well as perspective directions for further development.

The GW Compendium: A Practical Guide to Theoretical Photoemission Spectroscopy

The GW approximation in electronic structure theory has become a widespread tool for predicting electronic excitations in chemical compounds and materials. In the realm of theoretical spectroscopy, the GW method provides access to charged excitations as measured in direct or inverse photoemission spectroscopy. The number of GW calculations in the past two decades has exploded with increased computing power and modern codes. The success of GW can be attributed to many factors: favorable scaling with respect to system size, a formal interpretation for charged excitation energies, the importance of dynamical screening in real systems, and its practical combination with other theories. In this review, we provide an overview of these formal and practical considerations. We expand, in detail, on the choices presented to the scientist performing GW calculations for the first time. We also give an introduction to the many-body theory behind GW, a review of modern applications like molecules and surfaces, and a perspective on methods which go beyond conventional GW calculations. This review addresses chemists, physicists and material scientists with an interest in theoretical spectroscopy. It is intended for newcomers to GW calculations but can also serve as an alternative perspective for experts and an up-to-date source of computational techniques.

Correlated materials design: prospects and challenges

The design of correlated materials challenges researchers to combine the maturing, high throughput framework of DFT-based materials design with the rapidly-developing first-principles theory for correlated electron systems. We review the field of correlated materials, distinguishing two broad classes of correlation effects, static and dynamics, and describe methodologies to take them into account. We introduce a material design workflow, and illustrate it via examples in several materials classes, including superconductors, charge ordering materials and systems near an electronically driven metal to insulator transition, highlighting the interplay between theory and experiment with a view towards finding new materials. We review the statistical formulation of the errors of currently available methods to estimate formation energies. We formulate an approach for estimating a lower-bound for the probability of a new compound to form. Correlation effects have to be considered in all the material design steps. These include bridging between structure and property, obtaining the correct structure and predicting material stability. We introduce a post-processing strategy to take them into account.

Toward a predictive theory of correlated materials

Correlated electron materials display a rich variety of notable properties ranging from unconventional superconductivity to metal-insulator transitions. These properties are of interest from the point of view of applications but are hard to treat theoretically, as they result from multiple competing energy scales. Although possible in more weakly correlated materials, theoretical design and spectroscopy of strongly correlated electron materials have been a difficult challenge for many years. By treating all the relevant energy scales with sufficient accuracy, complementary advances in Green's functions and quantum Monte Carlo methods open a path to first-principles computational property predictions in this class of materials.

Effective Heisenberg Model and Exchange Interaction for Strongly Correlated Systems

We consider the extended Hubbard model and introduce a corresponding Heisenberg-like problem written in terms of spin operators. The derived formalism is reminiscent of Anderson's idea of the effective exchange interaction and takes into account nonlocal correlation effects. The results for the exchange interaction and spin susceptibility in the magnetic phase are expressed in terms of single-particle quantities. This fact not only can be used for realistic calculations of multiband systems but also allows us to reconsider a general description of many-body effects in the most interesting physical regimes, where the physical properties of the system are dominated by collective (bosonic) fluctuations. In the strongly spin-polarized limit, when the local magnetic moment is well defined, the exchange interaction reduces to a standard expression of the density functional theory that has been successfully used in practical calculations of magnetic properties of real materials.

Thermoelectricity in correlated narrow-gap semiconductors

We review many-body effects, their microscopic origin, as well as their impact on thermoelectricity in correlated narrow-gap semiconductors. Members of this class-such as FeSi and FeSb₂-display an unusual temperature dependence in various observables: insulating with large thermopowers at low temperatures, they turn bad metals at temperatures much smaller than the size of their gaps. This insulator-to-metal crossover is accompanied by spectral weight-transfers over large energies in the optical conductivity and by a gradual transition from activated to Curie-Weiss-like behaviour in the magnetic susceptibility. We show a retrospective of the understanding of these phenomena, discuss the relation to heavy-fermion Kondo insulators-such as Ce₃Bi₄Pt₃ for which we present new results-and propose a general classification of paramagnetic insulators. From the latter, FeSi emerges as an orbital-selective Kondo insulator. Focussing on intermetallics such as silicides, antimonides, skutterudites, and Heusler compounds we showcase successes and challenges for the realistic simulation of transport properties in the presence of electronic correlations. Further, we explore new avenues in which electronic correlations may contribute to the improvement of thermoelectric performance.

Variational structure of Luttinger-Ward formalism and bold diagrammatic expansion for Euclidean lattice field theory

The Luttinger-Ward functional was proposed more than five decades ago and has been used to formally justify most practically used Green's function methods for quantum many-body systems. Nonetheless, the very existence of the Luttinger-Ward functional has been challenged by recent theoretical and numerical evidence. We provide a rigorously justified Luttinger-Ward formalism, in the context of Euclidean lattice field theory. Using the Luttinger-Ward functional, the free energy can be variationally minimized with respect to Green's functions in its domain. We then derive the widely used bold diagrammatic expansion rigorously, without relying on formal arguments such as partial resummation of bare diagrams to infinite order.

Fierz Convergence Criterion: A Controlled Approach to Strongly Interacting Systems with Small Embedded Clusters

We present an embedded-cluster method, based on the triply irreducible local expansion formalism. It turns the Fierz ambiguity, inherent to approaches based on a bosonic decoupling of local fermionic interactions, into a convergence criterion. It is based on the approximation of the three-leg vertex by a coarse-grained vertex computed from a self-consistently determined cluster impurity model. The computed self-energies are, by construction, continuous functions of momentum. We show that, in three interaction and doping regimes of the two-dimensional Hubbard model, self-energies obtained with clusters of size four only are very close to numerically exact benchmark results. We show that the Fierz parameter, which parametrizes the freedom in the Hubbard-Stratonovich decoupling, can be used as a quality control parameter. By contrast, the GW+extended dynamical mean field theory approximation with four cluster sites is shown to yield good results only in the weak-coupling regime and for a particular decoupling. Finally, we show that the vertex has spatially nonlocal components only at low Matsubara frequencies.

Breakdown of Traditional Many-Body Theories for Correlated Electrons

Starting from the (Hubbard) model of an atom, we demonstrate that the uniqueness of the mapping from the interacting to the noninteracting Green function, G→G_{0}, is strongly violated, by providing numerous explicit examples of different G_{0} leading to the same physical G. We argue that there are indeed infinitely many such G_{0}, with numerous crossings with the physical solution. We show that this rich functional structure is directly related to the divergence of certain classes of (irreducible vertex) diagrams, with important consequences for traditional many-body physics based on diagrammatic expansions. Physically, we ascribe the onset of these highly nonperturbative manifestations to the progressive suppression of the charge susceptibility induced by the formation of local magnetic moments and/or resonating valence bond (RVB) states in strongly correlated electron systems.

Local magnetic moments in iron and nickel at ambient and Earth's core conditions

Some Bravais lattices have a particular geometry that can slow down the motion of Bloch electrons by pre-localization due to the band-structure properties. Another known source of electronic localization in solids is the Coulomb repulsion in partially filled d or f orbitals, which leads to the formation of local magnetic moments. The combination of these two effects is usually considered of little relevance to strongly correlated materials. Here we show that it represents, instead, the underlying physical mechanism in two of the most important ferromagnets: nickel and iron. In nickel, the van Hove singularity has an unexpected impact on the magnetism. As a result, the electron-electron scattering rate is linear in temperature, in violation of the conventional Landau theory of metals. This is true even at Earth's core pressures, at which iron is instead a good Fermi liquid. The importance of nickel in models of geomagnetism may have therefore to be reconsidered.

Coulomb correlations in 4d and 5d oxides from first principles-or how spin-orbit materials choose their effective orbital degeneracies

The interplay of spin-orbit coupling and Coulomb correlations has become a hot topic in condensed matter theory and is especially important in 4d and 5d transition metal oxides, like iridates or rhodates. Here, we review recent advances in dynamical mean-field theory (DMFT)-based electronic structure calculations for treating such compounds, introducing all necessary implementation details. We also discuss the evaluation of Hubbard interactions in spin-orbit materials. As an example, we perform DMFT calculations on insulating strontium iridate (Sr₂IrO₄) and its 4d metallic counterpart, strontium rhodate (Sr₂RhO₄). While a Mott-insulating state is obtained for Sr₂IrO₄ in its paramagnetic phase, the spectral properties and Fermi surfaces obtained for Sr₂RhO₄ show excellent agreement with available experimental data. Finally, we discuss the electronic structure of these two compounds by introducing the notion of effective spin-orbital degeneracy as the key quantity that determines the correlation strength. We stress that effective spin-orbital degeneracy introduces an additional axis into the conventional picture of a phase diagram based on filling and on the ratio of interactions to bandwidth, analogous to the degeneracy-controlled Mott transition in d¹ perovskites.

Nonequilibrium GW+EDMFT: Antiscreening and Inverted Populations from Nonlocal Correlations

We study the dynamics of screening in photodoped Mott insulators with long-ranged interactions using a nonequilibrium implementation of the GW plus extended dynamical mean-field theory formalism. Our study demonstrates that the complex interplay of the injected carriers with bosonic degrees of freedom (charge fluctuations) can result in long-lived transient states with properties that are distinctly different from those of thermal equilibrium states. Systems with strong nonlocal interactions are found to exhibit a self-sustained population inversion of the doublons and holes. This population inversion leads to low-energy antiscreening which can be detected in time-resolved electron-energy-loss spectra.

Electronic Phase Separation and Dramatic Inverse Band Renormalization in the Mixed-Valence Cuprate LiCu_{2}O_{2}

We measured, by angle-resolved photoemission spectroscopy, the electronic structure of LiCu_{2}O_{2}, a mixed-valence cuprate where planes of Cu(I) (3d^{10}) ions are sandwiched between layers containing one-dimensional edge-sharing Cu(II) (3d^{9}) chains. We find that the Cu(I)- and Cu(II)-derived electronic states form separate electronic subsystems, in spite of being coupled by bridging O ions. The valence band, of the Cu(I) character, disperses within the charge-transfer gap of the strongly correlated Cu(II) states, displaying an unprecedented 250% broadening of the bandwidth with respect to the predictions of density functional theory. Our observation is at odds with the widely accepted tenet of many-body theory that correlation effects generally yield narrower bands and larger electron masses and suggests that present-day electronic structure techniques provide an intrinsically inappropriate description of ligand-to-d hybridizations in late transition metal oxides.

Theoretical calculation and experiment study on the electronic structure, microstructures and photocatalytic activity of N–Al codoped TiO2

To more accurately demonstrate how the electronic structure was changed and to explain the origin of the enhanced photocatalytic activity by N–Al codoped TiO₂, we investigated N, Al monodoped, and codoped TiO₂by theoretical calculation and experiments.

Dynamical screening in correlated electron systems-from lattice models to realistic materials

Recent progress in treating the dynamical nature of the screened Coulomb interaction in strongly correlated lattice models and materials is reviewed with a focus on computational schemes based on the dynamical mean field approximation. We discuss approximate and exact methods for the solution of impurity models with retarded interactions, and explain how these models appear as auxiliary problems in various extensions of the dynamical mean field formalism. The current state of the field is illustrated with results from recent applications of these schemes to U-V Hubbard models and correlated materials.