Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2101.04136

Nonlinear terahertz electro-optical responses in centrosymmetric electronic systems

Indranil Paul

Published Jan 11, 2021• 1 citation
arXivHTMLPDF

Abstract

Motivated by the recent developments in terahertz spectroscopy using pump-probe setups to study correlated electronic materials, we review the field theoretical formalism to compute finite frequency nonlinear electro-optical responses in centrosymmetric systems starting from basic time dependent perturbation theory. We express the nonlinear current kernel as a sum of several causal response functions. These causal functions cannot be evaluated using perturbative field theory methods, since they are not contour ordered. Consequently, we associate each response function with a corresponding imaginary time ordered current correlation function, since the latter can be factorized using Wick's theorem. The mapping between the response functions and the correlation functions, suitably analytically continued to real frequencies, is proven exactly. We derive constraints satisfied by the nonlinear current kernel and we prove a generalized $f$-sum rule for the nonlinear conductivity, all of which are consequences of particle number conservation. The constraints guarantee that the nonlinear static responses are free from spurious divergences. We apply the theory to compute the gauge invariant nonlinear conductivity of a system of noninteracting electrons in the presence of weak disorder. As special cases of this generalized nonlinear response, we discuss its third harmonic and its instantaneous terahertz Kerr signals. The formalism can be used to compute the nonlinear conductivity in symmetry broken phases of electronic systems such as superconductors, density waves and nematic states.

DOI
10.48550/arxiv.2101.04136

Related papers

Field theoretic description of nonlinear electro-optical responses in centrosymmetric electronic systems

Indranil Paul

Abstract Motivated by the recent developments in terahertz spectroscopy using pump-probe setups to study correlated electronic materials, we review the field theoretical formalism to compute finite frequency nonlinear electro-optical responses in centrosymmetric systems starting from basic time dependent perturbation theory. We express the nonlinear current kernel as a sum of several causal response functions. These causal functions cannot be evaluated using perturbative field theory methods, since they are not contour ordered. Consequently, we associate each response function with a corresponding imaginary time ordered current correlation function, since the latter can be factorized using Wick’s theorem. The mapping between the response functions and the correlation functions, suitably analytically continued to real frequencies, is proven exactly. We derive constraints satisfied by the nonlinear current kernel and we prove a generalized f-sum rule for the nonlinear conductivity, all of which are consequences of particle number conservation. The constraints guarantee that the nonlinear static responses are free from spurious divergences. We apply the theory to compute the gauge invariant nonlinear conductivity of a system of noninteracting electrons in the presence of weak disorder. As special cases of this generalized nonlinear response, we discuss its third harmonic and its instantaneous terahertz Kerr signals. The formalism can be used to compute the nonlinear conductivity in symmetry broken phases of electronic systems such as superconductors, density waves and nematic states.

Journal of Physics Condensed Matter•Jul 19, 2024•2 citations

Nonlinear terahertz electro-optical responses in metals

Indranil Paul

Motivated by the recent developments in terahertz spectroscopy using pump-probe setups, we develop the theory of finite frequency nonlinear electro-optical responses in centrosymmetric metals starting from basic time dependent perturbation theory. We express the nonlinear current kernel as a sum of several causal response functions. These functions cannot be evaluated using perturbative field theory methods. Consequently, we associate each response function with an imaginary time ordered current correlation function, which can be factorized using Wick's theorem. The mapping between the response functions and the correlation functions, suitably analytically continued to real frequencies, is proven exactly. We derive constraints satisfied by the nonlinear current kernel and we prove a generalized $f$-sum rule for the nonlinear conductivity, all of which are consequences of particle number conservation. We apply the theory to compute the gauge invariant nonlinear conductivity of a system of noninteracting electrons in the presence of weak disorder. As special cases of this generalized nonlinear response, we discuss its third harmonic and its instantaneous terahertz Kerr signals. Our formalism can be used to compute the nonlinear conductivity in symmetry broken phases such as density waves and nematic states.

•Jan 11, 2021•1 citation

On the General Properties of Non-linear Optical Conductivities

Haruki Watanabe, Yankang Liu, Masaki Oshikawa

Abstract The optical conductivity is the basic defining property of materials characterizing the current response toward time-dependent electric fields. In this work, following the approach of Kubo’s response theory, we study the general properties of the nonlinear optical conductivities of quantum many-body systems both in equilibrium and non-equilibrium. We obtain an expression of the second- and the third-order optical conductivity in terms of correlation functions and present a perturbative proof of the generalized Kohn formula proposed recently. We also discuss a generalization of the f-sum rule to a non-equilibrium setting by focusing on the instantaneous response.

