Class eos_leptons (o2scl)

O2scl : Class List

class eos_leptons

Lepton and photon EOS.

Public Types

typedef boost::numeric::ublas::vector<double> ubvector

Public Functions

eos_leptons()
int pair_mu(double T)

Thermodynamics from the electron and muon chemical potentials.

The temperature should be in units of \( 1/\mathrm{fm} \) .

int pair_mu_eq(double T)

Thermodynamics from the electron chemical potential in weak equilibrium.

When include_muons is false, this function is essentially equivalent to pair_mu().

The temperature should be in units of \( 1/\mathrm{fm} \) .

inline void high_acc()
inline void standard_acc()
int pair_density(double T)

Thermodynamics from the electron and muon densities.

The temperature should be in units of \( 1/\mathrm{fm} \) .

The current values of the electron and muon chemical potentials are used as initial guesses to fermion_rel::pair_density() for both the electrons and (if included) muons.

int pair_density_eq(double nq, double T)

Thermodynamics from the charge density in weak equilibrium.

The first argument nq, is the total negative charge density including electrons (and muons if include_muons is true) and T is the temperature.

When include_muons is false, this function is essentially equivalent to pair_density() using nq for the electron density.

The charge density should be in units of \( 1/\mathrm{fm}^{-3} \) and the temperature should be in units of \( 1/\mathrm{fm} \).

The current values of the electron chemical potential potentials is used as initial guess. If pde_from_density is true, then the current value of the electron density is also used as an initial guess.

Public Members

mroot_hybrids mh

Solver for pair_density_eq()

fermion_rel frel

Relativistic fermion thermodyanmics.

fermion e

Electron.

bool pde_from_density

If true, use the electron density for pair_density_eq_fun().

fermion mu

Muon.

boson ph

Photon.

bool err_nonconv

If true, call the error handler if msolve() does not converge (default true)

bool include_muons

If true, include muons (default true)

bool include_photons

If true, include photons (default false)

thermo th

Thermodynamic quantities.

int verbose

Protected Functions

int electron_density(double T)

Electron thermodynamics from the electron density.

The temperature should be in units of \( 1/\mathrm{fm} \) .

Because this function uses fermion_rel::pair_density(), the current electron chemical potential is used as an initial guess.

Note

This internal function is designed to be used when fermion_rel::err_nonconv is false .

Note

If this function fails, it does not necessarily call the error handler. This allows it to be used inside pair_density_eq_fun() and is the reason why this function is protected instead of public.

int pair_density_eq_fun(size_t nv, const ubvector &x, ubvector &y, double T, double nq)

Function to solve for pair_density_eq()

The temperature should be in units of \( 1/\mathrm{fm} \) .

Note

This internal function presumes that include_muons is true (otherwise there’s nothing to solve) and that fermion_rel::err_nonconv is false .