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

Physical review. B, Condensed matter

Quantum Monte Carlo simulations and maximum entropy: Dynamics from imaginary-time data

James E. Gubernatis, Mark Jarrell, Richard N. Silver, Devinderjit Singh Sivia

Published Sep 15, 1991• 389 citations
Publisher

Abstract

We report the details of an application of the method of maximum entropy to the extraction of spectral and transport properties from the imaginary-time correlation functions generated from quantum Monte Carlo simulations of the nondegenerate, symmetric, single-impurity Anderson model. We find that these physical properties are approximately universal functions of temperature and frequency when these parameters are scaled by the Kondo temperature. We also found that important details for successful extractions included the generation of statistically independent, Gaussian-distributed data, and a good choice of a default model to represent the state of our prior knowledge about the result in the absence of data. We suggest that our techniques are not restricted to the Hamiltonian and quantum Monte Carlo algorithm used here, but that maximum entropy and these techniques lay the general groundwork for the extraction of dynamical information from imaginary-time data generated by other quantum Monte Carlo simulations.

DOI
10.1103/physrevb.44.6011

Related papers

Dynamical properties from quantum Monte Carlo by the Maximum Entropy Method

Richard N. Silver, Devinderjit Singh Sivia, James E. Gubernatis

An outstanding problem in the simulation of condensed matter phenomena is how to obtain dynamical information. We consider the numerical analytic continuation of imaginary time Quantum Monte Carlo data to obtain real frequency spectral functions. We suggest an image reconstruction approach which has been widely applied to data analysis in experimental research, the Maximum Entropy Method (MaxEnt). We report encouraging preliminary results for the Fano-Anderson model of an impurity state in a continuum. The incorporation of additional prior information, such as sum rules and asymptotic behavior, can be expected to significantly improve results. We also compare MaxEnt to alternative methods. 17 refs., 4 figs.

•Jan 1, 1989•1 citation

Maximum-entropy method for analytic continuation of quantum Monte Carlo data

Richard N. Silver, Devinderjit Singh Sivia, James E. Gubernatis

An outstanding problem in the simulation of condensed-matter phenomena is how to obtain dynamical information. We consider the numerical analytic continuation of imaginary-time quantum Monte Carlo data to obtain real-frequency spectral functions. This is an extremely ill-posed problem similar to the inversion of a Laplace transform. We suggest an image-reconstruction approach, which has been widely applied to data analysis in experimental research. Specifically, we apply the maximum-entropy method (ME) to the analytic continuation of quantum Monte Carlo data. We report encouraging preliminary results for the Fano-Anderson model of an impurity state in a continuum. The incorporation of additional prior information, such as sum rules and asymptotic behavior, can be expected to significantly improve results. We compare (ME) to alternative methods. We also discuss statistical error propagation for the analytic continuation problem via the likelihood function, which is independent of the choice of image-reconstruction method. This includes the sensitivity of the data to structure in the spectral function, the optimization of Monte Carlo simulations, and how to incorporate covariance in the statistical errors of the Monte Carlo method.

Physical review. B, Condensed matter•Feb 1, 1990•306 citations

Continuous-time quantum Monte Carlo and maximum entropy approach to an imaginary-time formulation of strongly correlated steady-state transport

Andreas Dirks, Philipp Werner, Mark Jarrell +1 more

Recently, Han and Heary [Phys. Rev. Lett. 99, 236808 (2007)] proposed an approach to steady-state quantum transport through mesoscopic structures, which maps the nonequilibrium problem onto a family of auxiliary quantum impurity systems subject to imaginary voltages. We employ continuous-time quantum Monte-Carlo solvers to calculate accurate imaginary time data for the auxiliary models. The spectral function is obtained from a maximum entropy analytical continuation in both Matsubara frequency and complexified voltage. To enable the analytical continuation we construct a kernel which is compatible with the analytical structure of the theory. While it remains a formidable task to extract reliable spectral functions from this unbiased procedure, particularly for large voltages, our results indicate that the method in principle yields results in agreement with those obtained by other methods.

arXiv:1002.4081•Physical Review E•Aug 5, 2010•30 citations

Applications of Maxent to Quantum Monte Carlo

Richard N. Silver, Devinderjit Singh Sivia, James E. Gubernatis +1 more

We consider the application of maximum entropy methods to the analysis of data produced by computer simulations. The focus is the calculation of the dynamical properties of quantum many-body systems by Monte Carlo methods, which is termed the "Analytic Continuation Problem." For the Anderson model of dilute magnetic impurities in metals, we obtain spectral functions and transport coefficients which obey "Kondo Universality."

•Jan 1, 1991•0 citations

A Bayesian analysis of Green’s function Monte Carlo correlation functions

Michel Caffarel, David M. Ceperley

