Showing 1–38 of 38 results for author: van Loon, E G C P

Second-order phase transitions and divergent linear response in dynamical mean-field theory

Authors: Erik G. C. P. van Loon

Abstract: Second-order phase transitions appear as a divergence in one of the linear response functions. For a system of correlated electrons, the relevant divergent response can and does involve many-particle observables, most famously the double occupancy. Generally, evaluating the linear response function of many-particle observables requires a many-particle generalization of the Bethe-Salpeter equation.… ▽ More Second-order phase transitions appear as a divergence in one of the linear response functions. For a system of correlated electrons, the relevant divergent response can and does involve many-particle observables, most famously the double occupancy. Generally, evaluating the linear response function of many-particle observables requires a many-particle generalization of the Bethe-Salpeter equation. However, here I show that the divergence of linear response functions in dynamical mean-field theory is governed by a two-particle Bethe-Salpeter equation, even for many-particle observables. The reason for this is that the divergence at the second-order phase transition is produced by the self-consistent feedback of the dynamical mean-field. △ Less

Submitted 14 June, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Journal ref: Phys. Rev. B 109, L241110 (2024)

Downfolding from Ab Initio to Interacting Model Hamiltonians: Comprehensive Analysis and Benchmarking of the DFT+cRPA Approach

Authors: Yueqing Chang, Erik G. C. P. van Loon, Brandon Eskridge, Brian Busemeyer, Miguel A. Morales, Cyrus E. Dreyer, Andrew J. Millis, Shiwei Zhang, Tim O. Wehling, Lucas K. Wagner, Malte Rösner

Abstract: Model Hamiltonians are regularly derived from first-principles data to describe correlated matter. However, the standard methods for this contain a number of largely unexplored approximations. For a strongly correlated impurity model system, here we carefully compare a standard downfolding technique with the best possible ground-truth estimates for charge-neutral excited state energies and wavefun… ▽ More Model Hamiltonians are regularly derived from first-principles data to describe correlated matter. However, the standard methods for this contain a number of largely unexplored approximations. For a strongly correlated impurity model system, here we carefully compare a standard downfolding technique with the best possible ground-truth estimates for charge-neutral excited state energies and wavefunctions using state-of-the-art first-principles many-body wave function approaches. To this end, we use the vanadocene molecule and analyze all downfolding aspects, including the Hamiltonian form, target basis, double counting correction, and Coulomb interaction screening models. We find that the choice of target-space basis functions emerges as a key factor for the quality of the downfolded results, while orbital-dependent double counting correction diminishes the quality. Background screening to the Coulomb interaction matrix elements primarily affects crystal-field excitations. Our benchmark uncovers the relative importance of each downfolding step and offers insights into the potential accuracy of minimal downfolded model Hamiltonians △ Less

Submitted 8 July, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: 15 pages (+8 pages Supplemental Material), 10 figures

Journal ref: npj Computational Materials volume 10, Article number: 129 (2024)

Unconventional charge-density-wave gap in monolayer NbS$_2$

Authors: Timo Knispel, Jan Berges, Arne Schobert, Erik G. C. P. van Loon, Wouter Jolie, Tim Wehling, Thomas Michely, Jeison Fischer

Abstract: Using scanning tunneling microscopy and spectroscopy, for a monolayer of transition metal dichalcogenide H-NbS$_2$ grown by molecular beam epitaxy on graphene, we provide unambiguous evidence for a charge density wave (CDW) with a 3$\times$3 superstructure, which is not present in bulk NbS$_2$. Local spectroscopy displays a pronounced gap of the order of 20 meV at the Fermi level. Within the gap l… ▽ More Using scanning tunneling microscopy and spectroscopy, for a monolayer of transition metal dichalcogenide H-NbS$_2$ grown by molecular beam epitaxy on graphene, we provide unambiguous evidence for a charge density wave (CDW) with a 3$\times$3 superstructure, which is not present in bulk NbS$_2$. Local spectroscopy displays a pronounced gap of the order of 20 meV at the Fermi level. Within the gap low energy features are present. The gap structure with its low energy features is at variance with the expectation for a gap opening in the electronic band structure due to a CDW. Instead, comparison with \it{ab initio} calculations indicates that the observed gap structure must be attributed to combined electron-phonon quasiparticles. The phonons in question are the elusive amplitude and phase collective modes of the CDW transition. Our findings advance the understanding of CDW mechanisms in two dimensional materials and their spectroscopic signatures. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Dual Bethe-Salpeter equation for the multi-orbital lattice susceptibility within dynamical mean-field theory

Authors: Erik G. C. P. van Loon, Hugo U. R. Strand

Abstract: Dynamical mean-field theory describes the impact of strong local correlation effects in many-electron systems. While the single-particle spectral function is directly obtained within the formalism, two-particle susceptibilities can also be obtained by solving the Bethe-Salpeter equation. The solution requires handling infinite matrices in Matsubara frequency space. This is commonly treated using a… ▽ More Dynamical mean-field theory describes the impact of strong local correlation effects in many-electron systems. While the single-particle spectral function is directly obtained within the formalism, two-particle susceptibilities can also be obtained by solving the Bethe-Salpeter equation. The solution requires handling infinite matrices in Matsubara frequency space. This is commonly treated using a finite frequency cut-off, resulting in slow linear convergence. A decomposition of the two-particle response in local and non-local contributions enables a reformulation of the Bethe-Salpeter equation inspired by the dual boson formalism. The re-formulation has a drastically improved cubic convergence with respect to the frequency cut-off, facilitating the calculation of susceptibilities in multi-orbital systems considerably. This improved convergence arises from the fact that local contributions can be measured in the impurity solver. The dual Bethe-Salpeter equation uses the fully reducible vertex which is free from vertex divergences. We benchmark the approach on several systems including the spin susceptibility of strontium ruthenate Sr$_2$RuO$_4$, a strongly correlated Hund's metal with three active orbitals. △ Less

Submitted 15 May, 2024; v1 submitted 8 June, 2023; originally announced June 2023.

Journal ref: Phys. Rev. B 109, 155157 (2024)

Nb$_3$Cl$_8$: A Prototypical Layered Mott-Hubbard Insulator

Authors: Sergii Grytsiuk, Mikhail I. Katsnelson, Erik G. C. P. van Loon, Malte Rösner

