No Arabic abstract
We present $texttt{Maxent}$, a tool for performing analytic continuation of spectral functions using the maximum entropy method. The code operates on discrete imaginary axis datasets (values with uncertainties) and transforms this input to the real axis. The code works for imaginary time and Matsubara frequency data and implements the Legendre representation of finite temperature Greens functions. It implements a variety of kernels, default models, and grids for continuing bosonic, fermionic, anomalous, and other data. Our implementation is licensed under GPLv2 and extensively documented. This paper shows the use of the programs in detail.
We present the $texttt{TRIQS}/texttt{SOM}$ analytic continuation package, an efficient implementation of the Stochastic Optimization Method proposed by A. Mishchenko et al [Phys. Rev. B $textbf{62}$, 6317 (2000)]. $texttt{TRIQS}/texttt{SOM}$ strives to provide a high quality open source (distributed under the GNU General Public License version 3) alternative to the more widely adopted Maximum Entropy continuation programs. It supports a variety of analytic continuation problems encountered in the field of computational condensed matter physics. Those problems can be formulated in terms of response functions of imaginary time, Matsubara frequencies or in the Legendre polynomial basis representation. The application is based on the $texttt{TRIQS}$ C++/Python framework, which allows for easy interoperability with $texttt{TRIQS}$-based quantum impurity solvers, electronic band structure codes and visualization tools. Similar to other $texttt{TRIQS}$ packages, it comes with a convenient Python interface.
A method for analytic continuation of imaginary-time correlation functions (here obtained in quantum Monte Carlo simulations) to real-frequency spectral functions is proposed. Stochastically sampling a spectrum parametrized by a large number of delta-functions, treated as a statistical-mechanics problem, it avoids distortions caused by (as demonstrated here) configurational entropy in previous sampling methods. The key development is the suppression of entropy by constraining the spectral weight to within identifiable optimal bounds and imposing a set number of peaks. As a test case, the dynamic structure factor of the S=1/2 Heisenberg chain is computed. Very good agreement is found with Bethe Ansatz results in the ground state (including a sharp edge) and with exact diagonalization of small systems at elevated temperatures.
We present a Fourier Continuation-based parallel pseudospectral method for incompressible fluids in cuboid non-periodic domains. The method produces dispersionless and dissipationless derivatives with fast spectral convergence inside the domain, and with very high order convergence at the boundaries. Incompressibility is imposed by solving a Poisson equation for the pressure. Being Fourier-based, the method allows for fast computation of spectral transforms. It is compatible with uniform grids (although refined or nested meshes can also be implemented), which in turn allows for explicit time integration at sufficiently high Reynolds numbers. Using a new parallel code named SPECTER we illustrate the method with two problems: channel flow, and plane Rayleigh-Benard convection under the Boussinesq approximation. In both cases the method yields results compatible with previous studies using other high-order numerical methods, with mild requirements on the time step for stability.
Traditional maximum entropy and sparsity-based algorithms for analytic continuation often suffer from the ill-posed kernel matrix or demand tremendous computation time for parameter tuning. Here we propose a neural network method by convex optimization and replace the ill-posed inverse problem by a sequence of well-conditioned surrogate problems. After training, the learned optimizers are able to give a solution of high quality with low time cost and achieve higher parameter efficiency than heuristic full-connected networks. The output can also be used as a neural default model to improve the maximum entropy for better performance. Our methods may be easily extended to other high-dimensional inverse problems via large-scale pretraining.
We define the analytic continuation of the number of black hole microstates in Loop Quantum Gravity to complex values of the Barbero-Immirzi parameter $gamma$. This construction deeply relies on the link between black holes and Chern-Simons theory. Technically, the key point consists in writing the number of microstates as an integral in the complex plane of a holomorphic function, and to make use of complex analysis techniques to perform the analytic continuation. Then, we study the thermodynamical properties of the corresponding system (the black hole is viewed as a gas of indistinguishable punctures) in the framework of the grand canonical ensemble where the energy is defined a la Frodden-Gosh-Perez from the point of view of an observer located close to the horizon. The semi-classical limit occurs at the Unruh temperature $T_U$ associated to this local observer. When $gamma=pm i$, the entropy reproduces at the semi-classical limit the area law with quantum corrections. Furthermore, the quantum corrections are logarithmic provided that the chemical potential is fixed to the simple value $mu=2T_U$.