Imaginary-time correlation functions calculated by quantum Monte Carlo (QMC) are analyzed using the maximum entropy method (MaxEnt) to determine the ground-state energy and spectral overlap function. In contrast to earlier applications of MaxEnt, the data is obtained from importanced-sampled zero-temperature quantum Monte Carlo simulations. The analysis includes two steps. First, that spectral overlap function and ground state energy which maximizes the entropy and agrees with the QMC correlation functions is obtained. Then the errors in the energy are evaluated by averaging over all the possible images (average MaxEnt method), the multidimensional integrals being computed using the Metropolis algorithm. The central feature of this approach is that all the information present in the correlation functions is used in the only way consistent with fundamental probabilistic hypotheses. This allows us to fully exploit the information contained in the correlation functions at small imaginary times, thus avoiding large statistical fluctuations associated with large imaginary times. In addition, the computed errors include both the statistical errors and systematic extrapolation errors. The method is illustrated with a harmonic oscillator and the four-electron LiH molecule.

The Journal of Chemical Physics•Dec 1, 1992•43 citations

Bayesian methods, maximum entropy, and quantum Monte Carlo

James E. Gubernatis, Richard N. Silver, Mark Jarrell

We heuristically discuss the application of the method of maximum entropy to the extraction of dynamical information from imaginary-time, quantum Monte Carlo data. The discussion emphasizes the utility of a Bayesian approach to statistical inference and the importance of statistically well-characterized data. 14 refs.

•Jan 1, 1991•0 citations
We report the details of an application of the method of maximum entropy to the extraction of spectral and transport properties from the imaginary-time correlation functions generated from quantum Monte Carlo simulations of the nondegenerate, symmetric, single-impurity Anderson model. We find that these physical properties are approximately universal functions of temperature and frequency when these parameters are scaled by the Kondo temperature. We also found that important details for successful extractions included the generation of statistically independent, Gaussian-distributed data, and a good choice of a default model to represent the state of our prior knowledge about the result in the absence of data. We suggest that our techniques are not restricted to the Hamiltonian and quantum Monte Carlo algorithm used here, but that maximum entropy and these techniques lay the general groundwork for the extraction of dynamical information from imaginary-time data generated by other quantum Monte Carlo simulations.
An outstanding problem in the simulation of condensed matter phenomena is how to obtain dynamical information. We consider the numerical analytic continuation of imaginary time Quantum Monte Carlo data to obtain real frequency spectral functions. We suggest an image reconstruction approach which has been widely applied to data analysis in experimental research, the Maximum Entropy Method (MaxEnt). We report encouraging preliminary results for the Fano-Anderson model of an impurity state in a continuum. The incorporation of additional prior information, such as sum rules and asymptotic behavior, can be expected to significantly improve results. We also compare MaxEnt to alternative methods. 17 refs., 4 figs.
An outstanding problem in the simulation of condensed-matter phenomena is how to obtain dynamical information. We consider the numerical analytic continuation of imaginary-time quantum Monte Carlo data to obtain real-frequency spectral functions. This is an extremely ill-posed problem similar to the inversion of a Laplace transform. We suggest an image-reconstruction approach, which has been widely applied to data analysis in experimental research. Specifically, we apply the maximum-entropy method (ME) to the analytic continuation of quantum Monte Carlo data. We report encouraging preliminary results for the Fano-Anderson model of an impurity state in a continuum. The incorporation of additional prior information, such as sum rules and asymptotic behavior, can be expected to significantly improve results. We compare (ME) to alternative methods. We also discuss statistical error propagation for the analytic continuation problem via the likelihood function, which is independent of the choice of image-reconstruction method. This includes the sensitivity of the data to structure in the spectral function, the optimization of Monte Carlo simulations, and how to incorporate covariance in the statistical errors of the Monte Carlo method.
Recently, Han and Heary [Phys. Rev. Lett. 99, 236808 (2007)] proposed an approach to steady-state quantum transport through mesoscopic structures, which maps the nonequilibrium problem onto a family of auxiliary quantum impurity systems subject to imaginary voltages. We employ continuous-time quantum Monte-Carlo solvers to calculate accurate imaginary time data for the auxiliary models. The spectral function is obtained from a maximum entropy analytical continuation in both Matsubara frequency and complexified voltage. To enable the analytical continuation we construct a kernel which is compatible with the analytical structure of the theory. While it remains a formidable task to extract reliable spectral functions from this unbiased procedure, particularly for large voltages, our results indicate that the method in principle yields results in agreement with those obtained by other methods.
Imaginary-time correlation functions calculated by quantum Monte Carlo (QMC) are analyzed using the maximum entropy method (MaxEnt) to determine the ground-state energy and spectral overlap function. In contrast to earlier applications of MaxEnt, the data is obtained from importanced-sampled zero-temperature quantum Monte Carlo simulations. The analysis includes two steps. First, that spectral overlap function and ground state energy which maximizes the entropy and agrees with the QMC correlation functions is obtained. Then the errors in the energy are evaluated by averaging over all the possible images (average MaxEnt method), the multidimensional integrals being computed using the Metropolis algorithm. The central feature of this approach is that all the information present in the correlation functions is used in the only way consistent with fundamental probabilistic hypotheses. This allows us to fully exploit the information contained in the correlation functions at small imaginary times, thus avoiding large statistical fluctuations associated with large imaginary times. In addition, the computed errors include both the statistical errors and systematic extrapolation errors. The method is illustrated with a harmonic oscillator and the four-electron LiH molecule.