Coulomb_logarithm¶

plasmapy.formulary.collisions.
Coulomb_logarithm
(T: Unit("K"), n_e: Unit("1 / m3"), species: (<class 'plasmapy.particles.particle_class.Particle'>, <class 'plasmapy.particles.particle_class.Particle'>), z_mean: Unit(dimensionless) = <Quantity nan>, V: Unit("m / s") = <Quantity nan m / s>, method='classical')¶ Computes the Coulomb logarithm.
 Parameters
T (Quantity) – Temperature in units of temperature or energy per particle, which is assumed to be equal for both the test particle and the target particle.
n_e (Quantity) – The electron number density in units convertible to per cubic meter.
species (tuple) – A tuple containing string representations of the test particle (listed first) and the target particle (listed second).
z_mean (Quantity, optional) – The average ionization (arithmetic mean) of a plasma for which a macroscopic description is valid. This parameter is used to compute the average ion density (given the average ionization and electron density) for calculating the ion sphere radius for nonclassical impact parameters.
z_mean
is a required parameter ifmethod
is"ls_full_interp"
,"hls_max_interp"
, or"hls_full_interp"
.V (Quantity, optional) – The relative velocity between particles. If not provided, thermal velocity is assumed: \(\mu V^2 \sim 2 k_B T\) where
mu
is the reduced mass.method (str, optional) – The method by which to compute the Coulomb logarithm. The default method is the classical straightline LandauSpitzer method (
"classical"
or"ls"
). The other 6 supported methods are"ls_min_interp"
,"ls_full_interp"
,"ls_clamp_mininterp"
,"hls_min_interp"
,"hls_max_interp"
, and"hls_full_interp"
. Please refer to the Notes section of this docstring for more information, including about abbreviated aliases of these names.
 Returns
ln_Lambda – The dimensionless Coulomb logarithm.
 Return type
 Raises
ValueError – If the mass or charge of either particle cannot be found, or any of the inputs contain incorrect values.
UnitConversionError – If the units on any of the inputs are incorrect. If the n_e, T, or V are not Quantities.
PhysicsError – If the result is smaller than 1.
RelativityError – If the input velocity is same or greater than the speed of light.
 Warns
UnitsWarning
– If units are not provided, SI units are assumed.RelativityWarning
– If the input velocity is greater than 5% of the speed of light.
Notes
Summary of Supported Methods of Computing the Coulomb Logarithm
PlasmaPy supports 7 methods of computing the Coulomb logarithm:
"classical"
or"ls"
"ls_min_interp"
or"GMS1"
"ls_full_interp"
or"GMS2"
"ls_clamp_mininterp"
or"GMS3"
"hls_min_interp"
or"GMS4"
"hls_max_interp"
or"GMS5"
"hls_full_interp"
or"GMS6"
Options 1–4 are straightline LandauSpitzer (
"ls..."
) methods in which the trajectory of a Coulomb collision is modeled as a straight line. For the straightline LandauSpitzer methods, the Coulomb logarithm (\(\ln{\Lambda}\)) is defined to be:\[\ln{\Lambda} \equiv \ln\left( \frac{b_{max}}{b_{min}} \right)\]Options 5–7 are hyperbolic LandauSpitzer (
"hls..."
) methods in which the trajectory of a Coulomb collision is modeled as a hyperbola. For the hyperbolic LandauSpitzer methods, the Coulomb logarithm (\(\ln{\Lambda}\)) is defined to be:\[\ln{\Lambda} \equiv \frac{1}{2} \ln\left(1 + \frac{b_{max}^2}{b_{min}^2} \right)\]For all 7 methods, \(b_{min}\) and \(b_{max}\) are the inner impact parameter and the outer impact parameter, respectively, for Coulomb collisions 1; \(b_{min}\) and \(b_{max}\) are each computed by
impact_parameter
, another function.The abbreviations of Options 27 (
"GMS..."
) refer to the first initials of the three authors of Reference 4.Note
For stronglycoupled plasma, PlasmaPy recommends Option 7,
"hls_full_interp"
or"GMS6"
, because of its high accuracy regardless of a plasma’s strength of coupling.Explanation of Supported Methods of Computing the Coulomb Logarithm
In this section, further information about each method, such as about interpolation and other special features, is documented. Please refer to Reference 1 and Reference 4 for additional information about these methods.
 Option 1:
