Thomson scattering in inhomogeneous plasmas: The Role of the Fluctuation-Dissipation Theorem.

A self-consistent kinetic theory of Thomson scattering of an electromagnetic field by a non-uniform plasma is derived. We draw the readers' attention to the inconsistency in recent results on the Thomson scattering in inhomogeneous plasma, which leads to violation of the Fluctuation-Dissipation Theorem. We show, that not only the imaginary part, but also the derivatives of the real part of the dielectric susceptibility determine the amplitude and the width of the Thomson scattering spectral lines. As a result of inhomogeneity, these properties become asymmetric with respect to inversion of the sign of the frequency. A method is proposed for measuring local gradients of the electron density with the aid of Thomson scattering.Arising from: P. Kozlowski, et al. Sci. Rep. 6, 24283 (2016); https://doi.org/10.1038/srep24283 .

Introduction

When an electromagnetic wave propagates in a plasma, its interaction with fluctuational oscillations of the plasma may result in scattering of the wave, which can be accompanied by a change in its frequency and wave vector. The intensity of scattered waves depends on both the intensity of the incident wave and the level of plasma fluctuations. Since the spectrum of plasma fluctuations exhibits sharp maxima at proper plasma frequencies, the spectrum of scattered waves will also exhibit sharp maxima at frequencies differing from the frequency of the incident wave by the according frequencies of the plasma fluctuations. The shift, width and shape of spectral lines carry information on such parameters of the plasma as its density, temperature, mean velocity, ion composition etc. A method of remote probing of a plasma, termed Thomson scattering, is a powerful plasma diagnostic tool that is widely employed in measurements of plasma parameters over a fairly broad range of plasma densities from the ionosphere to strongly coupled plasma. In such measurements the plasma must be transparent to the probe electromagnetic radiation. This may be microwave[1], laser[2] or X-ray radiation. A comprehensive exposition of the state-of-the-art of X -ray Thomson scattering is presented in Review[3].

The differential Thomson scattering cross section, within an elementary solid angle and for a frequency interval is described by the expression[4,5]:where , are the wave vectors and the frequencies of the incident and scattered electromagnetic fields. Thus, the problem reduces to finding the spectral characteristics of electron density fluctuations - the dynamic electron structure factor. The theory of equilibrium and nonequilibrium plasma fluctuations was successfully developed in the second half of the past century[5-9]. In accordance with the Poisson equation, the dynamic electron structure factor in a spatially homogeneous system is directly linked to the electrostatic field fluctuations. In thermodynamic equilibrium, the electrostatic field fluctuations satisfy the famous Callen-Welton Fluctuation-Dissipation Theorem (FDT)[10], linking their intensity to the imaginary part of the dielectric function ε(ω, k) and to the temperature TEq. (2) refers to the steady state, for a space uniform system. However, it is not evident that the plasma parameters can be kept constant in both space and time. Inhomogeneities in space and time of these quantities will certainly also contribute to the fluctuations. Hence it is challenging to formulate the generalization of the FDT for inhomogeneous plasma and reformulate accordingly the results for the Thomson scattering.

An attempt to solve this important problem of describing Thomson scattering in an inhomogeneous plasma has been made recently in[11]. The authors proposed the following ad hoc generalization for the quantity S(k, ω):“where is the dynamic structure factor for an ideal (noniteracting gas), and the dielectric (screening) function ε(k, ω) in the denominator of Eq. (3) in a first order gradient expansion in microscopic variable is:χ is the susceptibility of the ideal Coulomb plasma. The index “eq” labels the susceptibility for a homogeneous system in thermodynamic equilibrium”[11]. Noteworthily, while the authors applied the expansion for the denominator, the numerator in Eq. (3) has not been correspondingly expanded. Although, this approximation based on the “physical intuition” reflects some properties of the system, it fails, unfortunately, to satisfy the basic principles. This entails a dramatic inaccuracy of this approach. Namely, this resulted in two consequences: Firstly, the obtained result is nonphysical, since it contradicts FDT in the local equilibrium state. The FDT for a local equilibrium state was proved by Balescu[12]. The parameters of a system in a local equilibrium state can be changed adiabatically on a scale greater than the particle mean free path. Inhomogeneity and nonstationarity of plasma fluctuations are manifested via a nonlocal dependence upon time[13] and coordinates[14]. The FDT for a nonlocal plasma was given in our paper[15]. Secondly, the obtained correction due to the inhomogeneity in the denominator Eq. (3) is erroneous for Langmuir oscillations, especially in the case of small wave numbers k < k, which usually occurs in experiments. And last but not least: the rigorous kinetic theory predicts asymmetry of spectral lines in an inhomogeneous plasma. Such asymmetry has been indeed detected in spectroscopic studies of plasma flows in magnetic traps[16,17].

