
Investigating the cluster production mechanism with isospin triggering:
Thermal models versus coalescence models
Apiwit Kittiratpattana1,3, Tom Reichert1,2, Pengcheng Li4,5, Ayut
Limphirat3, Christoph Herold3,∗, Jan Steinheimer6, Marcus Bleicher1,2
1Institut f¨ur Theoretische Physik, Goethe Universit¨at Frankfurt,
Max-von-Laue-Strasse 1, D-60438 Frankfurt am Main, Germany
2Helmholtz Research Academy Hesse for FAIR (HFHF),
GSI Helmholtz Center for Heavy Ion Physics, Campus Frankfurt,
Max-von-Laue-Str. 12, 60438 Frankfurt, Germany
3Center of Excellence in High Energy Physics & Astrophysics,
School of Physics, Suranaree University of Technology,
University Avenue 111, Nakhon Ratchasima 30000, Thailand
4School of Nuclear Science and Technology, Lanzhou University, Lanzhou 730000, China
5School of Science, Huzhou University, Huzhou 313000, China and
6Frankfurt Institute for Advanced Studies (FIAS),
Ruth-Moufang-Str.1, D-60438 Frankfurt am Main, Germany∗
Isospin triggering allows to distinguish coalescence from thermal production of light clusters in
heavy ion collisions. Triggering on Y(π−)−Y(π+) allows to select very neutron or proton rich
final states. The deuteron (cluster) production with coalescence (d∝n·p) leads then to an inverse
parabolic dependence of the deuteron yield on ∆Yπ. In contrast, in a thermal model, cluster
production is independent on ∆Yπ. The observation of a maximum deuteron (cluster) yield as
function of ∆Yπprovides confirmation of the coalescence mechanism.
I. INTRODUCTION
The exploration of the properties of matter governed
by the theory of strong interaction (Quantum-Chromo-
Dynamics, QCD) is a topic of highest interest. Ab-initio
calculations based on lattice QCD methods have shown
that such matter undergoes a transition at sufficiently
high temperatures and/or baryonic densities [1, 2]. Such
temperatures are e.g. reached in accelerator facilities like
the CERN-LHC, BNL-RHIC or CERN-SPS. At the high
density frontier, laboratory experiments are performed at
RHIC in the beam energy scan program, at GSI’s SIS18
accelerator or at the future FAIR facility. In nature the
high temperature transition from the deconfined Quark-
Gluon-Plasma state to a hadronic system happened ap-
proximately a few microseconds after the Big Bang, while
the high density regime is probed by neutron stars and
neutron star mergers. Especially neutron star merg-
ers have renewed the interest in the equation-of-state of
nuclear matter at highest densities [3] because gravita-
tional wave measurements might allow to pin down the
equation-of-state of QCD matter with very high precision
[4].
A central tool that is often used to infer the properties
of the created matter are light clusters, e.g. deuterons,
tritons and helium. For the production of such states
one uses generally two complementary approaches: The
statistical (thermal) model [5–13] or coalescence [14–31].
While both models provide similar results [25, 32] over a
wide range of collision energies they are very different in
their physics assumptions:
∗herold@g.sut.ac.th
(I) The thermal model assumes the creation of a fully
thermalized (mostly assumed grand canonical) fire-
ball, which means that the clusters are produced
at the chemical freeze-out at a temperature of
60-150 MeV (depending on the collision energies
probed in large systems like Au+Au or Pb+Pb)
from √sNN = 2.4−13000 GeV. An often discussed
problem with this model is the fact that lightly
bound clusters may not form or survive in such a
hot environment. This is also known from Big Bang
nucleosynthesis under the term of deuteron bottle-
neck, which means that deuterons and higher mass
light elements can only be formed if the tempera-
ture is on the order or below the binding energy of
a few MeV.
Even within the thermal model itself such a ten-
sion is visible in certain energy ranges [33], and the
clusters are often removed from the thermal fitting
as they worsen the quality of the thermal fit signif-
icantly [34].
(II) In contrast, the coalescence model assumes that
light clusters are formed at kinetic freeze-out, i.e.
after the last collisions/decays have ceased and the
system reaches the free-streaming regime. Here the
formation is possible due to lower temperatures and
due to the fact that no further collisions will destroy
the formed cluster. It is clear that at the earlier
chemical freeze-out, the temperature is higher and
the volume of the source is smaller than at the later
kinetic freeze-out where the temperature is lower
and the volume is larger [35].
Up to now it has not been possible to distinguish be-
tween both methods for cluster production, because the
arXiv:2210.11699v2 [nucl-th] 7 May 2023