## Abstract

We present a novel method to solve the spatially homogeneous and isotropic relativistic Boltzmann equation. We employ a basis set of orthogonal polynomials dynamically adapted to allow for emergence of chemical non-equilibrium. Two time dependent parameters characterize the set of orthogonal polynomials, the effective temperature T(t) and phase space occupation factor ϒ{hooked}(t). In this first paper we address (effectively) massless fermions and derive dynamical equations for T(t) and ϒ{hooked}(t) such that the zeroth order term of the basis alone captures the particle number density and energy density of each particle distribution. We validate our method and illustrate the reduced computational cost and the ability to easily represent final state chemical non-equilibrium by studying a model problem that is motivated by the physics of the neutrino freeze-out processes in the early Universe, where the essential physical characteristics include reheating from another disappearing particle component (e^{±}-annihilation).

Original language | English (US) |
---|---|

Pages (from-to) | 896-916 |

Number of pages | 21 |

Journal | Journal of Computational Physics |

Volume | 281 |

DOIs | |

State | Published - Jan 5 2015 |

## Keywords

- Chemical non-equilibrium
- Orthogonal polynomial spectral method
- Relativistic Boltzmann equation

## ASJC Scopus subject areas

- Numerical Analysis
- Modeling and Simulation
- Physics and Astronomy (miscellaneous)
- General Physics and Astronomy
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics