Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2507.10717 · Physical review. B./Physical review. B

Mixed-configuration approximation for multiorbital systems out of equilibrium

Tommaso Maria Mazzocchi, Daniel Werner, Markus Aichhorn, Enrico Arrigoni

Published Sep 30, 2025• 3 citations
arXivHTMLPDF

Abstract

We propose a mixed-configuration approximation based on single-band impurity solvers to efficiently study nonequilibrium multiorbital systems at moderate computational cost. In this work, we merge the approach with the so-called auxiliary master equation approach. As a benchmark, we first show that our approach reproduces the results of quantum Monte Carlo (QMC) for two-orbital impurity models at equilibrium with overall good accuracy, especially for nondegenerate orbitals. We then use our approach as an impurity solver for dynamical mean-field theory (DMFT) to address the case of a two-orbital, realistic layered structure, recovering the strong crystal-field-driven charge polarization observed by solving the DMFT self-consistent cycle with QMC, albeit slightly reduced. Finally, we address a prototype nonequilibrium setup by sandwiching this layer between metallic contacts subject to a bias voltage described by different chemical potentials. This simplified model demonstrates our method's potential to access nonequilibrium steady-state behavior of multiorbital, realistic materials. These findings provide a first-step basis for theoretical studies of nonequilibrium properties of multiorbital compounds directly in the real frequency domain.

DOI
10.1103/j3c8-cy15

Related papers

Transport through a correlated interface: Auxiliary master equation approach

Irakli Titvinidze, Antonius Dorda, Wolfgang von der Linden +1 more

We present improvements of a recently introduced numerical method [E. Arrigoni et al., Phys. Rev. Lett. 110, 086403 (2013)] to compute steady-state properties of strongly correlated electronic systems out of equilibrium. The method can be considered as a nonequilibrium generalization of exact diagonalization based dynamical mean-field theory (DMFT). The key modification for the nonequilibrium situation consists in addressing the DMFT impurity problem within an auxiliary system consisting of the correlated impurity, ${N}_{b}$ uncorrelated bath sites, and two Markovian environments (sink and reservoir). Algorithmic improvements in the impurity solver allow to treat efficiently larger values of ${N}_{b}$ than previously in DMFT. This increases the accuracy of the results and is crucial for a correct description of the physical behavior of the system in the relevant parameter range including a semiquantitative description of the Kondo regime. To illustrate the approach, we consider a monoatomic layer of correlated orbitals, described by the single-band Hubbard model, attached to two metallic leads. The nonequilibrium situation is driven by a bias voltage applied to the leads. For this system, we investigate the spectral function and the steady-state current-voltage characteristics in the weakly as well as in the strongly interacting limit. In particular, we investigate the nonequilibrium behavior of quasiparticle excitations within the Mott gap of the correlated layer. We find for low-bias voltage Kondo-type behavior in the vicinity of the insulating phase. In particular, we observe a splitting of the Kondo resonance as a function of the bias voltage.

arXiv:1508.02953•Physical Review B•Dec 21, 2015•24 citations

Nonequilibrium Dynamical Mean-Field Theory: An Auxiliary Quantum Master Equation Approach

Enrico Arrigoni, Michael Knap, Wolfgang von der Linden

We introduce a versatile method to compute electronic steady-state properties of strongly correlated extended quantum systems out of equilibrium. The approach is based on dynamical mean-field theory (DMFT), in which the original system is mapped onto an auxiliary nonequilibrium impurity problem imbedded in a Markovian environment. The steady-state Green's function of the auxiliary system is solved by full diagonalization of the corresponding Lindblad equation. The approach can be regarded as the nontrivial extension of the exact-diagonalization-based DMFT to the nonequilibrium case. As a first application, we consider an interacting Hubbard layer attached to two metallic leads and present results for the steady-state current and the nonequilibrium density of states.

arXiv:1210.4167•Physical Review Letters•Feb 20, 2013•141 citations

Configuration interaction based nonequilibrium steady state impurity solver

Daniel Werner, Jan Lotze, Enrico Arrigoni