Abstract: The Hubbard model provides an idealized description of electronic correlations in solids. Despite its simplicity, the model features a competition between several different phases that have made it one of the most studied systems in theoretical physics. Real materials usually deviate from the ideal of the Hubbard model in several ways, but the monolayer of Nb$_3$Cl$_8$ has recently appeared as a p… ▽ More The Hubbard model provides an idealized description of electronic correlations in solids. Despite its simplicity, the model features a competition between several different phases that have made it one of the most studied systems in theoretical physics. Real materials usually deviate from the ideal of the Hubbard model in several ways, but the monolayer of Nb$_3$Cl$_8$ has recently appeared as a potentially optimal candidate for the realization of such a single-orbital Hubbard model. Here we show how this single orbital Hubbard model can be indeed constructed within a "molecular" rather than atomic basis set using ab initio constrained random phase approximation calculations. This way, we provide the essential ingredients to connect experimental reality with ab initio material descriptions and correlated electron theory, which clarifies that monolayer Nb$_3$Cl$_8$ is a Mott insulator with a gap of about 1 to 1.2eV depending on its dielectric environment. By comparing with an atomistic three-orbital model, we show that the single molecular orbital description is indeed adequate. We furthermore comment on the expected electronic and magnetic structure of the compound and show that the Mott insulating state survives in the low-temperature and bulk phases of the material. △ Less

Submitted 8 May, 2023; originally announced May 2023.

Comments: 14 pages, 7 figures, 4 tables

Ab initio electron-lattice downfolding: potential energy landscapes, anharmonicity, and molecular dynamics in charge density wave materials

Authors: Arne Schobert, Jan Berges, Erik G. C. P. van Loon, Michael A. Sentef, Sergey Brener, Mariana Rossi, Tim O. Wehling

Abstract: The interplay of electronic and nuclear degrees of freedom presents an outstanding problem in condensed matter physics and chemistry. Computational challenges arise especially for large systems, long time scales, in nonequilibrium, or in systems with strong correlations. In this work, we show how downfolding approaches facilitate complexity reduction on the electronic side and thereby boost the si… ▽ More The interplay of electronic and nuclear degrees of freedom presents an outstanding problem in condensed matter physics and chemistry. Computational challenges arise especially for large systems, long time scales, in nonequilibrium, or in systems with strong correlations. In this work, we show how downfolding approaches facilitate complexity reduction on the electronic side and thereby boost the simulation of electronic properties and nuclear motion - in particular molecular dynamics (MD) simulations. Three different downfolding strategies based on constraining, unscreening, and combinations thereof are benchmarked against full density functional calculations for selected charge density wave (CDW) systems, namely 1H-TaS$_2$, 1T-TiSe$_2$, 1H-NbS$_2$, and a one-dimensional carbon chain. We find that the downfolded models can reproduce potential energy surfaces on supercells accurately and facilitate computational speedup in MD simulations by about five orders of magnitude in comparison to purely ab initio calculations. For monolayer 1H-TaS$_2$ we report classical replica exchange and quantum path integral MD simulations, revealing the impact of thermal and quantum fluctuations on the CDW transition. △ Less

Submitted 16 January, 2024; v1 submitted 13 March, 2023; originally announced March 2023.

Journal ref: SciPost Phys. 16, 046 (2024)

Larmor precession in strongly correlated itinerant electron systems

Authors: Erik G. C. P. van Loon, Hugo U. R. Strand

Abstract: Many-electron systems undergo a collective Larmor precession in the presence of a magnetic field. In a paramagnetic metal, the resulting spin wave provides insight into the correlation effects generated by the electron-electron interaction. Here, we use dynamical mean-field theory to investigate the collective Larmor precession in the strongly correlated regime, where dynamical correlation effects… ▽ More Many-electron systems undergo a collective Larmor precession in the presence of a magnetic field. In a paramagnetic metal, the resulting spin wave provides insight into the correlation effects generated by the electron-electron interaction. Here, we use dynamical mean-field theory to investigate the collective Larmor precession in the strongly correlated regime, where dynamical correlation effects such as quasiparticle lifetimes and non-quasiparticle states are essential. We study the spin excitation spectrum, which includes a dispersive Larmor mode as well as electron-hole excitations that lead to Stoner damping. We also extract the momentum-resolved damping of slow spin waves. The accurate theoretical description of these phenomena relies on the Ward identity, which guarantees a precise cancellation of self-energy and vertex corrections at long wavelengths. Our findings pave the way towards a better understanding of spin wave damping in correlated materials. △ Less

Submitted 16 October, 2023; v1 submitted 6 March, 2023; originally announced March 2023.

Journal ref: Commun Phys 6, 289 (2023)

Two-particle correlations and the metal-insulator transition: Iterated Perturbation Theory revisited

Authors: Erik G. C. P. van Loon

Abstract: Recent advances in many-body physics have made it possible to study correlated electron systems at the two-particle level. In Dynamical Mean-Field theory, it has been shown that the metal-insulator phase diagram is closely related to the eigenstructure of the susceptibility. So far, this situation has been studied using accurate but numerically expensive solvers. Here, the Iterated Perturbation Th… ▽ More Recent advances in many-body physics have made it possible to study correlated electron systems at the two-particle level. In Dynamical Mean-Field theory, it has been shown that the metal-insulator phase diagram is closely related to the eigenstructure of the susceptibility. So far, this situation has been studied using accurate but numerically expensive solvers. Here, the Iterated Perturbation Theory (IPT) approximation is used instead. Its simplicity makes it possible to obtain analytical results for the two-particle vertex and the DMFT Jacobian. The limited computational cost also enables a detailed comparison of analytical expressions for the response functions to results obtained using finite differences. At the same time, the approximate nature of IPT precludes an interpretation of the metal-insulator transition in terms of a Landau free energy functional. △ Less

Submitted 31 May, 2022; v1 submitted 21 October, 2021; originally announced October 2021.

Comments: Revised version

Journal ref: Phys. Rev. B 105, 245104 (2022)

Local Plaquette Physics as Key Ingredient of High-Temperature Superconductivity in Cuprates

Authors: M. Danilov, E. G. C. P. van Loon, S. Brener, S. Iskakov, M. I. Katsnelson, A. I. Lichtenstein

