sarkas.utilities.fdints.fermidirac_integral
sarkas.utilities.fdints.fermidirac_integral#
- sarkas.utilities.fdints.fermidirac_integral(p, eta)[source]#
Wrapper function to calculate the Fermi-Dirac integral using approximation found in Ref. [Fukushima, 2015].
- Parameters
p (float) – Order of the unnormalized Fermi-Dirac integral.
eta (float) – Value of the normalized chemical potential \(\mu/(k_BT)\).
- Return type
- Returns
fd (float) – Value of the Fermi-Dirac integral.
- Raises
SarkasError – If the order of the integral is not supported.
Examples
>>> from numpy import pi, sqrt >>> # Electron density >>> eta = -0.2860 >>> I = fermidirac_integral(p = 0.5, eta = eta) >>> lambda_deB = 1.957093e-11 # [m] >>> ne = 4.0/( sqrt(pi) * lambda_deB**3 ) * I >>> f"{ne:.4e}" '1.6201e+32'