PHSD

Parton-Hadron-String Dynamics

A microscopic off-shell transport approach for the description of strongly interacting hadronic and partonic matter in and out of equilibrium.

About PHSD

The Parton-Hadron-String Dynamics (PHSD) is a microscopic off-shell transport approach for the description of strongly interacting hadronic and partonic matter in and out-of-equilibrium. It is based on the solution of Kadanoff-Baym equations in first-order gradient expansion in phase space, allowing a causal description of the time evolution of nonperturbative interacting systems. PHSD consistently describes the full evolution of a relativistic heavy-ion collision from the initial hard scatterings and string formation through the dynamical deconfinement phase transition to the quark-gluon plasma, followed by hadronization and subsequent hadronic interactions. The PHSD has been developed by the Giessen/Frankfurt groups on the basis of the Hadron-String Dynamics transport approach (HSD).

Au+Au at 200 GeV
Au+Au at 200 GeV
Click to download the full-size video

In PHSD the transition from partonic (quarks and gluons) to hadronic degrees of freedom is described by covariant transition rates for quark-antiquark fusion to mesonic resonances or three quarks (antiquarks) to baryonic states, i.e. by a dynamical hadronization obeying flavor current-conservation, color neutrality and energy-momentum conservation. The two-particle correlations resulting from finite parton spectral widths are treated dynamically by generalized off-shell Cassing-Juchem transport equations that go beyond mean-field or Boltzmann approximations.

The transport-theoretical description of quarks and gluons in PHSD is based on the Dynamical Quasi-Particle Model (DQPM), constructed to reproduce lattice-QCD results for the quark-gluon plasma in thermodynamic equilibrium. DQPM provides mean fields and effective two-body interactions for gluons and quarks. Close to the phase transition, dynamically massive quarks and antiquarks form resonant pre-hadronic color-dipole states (q-qbar or qqq) with high invariant mass that sequentially decay to meson and baryon octets. The resulting hadronization increases total entropy and thus remains consistent with the second law of thermodynamics.

Click to download the full-size video

The PHSD approach has been applied to nucleus-nucleus collisions from low SPS to LHC energies to explore space-time regions of partonic matter. It provides a consistent description of bulk heavy-ion observables - rapidity spectra, transverse-mass distributions, and azimuthal asymmetries (v1, v2, v3, v4) - and has been successfully used for dilepton analyses from hadronic and partonic sources.

Pb+Pb at 158 AGeV
Pb+Pb at 158 AGeV

Equilibrium properties of QGP were also studied with PHSD in a finite box with periodic boundary conditions at fixed temperature T. The ratio of shear viscosity to entropy density shows a minimum close to Tc = 160 MeV and approaches the perturbative QCD limit at higher temperatures, in line with lattice-QCD results. This indicates a strongly interacting liquid-like QGP rather than a weakly interacting gas of partons.

How To Cite PHSD

Recommended primary references for PHSD method and implementation.

Parton transport and hadronization from the dynamical quasiparticle point of view

W. Cassing and E. L. Bratkovskaya · Phys. Rev. C 78 (2008) 034919

Show abstract

The hadronization of an expanding partonic fireball is studied within the Parton-Hadron-Strings Dynamics (PHSD) approach which is based on a dynamical quasiparticle model (DQPM) matched to reproduce lattice QCD results in thermodynamic equilibrium. Apart from strong parton interactions the expansion and development of collective flow is driven by strong gradients in the parton mean-fields. An analysis of the elliptic flow v2 demonstrates a linear correlation with the spatial eccentricity epsilon as in case of ideal hydrodynamics. The hadronization occurs by quark-antiquark fusion or 3 quark/3 antiquark recombination which is described by covariant transition rates. Since the dynamical quarks become very massive, the formed resonant pre-hadronic color-dipole states (q-qbar or qqq) are of high invariant mass, too, and sequentially decay to the groundstate meson and baryon octets increasing the total entropy. This solves the entropy problem in hadronization in a natural way. The resulting particle ratios turn out to be in line with those from a grandcanonical partition function at temperature T approximately 170 MeV, rather independent from the initial temperature, and indicate an approximate strangeness equilibration.

