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:1401.0526 · Physical Review B

Green's functions from real-time bold-line Monte Carlo

Guy Cohen, David R. Reichman, Andrew J. Millis, Emanuel Gull

Published Mar 31, 2014• 77 citations
arXivHTMLPDF

Abstract

We present two methods for computing two-time correlation functions or Green's functions from real-time bold-line continuous-time quantum Monte Carlo. One method is a formally exact generalized auxiliary lead formalism by which spectral properties may be obtained from single-time observables. The other involves the evaluation of diagrams contributing to two-time observables directly on the Keldysh contour. Additionally, we provide a detailed description of the bold-line Monte Carlo method. Our methods are general and numerically exact, and able to reliably resolve high-energy features such as band edges. We compare the spectral functions obtained from real-time methods to analytically continued spectral functions obtained from imaginary-time Monte Carlo, thus probing the limits of analytic continuation.

DOI
10.1103/physrevb.89.115139

Related papers

Green’s Functions from Real-Time Bold-Line Monte Carlo Calculations: Spectral Properties of the Nonequilibrium Anderson Impurity Model

Guy Cohen, Emanuel Gull, David R. Reichman +1 more

The nonequilibrium spectral properties of the Anderson impurity model with a chemical potential bias are investigated within a numerically exact real-time quantum Monte Carlo formalism. The two-time correlation function is computed in a form suitable for nonequilibrium dynamical mean field calculations. Additionally, the evolution of the model's spectral properties are simulated in an alternative representation, defined by a hypothetical but experimentally realizable weakly coupled auxiliary lead. The voltage splitting of the Kondo peak is confirmed and the dynamics of its formation after a coupling or gate quench are studied. This representation is shown to contain additional information about the dot's population dynamics. Further, we show that the voltage-dependent differential conductance gives a reasonable qualitative estimate of the equilibrium spectral function, but significant qualitative differences are found including incorrect trends and spurious temperature dependent effects.

arXiv:1310.4151•Physical Review Letters•Apr 9, 2014•108 citations

“Correlation Between Ct Qmc- (Continuous Time Quantum Monte Carlo) And Green's Function Theory”.

Piyali J. Chakravorty

Recent years have seen a rise in the popularity of the physics of strongly interacting electronic systems, and in most of the more interesting cases the perturbative regime is not applicable. The continuous-time quantum Monte Carlo (CT-QMC) has already established its leading stochastic simulation method to address quantum impurity models that have gained central focus of DMFT. Lugiato-Lefever function approach, on the one hand, can provide a rigorous account of the electronic structure, excitation spectrum and phenomenon of correlation to most of the many-body systems. In the paper the connection between the CT-QMC technique and theory of Green functions are stated with a specification in the connection of the two, complementing each other in the solving of quantum many-body problems. We present the approach of CT-QMC that is numerically inexpensive but computes Green functions directly in the imaginary time that you then analytically continue to the real frequencies to obtain physical observables. Further, we discuss the ability to investigate quasiparticle behaviour, spectral density and transport phenomenon through self-energies acquired as Green functions in CT-QMC.

Advances in Nonlinear Variational Inequalities•Jan 24, 2025•0 citations

Extracting spectral properties from Keldysh Green functions

Andreas Dirks, Martin Eckstein, Thomas Pruschke +1 more

We investigate the possibility to assist the numerically ill-posed calculation of spectral properties of interacting quantum systems in thermal equilibrium by extending the imaginary-time simulation to a finite Schwinger-Keldysh contour. The effect of this extension is tested within the standard maximum entropy approach to analytic continuation. We find that the inclusion of real-time data improves the resolution of structures at high energy, while the imaginary-time data are needed to correctly reproduce low-frequency features such as quasiparticle peaks. As a nonequilibrium application, we consider the calculation of time-dependent spectral functions from retarded Green function data on a finite time interval, and compare the maximum entropy approach to direct Fourier transformation and a method based on Padé approximants.

arXiv:1210.7867•Physical Review E•Feb 21, 2013•12 citations

Flat histogram quantum Monte Carlo for analytic continuation to real time

Diamantis, Nikolaos G., Efstratios Manousakis

The Quantum Monte Carlo (QMC) method can yield the imaginary-time dependence of a correlation function $C(τ)$ of an operator $\hat O$. The analytic continuation to real-time proceeds by means of a "numerical inversion" of these data to find the response function or spectral density $A(ω)$ corresponding to $\hat O$. Such a technique is very sensitive to the statistical errors in $C(τ)$ especially for large values of $τ$, when we are interested in the low-energy excitations. In this paper, we find that if we use the flat histogram technique in the QMC method, in such a way to make the {\it histogram of} $C(τ)$ flat, the results of the analytic continuation for low-energy excitations improve using the same amount of computational time. To demonstrate the idea we select an exactly soluble version of the single-hole motion in the $t-J$ model and the diagrammatic Monte Carlo technique.

