nse_eos Class Reference

Equation of state for nuclei in statistical equilibrium. More...

#include <nse_eos.h>


Detailed Description

Equation of state for nuclei in statistical equilibrium.

This class computes the composition of matter in nuclear statistical equilibrium. The chemical potential of a nucleus X with proton number $ Z_X $ and neutron number $ N_X $ is given by

\[ \mu_X = N \mu_n + Z \mu_p - E_{\mathrm{bind},X} \]

where $ \mu_n $ and $ \mu_p $ are the neutron and proton chemical potentials and $ E_{\mathrm{bind},X} $ is the binding energy of the nucleus.

The baryon number density and electron fraction are then given by

\[ n_B = n_{X} (N_X + Z_X) \qquad Y_e n_B = n_X Z_X \]

where $ n_X $ is the number density which is determined from the chemical potential above.

This implicitly assumes that the nuclei are non-interacting.

Idea for future:
Right now calc_density() needs a very good guess. This could be fixed, probably by solving for the log(mu/T) instead of mu.

Definition at line 62 of file nse_eos.h.


Public Member Functions

int calc_mu (double mun, double mup, double T, double &nb, double &Ye, thermo &th, nuclear_dist &nd)
 Calculate the equation of state as a function of the chemical potentials.
int calc_density (double nb, double Ye, double T, double &mun, double &mup, thermo &th, nuclear_dist &nd)
 Calculate the equation of state as a function of the densities.
int set_mroot (mroot< solve_parms, mm_funct< solve_parms > > &rp)
 Set the solver for use in computing the chemical potentials.

Data Fields

gsl_mroot_hybrids< solve_parms,
mm_funct< solve_parms > > 
def_root
 Default solver.

Member Function Documentation

int calc_density ( double  nb,
double  Ye,
double  T,
double &  mun,
double &  mup,
thermo th,
nuclear_dist nd 
)

Calculate the equation of state as a function of the densities.

Given the baryon number density nb, and the electron fraction Ye and the temperature T, this computes the composition (the individual densities are stored in the distribution nd) and the chemical potentials are given in mun and mup .

This function uses the solver to self-consistently compute the chemical potentials.

int calc_mu ( double  mun,
double  mup,
double  T,
double &  nb,
double &  Ye,
thermo th,
nuclear_dist nd 
)

Calculate the equation of state as a function of the chemical potentials.

Given mun, mup and T, this computes the composition (the individual densities are stored in the distribution nd) the baryon number density nb, and the electron fraction Ye.

This function does not use the solver.


The documentation for this class was generated from the following file:

Documentation generated with Doxygen and provided under the GNU Free Documentation License. See License Information for details.

Project hosting provided by SourceForge.net Logo, O2scl Sourceforge Project Page