Time-dependent Numerical Renormalization Group Research Articles

Dynamical quantum phase transitions in a mesoscopic superconducting system

We inspect signatures of dynamical quantum phase transitions driven by two types of quenches acting on a correlated quantum dot embedded between superconducting and metallic reservoirs. Under stationary conditions the proximity induced on-dot pairing combined with the strong Coulomb repulsion prefers the quantum dot to be either in the singly occupied (spinful) or BCS-type (spinless) ground state configuration. We study the time evolution upon traversing such a phase boundary due to quantum quenches by means of the time-dependent numerical renormalization group approach, revealing non-analytic features in the low-energy return rate. Quench protocols can be realized in a controllable manner and we are confident that detection of this dynamical singlet-doublet phase transition would be feasible by the charge tunnelling spectroscopy.

Self-energy method for time-dependent spectral functions of the Anderson impurity model within the time-dependent numerical renormalization group approach

The self-energy method for quantum impurity models expresses the correlation part of the self-energy in terms of the ratio of two Green's functions and allows for a more accurate calculation of equilibrium spectral functions than is possible directly from the one-particle Green's function [Bulla et al., J. Phys.: Condens. Matter 10, 8365 (1998)], for example, within the numerical renormalization group method. In addition, the self-energy itself is a central quantity required in the dynamical mean field theory of strongly correlated lattice models. Here, we show how to generalize the self-energy method to the time-dependent situation for the prototype model of strong correlations, the Anderson impurity model. We use the equation of motion method to obtain closed expressions for the local Green's function in terms of a time-dependent correlation self-energy, with the latter being given as a ratio of a one-particle time-dependent Green's function and a higher-order correlation function. We benchmark this self-energy approach to time-dependent spectral functions against the direct approach within the time-dependent numerical renormalization group method. The self-energy approach improves the accuracy of time-dependent spectral function calculations, and the closed-form expressions for the Green's function allow for a clear picture of the time-evolution of spectral features at the different characteristic time-scales. The self-energy approach is of potential interest also for other quantum impurity solvers for real-time evolution, including time-dependent density matrix renormalization group and continuous-time quantum Monte Carlo techniques.

Transient effects in a double quantum dot sandwiched laterally between superconducting and metallic leads

We study the transient phenomena appearing in a subgap region of the double quantum dot coupled in series between the superconducting and normal metallic leads, focusing on the development of the superconducting proximity effect. For the uncorrelated nanostructure we derive explicit expressions of the time-dependent occupancies in both quantum dots, charge currents, and electron pairing induced on individual dots and between them. We show that the initial configurations substantially affect the dynamical processes, in which the in-gap bound states emerge upon coupling the double quantum dot to superconducting reservoir. In particular, the superconducting proximity effect would be temporarily blocked whenever the quantum dots are initially singly occupied. Such {\it triplet}/{\it Andreev blockade} has been recently reported experimentally for double quantum dots embedded in the Josephson [D. Bouman et al., Phys. Rev. B 102, 220505 (2020)] and Andreev [P. Zhang et al., arXiv:2102.03283 (2021)] junctions. We also address the role of correlation effects within the lowest-order decoupling scheme and by the time-dependent numerical renormalization group calculations. Competition of the repulsive Coulomb interactions with the superconducting proximity effect leads to renormalization of the in-gap quasiparticles, speeding up the quantum oscillations and narrowing a region of transient phenomena, whereas the dynamical Andreev blockade is well pronounced in the weak inter-dot coupling limit. We propose feasible methods for detecting the characteristic time-scales that could be observable by the Andreev spectroscopy.

Restoring the continuum limit in the time-dependent numerical renormalization group approach

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 non-equilibrium 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 $\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 quantum-impurity systems. An analytical solution of the resonant-level model serves as a benchmark for the accuracy of the method which is then applied to non-trivial models, such as the interacting resonant-level model and the single impurity Anderson model.

