arXiv:1001.4351 · Physical Review B

Analytical continuation of imaginary axis data using maximum entropy

O. Gunnarsson, M. W. Haverkort, Giorgio Sangiovanni

Published Apr 9, 2010• 47 citations

Abstract

We study the maximum entropy (MaxEnt) approach for analytical continuation of spectral data from imaginary times to real frequencies. The total error is divided in a statistical error, due to the noise in the input data, and a systematic error, due to deviations of the default function, used in the MaxEnt approach, from the exact spectrum. We find that the MaxEnt approach in its classical formulation can lead to a nonoptimal balance between the two types of errors, leading to an unnecessary large statistical error. The statistical error can be reduced by splitting up the data in several batches, performing a MaxEnt calculation for each batch and averaging. This can outweigh an increase in the systematic error resulting from this approach. The output from the MaxEnt calculation can be used as a default function for a new MaxEnt calculation. Such iterations often lead to worse results due to an increase in the statistical error. By splitting up the data in batches, the statistical error is reduced and the increase resulting from iterations can be outweighed by a decrease in the systematic error. Finally we consider a linearized version to obtain a better understanding of the method.

DOI: 10.1103/physrevb.81.155107

Related papers

Analytical continuation of imaginary axis data for optical conductivity

O. Gunnarsson, M. W. Haverkort, Giorgio Sangiovanni

We compare different methods for performing analytical continuation of spectral data from the imaginary time or frequency axis to the real frequency axis for the optical conductivity $\ensuremath{\sigma}(\ensuremath{\omega})$. We compare the maximum entropy (MaxEnt), singular value decomposition (SVD), sampling, and Pad\'e methods for analytical continuation. We also study two direct methods for obtaining $\ensuremath{\sigma}(0)$. For the MaxEnt approach we focus on a recent modification. The data are split up in batches, a separate MaxEnt calculation is done for each batch and the results are averaged. For the problems studied here, we find that typically the SVD, sampling, and modified MaxEnt methods give comparable accuracy while the Pad\'e approximation is usually less reliable.

arXiv:1012.5934Physical Review BOct 28, 201066 citations

Analytical continuation of spectral data from imaginary time axis to real frequency axis using statistical sampling

Kiamars Vafayi, O. Gunnarsson

We present a method for performing analytical continuation of spectral data from imaginary time to real frequencies based on a statistical sampling method. Compared with the maximum entropy method (MEM), an advantage is that no default model needs to be introduced. For the problems studied here, the statistical sampling method gives comparable or slightly better results than MEM using quite accurate default models.

Physical Review BJul 19, 200746 citations

Analytic continuation of Green's functions with a neural network

Fakher F. Assaad, Johanna Erdmenger B, Anika Götz +3 more

An important problem in many-body physics is to reconstruct the spectral density from the imaginary-time domain Green's function. Typically, the imaginary-time Green's function is generated by Monte Carlo methods. As the one-point fermionic kernel diverges exponentially for large frequencies, numerical noise generically causes instabilities. We use a convolutional neural network to obtain the spectral density for a given imaginary time Green's function. The network is trained by data which we generate using random Gaussians. We improve the training data set available by including collision centers for the Gaussians rather than employing uniformly distributed Gaussians. Our network is constructed in such a way that its output fulfills positive semidefiniteness. We compare the results of our network with results of the Maximum Entropy method (MaxEnt), a standard method for the same reconstruction problem for the spectral density. This comparison is performed for three different cases, namely our Gaussian based test data as well as two physical models, the 1d Hubbard model showing spin-charge separation, and the two-dimensional SSH model in the self-consistent Born approximation. We find that the network outperforms MaxEnt when presented data close to the training set. For the physical models considered, MaxEnt recognizes physical features more precisely as compared to our network prediction. While it is hard to improve MaxEnt, the quality of the network depends on the training data set which can be systematically enhanced and improved.

arXiv:2602.18561Feb 20, 20260 citations

Suppression of Overfitting in Extraction of Spectral Data from Imaginary Frequency Green Function Using Maximum Entropy Method

Enzhi Li

Although maximum entropy method (maxEnt method) is currently the standard algorithm for extracting real frequency information from imaginary frequency Green function, still this method is beset with overfitting problem, which manifests itself as the spurious spikes in the resultant spectral functions. To address this issue and motivated by the regularization techniques widely used in machine learning and statistics, here we propose to add one more regularization term into the original maxEnt loss function to suppress these redundant spikes. The essence of this extra regularization term is to demand that the resultant spectral functions should pay a price for being spiky. We test our algorithm with both artificial and real data, and find that spurious spikes in the resultant spectral functions can be effectively suppressed by this method.

arXiv:1810.06147Oct 15, 20180 citations

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 matterFeb 1, 1990306 citations

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, 19891 citation