Fast and Efficient Stochastic Optimization for Analytic Continuation
Abstract
In this analytic continuation of imaginary-time quantum Monte Carlo data to extract real-frequency spectra remains a key problem in connecting theory with experiment. Here we present a fast and efficient stochastic optimization method (FESOM) as a more accessible variant of the stochastic optimization method introduced by Mishchenko et al. [Phys. Rev. B 62, 6317 (2000)], and we benchmark the resulting spectra with those obtained by the standard maximum entropy method for three representative test cases, including data taken from studies of the two-dimensional Hubbard model. Genearally, we find that our FESOM approach yields spectra similar to the maximum entropy results. In particular, while the maximum entropy method yields superior results when the quality of the data is strong, we find that FESOM is able to resolve fine structure with more detail when the quality of the data is poor. In addition, because of its stochastic nature, the method provides detailed information on the frequency-dependent uncertainty of the resulting spectra, while the maximum entropy method does so only for the spectral weight integrated over a finite frequency region. Therefore, we believe that this variant of the stochastic optimization approach provides a viable alternative to the routinely used maximum entropy method, especiallymore »
- Authors:
-
- Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
- Virginia Polytechnic Inst. and State Univ. (Virginia Tech), Blacksburg, VA (United States)
- Publication Date:
- Research Org.:
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
- Sponsoring Org.:
- USDOE Laboratory Directed Research and Development (LDRD) Program
- OSTI Identifier:
- 1336574
- Alternate Identifier(s):
- OSTI ID: 1327072
- Grant/Contract Number:
- AC05-00OR22725
- Resource Type:
- Journal Article: Accepted Manuscript
- Journal Name:
- Physical Review B
- Additional Journal Information:
- Journal Volume: 94; Journal Issue: 12; Journal ID: ISSN 2469-9950
- Publisher:
- American Physical Society (APS)
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS
Citation Formats
Bao, Feng, Zhang, Guannan, Webster, Clayton G, Tang, Yanfei, Scarola, Vito, Summers, Michael Stuart, and Maier, Thomas A. Fast and Efficient Stochastic Optimization for Analytic Continuation. United States: N. p., 2016.
Web. doi:10.1103/PhysRevB.94.125149.
Bao, Feng, Zhang, Guannan, Webster, Clayton G, Tang, Yanfei, Scarola, Vito, Summers, Michael Stuart, & Maier, Thomas A. Fast and Efficient Stochastic Optimization for Analytic Continuation. United States. https://doi.org/10.1103/PhysRevB.94.125149
Bao, Feng, Zhang, Guannan, Webster, Clayton G, Tang, Yanfei, Scarola, Vito, Summers, Michael Stuart, and Maier, Thomas A. 2016.
"Fast and Efficient Stochastic Optimization for Analytic Continuation". United States. https://doi.org/10.1103/PhysRevB.94.125149. https://www.osti.gov/servlets/purl/1336574.
@article{osti_1336574,
title = {Fast and Efficient Stochastic Optimization for Analytic Continuation},
author = {Bao, Feng and Zhang, Guannan and Webster, Clayton G and Tang, Yanfei and Scarola, Vito and Summers, Michael Stuart and Maier, Thomas A},
abstractNote = {In this analytic continuation of imaginary-time quantum Monte Carlo data to extract real-frequency spectra remains a key problem in connecting theory with experiment. Here we present a fast and efficient stochastic optimization method (FESOM) as a more accessible variant of the stochastic optimization method introduced by Mishchenko et al. [Phys. Rev. B 62, 6317 (2000)], and we benchmark the resulting spectra with those obtained by the standard maximum entropy method for three representative test cases, including data taken from studies of the two-dimensional Hubbard model. Genearally, we find that our FESOM approach yields spectra similar to the maximum entropy results. In particular, while the maximum entropy method yields superior results when the quality of the data is strong, we find that FESOM is able to resolve fine structure with more detail when the quality of the data is poor. In addition, because of its stochastic nature, the method provides detailed information on the frequency-dependent uncertainty of the resulting spectra, while the maximum entropy method does so only for the spectral weight integrated over a finite frequency region. Therefore, we believe that this variant of the stochastic optimization approach provides a viable alternative to the routinely used maximum entropy method, especially for data of poor quality.},
doi = {10.1103/PhysRevB.94.125149},
url = {https://www.osti.gov/biblio/1336574},
journal = {Physical Review B},
issn = {2469-9950},
number = 12,
volume = 94,
place = {United States},
year = {Wed Sep 28 00:00:00 EDT 2016},
month = {Wed Sep 28 00:00:00 EDT 2016}
}
Web of Science
Works referenced in this record:
Dynamics of a hole in the t - J model with local disorder: Exact results for high dimensions
journal, December 1992
- Strack, Rainer; Vollhardt, Dieter
- Physical Review B, Vol. 46, Issue 21
Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions
journal, January 1996
- Georges, Antoine; Kotliar, Gabriel; Krauth, Werner
- Reviews of Modern Physics, Vol. 68, Issue 1
Spectral densities of the symmetric Anderson model
journal, July 1990
- Silver, R. N.; Gubernatis, J. E.; Sivia, D. S.
- Physical Review Letters, Vol. 65, Issue 4
Quantum Monte Carlo simulations and maximum entropy: Dynamics from imaginary-time data
journal, September 1991
- Gubernatis, J. E.; Jarrell, Mark; Silver, R. N.
- Physical Review B, Vol. 44, Issue 12
Maximum-entropy method for analytic continuation of quantum Monte Carlo data
journal, February 1990
- Silver, R. N.; Sivia, D. S.; Gubernatis, J. E.
- Physical Review B, Vol. 41, Issue 4
Self-consistent large- expansion for normal-state properties of dilute magnetic alloys
journal, August 1987
- Bickers, N. E.; Cox, D. L.; Wilkins, J. W.
- Physical Review B, Vol. 36, Issue 4
Spectral analysis by the method of consistent constraints
journal, August 2013
- Prokof’ev, N. V.; Svistunov, B. V.
- JETP Letters, Vol. 97, Issue 11
Continuous-time auxiliary-field Monte Carlo for quantum impurity models
journal, May 2008
- Gull, E.; Werner, P.; Parcollet, O.
- EPL (Europhysics Letters), Vol. 82, Issue 5
Quantum cluster theories
journal, October 2005
- Maier, Thomas; Jarrell, Mark; Pruschke, Thomas
- Reviews of Modern Physics, Vol. 77, Issue 3
The pseudogap: friend or foe of high T c ?
journal, December 2005
- Norman, M. R.; Pines, D.; Kallin, C.
- Advances in Physics, Vol. 54, Issue 8
Analytic continuation of quantum Monte Carlo data by stochastic analytical inference
journal, May 2010
- Fuchs, Sebastian; Pruschke, Thomas; Jarrell, Mark
- Physical Review E, Vol. 81, Issue 5
Bayesian inference and the analytic continuation of imaginary-time quantum Monte Carlo data
journal, May 1996
- Jarrell, Mark; Gubernatis, J. E.
- Physics Reports, Vol. 269, Issue 3
Stochastic method for analytic continuation of quantum Monte Carlo data
journal, May 1998
- Sandvik, Anders W.
- Physical Review B, Vol. 57, Issue 17
Phase diagram of the Hubbard model: Beyond the dynamical mean field
journal, November 2001
- Jarrell, M.; Maier, Th; Hettler, M. H.
- Europhysics Letters (EPL), Vol. 56, Issue 4
Diagrammatic quantum Monte Carlo study of the Fröhlich polaron
journal, September 2000
- Mishchenko, A.; Prokof’ev, N.; Sakamoto, A.
- Physical Review B, Vol. 62, Issue 10
Maximum entropy analysis of oversampled data problems
journal, April 1990
- Bryan, R. K.
- European Biophysics Journal, Vol. 18, Issue 3
Works referencing / citing this record:
Numerical analytic continuation: Answers to well-posed questions
text, January 2016
- Goulko, Olga; Mishchenko, Andrey S.; Pollet, Lode
- arXiv
Bayesian parametric analytic continuation of Green's functions
text, January 2019
- Rumetshofer, Michael; Bauernfeind, Daniel; von der Linden, Wolfgang
- arXiv