pynucastro.neutrino_cooling.sneut5_mod module

pynucastro.neutrino_cooling.sneut5_mod module#

Classes and methods that implement the Itoh et al. 1996 thermal neutron loss terms.

class pynucastro.neutrino_cooling.sneut5_mod.NeutrinoComponents[source]#

Bases: object

A simple container that holds the individual components to the neutrino cooling.

splas#

Plasma neutrino cooling rate

Type:

float

spair#

Pair neutrino cooling rate

Type:

float

sphot#

Photo-ionization neutrino cooling rate

Type:

float

sbrem#

Bremmstrahlung neutrino cooling rate

Type:

float

sreco#

Recombination neutrino cooling rate

Type:

float

pynucastro.neutrino_cooling.sneut5_mod.ifermi12(f)[source]#

Apply a rational function expansion to get the inverse fermi-dirac integral of order 1/2 when it is equal to f.

The maximum error is 4.19e-9_rt (see Antia ApJS 84,101 1993)

Parameters:

f (float) – Value of the Fermi integral

Return type:

float

pynucastro.neutrino_cooling.sneut5_mod.sneut5(rho, T, comp=None, *, abar=None, zbar=None, full_output=False)[source]#

Compute thermal neutrino losses (in erg/g/s) from the analytic fits of Itoh et al. ApJS 102, 411, 1996. Note that either a Composition object of abar/zbar need to be provided.

Parameters:

  • rho (float) – density

  • T (float) – temperature

  • comp (Composition) – composition of the plasma

  • abar (float) – mean molecular weight of the plasma (1 / sum{X_k/A_k})

  • zbar (float) – average charge of the plasma (abar * sum{X_k Z_k / A_k})

  • full_output (bool) – do we return the individual contributions to the cooling in addition to the total cooling rate?

Returns:

  • float – The cooling rate if full_output is False

  • tuple of (float, NeutrinoComponents) – The cooling rate and components if full_output is True

Contents