We present further developments of the auxiliary master equation approach (AMEA), a numerical method to simulate many-body quantum systems in as well as out of equilibrium and apply it to the interacting resonant level model to benchmark the new developments. In particular, our results are obtained by employing the stochastic wave functions method to solve the auxiliary open quantum system arising within AMEA. This development allows us to reach extremely low wall times for the calculation of correlation functions with respect to previous implementations of AMEA. An additional significant improvement is obtained by extrapolating a series of results obtained by increasing the number of auxiliary bath sites, N_{B}, used within the auxiliary open quantum system formally to the limit of N_{B}→∞. Results for the current-voltage characteristics and for equilibrium correlation functions are compared with the one obtained by exact and matrix-product states-based approaches. Further, we complement this benchmark by the presentation of spectral functions for higher temperatures where we find different behaviors around zero frequency depending on the hybridization strength.
We present further developments of the auxiliary master equation approach (AMEA), a numerical method to simulate many-body quantum systems in as well as out of equilibrium and apply it to the interacting resonant level model to benchmark the new developments. In particular, our results are obtained by employing the stochastic wave functions method to solve the auxiliary open quantum system arising within AMEA. This development allows us to reach extremely low wall times for the calculation of correlation functions with respect to previous implementations of AMEA. An additional significant improvement is obtained by extrapolating a series of results obtained by increasing the number of auxiliary bath sites, N_{B}, used within the auxiliary open quantum system formally to the limit of N_{B}→∞. Results for the current-voltage characteristics and for equilibrium correlation functions are compared with the one obtained by exact and matrix-product states-based approaches. Further, we complement this benchmark by the presentation of spectral functions for higher temperatures where we find different behaviors around zero frequency depending on the hybridization strength.
We present a numerical method for the study of correlated quantum impurity problems out of equilibrium, which is particularly suited to address steady-state properties within dynamical mean field theory. The approach, recently introduced by Arrigoni et al. [Phys. Rev. Lett. 110, 086403 (2013)], is based upon a mapping of the original impurity problem onto an auxiliary open quantum system, consisting of the interacting impurity coupled to bath sites as well as to a Markovian environment. The dynamics of the auxiliary system is governed by a Lindblad master equation whose parameters are used to optimize the mapping. The accuracy of the results can be readily estimated and systematically improved by increasing the number of auxiliary bath sites, or by introducing a linear correction. Here, we focus on a detailed discussion of the proposed approach including technical remarks. To solve for the Green's functions of the auxiliary impurity problem, a non-Hermitian Lanczos diagonalization is applied. As a benchmark, results for the steady-state current-voltage characteristics of the single-impurity Anderson model are presented. Furthermore, the bias dependence of the single-particle spectral function and the splitting of the Kondo resonance are discussed. In its present form, the method is fast, efficient, and features a controlled accuracy.
Within the recently introduced auxiliary master equation approach it is possible to address steady state properties of strongly correlated impurity models, small molecules or clusters efficiently and with high accuracy. It is particularly suited for dynamical mean field theory in the nonequilibrium as well as in the equilibrium case. The method is based on the solution of an auxiliary open quantum system, which can be made quickly equivalent to the original impurity problem. In its first implementation a Krylov space method was employed. Here, we aim at extending the capabilities of the approach by adopting matrix product states for the solution of the corresponding auxiliary quantum master equation. This allows for a drastic increase in accuracy and permits us to access the Kondo regime for large values of the interaction. In particular, we investigate the nonequilibrium steady state of a single impurity Anderson model and focus on the spectral properties for temperatures T below the Kondo temperature T_K and for small bias voltages phi. For the two cases considered, with T=T_K/4 and T=T_K/10 we find a clear splitting of the Kondo resonance into a two-peak structure for phi close above T_K. In the equilibrium case (phi=0) and for T=T_K/4, the obtained spectral function essentially coincides with the one from numerical renormalization group.
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 uncorrelated bath sites, and two Markovian environments (sink and reservoir). Algorithmic improvements in the impurity solver allow to treat efficiently larger values of 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.
This paper provides a detailed account of the numerical implementation of the stochastic equation of motion (SEOM) method for the dissipative dynamics of fermionic open quantum systems. To enable direct stochastic calculations, a minimal auxiliary space (MAS) mapping scheme is adopted, with which the time-dependent Grassmann fields are represented by c-number noises and a set of pseudo-operators. We elaborate on the construction of the system operators and pseudo-operators involved in the MAS-SEOM, along with the analytic expression for the particle current. The MAS-SEOM is applied to study the relaxation and voltage-driven dynamics of quantum impurity systems described by the single-level Anderson impurity model, and the numerical results are benchmarked against those of the highly accurate hierarchical equations of motion method. The advantages and limitations of the present MAS-SEOM approach are discussed extensively.
This thesis describes the development of practical and efficient computational approaches to the quantum dynamics of complex systems. Most of the work presented here relies on the generalized quantum master equation (GQME) formalism, which provides a simple equation of motion of reduced dimensionality for a set of dynamical quantities, e.g., nonequilibrium averages and equilibrium time correlation functions. The reduced dimensionality of the GQME comes at a cost: the introduction of the memory kernel, which accounts for the influence of all ``excluded'' degrees of freedom. Focusing first on the second-order perturbative treatment of the memory kernel known as Redfield theory, I present a collaborative effort to extend its applicability into highly non-Markovian regions via a mode freezing approach. In this method, a portion of bath modes characterized by low frequencies are treated as sources of static disorder and used to calculate modified Redfield dynamics. Application of the method to the spin-boson and FMO complex models indicates that the Redfield+frozen modes scheme consistently produces dynamics that are as good or better than bare Redfield dynamics. Next, we explore GQME approach coupled to the self-consistent solution of the memory kernel, which requires the calculation of auxiliary kernels. Previous implementations of the method had shown impressive boosts in efficiency and, when approximate methods were used to calculate the auxiliary kernels, accuracy over direct calculation of nonequilibrium averages. We show that this method, when formulated from the Mori perspective, is equally applicable to nonequilibrium averages and equilibrium correlation functions. In addition, we examine the dependence of the improvements afforded by the GQME framework on the choice projection operator and kernel closure. In particular, we demonstrate that improvements in efficiency, which rely on short memory lifetimes, are sensitively dependent on the choice of projection operator, and that the choice of kernel closure directly dictates the improvements in accuracy. In addition, we present evidence that indicates that the success of the GQME formalism when the auxiliary kernels are calculated via semi- and quasi-classical methods is largely due to the exact sampling of bath operators at t=0 required by the calculation of specific kernel closures. Next, we provide analytical arguments that delineate when the GQME framework coupled to the self-consistent solution of the memory kernel is likely to provide improvements in efficiency and accuracy. Finally, we present a path integral framework that can efficiently render the partially Wigner-transformed canonical density operator for systems coupled linearly to harmonic baths. This approach permits the direct calculation of any thermodynamic quantity and can be integrated into dynamical schemes like the Ehrenfest, surface hopping, or linearized semi-classical initial value representation methods to calculate equilibrium correlation functions.
We extend the twin-space formulation of the hierarchical equations of motion approach in combination with the matrix product state representation [R. Borrelli, J. Chem. Phys. 150, 234102 (2019)] to nonequilibrium scenarios where the open quantum system is coupled to a hybrid fermionic and bosonic environment. The key ideas used in the extension are a reformulation of the hierarchical equations of motion for the auxiliary density matrices into a time-dependent Schrödinger-like equation for an augmented multi-dimensional wave function as well as a tensor decomposition into a product of low-rank matrices. The new approach facilitates accurate simulations of non-equilibrium quantum dynamics in larger and more complex open quantum systems. The performance of the method is demonstrated for a model of a molecular junction exhibiting current-induced mode-selective vibrational excitation.