Stochastic Optimization Method for Analytic Continuation

The TRIQS-based Stochastic Optimization Method for Analytic Continuation allows to solve a family of Fredholm integral equations of the first kind. Numerical solution of such equations is required to reconstruct the spectral representation of physical observables (Green’s functions, dynamical susceptibilities) from noisy results of Quantum Monte Carlo simulations.

The integral equations have the following general form,

\[\int_{\epsilon_\mathrm{min}}^{\epsilon_\mathrm{max}} K(m,\epsilon) A(\epsilon) d\epsilon = G(m),\]

where \(G(m)\) is an input quantity, known with some uncertainty, and \(A(\epsilon)\) is a spectral function to be found. Meaning of the index \(m\), limits of the integral, and the form of the kernel \(K(m,\epsilon)\) depend on the problem at hand. \(A(\epsilon)\) is assumed to be non-negative and normalized to a known constant. A full list of supported observables with respective integral kernels is given here.

The family of Fredholm integral equations of the first kind shares a common property of being ill-posed problems; tiny variations of the right hand part \(S(m)\) may result in huge changes of the solution.

The Stochastic Optimization Method uses Markov chain sampling to minimize an objective function

\[\mathcal{D} = \sum_m \frac{1}{S(m)}\left| \int_{\epsilon_\mathrm{min}}^{\epsilon_\mathrm{max}} K(m,\epsilon) A(\epsilon) d\epsilon - G(m) \right|.\]

One can use a weight factor \(S(m)\) to stress importance of some data points \(G(m)\).

An important feature of the algorithm is a “continuous-energy” representation of the solutions. Function \(A(\epsilon)\) is parameterized as a sum of rectangles with certain positions, widths and heights.

Despite a number of optimizations, this code is considerably more CPU-intensive than algorithms based on the Bayesian statistical inference (MaxEnt) or Padé approximants. Its primary design goal is the ability to resolve fine spectral features, rather than giving a quick insight into general shape of the spectrum.