Solving the Transient Dyson Equation with Quasilinear Complexity via Matrix Compression

Abstract

We introduce a numerical strategy to efficiently solve the out-of-equilibrium Dyson equation in the transient regime. By discretizing the equation into a compact matrix form and applying state-of-the-art matrix compression techniques, we achieve significant improvements in computational efficiency, which result in quasi-linear scaling of both time and space complexity with propagation time. This enables to compute accurate solutions even for systems with multiple and disparate time scales. We benchmark our solver by simulating a voltage-biased Josephson junction formed by a quantum dot connected to two superconducting leads.

References in corpus (16)

• Functional renormalization group approach to correlated fermion systems
• Resonant tunneling through a small quantum dot coupled to superconducting leads
• Diagrammatic Monte Carlo for Correlated Fermions
• NESSi: The Non-Equilibrium Systems Simulation package
• Quantum Monte-Carlo for correlated out-of-equilibrium nanoelectronics devices
• Low rank compression in the numerical solution of the nonequilibrium Dyson equation
• Non-equilibrium dynamics of Andreev states in the Kondo regime
• Tkwant: a software package for time-dependent quantum transport
• Legendre-spectral Dyson equation solver with super-exponential convergence
• Time-linear quantum transport simulations with correlated nonequilibrium Green's functions
• Anomalous Josephson current through a driven double quantum dot
• A fast time domain solver for the equilibrium Dyson equation
• A scalable numerical approach to the solution of the Dyson equation for the non-equilibrium single-particle Green's function
• Nonadiabatic Josephson current pumping by microwave irradiation
• Theory of non-equilibrium noise in general multi-terminal superconducting hydrid devices: application to multiple Cooper pair resonances
• Josephson radiation in a superconductor-quantum dot-superconductor junction