Neutrino Quantum Kinetics in Compact Objects

This contains the codebase and data from "Neutrino Quantum Kinetics in
Compact Objects" by Sherwood Richers, Gail McLaughlin, Alexey
Vlasenko, and James P. Kneller. (2019, submitted to PrD)

======================
All f.h5 file contents
======================
r(cm)  - proxy for the time (r=ct)
dr_block(cm) - proxy for the block timestep (dr=cdt). Stored for
               recovery purposes
dr_int(cm)   - proxy for the interaction integration timestep
               (dr=cdt). Stored for recovery purposes.
dr_osc(cm)   - proxy for the oscillation integration timestep
               (dr=cdt). Stored for recovery purposes.
fmatrixf     - dimensionless distribution function. Accessed with
               [time][helicity][energy][flavor1][flavor2][0=real,1=imaginary]

==================
Directory contents
==================

fiducial - simulation output from the fiducial case of a background
with rho=1e10g/ccm, T=10MeV, Ye=0.3, the HShen EOS, including
oscillations, neutrino absorption on nucleons and nuclei, inelastic
electron scattering, elastic nucleon scattering, electron-positron
pair production, effective-absorption nucleon-nucleon bremsstrahlung,
and neutrino-neutrino scattering and pair processes. Using 25 energy
bins covering up to 100 MeV.

noosc_full - simulation output from the fiducial case with the same
circumstanses as above, except without oscillations and using 50
energy bins up to 100 MeV.

noosc_supernova_background.dat - snapshot from a 1D core-collapse
                         supernova simulation. The only relevant
                         quantities are the radius (col. 2), density
                         (col. 3), electron fraction (col. 4),
                         temperature (col. 5), electron chemical
                         potential (col. 9), proton chemical potential
                         (col. 10), and neutron chemical potential
                         (col. 11)
            
noosc_*_tdecohere.dat - decoherence time (col. 4) for neutrinos at
                every radial point (col. 1), each helicity
            (col. 2), and every energy (col. 3), using 50
            energy bins at integer multiples of 2 MeV.

   full - calculations done including the full set of
          collision interactions listed in fiducial.

   nonu4 - calculations done including the full set of
           collision interactions listed in fiducial
           except neutrino-neutrino scattering/pair
           processes.