In the present paper, applying the Klimontovich-Langevin approach[18] and the time-space multiscale technique, we show that not only the imaginary part but also the derivatives of the real part of the dielectric susceptibility determine the amplitude and width of spectral lines of the electrostatic field fluctuations and of the dynamic electron structure factor, as well. As a result of the inhomogeneity, these properties become asymmetric with respect to inversion of the sign of the frequency. In the kinetic regime the structure factor is more sensitive to space gradients than the spectral function of the electrostatic field fluctuations. Note that for simple fluids and gases a general theory of hydrodynamic fluctuations for nonequilibrium stationary inhomogeneous states has been developed in[19,20]. In particular, it has been found that there exists an asymmetry of the spectrum for Brillouin scattering from a fluid in a shear flow or in a temperature gradient. The situation for the plasma problem we are considering is, however, quite different.

Results

To treat the problem, a kinetic approach is required, especially when the wavelength of the fluctuations is larger than the Debye wavelength. To derive nonlocal expressions for the spectral function of the electrostatic field fluctuation and for the dynamic electron structure factor, we adopt the Klimontovich-Langevin approach to describe kinetic fluctuations[18]. A kinetic equation for the fluctuation δ f of the one-particle distribution function (DF) with respect to the reference state f is considered. In the general case the reference state is a nonequilibrium DF which varies in space and time both on the kinetic scale (mean free path l and interparticle collision time ) and, also, on the larger hydrodynamic scales. These scales are much larger than the characteristic fluctuation time ω−1. In the nonequilibrium case we can, therefore, introduce a small parameter , which allows us to describe fluctuations on the basis of a multiple space and time scale analysis. Obviously, the fluctuations vary on both the “ fast” (r, t) and the “slow” (μr, μt) space and time scales: and Here x stands for the phase-space coordinates (r, p). The Langevin kinetic equation for has the form[18]where e is the charge of the particle of species a, δE is the electrostatic field fluctuation, and the operator is defined byand is the linearized collision operator. A model collision operator for plasma is presented in[21].

The Langevin source in Eq. (5) is determined by the following equation:The solution of Eq. (5) has the formwhere the Green function of the operator is determined bywith the causality condition:when .

Thus, and are connected by the relation:

For stationary and spatially uniform systems the DF f and the operator do not depend on time and space. In this case, the dependence on time and space of the Green function is manifested only through the difference and . However, when the DF and are slowly varying quantities in time and space, and when nonlocal effects are considered, the time and space dependence of is more subtle:For the homogeneous case this non-trivial result was obtained for the first time in our previous work[13]. This result was extended to inhomogeneous systems[14]. Here we want to stress that the nonlocal effects appear due to the slow time and space dependencies and .

Relationship (12) is directly linked with the constitutive relation between the electric displacement and the electric fieldPreviously two kinds of constitutive relations were proposed phenomenologically for a weakly inhomogeneous and slowly time-varying medium. Kadomtsev[22] formulated the so-called symmetrized constitutive relation.Rukhadze and Silin[23] proposed a nonsymmetrized constitutive relationBoth phenomenological formulations are unsatisfactory. The correct expression should beIn the first order, expansion with respect to μ, Eq. (8) leads towith and .

Using the Poisson equation and performing the Fourier-Laplace transformation for the fast variableswe have that the spectral function of the nonequilibrium electrostatic field fluctuations assumes the form:where we introduced the effective dielectric function as:whereis the susceptibility for a nonequilibrium plasma.

If the authors of paper[11] deal with the nonequilibrium case, they should use an expression similar to Eq. (19) in Eq. (3) as well as the nonequilibrium susceptibility (21).

