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.
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.
Introduction There are countless situations in physics when one is interested in calculating the response of a system to a small time-dependent perturbation acting on it. With some luck the response can be expanded in a power series of the strength of the perturbation, so that, to first order, it is a linear function of the latter. To compute this function is the objective of the linear response theory (LRT). Linear response theory has many important applications to the study of electronic matter. Virtually all interactions of electrons with experimental probes (electromagnetic fields, beams of particles) can be regarded as small perturbations to the system: if they were not, one would not be probing the system, but the system modified by the probe. Consequently, the results of these experiments can be expressed in terms of linear response functions, which are properties of the unperturbed system. In particular it will turn out that the analytic structure of these functions is entirely determined by the eigenvalues and eigenfunctions of the unperturbed system. Conversely, a measure of the linear response as a function of the frequency of the perturbation enables us to determine the excitation energies of the system. Beside being a cornerstone for the theory of single-particle properties to be developed in Chapter 8, the linear response functions also provide invaluable information in their own right. For example, as we will discuss in Chapter 5, the extent to which an external electrostatic potential is reduced by screening is controlled by the dynamical dielectric function which, in turn, is determined by the density–density response function.
There exists a regime of overlap between the equilibrium and non-equilibrium behavior of a system, the non-equilibrium behavior of weakly perturbed states. When a system is perturbed ever so slightly, its response will be linear in the perturbation, say the current of the conduction electrons in a metal will be proportional to the strength of the applied electric field. This regime is called the linear response regime, and though the system is in a non-equilibrium state all its characteristics can be inferred from the properties of its equilibrium state. In the next chapter we shall go beyond the linear regime by showing how to obtain quantum kinetic equations. The kinetic-equation approach to transport is a general method, and allows in principle nonlinear effects to be considered. However, in many practical situations one is interested only in the linear response of a system to an external force. The linear response limit is a tremendous simplification in comparison with general non-equilibrium conditions, and is the subject matter of this chapter. In particular the linear response of the density and current of an electron gas are discussed. The symmetry properties of response functions, and the fluctuation–dissipation theorem are established. Lastly we demonstrate how correlation functions can be measured in scattering experiments, as illustrated by considering neutron scattering from matter. Needless to say, in measurements of (say) the current in a macroscopic body, far less information in the current correlation function is probed. Linear response In this section we consider the response of an arbitrary property of a system to a general perturbation.
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 A3 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.
The frequency-dependent exchange-correlation potential, which appears in the usual Kohn-Sham formulation of a time-dependent linear response problem, is a strongly nonlocal functional of the density, so that a consistent local density approximation generally does not exist. This problem can be avoided by choosing the current density as the basic variable in a generalized Kohn-Sham theory. This theory admits a local approximation which, for fixed frequency, is exact in the limit of slowly varying densities and perturbing potentials.
Time-Independent and Periodic Time-Dependent Potentials ...................... 73 6.3 Formal Foundations of Density Functional Theory for Time-Dependent Electric and Magnetic Fields ........................................... 74 6.4 Practical Schemes for the Calculation of Density and Current Density Variables in TDDFT ......................................................... 77 6.5 Linear Response within TDDFT ................................................................. 79 6.6 Concluding Remarks.................................................................................... 80 Acknowledgments................................................................................................... 80 References ............................................................................................................... 81 The study of behavior of many-electron systems such as atoms, molecules, and solids under the action of time-dependent (TD) external fields, which includes interaction with radiation, has been an important area of research. In the linear response regime, where one considers the external field to cause a small perturbation to the initial ground state of the system, one can obtain many important physical quantities such as polarizabilities, dielectric functions, excitation energies, photoabsorption spectra, van der Waals coefficients, etc. In many situations, for example, in the case of interaction of many-electron systems with strong laser field, however, it is necessary to go beyond linear response for investigation of the properties. Since a full theoretical description based on accurate solution of TD Schrodinger equation is not yet within the reach of computational capabilities, new methods which can efficiently handle the TD many-electron correlations need to be explored, and time-dependent density functional theory (TDDFT) is one such valuable approach.
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.