Maxwellian_1D

plasmapy.physics.distribution.Maxwellian_1D(v, T, particle='e', v_drift=0, vTh=nan, units='units')

Probability distribution function of velocity for a Maxwellian distribution in 1D.

Returns the probability density function at the velocity v in m/s to find a particle particle in a plasma of temperature T following the Maxwellian distribution function.

Parameters:
  • v (Quantity) – The velocity in units convertible to m/s.
  • T (Quantity) – The temperature in Kelvin.
  • particle (str, optional) – Representation of the particle species(e.g., 'p' for protons, 'D+' for deuterium, or 'He-4 +1' for \(He_4^{+1}\) (singly ionized helium-4)), which defaults to electrons.
  • v_drift (Quantity, optional) – The drift velocity in units convertible to m/s.
  • vTh (Quantity, optional) – Thermal velocity (most probable velocity) in m/s. This is used for optimization purposes to avoid re-calculating vTh, for example when integrating over velocity-space.
  • units (str, optional) – Selects whether to run function with units and unit checks (when equal to “units”) or to run as unitless (when equal to “unitless”). The unitless version is substantially faster for intensive computations.
Returns:

f – Probability density in units of Velocity^-1, normalized so that \(\int_{-\infty}^{+\infty} f(v) dv = 1\).

Return type:

Quantity

Raises:
  • TypeError – The parameter arguments are not Quantities and cannot be converted into Quantities.
  • UnitConversionError – If the parameters are not in appropriate units.
  • ValueError – If the temperature is negative, or the particle mass or charge state cannot be found.

Notes

In one dimension, the Maxwellian distribution function for a particle of mass m, velocity v, a drift velocity V and with temperature T is:

\[f = \sqrt{\frac{m}{2 \pi k_B T}} e^{-\frac{m}{2 k_B T} (v-V)^2} \equiv \frac{1}{\sqrt{\pi v_{Th}^2}} e^{-(v - v_{drift})^2 / v_{Th}^2}\]

where \(v_{Th} = \sqrt{2 k_B T / m}\) is the thermal speed

Examples

>>> from plasmapy.physics import Maxwellian_1D
>>> from astropy import units as u
>>> v=1*u.m/u.s
>>> Maxwellian_1D(v=v, T=30000 * u.K, particle='e', v_drift=0 * u.m / u.s)
<Quantity 5.91632969e-07 s / m>