Time evolution of the elliptic flow v2.
Time evolution in the parton, meson and baryon number.

Parton-Hadron-String Dynamics: an off-shell transport approach for relativistic energies

W. Cassingand E. L. Bratkovskaya · Nucl. Phys. A 831 (2009) 215-242

Show abstract

The dynamics of partons, hadrons and strings in relativistic nucleus-nucleus collisions is analyzed within the Parton-Hadron-String Dynamics (PHSD) transport approach, which is based on a Dynamical QuasiParticle Model (DQPM) for partons matched to reproduce lattice-QCD results in thermodynamic equilibrium. Scalar- and vector-interaction densities are extracted from the DQPM as well as effective scalar- and vector-mean fields for the partons. The transition from partonic to hadronic degrees of freedom is described by covariant transition rates for the fusion of quark-antiquark pairs or three quarks (antiquarks), obeying flavor current-conservation, color neutrality and energy-momentum conservation. Since the dynamical quarks and antiquarks become very massive close to the phase transition, the formed resonant pre-hadronic color-dipole states (q-qbar or qqq) are of high invariant mass, too, and sequentially decay to the groundstate meson and baryon octets increasing the total entropy. The PHSD approach is applied to nucleus-nucleus collisions from 20 to 160 A GeV to explore partonic matter and shows a sizeable influence on transverse-mass distributions of kaons and on multi-strange antibaryon production.

partonic energy fraction with time
partonic energy fraction with time

Exploring the partonic phase at finite chemical potential within an extended off-shell transport approach

P. Moreau, O. Soloveva, L. Oliva, T. Song, W. Cassing, E. Bratkovskaya· Phys. Rev. C, 100 (2019) 1, 014911

Show abstract

We extend the Parton-Hadron-String Dynamics (PHSD) transport approach in the partonic sector by explicitly calculating the total and differential partonic scattering cross sections as a function of temperature T and baryon chemical potential \mu_B on the basis of the effective propagators and couplings from the Dynamical QuasiParticle Model (DQPM) that is matched to reproduce the equation of state of the partonic system above the deconfinement temperature T_c from lattice QCD. The ratio of shear viscosity \eta over entropy density s, i.e. \eta/s, is evaluated using the collisional widths and compared to lQCD calculations for \mu_B = 0 as well. We find only a very modest change of \eta/s with the baryon chemical \mu_B. This also holds for a variety of hadronic observables from central A+A and C+Au collisions in the energy range 5 GeV \leq \sqrt{s_{NN}} \leq 200 GeV when implementing the differential cross sections into the PHSD approach. We only observe small differences in the strangeness and antibaryon sector with practically no sensitivity of rapidity and p_T distributions to the \mu_B dependence of the partonic cross sections. Since we find only small traces of a \mu_B-dependence in heavy-ion observables - although the effective partonic masses and widths as well as their partonic cross sections clearly depend on \mu_B - this implies that one needs a sizable partonic density and large space-time QGP volume to explore the dynamics in the partonic phase. These conditions are only fulfilled at high bombarding energies where \mu_B is, however, rather low. On the other hand, when decreasing the bombarding energy and thus increasing \mu_B, the hadronic phase becomes dominant and accordingly, it will be difficult to extract signals from the partonic dynamics based on "bulk" observables.

The energy density for central Pb+Pb collisions
The energy density for central Pb+Pb collisions

Basic References To The PHSD Model

From Kadanoff-Baym dynamics to off-shell parton transport

W. Cassing · Eur. Phys. J. ST 168 (2009) 3-87

Show abstract

This review provides a written version of lectures on non-equilibrium quantum field theory and shows the way from quantum-field theory in two-particle irreducible approximation to the Kadanoff-Baym equations and practical approximation schemes in phase space. This leads to an off-shell transport theory that incorporates the underlying quantum physics and can be solved in a testparticle representation for weak and strong coupling regimes. Applications to dilepton production in heavy-ion reactions and to hadronization from partonic to hadronic degrees of freedom are discussed.

dilepton emission rate DQPM vs Born
dilepton emission rate DQPM vs Born

Parton-Hadron-String Dynamics at relativistic collider energies