For the local equilibrium case, where the reference state is Maxwellian, we have the identity:and Eq. (19) takes the formwhere . For the case of equal temperatures and Va = 0, Eq. (23) satisfies the FDT. In this case the small parameter μ is determined on the slower hydrodynamic scale. Imaginary part of the susceptibility determines the width of the spectral line near the resonance:In Eq. (24) there appear additional first-order terms of the small parameter μ. It is important to note that the imaginary part of the dielectric susceptibility is now replaced by the real part, which in the plasma resonance may be greater than the imaginary part by the same factor μ−1. Therefore, the second and third terms in Eq. (24) in the kinetic regime have an effect comparable to that of the first term. Second-order corrections in the expansion in μ only appear in the imaginary part of the susceptibility, and they can be reasonably neglected. It is therefore sufficient to retain the first-order corrections to resolve the problem. The width of the spectral lines Eq. (24) is affected by new nonlocal terms. They are not related to Joule dissipation and appear because of an additional phase shift between the induction vector and the electric field. This phase shift results from the finite time needed to set the polarization in the plasma with dispersion[24]. Such a phase shift in the plasma with space dispersion appears due to the medium inhomogeneity.

For the case where the system parameters are homogeneous in space but vary in time, the correction to the width of the spectral lines in Eq. (24) is still symmetric with respect to the change in sign of ω. However, when the plasma parameters are space dependent this symmetry is lost. The real part of the susceptibility in Eq. (24) is an even function of ω. This property implies that the contribution of the space derivative to the expression for the width of the spectral lines is odd function of ω. Moreover, this term gives rise to an anisotropy in k space.

Let us estimate this correction for the plasma mode andOn the hydrodynamic scale , and γ > 0.

For the spatially homogeneous case there is no difference between the spectral properties of the longitudinal electric field and of the electron density, because they are related by the Poisson equation. This statement is no longer valid when an inhomogeneous plasma is considered. Indeed the longitudinal electric field is linked to the particle density by the nonlocal relation:In the same approximation as in Eq. (23) the expression for the dynamic electron structure factor for a two-component local equilibrium plasma has the form[15]:where k is the inverse Debye length,The inhomogeneous correction in Eq. (30) is not the same as in Eq. (24) . The origin of this difference is that the Green functions for the electrostatic field fluctuations and particle density fluctuations are not the same in an inhomogeneous situation. As above, we can expand near the plasma resonance . Thus, for the Langmuir linewhereis the width of the dynamic electron structure factor. An estimate for the plasma mode is then:From this equation we see that the inhomogeneous correction in Eq. (33) is greater than the one in Eq. (26) by the factor . For the same inhomogeneity, i.e., the same gradient of the density, we plot the dynamic electron structure factor together with the as functions of frequency (Fig. 1). This figure shows that the asymmetry of the spectral lines is present both for and . However, this effect is more pronounced in than in . Such asymmetry has been indeed detected in inhomogeneous plasma[16,17]. The asymmetry of lines can be used as a new diagnostic tool for measuring local gradients in plasmas by Thomson scattering.

Figure 1

The dynamic electron structure factor (solid line) and the spectral function of electrostatic field fluctuations (dashed line) as a function of frequency. ; .

The Langmuir line (31) takes the Lorentz form. The amplitude of the spectral line A is inversely proportional to its widthFrom Eq. (34) and Eq. (33) quite a simple formula for calculation of the electron density gradient from the Thomson scattering spectrum follows:here and are the amplitudes and the half-widths of the red and blue Langmuir satellites, respectively (Fig. 1).

Thus, intensity and width measurements of the red and blue lines of the spectrum allow to determine the scalar product of the electron density gradient and the scattering vector at a given point. To determine the vector it is sufficient to measure the radiation scattered in three directions simultaneously.

Similar calculations can also be performed for degenerate high density plasma[3].

Conclusion

A first-principle kinetic theory of Thomson scattering in a non-uniform plasma is constructed, which agrees with the basic FDT and provides quantitatively correct results, that have been confirmed experimentally[16,17]. Moreover, our theory provides a novel and unique method for the remote probing and measurement of electron density gradients in a plasma; this is based on the demonstrated asymmetry of the Thomson scattering lines. The latter may be important for numerous technological applications, e.g. for the tokamak[25], for the high energy density plasma[3] etc. Our findings are in a sharp contrast with the results of the recent publication[11], where the suggested ad hoc theory did not agree with the local FDT (which proved to hold valid) and led to quantitatively (and qualitatively) incorrect predictions.