Information about the pairing mechanism for superconductivity is contained in the spectral weight for the anomalous (Gorkov) Green's function. In the most general case, this spectral weight can change sign on the positive real axis or even be complex in the presence of broken time-reversal symmetry. This renders impossible the direct analytic continuation with maximum entropy methods of numerical results obtained for the anomalous Green's function either in Matsubara frequency or in imaginary time. Here we show that one can define auxiliary Green's functions that allow one to extract the anomalous spectral weight in the most general multiorbital spin-singlet or triplet cases with no particular symmetry. As a simple example, we treat the case of lead in Eliashberg theory: Maximum entropy analytic continuation with auxiliary Green's functions (MaxEntAux) allows us to recover spectral functions that change sign on the positive real axis. The approach can be extended to transport quantities.
Information about the pairing mechanism for superconductivity is contained in the spectral weight for the anomalous (Gorkov) Green's function. In the most general case, this spectral weight can change sign on the positive real axis or even be complex in the presence of broken time-reversal symmetry. This renders impossible the direct analytic continuation with maximum entropy methods of numerical results obtained for the anomalous Green's function either in Matsubara frequency or in imaginary time. Here we show that one can define auxiliary Green's functions that allow one to extract the anomalous spectral weight in the most general multiorbital spin-singlet or triplet cases with no particular symmetry. As a simple example, we treat the case of lead in Eliashberg theory: Maximum entropy analytic continuation with auxiliary Green's functions (MaxEntAux) allows us to recover spectral functions that change sign on the positive real axis. The approach can be extended to transport quantities.
The computation of transport coefficients, even in linear response, is a major challenge for theoretical methods that rely on analytic continuation of correlation functions obtained numerically in Matsubara space. While maximum entropy methods can be used for certain correlation functions, this is not possible in general, important examples being the Seebeck, Hall, Nernst, and Reggi-Leduc coefficients. Indeed, positivity of the spectral weight on the positive real-frequency axis is not guaranteed in these cases. The spectral weight can even be complex in the presence of broken time-reversal symmetry. Various workarounds, such as the neglect of vertex corrections or the study of the infinite frequency or Kelvin limits, have been proposed. Here, we show that one can define auxiliary response functions that allow one to extract the desired real-frequency susceptibilities from maximum entropy methods in the most general multiorbital cases with no particular symmetry. As a benchmark case, we study the longitudinal thermoelectric response and corresponding Onsager coefficient in the single-band two-dimensional Hubbard model treated with dynamical mean-field theory and continuous-time quantum Monte Carlo. We thereby extend the maximum entropy analytic continuation with auxiliary functions (MaxEntAux method), developed for the study of the superconducting pairing dynamics of correlated materials, to transport coefficients.
This thesis is mainly a computational work studying the analytical continuation
of Green’s functions using the Maximum Entropy method. A strongly
correlated electron system is described with the single-band Hubbard model
and paramagnetic solutions are studied using Dynamic Mean Field Theory
on a Bethe lattice. Continuous Time Quantum Monte Carlo is used as Impurity
solver, for the infinite Anderson model at a finite temperature, to obtain
the Matsubara single-particle Green’s function propagator. Both metallic
and insulating spectral functions are obtained using the Maximum Entropy
Method. General properties of the Maximum Entropy Method as an analytic
continuation method from imaginary to real time are also discussed.
We present a generalization of the maximum entropy method to the analytic continuation of matrix-valued Green's functions. To treat off-diagonal elements correctly based on Bayesian probability theory, the entropy term has to be extended for spectral functions that are possibly negative in some frequency ranges. In that way, all matrix elements of the Green's function matrix can be analytically continued; we introduce a computationally cheap element-wise method for this purpose. However, this method cannot ensure important constraints on the mathematical properties of the resulting spectral functions, namely positive semidefiniteness and Hermiticity. To improve on this, we present a full matrix formalism, where all matrix elements are treated simultaneously. We show the capabilities of these methods using insulating and metallic dynamical mean-field theory (DMFT) Green's functions as test cases. Finally, we apply the methods to realistic material calculations for LaTiO3, where off-diagonal matrix elements in the Green's function appear due to the distorted crystal structure.
Analytic continuation is a critical step in quantum many-body computations, connecting imaginary-time or Matsubara Green's functions with real-frequency spectral functions, which can be directly compared to experimental results. However, because of the ill-posed nature of the analytic continuation problems, they have not been completely solved so far. In this paper, we suggest a simple, yet highly efficient method for analytic continuations of Matsubara Green's functions. This method takes advantage of barycentric rational functions to directly interpolate Matsubara Green's functions. At first, the nodes and weights of the barycentric rational functions are determined by the adaptive Antoulas-Anderson algorithm, avoiding reliance on nonconvex optimization. Next, the retarded Green's functions and the related spectral functions are evaluated by the resulting interpolants. We systematically explore the performance of this method through a series of toy models and realistic examples, comparing its accuracy and efficiency with other popular methods, such as the maximum entropy method. The benchmark results demonstrate that the new method can accurately reproduce not only continuous but also discrete spectral functions, irrespective of their positive definiteness. It works well even in the presence of intermediate noise, and outperforms traditional analytic continuation methods in computational speed. We believe that this method should stand out for its robustness against noise, broad applicability, high precision, and ultra-efficiency, offering a promising alternative to the maximum entropy method.
The presence of dynamical information is a feature distinguishing a finite-temperature quantum Monte Carlo simulation from a classical one. We now discuss numerical methods for extracting this information that use techniques and concepts borrowed from an area of probability theory called Bayesian statistical inference. The use of these techniques and concepts provided a solution to the very difficult problem of analytically continuing imaginary-time Green's functions, estimated by a quantum Monte Carlo simulation, to the real-time axis. Baym and Mermin (1961) proved that a unique mapping between these functions exists. However, executing this mapping numerically, with a simulation's incomplete and noisy data, transforms the problem into one without a unique solution and thus into a problem of finding a “best” solution according to some reasonable criterion. Instead of executing the analytic continuation between imaginary- and real-time Green's functions, thereby obtaining real-time dynamics, we instead estimate the experimentally relevant spectral density function these Green's functions share. We present three “best” solutions and emphasize that making the simulation data consistent with the assumptions of the numerical approach is a key step toward finding any of these best solutions. Preliminary comments The title of this chapter, “Analytic Continuation,” is unusual in the sense that it describes the task we wish to accomplish instead of the method we use to accomplish it. If we used the name of the method, the title would be something like “Bayesian Statistical Inference Using an Entropic Prior.” A shorter title would be “The Maximum Entropy Method.”We hope by the end of the chapter the reader will agree that using the short title is perhaps too glib and the longer one has meaningful content.
In this paper, we propose an analytic continuation method to extract real-frequency spectral functions from imaginary-frequency Green's functions of quantum many-body systems. This method is based on the pole representation of Matsubara Green's function and a stochastic sampling procedure is utilized to optimize the amplitudes and locations of poles. In order to capture narrow peaks and sharp band edges in the spectral functions, a constrained sampling algorithm and a self-adaptive sampling algorithm are developed. To demonstrate the usefulness and performance of this method, we at first apply it to study the spectral functions of representative fermionic and bosonic correlators. Then we employ this method to tackle the analytic continuation problems of matrix-valued Green's functions. The synthetic Green's functions, as well as realistic correlation functions from finite-temperature quantum many-body calculations, are used as input. The benchmark results demonstrate that this method is capable of reproducing most of the key characteristics in the spectral functions. The sharp, smooth, and multipeak features in both low- and high-frequency regions of spectral functions could be accurately resolved, which overcomes one of the main limitations of the traditional maximum entropy method. More importantly, it exhibits excellent robustness with respect to noisy and incomplete input data. The causality of spectral function is always satisfied even in the presence of sizable noises. As a by-product, this method could derive a fitting formula for the Matsubara data, which provides a compact approximation to the many-body Green's functions. Hence, we expect that this method could become a pivotal workhorse for numerically analytic continuation and be broadly useful in many applications.