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

Physics Procedia

Application of Flat Histogram Methods to Extract the Long Imaginary-time Behavior in Diagrammatic Monte Carlo

Nikolaos G. Diamantis, Efstratios Manousakis

Published Jan 1, 2014• 0 citations
Publisher

Abstract

We demonstrate that using flat histogram methods we can extract the imaginary-time behavior of the Green's function G from diagrammatic Monte Carlo simulations very accurately even when G changes by many orders of magnitude.
DOI
10.1016/j.phpro.2014.08.130

Related papers

Flat histogram diagrammatic Monte Carlo method: Calculation of the Green's function in imaginary time

Nikolaos G. Diamantis, Efstratios Manousakis

The diagrammatic Monte Carlo (DiagMC) method is a numerical technique which samples the entire diagrammatic series of the Green's function in quantum many-body systems. In this work, we incorporate the flat histogram principle in the diagrammatic Monte Carlo method, and we term the improved version the "flat histogram diagrammatic Monte Carlo" method. We demonstrate the superiority of this method over the standard DiagMC in extracting the long-imaginary-time behavior of the Green's function, without incorporating any a priori knowledge about this function, by applying the technique to the polaron problem.

arXiv:1306.6320•Physical Review E•Oct 9, 2013•5 citations

Flat-histogram methods in quantum Monte Carlo simulations: Application to the t-J model

Nikolaos G. Diamantis, Efstratios Manousakis

We discuss that flat-histogram techniques can be appropriately applied in the sampling of quantum Monte Carlo simulation in order to improve the statistical quality of the results at long imaginary time or low excitation energy. Typical imaginary-time correlation functions calculated in quantum Monte Carlo are subject to exponentially growing errors as the range of imaginary time grows and this smears the information on the low energy excitations. We show that we can extract the low energy physics by modifying the Monte Carlo sampling technique to one in which configurations which contribute to making the histogram of certain quantities flat are promoted. We apply the diagrammatic Monte Carlo (diag-MC) method to the motion of a single hole in the t-J model and we show that the implementation of flat-histogram techniques allows us to calculate the Green's function in a wide range of imaginary-time. In addition, we show that applying the flat-histogram technique alleviates the "sign"-problem associated with the simulation of the single-hole Green's function at long imaginary time.

Journal of Physics Conference Series•Sep 1, 2016•0 citations

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

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

Analytical solution for time integrals in diagrammatic expansions: Application to real-frequency diagrammatic Monte Carlo

J. Vučičević, P. Stipsić, Michel Ferrero

The authors find an analytical solution to the imaginary time integrals that arise in quantum many-body perturbation theory, and use it to investigate the two-dimensional Hubbard model with diagrammatic Monte Carlo algorithms.

arXiv:2011.08226•Physical Review Research•Apr 29, 2021•23 citations

Real-time dynamics from imaginary-time quantum Monte Carlo simulations: Tests on oscillator chains

J. Bonča, James E. Gubernatis

We used methods of Bayesian statistical inference and the principle of maximum entropy to analytically continue imaginary-time Green's functions generated in quantum Monte Carlo simulations to obtain the real-time Green's functions. For test problems, we considered chains of harmonic and anharmonic oscillators whose properties we simulated by a hybrid path-integral quantum Monte Carlo method. From the imaginary-time displacement-displacement Green's function, we first obtained its spectral density. For harmonic oscillators, we demonstrated the peaks of this function were in the correct position and their areas satisfied a sum rule. Additionally, as a function of wave number, the peak positions followed the correct dispersion relation. For a double-well oscillator, we demonstrated that the peak location correctly predicted the tunnel splitting. Transforming the spectral densities to real-time Green's functions, we conclude that we can predict the real-time dynamics for length of times corresponding to five to ten times the natural period of the model. The length of time was limited by an overbroadening of the peaks in the spectral density caused by the simulation algorithm. \textcopyright{} 1996 The American Physical Society.

arXiv:cond-mat/9510098•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jun 1, 1996•26 citations
The diagrammatic Monte Carlo (DiagMC) method is a numerical technique which samples the entire diagrammatic series of the Green's function in quantum many-body systems. In this work, we incorporate the flat histogram principle in the diagrammatic Monte Carlo method, and we term the improved version the "flat histogram diagrammatic Monte Carlo" method. We demonstrate the superiority of this method over the standard DiagMC in extracting the long-imaginary-time behavior of the Green's function, without incorporating any a priori knowledge about this function, by applying the technique to the polaron problem.
We discuss that flat-histogram techniques can be appropriately applied in the sampling of quantum Monte Carlo simulation in order to improve the statistical quality of the results at long imaginary time or low excitation energy. Typical imaginary-time correlation functions calculated in quantum Monte Carlo are subject to exponentially growing errors as the range of imaginary time grows and this smears the information on the low energy excitations. We show that we can extract the low energy physics by modifying the Monte Carlo sampling technique to one in which configurations which contribute to making the histogram of certain quantities flat are promoted. We apply the diagrammatic Monte Carlo (diag-MC) method to the motion of a single hole in the t-J model and we show that the implementation of flat-histogram techniques allows us to calculate the Green's function in a wide range of imaginary-time. In addition, we show that applying the flat-histogram technique alleviates the "sign"-problem associated with the simulation of the single-hole Green's function at long imaginary time.
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.
We used methods of Bayesian statistical inference and the principle of maximum entropy to analytically continue imaginary-time Green's functions generated in quantum Monte Carlo simulations to obtain the real-time Green's functions. For test problems, we considered chains of harmonic and anharmonic oscillators whose properties we simulated by a hybrid path-integral quantum Monte Carlo method. From the imaginary-time displacement-displacement Green's function, we first obtained its spectral density. For harmonic oscillators, we demonstrated the peaks of this function were in the correct position and their areas satisfied a sum rule. Additionally, as a function of wave number, the peak positions followed the correct dispersion relation. For a double-well oscillator, we demonstrated that the peak location correctly predicted the tunnel splitting. Transforming the spectral densities to real-time Green's functions, we conclude that we can predict the real-time dynamics for length of times corresponding to five to ten times the natural period of the model. The length of time was limited by an overbroadening of the peaks in the spectral density caused by the simulation algorithm. \textcopyright{} 1996 The American Physical Society.