Voltage Quench Dynamics of a Kondo System

Quantum Monte Carlo Method in the Steady State

We present a numerically exact steady-state inchworm Monte Carlo method for nonequilibrium quantum impurity models. Rather than propagating an initial state to long times, the method is directly formulated in the steady state. This eliminates any need to traverse the transient dynamics and grants access to a much larger range of parameter regimes at vastly reduced computational costs. We benchmark the method on equilibrium Green's functions of quantum dots in the noninteracting limit and in the unitary limit of the Kondo regime. We then consider correlated materials described with dynamical mean field theory and driven away from equilibrium by a bias voltage. We show that the response of a correlated material to a bias voltage differs qualitatively from the splitting of the Kondo resonance observed in bias-driven quantum dots.

Generalized master equation for charge transport in a molecular junction: Exact memory kernels and their high order expansion

We derive a set of generalized master equations (GMEs) to study charge transport dynamics in molecular junctions using the Nakajima-Zwanzig-Mori projection operator approach. In the new GME, time derivatives of population on each quantum state of the molecule, as well as the tunneling current, are calculated as the convolution of time non-local memory kernels with populations on all system states. The non-Markovian memory kernels are obtained by combining the hierarchical equations of motion (HEOM) method and a previous derived Dyson relation for the exact kernel. A perturbative expansion of these memory kernels is then calculated using the extended HEOM developed in our previous work [M. Xu et al., J. Chem. Phys. 146, 064102 (2017)]. By using the resonant level model and the Anderson impurity model, we study properties of the exact memory kernels and analyze convergence properties of their perturbative expansions with respect to the system-bath coupling strength and the electron-electron repulsive energy. It is found that exact memory kernels calculated from HEOM exhibit short memory times and decay faster than the population and current dynamics. The high order perturbation expansion of the memory kernels can give converged results in certain parameter regimes. The Padé and Landau-Zener resummation schemes are also found to give improved results over low order perturbation theory.

Hierarchical equations of motion method based on Fano spectrum decomposition for low temperature environments

The hierarchical equations of motion (HEOM) method has become one of the most popular methods for the studies of the open quantum system. However, its applicability to systems at ultra-low temperatures is largely restrained by the enormous computational cost, which is caused by the numerous exponential functions required to accurately characterize the non-Markovian memory of the reservoir environment. To overcome this problem, a Fano spectrum decomposition (FSD) scheme has been proposed recently [Cui et al., J. Chem. Phys. 151, 024110 (2019)], which expands the reservoir correlation functions using polynomial-exponential functions and hence greatly reduces the size of the memory basis set. In this work, we explicitly establish the FSD-based HEOM formalisms for both bosonic and fermionic environments. The accuracy and efficiency of the FSD-based HEOM are exemplified by the calculated low-temperature dissipative dynamics of a spin-boson model and the dynamic and static properties of a single-orbital Anderson impurity model in the Kondo regime. The encouraging numerical results highlight the practicality and usefulness of the FSD-based HEOM method for general open systems at ultra-low temperatures.

Nonequilibrium Steady-State Transport in Quantum Impurity Models: A Thermofield and Quantum Quench Approach Using Matrix Product States

The numerical renormalization group (NRG) is tailored to describe interacting impurity models in equilibrium, but it faces limitations for steady-state nonequilibrium, arising, e.g., due to an applied bias voltage. We show that these limitations can be overcome by describing the thermal leads using a thermofield approach, integrating out high energy modes using NRG, and then treating the nonequilibrium dynamics at low energies using a quench protocol, implemented using the time-dependent density matrix renormalization group. This yields quantitatively reliable results for the current (with errors ≲3%) down to the exponentially small energy scales characteristic of impurity models. We present results of benchmark quality for the temperature and magnetic field dependence of the zero-bias conductance peak for the single-impurity Anderson model.

Time Evolution of the Kondo Resonance in Response to a Quench

We investigate the time evolution of the Kondo resonance in response to a quench by applying the time-dependent numerical renormalization group (TDNRG) approach to the Anderson impurity model in the strong correlation limit. For this purpose, we derive within the TDNRG approach a numerically tractable expression for the retarded two-time nonequilibrium Green function G(t+t^{'},t), and its associated time-dependent spectral function, A(ω,t), for times t both before and after the quench. Quenches from both mixed valence and Kondo correlated initial states to Kondo correlated final states are considered. For both cases, we find that the Kondo resonance in the zero temperature spectral function, a preformed version of which is evident at very short times t→0^{+}, only fully develops at very long times t≳1/T_{K}, where T_{K} is the Kondo temperature of the final state. In contrast, the final state satellite peaks develop on a fast time scale 1/Γ during the time interval -1/Γ≲t≲+1/Γ, where Γ is the hybridization strength. Initial and final state spectral functions are recovered in the limits t→-∞ and t→+∞, respectively. Our formulation of two-time nonequilibrium Green functions within the TDNRG approach provides a first step towards using this method as an impurity solver within nonequilibrium dynamical mean field theory.

Slowdown of the Electronic Relaxation Close to the Mott Transition

We investigate the time-dependent reformation of the quasiparticle peak in a correlated metal near the Mott transition, after the system is quenched into a hot electron state and equilibrates with an environment which is colder than the Fermi-liquid crossover temperature. Close to the transition, we identify a purely electronic bottleneck time scale, which depends on the spectral weight around the Fermi energy in the bad metallic phase in a nonlinear way. This time scale can be orders of magnitude larger than the bare and renormalized electronic hopping time, so that a separation of electronic and lattice time scales may break down. The results are obtained using nonequilibrium dynamical mean-field theory and a slave-rotor representation of the Anderson impurity model.