Boltzmann equation solver adapted to emergent chemical non-equilibrium

Jeremiah Birrell, Jon Wilkening, Johann Rafelski

Research output: Contribution to journalArticle

3 Citations (Scopus)

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 languageEnglish (US)
Pages (from-to)896-916
Number of pages21
JournalJournal of Computational Physics
Volume281
DOIs
StatePublished - Jan 5 2015

Fingerprint

Boltzmann equation
Polynomials
Fermions
polynomials
Physics
occupation
neutrinos
flux density
universe
fermions
costs
Costs
physics
heating
Temperature
temperature

Keywords

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

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy (miscellaneous)

Cite this

Boltzmann equation solver adapted to emergent chemical non-equilibrium. / Birrell, Jeremiah; Wilkening, Jon; Rafelski, Johann.

In: Journal of Computational Physics, Vol. 281, 05.01.2015, p. 896-916.

Research output: Contribution to journalArticle

@article{c44b64d663b34d14a40c2cea2693bd9d,
title = "Boltzmann equation solver adapted to emergent chemical non-equilibrium",
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).",
keywords = "Chemical non-equilibrium, Orthogonal polynomial spectral method, Relativistic Boltzmann equation",
author = "Jeremiah Birrell and Jon Wilkening and Johann Rafelski",
year = "2015",
month = "1",
day = "5",
doi = "10.1016/j.jcp.2014.10.056",
language = "English (US)",
volume = "281",
pages = "896--916",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",

}

TY - JOUR

T1 - Boltzmann equation solver adapted to emergent chemical non-equilibrium

AU - Birrell, Jeremiah

AU - Wilkening, Jon

AU - Rafelski, Johann

PY - 2015/1/5

Y1 - 2015/1/5

N2 - 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).

AB - 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).

KW - Chemical non-equilibrium

KW - Orthogonal polynomial spectral method

KW - Relativistic Boltzmann equation

UR - http://www.scopus.com/inward/record.url?scp=84916893280&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84916893280&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2014.10.056

DO - 10.1016/j.jcp.2014.10.056

M3 - Article

AN - SCOPUS:84916893280

VL - 281

SP - 896

EP - 916

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -