spectral_density¶

plasmapy.diagnostics.thomson.
spectral_density
(wavelengths: Unit("nm"), probe_wavelength: Unit("nm"), n: Unit("1 / m3"), Te: Unit("K"), Ti: Unit("K"), efract: numpy.ndarray = None, ifract: numpy.ndarray = None, ion_species: Union[str, List[str], plasmapy.particles.particle_class.Particle, List[plasmapy.particles.particle_class.Particle]] = 'H+', electron_vel: Unit("m / s") = None, ion_vel: Unit("m / s") = None, probe_vec=array([1, 0, 0]), scatter_vec=array([0, 1, 0])) → Tuple[Union[numpy.floating, numpy.ndarray], numpy.ndarray]¶ Calculate the spectral density function for Thomson scattering of a probe laser beam by a multispecies Maxwellian plasma.
This function calculates the spectral density function for Thomson scattering of a probe laser beam by a plasma consisting of one or more ion species and a one or more thermal electron populations (the entire plasma is assumed to be quasineutral)
\[S(k,\omega) = \sum_e \frac{2\pi}{k} \bigg 1  \frac{\chi_e}{\epsilon} \bigg ^2 f_{e0,e} \bigg (\frac{\omega}{k} \bigg ) + \sum_i \frac{2\pi Z_i}{k} \bigg \frac{\chi_e}{\epsilon} \bigg ^2 f_{i0,i} \bigg ( \frac{\omega}{k} \bigg )\]where \(\chi_e\) is the electron component susceptibility of the plasma and \(\epsilon = 1 + \sum_e \chi_e + \sum_i \chi_i\) is the total plasma dielectric function (with \(\chi_i\) being the ion component of the susceptibility), \(Z_i\) is the charge of each ion, \(k\) is the scattering wavenumber, \(\omega\) is the scattering frequency, and \(f_{e0,e}\) and \(f_{i0,i}\) are the electron and ion velocity distribution functions respectively. In this function the electron and ion velocity distribution functions are assumed to be Maxwellian, making this function equivalent to Eq. 3.4.6 in Sheffield.
Parameters:  wavelengths (
Quantity
) – Array of wavelengths over which the spectral density function will be calculated. (convertible to nm)  probe_wavelength (
Quantity
) – Wavelength of the probe laser. (convertible to nm)  n (
Quantity
) – Mean (0th order) density of all plasma components combined. (convertible to cm^3.)  Te (
Quantity
, shape (Ne, )) – Temperature of each electron component. Shape (Ne, ) must be equal to the number of electron components Ne. (in K or convertible to eV)  Ti (
Quantity
, shape (Ni, )) – Temperature of each ion component. Shape (Ni, ) must be equal to the number of ion components Ni. (in K or convertible to eV)  efract (array_like, shape (Ne, ), optional) – An arraylike object where each element represents the fraction (or ratio) of the electron component number density to the total electron number density. Must sum to 1.0. Default is a single electron component.
 ifract (array_like, shape (Ni, ), optional) – An arraylike object where each element represents the fraction (or ratio) of the ion component number density to the total ion number density. Must sum to 1.0. Default is a single ion species.
 ion_species (str or
Particle
, shape (Ni, ), optional) – A list or single instance ofParticle
, or strings convertible toParticle
. Default is'H+'
corresponding to a single species of hydrogen ions.  electron_vel (
Quantity
, shape (Ne, 3), optional) – Velocity of each electron component in the rest frame. (convertible to m/s) Defaults to a stationary plasma [0, 0, 0] m/s.  ion_vel (
Quantity
, shape (Ni, 3), optional) – Velocity vectors for each electron population in the rest frame (convertible to m/s) Defaults zero drift for all specified ion species.  probe_vec (float
ndarray
, shape (3, )) – Unit vector in the direction of the probe laser. Defaults to [1, 0, 0].  scatter_vec (float
ndarray
, shape (3, )) – Unit vector pointing from the scattering volume to the detector. Defaults to [0, 1, 0] which, along with the defaultprobe_vec
, corresponds to a 90 degree scattering angle geometry.
Returns:  alpha (float) – Mean scattering parameter, where
alpha
> 1 corresponds to collective scattering andalpha
< 1 indicates noncollective scattering. The scattering parameter is calculated based on the total plasma density n.  Skw (
Quantity
) – Computed spectral density function over the inputwavelengths
array with units of s/rad.
Notes
For details, see “Plasma Scattering of Electromagnetic Radiation” by Sheffield et al. ISBN 9780123748775. This code is a modified version of the program described therein.
For a concise summary of the relevant physics, see Chapter 5 of Derek Schaeffer’s thesis, DOI: 10.5281/zenodo.3766933.
 wavelengths (