This will be a toy Boltzmann Quantum Kinetic solver for a very small number of neutrinos in a two-beam scenario. The purpose of this project is NOT to solve all of neutrino quantum kinetics!

Instead, it is to use neutrino quantum kinetics as a vehicle to explore numerical integration techniques and explore the difficulties that even overly simplified physical problems can pose.

Please reach out if you would like to collaborate, play yourself, or expand on what the code does.

Compile:

make clean make

Adjust runtime parameters: parameters.h

Execute: ./boltzmann_solver

evaluate results: python3 plot_results.py --save-gif