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

Physical Review A

Theory of nonlinear response

Akio Morita

Published Aug 1, 1986• 41 citations
Publisher

Abstract

The problem of nonlinear response is considered by employing a general time-evolution equation, and a Green's function which is the transition or conditional probability density for an unperturbed system. Expansion of the Green's function in terms of orthonormal functions enables us to express the distribution function describing the nonlinear behavior by means of matrix products whose elements are composed of correlation functions in the absence of the perturbation. In other words, it is shown how the distribution function induced by a strong perturbation may be calculated by knowing the Green's function without the perturbation. As the special case of the linear response, we have obtained Kubo's relation. The Laplace-transform technique with respect to time is found quite useful in developing the present theory in which the transient effect is also taken into account. As an application of the theory, a new relation valid in the region of the second-order perturbation connecting the transient rise and decay with the stationary alternating perturbations has been obtained.

DOI
10.1103/physreva.34.1499

Related papers

Linear response theory

Helmut Hofmann

Abstract This chapter presents a universal definition of linear response, which at first is studied for the example of the damped oscillator. From the quantal time evolution operator lowest order transition rates are deduced, together with microscopic response functions. Their fundamental properties are studied, in particular, the meaning of their reactive and dissipative parts for energy transfer. Correlation functions are introduced and related to response functions by way of the fluctuation dissipation theorem. For a random matrix model, the latter is shown to also hold true for microcanonical ensembles. Connections to thermal Green functions are established and an extension of response theory to unstable modes is discussed. With the help of the Mori product, relations among the static response, the adiabatic and the isothermal susceptibilities are given. Time dependent irreversible processes are studied by introducing relaxation functions and by following the variation of the density operator and of entropy. Onsager relations are presented and Kubo's expressions for transport coefficients are given.

•Apr 17, 2008•1 citation

Beyond Linear-Response Theory

Giancarlo Calvanese Strinati

The general expressions for the number density and current are first considered under the action of an external time-dependent perturbing potential of arbitrary strength and cast in terms of the single-particle lesser Green’s function. The expansion of the number density up to linear order in the perturbing potential is then considered, yielding the density–density correlation function of linear-response theory. A connection is also considered with the temperature correlation function of the Matsubara formalism via an analytic continuation in frequency space.

Cambridge University Press eBooks•Feb 6, 2026•0 citations

Dynamics of Nonlinear Stochastic Systems

Robert H. Kraichnan

A method for treating nonlinear stochastic systems is described which it is hoped will be useful in both the quantum-mechanical many-body problem and the theory of turbulence. In this method the true problem is replaced by models that lead to closed equations for correlation functions and averaged Green's functions. The model solutions are exact descriptions of possible dynamical systems, and, as a result, they display certain consistency properties. For example, spectral components of Green's functions which must be positive-definite in the true problem automatically are so for the models. The models involve a new stochastic element: Random couplings are introduced among an infinite collection of similar systems, the true problem corresponding to the limit where these couplings vanish. The method is first applied to a linear oscillator with random frequency parameter. The mean impulse-response function of the oscillator is obtained explicitly for two successive models. The results suggest the existence of a sequence of model solutions which converges rapidly to the exact solution of the true problem. Applications then are made to the Schrödinger equation of a particle in a random potential and to Burgers' analog for turbulence dynamics. For both problems, closed model equations are obtained which determine the average Green's function, the amplitude of the mean field, and the covariance of the fluctuating field. The model solutions can be expressed as sums of infinite classes of terms from the formal perturbation expansions of the solutions to the true problems. It is suggested that correspondence to stochastic models may be a useful criterion to help judge the validity of partial summations of perturbation series.

Journal of Mathematical Physics•Jan 1, 1961•367 citations

General linear response theory

No¨elle Pottier

Abstract Linear response theory is a general formalism applicable to any physical system slightly departing from equilibrium. It allows us to express the linear response functions and the generalized susceptibilities in terms of the equilibrium correlation functions of the relevant dynamical variables. For instance, in the case of the response to an applied external field, one of the two relevant dynamical variables represents the physical quantity conjugate to the field whereas the other represents the physical quantity whose out-of-equilibrium average is measured. The computation of the linear response functions relies on a first-order perturbation expansion of the density operator (or of the phase space distribution function) of the system coupled to the field. The expressions thus obtained constitute the Kubo formulas of linear response theory.

•Sep 17, 2009•0 citations

Perturbation expansions of the response function and the asymptotic regime

J. Besprosvany

Perturbation expansions of the response function S(q,\ensuremath{\omega}) of a many-body system in orders of any chosen component of the Hamiltonian are obtained. The response S(q,\ensuremath{\omega}) describes the inclusive scattering of a projectile interacting weakly with a target, where q and \ensuremath{\omega} are the momentum and energy transfer. Each term of a given order is written as an expectation value of products of paired Green's functions representing each initial- and final-state effects. For several choices of the perturbing operator, close expressions are obtained for the first orders of these expansions. The resulting terms represent partial summations over the Gersch series, whose relation to Green's-function expansions is established. Also, they are finite for singular potentials and are useful in the asymptotic regime of large q. In this regime they can give the dominant contribution to the response and initial- and final-state interactions of the hit particle and other particles in the system may be fully accounted for. Contributions of the n-particle density matrix can be separated for some expansions. In particular, an expression for the complete two-body contribution involving the hit particle is given.

