aiqbench analysis
This paper establishes a rigorous mathematical proof that causal nonlinear response functions can be derived from imaginary-time-ordered Matsubara functions across all orders of perturbation theory, extending techniques previously limited to linear regimes. By solving equations of motion via induction, the authors provide explicit spectral density representations and derive generalized sum rules applicable to complex condensed matter systems like spin Hall insulators. This theoretical advancement offers a more efficient computational pathway for analyzing nonlinear dynamics in interacting quantum materials without relying solely on diagrammatic methods.
Imaginary Time Formalism for Causal Nonlinear Response Functions
arXiv:2506.21428v2 Announce Type: replace-cross Abstract: 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.
Abstract: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.
v2: accepted version. 27pgs
Mesoscale and Nanoscale Physics (cond-mat.mes-hall); Statistical Mechanics (cond-mat.stat-mech); Quantum Physics (quant-ph)
arXiv:2506.21428 [cond-mat.mes-hall]
(or arXiv:2506.21428v2 [cond-mat.mes-hall] for this version)
https://doi.org/10.48550/arXiv.2506.21428 Focus to learn more arXiv-issued DOI via DataCite
Phys. Rev. B 114, 045113 (2026)
https://doi.org/10.1103/7sn8-t9ds Focus to learn more DOI(s) linking to related resources
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.
Full text under the publisher's open licence · 390 words
Related coverage
Function-Valued Causal Influence in Nonlinear Time Series
arXiv · Oct 6, 2026

The Entangled Dimensions of Decision Theory
LessWrong · Jul 30, 2026
When are safe Pareto improvements individually rational?
LessWrong · Oct 13, 2026
Stop doing decision theory without metaphysics
LessWrong · Jul 20, 2026