Abstract: A major pathway towards understanding complex systems is given by exactly solvable reference systems that contain the essential physics of the system. For the $t-t'-U$ Hubbard model, the four-site plaquette is known to have a quantum critical point in the $U-μ$ space where states with electron occupations $N=2, 3, 4$ per plaquette are degenerate [Phys. Rev. B {\bf 94}, 125133 (2016)]. We show that… ▽ More A major pathway towards understanding complex systems is given by exactly solvable reference systems that contain the essential physics of the system. For the $t-t'-U$ Hubbard model, the four-site plaquette is known to have a quantum critical point in the $U-μ$ space where states with electron occupations $N=2, 3, 4$ per plaquette are degenerate [Phys. Rev. B {\bf 94}, 125133 (2016)]. We show that such a critical point in the lattice causes an instability in the particle-particle singlet d-wave channel and manifests some of the essential elements of the cuprate superconductivity. For this purpose we design an efficient superperturbation theory -- based on the dual fermion approach -- with the critical plaquette as the reference system. Thus, the perturbation theory already contains the relevant d-wave fluctuations from the beginning via the two-particle correlations of the plaquette. We find that d-wave superconductivity remains a leading instability channel under reasonably broad range of parameters. The next-nearest-neighbour hopping $t'$ is shown to play a crucial role in a formation of strongly bound electronic bipolarons whose coherence at lower temperature results in superconductivity. The physics of the pseudogap within the developed picture is also discussed. △ Less

Submitted 27 August, 2021; v1 submitted 23 July, 2021; originally announced July 2021.

Comments: 25 pages, 28 figures

Journal ref: npj Quantum Materials 7, 50 (2022)

Random Phase Approximation for gapped systems: role of vertex corrections and applicability of the constrained random phase approximation

Authors: Erik G. C. P. van Loon, Malte Rösner, Mikhail I. Katsnelson, Tim O. Wehling

Abstract: The many-body theory of interacting electrons poses an intrinsically difficult problem that requires simplifying assumptions. For the determination of electronic screening properties of the Coulomb interaction, the Random Phase Approximation (RPA) provides such a simplification. Here, we explicitly show that this approximation is justified for band structures with sizeable band gaps. This is when… ▽ More The many-body theory of interacting electrons poses an intrinsically difficult problem that requires simplifying assumptions. For the determination of electronic screening properties of the Coulomb interaction, the Random Phase Approximation (RPA) provides such a simplification. Here, we explicitly show that this approximation is justified for band structures with sizeable band gaps. This is when the electronic states responsible for the screening are energetically far away from the Fermi level, which is equivalent to a short electronic propagation length of these states. The RPA contains exactly those diagrams in which the classical Coulomb interaction covers all distances, whereas neglected vertex corrections involve quantum tunneling through the barrier formed by the band gap. Our analysis of electron-electron interactions provides a real-space analogy to Migdal's theorem on the smallness of vertex corrections in electron-phonon problems. An important application is the increasing use of constrained Random Phase Approximation (cRPA) calculations of effective interactions. We find that their usage of Kohn-Sham energies already accounts for the leading local (excitonic) vertex correction in insulators. △ Less

Submitted 23 July, 2021; v1 submitted 7 March, 2021; originally announced March 2021.

Journal ref: Phys. Rev. B 104, 045134 (2021)

Downfolding approaches to electron-ion coupling: Constrained density-functional perturbation theory for molecules

Authors: Erik G. C. P. van Loon, Jan Berges, Tim O. Wehling

Abstract: Constrained electronic-structure theories enable the construction of effective low-energy models consisting of partially dressed particles. However, the interpretation and physical content of these theories is not straightforward. Here, we carefully explore the properties of downfolding theories for electron-ion problems, in particular constrained density-functional perturbation theory (cDFPT). We… ▽ More Constrained electronic-structure theories enable the construction of effective low-energy models consisting of partially dressed particles. However, the interpretation and physical content of these theories is not straightforward. Here, we carefully explore the properties of downfolding theories for electron-ion problems, in particular constrained density-functional perturbation theory (cDFPT). We show that the dipole selection rules determine whether the partially dressed phonons satisfy Goldstone's theorem, and we prove that electronic screening always lowers the phonon frequencies. We illustrate the theory with cDFPT calculations for minimal example systems: the nitrogen and benzene molecule as well as graphene. △ Less

Submitted 10 May, 2021; v1 submitted 19 February, 2021; originally announced February 2021.

Comments: Revised manuscript with new appendix and ancillary scripts and data

Journal ref: Phys. Rev. B 103, 205103 (2021)

Efficient fluctuation exchange approach to low-temperature spin fluctuations and superconductivity: from the Hubbard model to Na$_x$CoO$_2\cdot y$H$_2$O

Authors: Niklas Witt, Erik G. C. P. van Loon, Takuya Nomoto, Ryotaro Arita, Tim Wehling

Abstract: Superconductivity arises mostly at energy and temperature scales that are much smaller than the typical bare electronic energies. Since the computational effort of diagrammatic many-body techniques increases with the number of required Matsubara frequencies and thus with the inverse temperature, phase transitions that occur at low temperatures are typically hard to address numerically. In this wor… ▽ More Superconductivity arises mostly at energy and temperature scales that are much smaller than the typical bare electronic energies. Since the computational effort of diagrammatic many-body techniques increases with the number of required Matsubara frequencies and thus with the inverse temperature, phase transitions that occur at low temperatures are typically hard to address numerically. In this work, we implement a fluctuation exchange (FLEX) approach to spin fluctuations and superconductivity using the "intermediate representation basis" (IR) [Shinaoka et al., PRB 96, 2017] for Matsubara Green functions. This FLEX+IR approach is numerically very efficient and enables us to reach temperatures on the order of $10^{-4}$ in units of the electronic band width in multi-orbital systems. After benchmarking the method in the doped repulsive Hubbard model on the square lattice, we study the possibility of spin-fluctuation-mediated superconductivity in the hydrated sodium cobalt material Na$_x$CoO$_2\cdot y$H$_2$O reaching the scale of the experimental transition temperature $T_{\mathrm{c}}=4.5$ K and below. △ Less

Submitted 26 May, 2021; v1 submitted 8 December, 2020; originally announced December 2020.

Comments: 13 pages, 8 figures, includes last revisions to the manuscript

Journal ref: Phys. Rev. B 103, 205148 (2021)

Second-order dual fermion for multi-orbital systems

Authors: Erik G. C. P. van Loon

Abstract: In dynamical mean-field theory, the correlations between electrons are assumed to be purely local. The dual fermion approach provides a systematic way of adding non-local corrections to the dynamical mean-field theory starting point. Initial applications of this method were largely restricted to the single-orbital Hubbard model. Here, we present an implementation of second-order dual fermion for g… ▽ More In dynamical mean-field theory, the correlations between electrons are assumed to be purely local. The dual fermion approach provides a systematic way of adding non-local corrections to the dynamical mean-field theory starting point. Initial applications of this method were largely restricted to the single-orbital Hubbard model. Here, we present an implementation of second-order dual fermion for general multi-orbital systems and use this approach to investigate spatial correlations in SrVO$_3$. In addition, the approach is benchmarked in several exactly solvable small systems. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: Example scripts are included in the arxiv source

