Imaginary-time formulation of steady-state nonequilibrium: application to strongly correlated transport

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.

Numerical study of non-adiabatic quantum thermodynamics of the driven resonant level model: non-equilibrium entropy production and higher order corrections

We present our numerical study on quantum thermodynamics of the resonant level model subjected to non-equilibrium condition as well as external driving. Following our previous work on non-equilibrium quantum thermodynamics (Douet al2020Phys. Rev.B101184304), we expand the density operator into a series of power in the driving speed, where we can determine the non-adiabatic thermodynamic quantities. Particularly, we calculate the non-equilibrium entropy production rate as well as higher order non-adiabatic corrections to the energy and/or population, which is not determined previously in Douet al(2020Phys. Rev.B101184304). In the limit of weak system-bath coupling, our results reduce to the one from the quantum master equation.

Analyzing scanning tunneling spectroscopy for Fe-based superconductors and extracting sample density of states

We extract the density of states (DOS) from the scanning tunneling spectroscopy data for Ba_1-x K _x Fe₂As₂ superconductor. The obtained sample DOS is composed of two ordinary s-wave types from the band at [Formula: see text] point and a linear-like DOS within the s-wave gap from the band at M point in the Brillouin zone, and is consistent with the corresponding data from angle-resolved photoemission spectroscopy. We clarify that the major peak of the tunneling conductance is not related to the DOS but is rather the effect of nonequilibrium coherent tunneling including all coherent spins in the tip and sample.

Auxiliary master equation approach within stochastic wave functions: Application to the interacting resonant level model

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.

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.

Nonequilibrium density matrix in quantum open systems: generalization for simultaneous heat and charge steady-state transport

We suggest a generalization of the expression of the nonequilibrium (NE) density matrix obtained by Hershfield's method for the cases where both heat and charge steady-state currents are present in a quantum open system. The finite-size quantum system, connected to two temperature and particle reservoirs, is driven out of equilibrium by the presence of both a temperature gradient and a chemical potential gradient between the two reservoirs. We show that the NE density matrix is given by a generalized Gibbs-like ensemble and is in full agreement with the general results of the McLennan-Zubarev nonequilibrium ensembles. The extra nonequilibrium terms are related to the entropy production in the system and characterize the fluxes of heat and particle. An explicit example, for the lowest-order expansion, is provide for a model system of noninteracting fermions.

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

The nonequilibrium spectral properties of the Anderson impurity model with a chemical potential bias are investigated within a numerically exact real-time quantum Monte Carlo formalism. The two-time correlation function is computed in a form suitable for nonequilibrium dynamical mean field calculations. Additionally, the evolution of the model's spectral properties are simulated in an alternative representation, defined by a hypothetical but experimentally realizable weakly coupled auxiliary lead. The voltage splitting of the Kondo peak is confirmed and the dynamics of its formation after a coupling or gate quench are studied. This representation is shown to contain additional information about the dot's population dynamics. Further, we show that the voltage-dependent differential conductance gives a reasonable qualitative estimate of the equilibrium spectral function, but significant qualitative differences are found including incorrect trends and spurious temperature dependent effects.

Kondo memory in driven strongly correlated quantum dots

We investigate the real-time current response of strongly correlated quantum dot systems under sinusoidal driving voltages. By means of an accurate hierarchical equations of motion approach, we demonstrate the presence of prominent memory effects induced by the Kondo resonance on the real-time current response. These memory effects appear as distinctive hysteresis line shapes and self-crossing features in the dynamic current-voltage characteristics, with concomitant excitation of odd-number overtones. They emerge as a cooperative effect of quantum coherence-due to inductive behavior-and electron correlations-due to the Kondo resonance. We also show the suppression of memory effects and the transition to classical behavior as a function of temperature. All these phenomena can be observed in experiments and may lead to novel quantum memory applications.

Nonequilibrium density matrix for quantum transport: Hershfield approach as a McLennan-Zubarev form of the statistical operator

