This content will become publicly available on June 29, 2027

Title: Imaginary Time Formalism for Causal Nonlinear Response Functions
It is well established that causal linear response functions can be found by computing the much simpler imaginary time-ordered Matsubara functions and performing an analytic continuation. This principle is the basis for much of our understanding of linear response for interacting and disordered systems, via diagrammatic perturbation theory. Similar imaginary-time approaches have recently been introduced for computing nonlinear response functions as well, for example in [Annalen der Physik 536, 2300504 (2024); Physical Review X 11, 041006 (2021)], where the authors analytically continue the Matsubara functions to obtain the Keldysh response functions. In this work, we provide a proof of this connection to all orders in perturbation theory using an equation of motion based approach. We show by induction that causal nonlinear response functions at every order can be obtained from an analytic continuation of an appropriate time-ordered Matsubara function. We demonstrate this connection explicitly for second order response functions in the Lehmann representation. As a byproduct of our approach, we derive an explicit expression for the Lehmann representation of $$n$$-th order response functions by solving the equations of motion. We also use our result to find an analytic spectral density representation for both causal response functions and Matsubara functions. As an example, we apply our method to derive the non-linear $A^3$ term in the $SU(2)$ spin Hall response of an insulator with spin rotation symmetry. Finally, we show how our results lead to a family of generalized sum rules, focusing explicitly on the asymptotic expression for $$n$$-th harmonic generation rate. Our work opens the door to using imaginary time approaches to study nonlinear response functions in general condensed matter systems.  more » « less
Award ID(s):
2510219
PAR ID:
10693815
Author(s) / Creator(s):
;
Publisher / Repository:
American Physical Society
Date Published:
Journal Name:
Physical review B
ISSN:
2469-9969
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Analytic continuation is a central step in the simulation of finite-temperature field theories in which numerically obtained Matsubara data are continued to the real frequency axis for a physical interpretation. Numerical analytic continuation is considered to be an ill-posed problem where uncertainties on the Matsubara axis are amplified exponentially. Here, we present a systematic and controlled procedure that approximates any Matsubara function by a minimal pole representation to within a predefined precision. We then show a systematic convergence to the exact spectral function on the real axis as a function of our control parameter for a range of physically relevant setups. Our methodology is robust to noise and paves the way towards reliable analytic continuation in many-body theory and, by providing access to the analytic structure of the functions, a direct theoretical interpretation of physical properties. 
    more » « less
  3. A<sc>bstract</sc> We study two-point functions of symmetric traceless local operators in the bulk of de Sitter spacetime. We derive the Källén-Lehmann spectral decomposition for any spin and show that unitarity implies its spectral densities are nonnegative. In addition, we recover the Källén-Lehmann decomposition in Minkowski space by taking the flat space limit. Using harmonic analysis and the Wick rotation to Euclidean Anti de Sitter, we derive an inversion formula to compute the spectral densities. Using the inversion formula, we relate the analytic structure of the spectral densities to the late-time boundary operator content. We apply our technical tools to study two-point functions of composite operators in free and weakly coupled theories. In the weakly coupled case, we show how the Källén-Lehmann decomposition is useful to find the anomalous dimensions of the late-time boundary operators. We also derive the Källén-Lehmann representation of two-point functions of spinning primary operators of a Conformal Field Theory on de Sitter. 
    more » « less
  4. 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
  5. 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