skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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 » for data of poor quality.« less

Authors:
 [1];  [1];  [1];  [1];  [2];  [1];  [1]
  1. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
  2. 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}
}

Journal Article:

Citation Metrics:
Cited by: 12 works
Citation information provided by
Web of Science

Save / Share:

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


Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions
journal, January 1996


Spectral densities of the symmetric Anderson model
journal, July 1990


Quantum Monte Carlo simulations and maximum entropy: Dynamics from imaginary-time data
journal, September 1991


Maximum-entropy method for analytic continuation of quantum Monte Carlo data
journal, February 1990


Self-consistent large- N expansion for normal-state properties of dilute magnetic alloys
journal, August 1987


Spectral analysis by the method of consistent constraints
journal, August 2013


Continuous-time auxiliary-field Monte Carlo for quantum impurity models
journal, May 2008


Quantum cluster theories
journal, October 2005


The pseudogap: friend or foe of high T c ?
journal, December 2005


Analytic continuation of quantum Monte Carlo data by stochastic analytical inference
journal, May 2010


Bayesian inference and the analytic continuation of imaginary-time quantum Monte Carlo data
journal, May 1996


Stochastic method for analytic continuation of quantum Monte Carlo data
journal, May 1998


Phase diagram of the Hubbard model: Beyond the dynamical mean field
journal, November 2001


Diagrammatic quantum Monte Carlo study of the Fröhlich polaron
journal, September 2000


Markov Chain Monte Carlo
book, January 2005


Maximum entropy analysis of oversampled data problems
journal, April 1990


Works referencing / citing this record:

Numerical analytic continuation: Answers to well-posed questions
text, January 2016


Bayesian parametric analytic continuation of Green's functions
text, January 2019