In this paper, we formally demonstrate that the nonequilibrium density matrix developed by Hershfield for the steady state has the form of a McLennan-Zubarev nonequilibrium ensemble. The correction term in this pseudoequilibrium Gibbs-like ensemble is directly related to the entropy production in the quantum open system. The fact that both methods state that a nonequilibrium steady state can be mapped onto a pseudoequilibrium, permits us to develop nonequilibrium quantities from formal expressions equivalent to the equilibrium case. We provide an example: the derivation of a nonequilibrium distribution function for the electron population in a scattering region in the context of quantum transport.

Nonequilibrium dynamical mean-field theory: an auxiliary quantum master equation approach

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.

Effective projection theory in a correlated electron system

We propose an efficient method for nonperturbative calculation of Green's function in a correlated electron system. The key idea of the method is to project out irrelevant operators having zero norm in the ground state, which we refer to as effective projection theory. We apply the method to a mesoscopic Anderson model and show that for a given wavefunction ansatz, equations of motion can be closed only by relevant operators, allowing easy calculation of the zero-temperature Green's function. It turns out that the resulting Green's functions reproduce exact limits of both weak and strong interactions. The accuracy is also verified for small systems by comparison with exact diagonalization results, revealing that effective projection theory captures the essential correlated features in the entire regime of interactions.

Continuous-time quantum Monte Carlo and maximum entropy approach to an imaginary-time formulation of strongly correlated steady-state transport

Recently, Han and Heary [Phys. Rev. Lett. 99, 236808 (2007)] proposed an approach to steady-state quantum transport through mesoscopic structures, which maps the nonequilibrium problem onto a family of auxiliary quantum impurity systems subject to imaginary voltages. We employ continuous-time quantum Monte-Carlo solvers to calculate accurate imaginary time data for the auxiliary models. The spectral function is obtained from a maximum entropy analytical continuation in both Matsubara frequency and complexified voltage. To enable the analytical continuation we construct a kernel which is compatible with the analytical structure of the theory. While it remains a formidable task to extract reliable spectral functions from this unbiased procedure, particularly for large voltages, our results indicate that the method in principle yields results in agreement with those obtained by other methods.

Charge transport through single molecules, quantum dots and quantum wires

We review recent progress in the theoretical description of correlation and quantum fluctuation phenomena in charge transport through single molecules, quantum dots and quantum wires. Various physical phenomena are addressed, relating to cotunneling, pair-tunneling, adiabatic quantum pumping, charge and spin fluctuations, and inhomogeneous Luttinger liquids. We review theoretical many-body methods to treat correlation effects, quantum fluctuations, non-equilibrium physics, and the time evolution into the stationary state of complex nanoelectronic systems.

Steady-state nonequilibrium density of States of driven strongly correlated lattice models in infinite dimensions

An exact formalism for calculating the retarded and advanced Green's functions of strongly correlated lattice models in a uniform electric field is derived within dynamical mean-field theory. To illustrate the method, we solve for the nonequilibrium density of states of the Hubbard model in both the metallic and Mott-insulating phases at half-filling (with an arbitrary strength electric field) by employing the approximate numerical renormalization group as the impurity solver. This general approach can be applied to any strongly correlated lattice model in the limit of large dimensions.

Steady-state currents through nanodevices: a scattering-states numerical renormalization-group approach to open quantum systems

We propose a numerical renormalization group (NRG) approach to steady-state currents through nanodevices. A discretization of the scattering-states continuum ensures the correct boundary condition for an open quantum system. We introduce two degenerate Wilson chains for current carrying left- and right-moving electrons reflecting time-reversal symmetry in the absence of a finite bias V. We employ the time-dependent NRG to evolve the known steady-state density operator for a noninteracting junction into the density operator of the fully interacting nanodevice at finite bias. We calculate the differential conductance as function of V, T, and the external magnetic field.