Analytical and numerical study of the out-of-equilibrium current through a helical edge coupled to a magnetic impurity

We study the conductance of a time-reversal symmetric helical electronic edge coupled antiferromagnetically to a magnetic impurity, employing analytical and numerical approaches. The impurity can reduce the perfect conductance $G_0$ of a noninteracting helical edge by generating a backscattered current. The backscattered steady-state current tends to vanish below the Kondo temperature $T_K$ for time-reversal symmetric setups. We show that the central role in maintaining the perfect conductance is played by a global $U(1)$ symmetry. This symmetry can be broken by an anisotropic exchange coupling of the helical modes to the local impurity. Such anisotropy, in general, dynamically vanishes during the renormalization group (RG) flow to the strong coupling limit at low-temperatures. The role of the anisotropic exchange coupling is further studied using the time-dependent Numerical Renormalization Group (TD-NRG) method, uniquely suitable for calculating out-of-equilibrium observables of strongly correlated setups. We investigate the role of finite bias voltage and temperature in cutting the RG flow before the isotropic strong-coupling fixed point is reached, extract the relevant energy scales and the manner in which the crossover from the weakly interacting regime to the strong-coupling backscattering-free screened regime is manifested. Most notably, we find that at low temperatures the conductance of the backscattering current follows a power-law behavior $G\sim (T/T_K)^2$, which we understand as a strong nonlinear effect due to time-reversal symmetry breaking by the finite-bias.

Time-dependent spectral functions of the Anderson impurity model in response to a quench with application to time-resolved photoemission spectroscopy

We investigate several definitions of the time-dependent spectral function $A(\omega,t)$ of the Anderson impurity model following a quench and within the time-dependent numerical renormalization group method. In terms of the two-time retarded Green function $G^r(t_1,t_2)$, the definitions differ in the choice of the time variable $t$ with respect to $t_1$ and/or $t_2$. In a previous study [Nghiem {\it et al.} Phys. Rev. Lett. 119, 156601 (2017)], we investigated the spectral function, obtained from the Fourier transform of ${\rm Im}[G^r(t_1,t_2)]$ w.r.t. the time difference $t'=t_1-t_2$, with $t=t_2$. Here, we derivie expressions for the retarded Green function for the choices $t=t_1$ and the average time $t=(t_1+t_2)/2$, within the TDNRG approach. We compare and contrast the resulting $A(\omega,t)$ for the different choices of time reference. Expressions for the lesser, greater and advanced Green functions are also derived within TDNRG for all choices of time reference. The average time lesser Green function $G^<(\omega,t)$ is particularly interesting, as it determines the time-dependent occupied density of states $N(\omega,t)=G^<(\omega,t)/(2\pi i)$, a quantity that determines the photoemission current in time-resolved pump-probe photoemission spectroscopy. We present calculations for $N(\omega,t)$ for the Anderson model following a quench, and discuss the resulting time evolution of the spectral features, such as the Kondo resonance and high-energy peaks. We also discuss the issue of thermalization at long times for $N(\omega,t)$. Finally, we use the results for $N(\omega,t)$ to calculate the time-resolved photoemission current for the Anderson model following a quench (acting as the pump) and study the different behaviors that can be observed for different resolution times of a Gaussian probe pulse.

Time-dependent numerical renormalization group method for multiple quenches: Towards exact results for the long-time limit of thermodynamic observables and spectral functions