Journal ref: J. Phys.: Condens. Matter 33 135601 (2021)

The Bethe-Salpeter equation at the critical end-point of the Mott transition

Authors: Erik G. C. P. van Loon, Friedrich Krien, Andrey A. Katanin

Abstract: Strong repulsive interactions between electrons can lead to a Mott metal-insulator transition. The Dynamical Mean-Field Theory (DMFT) explains the critical end-point and hysteresis region with single-particle concepts such as the spectral function and the quasiparticle weight. In this work, we reconsider the critical end point of the metal-insulator transition on the two-particle level. We show th… ▽ More Strong repulsive interactions between electrons can lead to a Mott metal-insulator transition. The Dynamical Mean-Field Theory (DMFT) explains the critical end-point and hysteresis region with single-particle concepts such as the spectral function and the quasiparticle weight. In this work, we reconsider the critical end point of the metal-insulator transition on the two-particle level. We show that the relevant eigenvalue and eigenvector of the non-local Bethe-Salpeter kernel provide a unified picture of the hysteresis region and of the critical end point of the Mott transition. In particular, they simultaneously explain the thermodynamics of the hysteresis region and the iterative stability of the DMFT equations. This analysis paves the way for a deeper understanding of phase transitions in correlated materials. △ Less

Submitted 24 September, 2020; v1 submitted 28 February, 2020; originally announced February 2020.

Comments: 5 pages + supplement

Journal ref: Phys. Rev. Lett. 125, 136402 (2020)

Coulomb Engineering of two-dimensional Mott materials

Authors: Erik G. C. P. van Loon, Malte Schüler, Daniel Springer, Giorgio Sangiovanni, Jan M. Tomczak, Tim O. Wehling

Abstract: 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… ▽ More 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. △ Less

Submitted 12 October, 2022; v1 submitted 6 January, 2020; originally announced January 2020.

Journal ref: npj 2D Mater Appl 7, 47 (2023)

Ab initio phonon self-energies and fluctuation diagnostics of phonon anomalies: Lattice instabilities from Dirac pseudospin physics in transition metal dichalcogenides

Authors: Jan Berges, Erik G. C. P. van Loon, Arne Schobert, Malte Rösner, Tim O. Wehling

Abstract: We present an ab initio approach for the calculation of phonon self-energies and their fluctuation diagnostics, which allows us to identify the electronic processes behind phonon anomalies. Application to the transition-metal-dichalcogenide monolayer 1H-TaS$_2$ reveals that coupling between the longitudinal-acoustic phonons and the electrons from an isolated low-energy metallic band is entirely re… ▽ More We present an ab initio approach for the calculation of phonon self-energies and their fluctuation diagnostics, which allows us to identify the electronic processes behind phonon anomalies. Application to the transition-metal-dichalcogenide monolayer 1H-TaS$_2$ reveals that coupling between the longitudinal-acoustic phonons and the electrons from an isolated low-energy metallic band is entirely responsible for phonon anomalies such as the mode softening and associated charge-density waves observed in this material. Our analysis allows us to distinguish between different mode-softening mechanisms including matrix-element effects, Fermi-surface nesting, and Van Hove scenarios. We find that matrix-element effects originating from a peculiar type of Dirac pseudospin textures control the charge-density-wave physics in 1H-TaS$_2$ and similar transition metal dichalcogenides. △ Less

Submitted 3 November, 2021; v1 submitted 6 November, 2019; originally announced November 2019.

Comments: Revised manuscript with ancillary scripts and data, 15 pages, 9 figures

Journal ref: Phys. Rev. B 101, 155107 (2020)

Thermodynamics of the metal-insulator transition in the extended Hubbard model

Authors: M. Schüler, E. G. C. P. van Loon, M. I. Katsnelson, T. O. Wehling

Abstract: In contrast to the Hubbard model, the extended Hubbard model, which additionally accounts for non-local interactions, lacks systemic studies of thermodynamic properties especially across the metal-insulator transition. Using a variational principle, we perform such a systematic study and describe how non-local interactions screen local correlations differently in the Fermi-liquid and in the insula… ▽ More In contrast to the Hubbard model, the extended Hubbard model, which additionally accounts for non-local interactions, lacks systemic studies of thermodynamic properties especially across the metal-insulator transition. Using a variational principle, we perform such a systematic study and describe how non-local interactions screen local correlations differently in the Fermi-liquid and in the insulator. The thermodynamics reveal that non-local interactions are at least in parts responsible for first-order metal-insulator transitions in real materials. △ Less

Submitted 29 May, 2019; v1 submitted 24 March, 2019; originally announced March 2019.

Comments: 19 pages, 9 figures. Resubmission to SciPost

Journal ref: SciPost Phys. 6, 067 (2019)

Dual Boson approach with instantaneous interaction

Authors: L. Peters, E. G. C. P. van Loon, A. N. Rubtsov, A. I. Lichtenstein, M. I. Katsnelson, E. A. Stepanov

Abstract: The Dual Boson approach to strongly correlated systems generally involves a dynamic (frequency-dependent) interaction in the auxiliary impurity model. In this work, we explore the consequences of forcing this interaction to be instantaneous (frequency-independent) via the use of a self-consistency condition on the instantaneous susceptibility. The result is a substantial simplification of the impu… ▽ More The Dual Boson approach to strongly correlated systems generally involves a dynamic (frequency-dependent) interaction in the auxiliary impurity model. In this work, we explore the consequences of forcing this interaction to be instantaneous (frequency-independent) via the use of a self-consistency condition on the instantaneous susceptibility. The result is a substantial simplification of the impurity model, especially with an eye on realistic multiband implementations, while keeping desireable properties of the Dual Boson approach, such as the charge conservation law, intact. We show and illustrate numerically that this condition enforces the absence of phase transitions in finite systems, as should be expected from general physical considerations, and respects the Mermin-Wagner theorem. In particular, the theory does not allow the metal to insulator phase transition associated with the formation of the magnetic order in the two-dimensional system. At the same time, the metal to charge ordered phase transition is allowed, as it is not associated with the spontaneous breaking of a continuous symmetry, and is accurately captured by the introduced approach. △ Less

Submitted 25 September, 2019; v1 submitted 18 February, 2019; originally announced February 2019.

Journal ref: Phys. Rev. B 100, 165128 (2019)

Bandwidth renormalization due to the intersite Coulomb interaction

Authors: Yann in 't Veld, Malte Schüler, Tim Wehling, Mikhail I. Katsnelson, Erik G. C. P. van Loon

