The conventional, single band Hubbard model captures the dynamics that arise in strongly correlated electron systems in the competition between electron transport and local Coulomb interactions. However, some of the important physics of real lattices is overlooked by its simplicity. The dynamic Hubbard model (DHM) provides additional degrees of freedom for electron states in a lattice by including a higher energy orbital and hybridization parameters for electrons to take advantage of a reduced Coulomb repulsion in the upper orbital when the lattice is at half-filling. This ``relaxation'' of the electronic state is accompanied by an asymmetry between electron and hole quasiparticles in the model. In this thesis we employ the computational tool of multiorbital dynamical mean field theory (MODMFT) to study electron-hole asymmetry in the DHM as evidenced in the behaviour of the quasiparticle weight and optical conductivity of electrons. In agreement with previous exact diagonalization work on the DHM, we find that the relaxation of electronic states is associated with quasiparticle dressing in the hole regime of the lower orbital, and a transfer of optical conductivity spectral weight from low to high frequencies. These observations are evidences of the electron-hole asymmetry in the two-orbital system. In this context, we evaluate the role of two critical variables in the DHM: the local and nonlocal hybridizations which parameterize interorbital electron transitions in the Hamiltonian. We demonstrate that the electron-hole asymmetry of the model has a more complex relationship to these hybridization parameters than previously identified. Mott physics plays an essential role in the quasiparticle dressing due to the local hybridization, while the nonlocal hybridization promotes electron-hole asymmetry by providing a separate channel for electron transport which serves to reduce the quasiparticle weight of holes in the lower band. We also briefly examine the phenomenon of orbital selective Mott transitions in the DHM within the parameter regime of electron state relaxation, and evaluate the role played by the hybridization parameters in this context.The conventional, single band Hubbard model captures the dynamics that arise in strongly correlated electron systems in the competition between electron transport and local Coulomb interactions. However, some of the important physics of real lattices is overlooked by its simplicity. The dynamic Hubbard model (DHM) provides additional degrees of freedom for electron states in a lattice by including a higher energy orbital and hybridization parameters for electrons to take advantage of a reduced Coulomb repulsion in the upper orbital when the lattice is at half-filling. This ``relaxation'' of the electronic state is accompanied by an asymmetry between electron and hole quasiparticles in the model. In this thesis we employ the computational tool of multiorbital dynamical mean field theory (MODMFT) to study electron-hole asymmetry in the DHM as evidenced in the behaviour of the quasiparticle weight and optical conductivity of electrons. In agreement with previous exact diagonalization work on the DHM, we find that the relaxation of electronic states is associated with quasiparticle dressing in the hole regime of the lower orbital, and a transfer of optical conductivity spectral weight from low to high frequencies. These observations are evidences of the electron-hole asymmetry in the two-orbital system. In this context, we evaluate the role of two critical variables in the DHM: the local and nonlocal hybridizations which parameterize interorbital electron transitions in the Hamiltonian. We demonstrate that the electron-hole asymmetry of the model has a more complex relationship to these hybridization parameters than previously identified. Mott physics plays an essential role in the quasiparticle dressing due to the local hybridization, while the nonlocal hybridization promotes electron-hole asymmetry by providing a separate channel for electron transport which serves to reduce the quasiparticle weight of holes in the lower band. We also briefly examine the phenomenon of orbital selective Mott transitions in the DHM within the parameter regime of electron state relaxation, and evaluate the role played by the hybridization parameters in this context.

Jan 1, 20150 citations

Two-particle self-consistent approach for multiorbital models: Application to the Emery model

C. Gauvin-Ndiaye, J. Leblanc, S. Marin +3 more