"classical"
or"ls"
(LandauSpitzer) The classical straightline LandauSpitzer method in which \(b_{min}\) is defined to be the higher of the de Broglie wavelength (\(\lambda_{de Broglie}\)) and the distance of closest approach (\(\rho_{\perp}\)) if they are not equal (and either of the two if they are equal) and \(b_{max}\) is defined to be the Debye length (\(\lambda_{Debye}\)).
\[\ln{\Lambda} \equiv \ln\left( \frac{b_{max}}{b_{min}} \right)\]\[\begin{split}b_{min} \equiv \left\{ \begin{array}{ll} \lambda_{de Broglie} & \mbox{if } \lambda_{de Broglie} \geq \rho_{\perp} \\ \rho_{\perp} & \mbox{if } \rho_{\perp} \geq \lambda_{de Broglie} \end{array} \right.\end{split}\]\[b_{max} \equiv \lambda_{Debye}\]The inner impact parameter (\(b_{min}\)) is the higher of \(\lambda_{de Broglie}\) and \(\rho_{\perp}\) because for impact parameters lower than \(\lambda_{de Broglie}\), quantum effects cause the collision to be nonCoulombic 2 3.
The outer impact parameter (\(b_{max}\)) is defined to be the Debye length (\(\lambda_{Debye}\)) because at distances higher than the Debye length, the electric fields created by other particles are screened out by the electrons rearranging themselves.
The uncertainty of the classical straightline LandauSpitzer method is on the order of its reciprocal.
This method is invalid if \(\ln{\Lambda} < 2\) because of the uncertainty of this method and is invalid if \(\ln{\Lambda} < 0\), which may be true if the coupling parameter is high (such as for nonideal, dense, cold plasmas).
Please refer to Reference 1 for additional information about this method.
 Option 2:
"ls_min_interp"
or"GMS1"
(LandauSpitzer, interpolation of \(b_{min}\)) A straightline LandauSpitzer method in which \(b_{min}\) is interpolated between the de Broglie wavelength (\(\lambda_{de Broglie}\)) and the distance of closest approach (\(\rho_{\perp}\)) and \(b_{max}\) is defined to be the Debye length (\(\lambda_{Debye}\)).
\[\ln{\Lambda} \equiv \ln\left( \frac{b_{max}}{b_{min}} \right)\]\[b_{min} \equiv \sqrt{\lambda_{de Broglie}^2 + \rho_{\perp}^2}\]\[b_{max} \equiv \lambda_{Debye}\]This method is invalid if \(\ln{\Lambda} < 0\), which may be true if the coupling parameter is high (such as for nonideal, dense, cold plasmas).
Note: This is the first method in Table 1 of Reference 4.
 Option 3:
"ls_full_interp"
or"GMS2"
(LandauSpitzer, interpolation of \(b_{min}\) and \(b_{max}\)) A straightline LandauSpitzer method in which \(b_{min}\) and \(b_{max}\) are each interpolated. \(b_{min}\) is interpolated between the de Broglie wavelength (\(\lambda_{de Broglie}\)) and the distance of closest approach (\(\rho_{\perp}\)). \(b_{max}\) is interpolated between the Debye length (\(\lambda_{Debye}\)) and the ion sphere radius (\(a_i\)).
\[\ln{\Lambda} \equiv \ln\left( \frac{b_{max}}{b_{min}} \right)\]\[b_{min} \equiv \sqrt{\lambda_{de Broglie}^2 + \rho_{\perp}^2}\]\[b_{max} \equiv \sqrt{\lambda_{Debye}^2 + a_i^2}\]This method is invalid if \(\ln{\Lambda} < 0\), which may be true if the coupling parameter is high (such as for nonideal, dense, cold plasmas).
Note: This is the second method in Table 1 of Reference 4.
 Option 4:
"ls_clamp_mininterp"
or"GMS3"
(LandauSpitzer with a clamp, interpolation of \(b_{min}\)) A straightline LandauSpitzer method in which the value of \(\ln{\Lambda}\) is clamped at a minimum of \(2\) so that it is impossible for \(\ln{\Lambda} < 0\) (which is possible by the classical LandauSpitzer method). \(b_{min}\) is interpolated between the de Broglie wavelength (\(\lambda_{de Broglie}\)) and the distance of closest approach (\(\rho_{\perp}\)). \(b_{max}\) is defined to be the Debye length (\(\lambda_{Debye}\)).
\[\begin{split}\ln{\Lambda} \equiv \left\{ \begin{array}{ll} \ln\left( \frac{b_{max}}{b_{min}} \right) & \mbox{if } \ln\left( \frac{b_{max}}{b_{min}} \right) \geq 2 \\ 2 & \mbox{if } \ln\left( \frac{b_{max}}{b_{min}} \right) \leq 2 \end{array} \right.\end{split}\]\[b_{min} \equiv \sqrt{\lambda_{de Broglie}^2 + \rho_{\perp}^2}\]\[b_{max} \equiv \lambda_{Debye}\]This method is valid for any plasma because it is impossible for \(\ln{\Lambda} < 0\) by this method, even if the coupling parameter is high (such as for nonideal, dense, cold plasmas).
Note: This is the third method in Table 1 of Reference 4.
 Option 5:
"hls_min_interp"
or"GMS4"
(Hyperbolic LandauSpitzer, interpolation of \(b_{min}\)) A hyperbolic LandauSpitzer method in which \(b_{min}\) is interpolated between the de Broglie wavelength (\(\lambda_{de Broglie}\)) and the distance of closest approach (\(\rho_{\perp}\)) and \(b_{max}\) is defined to be the Debye length (\(\lambda_{Debye}\)).
\[\ln{\Lambda} \equiv \frac{1}{2} \ln\left(1 + \frac{b_{max}^2}{b_{min}^2} \right)\]\[b_{min} \equiv \sqrt{\lambda_{de Broglie}^2 + \rho_{\perp}^2}\]\[b_{max} \equiv \lambda_{Debye}\]This method is valid for any plasma because it is impossible for \(\ln{\Lambda} < 0\) by this method, even if the coupling parameter is high (such as for nonideal, dense, cold plasmas).
Note: This is the fourth method in Table 1 of Reference 4.
 Option 6:
"hls_max_interp"
or"GMS5"
(Hyperbolic LandauSpitzer, interpolation of \(b_{max}\)) A hyperbolic LandauSpitzer method in which \(b_{max}\) is interpolated between the Debye length (\(\lambda_{Debye}\)) and the ion sphere radius (\(a_i\)) and \(b_{min}\) is defined to be the distance of closest approach (\(\rho_{\perp}\)).
\[\ln{\Lambda} \equiv \frac{1}{2} \ln\left(1 + \frac{b_{max}^2}{b_{min}^2} \right)\]\[b_{min} \equiv \rho_{\perp}\]\[b_{max} \equiv \sqrt{\lambda_{Debye}^2 + a_i^2}\]This method is valid for any plasma because it is impossible for \(\ln{\Lambda} < 0\) by this method, even if the coupling parameter is high (such as for nonideal, dense, cold plasmas).
Caution: This method overestimates \(\ln{\Lambda}\) at high temperatures.
Note: This is the fifth method in Table 1 of Reference 4.
 Option 7:
"hls_full_interp"
or"GMS6"
(Hyperbolic LandauSpitzer, interpolation of \(b_{min}\) and \(b_{max}\)) A hyperbolic LandauSpitzer method in which \(b_{min}\) and \(b_{max}\) are each interpolated. \(b_{min}\) is interpolated between the de Broglie wavelength (\(\lambda_{de Broglie}\)) and the distance of closest approach (\(\rho_{\perp}\)). \(b_{max}\) is interpolated between the Debye length (\(\lambda_{Debye}\)) and the ion sphere radius (\(a_i\)).
\[\ln{\Lambda} \equiv \frac{1}{2} \ln\left(1 + \frac{b_{max}^2}{b_{min}^2} \right)\]\[b_{min} \equiv \sqrt{\lambda_{de Broglie}^2 + \rho_{\perp}^2}\]\[b_{max} \equiv \sqrt{\lambda_{Debye}^2 + a_i^2}\]This method is valid for any plasma because it is impossible for \(\ln{\Lambda} < 0\) by this method, even if the coupling parameter is high (such as for nonideal, dense, cold plasmas).
Note: This is the sixth method in Table 1 of Reference 4.
Examples
>>> from astropy import units as u >>> n = 1e19 * u.m**3 >>> T = 1e6 * u.K >>> Coulomb_logarithm(T, n, ('e', 'p+')) 14.545527... >>> Coulomb_logarithm(T, n, ('e', 'p+'), V = 1e6 * u.m / u.s) 11.363478...
References
 1(1,2,3)
Physics of Fully Ionized Gases, L. Spitzer (1962)
 2
Francis, F. Chen. Introduction to plasma physics and controlled fusion 3rd edition. Ch 5 (Springer 2015).
 3
Comparison of Coulomb Collision Rates in the Plasma Physics and Magnetically Confined Fusion Literature, W. Fundamenski and O.E. Garcia, EFDA–JET–R(07)01 (http://www.eurofusionscipub.org/wpcontent/uploads/2014/11/EFDR07001.pdf)
 4(1,2,3,4,5,6,7,8)
Dense plasma temperature equilibration in the binary collision approximation. D. O. Gericke et. al. PRE, 65, 036418 (2002). DOI: 10.1103/PhysRevE.65.036418
See also
impact_parameter
Computes \(b_{min}\) and \(b_{max}\).