We develop an alternative time-dependent numerical renormalization group (TDNRG) formalism for multiple quenches and implement it to study the response of a quantum impurity system to a general pulse. Within this approach, we reduce the contribution of the NRG approximation to numerical errors in the time evolution of observables by a formulation that avoids the use of the generalized overlap matrix elements in our previous multiple-quench TDNRG formalism [Nghiem {\em et al.,} Phys. Rev. B {\bf 89}, 075118 (2014); Phys. Rev. B {\bf 90}, 035129 (2014)]. We demonstrate that the formalism yields a smaller cumulative error in the trace of the projected density matrix as a function of time and a smaller discontinuity of local observables between quenches than in our previous approach. Moreover, by increasing the switch-on time, the time between the first and last quench of the discretized pulse, the long-time limit of observables systematically converges to its expected value in the final state, i.e., the more adiabatic the switching, the more accurately is the long-time limit recovered. The present formalism can be straightforwardly extended to infinite switch-on times. We show that this yields highly accurate results for the long-time limit of both thermodynamic observables and spectral functions, and overcomes the significant errors within the single quench formalism [Anders {\em et al.}, Phys. Rev. Lett. {\bf 95}, 196801 (2005); Nghiem {\em et al.}, Phys. Rev. Lett. {\bf 119}, 156601 (2017)]. This improvement provides a first step towards an accurate description of nonequilibrium steady states of quantum impurity systems, e.g., within the scattering states NRG approach [Anders, Phys. Rev. Lett. {\bf 101}, 066804 (2008)].

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.

Dynamics of the two-spin spin-boson model with a common bath.

Dynamics of the two-spin spin-boson model in the presence of Ohmic and sub-Ohmic baths is investigated by employing a multitude of the Davydov D1 trial states, also known as the multi-D1 Ansatz. Its accuracy in dynamics simulations of the two-spin SBM is improved significantly over the single D1 Ansatz, especially in the weak to moderately strong coupling regime. Validity of the multi-D1 Ansatz for various coupling strengths is also systematically examined by making use of the deviation vector which quantifies how faithfully the trial state obeys the Schrödinger equation. The time evolution of population difference and entanglement has been studied for various initial conditions and coupling strengths. Careful comparisons are carried out between our approach and three other methods, i.e., the time-dependent numerical renormalization group (TD-NRG) approach, the Bloch-Redfield theory, and a method based on a variational master equation. For strong coupling, the multi-D1 trial state yields consistent results as the TD-NRG approach in the Ohmic regime while the two disagree in the sub-Ohmic regime, where the multi-D1 trial state is shown to be more accurate. For weak coupling, the multi-D1 trial state agrees with the two master-equation methods in the presence of both Ohmic and sub-Ohmic baths, but shows considerable differences with the TD-NRG approach in the presence of a sub-Ohmic bath, calling into question the validity of the TD-NRG results at long times in the literature.

Real-time dynamics induced by quenches across the quantum critical points in gapless Fermi systems with a magnetic impurity

The energy-dependent scattering of fermions from a localized orbital at an energy-dependent rate $\ensuremath{\Gamma}(\ensuremath{\varepsilon})\ensuremath{\propto}{|\ensuremath{\varepsilon}|}^{r}$ gives rise to quantum critical points (QCPs) in the pseudogap single-impurity Anderson model separating a local moment phase with an unscreened spin moment from a strong-coupling phase which slightly deviates from the screened phase of standard Kondo problem. Using the time-dependent numerical renormalization group (TD-NRG) approach we show that local dynamic properties always equilibrate towards a steady-state value even for quenches across the QCP but with systematic deviations from the thermal equilibrium depending on the distance to the critical coupling. Local nonequilibrium properties are presented for interaction quenches and hybridization quenches. We augment our numerical data by an analytical calculation that becomes exact at short times and find excellent agreement between the numerics and the analytical theory. For interaction quenches within the screened phase we find a universal function for the time-dependent local double occupancy. We trace back the discrepancy between our results and the data obtained by a time-dependent Gutzwiller variational approach to restrictions of the wave-function ansatz in the Gutzwiller theory: while the NRG ground states properly account for the formation of an extended spin moment which decouples from the system in the unscreened phase, the Gutzwiller ansatz only allows the formation of the spin moment on the local impurity orbital.

Time-dependent numerical renormalization group method for multiple quenches: Application to general pulses and periodic driving