E. L. Bratkovskaya, W. Cassing, V. P. Konchakovski, O. Linnyk · Nucl. Phys. A 856 (2011) 162-182

Show abstract

The Parton-Hadron-String Dynamics (PHSD) transport approach is applied to nucleus-nucleus collisions at RHIC energies with respect to differential hadronic spectra in comparison to available data. PHSD is based on DQPM partons matched to lattice-QCD results and includes covariant transition rates for hadronization with all conservation laws. Dynamical studies for Au+Au collisions at top RHIC energy are compared to HSD results and to STAR, PHENIX, BRAHMS and PHOBOS data. The approach gives a reasonable description of rapidity distributions, transverse mass spectra and elliptic flow, and shows approximate quark-number scaling of v2 for hadrons.

The hadron elliptic flow v2
Scaled elliptic flow v2/nq versus transverse kinetic energy per constituent quark.

Azimuthal anisotropies for Au+Au collisions in the parton-hadron transient energy range

V. P. Konchakovski, E. L. Bratkovskaya, W. Cassing, V. D. Toneev, V. Voronyuk · Phys. Rev. C 85 (2012) 044922

Show abstract

The azimuthal anisotropies of the collective transverse flow of charged hadrons are investigated in a wide range of heavy-ion collision energies within the microscopic parton-hadron-string dynamics (PHSD) transport approach which incorporates explicit partonic degrees of freedom in terms of strongly interacting quasiparticles (quarks and gluons) in line with an equation of state from lattice QCD as well as the dynamical hadronization and hadronic collision dynamics in the final reaction phase. The experimentally observed increase of the elliptic flow v_2 of charged hadrons with collision energy is successfully described in terms of the PHSD approach. The PHSD scaling properties of various collective observables are confronted with experimental data as well as with hydrodynamic predictions. The analysis of higher-order harmonics v_3 and v_4 in the azimuthal angular distribution shows a similar tendency of growing deviations between partonic and purely hadronic models with increasing collision energy. This demonstrates that the excitation functions of azimuthal anisotropies reflect the increasing role of quark-gluon degrees of freedom in the early phase of relativistic heavy-ion collisions. Furthermore, the specific variation of the ratio v_4/(v_2)^2 with respect to bombarding energy, centrality and transverse momentum is found to provide valuable information on the underlying dynamics.

The evolution of the parton fraction of the total energy density.
Impact parameter dependence of anisotropic flows of charged particles.

Tomography of the Quark-Gluon-Plasma by Charm Quarks

T. Song, H. Berrehrah, D. Cabrera, J. M. Torres-Rincon, L. Tolos · Phys.Rev.C 92 (2015) 014910

Show abstract

We study charm production in ultra-relativistic heavy-ion collisions by using the Parton-Hadron-String Dynamics (PHSD) transport approach. The initial charm quarks are produced by the Pythia event generator tuned to fit the transverse momentum spectrum and rapidity distribution of charm quarks from Fixed-Order Next-to-Leading Logarithm (FONLL) calculations. The produced charm quarks scatter in the quark-gluon plasma (QGP) with the off-shell partons whose masses and widths are given by the Dynamical Quasi-Particle Model (DQPM) which reproduces the lattice QCD equation-of-state in thermal equilibrium. The relevant cross section are calculated in a consistent way by employing the effective propagators and couplings from the DQPM. Close to the critical energy density of the phase transition, the charm quarks are hadronized into D mesons through coalescence and/or fragmentation depending on transverse momentum. The hadronized D mesons then interact with the various hadrons in the hadronic phase with cross sections calculated in an effective lagrangian approach with heavy-quark spin symmetry. Finally, the nuclear modification factor \rm R_{AA} and the elliptic flow v_2 of D^0 mesons from PHSD are compared with the experimental data from the STAR Collaboration for Au+Au collisions at \sqrt{s_{\rm NN}} =200 GeV. We find that in the PHSD the energy loss of D mesons at high p_T can be dominantly attributed to partonic scattering while the actual shape of \rm R_{AA} versus p_T reflects the heavy quark hadronization scenario, i.e. coalescence versus fragmentation. Also the hadronic rescattering is important for the \rm R_{AA} at low p_T and enhances the D-meson elliptic flow v_2.

