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:1409.3374

A Davidson-Lanczos iteration method for computation of continued-fraction expansion of the Green's function at very low temperatures: Applications to the dynamical mean field theory

Medha Sharma, M. A. H. Ahsan

Published Sep 11, 2014• 0 citations
arXivHTMLPDF

Abstract

We present a combination method based on orignal version of Davidson algorithm for extracting few of the lowest eigenvalues and eigenvectors of a sparse symmetric Hamiltonian matrix and the simplest version of Lanczos technique for obtaining a tridiagonal representation of the Hamiltonian to compute the continued fraction expansion of the Green's function at a very low temperature. We compare the Davidson$+$Lanczos method with the full diagonalization on a one-band Hubbard model on a Bethe lattice of infinite-coordination using dynamical mean field theory.

DOI
10.48550/arxiv.1409.3374

Related papers

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

Augmented hybrid exact-diagonalization solver for dynamical mean field theory

Cédric Weber, A. Amaricci, Massimo Capone +1 more

We present a methodology to solve the Anderson impurity model in the context of dynamical mean field theory based on the exact-diagonalization method. We propose a strategy to effectively refine the exact-diagonalization solver by combining a finite-temperature Lanczos algorithm with an adapted version of the cluster perturbation theory. We show that the augmented diagonalization yields an improved accuracy in the description of the spectral function of the single-band Hubbard model and is a reliable approach for a full $d$-orbital manifold calculation.

arXiv:1207.5610•Physical Review B•Sep 26, 2012•30 citations

Solving the dynamical mean-field theory at very low temperatures using the Lanczos exact diagonalization

Massimo Capone, Luca de’ Medici, Antoine Georges

We present an efficient method to solve the impurity Hamiltonians involved in dynamical mean-field theory at low but finite temperature based on the extension of the Lanczos algorithm from ground state properties alone to excited states. We test the approach on the prototypical Hubbard model and find extremely accurate results from $T=0$ up to relatively high temperatures up to the scale of the critical temperature for the Mott transition. The algorithm substantially decreases the computational effort involved in finite temperature calculations.

arXiv:cond-mat/0512484•Physical Review B•Dec 17, 2007•119 citations

Quantum subspace expansion algorithm for Green's functions

François Jamet, Abhishek Agarwal, Ivan Rungger

We present an algorithm to compute Green's functions on quantum computers for interacting electron systems, which is a challenging task on conventional computers. It uses a continued fraction representation based on the Lanczos method, where the wave functions are expanded as linear combination of basis states within a quantum subspace. While on conventional computers the cost of the computation grows exponentially with system size, limiting the method to small systems, by representing the basis states on a quantum computer one may overcome this exponential scaling barrier. We propose a two-level multigrid Trotter time evolution for an efficient preparation of the basis states in a quantum circuit, which takes advantage of the robustness of the subspace expansion against Trotter errors. Using a quantum emulator we demonstrate the algorithm for the Hubbard model on a Bethe lattice with infinite coordination, which we map to a 16 qubit Anderson impurity model within the dynamical mean field theory. Our algorithm computes the Green's function accurately for both the metallic and Mott insulating regimes, with a circuit depth several orders of magnitude below what has been proposed using time evolution. The two-level multigrid time evolution reduces the number of Trotter steps required to compute the Green's function to about four to six. We therefore expect that the method can be used on near term quantum computers for moderate system sizes, while allowing for scalability to larger circuit depths and qubit numbers on future fault tolerant quantum computers.

arXiv:2205.00094•Apr 29, 2022•8 citations

Spectral density of the Hubbard model by the continued fraction method

R. Hayn, P. Lombardo, Konrad Matho

We present the continued fraction method (CFM) as a microscopic approximation to the spectral density of the Hubbard model in the correlated metal phase away from half filling. The quantity expanded as a continued fraction is the single-particle Green's function. Leading spectral moments are taken into account through a set of real expansion coefficients, as known from the projection technique. Further stages are added to the continued fraction, with complex coefficients, thus defining a terminator function. This enables us to treat the entire spectral range of the Green's function on equal footing and determine the energy scale of the Fermi liquid quasiparticles by minimizing the total energy. The solution is free of phenomenological parameters and remains well defined in the strong-coupling limit, near the doping-controlled metal-insulator transition. Our results for the density of states agree reasonably with those from several variants of the dynamical mean-field theory. The CFM requires minimal numerical effort and can be generalized in several ways that are interesting for applications to real materials.

arXiv:cond-mat/0605109•Physical Review B•Nov 30, 2006•7 citations

Adaptive Lanczos-vector method for dynamic properties within the density matrix renormalization group

