An official website of the United States government
We introduce a Julia implementation of the recently proposed Nevanlinna analytic continuation method. The method is based on Nevanlinna interpolants and inherently preserves the causality of a response function due to its construction. For theoretical calculations without statistical noise, this continuation method is a powerful tool to extract real-frequency information from numerical input data on the Matsubara axis. This method has been applied to first-principles calculations of correlated materials. This paper presents its efficient and full-featured open-source implementation of the method including the Hamburger moment problem and smoothing.
more »
« less
TRIQS/Nevanlinna: Implementation of the Nevanlinna Analytic Continuation method for noise-free data
We present the TRIQS/Nevanlinna analytic continuation package, an efficient implementation of the methods proposed by J. Fei et al. (2021) [53] and (2021) [55]. TRIQS/Nevanlinna 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 based analytic continuation programs. With the additional Hardy functions optimization procedure, it allows for an accurate resolution of wide band and sharp features in the spectral function. Those problems can be formulated in terms of imaginary time or Matsubara frequency response functions. The application is based on the TRIQS C++/Python framework, which allows for easy interoperability with other TRIQS-based applications, electronic band structure codes and visualization tools. Similar to other TRIQS packages, it comes with a convenient Python interface.
more »
« less
- Award ID(s):
- 2001465
- PAR ID:
- 10616825
- Publisher / Repository:
- Mendeley Data
- Date Published:
- Subject(s) / Keyword(s):
- Condensed Matter Physics FOS: Physical sciences Computational Physics Analytic Continuation
- Format(s):
- Medium: X
- Right(s):
- GNU Public License Version 3
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract In this paper, we introduce the notion of an Evertse–Ferretti Nevanlinna constant and compare it with the birational Nevanlinna constant introduced by the authors in a recent joint paper. We then use it to recover several previously known results. This includes a 1999 example of Faltings from hisBaker’s Gardenarticle. We also extend the theory of these Nevanlinna constants to what we call “multidivisor Nevanlinna constants,” which allow the proximity function to involve the maximum of Weil functions for finitely many divisors.more » « less
-
We present a minimal pole method for analytically continuing matrix-valued imaginary frequency correlation functions to the real axis, enabling precise access to off-diagonal elements and thus improving the interpretation of self-energies and susceptibilities in quantum simulations. Traditional methods for matrix-valued analytic continuation tend to be either noise sensitive or make ad hoc positivity assumptions. Our approach avoids these issues via the construction of a compact pole representation with shared poles through exponential fits, expanding upon prior work focused on scalar functions. We test our method across various scenarios, including fermionic and bosonic response functions, with and without noise, and for both continuous and discrete spectra of real materials and model systems. Our findings demonstrate that this technique addresses the shortcomings of existing methodologies, such as artificial broadening and positivity violations. The paper is supplemented with a sample implementation in PYTHON.more » « less
-
Quantum mechanics has about a dozen exactly solvable potentials. Normally, the time-independent Schrödinger equation for them is solved by using a generalized series solution for the bound states (using the Fröbenius method) and then an analytic continuation for the continuum states (if present). In this work, we present an alternative way to solve these problems, based on the Laplace method. This technique uses a similar procedure for the bound states and for the continuum states. It was originally used by Schrödinger when he solved the wave functions of hydrogen. Dirac advocated using this method too. We discuss why it is a powerful approach to solve all problems whose wave functions are represented in terms of confluent hypergeometric functions, especially for the continuum solutions, which can be determined by an easy-to-program contour integral.more » « less
-
We propose a data-driven learning framework for the analytic continuation problem in numerical quantum many-body physics. Designing an accurate and efficient framework for the analytic continuation of imaginary time using computational data is a grand challenge that has hindered meaningful links with experimental data. The standard Maximum Entropy (MaxEnt)-based method is limited by the quality of the computational data and the availability of prior information. Also, the MaxEnt is not able to solve the inversion problem under high level of noise in the data. Here we introduce a novel learning model for the analytic continuation problem using a Adams-Bashforth residual neural network (AB-ResNet). The advantage of this deep learning network is that it is model independent and, therefore, does not require prior information concerning the quantity of interest given by the spectral function. More importantly, the ResNet-based model achieves higher accuracy than MaxEnt for data with higher level of noise. Finally, numerical examples show that the developed AB-ResNet is able to recover the spectral function with accuracy comparable to MaxEnt where the noise level is relatively small.more » « less
