PHQMD

Parton-Hadron-Quantum-Molecular Dynamics

A microscopic N-body transport approach for heavy-ion collisions, cluster formation, and hypernuclei production.

About PHQMD

The Parton-Hadron-Quantum-Molecular Dynamics (PHQMD) transport approach is an n-body microscopic framework designed to describe nuclear cluster and hypernucleus formation, as well as general particle production, in relativistic heavy-ion reactions.

PHQMD combines baryon propagation from the Quantum Molecular Dynamics (QMD) model with the dynamical properties and interactions of hadronic and partonic degrees of freedom from the Parton-Hadron-String Dynamics (PHSD) approach.

A recent development in PHQMD is the inclusion of a momentum-dependent nucleon potential in addition to the static, density-dependent Skyrme interaction. This extension enables three distinct equation-of-state (EoS) scenarios: two static cases (“soft” and “hard”), which differ in compressibility, and a soft momentum-dependent EoS calibrated to pA elastic scattering data.

Cluster Production in PHQMD

In contrast to coalescence or statistical approaches, PHQMD forms clusters dynamically through baryon interactions described within QMD. This allows the propagation of the n-body Wigner density and n-body phase-space correlations, which are essential for cluster formation.

Clusters can be identified in PHQMD using three different algorithms:

  1. Potential mechanism: The attractive potential between baryons with small relative momentum can keep them close together and lead to bound groups of nucleons. These co-moving nucleons can be identified during the dynamical evolution using the advanced Minimum Spanning Tree (aMST) method or the SACA (Simulated Annealing Cluster Algorithm), which searches for the most bound configuration of nucleons and hyperons.

    In the MST algorithm, two nucleons i and j are considered bound if |ri* - rj*| < rclus, where the positions are evaluated in the center-of-mass frame of the pair and rclus = 4 fm, roughly corresponding to the range of the attractive NN potential. In addition, aMST requires a negative binding energy, EB < 0, and includes a stabilization procedure to prevent the spontaneous decay of clusters caused by the semiclassical treatment. It should be noted that (a)MST acts as a cluster recognition tool, applied perturbatively, rather than as the actual formation mechanism, since QMD propagates baryons and not pre-formed clusters.
  2. Kinetic mechanism: Deuterons are produced through catalytic hadronic reactions πNN ↔ πd and NNN ↔ Nd in different isospin channels. Quantum effects are incorporated through excluded-volume corrections and projection onto the deuteron wave function in momentum space, which reduces the production probability, especially at target and projectile rapidities.
  3. Coalescence mechanism: A proton and neutron form a deuteron if their phase-space distance at freeze-out satisfies |r1 - r2| ≤ 3.575 fm and |p1 - p2| ≤ 285 MeV/c. In PHQMD, this perturbative method is used only for model studies and comparison purposes.

How To Cite PHQMD

Recommended primary references for PHQMD mode and usage.

Parton-Hadron-Quantum-Molecular Dynamics (PHQMD) - A Novel Microscopic N-Body Transport Approach for Heavy-Ion Collisions, Dynamical Cluster Formation and Hypernuclei Production

J. Aichelin, E. Bratkovskaya, A. Le Fevre, V. Kireyeu, V. Kolesnikov, Y. Leifels, V. Voronyuk, G. Coci · Phys. Rev. C 101 (2020) 044905

Show abstract

Cluster and hypernuclei production in heavy-ion collisions is presently under active experimental and theoretical investigation. Since clusters are weekly bound objects, their production is very sensitive to the dynamical evolution of the system and its interactions. The theoretical description of cluster formation is related to the n-body problem. Here we present the novel n-body dynamical transport approach PHQMD (Parton-Hadron-Quantum-Molecular Dynamics) which is designed to provide a microscopic description of nuclear cluster and hypernucleus formation as well as of general particle production in heavy-ion reactions at relativistic energies. In difference to the coalescence or statistical models, often used for the cluster formation, in PHQMD clusters are formed dynamically due to the interactions between baryons described on a basis of Quantum Molecular Dynamics (QMD)which allows to propagate the n-body Wigner density and n-body correlations in phase-space, essential for the cluster formation. The clusters are identified by the MST (Minimum Spanning Tree) or the SACA ('Simulated Annealing Cluster Algorithm') algorithm which finds the most bound configuration of nucleons and clusters. Collisions among hadrons as well as Quark-Gluon-Plasma formation and parton dynamics in PHQMD are treated in the same way as in the established PHSD (Parton-Hadron-String Dynamics)transport approach. In order to verify our approach with respect to the general dynamics we present here the first PHQMD results for general 'bulk' observables such as rapidity distributions and transverse mass spectra for hadrons (\pi, K, \bar K, p, \bar p, \Lambda, \bar \Lambda) from SIS to RHIC energies. We find a good description of the 'bulk' dynamics which allows us to proceed with the results on cluster production, including hypernuclei.

