Track papers across labs, topics, and trends.
Search papers, labs, and topics across Lattice.
The authors derive the relativistic Anderson-Witting relaxation-time Boltzmann equation as the mean-field limit of a stochastic $N$-body system with exact microcanonical momentum resets, implementing it within a 3D continuous-momentum Monte Carlo solver called ShARK. By formulating the relaxation kernel in continuous phase space rather than on a discrete momentum grid, the framework circumvents the severe truncation errors that degrade standard lattice Boltzmann methods far from local thermodynamic equilibrium. Numerical validations across Bjorken and Gubser conformal expansions show that finite-$N$ phase-space sampling converges robustly to analytical solutions and the Jüttner-Boltzmann equilibrium.
Discretization errors in far-from-equilibrium relativistic transport vanish when relaxation kinetics are reformulated as exact, continuous-momentum $N$-body microcanonical phase-space redraws.
We show that the Anderson-Witting relaxation-time approximation (relativistic Bhatnagar-Gross-Krook (BGK) equation) emerges as the $N\rightarrow \infty$ limit of a stochastic $N$-body particle gas, in which collisions perform a full microcanonical momentum redraw that enforces local energy-momentum conservation at every event. Using the RAMBO algorithm to sample the Lorentz-invariant phase space, we derive the finite-$N$ single-particle momentum spectrum and show that it converges to the Jüttner-Boltzmann equilibrium distribution as the local particle number grows. We couple this relaxation kernel to an advection-relaxation splitting scheme to construct ShARK (Stochastic Advection Relaxation Kinetics), a 3D Monte Carlo relativistic solver for the relaxation time approximation Boltzmann equation. The framework is conceptually related to lattice Boltzmann methods, but formulated in continuous momentum space to avoid the accuracy loss caused by finite-order momentum discretizations far from local equilibrium. We validate the numerical results against analytical solutions for two highly symmetric conformal expansions, the Bjorken and Gubser flows. The framework provides a first-principles route to far-from-equilibrium relativistic transport, with applications ranging from heavy-ion collisions to expanding astrophysical plasmas.