Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2408.06554 · Physical review. E

Application of the shift-invert Lanczos algorithm to a nonequilibrium Green's function for transport problems

Kunihito Uzawa, K. Hagino

Published Nov 4, 2024• 2 citations
arXivHTMLPDF

Abstract

Nonequilibrium Green's function theory and related methods are widely used to describe transport phenomena in many-body systems, but they often require a costly inversion of a large matrix. We show here that the shift-invert Lanczos method can dramatically reduce the computational effort. We apply the method to two test problems, namely a simple model Hamiltonian and to a more realistic Hamiltonian for nuclear fission. For a Hamiltonian of dimension 66 103 we find that the computation time is reduced by a factor of 33 compared to the direct calculation of the Green's function.

DOI
10.1103/physreve.110.055302

Related papers

Krylov-space approach to the equilibrium and nonequilibrium single-particle Green’s function

Matthias Balzer, Nadine Gdaniec, Michael Potthoff

The zero-temperature single-particle Green's function of correlated fermion models with moderately large Hilbert-space dimensions can be calculated by means of Krylov-space techniques. The conventional Lanczos approach consists of finding the ground state in a first step, followed by an approximation for the resolvent of the Hamiltonian in a second step. We analyze the character of this approximation and discuss a numerically exact variant of the Lanczos method which is formulated in the time domain. This method is extended to obtain the nonequilibrium single-particle Green's function defined on the Keldysh-Matsubara contour in the complex time plane which describes the system's nonperturbative response to a sudden parameter switch in the Hamiltonian. The proposed method will be important as an exact-diagonalization solver in the context of self-consistent or variational cluster-embedding schemes. For the recently developed nonequilibrium cluster-perturbation theory, we discuss its efficient implementation and demonstrate the feasibility of the Krylov-based solver. The dissipation of a strong local magnetic excitation into a non-interacting bath is considered as an example for applications.

arXiv:1109.1205•Journal of Physics Condensed Matter•Dec 20, 2011•34 citations

Lanczos-based Low-Rank Correction Method for Solving the Dyson Equation in Inhomogenous Dynamical Mean-Field Theory

Pierre Carrier, Jok Man Tang, Yousef Saad +1 more

Inhomogeneous dynamical mean-field theory has been employed to solve many interesting strongly interacting problems from transport in multilayered devices to the properties of ultracold atoms in a trap. The main computational step, especially for large systems, is the problem of calculating the inverse of a large sparse matrix to solve Dyson's equation and determine the local Green's function at each lattice site from the corresponding local self-energy. We present a new e_cient algorithm, the Lanczos-based low-rank algorithm, for the calculation of the inverse of a large sparse matrix which yields this local (imaginary time) Green's function. The Lanczos-based low-rank algorithm is based on a domain decomposition viewpoint, but avoids explicit calculation of Schur complements and relies instead on low-rank matrix approximations derived from the Lanczos algorithm, for solving the Dyson equation. We report at least a 25-fold improvement of performance compared to explicit decomposition (such as sparse LU) of the matrix inverse. We also report that scaling relative to matrix sizes, of the low-rank correction method on the one hand and domain decomposition methods on the other, are comparable.

arXiv:1102.3137•Physics Procedia•Jan 1, 2011•2 citations

Time-dependent Green's Functions Approach to Nuclear Reactions

Arnau Rios Huguet, Paweł Danielewicz, Brent Barker

Nonequilibrium Green's functions represent underutilized means of studying the time evolution of quantum many-body systems. In view of a rising computer power, an effort is underway to apply the Green's functions formalism to the dynamics of central nuclear reactions. As the first step, mean-field evolution for the density matrix for colliding slabs is studied in one dimension. The strategy to extend the dynamics to correlations is described.

Bulletin of the American Physical Society•Oct 24, 2008•0 citations

Quantum Computed Green's Functions using a Cumulant Expansion of the Lanczos Method

Gabriel Greene‐Diniz, David Zsolt Manrique, Kentaro Yamamoto +5 more

In this paper, we present a quantum computational method to calculate the many-body Green's function matrix in a spin orbital basis. We apply our approach to finite-sized fermionic Hubbard models and related impurity models within Dynamical Mean Field Theory, and demonstrate the calculation of Green's functions on Quantinuum's H1-1 trapped-ion quantum computer. Our approach involves a cumulant expansion of the Lanczos method, using Hamiltonian moments as measurable expectation values. This bypasses the need for a large overhead in the number of measurements due to repeated applications of the variational quantum eigensolver (VQE), and instead measures the expectation value of the moments with one set of measurement circuits. From the measured moments, the tridiagonalised Hamiltonian matrix can be computed, which in turn yields the Green's function via continued fractions. While we use a variational algorithm to prepare the ground state in this work, we note that the modularity of our implementation allows for other (non-variational) approaches to be used for the ground state.