PHQMD result q1
’Rise and fall’ of the multiplicity of clusters with Z ∈ [3, 30] as a function of the total bound charge

Dynamical mechanisms for deuteron production at mid-rapidity in relativistic heavy-ion collisions from energies available at the GSI Schwerionensynchrotron to those at the BNL Relativistic Heavy Ion Collider

G. Coci, S. Gläßel, V. Kireyeu, J. Aichelin, C. Blume, E. Bratkovskaya, V. Kolesnikov, V. Voronyuk · Phys. Rev. C 108 (2023) 014902

Show abstract

The understanding of the mechanisms for the production of weakly bound clusters, such as a deuteron d, in heavy-ion reactions at mid-rapidity is presently one of the challenging problems which is also known as the "ice in a fire" puzzle. In this study we investigate the dynamical formation of deuterons within the Parton-Hadron-Quantum-Molecular Dynamics (PHQMD) microscopic transport approach and advance two microscopic production mechanisms to describe deuterons in heavy-ion collisions from SIS to RHIC energies: kinetic production by hadronic reactions and potential production by the attractive potential between nucleons. Differently to other studies, for the "kinetic" deuterons we employ the full isospin decomposition of the various \pi NN\leftrightarrow \pi d, NNN\leftrightarrow N d channels and take into account the finite size properties of the deuteron by means of an excluded volume condition in coordinate space and by the projection onto the deuteron wave function in momentum space. We find that considering the quantum nature of the deuteron in coordinate and momentum space reduces substantially the kinetic deuteron production in a dense medium as encountered in heavy-ion collisions. If we add the "potential" deuterons by applying an advanced Minimum Spanning Tree (aMST) procedure, we obtain good agreement with the available experimental data from SIS energies up to the top RHIC energy.

PHQMD result q2
The mid-rapidity |y| < 0.3 excitation function for dN/dy of deuterons for Au+Au central collisions

Basic References To The PHQMD Model

Selected PHQMD references on cluster formation and dynamical mechanisms.

Cluster and hyper-cluster production in relativistic heavy-ion collisions within the PHQMD approach

S. Gläßel, V. Kireyeu, V. Voronyuk, J. Aichelin, C. Blume, E. Bratkovskaya, G. Coci, V. Kolesnikov, M. Winn · Phys. Rev. C 105 (2022) 014908

Show abstract

We study cluster and hypernuclei production in heavy-ion collisions at relativistic energies employing the Parton-Hadron-Quantum-Molecular-Dynamics (PHQMD) approach, a microscopic n-body transport model based on the QMD propagation of the baryonic degrees of freedom with density dependent 2-body potential interactions. All other ingredients of PHQMD, including the collision integral and the treatment of the quark-gluon plasma (QGP) phase, are adopted from the Parton-Hadron-String Dynamics (PHSD) approach. In PHQMD the cluster formation occurs dynamically, caused by the interactions. The clusters are recognized by the Minimum Spanning Tree (MST) algorithm. We present the PHQMD results for cluster and hypernuclei formation in comparison with the available experimental data at AGS, SPS, RHIC-BES and RHIC fixed target energies. We also provide predictions on cluster production for the upcoming FAIR and NICA experiments. PHQMD allows to study the time evolution of formed clusters and the origin of their production, which helps to understand how such weakly bound objects are formed and survive in the rather dense and hot environment created in heavy-ion collisions. It offers therefore an explanation of the 'ice in the fire' puzzle.

PHQMD result q3
Transverse momentum distribution of 4ΛH for different rapidity intervals for central Au+Au collisions

Deuteron production in ultrarelativistic heavy-ion collisions: A comparison of the coalescence and the minimum spanning tree procedure

V. Kireyeu, J. Steinheimer, J. Aichelin, M. Bleicher, E. Bratkovskaya · Phys. Rev. C 105 (2022) 044909

Show abstract