The Emery model, or three-band Hubbard model, is a Hamiltonian that is thought to contain much of the physics of cuprate superconductors. This model includes two noninteracting $p$ orbitals and one interacting $d$ orbital per unit cell. Few methods that can solve multiorbital interacting Hamiltonians reliably and efficiently exist. Here, we introduce an application of the two-particle self-consistent (TPSC) approach to the Emery model. We construct this method within the framework of the TPSC $+$ DMFT method, which can be seen as a way to introduce nonlocal corrections to dynamical mean-field theory (DMFT). We show that interacting orbital densities, rather than the noninteracting ones, must be used in the calculations. For the Emery model, we find that at constant bare interaction $U$ , the vertex for spin fluctuations, ${U}_{sp}$ , decreases rapidly with filling. This may be one of the factors that contributes to electron-doped cuprates appearing less correlated than hole-doped ones. More generally, our work opens the road to the application of the TPSC approach to spin fluctuations in multiorbital models.The Emery model, or three-band Hubbard model, is a Hamiltonian that is thought to contain much of the physics of cuprate superconductors. This model includes two noninteracting $p$ orbitals and one interacting $d$ orbital per unit cell. Few methods that can solve multiorbital interacting Hamiltonians reliably and efficiently exist. Here, we introduce an application of the two-particle self-consistent (TPSC) approach to the Emery model. We construct this method within the framework of the TPSC $+$ DMFT method, which can be seen as a way to introduce nonlocal corrections to dynamical mean-field theory (DMFT). We show that interacting orbital densities, rather than the noninteracting ones, must be used in the calculations. For the Emery model, we find that at constant bare interaction $U$, the vertex for spin fluctuations, ${U}_{sp}$, decreases rapidly with filling. This may be one of the factors that contributes to electron-doped cuprates appearing less correlated than hole-doped ones. More generally, our work opens the road to the application of the TPSC approach to spin fluctuations in multiorbital models.

arXiv:2308.14091Physical review. B./Physical review. BApr 2, 20245 citations

Application of Multi-Orbital DMFT to the Dynamic Hubbard Model

Christopher Polachic, F. Marsiglio

Bulletin of the American Physical SocietyMar 19, 20130 citations

Precise ground state of multiorbital Mott systems via the variational discrete action theory

Zhengqian Cheng, Chris A. Marianetti

Determining the ground state of multiorbital Hubbard models is critical for understanding strongly correlated electron materials, yet, existing methods struggle to simultaneously reach zero temperature and infinite system size. The de facto standard is to approximate a finite dimension multiorbital Hubbard model with a $d=\ensuremath{\infty}$ version, which can then be formally solved via the dynamical mean-field theory (DMFT), although the DMFT solution is limited by the state of unbiased impurity solvers for zero temperature and multiple orbitals. The recently developed variational discrete action theory (VDAT) offers a new approach to solve the $d=\ensuremath{\infty}$ Hubbard model, with a variational ansatz that is controlled by an integer $\mathcal{N}$ , and monotonically approaches the exact solution at an increasing computational cost. Here we propose a decoupled minimization algorithm to implement VDAT for the multiorbital Hubbard model in $d=\ensuremath{\infty}$ and study $\mathcal{N}=2\ensuremath{-}4$. At $\mathcal{N}=2$ , VDAT rigorously recovers the multiorbital Gutzwiller approximation, reproducing known results. At $\mathcal{N}=3$ , VDAT precisely captures the competition between the Hubbard $U$ , Hund $J$ , and crystal field $\mathrm{\ensuremath{\Delta}}$ in the two-orbital Hubbard model over all parameter space, with a negligible computational cost. For sufficiently large $U/t$ and $J/U$ , we show that $\mathrm{\ensuremath{\Delta}}$ drives a first-order transition within the Mott insulating regime. In the large orbital polarization limit with finite $J/U$ , we find that interactions have a nontrivial effect even for small $U/t$ . VDAT will have far ranging implications for understanding multiorbital model Hamiltonians and strongly correlated electron materials.Determining the ground state of multiorbital Hubbard models is critical for understanding strongly correlated electron materials, yet, existing methods struggle to simultaneously reach zero temperature and infinite system size. The de facto standard is to approximate a finite dimension multiorbital Hubbard model with a $d=\ensuremath{\infty}$ version, which can then be formally solved via the dynamical mean-field theory (DMFT), although the DMFT solution is limited by the state of unbiased impurity solvers for zero temperature and multiple orbitals. The recently developed variational discrete action theory (VDAT) offers a new approach to solve the $d=\ensuremath{\infty}$ Hubbard model, with a variational ansatz that is controlled by an integer $\mathcal{N}$, and monotonically approaches the exact solution at an increasing computational cost. Here we propose a decoupled minimization algorithm to implement VDAT for the multiorbital Hubbard model in $d=\ensuremath{\infty}$ and study $\mathcal{N}=2\ensuremath{-}4$. At $\mathcal{N}=2$, VDAT rigorously recovers the multiorbital Gutzwiller approximation, reproducing known results. At $\mathcal{N}=3$, VDAT precisely captures the competition between the Hubbard $U$, Hund $J$, and crystal field $\mathrm{\ensuremath{\Delta}}$ in the two-orbital Hubbard model over all parameter space, with a negligible computational cost. For sufficiently large $U/t$ and $J/U$, we show that $\mathrm{\ensuremath{\Delta}}$ drives a first-order transition within the Mott insulating regime. In the large orbital polarization limit with finite $J/U$, we find that interactions have a nontrivial effect even for small $U/t$. VDAT will have far ranging implications for understanding multiorbital model Hamiltonians and strongly correlated electron materials.

