Search for: All records

Thomson scatter radars have successfully measured plasma parameters in the ionosphere for over 60 years. Fundamentally, the radars measure increased power returns when the Bragg scattering condition is met by a source of density fluctuations in the plasma. Typically, wave modes of the plasma provide the source of structuring, and the radars measure strong power returns at the ion line which is associated with the ion-acoustic mode, the gyro line which is associated with the electrostatic whistler mode, and the plasma line that comes from the Langmuir mode. However, the existence of an ion-acoustic mode or electrostatic whistler mode is not guaranteed in the ionosphere. In this study, a formalism is developed to explain non-resonant wave modes as features occurring at frequencies where the dielectric function has a local minimum as opposed to a root corresponding to the typical resonant wave mode. With this formalism, the frequency of non-resonant waves is numerically solved as a function of basic plasma parameters. By solving for minima of the dielectric function, the frequency and intensity of gyro lines is determined for a wide range of plasma temperatures and densities. This analysis explains why Arecibo gyro lines are typically weak in intensity and result from non-resonant waves. For VHF systems like EISCAT, gyro lines are shown to be strong spectral peaks corresponding to standard resonant solutions for electrostatic whistler waves. 

Kinetic plasma studies often require computing integrals of the velocity distribution over a complex-valued pole. The standard method is to solve the integral in the complex plane using the Plemelj theorem, resulting in the standard plasma dispersion function for Maxwellian plasmas. For non-Maxwellian plasmas, the Plemelj theorem does not generalize to an analytic form, and computational methods must be used. In this paper, a new computational method is developed to accurately integrate a non-Maxwellian velocity distribution over an arbitrary set of complex valued poles. This method works by keeping the integration contour on the real line, and applying a trapezoid rule-like integration scheme over all discretized intervals. In intervals containing a pole, the velocity distribution is linearly interpolated, and the analytic result for the integral over a linear function is used. The integration scheme is validated by comparing its results to the analytic plasma dispersion function for Maxwellian distributions. We then show the utility of this method by computing the Thomson scattering spectra for several non-Maxwellian distributions: the kappa, super-Gaussian, and toroidal distributions. Thomson scattering is a valuable plasma diagnostic tool for both laboratory and space plasmas, but the technique relies on fitting measured wave spectra to a forward model, which typically assumes Maxwellian plasmas. Therefore, this integration method can expand the capabilities of Thomson scatter diagnostics to regimes where the plasma is non-Maxwellian, including high energy density plasmas, frictionally heated plasmas in the ionosphere, and plasmas with a substantial suprathermal electron tail.