The formation of deuterons in heavy-ion collisions at relativistic energies is investigated by employing two recently advanced models -- the Minimum Spanning Tree (MST) method and the coalescence model by embedding them in the PHQMD and the UrQMD transport approaches. While the coalescence mechanism combines nucleons into deuterons at the kinetic freeze-out hypersurface, the MST identifies the clusters during the different stages of time evolution. We find that both clustering procedures give very similar results for the deuteron observables in the UrQMD as well as in the PHQMD environment. Moreover, the results agree well with the experimental data on deuteron production in Pb+Pb collisions at \sqrt{s_{NN}} = 8.8 GeV (selected for the comparison of the methods and models in this study). A detailed investigation shows that the coordinate space distribution of the produced deuterons differs from that of the free nucleons and other hadrons. Thus, deuterons are not destroyed by additional rescattering.

PHQMD result q4
Rapidity distribution of protons at the "freeze-out" time before the coalescence

Cluster formation near midrapidity - can the mechanism be identified experimentally?

V. Kireyeu, G. Coci, S. Gläßel, J. Aichelin, C. Blume, E. Bratkovskaya · Phys. Rev. C 109 (2024) 044906

Show abstract

The formation of weakly bound clusters in the hot and dense environment at midrapidity is one of the surprising phenomena observed experimentally in heavy-ion collisions from a low center of mass energy of a few GeV up to a ultra-relativistic energy of several TeV. Three approaches have been advanced to describe the cluster formation: coalescence at kinetic freeze-out, cluster formation during the entire heavy-ion collision by potential interaction between nucleons and deuteron production by hadronic reactions. We identify experimental observables, which can discriminate these production mechanisms for deuterons.

PHQMD result q5
The transverse momentum distributions of deuterons for central Au+Au collisions

Cluster dynamics studied with the phase-space minimum spanning tree approach

V. Kireyeu · Phys. Rev. C 103 (2021) 054905

Show abstract

The origin of weakly bound objects like clusters and hypernuclei, observed in heavy-ion collisions, is of theoretical and experimental interest. It is in the focus of the experiments at RHIC and LHC since it is not evident how such weakly bound objects can survive in an environment whose hadronic decay products point to a temperature of the order of 150 MeV. It is as well one of the key research topics in the future facilities of FAIR and NICA which are under construction in Darmstadt (Germany) and Dubna (Russia), respectively. We present here first results on the cluster dynamics within the model-independent cluster recognition library "phase-space Minimum Spanning Tree" (psMST) applied to different transport approaches: PHQMD, PHSD, SMASH and UrQMD. The psMST is based on the "Minimum Spanning Tree" (MST) algorithm for the cluster recognition which exploits correlations in coordinate space, and it is extended to correlations of baryons in the clusters in momentum space. We show the sensitivity of the cluster formation on the microscopic realization of the n-body dynamics and on the potential interactions in heavy-ion collisions.

PHQMD result psMST
The rapidity distributions of clusters with the mass number 4 ≤ A ≤ 20

Constraints on the equation-of-state from low energy heavy-ion collisions within the PHQMD microscopic approach with momentum-dependent potential

V. Kireyeu, V. Voronyuk, M. Winn, S. Gläßel, J. Aichelin · e-Print: 2411.04969 [nucl-th], 2024

Show abstract

We investigate the influence of the equation-of-state (EoS) of nuclear matter on collective observables, the directed (v_1) and the elliptic flow (v_2) of nucleons and light clusters in heavy-ion collisions at GeV beam energies employing the Parton-Hadron-Quantum-Molecular Dynamics (PHQMD) microscopic transport approach. Here the clusters are formed dynamically during the entire heavy-ion collision by potential interaction between nucleons, including additionally deuteron production by hadronic kinetic reactions. We employ three different EoS - realized via potential interactions: two static EoS, dubbed 'soft' and 'hard', which differ in the compressibility modulus, as well as a soft momentum dependent EoS, adjusted to pA elastic scattering data. We find that the momentum dependent potential has different consequences for rapidity and transverse momentum spectra than for flow coefficients. We obtain the best description of the HADES and FOPI data on the directed and elliptic flow coefficients of protons and light clusters applying a momentum dependent EoS. Moreover, we observe a scaling behavior of v_2 versus p_T with atomic number A. Finally we demonstrate that flow observables can help to identify the cluster production mechanisms.

PHQMD result psMST
v2 of deuterons as a function of rapidity for 20-30% central Au+Au collisions

Probing the nuclear equation of state with clusters and hypernuclei

Y. Zhou, S. Glässel, L. Yue-Hang, V. Kireyeu, J. Zhao, H. Liu, C. Blume, I. Vassiliev, V. Voronyuk, M. Winn, N. Herrmann, Y. Wang, N. Xu, J. Aichelin, E. Bratkovskaya · Phys.Rev.C 113 (2026) 014909

Show abstract