Abstract: 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 renormal… ▽ More 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 observe that the most interesting charge fluctuation phenomena actually occur at elevated temperatures, substantially higher than studied in previous investigations. △ Less

Submitted 28 March, 2019; v1 submitted 31 January, 2019; originally announced January 2019.

Comments: New appendix on benzene

Journal ref: J. Phys.: Condens. Matter 31 465603 (2019)

Two-particle Fermi liquid parameters at the Mott transition: Vertex divergences, Landau parameters, and incoherent response in dynamical mean-field theory

Authors: Friedrich Krien, Erik G. C. P. van Loon, Mikhail I. Katsnelson, Alexander I. Lichtenstein, Massimo Capone

Abstract: We consider the interaction-driven Mott transition at zero temperature from the viewpoint of microscopic Fermi liquid theory. To this end, we derive an exact expression for the Landau parameters within the dynamical mean-field theory (DMFT) approximation to the single-band Hubbard model. At the Mott transition the symmetric and the anti-symmetric Landau parameter diverge. The vanishing compressibi… ▽ More We consider the interaction-driven Mott transition at zero temperature from the viewpoint of microscopic Fermi liquid theory. To this end, we derive an exact expression for the Landau parameters within the dynamical mean-field theory (DMFT) approximation to the single-band Hubbard model. At the Mott transition the symmetric and the anti-symmetric Landau parameter diverge. The vanishing compressibility at the Mott transition directly implies the divergence of the forward scattering amplitude in the charge sector, which connects the proximity of the Mott phase to a tendency towards phase separation. We verify the expected behavior of the Landau parameters in a DMFT application to the Hubbard model on the triangular lattice at finite temperature. Exact conservation laws and the Ward identity are crucial to capture vertex divergences related to the Mott transition. We furthermore generalize Leggett's formula for the static susceptibility of the Fermi liquid to the static fermion-boson response function. In the charge sector the limits of small transferred momentum and frequency of this response function commute at the Mott transition. △ Less

Submitted 16 June, 2019; v1 submitted 1 November, 2018; originally announced November 2018.

Comments: Only 17 pages (+ 6 pages Appendix)

Journal ref: Phys. Rev. B 99, 245128 (2019)

Fermion-boson vertex within Dynamical Mean-Field Theory

Authors: Erik G. C. P. van Loon, Friedrich Krien, Hartmut Hafermann, Alexander I. Lichtenstein, Mikhail I. Katsnelson

Abstract: In the study of strongly interacting systems, correlations on the two-particle level are receiving more and more attention. In this work, we study a particular two-particle correlation function: the fermion-boson vertex. It describes the response of the Green's function when an external field is applied and is an important ingredient of diagrammatic extensions of Dynamical Mean-Field Theory. We pr… ▽ More In the study of strongly interacting systems, correlations on the two-particle level are receiving more and more attention. In this work, we study a particular two-particle correlation function: the fermion-boson vertex. It describes the response of the Green's function when an external field is applied and is an important ingredient of diagrammatic extensions of Dynamical Mean-Field Theory. We provide several perspectives on this object, using Ward identities, sum rules, perturbative analysis and asymptotic relations. We then use these tools to study the vertex across the doping-driven metal-insulator transition and find a divergence in the imaginary part. △ Less

Submitted 5 November, 2018; v1 submitted 27 June, 2018; originally announced June 2018.

Comments: 21 pages

Journal ref: Phys. Rev. B 98, 205148 (2018)

Second-order dual fermion approach to the Mott transition in the two-dimensional Hubbard model

Authors: Erik G. C. P. van Loon, Mikhail I. Katsnelson, Hartmut Hafermann

Abstract: We apply the dual fermion approach with a second-order approximation to the self-energy to the Mott transition in the two-dimensional Hubbard model. The approximation captures nonlocal dynamical short-range correlations as well as several features observed in studies using cluster dynamical mean-field theory. This includes a strong reduction of the critical interaction and inversion of the slope o… ▽ More We apply the dual fermion approach with a second-order approximation to the self-energy to the Mott transition in the two-dimensional Hubbard model. The approximation captures nonlocal dynamical short-range correlations as well as several features observed in studies using cluster dynamical mean-field theory. This includes a strong reduction of the critical interaction and inversion of the slope of the transition lines with respect to single-site dynamical mean-field theory. We show that these effects coincide with a much smaller momentum differentiation compared to cluster methods. We further discuss the role of the self-consistency condition and show that the approximation behaves as an asymptotic series at low temperature. △ Less

Submitted 8 October, 2018; v1 submitted 22 May, 2018; originally announced May 2018.

Comments: 11 pages, 12 figures

Journal ref: Phys. Rev. B 98, 155117 (2018)

Precursors of the insulating state in the square lattice Hubbard model

Authors: Erik G. C. P. van Loon, Hartmut Hafermann, Mikhail Katsnelson

Abstract: We study the two-dimensional square lattice Hubbard model for small to moderate interaction strengths $1\leq U/t\leq 4$ by means of the ladder dual fermion approach. The non-local correlations beyond dynamical mean-field theory lower the potential energy, lead to a maximum in the uniform susceptibility and induce a pseudogap in the density of states. While the self-energy exhibits precursors of a… ▽ More We study the two-dimensional square lattice Hubbard model for small to moderate interaction strengths $1\leq U/t\leq 4$ by means of the ladder dual fermion approach. The non-local correlations beyond dynamical mean-field theory lower the potential energy, lead to a maximum in the uniform susceptibility and induce a pseudogap in the density of states. While the self-energy exhibits precursors of a possible insulating phase linked to the appearance of long-range fluctuations, the metallic phase persists within the accessible temperature range. Finite-size effects affect results qualitatively. Upper bounds on the crossover temperature are found to be significantly lower than previously reported dynamical vertex approximation results at $U/t=1$. △ Less

Submitted 22 December, 2017; originally announced December 2017.

Journal ref: Phys. Rev. B 97, 085125 (2018)

The extended Hubbard model with attractive interactions

Authors: E. G. C. P. van Loon, M. I. Katsnelson

Abstract: We discuss the phase diagram of the extended Hubbard model with both attractive and repulsive local and nonlocal interactions. The extended dynamical mean-field theory (EDMFT) and the dual boson method (DB) are compared. The latter contains additional nonlocal correlation effects that are not incorporated in EDMFT. We find that EDMFT and DB give almost identical results in the attractive $V$ regim… ▽ More We discuss the phase diagram of the extended Hubbard model with both attractive and repulsive local and nonlocal interactions. The extended dynamical mean-field theory (EDMFT) and the dual boson method (DB) are compared. The latter contains additional nonlocal correlation effects that are not incorporated in EDMFT. We find that EDMFT and DB give almost identical results in the attractive $V$ regime, where phase separation occurs. This is quite a difference with the previously studied repulsive $V$ regime, where EDMFT and DB give very different phase boundaries for the checkerboard order phase, especially at small $U$. △ Less

