Comparison of heavy-ion transport simulations: Collision integral in a box

TL;DR

This study compares 15 heavy-ion transport codes in a controlled box setup to evaluate their collision integral implementations, revealing phase-space fluctuation issues affecting Pauli blocking accuracy and guiding improvements in simulation strategies.

Contribution

First comprehensive comparison of collision integral handling in multiple transport codes using a controlled box setup, highlighting phase-space fluctuation impacts on Pauli blocking accuracy.

Findings

Collision rates match analytical results when Pauli blocking is suppressed.

Pauli blocking probabilities deviate due to phase-space fluctuations.

Most codes show a reduction in blocking probability, affecting system evolution.

Abstract

Simulations by transport codes are indispensable to extract valuable physics information from heavy ion collisions. In order to understand the origins of discrepancies between different widely used transport codes, we compare 15 such codes under controlled conditions of a system confined to a box with periodic boundary, initialized with Fermi-Dirac distributions at saturation density and temperatures of either 0 or 5 MeV. In such calculations, one is able to check separately the different ingredients of a transport code. In this second publication of the code evaluation project, we only consider the two-body collision term, i.e. we perform cascade calculations. When the Pauli blocking is artificially suppressed, the collision rates are found to be consistent for most codes (to within $1\%$ or better) with analytical results, or completely controlled results of a basic cascade code after eliminating the correlations within the same pair of colliding particles. In calculations with active Pauli blocking, the blocking probability was found to deviate from the expected reference values. The reason is found in substantial phase-space fluctuations and smearing tied to numerical algorithms and model assumptions in the representation of phase space. This results in the reduction of the blocking probability in most transport codes, so that the simulated system gradually evolves away from the Fermi-Dirac towards a Boltzmann distribution. As a result of this investigation, we are able to make judgements about the most effective strategies in transport simulations for determining the collision probabilities and the Pauli blocking. Investigation in a similar vein of other ingredients in transport calculations, like the mean field propagation or the production of nucleon resonances and mesons, will be discussed in the future publications.