Physical review. B, Condensed matter•May 1, 1991•12 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
The problem of nonlinear response is considered by employing a general time-evolution equation, and a Green's function which is the transition or conditional probability density for an unperturbed system. Expansion of the Green's function in terms of orthonormal functions enables us to express the distribution function describing the nonlinear behavior by means of matrix products whose elements are composed of correlation functions in the absence of the perturbation. In other words, it is shown how the distribution function induced by a strong perturbation may be calculated by knowing the Green's function without the perturbation. As the special case of the linear response, we have obtained Kubo's relation. The Laplace-transform technique with respect to time is found quite useful in developing the present theory in which the transient effect is also taken into account. As an application of the theory, a new relation valid in the region of the second-order perturbation connecting the transient rise and decay with the stationary alternating perturbations has been obtained.
Abstract This chapter presents a universal definition of linear response, which at first is studied for the example of the damped oscillator. From the quantal time evolution operator lowest order transition rates are deduced, together with microscopic response functions. Their fundamental properties are studied, in particular, the meaning of their reactive and dissipative parts for energy transfer. Correlation functions are introduced and related to response functions by way of the fluctuation dissipation theorem. For a random matrix model, the latter is shown to also hold true for microcanonical ensembles. Connections to thermal Green functions are established and an extension of response theory to unstable modes is discussed. With the help of the Mori product, relations among the static response, the adiabatic and the isothermal susceptibilities are given. Time dependent irreversible processes are studied by introducing relaxation functions and by following the variation of the density operator and of entropy. Onsager relations are presented and Kubo's expressions for transport coefficients are given.
The general expressions for the number density and current are first considered under the action of an external time-dependent perturbing potential of arbitrary strength and cast in terms of the single-particle lesser Green’s function. The expansion of the number density up to linear order in the perturbing potential is then considered, yielding the density–density correlation function of linear-response theory. A connection is also considered with the temperature correlation function of the Matsubara formalism via an analytic continuation in frequency space.
A method for treating nonlinear stochastic systems is described which it is hoped will be useful in both the quantum-mechanical many-body problem and the theory of turbulence. In this method the true problem is replaced by models that lead to closed equations for correlation functions and averaged Green's functions. The model solutions are exact descriptions of possible dynamical systems, and, as a result, they display certain consistency properties. For example, spectral components of Green's functions which must be positive-definite in the true problem automatically are so for the models. The models involve a new stochastic element: Random couplings are introduced among an infinite collection of similar systems, the true problem corresponding to the limit where these couplings vanish. The method is first applied to a linear oscillator with random frequency parameter. The mean impulse-response function of the oscillator is obtained explicitly for two successive models. The results suggest the existence of a sequence of model solutions which converges rapidly to the exact solution of the true problem. Applications then are made to the Schrödinger equation of a particle in a random potential and to Burgers' analog for turbulence dynamics. For both problems, closed model equations are obtained which determine the average Green's function, the amplitude of the mean field, and the covariance of the fluctuating field. The model solutions can be expressed as sums of infinite classes of terms from the formal perturbation expansions of the solutions to the true problems. It is suggested that correspondence to stochastic models may be a useful criterion to help judge the validity of partial summations of perturbation series.
Abstract Linear response theory is a general formalism applicable to any physical system slightly departing from equilibrium. It allows us to express the linear response functions and the generalized susceptibilities in terms of the equilibrium correlation functions of the relevant dynamical variables. For instance, in the case of the response to an applied external field, one of the two relevant dynamical variables represents the physical quantity conjugate to the field whereas the other represents the physical quantity whose out-of-equilibrium average is measured. The computation of the linear response functions relies on a first-order perturbation expansion of the density operator (or of the phase space distribution function) of the system coupled to the field. The expressions thus obtained constitute the Kubo formulas of linear response theory.
Perturbation expansions of the response function S(q,\ensuremath{\omega}) of a many-body system in orders of any chosen component of the Hamiltonian are obtained. The response S(q,\ensuremath{\omega}) describes the inclusive scattering of a projectile interacting weakly with a target, where q and \ensuremath{\omega} are the momentum and energy transfer. Each term of a given order is written as an expectation value of products of paired Green's functions representing each initial- and final-state effects. For several choices of the perturbing operator, close expressions are obtained for the first orders of these expansions. The resulting terms represent partial summations over the Gersch series, whose relation to Green's-function expansions is established. Also, they are finite for singular potentials and are useful in the asymptotic regime of large q. In this regime they can give the dominant contribution to the response and initial- and final-state interactions of the hit particle and other particles in the system may be fully accounted for. Contributions of the n-particle density matrix can be separated for some expansions. In particular, an expression for the complete two-body contribution involving the hit particle is given.
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.