arXiv:2004.04561•Journal of Statistical Physics•Oct 18, 2020•22 citations

Tangent equations of motion for nonlinear response functions

Atsushi Ono

Nonlinear response functions, formulated as multipoint correlation functions or Volterra kernels, encode the dynamical and spectroscopic properties of physical systems and underpin a wide range of nonlinear transport and optical phenomena. However, their evaluation rapidly becomes prohibitive at high orders because of combinatorial (often factorial) scaling or severe numerical errors. Here, we establish a systematic and efficient framework to compute nonlinear response functions directly from real-time dynamics, without explicitly constructing multipoint correlators or relying on numerically unstable finite-difference methods for order-resolved extraction. Our approach is based on the Gateaux derivative with respect to the external field in function space, which yields a closed hierarchy of tangent equations of motion (TEOM). Propagating the TEOM alongside the original dynamics isolates each perturbative order with high accuracy, providing a term-by-term decomposition of physical contributions. The computational cost scales exponentially with response order in the fully general setting and reduces to polynomial complexity when all perturbation directions are identical; both regimes avoid the factorial scaling of explicit multipoint-correlator evaluations. We demonstrate the power of TEOM by computing frequency-resolved fifth-order response functions for a solid-state electron model and by obtaining nonlinear response functions up to the 49th order with controlled accuracy in a classical Duffing oscillator. We further show that our time-evolution formulation allows optical conductivities to be evaluated directly while remaining numerically stable even near zero frequency. TEOM can be incorporated seamlessly into existing real-time evolution methods, yielding a general framework for computing nonlinear response functions in quantum and classical dynamical systems.

cond-mat.str-elCross-listedarXiv:2603.22137•Physical Review Research•Mar 24, 2026•0 citations

Gauge invariance and Ward identities in nonlinear response theory

Habib Rostami, Mikhail I. Katsnelson, Giovanni Vignale +1 more

We present a formal analysis of nonlinear response functions in terms of correlation functions in real- and imaginary-time domains. In particular, we show that causal nonlinear response functions, expressed in terms of nested commutators in real time, can be obtained from the analytic continuation of time-ordered response functions, which are more easily amenable to diagrammatic calculation. This generalizes the well-known result of linear response theory. We then use gauge invariance arguments to derive exact relations between second-order response functions in density and current channels. These identities, which are non-perturbative in the strength of inter-particle interactions, allow us to establish exact connections between nonlinear optics calculations done in different electromagnetic gauges.

arXiv:2102.04425•Annals of Physics•May 26, 2021•40 citations

Imaginary time formalism for causal nonlinear response functions

Sounak Sinha, Barry Bradlyn

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.

