Two-dimensional PIC simulation of collective Thomson scattering in a beam-plasma system

Collective Thomson scattering (CTS) in a beam-plasma system is reproduced by using 2D PIC simulations and the characteristics of the scattered wave spectrum are examined. By formulating the geometric shape of the scattered wave spectrum in wave number space, where the velocity vector of the beam component and the wave vectors of the incident and scattered waves are arbitrary, it is demonstrated that the spectrum in 2D wave number space becomes asymmetric. The spectrum of scattered waves propagating in a specific direction is presented as a function of wavelength to show that the electron (ion) feature is amplified and becomes asymmetric or distorted when Buneman (ion acoustic) instability occurs. An additional simulation is conducted for a weak, linearly stable beam-plasma system with a hot beam, and confirmed that the obtained scattered wave spectrum shows asymmetric feature. The results are expected to be applicable to the interpretation of radar observations of ionospheric plasmas as well as CTS measurements in laboratory plasmas.

💡 Research Summary

This paper presents a comprehensive investigation of Collective Thomson Scattering (CTS) in non-equilibrium beam-plasma systems using two-dimensional, self-consistent electromagnetic Particle-in-Cell (PIC) simulations. The primary goal is to understand how beam-driven instabilities and non-Maxwellian velocity distributions manifest in the multidimensional spectrum of scattered electromagnetic waves, with implications for diagnosing laboratory and space plasmas.

The study employs a 2D periodic simulation domain with three particle species: background electrons, background ions, and beam ions. A finite packet of monochromatic, linearly polarized electromagnetic wave is injected into the system to act as the incident wave for CTS. The research compares three primary scenarios: a thermal equilibrium plasma without a beam (Run 1), a system with a strong, warm beam that excites the Buneman instability (Run 2), and a system with a weak, cold beam that excites the ion acoustic instability (Run 3). An additional simulation models a linearly stable, weak beam system with a hot beam component.

The key findings are multi-faceted. First, the study successfully reproduces the well-known CTS spectrum for thermal plasma, showing symmetric “ion” and “electron” features corresponding to scattering off ion acoustic and Langmuir waves, respectively. The novel contribution lies in the analysis of the scattered wave spectrum in 2D wavenumber space (kx-ky). The authors derive geometric formulations showing that the resonance condition involving a beam mode (with dispersion ω = k·ub) leads to an asymmetric circle equation in (kx, ky) space, explaining the distorted, non-concentric ring structures observed in the simulations when a beam is present.

Second, by examining the spectrum of waves scattered into a specific direction (e.g., the negative y-direction), the impact of instabilities becomes starkly clear. In the case of Buneman instability (Run 2), the high-frequency “electron feature” peak on the short-wavelength side is dramatically amplified and asymmetric. In the case of ion acoustic instability (Run 3), the low-frequency “ion feature” is amplified, also exhibiting a pronounced asymmetry favoring the short-wavelength side. This selective amplification is directly linked to the growth of unstable wave modes in a specific direction dictated by the beam velocity.

Third, and crucially, the simulation of a stable, weak beam-plasma system (with a negative slope everywhere in the ion distribution function) still yields an asymmetric ion feature in the CTS spectrum. This demonstrates that asymmetry in the scattered spectrum is not solely a signature of an ongoing instability but can also be a direct consequence of a non-thermal, anisotropic velocity distribution itself. The beam’s presence modifies the plasma wave fluctuations, which in turn imprints its signature onto the CTS spectrum.

In conclusion, this work provides a significant advancement in simulating CTS self-consistently in multidimensional, non-equilibrium plasmas. It establishes a clear link between beam parameters (density, temperature, drift velocity and direction), the resulting plasma wave spectra (whether unstable or stable), and the observable characteristics of the CTS spectrum. The results offer a powerful interpretive framework for deciphering complex, often asymmetric, CTS spectra observed in laboratory experiments (e.g., high-power laser-plasma interactions) and in natural settings like the ionosphere, where non-thermal ion beams are common. The 2D PIC approach proves to be a robust tool for isolating and understanding the geometric and kinetic origins of spectral distortions in collective scattering diagnostics.