No Arabic abstract
We present here the details of a method [A. B. Culver and N. Andrei, Phys. Rev. B 103, L201103 (2021)] for calculating the time-dependent many-body wavefunction that follows a local quench. We apply the method to the voltage-driven nonequilibrium Kondo model to find the exact time-evolving wavefunction following a quench where the dot is suddenly attached to the leads at $t=0$. The method, which does not use Bethe ansatz, also works in other quantum impurity models and may be of wider applicability. We show that the long-time limit (with the system size taken to infinity first) of the time-evolving wavefunction of the Kondo model is a current-carrying nonequilibrium steady state that satisfies the Lippmann-Schwinger equation. We show that the electric current in the time-evolving wavefunction is given by a series expression that can be expanded either in weak coupling or in strong coupling, converging to all orders in the steady-state limit in either case. The series agrees to leading order with known results in the well-studied regime of weak antiferromagnetic coupling and also reveals a universal regime of strong ferromagnetic coupling with Kondo temperature $T_K^{(F)} = D e^{-frac{3pi^2}{8} rho |J|}$ ($J<0$, $rho|J|toinfty$). In this regime, the differential conductance $dI/dV$ reaches the unitarity limit $2e^2/h$ asymptotically at large voltage or temperature.
We present a method for calculating the time-dependent many-body wavefunction that follows a local quench. We apply the method to the voltage-driven nonequilibrium Kondo model to find the exact time-evolving wavefunction following a quench where the dot is suddenly attached to the leads at $t=0$. The method, which does not use Bethe ansatz, also works in other quantum impurity models (we include results for the interacting resonant level and the Anderson impurity model) and may be of wider applicability. In the particular case of the Kondo model, we show that the long-time limit (with the system size taken to infinity first) of the time-evolving wavefunction is a current-carrying nonequilibrium steady state that satisfies the Lippmann-Schwinger equation. We show that the electric current in the time-evolving wavefunction is given by a series expression that can be expanded either in weak coupling or in strong coupling, converging to all orders in the steady-state limit in either case. The series agrees to leading order with known results in the well-studied regime of weak antiferromagnetic coupling and also reveals another universal regime of strong ferromagnetic coupling, with Kondo temperature $T_K^{(F)} = D e^{-frac{3pi^2}{8} rho |J|}$ ($J<0$, $rho|J|toinfty$). In this regime, the differential conductance $dI/dV$ reaches the unitarity limit $2e^2/h$ asymptotically at large voltage or temperature.
We extend the general formalism discussed in the previous paper [A. B. Culver and N. Andrei, Phys. Rev. B 103, 195106 (2021)] to two models with charge fluctuations: the interacting resonant level model and the Anderson impurity model. In the interacting resonant level model, we find the exact time-evolving wavefunction and calculate the steady state impurity occupancy to leading order in the interaction. In the Anderson impurity model, we find the nonequilibrium steady state for small or large Coulomb repulsion $U$, and we find that the steady state current to leading order in $U$ agrees with a Keldysh perturbation theory calculation.
We theoretically investigate the non-equilibrium quantum phase transition in a generic setup: the pseudogap Kondo model where a quantum dot couples to two-left (L) and right (R)-voltage-biased fermionic leads with power-law density of states (DOS) with respect to their Fermi levels {mu}_L/R, {rho}_c,L(R) ({omega}) propto |{omega} - {mu}_L(R) |r, and 0 < r < 1. In equilibrium (zero bias voltage) and for 0 < r < 1/2, with increasing Kondo correlations, in the presence of particle-hole symmetry this model exhibits a quantum phase transition from a unscreened local moment (LM) phase to the Kondo phase. Via a controlled frequency-dependent renormalization group (RG) approach, we compute analytically and numerically the non-equilibrium conductance, conduction electron T-matrix and local spin susceptibility at finite bias voltages near criticality. The current-induced decoherence shows distinct nonequilibrium scaling, leading to new universal non-equilibrium quantum critical behaviors in the above observables. Relevance of our results for the experiments is discussed.
Closed quantum many-body systems out of equilibrium pose several long-standing problems in physics. Recent years have seen a tremendous progress in approaching these questions, not least due to experiments with cold atoms and trapped ions in instances of quantum simulations. This article provides an overview on the progress in understanding dynamical equilibration and thermalisation of closed quantum many-body systems out of equilibrium due to quenches, ramps and periodic driving. It also addresses topics such as the eigenstate thermalisation hypothesis, typicality, transport, many-body localisation, universality near phase transitions, and prospects for quantum simulations.
A numerical approach is presented that allows to compute nonequilibrium steady state properties of strongly correlated quantum many-body systems. The method is imbedded in the Keldysh Greens function formalism and is based upon the idea of the variational cluster approach as far as the treatment of strong correlations is concerned. It appears that the variational aspect is crucial as it allows for a suitable optimization of a reference system to the nonequilibrium target state. The approach is neither perturbative in the many-body interaction nor in the field, that drives the system out of equilibrium, and it allows to study strong perturbations and nonlinear responses of systems in which also the correlated region is spatially extended. We apply the presented approach to non-linear transport across a strongly correlated quantum wire described by the fermionic Hubbard model. We illustrate how the method bridges to cluster dynamical mean-field theory upon coupling two baths containing and increasing number of uncorrelated sites.