arXiv:2309.09685•Quantum•Jun 20, 2024•10 citations

NESSi 2.0: The Non-Equilibrium Systems Simulation package version 2.0

Fabian Künzel, Michael Schüler, Denis Golež +7 more

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.

arXiv:2605.02606•Computer Physics Communications•Oct 9, 2026•0 citations

Computing the real-time Green's Functions of large Hamiltonian matrices

T. Iitaka

A numerical method is developed for calculating the real time Green's functions of very large sparse Hamiltonian matrices, which exploits the numerical solution of the inhomogeneous time-dependent Schroedinger equation. The method has a clear-cut structure reflecting the most naive definition of the Green's functions, and is very suitable to parallel and vector supercomputers. The effectiveness of the method is illustrated by applying it to simple lattice models. An application of this method to condensed matter physics will be found in H. Tanaka, Phys. PRB 57, 2168 (1998).

arXiv:physics/9802021•Feb 11, 1998•3 citations
Nonequilibrium Green's function theory and related methods are widely used to describe transport phenomena in many-body systems, but they often require a costly inversion of a large matrix. We show here that the shift-invert Lanczos method can dramatically reduce the computational effort. We apply the method to two test problems, namely a simple model Hamiltonian and to a more realistic Hamiltonian for nuclear fission. For a Hamiltonian of dimension 66 103 we find that the computation time is reduced by a factor of 33 compared to the direct calculation of the Green's function.
The zero-temperature single-particle Green's function of correlated fermion models with moderately large Hilbert-space dimensions can be calculated by means of Krylov-space techniques. The conventional Lanczos approach consists of finding the ground state in a first step, followed by an approximation for the resolvent of the Hamiltonian in a second step. We analyze the character of this approximation and discuss a numerically exact variant of the Lanczos method which is formulated in the time domain. This method is extended to obtain the nonequilibrium single-particle Green's function defined on the Keldysh-Matsubara contour in the complex time plane which describes the system's nonperturbative response to a sudden parameter switch in the Hamiltonian. The proposed method will be important as an exact-diagonalization solver in the context of self-consistent or variational cluster-embedding schemes. For the recently developed nonequilibrium cluster-perturbation theory, we discuss its efficient implementation and demonstrate the feasibility of the Krylov-based solver. The dissipation of a strong local magnetic excitation into a non-interacting bath is considered as an example for applications.
Inhomogeneous dynamical mean-field theory has been employed to solve many interesting strongly interacting problems from transport in multilayered devices to the properties of ultracold atoms in a trap. The main computational step, especially for large systems, is the problem of calculating the inverse of a large sparse matrix to solve Dyson's equation and determine the local Green's function at each lattice site from the corresponding local self-energy. We present a new e_cient algorithm, the Lanczos-based low-rank algorithm, for the calculation of the inverse of a large sparse matrix which yields this local (imaginary time) Green's function. The Lanczos-based low-rank algorithm is based on a domain decomposition viewpoint, but avoids explicit calculation of Schur complements and relies instead on low-rank matrix approximations derived from the Lanczos algorithm, for solving the Dyson equation. We report at least a 25-fold improvement of performance compared to explicit decomposition (such as sparse LU) of the matrix inverse. We also report that scaling relative to matrix sizes, of the low-rank correction method on the one hand and domain decomposition methods on the other, are comparable.
In this paper, we present a quantum computational method to calculate the many-body Green's function matrix in a spin orbital basis. We apply our approach to finite-sized fermionic Hubbard models and related impurity models within Dynamical Mean Field Theory, and demonstrate the calculation of Green's functions on Quantinuum's H1-1 trapped-ion quantum computer. Our approach involves a cumulant expansion of the Lanczos method, using Hamiltonian moments as measurable expectation values. This bypasses the need for a large overhead in the number of measurements due to repeated applications of the variational quantum eigensolver (VQE), and instead measures the expectation value of the moments with one set of measurement circuits. From the measured moments, the tridiagonalised Hamiltonian matrix can be computed, which in turn yields the Green's function via continued fractions. While we use a variational algorithm to prepare the ground state in this work, we note that the modularity of our implementation allows for other (non-variational) approaches to be used for the ground state.
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.
A numerical method is developed for calculating the real time Green's functions of very large sparse Hamiltonian matrices, which exploits the numerical solution of the inhomogeneous time-dependent Schroedinger equation. The method has a clear-cut structure reflecting the most naive definition of the Green's functions, and is very suitable to parallel and vector supercomputers. The effectiveness of the method is illustrated by applying it to simple lattice models. An application of this method to condensed matter physics will be found in H. Tanaka, Phys. PRB 57, 2168 (1998).