Green's functions from real-time bold-line Monte Carlo calculations: spectral properties of the nonequilibrium Anderson impurity model

Nonequilibrium steady state full counting statistics in the noncrossing approximation

Quantum transport is often characterized not just by mean observables like the particle or energy current but by their fluctuations and higher moments, which can act as detailed probes of the physical mechanisms at play. However, relatively few theoretical methods are able to access the full counting statistics (FCS) of transport processes through electronic junctions in strongly correlated regimes. While most experiments are concerned with steady state properties, most accurate theoretical methods rely on computationally expensive propagation from a tractable initial state. Here, we propose a simple approach for computing the FCS through a junction directly at the steady state, utilizing the propagator noncrossing approximation. Compared to time propagation, our method offers reduced computational cost at the same level of approximation, but the idea can also be used within other approximations or as a basis for numerically exact techniques. We demonstrate the method's capabilities by investigating the impact of lead dimensionality on electronic transport in the nonequilibrium Anderson impurity model at the onset of Kondo physics. Our results reveal a distinct signature of one dimensional leads in the noise and Fano factor not present for other dimensionalities, showing the potential of FCS measurements as a probe of the environment surrounding a quantum dot.

Steady-state properties of multi-orbital systems using quantum Monte Carlo

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

Shaping Electronic Flows with Strongly Correlated Physics

Nonequilibrium quantum transport is of central importance in nanotechnology. Its description requires the understanding of strong electronic correlations that couple atomic-scale phenomena to the nanoscale. So far, research in correlated transport has focused predominantly on few-channel transport, precluding the investigation of cross-scale effects. Recent theoretical advances enable the solution of models that capture the interplay between quantum correlations and confinement beyond a few channels. This problem is the focus of this study. We consider an atomic impurity embedded in a metallic nanosheet spanning two leads, showing that transport is significantly altered by tuning only the phase of a single local hopping parameter. Furthermore─depending on this phase─correlations reshape the electronic flow throughout the sheet, either funneling it through the impurity or scattering it away from a much larger region. This demonstrates the potential for quantum correlations to bridge length scales in the design of nanoelectronic devices and sensors.

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.

A multi-fragment real-time extension of projected density matrix embedding theory: Non-equilibrium electron dynamics in extended systems

In this work, we derive a multi-fragment real-time extension of the projected density matrix embedding theory (pDMET) designed to treat non-equilibrium electron dynamics in strongly correlated systems. As in the previously developed static pDMET, the real time pDMET partitions the total system into many fragments; the coupling between each fragment and the rest of the system is treated through a compact representation of the environment in terms of a quantum bath. The real-time pDMET involves simultaneously propagating the wavefunctions for each separate fragment–bath embedding system along with an auxiliary mean-field wavefunction of the total system. The equations of motion are derived by (i) projecting the time-dependent Schrödinger equation in the fragment and bath space associated with each separate fragment and by (ii) enforcing the pDMET matching conditions between the global 1-particle reduced density matrix (1-RDM) obtained from the fragment calculations and the mean-field 1-RDM at all points in time. The accuracy of the method is benchmarked through comparisons to time-dependent density-matrix renormalization group and time-dependent Hartree–Fock (TDHF) theory; the methods were applied to a one- and two-dimensional single-impurity Anderson model and multi-impurity Anderson models with ordered and disordered distributions of the impurities. The results demonstrate a large improvement over TDHF and rapid convergence to the exact dynamics with an increase in fragment size. Our results demonstrate that the real-time pDMET is a promising and flexible method that balances accuracy and efficiency to simulate the non-equilibrium electron dynamics in heterogeneous systems of large size.

Auxiliary Master Equation for Nonequilibrium Dual-Fermion Approach

We introduce an auxiliary quantum master equation dual fermion method and argue that it presents a convenient way to describe steady states of correlated impurity models. The scheme yields an expansion around a reference that is much closer to the true nonequilibrium state than that in the original dual fermion formulation. In steady-state situations, the scheme is numerically inexpensive and avoids time propagation. The Anderson impurity model is used to test the approach against numerically exact benchmarks.

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.

Kinetic Schemes in Open Interacting Systems

We discuss utilization of kinetic schemes for description of open interacting systems, focusing on vibrational energy relaxation for an oscillator coupled to a nonequilibirum electronic bath. Standard kinetic equations with constant rate coefficients are obtained under the assumption of time scale separation between the system and bath, with the bath dynamics much faster than that of the system of interest. This assumption may break down in certain limits, and we show that ignoring this may lead to qualitatively wrong predictions. Connection with more general, nonequilibrium Green's function (NEGF) analysis is demonstrated. Our considerations are illustrated within generic molecular junction models with electron-vibration coupling.

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.

Voltage Quench Dynamics of a Kondo System

We examine the dynamics of a correlated quantum dot in the mixed valence regime. We perform numerically exact calculations 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 equilibration times and saturates upon the decrease of temperature at all times, indicating Kondo behavior both in the transient regime and in the steady state. The time-dependent current saturation temperature connects the equilibrium Kondo temperature to a substantially increased value at voltages outside of the linear response. These signatures are directly observable by experiments in the time domain.

Taming the Dynamical Sign Problem in Real-Time Evolution of Quantum Many-Body Problems

Current nonequilibrium Monte Carlo methods suffer from a dynamical sign problem that makes simulating real-time dynamics for long times exponentially hard. We propose a new "inchworm algorithm," based on iteratively reusing information obtained in previous steps to extend the propagation to longer times. The algorithm largely overcomes the dynamical sign problem, changing the scaling from exponential to quadratic. We use the method to solve the Anderson impurity model in the Kondo and mixed valence regimes, obtaining results both for quenches and for spin dynamics in the presence of an oscillatory magnetic field.

Towards a full ab initio theory of strong electronic correlations in nanoscale devices

In this paper I give a detailed account of an ab initio methodology for describing strong electronic correlations in nanoscale devices hosting transition metal atoms with open d- or f-shells. The method combines Kohn-Sham density functional theory for treating the weakly interacting electrons on a static mean-field level with non-perturbative many-body methods for the strongly interacting electrons in the open d- and f-shells. An effective description of the strongly interacting electrons in terms of a multi-orbital Anderson impurity model is obtained by projection onto the strongly correlated subspace properly taking into account the non-orthogonality of the atomic basis set. A special focus lies on the ab initio calculation of the effective screened interaction matrix U for the Anderson model. Solution of the effective Anderson model with the one-crossing approximation or other impurity solver techniques yields the dynamic correlations within the strongly correlated subspace giving rise e.g. to the Kondo effect. As an example the method is applied to the case of a Co adatom on the Cu(0 0 1) surface. The calculated low-bias tunnel spectra show Fano-Kondo lineshapes similar to those measured in experiments. The exact shape of the Fano-Kondo feature as well as its width depend quite strongly on the filling of the Co 3d-shell. Although this somewhat hampers accurate quantitative predictions regarding lineshapes and Kondo temperatures, the overall physical situation can be predicted quite reliably.