Submitted 26 September, 2017; v1 submitted 19 September, 2017; originally announced September 2017.

Comments: To be submitted to the proceedings of CCP2017, to appear in Journal of Physics Conference Series (IOP); comments are welcome

Journal ref: J. Phys.: Conf. Ser., 1136, 012006 (2018)

Competing Coulomb and electron-phonon interactions in NbS$_2$

Authors: E. G. C. P. van Loon, M. Rösner, G. Schönhoff, M. I. Katsnelson, T. O. Wehling

Abstract: The interplay of Coulomb and electron-phonon interactions with thermal and quantum fluctuations facilitates rich phase diagrams in two-dimensional electron systems. Layered transition metal dichalcogenides hosting charge, excitonic, spin and superconducting order form an epitomic material class in this respect. Theoretical studies of materials like NbS$_2$ have focused on the electron-phonon coupl… ▽ More The interplay of Coulomb and electron-phonon interactions with thermal and quantum fluctuations facilitates rich phase diagrams in two-dimensional electron systems. Layered transition metal dichalcogenides hosting charge, excitonic, spin and superconducting order form an epitomic material class in this respect. Theoretical studies of materials like NbS$_2$ have focused on the electron-phonon coupling whereas the Coulomb interaction, particularly strong in the monolayer limit, remained essentially untouched. Here, we analyze the interplay of short- and long-range Coulomb as well as electron-phonon interactions in NbS$_2$ monolayers. The combination of these interactions causes electronic correlations that are fundamentally different to what would be expected from the interaction terms separately. The fully interacting electronic spectral function resembles the non-interacting band structure but with appreciable broadening. An unexpected coexistence of strong charge and spin fluctuations puts NbS$_2$ close to spin and charge order, suggesting monolayer NbS$_2$ as a platform for atomic scale engineering of electronic quantum phases. △ Less

Submitted 2 July, 2018; v1 submitted 18 July, 2017; originally announced July 2017.

Journal ref: npj Quantum Materials 3, 32 (2018)

Conservation in two-particle self-consistent extensions of dynamical-mean-field-theory

Authors: F. Krien, E. G. C. P. van Loon, H. Hafermann, J. Otsuki, M. I. Katsnelson, A. I. Lichtenstein

Abstract: Extensions of dynamical-mean-field-theory (DMFT) make use of quantum impurity models as non-perturbative and exactly solvable reference systems which are essential to treat the strong electronic correlations. Through the introduction of retarded interactions on the impurity, these approximations can be made two-particle self-consistent. This is of interest for the Hubbard model, because it allows… ▽ More Extensions of dynamical-mean-field-theory (DMFT) make use of quantum impurity models as non-perturbative and exactly solvable reference systems which are essential to treat the strong electronic correlations. Through the introduction of retarded interactions on the impurity, these approximations can be made two-particle self-consistent. This is of interest for the Hubbard model, because it allows to suppress the antiferromagnetic phase transition in two-dimensions in accordance with the Mermin-Wagner theorem, and to include the effects of bosonic fluctuations. For a physically sound description of the latter, the approximation should be conserving. In this paper we show that the mutual requirements of two-particle self-consistency and conservation lead to fundamental problems. For an approximation that is two-particle self-consistent in the charge- and longitudinal spin channel, the double occupancy of the lattice and the impurity are no longer consistent when computed from single-particle properties. For the case of self-consistency in the charge- and longitudinal as well as transversal spin channels, these requirements are even mutually exclusive so that no conserving approximation can exist. We illustrate these findings for a two-particle self-consistent and conserving DMFT approximation. △ Less

Submitted 22 September, 2017; v1 submitted 30 June, 2017; originally announced June 2017.

Comments: 17 pages, 9 figures

Journal ref: Phys. Rev. B 96, 075155 (2017)

First-order metal-insulator transitions in the extended Hubbard model due to self-consistent screening of the effective interaction

Authors: M. Schüler, E. G. C. P. van Loon, M. I. Katsnelson, T. O. Wehling

Abstract: 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 lattice… ▽ More 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. △ Less

Submitted 24 April, 2018; v1 submitted 29 June, 2017; originally announced June 2017.

Comments: 10 pages, 11 figures

Journal ref: Phys. Rev. B 97, 165135 (2018)

A comparison between methods of analytical continuation for bosonic functions

Authors: Johan Schött, Erik G. C. P. van Loon, Inka L. M. Locht, Mikhail Katsnelson, Igor Di Marco

Abstract: In this article we perform a critical assessment of different known methods for the analytical continuation of bosonic functions, namely the maximum entropy method, the non-negative least-square method, the non-negative Tikhonov method, the Padé approximant method, and a stochastic sampling method. Three functions of different shape are investigated, corresponding to three physically relevant scen… ▽ More In this article we perform a critical assessment of different known methods for the analytical continuation of bosonic functions, namely the maximum entropy method, the non-negative least-square method, the non-negative Tikhonov method, the Padé approximant method, and a stochastic sampling method. Three functions of different shape are investigated, corresponding to three physically relevant scenarios. They include a simple two-pole model function and two flavours of the non-interacting Hubbard model on a square lattice, i.e. a single-orbital metallic system and a two-orbitals insulating system. The effect of numerical noise in the input data on the analytical continuation is discussed in detail. Overall, the stochastic method by Mishchenko et al. [Phys. Rev. B \textbf{62}, 6317 (2000)] is shown to be the most reliable tool for input data whose numerical precision is not known. For high precision input data, this approach is slightly outperformed by the Padé approximant method, which combines a good resolution power with a good numerical stability. Although none of the methods retrieves all features in the spectra in the presence of noise, our analysis provides a useful guideline for obtaining reliable information of the spectral function in cases of practical interest. △ Less

Submitted 14 July, 2016; originally announced July 2016.

Comments: 13 pages, 9 figures

Journal ref: Phys. Rev. B 94, 245140 (2016)

Capturing non-local interaction effects in the Hubbard model: optimal mappings and limits of applicability

Authors: E. G. C. P. van Loon, M. Schüler, M. I. Katsnelson, T. O. Wehling

Abstract: 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… ▽ More 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. △ Less

Submitted 18 October, 2016; v1 submitted 30 May, 2016; originally announced May 2016.

Comments: 12 pages, 10 figures

Journal ref: Phys. Rev. B 94, 165141 (2016)

