TThermalNuBackground Class

module: MassiveNu file: massive_neutrinos

Description

Quantities for the neutrino background momentum distribution assuming thermal

Subroutines

init(this)
- TThermalNuBackground :: this
rho(this, am, rhonu)
- TThermalNuBackground :: this
- real(dl) intent(in) :: am
- real(dl) intent(out) :: rhonu
rho_P(this, am, rhonu, pnu)
- TThermalNuBackground :: this
- real(dl) intent(in) :: am
- real(dl) intent(out) :: rhonu
- real(dl) intent(out) :: pnu

Functions

real(dl)
drho(this, adotoa, am)

Compute the time derivative of the mean density in massive neutrinos
- TThermalNuBackground :: this
- real(dl) :: adotoa
- real(dl) intent(IN) :: am
real(dl)
find_nu_mass_for_rho(this, rho)

Get eigenstate mass given input density (rho is neutrino density in units of one massless) nu_mass=m_n*c**2/(k_B*T_nu0). Get number density n of neutrinos from rho_massless/n = int q^3/(1+e^q) / int q^2/(1+e^q)=7/180 pi^4/Zeta(3) then m = Omega_nu/N_nu rho_crit /n if non-relativistic
- TThermalNuBackground :: this
- real(dl) intent(in) :: rho

Properties

real(dl) :: dam

step in a*m
real(dl) :: r1 dimension(:), allocatable
real(dl) :: p1 dimension(:), allocatable
real(dl) :: dr1 dimension(:), allocatable
real(dl) :: dp1 dimension(:), allocatable
real(dl) :: ddr1 dimension(:), allocatable
real(dl) :: target_rho private