The description of interacting quantum impurity models in steady-state nonequilibrium is an open challenge for computational many-particle methods: the numerical requirement of using a finite number of lead levels and the physical requirement of describing a truly open quantum system are seemingly incompatible. One possibility to bridge this gap is the use of Lindblad-driven discretized leads (LDDL): one couples auxiliary continuous reservoirs to the discretized lead levels and represents these additional reservoirs by Lindblad terms in the Liouville equation. For quadratic models governed by Lindbladian dynamics, we present an elementary approach for obtaining correlation functions analytically. In a second part, we use this approach to explicitly discuss the conditions under which the continuum limit of the LDDL approach recovers the correct representation of thermal reservoirs. As an analytically solvable example, the nonequilibrium resonant level model is studied in greater detail. Lastly, we present ideas towards a numerical evaluation of the suggested Lindblad equation for interacting impurities based on matrix product states. In particular, we present a reformulation of the Lindblad equation, which has the useful property that the leads can be mapped onto a chain where both the Hamiltonian dynamics and the Lindblad driving are local at the same time. Moreover, we discuss the possibility to combine the Lindblad approach with a logarithmic discretization needed for the exploration of exponentially small energy scales.
The description of interacting quantum impurity models in steady-state nonequilibrium is an open challenge for computational many-particle methods: the numerical requirement of using a finite number of lead levels and the physical requirement of describing a truly open quantum system are seemingly incompatible. One possibility to bridge this gap is the use of Lindblad-driven discretized leads (LDDL): one couples auxiliary continuous reservoirs to the discretized lead levels and represents these additional reservoirs by Lindblad terms in the Liouville equation. For quadratic models governed by Lindbladian dynamics, we present an elementary approach for obtaining correlation functions analytically. In a second part, we use this approach to explicitly discuss the conditions under which the continuum limit of the LDDL approach recovers the correct representation of thermal reservoirs. As an analytically solvable example, the nonequilibrium resonant level model is studied in greater detail. Lastly, we present ideas towards a numerical evaluation of the suggested Lindblad equation for interacting impurities based on matrix product states. In particular, we present a reformulation of the Lindblad equation, which has the useful property that the leads can be mapped onto a chain where both the Hamiltonian dynamics and the Lindblad driving are local at the same time. Moreover, we discuss the possibility to combine the Lindblad approach with a logarithmic discretization needed for the exploration of exponentially small energy scales.
We present a general scheme to address correlated nonequilibrium quantum impurity problems based on a mapping onto an auxiliary open quantum system of small size. The infinite fermionic reservoirs of the original system are thereby replaced by a small number N B of noninteracting auxiliary bath sites whose dynamics are described by a Lindblad equation, which can then be exactly solved by numerical methods such as Lanczos or matrix-product states. The mapping becomes exponentially exact with increasing N B , and is already quite accurate for small N B . Due to the presence of the intermediate bath sites, the overall dynamics acting on the impurity site is non-Markovian. While in previous work we put the focus on the manybody solution of the associated Lindblad problem, here we discuss the mapping scheme itself, which is an essential part of the overall approach. On the one hand, we provide technical details together with an in-depth discussion of the employed algorithms, and on the other hand, we present a detailed convergence study. The latter clearly demonstrates the above-mentioned exponential convergence of the procedure with increasing N B . Furthermore, the influence of temperature and an external bias voltage on the reservoirs is investigated. The knowledge of the particular convergence behavior is of great value to assess the applicability of the scheme to certain physical situations. Moreover, we study different geometries for the auxiliary system. On the one hand, this is of importance for advanced manybody solution techniques such as matrix product states which work well for short-ranged couplings, and on the other hand, it allows us to gain more insights into the underlying mechanisms when mapping non-Markovian reservoirs onto Lindblad-type impurity problems. Finally, we present results for the spectral function of the Anderson impurity model in and out of equilibrium and discuss the accuracy obtained with the different geometries of the auxiliary system. In particular, we show that allowing for complex Lindblad couplings produces a drastic improvement in the description of the Kondo resonance.
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.
The continuous coupling function in quantum impurity problems is exactly partitioned into a part represented by a finite-size Wilson chain and a part represented by a set of additional reservoirs, each coupled to one Wilson chain site. These additional reservoirs represent high-energy modes of the environment neglected by the numerical renormalization group and are required to restore the continuum limit of the original problem. We present a hybrid time-dependent numerical renormalization group approach which combines an accurate numerical renormalization group treatment of the nonequilibrium dynamics on the finite-size Wilson chain with a Bloch-Redfield formalism to include the effect of these additional reservoirs. Our approach overcomes the intrinsic shortcoming of the time-dependent numerical renormalization group approach induced by the bath discretization with a Wilson parameter $\mathrm{\ensuremath{\Lambda}}>1$. We analytically prove that for a system with a single chemical potential, the thermal equilibrium reduced density operator is the steady-state solution of the Bloch-Redfield master equation. For the numerical solution of this master equation, a Lanczos method is employed which couples all energy shells of the numerical renormalization group. The presented hybrid approach is applied to the real-time dynamics in correlated fermionic QISs. An analytical solution of the resonant-level model serves as a benchmark for the accuracy of the method which is then applied to nontrivial models, such as the interacting resonant-level model and the single-impurity Anderson model.
Theoretical studies of electron and energy transport in nanostructures usually involve nonequilibrium quantum impurity models, which feature a locally interacting quantum system (impurity) coupled to noninteracting bosonic and/or electronic baths each at its own equilibrium. In these models, the interplay of strong correlation, nonequilibrium effects and dissipation leads to rich and complex phenomena. As a nonequilibrium strongly correlated many-body problem, these models, except at some rare cases, have no exact analytical solution, and hence their treatment heavily relies upon numerical techniques. The goal of this thesis is to develop numerically accurate and relatively low-cost techniques for treatment of strongly correlated open nonequilibirum quantum systems. Such impurity solvers are in high demand both as standalone methods for simulation of response in nanoscale junctions and as part of divide-and-conquer schemes such as dynamical mean field theory (DMFT).For nonequilibirum open quantum systems, two characteristic energy scales governing the physics are intra-system interactions and strength of the system-baths couplings. Thus, systems strongly coupled to their baths with weak intra-system interactions or weakly coupled systems with strong local interactions can be treated within developed perturbation theory approaches. My research focuses on development of nonequilibirum dual techniques, which are capable to accurately treat situations where perturbative techniques are inapplicable due to absence of a small parameter in the system and which reduce to standard perturbative expansions in the limits of weak coupling or weak intra-system interaction. \tDual techniques rely on ability to solve exactly a simplified reference system. The latter includes exactly all intra-system interactions and accounts approximately for system-bath couplings. Dual techniques (described below) allow to account for deviation from the correct hybridization function in a organized and controlled manner of superperturbative expansion. Here, I use Markovian Lindbald quantum master equation (QME) as a solver for the reference system. In the thesis I discuss possibility of mapping between true non-Markov dynamics and simpler Markov Lindblad QME (mapping between physical and reference systems), numerical methods for the mapping and for solution of the resulting Lindblad QME, and application of the mapping in novel dual formulations: auxiliary-QME dual fermion and dual boson approaches. \tFirst I analytically prove that the continuous fermionic baths that the impurity is coupled to can be mapped onto a set of dissipative auxiliary modes under Lindblad QME. Based on this mapping, complicated non-Markovian dynamics of the open fermionic system can be treated by a much simpler Markovian Lindblad QME. However, nice mapping generally requires many auxiliary modes, and numerically only a small number of modes can be exactly diagonalized due to exponential increase of the Hilbert space. To perturbatively treat the difference between the true bath and the small-size auxiliary one, we introduce dual-fermion method. It can directly target the steady state and avoid long-time propagation. Besides current-voltage relation, it can also provide spectral functions, indicating itself as a potential impurity solver for DMFT. Furthermore, we generalize dual-fermion method to consider additionally the effects of bosonic environment; hence it is called dual-boson method. It can be used for studying electron-phonon interaction and light-matter interaction in nanojunctions. By comparing with numerically exact results we argue that both methods are very accurate and relatively cheap. \tFinally, matrix product state representation of the Lindblad QME enables us to treat a large number of auxiliary modes, which provide an exponentially improved mapping. Adding counting field to the auxiliary QME, we can compute the generating function of time-dependent full-counting statistics of electron transport, which yields not only electronic current but also higher-order cumulants (e.g., zero-frequency noise).
State-of-the-art approaches to extract transport coefficients of many-body quantum systems broadly fall into two categories: (i) they target the linear-response regime in terms of equilibrium correlation functions of the closed system; or (ii) they consider an open-system situation typically modeled by a Lindblad equation, where a nonequilibrium steady state emerges from driving the system at its boundaries. While quantitative agreement between (i) and (ii) has been found for selected model and parameter choices, also disagreement has been pointed out in the literature. Studying magnetization transport in the spin-1/2 XXZ chain, we here demonstrate that at weak driving, the nonequilibrium steady state in an open system, including its buildup in time, can remarkably be constructed just on the basis of correlation functions in the closed system. We numerically illustrate this direct correspondence of closed-system and open-system dynamics, and show that it allows the treatment of comparatively large open systems, usually only accessible to matrix product state simulations. We also point out potential pitfalls when extracting transport coefficients from nonequilibrium steady states in finite systems.
In the thesis we present an analytic approach towards exact description for steady state density operators of nonequilibrium quantum dynamics in the framework of open systems. We employ the so-called quantum Markovian semi-group evolution, i.e. a general form of time-autonomous positivity and trace-preserving dynamical equation for reduced density operators, by only allowing unitarity-breaking dissipative terms acting at the boundaries of a system. Such setup enables to simulate macroscopic reservoirs for different values of effective thermodynamic potentials, causing incoherent transitions between quantum states which are modeled with aid of the Lindblad operators. This serves as a simple minimalistic model for studying quantum transport properties, either in the linear response domain or in more general regimes far from canonical equilibrium. We are mainly exploring possibilities of identifying nonequilibrium situations which are amenable to exact description within matrix product state representation, by exclusively focusing on steady states, i.e. fixed points of the Lindblad equation, of certain prototypic interacting integrable spin chains driven by incoherent polarizing processes. Finally, we define a concept of pseudo-local extensive almost-conserved quantities by allowing a violation of time-invariance up to boundary-localized terms. We elucidate the role of such quantities on non-ergodic behaviour of temporal correlation functions rendering anomalous transport properties. It turns out that such conservation laws can be generated by means of boundary universal quantum transfer operators of the fundamental integrable models.