Piet E. Dargel, A. Honecker, Robert Peters +2 more

Current widely used approaches to calculate spectral functions using the density-matrix renormalization group in frequency space either necessarily include an artificial broadening (correction-vector method), have limited resolution (time-domain density-matrix renormalization group with Fourier transform method), or are limited to low-energy properties or single dominant modes (original continued fraction method). Here we propose an adaptive Lanczos-vector method to calculate the coefficients of a continued fraction expansion of the spectral function iteratively. We show that one can obtain a very accurate representation of the spectral function very efficiently, and that one can also directly extract the spectral weights and poles for the discrete system.

arXiv:1012.5543•Physical Review B•Apr 22, 2011•37 citations
We present a combination method based on orignal version of Davidson algorithm for extracting few of the lowest eigenvalues and eigenvectors of a sparse symmetric Hamiltonian matrix and the simplest version of Lanczos technique for obtaining a tridiagonal representation of the Hamiltonian to compute the continued fraction expansion of the Green's function at a very low temperature. We compare the Davidson+++Lanczos method with the full diagonalization on a one-band Hubbard model on a Bethe lattice of infinite-coordination using dynamical mean field theory.
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.
We present a methodology to solve the Anderson impurity model in the context of dynamical mean field theory based on the exact-diagonalization method. We propose a strategy to effectively refine the exact-diagonalization solver by combining a finite-temperature Lanczos algorithm with an adapted version of the cluster perturbation theory. We show that the augmented diagonalization yields an improved accuracy in the description of the spectral function of the single-band Hubbard model and is a reliable approach for a full ddd-orbital manifold calculation.
We present an efficient method to solve the impurity Hamiltonians involved in dynamical mean-field theory at low but finite temperature based on the extension of the Lanczos algorithm from ground state properties alone to excited states. We test the approach on the prototypical Hubbard model and find extremely accurate results from T=0T=0T=0 up to relatively high temperatures up to the scale of the critical temperature for the Mott transition. The algorithm substantially decreases the computational effort involved in finite temperature calculations.
We present an algorithm to compute Green's functions on quantum computers for interacting electron systems, which is a challenging task on conventional computers. It uses a continued fraction representation based on the Lanczos method, where the wave functions are expanded as linear combination of basis states within a quantum subspace. While on conventional computers the cost of the computation grows exponentially with system size, limiting the method to small systems, by representing the basis states on a quantum computer one may overcome this exponential scaling barrier. We propose a two-level multigrid Trotter time evolution for an efficient preparation of the basis states in a quantum circuit, which takes advantage of the robustness of the subspace expansion against Trotter errors. Using a quantum emulator we demonstrate the algorithm for the Hubbard model on a Bethe lattice with infinite coordination, which we map to a 16 qubit Anderson impurity model within the dynamical mean field theory. Our algorithm computes the Green's function accurately for both the metallic and Mott insulating regimes, with a circuit depth several orders of magnitude below what has been proposed using time evolution. The two-level multigrid time evolution reduces the number of Trotter steps required to compute the Green's function to about four to six. We therefore expect that the method can be used on near term quantum computers for moderate system sizes, while allowing for scalability to larger circuit depths and qubit numbers on future fault tolerant quantum computers.
We present the continued fraction method (CFM) as a microscopic approximation to the spectral density of the Hubbard model in the correlated metal phase away from half filling. The quantity expanded as a continued fraction is the single-particle Green's function. Leading spectral moments are taken into account through a set of real expansion coefficients, as known from the projection technique. Further stages are added to the continued fraction, with complex coefficients, thus defining a terminator function. This enables us to treat the entire spectral range of the Green's function on equal footing and determine the energy scale of the Fermi liquid quasiparticles by minimizing the total energy. The solution is free of phenomenological parameters and remains well defined in the strong-coupling limit, near the doping-controlled metal-insulator transition. Our results for the density of states agree reasonably with those from several variants of the dynamical mean-field theory. The CFM requires minimal numerical effort and can be generalized in several ways that are interesting for applications to real materials.
Current widely used approaches to calculate spectral functions using the density-matrix renormalization group in frequency space either necessarily include an artificial broadening (correction-vector method), have limited resolution (time-domain density-matrix renormalization group with Fourier transform method), or are limited to low-energy properties or single dominant modes (original continued fraction method). Here we propose an adaptive Lanczos-vector method to calculate the coefficients of a continued fraction expansion of the spectral function iteratively. We show that one can obtain a very accurate representation of the spectral function very efficiently, and that one can also directly extract the spectral weights and poles for the discrete system.