An official website of the United States government
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.
more »
« less
Thomson scattering diagnostics of nonthermal plasma from particle-in-cell simulations
Optical Thomson scattering is now a mature diagnostic tool for precisely measuring local plasma density and temperature. These measurements typically take advantage of a simplified analytical model of the scattered spectrum, which is built upon the assumption that each plasma species is in thermal equilibrium. However, this assumption fails for most laboratory plasmas of interest, which are often produced through high field ionization of atoms via ultrashort laser pulses and vulnerable to several kinetic instabilities. While it is possible to analytically model the Thomson scattered spectrum for some non-Maxwellian distribution functions, it is often not practical to do so for laboratory plasmas with highly complex and unstable distribution functions. We present a new method for predicting the Thomson scattered spectrum from any plasma directly from fully kinetic particle-in-cell simulations. This approach allows us to model the contributions of kinetic instabilities to the Thomson spectrum that aren’t taken into account in Maxwellian theory. We demonstrate this method’s capability to capture nonthermal features in the Thomson spectrum by simulating a simple bumpon- tail plasma as well as a more complex laser-ionized plasma. The versatility of this approach makes it an effective aid in the experimental design of Thomson diagnostics to directly characterize kinetic instabilities in laboratory plasmas. Index Terms—plasma measurement, low-temperature plasmas, plasma diagnostics, plasma simulation, plasma stability, plasma density, plasma temperature
more »
« less
- Award ID(s):
- 2003354
- PAR ID:
- 10462344
- Date Published:
- Journal Name:
- Advanced Accelerator Concepts Workshop 2022
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)We investigate kinetic entropy-based measures of the non-Maxwellianity of distribution functions in plasmas, i.e. entropy-based measures of the departure of a local distribution function from an associated Maxwellian distribution function with the same density, bulk flow and temperature as the local distribution. First, we consider a form previously employed by Kaufmann & Paterson ( J. Geophys. Res. , vol. 114, 2009, A00D04), assessing its properties and deriving equivalent forms. To provide a quantitative understanding of it, we derive analytical expressions for three common non-Maxwellian plasma distribution functions. We show that there are undesirable features of this non-Maxwellianity measure including that it can diverge in various physical limits and elucidate the reason for the divergence. We then introduce a new kinetic entropy-based non-Maxwellianity measure based on the velocity-space kinetic entropy density, which has a meaningful physical interpretation and does not diverge. We use collisionless particle-in-cell simulations of two-dimensional anti-parallel magnetic reconnection to assess the kinetic entropy-based non-Maxwellianity measures. We show that regions of non-zero non-Maxwellianity are linked to kinetic processes occurring during magnetic reconnection. We also show the simulated non-Maxwellianity agrees reasonably well with predictions for distributions resembling those calculated analytically. These results can be important for applications, as non-Maxwellianity can be used to identify regions of kinetic-scale physics or increased dissipation in plasmas.more » « less
-
Understanding how plasmas thermalize when density gradients are steep remains a fundamental challenge in plasma physics, with direct implications for fusion experiments and astrophysical phenomena. Standard hydrodynamic models break down in these regimes, and kinetic theories make predictions that have never been directly tested. Here, we present the first detailed phase-space measurements of a strongly coupled plasma as it evolves from sharp density gradients to thermal equilibrium. Using laser-induced fluorescence imaging of an ultracold calcium plasma, we track the complete ion distribution function f(x,v,t). We discover that commonly used kinetic models (Bhatnagar–Gross–Krook and Lenard–Bernstein) overpredict thermalization rates, even while correctly capturing the initial counterstreaming plasma formation. Our measurements reveal that the initial ion acceleration response scales linearly with electron temperature, and that the simulations underpredict the initial ion response. In our geometry we demonstrate the formation of well-controlled counterpropagating plasma beams. This experimental platform enables precision tests of kinetic theories and opens new possibilities for studying plasma stopping power and flow-induced instabilities in strongly coupled systems.more » « less
-
null (Ed.)ABSTRACT Various plasma waves and instabilities are abundantly present in the solar wind plasma, as evidenced by spacecraft observations. Among these, propagating modes and instabilities driven by temperature anisotropies are known to play a significant role in the solar wind dynamics. In situ measurements reveal that the threshold conditions for these instabilities adequately explain the solar wind conditions at large heliocentric distances. This paper pays attention to the combined effects of electron firehose instability driven by excessive parallel electron temperature anisotropy (T⊥e < T∥e) at high beta conditions, and electromagnetic ion cyclotron instability driven by excessive perpendicular proton temperature anisotropy (T⊥i > T∥i). By employing quasilinear kinetic theory based upon the assumption of bi-Maxwellian velocity distribution functions for protons and electrons, the dynamical evolution of the combined instabilities and their mutual interactions mediated by the particles is explored in depth. It is found that while in some cases, the two unstable modes are excited and saturated at distinct spatial and temporal scales, in other cases, the two unstable modes are intermingled such that a straightforward interpretation is not so easy. This shows that when the dynamics of protons and electrons are mutually coupled and when multiple unstable modes are excited in the system, the dynamical consequences can be quite complex.more » « less
-
Abstract The quasi-steady states of collisionless plasmas in space (e.g., in the solar wind and planetary environments) are governed by the interactions of charged particles with wave fluctuations. These interactions are responsible not only for the dissipation of plasma waves but also for their excitation. The present analysis focuses on two instabilities, mirror and electromagnetic ion cyclotron instabilities, associated with the same proton temperature anisotropyT⊥>T∥(where ⊥, ∥ are directions defined with respect to the local magnetic field vector). Theories relying on standard Maxwellian models fail to link these two instabilities (i.e., predicted thresholds) to the proton quasi-stable anisotropies measured in situ in a completely satisfactory manner. Here we revisit these instabilities by modeling protons with the generalized bi-Kappa (bi-κpower-law) distribution, and by a comparative analysis of a 2D hybrid simulation with the velocity-moment-based quasi-linear (QL) theory. It is shown that the two methods feature qualitative and, even to some extent, quantitative agreement. The reduced QL analysis based upon the assumption of a time-dependent bi-Kappa model thus becomes a valuable theoretical approach that can be incorporated into the present studies of solar wind dynamics.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Conference Paper:
- The DOI is not currently available.
