We develop a time-local nonequilibrium Green's function formulation for real-time dynamics in quantum systems coupled to superconducting leads. The superconducting lead self-energy is a strongly frequency-dependent matrix in Nambu space, giving rise to nonlocal memory kernels in the time domain. This makes direct propagation of the Kadanoff-Baym (KB) equations computationally demanding. To overcome this difficulty, we extend the auxiliary-mode expansion, originally developed for normal-metal leads, to Nambu-space self-energies. This allows us to decompose the superconducting lead self-energy into a finite number of exponential modes and to transform the KB equations with memory integrals into a closed set of ordinary differential equations. The resulting time-local equations enable efficient real-time simulations under general time-dependent bias voltages, superconducting phases, and one-body Hamiltonians of the central system, while retaining the memory effects induced by superconducting leads. As an application, we analyze voltage-quench dynamics in a superconductor-quantum-dot-superconductor junction and show that, after a dc bias is suddenly applied, the system evolves through a transient regime and relaxes to an ac Josephson periodic steady state. The resulting periodic steady-state current agrees with the Floquet Green's function solution, validating the present real-time formulation.
We develop a time-local nonequilibrium Green's function formulation for real-time dynamics in quantum systems coupled to superconducting leads. The superconducting lead self-energy is a strongly frequency-dependent matrix in Nambu space, giving rise to nonlocal memory kernels in the time domain. This makes direct propagation of the Kadanoff-Baym (KB) equations computationally demanding. To overcome this difficulty, we extend the auxiliary-mode expansion, originally developed for normal-metal leads, to Nambu-space self-energies. This allows us to decompose the superconducting lead self-energy into a finite number of exponential modes and to transform the KB equations with memory integrals into a closed set of ordinary differential equations. The resulting time-local equations enable efficient real-time simulations under general time-dependent bias voltages, superconducting phases, and one-body Hamiltonians of the central system, while retaining the memory effects induced by superconducting leads. As an application, we analyze voltage-quench dynamics in a superconductor-quantum-dot-superconductor junction and show that, after a dc bias is suddenly applied, the system evolves through a transient regime and relaxes to an ac Josephson periodic steady state. The resulting periodic steady-state current agrees with the Floquet Green's function solution, validating the present real-time formulation.
We propose a time-dependent many-body approach to study the short-time dynamics of correlated electrons in quantum transport through nanoscale systems contacted to metallic leads. This approach is based on the time propagation of the Kadanoff-Baym equations for the nonequilibrium many-body Green's function of open and interacting systems out of equilibrium. An important feature of the method is that it takes full account of electronic correlations and embedding effects in the presence of time-dependent external fields, while at the same time satisfying the charge conservation law. The method further extends the Meir-Wingreen formula to the time domain for initially correlated states. We study the electron dynamics of a correlated quantum wire attached to two-dimensional leads exposed to a sudden switch on of a bias voltage using conserving many-body approximations at Hartree-Fock, second Born and GW level. We obtain detailed results for the transient currents, dipole moments, spectral functions, charging times, and the many-body screening of the quantum wire as well as for the time-dependent density pattern in the leads, and we show how the time dependence of these observables provides a wealth of information on the energy level structure of the quantum wire out of equilibrium. For moderate interaction strengths the second Born and GW results are in excellent agreement at all times. We find that many-body effects beyond the Hartree-Fock approximation have a large effect on the qualitative behavior of the system and lead to a bias-dependent gap closing and quasiparticle broadening, shortening of the transient times and washing out of the step features in the current-voltage curves.
The method of nonequilibrium Green functions (NEGFs) constitutes a solid framework to describe closed and open quantum systems. It captures situations far from equilibrium including arbitrary perturbations and correlations - effects higher than first order in the binary interactions. On the contrary, in the absence of external fields, it reduces to the imaginary-time Matsubara formalism. The abilities (and successes) of the method originate from its self-consistency and the conserving character of many-body approximations that are introduced by diagram technique. On this account, NEGFs have, now for a long time, brought insight into the dynamics of nuclear matter, non-ideal quantum gases and (semiconductor) plasmas. However, applications to localized, finite and strongly interacting systems have emerged only within the recent decade. The dynamical properties follow from quantum kinetic equations - the two-time Kadanoff-Baym equations for the one-particle NEGF - the numerical solution of which is the central topic of this thesis. The non-Markovian structure of these equations thereby inhibits an unlimited, gradual propagation of the nonequilibrium Green function: Any self-consistent, future dynamics of the system is tightly linked to the (complete) past by non-local self-energies. Further, the success in obtaining solutions is being diminished if the system under investigation is not isotropic or inhomogeneous and, in addition, the spatial topology (e.g., the charge carrier or spin density) needs to be resolved in time. On top of this, small isolated systems are generally not candidates showing fast decoherence which offers little hope for smooth and well-damped integration kernels. Despite these difficulties, the thesis at hand is devoted to the computation of the two-time NEGF for finite, closed systems of which ensembles of electrons bound by nuclei forming atoms or molecules are the most natural representatives. Here, in the spirit of the well-studied one-dimensional helium atom, the electron dynamics in model atoms and molecules is addressed in ab initio fashion. To render calculations possible in the first place, (i) a novel hybrid representation of the NEGF is introduced based on the conception of finite elements (FE) and the discrete variable representation (DVR), and (ii) a parallel algorithm is developed that allows for an efficient, well scalable distributed-memory computation of the full two-time Green function. The first point attributes plenty of flexibility and makes the method readily adjustable to different Hamiltonians avoiding the need to, by hand or in intricate numerical manner, calculate the whole bunch of matrix elements of one- and two-particle operators. Also, in consequence, it drastically simplifies the evaluation of self-energies. The second point makes the NEGF approach ready for large-scale, high-performance computing facilities, the resources (of computer power: mainly memory!) of which are indispensable for (some of) the calculations presented. The applicability of the developed methods is demonstrated inter alia on the basis of helium, beryllium, hydrogen and lithium hydride modeled in one spatial dimension whereby electron-electron correlations are treated in the second(-order) Born approximation. After a thorough discussion on the self-consistent ground states, the thesis focuses on the electron dynamics in the linear response regime as well as in the presence of strong external fields. Overall, to assess the the performance of the second Born approximation, the results are compared to the Hartree-Fock approximation, on the one hand, and to the exact solution of the time-(in)dependent Schrödinger equation, on the other hand. In addition, aiming at the possibility of extended time-dependent calculations, the behavior of the generalized Kadanoff-Baym ansatz is discussed.
The time evolution in quantum many-body systems after external excitations is attracting high interest in many fields, including dense plasmas, correlated solids, laser-excited materials, or fermionic and bosonic atoms in optical lattices. The theoretical modeling of these processes is challenging, and the only rigorous quantum-dynamics approach that can treat correlated fermions in two and three dimensions is nonequilibrium Green functions (NEGF). However, NEGF simulations are computationally expensive due to their T3 scaling with the simulation duration T. Recently, T2 scaling was achieved with the generalized Kadanoff-Baym ansatz (GKBA), for the second-order Born (SOA) self energy, which has substantially extended the scope of NEGF simulations. In a recent Letter [Schl\"unzen et al., Phys. Rev. Lett. 124, 076601 (2020).] we demonstrated that GKBA-NEGF simulations can be efficiently mapped onto coupled time-local equations for the single-particle and two-particle Green functions on the time diagonal, hence the method has been called the G1-G2 scheme. This allows one to perform the same simulations with order T1 scaling, both for SOA and GW self energies giving rise to a dramatic speedup. Here we present more details on the G1-G2 scheme, including derivations of the basic equations including results for a general basis, for Hubbard systems, and for jellium. Also, we demonstrate how to incorporate initial correlations into the G1-G2 scheme. Further, the derivations are extended to a broader class of self energies, including the T matrix in the particle-particle and particle-hole channels and the dynamically-screened-ladder approximation. Finally, we demonstrate that, for all self energies, the CPU-time scaling of the G1-G2 scheme with the basis dimension Nb can be improved compared to our first report: The overhead compared to the original GKBA is not more than an additional factor Nb, even for Hubbard systems.
Nonequilibrium Green’s functions provide a powerful framework for studying quantum many-body dynamics including the laser-induced dynamics in solids. The N on- E quilibrium S ystems Si mulation package ( NESSi ) offers an efficient platform for such simulations, ranging from perturbative approaches like nonequilibrium GW to nonequilibrium dynamical mean-field theory. However, simulations based on nonequilibrium Green’s functions become computationally demanding when the dynamics span a large temporal range, such as from sub-femtosecond electron dynamics to the picosecond dynamics of collective modes. Due to the memory integral in the Kadanoff-Baym equations, which serve as equations of motion for nonequilibrium Green’s functions, the computational cost scales as 𝒪 ( 𝑁 3 𝑡 ) with the number of timesteps N t , and the memory requirement scales as 𝒪 ( 𝑁 2 𝑡 ) . In this work, we extend NESSi by incorporating techniques that aim to overcome this bottleneck: (i) By truncating the memory integrals in the KBE to a maximum of N c timesteps, the computational complexity is reduced to 𝒪 ( 𝑁 𝑡 𝑁 2 𝑐 ) , and the memory requirement to 𝒪 ( 𝑁 2 𝑐 ) . Provided that the results converge with respect to the cutoff N c , memory truncation allows to extend the simulations to significantly longer times. (ii) We introduce functionalities to describe nonequilibrium steady states, i.e. time-translationally invariant nonequilibrium states. Such states are relevant for transport settings, and they provide an approximate description of slowly evolving (prethermal) nonequilibrium states.
The nonequilibrium Dyson (or Kadanoff-Baym) equation, which is an equation of motion with a long-range memory kernel for real-time Green functions, underlies many numerical approaches based on the Keldysh formalism. In this paper we map the problem of solving the Dyson equation in real time onto a noninteracting auxiliary Hamiltonian with additional bath degrees of freedom. The solution of the auxiliary model does not require the evaluation of a memory kernel and can thus be implemented in a very memory efficient way. The mapping is derived for a self-energy which is local in space and is thus directly applicable within nonequilibrium dynamical mean-field theory (DMFT). We apply the method to study the interaction quench in the Hubbard model for an optical lattice with a narrow confinement, using inhomogeneous DMFT in combination with second-order weak-coupling perturbation theory. We find that, although the quench excites pronounced density oscillations, signatures of the two-stage relaxation similar to the homogeneous system can be observed by looking at the time-dependent occupations of natural orbitals.
We have developed a method based on the embedded Kadanoff-Baym equations to study the time evolution of open and inhomogeneous systems. The equation of motion for the Green's function on the Keldysh contour is solved using different conserving many-body approximations for the self-energy. Our formulation incorporates basic conservation laws, such as particle conservation, and includes both initial correlations and initial embedding effects, without restrictions on the time-dependence of the external driving field. We present results for the time-dependent density, current and dipole moment for a correlated tight binding chain connected to one-dimensional non-interacting leads exposed to DC and AC biases of various forms. We find that the self-consistent 2B and GW approximations are in extremely good agreement with each other at all times, for the long-range interactions that we consider. In the DC case we show that the oscillations in the transients can be understood from interchain and lead-chain transitions in the system and find that the dominant frequency corresponds to the HOMO-LUMO transition of the central wire. For AC biases with odd inversion symmetry odd harmonics to high harmonic order in the driving frequency are observed in the dipole moment, whereas for asymmetric applied bias also even harmonics have considerable intensity. In both cases we find that the HOMO-LUMO transition strongly mixes with the harmonics leading to harmonic peaks with enhanced intensity at the HOMO-LUMO transition energy.