Plasma Dispersion Function

If the unperturbed distribution function, $F_0$, appearing in Eq. (1007), is a Maxwellian then it is readily seen that, with a suitable scaling of the variables, the dispersion relation for electrostatic plasma waves can be expressed in terms of the function

$$Z(\zeta) = \pi^{-1/2}\int_{-\infty}^{\infty} \frac{{\rm e}^{-t^2}}{t-\zeta}\,dt,$$ (1022)

which is defined as it is written for ${\rm Im}(\zeta)>0$, and is analytically continued for ${\rm Im}(\zeta)\leq 0$. This function is known as the plasma dispersion function, and very often crops up in problems involving small-amplitude waves propagating through warm plasmas. Incidentally, $Z(\zeta)$ is the Hilbert transform of a Gaussian.

In view of the importance of the plasma dispersion function, and its regular occurrence in the literature of plasma physics, let us briefly examine its main properties. We first of all note that if we differentiate $Z(\zeta)$ with respect to $\zeta$ we obtain

$$Z'(\zeta) = \pi^{-1/2}\int_{-\infty}^{\infty} \frac{{\rm e}^{-t^2}}{(t-\zeta)^2}\, dt,$$ (1023)

which yields, on integration by parts,

$$Z'(\zeta) =-\pi^{-1/2}\int_{-\infty}^{\infty} \frac{2\,t}{t-\zeta}\, {\rm e}^{-t^2}\,dt = -2\,[1+\zeta\,Z].$$ (1024)

If we let $\zeta$ tend to zero from the upper half of the complex plane, we get

$$Z(0) = \pi^{-1/2}\,{\rm P}\!\int_{-\infty}^{\infty} \frac{{\rm e}^{-t^2}}{t}\,dt + {\rm i}\,\pi^{1/2} = {\rm i}\,\pi^{1/2}.$$ (1025)

Note that the principle part integral is zero because its integrand is an odd function of $t$.

Integrating the linear differential equation (1024), which possesses an integrating factor ${\rm e}^{\zeta^2}$, and using the boundary condition (1025), we obtain an alternative expression for the plasma dispersion function:

$$Z(\zeta) = {\rm e}^{-\zeta^2}\left({\rm i}\,\pi^{1/2} -2\!\int_0^\zeta {\rm e}^{x^2}\,dx\right).$$ (1026)

Making the substitution $t={\rm i}\,x$ in the integral, and noting that

$$\int_{-\infty}^0 {\rm e}^{-t^2}\,dt = \frac{\pi^{1/2}}{2},$$ (1027)

we finally arrive at the expression

$$Z(\zeta) = 2\,{\rm i}\, {\rm e}^{-\zeta^2}\int_{-\infty}^{{\rm i}\,\zeta} {\rm e}^{-t^2}\,dt.$$ (1028)

This formula, which relates the plasma dispersion function to an error function of imaginary argument, is valid for all values of $\zeta$.

For small $\zeta$ we have the expansion

$$Z(\zeta) = {\rm i}\,\pi^{1/2}\,{\rm e}^{-\zeta^2}-2\,\zeta\l... ...\frac{4\,\zeta^4}{15} - \frac{8\,\zeta^6}{105} +\cdots\right].$$ (1029)

For large $\zeta$, where $\zeta=x+{\rm i}\,y$, the asymptotic expansion for $x>0$ is written

$$Z(\zeta) \sim {\rm i}\,\pi^{1/2}\,\sigma\,{\rm e}^{-\zeta^2}... ...} + \frac{3}{4\,\zeta^4}+\frac{15}{8\,\zeta^6} +\cdots\right].$$ (1030)

Here,

$$\sigma = \left\{ \begin{array}{lll} 0&\mbox{\hspace{1cm}}& y... ...rt x\vert\\ [0.5ex] 2&&y< -1/\vert x\vert \end{array} \right..$$ (1031)

In deriving our expression for the Landau damping rate we have, in effect, used the first few terms of the above asymptotic expansion.

The properties of the plasma dispersion function are specified in exhaustive detail in a well-known book by Fried and Conte.