Holtsmark distribution explained

The (one-dimensional) Holtsmark distribution is a continuous probability distribution. The Holtsmark distribution is a special case of a stable distribution with the index of stability or shape parameter

\alpha equal to 3/2 and the skewness parameter \beta of zero. Since \beta equals zero, the distribution is symmetric, and thus an example of a symmetric alpha-stable distribution. The Holtsmark distribution is one of the few examples of a stable distribution for which a closed form expression of the probability density function is known. However, its probability density function is not expressible in terms of elementary functions; rather, the probability density function is expressed in terms of hypergeometric functions.

The Holtsmark distribution has applications in plasma physics and astrophysics. In 1919, Norwegian physicist Johan Peter Holtsmark proposed the distribution as a model for the fluctuating fields in plasma due to the motion of charged particles.^[1] It is also applicable to other types of Coulomb forces, in particular to modeling of gravitating bodies, and thus is important in astrophysics.^{[2] [3]}

Characteristic function

The characteristic function of a symmetric stable distribution is:

\varphi(t;\mu,c) = \exp\left[~it\mu - |c t| ^{\alpha~\right],}

where

\alpha is the shape parameter, or index of stability, \mu is the location parameter, and c is the scale parameter.

Since the Holtsmark distribution has

\alpha =3 / 2 , its characteristic function is:^[4] \varphi(t;\mu,c) = \exp\left[~it\mu - |c t| ^{3/2} ~\right] .

Since the Holtsmark distribution is a stable distribution with,

\mu represents the mean of the distribution.^{[5] [6]} Since, \mu also represents the median and mode of the distribution. And since, the variance of the Holtsmark distribution is infinite.^[5] All higher moments of the distribution are also infinite.^[5] Like other stable distributions (other than the normal distribution), since the variance is infinite the dispersion in the distribution is reflected by the scale parameter, c. An alternate approach to describing the dispersion of the distribution is through fractional moments.^[5]

Probability density function

In general, the probability density function, f(x), of a continuous probability distribution can be derived from its characteristic function by:

f(x) = (1)/(2 \pi) \int ^{inf ty} _{-inf ty} \varphi(t)e ^{-ixt} dt .

Most stable distributions do not have a known closed form expression for their probability density functions. Only the normal, Cauchy and Lévy distributions have known closed form expressions in terms of elementary functions. The Holtsmark distribution is one of two symmetric stable distributions to have a known closed form expression in terms of hypergeometric functions. When

\mu is equal to 0 and the scale parameter is equal to 1, the Holtsmark distribution has the probability density function:\begin{align} f(x; 0 , 1 ) & = { 1 \over \pi } \Gamma\left({ 5 \over 3 }\right) { _{2F} _{3}\left({5} \over 12 }, { 11 \over 12 }; { 1 \over 3 }, { 1 \over 2 }, { 5 \over 6 }; - { 4 x ^{6} \over 729 }\right) \\ & {} {} - { x ^{2} \over 3 \pi } { _{3F} _{4}\left({3} \over 4 }, { 1 }, { 5 \over 4 }; { 2 \over 3 }, { 5 \over 6 }, { 7 \over 6 }, { 4 \over 3 }; - { 4 x ^{6} \over 729 }\right) \\ & {} {} + { 7 x ^{4} \over 81 \pi } \Gamma\left({ 4 \over 3 }\right) { _{2F} _{3}\left({13} \over 12 }, { 19 \over 12 }; { 7 \over 6 }, { 3 \over 2 }, { 5 \over 3 }; - { 4 x ^{6} \over 729 }\right), \end{align}

where

{\Gamma(x)} is the gamma function and _{mF} _{n} is a hypergeometric function.^[7] One has also^[8] \begin{align} f(x; 0 , 1 )& = (- x ^{2})/(6 \pi) \left[~ _{2F} _{2\left(1,} (3)/(2) ; (4)/(3) , (5)/(3) ; - (4 ix ^{3})/(27) \right) + ~ _{2F} _{2\left(1,} (3)/(2) ; (4)/(3) , (5)/(3) ; (4 ix ^{3})/(\right)\right]\\ & + 27) (4)/(\left[ Bi '\left( - 3 x 3 ^{2/3}) (x ^{2})/(\right)\cos\left( 3 x 3 ^{1/3}) (2 x ^{3})/(\right) + 27) (x)/(~ Bi \left( - 3 ^{2/3}) (x ^{2})/(\right)\sin\left( 3 x 3 ^{1/3}) (2 x ^{3})/(27) \right)\right], \end{align}

where

Bi is the Airy function of the second kind and Bi ' its derivative. The arguments of the _{2F} _{2} functions are pure imaginary complex numbers, but the sum of the two functions is real. For x positive, the function Bi ( - x) is related to the Bessel functions of fractional order J _{-1/3} and J _{1/3} and its derivative to the Bessel functions of fractional order J _{-2/3} and J _{2/3} . Therefore, one can write^[8] \begin{align} f(x; 0 , 1 )& = (4 x ^{2})/(27 \sqrt{ 3) }\left\\\&-\frac\left[~_2F_2\left(1, \frac{3}{2}; \frac{4}{3}, \frac{5}{3}; -\frac{4ix^3}{27}\right)+~_2F_2\left(1, \frac{3}{2}; \frac{4}{3}, \frac{5}{3}; \frac{4ix^3}{27}\right)\right].\end

References

• Hummer . D. G. . 10.1016/0022-4073(86)90011-7 . Rational approximations for the holtsmark distribution, its cumulative and derivative . Journal of Quantitative Spectroscopy and Radiative Transfer . 36 . 1–5. 1986 . 1986JQSRT..36....1H .

Notes and References

1. 10.1002/andp.19193630702. 363. 7. 577–630. Holtsmark. J.. Uber die Verbreiterung von Spektrallinien. Annalen der Physik. 1919. 1919AnP...363..577H .
2. 10.1086/144420. 0004-637X. 95. 489. Chandrasekhar. S.. J. von Neumann. The Statistics of the Gravitational Field Arising from a Random Distribution of Stars. I. The Speed of Fluctuations. The Astrophysical Journal. 1942. 1942ApJ....95..489C.
3. 10.1103/RevModPhys.15.1. 15. 1. 1–89. Chandrasekhar. S.. Stochastic Problems in Physics and Astronomy. Reviews of Modern Physics. 1943-01-01. 1943RvMP...15....1C.
4. Book: One-Dimensional Stable Distributions. Zolotarev, V. M.. 1, 41. 1986. Providence, RI. American Mathematical Society. 978-0-8218-4519-6. registration. holtsmark..
5. Book: Basic Properties of Univariate Stable Distributions. Stable Distributions: Models for Heavy Tailed Data. Nolan, J. P.. 3, 15–16. http://lpmt-theory.wdfiles.com/local--files/michael-blog/stablePDF.pdf. 2008. 2011-02-06.
6. Book: Handbook of Heavy Tailed Distributions in Finance. limited. Rachev, S. T.. Modeling Financial Data. 111–112. Nolan, J. P.. 2003. Amsterdam. Elsevier. 978-0-444-50896-6.
7. Book: Continuous and Discrete Properties of Stochastic Processes. Lee, W. H.. 37–39. PhD thesis. 2010. University of Nottingham.
8. 10.1140/epjp/s13360-020-00248-4. 135. 236. Pain. Jean-Christophe. Expression of the Holtsmark function in terms of hypergeometric _{2F} _{2} and Airy Bi functions. Eur. Phys. J. Plus . 2020. 211030564. 2001.11893.