We present a solver for correlated impurity problems out of equilibrium based on a combination of the so-called auxiliary master equation approach (AMEA) and the configuration interaction expansion. Within AMEA one maps the original impurity model onto an auxiliary open quantum system with a restricted number of bath sites which can be addressed by numerical many-body approaches such as Lanczos/Arnoldi exact diagonalization (ED) or matrix product states (MPS). While the mapping becomes exponentially more accurate with increasing number of bath sites, ED implementations are severely limited due to the fast increase of the Hilbert space dimension for open systems, and the MPS solver typically requires rather long runtimes. Here, we propose to adopt a configuration interaction approach augmented by active space extension to solve numerically the correlated auxiliary open quantum system. This allows access to a larger number of bath sites at lower computational costs than for plain ED. We benchmark the approach with numerical renormalization group results in equilibrium and with MPS out of equilibrium. In particular, we evaluate the current, the conductance, as well as the Kondo peak and its splitting as a function of increasing bias voltage below the Kondo temperature ${T}_{\text{K}}$. We obtain a rather accurate scaling of the conductance as a function of the bias voltage and temperature rescaled by ${T}_{\text{K}}$ for moderate to strong interactions in a wide range of parameters. The approach combines the fast runtime of ED with an accuracy close to the one achieved by MPS making it an attractive solver for nonequilibrium dynamical mean field theory.

arXiv:2210.09623•Physical review. B./Physical review. B•Feb 8, 2023•33 citations

Non-Equilibrium Dynamical Mean Field Theory with Quantum Monte Carlo

Qiaoyuan Dong

This work is motivated by the fact that the investigation of non-equilibrium phenomena in strongly correlated electron systems has developed into one of the most active and exciting branches of condensed matter physics as it provides rich new insights that could not be obtained from the study of equilibrium situations. However, a theoretical description of those phenomena is missing. Therefore, in this thesis, we develop a numerical method that can be used to study two minimal models -- the Hubbard model and the Anderson impurity model with general parameter range and time dependence. We begin by introducing the theoretical framework and the general features of the Hubbard model. We then describe the dynamical mean field theory (DMFT), which was first invented by Georges in 1992. It provides a feasible way to approach strongly correlated electron systems and reduces the complexity of the calculations via a mapping of lattice models onto quantum impurity models subject to a self-consistency condition. We employ the non-equilibrium extension of DMFT and map the Hubbard model to the single impurity Anderson model (SIAM). Since the fundamental component of the DMFT method is a solver of the single impurity Anderson model, we continue with a description of the formalism to study the real-time dynamics of the impurity model staring at its thermal equilibrium state. We utilize the non-equilibrium strong-coupling perturbation theory and derive semi-analytical approximation methods such as the non-crossing approximation (NCA) and the one-crossing approximation (OCA). We then use the Quantum Monte-Carlo method (QMC) as a numerically exact method and present proper measurements of local observables, current and Green's functions. We perform simulations of the current after a quantum quench from equilibrium by rapidly applying a bias voltage in a wide range of initial temperatures. The current exhibits short equilibrium times and saturates upon the decrease of temperature at all times, indicating Kondo behavior both in the transient regime and in the steady state. However, this bare QMC solver suffers from a dynamical sign problem for long time propagations. To overcome the limitations of this bare treatment, we introduce the ``Inchworm algorithm'', based on iteratively reusing the information obtained in previous steps to extend the propagation to longer times and stabilize the calculations. We show that this algorithm greatly reduces the required order for each simulation and re-scales the exponential challenge to quadratic in time. We introduce a method to compute Green's functions, spectral functions, and currents for inchworm Monte Carlo and show how systematic error assessments in real time can be obtained. We illustrate the capabilities of the algorithm with a study of the behavior of quantum impurities after an instantaneous voltage quench from a thermal equilibrium state. We conclude with the applications of the unbiased inchworm impurity solver to DMFT calculations. We employ the methods for a study of the one-band paramagnetic Hubbard model on the Bethe lattice in equilibrium, where the DMFT approximation becomes exact. We begin with a brief introduction of the Mott metal insulator phase diagram. We present the results of both real time Green's functions and spectral functions from our nonequilibrium calculations. We observe the metal-insulator crossover as the on-site interaction is increased and the formation of a quasi-particle peak as the temperature is lowered. We also illustrate the convergence of our algorithms in different aspects.

