pyccoonStandard Model BBN test
This test checks that in the universe filled with photons, electrons and neutrinos:
-
is not conserved by a factor around
1.401and precise details of this process - neutrino non-equilibrium corrections reproduce the results of the Dolgov-Hansen-Semikoz papers
import os
from plotting import RadiationParticleMonitor
from particles import Particle
from library.SM import particles as SMP, interactions as SMI
from evolution import Universe
from common import UNITS, Params
folder = os.path.split(__file__)[0]
T_final = 0.0008 * UNITS.MeV
params = Params(T=10. * UNITS.MeV,
dy=0.025)
universe = Universe(params=params, folder=folder)
photon = Particle(**SMP.photon)
electron = Particle(**SMP.leptons.electron)
neutrino_e = Particle(**SMP.leptons.neutrino_e)
neutrino_mu = Particle(**SMP.leptons.neutrino_mu)
neutrino_tau = Particle(**SMP.leptons.neutrino_tau)
neutron = Particle(**SMP.hadrons.neutron)
proton = Particle(**SMP.hadrons.proton)
universe.add_particles([
photon,
electron,
neutrino_e,
neutrino_mu,
neutrino_tau,
neutron,
proton
])
universe.init_oscillations(SMP.leptons.oscillations_map(), (neutrino_e, neutrino_mu, neutrino_tau))
universe.interactions += ([SMI.baryons_interaction(neutron=neutron, proton=proton, electron=electron, neutrino=neutrino_e)] +
SMI.neutrino_interactions(leptons=[electron],
neutrinos=[neutrino_e, neutrino_mu, neutrino_tau]))
universe.init_kawano(electron=electron, neutrino=neutrino_e)
if universe.graphics:
universe.graphics.monitor([
(neutrino_e, RadiationParticleMonitor),
(neutrino_mu, RadiationParticleMonitor),
(neutrino_tau, RadiationParticleMonitor)
])
universe.evolve(T_final)