The evolution of the parton fraction of the total energy density.
The RAA of D0 mesons including partonic scattering with and wo. hadronic scattering

Effective QCD and transport description of dilepton and photon production in heavy-ion collisions and elementary processes

O. Linnyk, E. L. Bratkovskaya, W. Cassing · Prog.Part.Nucl.Phys. 87 (2016) 50-115

Show abstract

In this review we address the dynamics of relativistic heavy-ion reactions and in particular the information obtained from electromagnetic probes that stem from the partonic and hadronic phases. The out-of-equilibrium description of strongly interacting relativistic fields is based on the theory of Kadanoff and Baym. For the modeling of the partonic phase we introduce a dynamical quasiparticle model (DQPM) for QCD in equilibrium. The widths and masses of the quasiparticles are controlled by transport coefficients in comparison to lattice QCD results. The resulting off-shell transport approach - denoted by Parton-Hadron-String Dynamics (PHSD) - also includes covariant dynamical hadronization and keeps track of the hadronic interactions in the final phase. We show that PHSD captures the bulk dynamics of heavy-ion collisions from SPS to LHC energies and provides a basis for the evaluation of the electromagnetic emissivity, using the same dynamical parton propagators as for the system evolution. Direct photon production in elementary processes and heavy-ion reactions at RHIC and LHC energies is investigated and the status of the photon v2 puzzle - a large elliptic flow of the direct photons observed in A+A collisions - is addressed. We discuss the roles of hadronic and partonic sources for the photon spectra and the flow coefficients v2 and v3 and also the possibility to subtract the QGP signal from observables. Furthermore, the production of dilepton pairs is addressed from SIS to LHC energies. The low-mass dilepton yield is enhanced due to the in-medium modification of the rho-meson and at the lowest energy also due to a multiple regeneration of Delta-resonances. In addition, a signal of the partonic degrees-of-freedom is found in the intermediate dilepton mass regime (1.2GeV

The evolution of the parton fraction of the total energy density.
Spectrum of direct photons produced most central Au+Au collisions

In-medium effects in strangeness production in heavy-ion collisions at (sub)threshold energies

T. Song, L. Tolos, J. Wirth, J. Aichelin, and E. Bratkovskaya · Phys.Rev.C 103 (2021) 044901

Show abstract

We study the in-medium effects in strangeness production in heavy-ion collisions at (sub-)\-threshold energies of 1 - 2 A GeV based on the microscopic Parton-Hadron-String Dynamics (PHSD) transport approach. The in-medium modifications of the antikaon (\bar K = K^-, \bar K^0) properties are described via the self-consistent coupled-channel unitarized scheme based on a SU(3) chiral Lagrangian which incorporates explicitly the s- and p- waves of the kaon-nucleon interaction. This scheme provides the antikaon potential, spectral functions and reaction cross sections as well as their dependence on baryon density, temperature and antikaon momentum in the nuclear medium, which are incorporated in the off-shell dynamics of the PHSD. The in-medium modification of kaons (K = K^+, K^0) are accounted via the kaon-nuclear potential, which is assumed to be proportional to the local baryon density. The manifestation of the medium effects in observables is investigated for the K and {\bar K} rapidity distributions, p_T-spectra as well as the polar and azimuthal angular distributions, direct (v_1) and elliptic (v_2) flow in C+C, Ni+Ni, and Au+Au collisions. We find - by comparison to experimental data from the KaoS, FOPI and HADES Collaborations - that the modifications of (anti)kaon properties in nuclear matter are necessary to explain the data in a consistent manner. Moreover, we demonstrate the sensitivity of kaon observables to the equation-of-state of nuclear matter.

The evolution of the parton fraction of the total energy density.
The contour plot for the mass distributions of K− as a function of time

Electromagnetic emission from strongly interacting hadronic and partonic matter created in heavy-ion collisions

A. W. Romero Jorge, T. Song, Q. Zhou, E. Bratkovskaya · Phys. Rev. C 111 (2025) 064904

Show abstract

We investigate dilepton production in heavy-ion, proton-proton, and proton-nucleus collisions from low energies of 1 AGeV (SIS) to ultra-relativistic energies (LHC) using the Parton-Hadron-String Dynamics (PHSD) transport approach. PHSD is a microscopic, non-equilibrium approach that integrates hadronic and partonic degrees of freedom, providing a comprehensive description of relativistic heavy-ion collisions from initial nucleon-nucleon interactions to quark-gluon plasma (QGP) formation, hadronization, and final-state interactions. Key dilepton sources in PHSD include hadronic decays, bremsstrahlung, QGP radiation (q+\bar q \to e^+e^-, \ q+\bar q \to g+ e^+e^-, \ q+g \to q+ e^+e^-), primary Drell-Yan production, and semileptonic decays of correlated charm and bottom pairs. PHSD well describes dilepton data from HADES, STAR, and ALICE experiments. We examine in-medium effects, such as the vector meson spectral function broadening, and present the excitation function for the dilepton "excess" in the invariant mass range 0.4

Electromagnetic emission from strongly interacting hadronic and partonic matter created in heavy-ion collisions
Invariant mass spectra for Pb+Pb

Talks

Selected PHSD talks with direct slide links.

Electromagnetic emission from strongly interacting hadronic and partonic matter

Elena Bratkovskaya · The 11th International Symposium on Non-equilibrium Dynamics (NeD-2026), France · Feb-2026

Attenuation of jet partons and heavy quarks in strongly interacting QGP

Elena Bratkovskaya · Workshop High energy probes of the initial stages, CERN · Apr-2025

PHSD/PHQMD for FAIR

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

Constraining Dark Photon Kinetic Mixing in Heavy-Ion Collisions from SIS to LHC energies

A. W. Romero Jorge · Dark Matter Conference, Santander, Spain · Jun-2025

Combined constraints on dark photons from high-energy collisions, cosmology, and astrophysics

A. W. Romero Jorge · DMLab: DarkMatter@Bonn · Oct-2025

Charm dynamics in heavy-ion collisions

Elena Bratkovskaya · 2nd Workshop of Network NA7-HF-QGP (STRONG-2020) · Sep-2023

In-medium effects on hidden strangeness production in heavy-ion collisions

Taesoo Song · SQM, Korea · Jun-2022

The impact of electromagnetic and vortical fields in relativistic nuclear collisions

Lucia Oliva · Theoretical Physics Colloquium, Arizona State University (online) · 27 Jan 2021

Exploring the QGP at finite baryon chemical potential in and out of equilibrium

Elena Bratkovskaya · Theoretical Physics Colloquium, Arizona State University (online) · 20 Jan 2021

Transport coefficients of the dense QGP along the chiral phase transition

Olga Soloveva · Transport meeting (Frankfurt, Zoom) · 10 Dec 2020

Dynamical description of the partonic phase at finite chemical potential

Elena Bratkovskaya · NeD conference deck · 2019

Electromagnetic emissivity of hot and dense matter: QGP vs charm

Elena Bratkovskaya · 34th Winter Workshop on Nuclear Dynamics, Deshaies, Guadeloupe · 25-31 Mar 2018

Processes in collisions of heavy ions

Pierre Moreau · CBM Students Colloquium (29th CBM Week), GSI, Darmstadt · 19 Mar 2017

Observables of chiral symmetry restoration in heavy-ion collisions

Alessia Palmese · International School of Nuclear Physics, Erice, Italy · 21 Sep 2016

Chiral symmetry restoration versus deconfinement in heavy-ion collisions at high baryon density

Elena Bratkovskaya · CPOD 2016, Wroclaw, Poland · 30 May-4 Jun 2016

Tomography of the Quark-Gluon Plasma by charm quarks

Elena Bratkovskaya · Heavy Flavor Working Group Meeting, Berkeley Lab, USA · 18-22 Jan 2016

Electromagnetic emissivity of hot and dense matter

Elena Bratkovskaya · New Horizons in Fundamental Physics, Makutsi, South Africa · 23-28 Nov 2015

The QGP dynamics in relativistic heavy-ion collisions

Elena Bratkovskaya · Kruger2014 workshop, South Africa · 1-5 Dec 2014

From hadrons to partons and back

Elena Bratkovskaya · EMMI Nuclear and Quark Matter Seminar, GSI, Darmstadt · 9 Jul 2014