arXiv:2206.08433Physical review. B./Physical review. BNov 18, 20227 citations

DMFT+NRG study of orbital differentiation in three-orbital Hubbard models

Fabian B. Kugler, Seung‐Sup B. Lee, Manuel Zingl +4 more

Bulletin of the American Physical SocietyMar 5, 20190 citations

The Emery model, or three-band Hubbard model, is a Hamiltonian that is thought to contain much of the physics of cuprate superconductors. This model includes two noninteracting

p

orbitals and one interacting

d

orbital per unit cell. Few methods that can solve multiorbital interacting Hamiltonians reliably and efficiently exist. Here, we introduce an application of the two-particle self-consistent (TPSC) approach to the Emery model. We construct this method within the framework of the TPSC

+

DMFT method, which can be seen as a way to introduce nonlocal corrections to dynamical mean-field theory (DMFT). We show that interacting orbital densities, rather than the noninteracting ones, must be used in the calculations. For the Emery model, we find that at constant bare interaction

U

, the vertex for spin fluctuations,

{U}_{sp}

, decreases rapidly with filling. This may be one of the factors that contributes to electron-doped cuprates appearing less correlated than hole-doped ones. More generally, our work opens the road to the application of the TPSC approach to spin fluctuations in multiorbital models.

\mathcal{N}

, and monotonically approaches the exact solution at an increasing computational cost. Here we propose a decoupled minimization algorithm to implement VDAT for the multiorbital Hubbard model in $d=\ensuremath{\infty}$ and study $\mathcal{N}=2\ensuremath{-}4$. At

\mathcal{N}=2

, VDAT rigorously recovers the multiorbital Gutzwiller approximation, reproducing known results. At

\mathcal{N}=3

, VDAT precisely captures the competition between the Hubbard

U

, Hund

J

, and crystal field $\mathrm{\ensuremath{\Delta}}$ in the two-orbital Hubbard model over all parameter space, with a negligible computational cost. For sufficiently large

U/t

and

J/U

, we show that $\mathrm{\ensuremath{\Delta}}$ drives a first-order transition within the Mott insulating regime. In the large orbital polarization limit with finite

J/U

, we find that interactions have a nontrivial effect even for small

U/t

. VDAT will have far ranging implications for understanding multiorbital model Hamiltonians and strongly correlated electron materials.