From local to non-local correlations: the Dual Boson perspective

Authors: E. A. Stepanov, A. Huber, E. G. C. P. van Loon, A. I. Lichtenstein, M. I. Katsnelson

Abstract: Extended dynamical mean-field theory (EDMFT) is insufficient to describe non-local effects in strongly correlated systems, since corrections to the mean-field solution are generally large. We present an efficient scheme for the construction of diagrammatic extensions of EDMFT that avoids usual double counting problem by using an exact change of variables (the dual boson formalism) to distinguish t… ▽ More Extended dynamical mean-field theory (EDMFT) is insufficient to describe non-local effects in strongly correlated systems, since corrections to the mean-field solution are generally large. We present an efficient scheme for the construction of diagrammatic extensions of EDMFT that avoids usual double counting problem by using an exact change of variables (the dual boson formalism) to distinguish the correlations included in the mean-field solution and those beyond. With a computational efficiency comparable to EDMFT+GW approach, our scheme significantly improves on the charge order transition phase boundary in the extended Hubbard model. △ Less

Submitted 30 May, 2020; v1 submitted 26 April, 2016; originally announced April 2016.

Comments: Misprint in equation (A5) appeared in the previous version of the text has been corrected

Journal ref: Phys. Rev. B 94, 205110 (2016)

Interaction-driven Lifshitz transition with dipolar fermions in optical lattices

Authors: E. G. C. P. van Loon, M. I. Katsnelson, L. Chomaz, M. Lemeshko

Abstract: Anisotropic dipole-dipole interactions between ultracold dipolar fermions break the symmetry of the Fermi surface and thereby deform it. Here we demonstrate that such a Fermi surface deformation induces a topological phase transition -- so-called Lifshitz transition -- in the regime accessible to present-day experiments. We describe the impact of the Lifshitz transition on observable quantities su… ▽ More Anisotropic dipole-dipole interactions between ultracold dipolar fermions break the symmetry of the Fermi surface and thereby deform it. Here we demonstrate that such a Fermi surface deformation induces a topological phase transition -- so-called Lifshitz transition -- in the regime accessible to present-day experiments. We describe the impact of the Lifshitz transition on observable quantities such as the Fermi surface topology, the density-density correlation function, and the excitation spectrum of the system. The Lifshitz transition in ultracold atoms can be controlled by tuning the dipole orientation and -- in contrast to the transition studied in crystalline solids -- is completely interaction-driven. △ Less

Submitted 24 May, 2016; v1 submitted 30 March, 2016; originally announced March 2016.

Comments: 5 pages, 5 figures + 2 appendices

Journal ref: Phys. Rev. B 93, 195145 (2016)

Double occupancy in DMFT and the Dual Boson approach

Authors: Erik G. C. P. van Loon, Friedrich Krien, Hartmut Hafermann, Evgeny A. Stepanov, Alexander I. Lichtenstein, Mikhail I. Katsnelson

Abstract: We discuss the calculation of the double occupancy using Dynamical Mean-Field Theory (DMFT) in finite dimensions. The double occupancy can be determined from the susceptibility of the auxiliary impurity model or from the lattice susceptibility. The former method underestimates, whereas the latter overestimates the double occupancy. We illustrate this for the square-lattice Hubbard model. We propos… ▽ More We discuss the calculation of the double occupancy using Dynamical Mean-Field Theory (DMFT) in finite dimensions. The double occupancy can be determined from the susceptibility of the auxiliary impurity model or from the lattice susceptibility. The former method underestimates, whereas the latter overestimates the double occupancy. We illustrate this for the square-lattice Hubbard model. We propose an approach for which both methods lead to identical results by construction and which resolves this ambiguity. This self-consistent dual boson scheme results in a double occupancy that is numerically close to benchmarks available in the literature. △ Less

Submitted 28 April, 2016; v1 submitted 29 February, 2016; originally announced February 2016.

Comments: 5 pages, 5 figures

Journal ref: Phys. Rev. B 93, 155162 (2016)

Self-consistent Dual Boson approach to single-particle and collective excitations in correlated systems

Authors: E. A. Stepanov, E. G. C. P. van Loon, A. A. Katanin, A. I. Lichtenstein, M. I. Katsnelson, A. N. Rubtsov

Abstract: We propose an efficient dual boson scheme, which extends the DMFT paradigm to collective excitations in correlated systems. The theory is fully self-consistent both on the one- and on the two-particle level, thus describing the formation of collective modes as well as the renormalization of electronic and bosonic spectra on equal footing. The method employs an effective impurity model comprising b… ▽ More We propose an efficient dual boson scheme, which extends the DMFT paradigm to collective excitations in correlated systems. The theory is fully self-consistent both on the one- and on the two-particle level, thus describing the formation of collective modes as well as the renormalization of electronic and bosonic spectra on equal footing. The method employs an effective impurity model comprising both fermionic and bosonic hybridization functions. Only single- and two-electron Green's functions of the reference problem enter the theory, due to the optimal choice of the self-consistency condition for the effective bosonic bath. We show that the theory is naturally described by a dual Luttinger-Ward functional and obeys the relevant conservation laws. △ Less

Submitted 11 January, 2016; v1 submitted 28 August, 2015; originally announced August 2015.

Comments: 17 pages, 12 figures

Journal ref: Phys. Rev. B 93, 045107 (2016)

Ultralong-range order in the Fermi-Hubbard model with long-range interactions

Authors: Erik G. C. P. van Loon, Mikhail I. Katsnelson, Mikhail Lemeshko

Abstract: We use the dual boson approach to reveal the phase diagram of the Fermi-Hubbard model with long-range dipole-dipole interactions. By using a large-scale finite-temperature calculation on a $64 \times 64$ square lattice we demonstrate the existence of a novel phase, possessing an `ultralong-range' order. The fingerprint of this phase -- the density correlation function -- features a non-trivial beh… ▽ More We use the dual boson approach to reveal the phase diagram of the Fermi-Hubbard model with long-range dipole-dipole interactions. By using a large-scale finite-temperature calculation on a $64 \times 64$ square lattice we demonstrate the existence of a novel phase, possessing an `ultralong-range' order. The fingerprint of this phase -- the density correlation function -- features a non-trivial behavior on a scale of tens of the lattice sites. We study the properties and the stability of the ultralong-range ordered phase, and show that it is accessible in modern experiments with ultracold polar molecules and magnetic atoms. △ Less

Submitted 12 August, 2015; v1 submitted 19 June, 2015; originally announced June 2015.

Journal ref: Phys. Rev. B 92, 081106 (2015)