The time-dependent numerical renormalization group method (TDNRG) [Anders et al., Phys. Rev. Lett. {\bf 95}, 196801 (2005)] was recently generalized to multiple quenches and arbitrary finite temperatures [Nghiem et al., Phys. Rev. B {\bf 89}, 075118 (2014)] by using the full density matrix approach [Weichselbaum et al., Phys. Rev. Lett. {\bf 99}, 076402 (2007)]. In this paper, we numerically implement this formalism to study the response of a quantum impurity system to a general pulse and periodic driving which are approximated by a sufficient number of quenches. We show how the NRG approximation affects the trace of the projected density matrices and the continuity of the time-evolution of a local observable. For the general pulse case, the local observable in the long-time limit exhibits a dependence on the switch-on time, the time interval between the first and last quenches, as well as on the pulse shape. In particular, the long-time limit is improved for longer switch-on times and smoother pulses. This lends support to our earlier suggestion that the long-time limit of observables can be improved by replacing a sudden large quench by a sequence of smaller ones acting over a finite time-interval: longer switch-on times and smoother pulses, i.e., increased adiabaticity, favor relaxation of the system to its correct thermodynamic long-time limit. For the case of periodic driving, we compare the TDNRG results to exact analytic ones for the non-interacting resonant level model, finding better agreement at short to intermediate time scales in the case of smoother driving. Finally, we demonstrate the validity of the multiple-quench TDNRG formalism for arbitrary temperatures by studying the time-evolution of the occupation number in the Anderson impurity model in response to a periodic switching of the local level from the mixed valence to the Kondo regime at finite temperatures.

Generalization of the time-dependent numerical renormalization group method to finite temperatures and general pulses

The time-dependent numerical renormalization group (td-NRG) [Anders et al. Phys. Rev. Lett. {\bf 95}, 196801 (2006)] offers the prospect of investigating in a non-perturbative manner the time-dependence of local observables of interacting quantum impurity models at all time scales following a quantum quench. We present a generalization of this method to arbitrary finite temperature by making use of the full density matrix [Weichselbaum et al. Phys. Rev. Lett. {\bf 99}, 076402 (2007)]. We show that all terms in the projected density matrix $\rho^{i\to f} = \rho^{++} + \rho^{--} + \rho^{+-}+\rho^{-+}$ in the time-evolution of a local observable may be evaluated in closed form, with $\rho^{+-}=\rho^{-+}=0$. The expression for $\rho^{--}$ is shown to be finite at finite temperature, becoming negligible only in the limit of vanishing temperatures. We prove that this approach recovers the short-time limit for the expectation value of a local observable exactly at arbitrary temperatures. In contrast, the corresponding long-time limit is only recovered exactly in the limit $\Lambda\rightarrow 1^{+}$. This limit is impractical within NRG. We suggest how to overcome this problem by noting that the long-time limit becomes increasingly more accurate on reducing the size of the quantum quench. This suggests an improved generalized td-NRG approach in which the system is time-evolved between the initial and final states via a sequence of small quantum quenches within a finite time interval instead of a single large quantum quench. The formalism for this is provided, generalizing the td-NRG to multiple quantum quenches, periodic switching, and general pulses. This formalism, like our finite temperature generalization of the single-quench case, rests only on the NRG approximation. The results are illustrated by application to the Anderson impurity model.

Tensor networks and the numerical renormalization group

The full-density-matrix numerical renormalization group has evolved as a systematic and transparent setting for the calculation of thermodynamical quantities at arbitrary temperatures within the numerical renormalization group (NRG) framework. It directly evaluates the relevant Lehmann representations based on the complete basis sets introduced by Anders and Schiller [Phys. Rev. Lett. 95, 196801 (2005)]. In addition, specific attention is given to the possible feedback from low-energy physics to high energies by the explicit and careful construction of the full thermal density matrix, naturally generated over a distribution of energy shells. Specific examples are given in terms of spectral functions (fdmNRG), time-dependent NRG (tdmNRG), Fermi-golden-rule calculations (fgrNRG) as well as the calculation of plain thermodynamic expectation values. Furthermore, based on the very fact that, by its iterative nature, the NRG eigenstates are naturally described in terms of matrix product states, the language of tensor networks has proven enormously convenient in the description of the underlying algorithmic procedures. This paper therefore also provides a detailed introduction and discussion of the prototypical NRG calculations in terms of their corresponding tensor networks.

