Search

We study the influence of ferromagnetic nonlocal exchange on correlated electrons in terms of an SU(2)-Hubbard-Heisenberg model and address the interplay of on-site interaction induced local moment formation and the competition of ferromagnetic direct and antiferromagnetic kinetic exchange interactions. In order to simulate thermodynamic properties of the system in a way that largely accounts for the on-site interaction driven correlations in the system, we advance the correlated variational scheme introduced in [M. Sch\"uler et al., Phys. Rev. Lett. 111, 036601 (2013)] to account for explicitly symmetry broken electronic phases by introducing an auxiliary magnetic field. After benchmarking the method against exact solutions of a finite system, we study the SU(2) Hubbard-Heisenberg model on a square lattice. We obtain the $U\ensuremath{-}J$ finite temperature phase diagram of a SU(2)-Hubbard-Heisenberg model within the correlated variational approach and compare it to static mean-field theory. While the generalized variational principle and static mean-field theory yield transitions from dominant ferromagnetic to antiferromagnetic correlations in similar regions of the phase diagram, we find that the nature of the associated phase transitions differs between the two approaches. The fluctuations accounted for in the generalized variational approach render the transitions continuous, while static mean-field theory predicts discontinuous transitions between ferro- and antiferromagnetically ordered states.

While the Hubbard model is the standard model to study Mott metal-insulator transitions, it is still unclear to which extent it can describe metal-insulator transitions in real solids, where non-local Coulomb interactions are always present. By using a variational principle, we clarify this issue for short- and long-ranged non-local Coulomb interactions for half-filled systems on bipartite lattices. We find that repulsive non-local interactions generally stabilize the Fermi-liquid regime. The metal-insulator phase boundary is shifted to larger interaction strengths to leading order linearly with non-local interactions. Importantly, non-local interactions can raise the order of the metal-insulator transition. We present a detailed analysis of how the dimension and geometry of the lattice as well as the temperature determine the critical non-local interaction leading to a first-order transition: for systems in more than two dimensions with non-zero density of states at the Fermi energy the critical non-local interaction is arbitrarily small; otherwise it is finite., Comment: 10 pages, 11 figures

Using ab-initio calculations we reveal the nature of the insulating phase recently found experimentally in monolayer 1T-NbSe$_2$. We find soft phonon modes in a large parts of the Brillouin zone indicating the strong-coupling nature of a charge-density-wave instability. Structural relaxation of a $\sqrt{13}\times\sqrt{13}$ supercell reveals a Star-of-David reconstruction with an energy gain of 60 meV per primitive unit cell. The band structure of the distorted phase exhibits a half-filled flat band which is associated with orbitals that are delocalized over several atoms in each Star of David. By including many-body corrections through a combined GW, hybrid-functional, and DMFT treatment, we find the flat band to split into narrow Hubbard bands. The lowest energy excitation across the gap turns out to be between itinerant Se-$p$ states and the upper Hubbard band, determining the system to be a charge-transfer insulator. Combined hybrid-functional and GW calculations show that long-range interactions shift the Se-$p$ states to lower energies. Thus, a delicate interplay of local and long-range correlations determines the gap nature and its size in this distorted phase of the monolayer 1T-NbSe$_2$., Comment: 10 pages, 11 figures

In order for methods combining ab initio density-functional theory and many-body techniques to become routinely used, a flexible, fast, and easy-to-use implementation is crucial. We present an implementation of a general charge self-consistent scheme based on projected localized orbitals in the projector augmented wave framework in the Vienna Ab Initio Simulation Package (VASP). We give a detailed description on how the projectors are optimally chosen and how the total energy is calculated. We benchmark our implementation in combination with dynamical mean-field theory: first we study the charge-transfer insulator NiO using a Hartree-Fock approach to solve the many-body Hamiltonian. We address the advantages of the optimized against non-optimized projectors and furthermore find that charge self-consistency decreases the dependence of the spectral function - especially the gap - on the double counting. Second, using continuous-time quantum Monte Carlo we study a monolayer of SrVO$_3$, where strong orbital polarization occurs due to the reduced dimensionality. Using total-energy calculation for structure determination, we find that electronic correlations have a non-negligible influence on the position of the apical oxygens, and therefore on the thickness of the single SrVO$_3$ layer., 11 pages, 6 figures

Nanostructures with open shell transition metal or molecular constituents host often strong electronic correlations and are highly sensitive to atomistic material details. This tutorial review discusses method developments and applications of theoretical approaches for the realistic description of the electronic and magnetic properties of nanostructures with correlated electrons. First, the implementation of a flexible interface between density functional theory and a variant of dynamical mean field theory (DMFT) highly suitable for the simulation of complex correlated structures is explained and illustrated. On the DMFT side, this interface is largely based on recent developments of quantum Monte Carlo and exact diagonalization techniques allowing for efficient descriptions of general four fermion Coulomb interactions, reduced symmetries and spin-orbit coupling, which are explained here. With the examples of the Cr (001) surfaces, magnetic adatoms, and molecular systems it is shown how the interplay of Hubbard U and Hund's J determines charge and spin fluctuations and how these interactions drive different sorts of correlation effects in nanosystems. Non-local interactions and correlations present a particular challenge for the theory of low dimensional systems. We present our method developments addressing these two challenges, i.e., advancements of the dynamical vertex approximation and a combination of the constrained random phase approximation with continuum medium theories. We demonstrate how non-local interaction and correlation phenomena are controlled not only by dimensionality but also by coupling to the environment which is typically important for determining the physics of nanosystems., tutorial review submitted to EPJ-ST (scientific report of research unit FOR 1346); 14 figures, 26 pages

Graphene is an ideal platform to study many-body effects due to its semimetallic character and the possibility to dope it over a wide range. Here we study the width of graphene's occupied $\pi$-band as a function of doping using angle-resolved photoemission. Upon increasing electron doping, we observe the expected shift of the band to higher binding energies. However, this shift is not rigid and the bottom of the band moves less than the Dirac point. We show that the observed shift cannot be accounted for by band structure calculations in the local density approximation but that non-local exchange interactions must be taken into account., Comment: 5 pages, 3 figures

The electronic structure of the Cr(001) surface with its sharp resonance at the Fermi level is a subject of controversial debate of many experimental and theoretical works. To date, it is unclear whether the origin of this resonance is an orbital Kondo or an electron-phonon coupling effect. We have combined ab initio density functional calculations with dynamical mean-field simulations to calculate the orbitally resolved spectral function of the Cr(001) surface. The calculated orbital character and shape of the spectrum is in agreement with data from (inverse) photoemission experiments. We find that dynamic electron correlations crucially influence the surface electronic structure and lead to a low energy resonance in the $d_{z^2}$ and $d_{xz/yz}$ orbitals. Our results help to reconvene controversial experimental results from (I)PES and STM measurements., Comment: 8 pages, 5 figures

We investigate the Peierls-Feynman-Bogoliubov variational principle to map Hubbard models with nonlocal interactions to effective models with only local interactions. We study the renormalization of the local interaction induced by nearest-neighbor interaction and assess the quality of the effective Hubbard models in reproducing observables of the corresponding extended Hubbard models. We compare the renormalization of the local interactions as obtained from numerically exact determinant Quantum Monte Carlo to approximate but more generally applicable calculations using dual boson, dynamical mean field theory, and the random phase approximation. These more approximate approaches are crucial for any application with real materials in mind. Furthermore, we use the dual boson method to calculate observables of the extended Hubbard models directly and benchmark these against determinant Quantum Monte Carlo simulations of the effective Hubbard model., Comment: 12 pages, 10 figures

We describe a variational approach to solving Anderson impurity models by means of exact diagonalization. Optimized parameters of a discretized auxiliary model are obtained on the basis of the Peierls-Feynman-Bogoliubov principle. Thereby, the variational approach resolves ambiguities related with the bath discretization, which is generally necessary to make Anderson impurity models tractable by exact diagonalization. The choice of variational degrees of freedom made here allows systematic improvements of total energies over mean field decouplings like Hartree-Fock. Furthermore, our approach allows us to embed arbitrary bath discretization schemes in total energy calculations and to systematically optimize and improve on traditional routes to the discretization problem such as fitting of hybridization functions on Matsubara frequencies. Benchmarks in terms of a single orbital Anderson model demonstrate that the variational exact diagonalization method accurately reproduces free energies as well as several single- and two-particle observables obtained from an exact solution. Finally, we demonstrate the applicability of the variational exact diagonalization approach to realistic five orbital problems with the example system of Co impurities in bulk Cu and compare to continuous-time Monte Carlo calculations. The accuracy of established bath discretization schemes is assessed in the framework of the variational approach introduced here., Comment: 10 pages, 4 figures, accepted in PRB, v2: Added section on spectral information

A proper theoretical description of the electronic structure of the $3d$ orbitals in the metal centers of functional metalorganics is a challenging problem. We apply density functional theory and an exact diagonalization method in a many-body approach to study the ground-state electronic configuration of an iron porphyrin (FeP) molecule. Our study reveals that the consideration of multiple Slater determinants is important, and FeP is a potential candidate for realizing a spin crossover due to a subtle balance of crystal-field effects, on-site Coulomb repulsion, and hybridization between the Fe-$d$ orbitals and ligand N-$p$ states. The mechanism of switching between two close-lying electronic configurations of Fe-$d$ orbitals is shown. We discuss the generality of the suggested approach and the possibility to properly describe the electronic structure and related low-energy physics of the whole class of correlated metal-centered organometallic molecules.

To understand how nonlocal Coulomb interactions affect the phase diagram of correlated electron materials, we report on a method to approximate a correlated lattice model with nonlocal interactions by an effective Hubbard model with on-site interactions U* only. The effective model is defined by the Peierls-Feynman-Bogoliubov variational principle. We find that the local part of the interaction U is reduced according to U*=U-V', where V' is a weighted average of nonlocal interactions. For graphene, silicene and benzene we show that the nonlocal Coulomb interaction can decrease the effective local interaction by more than a factor of 2 in a wide doping range., Accepted for publication in PRL. Supplemental material added in this version