The pressure for non-interacting fermions with degeneracy is
where the second form is obtained with an integration by parts. We use units where . The variable substitutions from Johns96 are
,
, and
. (Presumably this choice of variables gives better results for non-relativistic fermions because the mass is separated from the chemical potential in the definition of
, but I haven't checked this). These replacements give
where . Re-expressing in terms of
, one obtains
The pressure is
where .
The Sommerfeld expansion for is
This is an asymptotic expansion, and must thus be used with care. Define . The first term in the Sommerfeld expansion for
depends only on
alone:
where . This expression cannot be used when
is small, but a Taylor series expansion can be used instead. A few terms are
The number density is
Note that because the density is a derivative, it is possible that the terms in the density fail before the terms in the pressure, thus we should use one less term for the density when using the expansion. The entropy is
The derivative of the number density with respect to the chemical potential is
The derivative of the number density with respect to the temperature is
and the derivative of the entropy density with respect to the temperature is
Finally, the derivative of the number density with respect to the mass is more involved because of the mass-dependent prefactor.
These expansions are used in o2scl::fermion_eval_thermo::calc_mu_deg() and o2scl::fermion_deriv_thermo::calc_mu_deg() .
There is a useful identity (Chandrasekhar10 and Tooper69)
which works well when . This result directly gives the sum in Johns96
The function is implemented in GSL as
gsl_sf_bessel_Kn_scaled
. In the case that one wants to include antiparticles, the result is similar
where the scaled Bessel function has been separated out. The density is
and
if antiparticles are included. The entropy density can also be written down similarly.
These expansions are used in o2scl::fermion_eval_thermo::calc_mu_ndeg() .
Documentation generated with Doxygen. Provided under the
GNU Free Documentation License (see License Information).