Thermodynamic consistency of the charge response in dynamical mean-field based approaches

Authors: Erik G. C. P. van Loon, Hartmut Hafermann, Alexander I. Lichtenstein, Mikhail I. Katsnelson

Abstract: We consider the thermodynamic consistency of the charge response function in the (extended) Hubbard model. In DMFT, thermodynamic consistency is preserved. We prove that the static, homogeneous DMFT susceptibility is consistent as long as vertex corrections obtained from the two-particle impurity correlation function are included. In presence of a nonlocal interaction, the problem may be treated w… ▽ More We consider the thermodynamic consistency of the charge response function in the (extended) Hubbard model. In DMFT, thermodynamic consistency is preserved. We prove that the static, homogeneous DMFT susceptibility is consistent as long as vertex corrections obtained from the two-particle impurity correlation function are included. In presence of a nonlocal interaction, the problem may be treated within extended DMFT (EDMFT), or its diagrammatic extension, the dual boson approach. We show that here, maintaining thermodynamic consistency requires knowledge of three- and four-particle impurity correlation functions, which are typically neglected. Nevertheless, the dual boson approximation to the response is remarkably close to consistency. This holds even when two-particle vertex corrections are neglected. EDMFT is consistent only in the strongly correlated regime and near half-filling, where the physics is predominantly local. △ Less

Submitted 5 August, 2015; v1 submitted 20 May, 2015; originally announced May 2015.

Comments: 11 pages (incl. appendix), 4 figures

Journal ref: Phys. Rev. B 92, 085106 (2015)

Beyond extended dynamical mean-field theory: Dual boson approach to the two-dimensional extended Hubbard model

Authors: Erik G. C. P. van Loon, Alexander I. Lichtenstein, Mikhail I. Katsnelson, Olivier Parcollet, Hartmut Hafermann

Abstract: The dual boson approach [Ann. Phys. 327, 1320 (2012)] provides a means to construct a diagrammatic expansion around the extended dynamical mean-field theory (EDMFT). In this paper, we present the numerical implementation of the approach and apply it to the extended Hubbard model with nearest-neighbor interaction $V$. We calculate the EDMFT phase diagram and study the effect of diagrams beyond EDMF… ▽ More The dual boson approach [Ann. Phys. 327, 1320 (2012)] provides a means to construct a diagrammatic expansion around the extended dynamical mean-field theory (EDMFT). In this paper, we present the numerical implementation of the approach and apply it to the extended Hubbard model with nearest-neighbor interaction $V$. We calculate the EDMFT phase diagram and study the effect of diagrams beyond EDMFT on the transition to the charge-ordered phase. Including diagrammatic corrections to the EDMFT polarization shifts the EDMFT phase boundary to lower values of $V$. The approach interpolates between the random phase approximation in the weak coupling limit and EDMFT for strong coupling. Neglecting vertex corrections leads to results reminiscent of the EDMFT+$GW$ approximation. We however find significant deviations from the dual boson results already for small values of the interaction, emphasizing the crucial importance of fermion-boson vertex corrections. △ Less

Submitted 29 December, 2014; v1 submitted 9 August, 2014; originally announced August 2014.

Comments: Published version; 24 pages, 32 figures

Journal ref: Phys. Rev. B 90, 235135 (2014)

Collective Charge Excitations of Strongly Correlated Electrons, Vertex Corrections and Gauge Invariance

Authors: Hartmut Hafermann, Erik G. C. P. van Loon, Mikhail I. Katsnelson, Alexander I. Lichtenstein, Olivier Parcollet

Abstract: We consider the collective, long-wavelength charge excitations in correlated media in presence of short- and long-range forces. As an example for the case of a short-range interaction, we examine the two-dimensional Hubbard model within dynamical mean-field theory (DMFT). It is shown explicitly that the DMFT susceptibility including vertex corrections respects the Ward identity and yields a manife… ▽ More We consider the collective, long-wavelength charge excitations in correlated media in presence of short- and long-range forces. As an example for the case of a short-range interaction, we examine the two-dimensional Hubbard model within dynamical mean-field theory (DMFT). It is shown explicitly that the DMFT susceptibility including vertex corrections respects the Ward identity and yields a manifestly gauge-invariant response in finite dimensions. For computing the susceptibility, we use a different expression and establish its formal equivalence to the standard DMFT formula. It allows for a more stable analytical continuation. We find a zero-sound mode expected for short-range forces. The relation between the vertex corrections, gauge invariance, and the appearance of the collective modes is discussed. Long-range forces are treated within extended dynamical mean-field theory. In order to obtain a gauge-invariant response, it is necessary to additionally incorporate some non local vertex corrections into the polarization. In doing so, we obtain plasmons in the three-dimensional Hubbard model. The plasma frequency is determined by the (single-particle) density distribution as a consequence of gauge invariance. We compare this result with the plasma frequency extracted from the analytical continuation of the susceptibility. It is in good agreement with the prediction from the gauge-invariance condition. △ Less

Submitted 7 December, 2014; v1 submitted 25 June, 2014; originally announced June 2014.

Comments: Published version; 21 pages, 22 figures

Journal ref: Phys. Rev. B 90, 235105 (2014)

Plasmons in strongly correlated systems: spectral weight transfer and renormalized dispersion

Authors: E. G. C. P. van Loon, H. Hafermann, A. I. Lichtenstein, A. N. Rubtsov, M. I. Katsnelson

Abstract: We study the charge-density dynamics within the two-dimensional extended Hubbard model in the presence of long-range Coulomb interaction across the metal-insulator transition point. To take into account strong correlations we start from self-consistent extended dynamical mean-field theory and include non-local dynamical vertex corrections through a ladder approximation to the polarization operator… ▽ More We study the charge-density dynamics within the two-dimensional extended Hubbard model in the presence of long-range Coulomb interaction across the metal-insulator transition point. To take into account strong correlations we start from self-consistent extended dynamical mean-field theory and include non-local dynamical vertex corrections through a ladder approximation to the polarization operator. This is necessary to fulfill charge conservation and to describe plasmons in the correlated state. The calculated plasmon spectra are qualitatively different from those in the random-phase approximation: they exhibit a spectral density transfer and a renormalized dispersion with enhanced deviation from the canonical $\sqrt{q}$-behavior. Both features are reminiscent of interaction induced changes found in single-electron spectra of strongly correlated systems. △ Less

Submitted 23 December, 2014; v1 submitted 24 June, 2014; originally announced June 2014.

Comments: 5 pages, 5 figures + appendix (3 pages, 1 figure)

Journal ref: Phys. Rev. Lett. 113, 246407 (2014)