arXiv:2506.21428•Physical review. B./Physical review. B•Jun 24, 2026•0 citations
Motivated by the recent developments in terahertz spectroscopy using pump-probe setups to study correlated electronic materials, we review the field theoretical formalism to compute finite frequency nonlinear electro-optical responses in centrosymmetric systems starting from basic time dependent perturbation theory. We express the nonlinear current kernel as a sum of several causal response functions. These causal functions cannot be evaluated using perturbative field theory methods, since they are not contour ordered. Consequently, we associate each response function with a corresponding imaginary time ordered current correlation function, since the latter can be factorized using Wick's theorem. The mapping between the response functions and the correlation functions, suitably analytically continued to real frequencies, is proven exactly. We derive constraints satisfied by the nonlinear current kernel and we prove a generalized fff-sum rule for the nonlinear conductivity, all of which are consequences of particle number conservation. The constraints guarantee that the nonlinear static responses are free from spurious divergences. We apply the theory to compute the gauge invariant nonlinear conductivity of a system of noninteracting electrons in the presence of weak disorder. As special cases of this generalized nonlinear response, we discuss its third harmonic and its instantaneous terahertz Kerr signals. The formalism can be used to compute the nonlinear conductivity in symmetry broken phases of electronic systems such as superconductors, density waves and nematic states.
Abstract Motivated by the recent developments in terahertz spectroscopy using pump-probe setups to study correlated electronic materials, we review the field theoretical formalism to compute finite frequency nonlinear electro-optical responses in centrosymmetric systems starting from basic time dependent perturbation theory. We express the nonlinear current kernel as a sum of several causal response functions. These causal functions cannot be evaluated using perturbative field theory methods, since they are not contour ordered. Consequently, we associate each response function with a corresponding imaginary time ordered current correlation function, since the latter can be factorized using Wick’s theorem. The mapping between the response functions and the correlation functions, suitably analytically continued to real frequencies, is proven exactly. We derive constraints satisfied by the nonlinear current kernel and we prove a generalized f-sum rule for the nonlinear conductivity, all of which are consequences of particle number conservation. The constraints guarantee that the nonlinear static responses are free from spurious divergences. We apply the theory to compute the gauge invariant nonlinear conductivity of a system of noninteracting electrons in the presence of weak disorder. As special cases of this generalized nonlinear response, we discuss its third harmonic and its instantaneous terahertz Kerr signals. The formalism can be used to compute the nonlinear conductivity in symmetry broken phases of electronic systems such as superconductors, density waves and nematic states.
Motivated by the recent developments in terahertz spectroscopy using pump-probe setups, we develop the theory of finite frequency nonlinear electro-optical responses in centrosymmetric metals starting from basic time dependent perturbation theory. We express the nonlinear current kernel as a sum of several causal response functions. These functions cannot be evaluated using perturbative field theory methods. Consequently, we associate each response function with an imaginary time ordered current correlation function, which can be factorized using Wick's theorem. The mapping between the response functions and the correlation functions, suitably analytically continued to real frequencies, is proven exactly. We derive constraints satisfied by the nonlinear current kernel and we prove a generalized fff-sum rule for the nonlinear conductivity, all of which are consequences of particle number conservation. We apply the theory to compute the gauge invariant nonlinear conductivity of a system of noninteracting electrons in the presence of weak disorder. As special cases of this generalized nonlinear response, we discuss its third harmonic and its instantaneous terahertz Kerr signals. Our formalism can be used to compute the nonlinear conductivity in symmetry broken phases such as density waves and nematic states.
Abstract The optical conductivity is the basic defining property of materials characterizing the current response toward time-dependent electric fields. In this work, following the approach of Kubo’s response theory, we study the general properties of the nonlinear optical conductivities of quantum many-body systems both in equilibrium and non-equilibrium. We obtain an expression of the second- and the third-order optical conductivity in terms of correlation functions and present a perturbative proof of the generalized Kohn formula proposed recently. We also discuss a generalization of the f-sum rule to a non-equilibrium setting by focusing on the instantaneous response.
Nonlinear response functions, formulated as multipoint correlation functions or Volterra kernels, encode the dynamical and spectroscopic properties of physical systems and underpin a wide range of nonlinear transport and optical phenomena. However, their evaluation rapidly becomes prohibitive at high orders because of combinatorial (often factorial) scaling or severe numerical errors. Here, we establish a systematic and efficient framework to compute nonlinear response functions directly from real-time dynamics, without explicitly constructing multipoint correlators or relying on numerically unstable finite-difference methods for order-resolved extraction. Our approach is based on the Gateaux derivative with respect to the external field in function space, which yields a closed hierarchy of tangent equations of motion (TEOM). Propagating the TEOM alongside the original dynamics isolates each perturbative order with high accuracy, providing a term-by-term decomposition of physical contributions. The computational cost scales exponentially with response order in the fully general setting and reduces to polynomial complexity when all perturbation directions are identical; both regimes avoid the factorial scaling of explicit multipoint-correlator evaluations. We demonstrate the power of TEOM by computing frequency-resolved fifth-order response functions for a solid-state electron model and by obtaining nonlinear response functions up to the 49th order with controlled accuracy in a classical Duffing oscillator. We further show that our time-evolution formulation allows optical conductivities to be evaluated directly while remaining numerically stable even near zero frequency. TEOM can be incorporated seamlessly into existing real-time evolution methods, yielding a general framework for computing nonlinear response functions in quantum and classical dynamical systems.
We present a formal analysis of nonlinear response functions in terms of correlation functions in real- and imaginary-time domains. In particular, we show that causal nonlinear response functions, expressed in terms of nested commutators in real time, can be obtained from the analytic continuation of time-ordered response functions, which are more easily amenable to diagrammatic calculation. This generalizes the well-known result of linear response theory. We then use gauge invariance arguments to derive exact relations between second-order response functions in density and current channels. These identities, which are non-perturbative in the strength of inter-particle interactions, allow us to establish exact connections between nonlinear optics calculations done in different electromagnetic gauges.
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 nnn-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 A3A^3A3 term in the SU(2)SU(2)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 nnn-th harmonic generation rate. Our work opens the door to using imaginary time approaches to study nonlinear response functions in general condensed matter systems.