•Jan 1, 2017•0 citations

Quantum Monte Carlo impurity solvers for multi-orbital problems and frequency-dependent interactions

Hiroshi Shinaoka, Fakher F. Assaad, Nils Blümer +1 more

The solution of an auxiliary quantum impurity system is the computationally expensive step in dynamical mean field theory simulations of lattice models and materials. In this review, we discuss Monte Carlo based impurity solvers, which are suitable for a wide range of applications. In particular, we present an efficient implementation of the hybridization expansion approach, which enables the simulation of multiorbital impurity problems with off-diagonal and complex hybridizations, and dynamically screened (retarded) density-density interactions. As a complementary approach, we discuss an impurity solver based on the determinant Monte Carlo method, which scales favorably with inverse temperature and hence provides access to the very low temperature regime. The usefulness of these state-of-the-art impurity solvers is demonstrated with applications to the downfolding problem, i.e., the systematic derivation of dynamically screened interactions for low-energy effective models, and to pyrochlore iridates, where the spin-orbit coupling leads to complex hybridization functions in a multi-orbital system. As a benchmark for cluster extensions of dynamical mean field theory, we also present results from lattice Monte Carlo simulations for the momentum dependence of the pseudo-gap in the half-filled two-dimensional Hubbard model.

The European Physical Journal Special Topics•Jul 1, 2017•11 citations

Weak-coupling tensor cross interpolation impurity solver for nonequilibrium dynamical mean-field theory

Shuta Matsuura, Hiroshi Shinaoka, Philipp Werner +1 more

Simulating nonequilibrium quantum many-body systems remains a major challenge due to the exponential growth of the computational complexity with real time. Here we implement a nonequilibrium impurity solver based on the weak-coupling expansion and the tensor cross interpolation (TCI), and apply it to nonequilibrium dynamical mean-field theory (DMFT). The method approximates the integrands of the high-dimensional integrals arising in the weak-coupling expansion in a tensor-train form, enabling efficient evaluations without stochastic sampling and thereby mitigating the sign problem affecting continuous-time quantum Monte Carlo (CT-QMC) methods. Benchmark calculations for an exactly solvable nonequilibrium impurity model agree well with the exact results and reveal a low-rank structure of the integrands. When applied to interaction-quench problems in the half-filled Hubbard model, the method reproduces fast thermalization at a critical interaction strength with accuracy comparable to CT-QMC. Away from half filling, where the sign problem becomes even more severe, the present approach remains well controlled, revealing a crossover instead of a sharply defined fast thermalization point in the 3/4-filled case. The solver can also be applied to steady-state DMFT problems, yielding accurate spectral functions in the metallic regime without analytic continuation.