Quantum transport through a quantum dot: Combining the scattering-states numerical renormalization group with nonequilibrium Green functions

Scattering states fulfill the correct boundary conditions of a current carrying open quantum system. Discretizing the energy continuum of these states allows for employing Wilson's numerical renormalization group approach without violating the boundary conditions by using a finite size system. We evolve the analytically known steady-state density operator for a non-interacting quantum-system at finite bias to the full interacting problem by the time-dependent numerical renormalization group after switching on the local charging energy. Using a newly developed algorithm for steady-state nonequilibrium Green functions, we can calculate the current I as function of bias voltage V for arbitrary temperature and magnetic field. A comparison with second-order and GW Kadanoff-Baym-Keldysh results shows excellent agreement for weak interaction strength U.

Comparative study of theoretical methods for non-equilibrium quantum transport

We present a detailed comparison of three different methods designed to tackle non-equilibrium quantum transport, namely the functional renormalization group (fRG), the time-dependent density matrix renormalization group (tDMRG) and the iterative summation of real-time path integrals (ISPI). For the non-equilibrium single-impurity Anderson model (including a Zeeman term at the impurity site), we demonstrate that the three methods are in quantitative agreement over a wide range of parameters at the particle–hole symmetric point as well as in the mixed-valence regime. We further compare these techniques with two quantum Monte Carlo approaches and the time-dependent numerical renormalization group method.

Nonequilibrium Spin Dynamics in the Ferromagnetic Kondo Model

Motivated by recent experiments on molecular quantum dots we investigate the relaxation of pure spin states when coupled to metallic leads. Under suitable conditions these systems are well described by a ferromagnetic Kondo model. Using two recently developed theoretical approaches, the time-dependent numerical renormalization group and an extended flow equation method, we calculate the real-time evolution of a Kondo spin into its partially screened steady state. We obtain exact analytical results which agree well with numerical implementations of both methods. Analytical expressions for the steady state magnetization and the dependence of the long-time relaxation on microscopic parameters are established. We find the long-time relaxation process to be much faster in the regime of anisotropic Kondo couplings. The steady state magnetization is found to deviate significantly from its thermal equilibrium value.

Spin precession and real-time dynamics in the Kondo model: Time-dependent numerical renormalization-group study

A detailed derivation of the recently proposed time-dependent numerical renormalization-group (TD-NRG) approach to nonequilibrium dynamics in quantum impurity systems is presented. We demonstrate that the method is suitable for fermionic as well as bosonic baths. A comparison with exact analytical results for the charge relaxation in the resonant-level model and for dephasing in the spin-boson model establishes the accuracy of the method. The real-time dynamics of a single spin coupled to both types of baths is investigated. We use the TD-NRG to calculate the spin relaxation and spin precession of a single Kondo impurity. The short- and long-time dynamics is studied as a function of temperature in the ferromagnetic and antiferromagnetic regimes. The short-time dynamics agrees very well with analytical results obtained at second order in the exchange coupling $J$. In the ferromagnetic regime, the long-time spin decay is described by the scaling variable $x = 2\rho_F J(T) T t$. In the antiferromagnetic regime it is governed for $T < T_K$ by the Kondo time scale $1/T_K$. Here $\rho_F$ is the conduction-electron density of states and $T_K$ is the Kondo temperature. Results for spin precession are obtained by rotating the external magnetic field from the x axis to the z axis.