Description: Stochastic Analytic Continuation: A Bayesian Approach

This title appears in the Scientific Report : 2017

Stochastic Analytic Continuation: A Bayesian Approach

The stochastic sampling method (StochS) is used for the analytic continuation of quantum Monte Carlo data from the imaginary axis to the real axis. Compared to the maximum entropy method, StochS does not have explicit parameters, and one would expect the results to be unbiased. We present a very eff...

Personal Name(s):	Ghanem, Khaldoon (Corresponding author)
Contributing Institute:	Jülich Supercomputing Center; JSC
Imprint:	2017
Physical Description:	183 p.
Dissertation Note:	Dissertation, RWTH Aachen University, 2017
Document Type:	Dissertation / PhD Thesis
Research Program:	AICES Aachen Institute for Advanced Study in Computational Engineering Science Computational Science and Mathematical Methods
Link:	Get full text OpenAccess
	Publikationsportal JuSER

Please use the identifier: http://hdl.handle.net/2128/15995 in citations.

The stochastic sampling method (StochS) is used for the analytic continuation of quantum Monte Carlo data from the imaginary axis to the real axis. Compared to the maximum entropy method, StochS does not have explicit parameters, and one would expect the results to be unbiased. We present a very efficient algorithm for performing StochS and use it to study the effect of the discretization grid. Surprisingly, we find that the grid affects the results of StochS acting as an implicit default model. We provide a recipe for choosing a reliable StochS grid.To reduce the effect of the grid, we extend StochS into a gridless method (gStochS) by sampling the grid points from a default model instead of having them fixed. The effect of the default model is much reduced in gStochS compared to StochS and depends mainly on its width rather than its shape. The proper width can then be chosen using a simple recipe like we did in StochS.Finally, to avoid fixing the width, we go one step further and extend gStochS to sample over a whole class of default models with different widths. The extended method (eStochS) is then able to automatically relocate the grid points and concentrate them in the important region. Test cases show that eStochS gives good results resolving sharp features in the spectrum without the need for fine tuning a default model.