arXiv:2607.00702•Physical review. B./Physical review. B•Oct 14, 2026•0 citations
We propose a mixed-configuration approximation based on single-band impurity solvers to efficiently study nonequilibrium multiorbital systems at moderate computational cost. In this work, we merge the approach with the so-called auxiliary master equation approach. As a benchmark, we first show that our approach reproduces the results of quantum Monte Carlo (QMC) for two-orbital impurity models at equilibrium with overall good accuracy, especially for nondegenerate orbitals. We then use our approach as an impurity solver for dynamical mean-field theory (DMFT) to address the case of a two-orbital, realistic layered structure, recovering the strong crystal-field-driven charge polarization observed by solving the DMFT self-consistent cycle with QMC, albeit slightly reduced. Finally, we address a prototype nonequilibrium setup by sandwiching this layer between metallic contacts subject to a bias voltage described by different chemical potentials. This simplified model demonstrates our method's potential to access nonequilibrium steady-state behavior of multiorbital, realistic materials. These findings provide a first-step basis for theoretical studies of nonequilibrium properties of multiorbital compounds directly in the real frequency domain.
We present improvements of a recently introduced numerical method [E. Arrigoni et al., Phys. Rev. Lett. 110, 086403 (2013)] to compute steady-state properties of strongly correlated electronic systems out of equilibrium. The method can be considered as a nonequilibrium generalization of exact diagonalization based dynamical mean-field theory (DMFT). The key modification for the nonequilibrium situation consists in addressing the DMFT impurity problem within an auxiliary system consisting of the correlated impurity, Nb{N}_{b}Nb​ uncorrelated bath sites, and two Markovian environments (sink and reservoir). Algorithmic improvements in the impurity solver allow to treat efficiently larger values of Nb{N}_{b}Nb​ than previously in DMFT. This increases the accuracy of the results and is crucial for a correct description of the physical behavior of the system in the relevant parameter range including a semiquantitative description of the Kondo regime. To illustrate the approach, we consider a monoatomic layer of correlated orbitals, described by the single-band Hubbard model, attached to two metallic leads. The nonequilibrium situation is driven by a bias voltage applied to the leads. For this system, we investigate the spectral function and the steady-state current-voltage characteristics in the weakly as well as in the strongly interacting limit. In particular, we investigate the nonequilibrium behavior of quasiparticle excitations within the Mott gap of the correlated layer. We find for low-bias voltage Kondo-type behavior in the vicinity of the insulating phase. In particular, we observe a splitting of the Kondo resonance as a function of the bias voltage.
We introduce a versatile method to compute electronic steady-state properties of strongly correlated extended quantum systems out of equilibrium. The approach is based on dynamical mean-field theory (DMFT), in which the original system is mapped onto an auxiliary nonequilibrium impurity problem imbedded in a Markovian environment. The steady-state Green's function of the auxiliary system is solved by full diagonalization of the corresponding Lindblad equation. The approach can be regarded as the nontrivial extension of the exact-diagonalization-based DMFT to the nonequilibrium case. As a first application, we consider an interacting Hubbard layer attached to two metallic leads and present results for the steady-state current and the nonequilibrium density of states.
We present a solver for correlated impurity problems out of equilibrium based on a combination of the so-called auxiliary master equation approach (AMEA) and the configuration interaction expansion. Within AMEA one maps the original impurity model onto an auxiliary open quantum system with a restricted number of bath sites which can be addressed by numerical many-body approaches such as Lanczos/Arnoldi exact diagonalization (ED) or matrix product states (MPS). While the mapping becomes exponentially more accurate with increasing number of bath sites, ED implementations are severely limited due to the fast increase of the Hilbert space dimension for open systems, and the MPS solver typically requires rather long runtimes. Here, we propose to adopt a configuration interaction approach augmented by active space extension to solve numerically the correlated auxiliary open quantum system. This allows access to a larger number of bath sites at lower computational costs than for plain ED. We benchmark the approach with numerical renormalization group results in equilibrium and with MPS out of equilibrium. In particular, we evaluate the current, the conductance, as well as the Kondo peak and its splitting as a function of increasing bias voltage below the Kondo temperature TK{T}_{\text{K}}TK​. We obtain a rather accurate scaling of the conductance as a function of the bias voltage and temperature rescaled by TK{T}_{\text{K}}TK​ for moderate to strong interactions in a wide range of parameters. The approach combines the fast runtime of ED with an accuracy close to the one achieved by MPS making it an attractive solver for nonequilibrium dynamical mean field theory.
This work is motivated by the fact that the investigation of non-equilibrium phenomena in strongly correlated electron systems has developed into one of the most active and exciting branches of condensed matter physics as it provides rich new insights that could not be obtained from the study of equilibrium situations. However, a theoretical description of those phenomena is missing. Therefore, in this thesis, we develop a numerical method that can be used to study two minimal models -- the Hubbard model and the Anderson impurity model with general parameter range and time dependence. We begin by introducing the theoretical framework and the general features of the Hubbard model. We then describe the dynamical mean field theory (DMFT), which was first invented by Georges in 1992. It provides a feasible way to approach strongly correlated electron systems and reduces the complexity of the calculations via a mapping of lattice models onto quantum impurity models subject to a self-consistency condition. We employ the non-equilibrium extension of DMFT and map the Hubbard model to the single impurity Anderson model (SIAM). Since the fundamental component of the DMFT method is a solver of the single impurity Anderson model, we continue with a description of the formalism to study the real-time dynamics of the impurity model staring at its thermal equilibrium state. We utilize the non-equilibrium strong-coupling perturbation theory and derive semi-analytical approximation methods such as the non-crossing approximation (NCA) and the one-crossing approximation (OCA). We then use the Quantum Monte-Carlo method (QMC) as a numerically exact method and present proper measurements of local observables, current and Green's functions. We perform simulations of the current after a quantum quench from equilibrium by rapidly applying a bias voltage in a wide range of initial temperatures. The current exhibits short equilibrium times and saturates upon the decrease of temperature at all times, indicating Kondo behavior both in the transient regime and in the steady state. However, this bare QMC solver suffers from a dynamical sign problem for long time propagations. To overcome the limitations of this bare treatment, we introduce the ``Inchworm algorithm'', based on iteratively reusing the information obtained in previous steps to extend the propagation to longer times and stabilize the calculations. We show that this algorithm greatly reduces the required order for each simulation and re-scales the exponential challenge to quadratic in time. We introduce a method to compute Green's functions, spectral functions, and currents for inchworm Monte Carlo and show how systematic error assessments in real time can be obtained. We illustrate the capabilities of the algorithm with a study of the behavior of quantum impurities after an instantaneous voltage quench from a thermal equilibrium state. We conclude with the applications of the unbiased inchworm impurity solver to DMFT calculations. We employ the methods for a study of the one-band paramagnetic Hubbard model on the Bethe lattice in equilibrium, where the DMFT approximation becomes exact. We begin with a brief introduction of the Mott metal insulator phase diagram. We present the results of both real time Green's functions and spectral functions from our nonequilibrium calculations. We observe the metal-insulator crossover as the on-site interaction is increased and the formation of a quasi-particle peak as the temperature is lowered. We also illustrate the convergence of our algorithms in different aspects.
The solution of an auxiliary quantum impurity system is the computationally expensive step in dynamical mean field theory simulations of lattice models and materials. In this review, we discuss Monte Carlo based impurity solvers, which are suitable for a wide range of applications. In particular, we present an efficient implementation of the hybridization expansion approach, which enables the simulation of multiorbital impurity problems with off-diagonal and complex hybridizations, and dynamically screened (retarded) density-density interactions. As a complementary approach, we discuss an impurity solver based on the determinant Monte Carlo method, which scales favorably with inverse temperature and hence provides access to the very low temperature regime. The usefulness of these state-of-the-art impurity solvers is demonstrated with applications to the downfolding problem, i.e., the systematic derivation of dynamically screened interactions for low-energy effective models, and to pyrochlore iridates, where the spin-orbit coupling leads to complex hybridization functions in a multi-orbital system. As a benchmark for cluster extensions of dynamical mean field theory, we also present results from lattice Monte Carlo simulations for the momentum dependence of the pseudo-gap in the half-filled two-dimensional Hubbard model.
Simulating nonequilibrium quantum many-body systems remains a major challenge due to the exponential growth of the computational complexity with real time. Here we implement a nonequilibrium impurity solver based on the weak-coupling expansion and the tensor cross interpolation (TCI), and apply it to nonequilibrium dynamical mean-field theory (DMFT). The method approximates the integrands of the high-dimensional integrals arising in the weak-coupling expansion in a tensor-train form, enabling efficient evaluations without stochastic sampling and thereby mitigating the sign problem affecting continuous-time quantum Monte Carlo (CT-QMC) methods. Benchmark calculations for an exactly solvable nonequilibrium impurity model agree well with the exact results and reveal a low-rank structure of the integrands. When applied to interaction-quench problems in the half-filled Hubbard model, the method reproduces fast thermalization at a critical interaction strength with accuracy comparable to CT-QMC. Away from half filling, where the sign problem becomes even more severe, the present approach remains well controlled, revealing a crossover instead of a sharply defined fast thermalization point in the 3/4-filled case. The solver can also be applied to steady-state DMFT problems, yielding accurate spectral functions in the metallic regime without analytic continuation.