arXiv:1411.5984•Nov 21, 2014•0 citations

Numerically exact Green functions from Hirsch-Fye quantum Monte Carlo simulations

Nils Blümer

We present a new method for extracting numerically exact imaginary-time Green functions from standard Hirsch-Fye quantum Monte Carlo (HF-QMC) simulations within dynamical mean-field theory. By analytic continuation, angular resolved spectra are obtained without the discretization bias previously associated with HF-QMC results. The method is shown to be accurate even at very low temperatures (T=W/800 for bandwidth W) in the strongly correlated regime.

arXiv:0712.1290•Dec 8, 2007•1 citation

Matsubara Green's Function and its Analytic Continuation using Flat Histogram Quantum Monte Carlo Simulation

Nikolaos G. Diamantis, Efstratios Manousakis

We demonstrate that by incorporating flat histogram techniques in the quantum Monte Carlo simulation, we can calculate the Matsubara Green's function for significantly larger values of the imaginary time. This allows us to obtain more accurate results for the spectral function, particularly for low energy excitations, using known techniques for analytic continuation to real time.

Physics Procedia•Jan 1, 2015•2 citations
We present two methods for computing two-time correlation functions or Green's functions from real-time bold-line continuous-time quantum Monte Carlo. One method is a formally exact generalized auxiliary lead formalism by which spectral properties may be obtained from single-time observables. The other involves the evaluation of diagrams contributing to two-time observables directly on the Keldysh contour. Additionally, we provide a detailed description of the bold-line Monte Carlo method. Our methods are general and numerically exact, and able to reliably resolve high-energy features such as band edges. We compare the spectral functions obtained from real-time methods to analytically continued spectral functions obtained from imaginary-time Monte Carlo, thus probing the limits of analytic continuation.
The nonequilibrium spectral properties of the Anderson impurity model with a chemical potential bias are investigated within a numerically exact real-time quantum Monte Carlo formalism. The two-time correlation function is computed in a form suitable for nonequilibrium dynamical mean field calculations. Additionally, the evolution of the model's spectral properties are simulated in an alternative representation, defined by a hypothetical but experimentally realizable weakly coupled auxiliary lead. The voltage splitting of the Kondo peak is confirmed and the dynamics of its formation after a coupling or gate quench are studied. This representation is shown to contain additional information about the dot's population dynamics. Further, we show that the voltage-dependent differential conductance gives a reasonable qualitative estimate of the equilibrium spectral function, but significant qualitative differences are found including incorrect trends and spurious temperature dependent effects.
Recent years have seen a rise in the popularity of the physics of strongly interacting electronic systems, and in most of the more interesting cases the perturbative regime is not applicable. The continuous-time quantum Monte Carlo (CT-QMC) has already established its leading stochastic simulation method to address quantum impurity models that have gained central focus of DMFT. Lugiato-Lefever function approach, on the one hand, can provide a rigorous account of the electronic structure, excitation spectrum and phenomenon of correlation to most of the many-body systems. In the paper the connection between the CT-QMC technique and theory of Green functions are stated with a specification in the connection of the two, complementing each other in the solving of quantum many-body problems. We present the approach of CT-QMC that is numerically inexpensive but computes Green functions directly in the imaginary time that you then analytically continue to the real frequencies to obtain physical observables. Further, we discuss the ability to investigate quasiparticle behaviour, spectral density and transport phenomenon through self-energies acquired as Green functions in CT-QMC.
We investigate the possibility to assist the numerically ill-posed calculation of spectral properties of interacting quantum systems in thermal equilibrium by extending the imaginary-time simulation to a finite Schwinger-Keldysh contour. The effect of this extension is tested within the standard maximum entropy approach to analytic continuation. We find that the inclusion of real-time data improves the resolution of structures at high energy, while the imaginary-time data are needed to correctly reproduce low-frequency features such as quasiparticle peaks. As a nonequilibrium application, we consider the calculation of time-dependent spectral functions from retarded Green function data on a finite time interval, and compare the maximum entropy approach to direct Fourier transformation and a method based on Padé approximants.
The Quantum Monte Carlo (QMC) method can yield the imaginary-time dependence of a correlation function C(τ)C(τ)C(τ) of an operator O^\hat OO^. The analytic continuation to real-time proceeds by means of a "numerical inversion" of these data to find the response function or spectral density A(ω)A(ω)A(ω) corresponding to O^\hat OO^. Such a technique is very sensitive to the statistical errors in C(τ)C(τ)C(τ) especially for large values of τττ, when we are interested in the low-energy excitations. In this paper, we find that if we use the flat histogram technique in the QMC method, in such a way to make the {\it histogram of} C(τ)C(τ)C(τ) flat, the results of the analytic continuation for low-energy excitations improve using the same amount of computational time. To demonstrate the idea we select an exactly soluble version of the single-hole motion in the t−Jt-Jt−J model and the diagrammatic Monte Carlo technique.