chemical_potential¶

plasmapy.formulary.
chemical_potential
(n_e: Unit("1 / m3"), T: Unit("K")) > Unit(dimensionless)¶ Calculate the ideal chemical potential.
Parameters: Returns: beta_mu – The dimensionless ideal chemical potential. That is the ratio of the ideal chemical potential to the thermal energy.
Return type: Raises: TypeError
– If argument is not aQuantity
.UnitConversionError
– If argument is in incorrect units.ValueError
– If argument contains invalid values.
Warns: ~astropy.units.UnitsWarning – If units are not provided, SI units are assumed.
Notes
The ideal chemical potential is given by [1]:
\[\chi_a = I_{1/2}(\beta \mu_a^{ideal})\]where \(\chi\) is the degeneracy parameter, \(I_{1/2}\) is the Fermi integral with order 1/2, \(\beta\) is the inverse thermal energy \(\beta = 1/(k_B T)\), and \(\mu_a^{ideal}\) is the ideal chemical potential.
The definition for the ideal chemical potential is implicit, so it must be obtained numerically by solving for the Fermi integral for values of chemical potential approaching the degeneracy parameter. Since values returned from the Fermi_integral are complex, a nonlinear LevenbergMarquardt least squares method is used to iteratively approach a value of \(\mu\) which minimizes \(I_{1/2}(\beta \mu_a^{ideal})  \chi_a\)
This function returns \(\beta \mu^{ideal}\) the dimensionless ideal chemical potential.
Warning: at present this function is limited to relatively small arguments due to limitations in the
mpmath
package’s implementation ofpolylog
, which PlasmaPy uses in calculating the Fermi integral.References
[1] Bonitz, Michael. Quantum kinetic theory. Stuttgart: Teubner, 1998. Example
>>> from astropy import units as u >>> chemical_potential(n_e=1e21*u.cm**3,T=11000*u.K) # doctest: +SKIP <Quantity 2.00039985e12>