The study of the nuclear equation-of-state (EoS) is a one of the primary goals of experimental and theoretical heavy-ion physics. The comparison of recent high statistics data from the STAR Collaboration with transport models provides a unique possibility to address this topic in a yet unexplored energy domain. Employing the microscopic N-body Parton-Hadron-Quantum-Molecular Dynamics (PHQMD) transport approach, which allows to describe the propagation and interactions of hadronic and partonic degrees of freedom including cluster and hyper-nucleus formation and dynamics, we investigate the influence of different EoS on bulk observables, the multiplicity, p_T and rapidity distributions of protons, \Lambdas and clusters up to A=4 as well as their influence on the collective flow. We explore three different EoS: two static EoS, dubbed 'soft' and 'hard', which differ in the compressibility modulus, as well as a soft momentum dependent EoS. We find that a soft momentum dependent EoS reproduces most baryon and cluster observables, including the flow observables, quantitatively, however, hard EOS show a similar trend.

PHQMD result psMST
v1 of deuterons as a function of pT for central Au+Au collisions

Results

Selected PHQMD comparison plots and links to references.

q1: rise and fall of cluster multiplicity
q1: Multiplicity of clusters versus Zbound2.

"Rise and fall" of the multiplicity of clusters with Z in [3, 30] as a function of total bound charge Zbound2.

q2: average charge of the largest cluster
q2: Average charge of the largest cluster versus bound charge.

PHQMD results with SACA at different times compared to ALADIN data.

q3: rapidity distribution from FOPI comparison
q3: Scaled rapidity distributions in central Au+Au collisions.

Comparison of PHQMD cluster-identification strategies (MST/SACA) to FOPI observables.

Compilation of experimental cluster data used for PHQMD validation

Experiment System Energy Species Paper
FOPI Au+Au E_lab = 1.5 A.GeV d, t 10.1016/j.nuclphysa.2010.09.008
E802 Au+Pb E_lab = 10.6 A.GeV d 10.1016/j.nuclphysa.2010.09.008
E864 Au+Pb / Au+Pt E_lab = 10.6 A.GeV d,t,3He, 4He PhysRevC.61.064908
PhysRevLett.83.5431
PhysRevC.70.024902
E878 + E886 Au+Au E_lab = 10.6 A.GeV d, t, 3He, 4He 10.1103/PhysRevC.58.1155
NA49 Pb+Pb P_lab = 20, 30, 40, 80, 158 GeV/c d,t,3He 10.1103/PhysRevC.94.044906
STAR BES Au+Au sqrt(sNN) = 7.7, 11.5, 19.6, 27, 39, 62.4, 200 GeV d 10.1103/PhysRevC.99.064905

Talks

Selected PHQMD talks and slide links.

A Microscopic Transport Framework for Heavy-Ion Collisions from High to Low Energies - PHSD/PHQMD

Jiaxing Zhao · Non-equilibrium QCD and Transport @ Huizhou · Dec-2025

Cluster production in PHQMD

Elena Bratkovskaya · EMMI Workshop, GSI-Darmstadt · Nov-2025

PHSD/PHQMD for FAIR

Elena Bratkovskaya · 45th CBM Collaboration Meeting, GSI-Darmstadt · Feb-2025

Kinetic and potential mechanisms for deuteron production in HICs

Gabriele Coci · STRONG, Italy · Oct-2023

Dynamical cluster and hypernuclei production in heavy ion collisions

Elena Bratkovskaya · SQM 2021

Prospects of studying the production of hypernuclei in heavy-ion interactions at the NICA collider at JINR

V. Kireyeu et al. · PANIC2021

Prospects for the study of the strangeness production within PHQMD model

V. Kireyeu · The 5th International Conference on Particle Physics and Astrophysics, MEPhI, Russia · Oct-2020

PHQMD: a microscopic transport approach to study cluster production

J. Aichelin · ICNFP 2020

Overview over possibilities to determine fragments dynamically from transport codes

Joerg Aichelin · EMMI Workshop 2019

Clusters and hypernuclei production within PHQMD

Elena Bratkovskaya · EMMI Workshop 2019

News from PHQMD

Joerg Aichelin · Dubna workshop 2019

PHQMD: A Novel Microscopic N-Body Transport Approach

Elena Bratkovskaya · SQM 2019

The PHQMD model for the formation of nuclear clusters and hypernuclei in heavy-ion collisions

Viktar Kireyeu · Nucleus-2019

Prospects for the study of hyperons and hypernuclei at NICA collider

Vadim Kolesnikov · Nucleus-2019