Researcher profile

Volodymyr Vovchenko

Volodymyr Vovchenko contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
25works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2026arXiv

Lattice-based equation of state with a critical point from constant entropy contours and its comparison to effective QCD approaches

In this work, we systematically assess the performance of a new method from [H. Shah et al., Phys. Rev. C 113, L012201] for locating the QCD critical point using constant-entropy contours by testing it against various effective QCD approaches. We demonstrate that, while the method yields spurious critical points in purely hadronic models (HRG) due to non-parabolic contour behavior at low temperatures ($T \lesssim 120$ MeV), it accurately reproduces the CP location in frameworks that feature a genuine phase transition and benchmarked against lattice QCD, such as Holographic Einstein-Maxwell-Dilaton, and Functional QCD approaches. Building on our previous determination of constant entropy contours using lattice data, we extend that analysis to construct a complete Lattice-based Equation of State (EoS) at finite density, which features a critical point at $(T, μ_B) \approx (114, 602)$ MeV. By integrating the extrapolated entropy density with respect to temperature, we reconstruct the pressure, baryon density, susceptibility, and speed of sound in the critical region, and analyze the focusing behavior of isentropic trajectories in the vicinity of the critical point.

preprint2023arXiv

Spinodal enhancement of fluctuations in nucleus-nucleus collisions

Subensemble Acceptance Method (SAM) [1,2] is an essential link between measured event-by-event fluctuations and their grand canonical theoretical predictions such as lattice QCD. The method allows quantifying the global conservation law effects in fluctuations. In its basic formulation, SAM requires a sufficiently large system such as created in central nucleus-nucleus collisions and sufficient space-momentum correlations. Directly in the spinodal region of the First Order Phase Transition (FOPT) different approximations should be used that account for finite size effects. Thus, we present the generalization of SAM applicable in both the pure phases, metastable and unstable regions of the phase diagram [3]. Obtained analytic formulas indicate the enhancement of fluctuations due to crossing the spinodal region of FOPT and are tested using molecular dynamics simulations. A rather good agreement is observed. Using transport model calculations with interaction potential we show that the spinodal enhancement of fluctuations survives till the later stages of collision via the memory effect [4]. However, at low collision energies the space-momentum correlation is not strong enough for this signal to be transferred to second and third order cumulants measured in momentum subspace. This result agrees well with recent HADES data on proton number fluctuations at $\sqrt{s_{NN}}=2.4$ GeV which are found to be consistent with the binomial momentum space acceptance [5].

preprint2022arXiv

Coalescence, the thermal model and multi-fragmentation: The energy and volume dependence of light nuclei production in heavy ion collisions

We present results of a phase space coalescence approach within the UrQMD transport and -hybrid model for a very wide range of beam energies from SIS to LHC. The coalescence model is able to qualitatively describe the whole range of experimental data with a fixed set of parameters. Some systematic deviations are observed for very low beam energies where the role of feed down from heavier nuclei and multi-fragmentation becomes relevant. The coalescence results are mostly very close to the thermal model fits. However, both the coalescence approach as well as thermal fits are struggling to simultaneously describe the triton multiplicities measured with the STAR and ALICE experiment. The double ratio of $tp/d^2$, in the coalescence approach, is found to be essentially energy and centrality independent for collisions of heavy nuclei at beam energies of $\mathrm{E_{lab}}> 10 A$ GeV. On the other hand the clear scaling of the $d/p^2$ and $t/p^3$ ratios with the systems volume is broken for peripheral collisions, where a canonical treatment and finite size effects become more important.

preprint2022arXiv

Constraining baryon annihilation in the hadronic phase of heavy-ion collisions via event-by-event fluctuations

We point out that the variance of net-baryon distribution normalized by the Skellam distribution baseline, $κ_2[B-\bar{B}]/\langle B+\bar{B}\rangle$, is sensitive to the possible modification of (anti)baryon yields due to $B\bar{B}$ annihilation in the hadronic phase. The corresponding measurements can thus place stringent limits on the magnitude of the $B\bar{B}$ annihilation and its inverse reaction. We perform Monte Carlo simulations of the hadronic phase in Pb-Pb collisions at the LHC via the recently developed subensemble sampler + UrQMD afterburner and show that the effect survives in net-proton fluctuations, which are directly accessible experimentally. The available experimental data of the ALICE Collaboration on net-proton fluctuations disfavors a notable suppression of (anti)baryon yields in $B\bar{B}$ annihilations predicted by the present version of UrQMD if only global baryon conservation is incorporated. On the other hand, the annihilations improve the data description when local baryon conservation is imposed. The two effects can be disentangled by measuring $κ_2[B+\bar{B}]/\langle B+\bar{B}\rangle$, which at the LHC is notably suppressed by annihilations but virtually unaffected by baryon number conservation.

preprint2022arXiv

Correcting event-by-event fluctuations in heavy-ion collisions for exact global conservation laws with the generalized subensemble acceptance method

We introduce the subensemble acceptance method 2.0 (SAM-2.0) -- a procedure to correct cumulants of a random number distribution inside a subsystem for the effect of exact global conservation of a conserved quantity to which this number is correlated, with applications to measurements of event-by-event fluctuations in heavy-ion collisions. The method expresses the corrected cumulants in terms of the cumulants inside and outside the subsystem that are not subject to the exact conservation. The derivation assumes that all probability distributions associated with the cumulants are peaked at the mean values but are otherwise of arbitrary shape. The formalism reduces to the original SAM [V. Vovchenko et al., Phys.Lett.B 811 (2020) 135868 [arXiv:2003.13905]] when applied to a coordinate space subvolume of a uniform thermal system. As the new method is restricted neither to the uniform systems nor to the coordinate space, it is applicable to fluctuations measured in heavy-ion collisions at various collision energies in different momentum space acceptances. The SAM-2.0 thus brings the experimental measurements and theoretical calculations of event-by-event fluctuations closer together, as the latter are typically performed without the account of exact global conservation.

preprint2022arXiv

Critical point particle number fluctuations from molecular dynamics

We study fluctuations of particle number in the presence of critical point by utilizing molecular dynamics simulations of the classical Lennard-Jones fluid in a periodic box. The numerical solution of the $N$-body problem naturally incorporates all correlations, exact conservation laws, and finite size effects, allowing us to study the fluctuation signatures of the critical point in a dynamical setup. We find that large fluctuations associated with the critical point are observed when measurements are performed in coordinate subspace, but, in the absence of collective flow and expansion, are essentially washed out when momentum cuts are imposed instead. We put our findings in the context of event-by-event fluctuations in heavy-ion collisions.

preprint2022arXiv

Proton number cumulants and correlation functions from hydrodynamics and the QCD phase diagram

We analyze the behavior of (net-)proton number cumulants in central collisions of heavy ions across a broad collision energy range by utilizing hydrodynamic simulations. The calculations incorporate essential non-critical contributions to proton fluctuations such as repulsive baryonic core and exact baryon number conservation. The experimental data are consistent with non-critical physics at collision energies $\sqrt{s_{\rm NN}} \geq 20$ GeV. The data from STAR and HADES Collaborations at lower collision energies indicate an excess of (multi-)proton correlations over the non-critical reference. This observation is discussed in the context of different mechanisms, including the possibility of a critical point in the baryon-rich region of the QCD phase diagram.

preprint2022arXiv

Proton number cumulants and correlation functions in Au-Au collisions at $\sqrt{s_{\rm NN}} = 7.7-200$ GeV from hydrodynamics

We present a dynamical description of (anti)proton number fluctuations cumulants and correlation functions in central Au-Au collisions at $\sqrt{s_{\rm NN}} = 7.7-200$ GeV by utilizing viscous hydrodynamics simulations. The cumulants of proton and baryon number are calculated in a given momentum acceptance analytically, via an appropriately extended Cooper-Frye procedure describing particlization of an interacting hadron resonance gas. The effects of global baryon number conservation are taken into account using a generalized subensemble acceptance method. The experimental data of the STAR collaboration are consistent at $\sqrt{s_{\rm NN}} \gtrsim 20$ GeV with simultaneous effects of global baryon number conservation and repulsive interactions in baryon sector, the latter being in line with the behavior of baryon number susceptibilities observed in lattice QCD. The data at lower collision energies show possible indications for sizable attractive interactions among baryons. The data also indicate sizable negative two-particle correlations between antiprotons that are not satisfactorily described by baryon conservation and excluded volume effects. We also discuss differences between cumulants and correlation functions (factorial cumulants) of (anti)proton number distribution, proton versus baryon number fluctuations, and effects of hadronic afterburner.

preprint2022arXiv

Thermal-model-based characterization of heavy-ion-collision systems at chemical freeze-out

We investigate the chemical freeze-out in heavy-ion collisions (HICs) and the impact of the hadronic spectrum on thermal model analyses. Detailed knowledge of the hadronic spectrum is still an open question, which has phenomenological consequences on the study of HICs. By varying the number of resonances included in Hadron Resonance Gas (HRG) Model calculations, we can shed light on which particles may be produced. Furthermore, we study the influence of the number of states on the so-called two flavor freeze-out scenario, in which strange and light particles can freeze-out separately. We consider results for the chemical freeze-out parameters obtained from thermal model fits and from calculating net-particle fluctuations. We will show the effect of using one global temperature to fit all particles and alternatively, allowing particles with and without strange quarks to freeze-out separately.

preprint2022arXiv

Thermodynamic approach to proton number fluctuations in baryon-rich heavy-ion matter created at moderate collision energies

We develop a framework to relate proton number cumulants measured in heavy-ion collisions within a momentum space acceptance to the susceptibilities of baryon number, assuming that particles are emitted from a fireball with uniform distribution of temperature and baryochemical potential, superimposed on a hydrodynamic flow velocity profile. The rapidity acceptance dependence of proton cumulants measured by the HADES Collaboration in $\sqrt{s_{\rm NN}} = 2.4$ GeV Au-Au appears to be consistent with thermal emission of nucleons from a grand-canonical heat bath, with the extracted baryon number susceptibilities exhibiting an hierarchy $χ_4^B \gg -χ_3^B \gg χ_2^B \gg χ_1^B$. Naively, this could indicate large non-Gaussian fluctuations that might point to the presence of the QCD critical point close to the chemical freeze-out at $T \sim 70$ MeV and $μ_B \sim 850-900$ MeV. However, the description of the experimental data at large rapidity acceptances becomes challenging once the effect of exact baryon number conservation is incorporated, suggesting that more theoretical and experimental studies are needed to reach a firm conclusion.

preprint2021arXiv

Higher order conserved charge fluctuations inside the mixed phase

General formulas are presented for higher order cumulants of the conserved charge statistical fluctuations inside the mixed phase. As a particular example the van der Waals model in the grand canonical ensemble is used. The higher order measures of the conserved charge fluctuations up to the hyperkurtosis are calculated in a vicinity of the critical point (CP). The analysis includes both the mixed phase region and the pure phases on the phase diagram. It is shown that even-order fluctuation measures, e.g. scaled variance, kurtosis, and hyperkurtosis, have only positive values in the mixed phase, and go to infinity at the CP. For odd-order measures, such as skewness and hyperskewness, the regions of positive and negative values are found near the left and right binodals, respectively. The obtained results are discussed in a context of the event-by-event fluctuation measurements in heavy-ion collisions.

preprint2021arXiv

Solving the puzzle of high temperature light (anti)-nuclei production in ultra-relativistic heavy ion collisions

The creation of loosely bound objects in heavy ion collisions, e.g. light clusters, near the phase transition temperature $\left(T_{\rm ch} \approx 155 \, \rm{MeV} \right)$ has been a puzzling observation that seems to be at odds with Big Bang nucleosynthesis suggesting that deuterons and other clusters are formed only below a temperature $T\approx 0.1-1 \, \text{MeV}$. We solve this puzzle by showing that the light cluster abundancies in heavy ion reactions stay approximately constant from chemical freeze-out to kinetic freeze-out. To this aim we develop an extensive network of coupled reaction rate equations including stable hadrons and hadronic resonances to describe the temporal evolution of the abundancies of light (anti-)(hyper-)nuclei in the late hadronic environment of an ultrarelativistic heavy ion collision. It is demonstrated that the chemical equilibration of the light nuclei occurs on a very short timescale as a consequence of the strong production and dissociation processes. However, because of the partial chemical equilibrium of the stable hadrons, including the nucleon feeding from $Δ$ resonances, the abundancies of the light nuclei stay nearly constant during the evolution and cooling of the hadronic phase. This solves the longstanding contradiction between the thermal fits and the late stage coalescence (and the Big Bang nucleosynthesis) and explains why the observed light cluster yields are compatible with both a high chemical production temperature and a late state emission as modelled by coalescence. We also note in passing that the abundancies of the light clusters in the present approach are in excellent agreement with those measured by ALICE at LHC.

preprint2020arXiv

Binomial acceptance corrections for particle number distributions in high-energy reactions

The binomial acceptance correction procedure is studied for particle number distributions detected in high energy reactions in finite regions of the momentum space. We present acceptance correction formulas for scaled variance, skewness, and kurtosis. Our considerations include various specific types of particles - positively or negatively charged, baryons and antibaryons - as well as conserved charges, namely, the net baryon number and electric charge. A simple model with effects of exact charge conservation, namely the Bessel distribution, is studied in some detail where effects of multi-particle correlations are present. The accuracy of the binomial filter is studied with UrQMD model simulations of inelastic proton-proton reactions. Binomial acceptance correction procedure works well when used inside a small region of phase space as well as for certain efficiency corrections, in particular for constructing net proton fluctuation from net baryon ones. Its performance is less accurate when applied to obtain UrQMD fluctuations in a finite rapidity window from fluctuations in the full $4π$ space.

preprint2020arXiv

Bose-Einstein condensation phenomenology in systems with repulsive interactions

The role of repulsive interactions in statistical systems of Bose particles is investigated. Three different phenomenological frameworks are considered: a mean field model, an excluded volume model, and a model with a medium dependent effective mass. All three models are tuned to yield similar equations of state, with only minor deviations from the ideal Bose gas at small chemical potentials. Our analysis indicates, however, that these models lead to qualitatively different results for the Bose-Einstein condensation phenomenon. We discuss the different aspects of this phenomenon, namely, an onset of the Bose-Einstein condensation, particle number fluctuations, and a behavior of the Bose condensate. The obtained results can be helpful for interpreting the lattice QCD data at small temperature and large isospin chemical potential and the data on multiple pion production in high energy nuclear collisions.

preprint2020arXiv

Critical point fluctuations: Finite size and global charge conservation effects

We investigate simultaneous effects of finite system size and global charge conservation on thermal fluctuations in the vicinity of a critical point. For that we consider a finite interacting system which exchanges particles with a finite reservoir (thermostat), comprising a statistical ensemble that is distinct from the common canonical and grand canonical ensembles. As a particular example the van der Waals model is used. The global charge conservation effects strongly influence the cumulants of particle number distribution when the system size is comparable to that of the reservoir. If the system size is large enough to capture all the physics associated with the interactions, the global charge conservation effects can be accurately described and corrected for analytically, within a recently developed subensemble acceptance method. The finite size effects start to play a significant role when the correlation length grows large due to proximity of the critical point or when the system is small enough to be comparable to an eigenvolume of an individual particle. We discuss our results in the context of fluctuation measurements in heavy-ion collisions.

preprint2020arXiv

Feeddown contributions from unstable nuclei in relativistic heavy-ion collisions

We estimate the feeddown contributions from decays of unstable $A = 4$ and $A = 5$ nuclei to the final yields of protons, deuterons, tritons, $^3$He, and $^4$He produced in relativistic heavy-ion collisions at $\sqrt{s_{\rm NN}} > 2.4$ GeV, using the statistical model. The feeddown contribution effects do not exceed 5% at LHC and top RHIC energies due to the large penalty factors involved, but are substantial at intermediate collision energies. We observe large feeddown contributions for tritons, $^3$He, and $^4$He at $\sqrt{s_{\rm NN}} \lesssim 10$ GeV, where they may account for as much as 70% of the final yield at the lower end of the collision energies considered. Sizable ($>10$%) effects for deuteron yields are observed at $\sqrt{s_{\rm NN}} \lesssim 4$ GeV. The results suggest that the excited nuclei feeddown cannot be neglected in the ongoing and future analysis of light nuclei production at intermediate collision energies, including HADES and CBM experiments at FAIR, NICA at JINR, RHIC beam energy scan and fixed-target programmes, and NA61/SHINE at CERN. We further show that the freeze-out curve in the $T$-$μ_B$ plane itself is affected significantly by the light nuclei at high baryochemical potential.

preprint2020arXiv

Hadron resonance gas with van der Waals interactions

An overview of a hadron resonance gas (HRG) model that includes van der Waals interactions between hadrons is presented. Applications of the excluded volume HRG model to heavy-ion collision data and lattice QCD equation of state are discussed. A recently developed quantum van der Waals hadron resonance gas model is covered as well. Applications of this model in the context of the QCD critical point are elaborated.

preprint2020arXiv

Kinetic freeze-out temperature from yields of short-lived resonances

A method to determine the kinetic freeze-out temperature in heavy-ion collisions from measured yields of short-lived resonances is presented. The resonance production is treated in the framework of thermal model with an evolution between chemical and kinetic freeze-outs. The yields of many short-lived resonances are suppressed at $T = T_{\rm kin} < T_{\rm ch}$. We determine the values of $T_{\rm kin}$ and $T_{\rm ch}$ for various centralities in Pb--Pb collisions at $\sqrt{s_{_{NN}}} = 2.76$ TeV by fitting the abundances of both the stable hadrons and the short-lived resonances such as $ρ^0$ and $ \text{K}^{*0}$, that were measured by the ALICE collaboration. This allows to extract the kinetic freeze-out temperature from the measured hadron and resonance yields alone, independent of assumptions about the flow velocity profile and the freeze-out hypersurface. The extracted $T_{\rm ch}$ values exhibit a moderate multiplicity dependence whereas $T_{\rm kin}$ drops, from $T_{\rm kin} \simeq T_{\rm ch} \simeq 155$ MeV in peripheral collisions to $T_{\rm kin} \simeq 110$ MeV in 0-20% central collisions. Predictions for other short-lived resonances are presented. A potential (non-)observation of a suppressed $f_0(980)$ meson yield will allow to constrain the lifetime of that meson.

preprint2020arXiv

Pion condensation in the early Universe at nonvanishing lepton flavor asymmetry and its gravitational wave signatures

We investigate the possible formation of a Bose-Einstein condensed phase of pions in the early Universe at nonvanishing values of lepton flavor asymmetries. A hadron resonance gas model with pion interactions, based on first-principle lattice QCD simulations at nonzero isospin density, is used to evaluate cosmic trajectories at various values of electron, muon, and tau lepton asymmetries that satisfy the available constraints on the total lepton asymmetry. The cosmic trajectory can pass through the pion condensed phase if the combined electron and muon asymmetry is sufficiently large: $|l_e + l_μ| \gtrsim 0.1$, with little sensitivity to the difference $l_e - l_μ$ between the individual flavor asymmetries. Future constraints on the values of the individual lepton flavor asymmetries will thus be able to either confirm or rule out the condensation of pions during the cosmic QCD epoch. We demonstrate that the pion condensed phase leaves an imprint both on the spectrum of primordial gravitational waves and on the mass distribution of primordial black holes at the QCD scale e.g. the black hole binary of recent LIGO event GW190521 can be formed in that phase.

preprint2020arXiv

QCD equation of state at vanishing and high baryon density: Chiral Mean Field model

The thermodynamic properties of high temperature and high density QCD-matter are studied using the Chiral SU(3)-flavor parity-doublet Polyakov-loop quark-hadron mean-field model, CMF. The CMF model provides a proper description of lattice QCD data, heavy-ions physics, and static neutron stars. The behavior of lines of constant pressure with increase of baryon density is discussed. The rapid change of pressure behavior at $μ_B/T\approx3$ suggests a strong contribution of baryons to thermodynamic properties at this region. The position of this region is very close to the radius of convergence for a Taylor expansion of the QCD pressure. The role of mesons and unstable hadrons in the hydrodynamic expansion of strongly interacting matter is also discussed.

preprint2020arXiv

Review of opportunities at the HESR-Collider with PANDA at FAIR

Exciting new scientific opportunities are presented for the PANDA detector at the High Energy Storage Ring in the redefined $\overline{\text{p}} \text{p}(A)$ collider mode, HESR-C, at the Facility for Antiproton and Ion Research (FAIR) in Europe. The high luminosity, $L \sim 10^{31}$ cm$^{-2}$ s$^{-1}$, and a wide range of intermediate and high energies, $\sqrt{s_{\text{NN}}}$ up to 30 GeV for $\overline{\text{p}} \text{p}(A)$ collisions will allow to explore a wide range of exciting topics in QCD, including the study of the production of excited open charm and bottom states, nuclear bound states containing heavy (anti)quarks, the interplay of hard and soft physics in the dilepton production, probing short-range correlations in nuclei, and the exploration of the early, complete $\overline{\text{p}}$-p-annihilation phase, where an intially pure Yang-Mills gluon plasma is formed.

preprint2020arXiv

Traces of the nuclear liquid-gas phase transition in the analytic properties of hot QCD

The nuclear liquid-gas transition at normal nuclear densities, $n \sim n_0 = 0.16$ fm$^{-3}$, and small temperatures, $T \sim 20$ MeV, has a large influence on analytic properties of the QCD grand-canonical thermodynamic potential. A classical van der Waals equation is used to determine these unexpected features due to dense cold matter qualitatively. The existence of the nuclear matter critical point results in thermodynamic branch points, which are located at complex chemical potential values, for $T > T_c \simeq 20$ MeV, and exhibit a moderate model dependence up to rather large temperatures $T \lesssim 100$ MeV. The behavior at higher temperatures is studied using the van der Waals hadron resonance gas (vdW-HRG) model. The baryon-baryon interactions have a decisive influence on the QCD thermodynamics close to $μ_B = 0$. In particular, nuclear matter singularities limit the radius of convergence $r_{μ_B/T}$ of the Taylor expansion in $μ_B/T$, with $r_{μ_B/T} \sim 2-3$ values at $T \sim 140-170$ MeV obtained in the vdW-HRG model.

preprint2019arXiv

Critical point signatures in the cluster expansion in fugacities

The QCD baryon number density can formally be expanded into a Laurent series in fugacity, which is a relativistic generalization of Mayer&#39;s cluster expansion. We determine properties of the cluster expansion in a model with a phase transition and a critical point at finite baryon density, in which the Fourier coefficients of the expansion can be determined explicitly and to arbitrary order. The asymptotic behavior of Fourier coefficients changes qualitatively as one traverses the critical temperature and it is connected to the branch points of a thermodynamic potential associated with the phase transition. The results are discussed in the context of lattice QCD simulations at imaginary chemical potential. We argue that the location of a branch point closest to the imaginary chemical potential axis can be extracted through an analysis of an exponential suppression of Fourier coefficients. This is illustrated using the four leading coefficients both in a toy model as well as by using recent lattice QCD data.

preprint2019arXiv

Equation of state for hot QCD and compact stars from a mean field approach

The thermodynamic properties of high temperature and high density QCD-matter are explored within the Chiral SU(3)-flavor parity-doublet Polyakov-loop quark-hadron mean-field model, CMF. The quark sector of the CMF model is tuned to describe the $μ_B=0$ thermodynamics data of lattice QCD. The resulting lines of constant physical variables as well as the baryon number susceptibilities are studied in some detail in the temperature/chemical potential plane. The CMF model predicts three consecutive transitions, the nuclear first-order liquid-vapor phase transition, chiral symmetry restoration, and the cross-over transition to a quark-dominated phase. All three phenomena are cross-over, for most of the $T-μ_B$-plane. The deviations from the free ideal hadron gas baseline at $μ_B=0$ and $T\approx 100-200$ MeV can be attributed to remnants of the liquid-vapor first order phase transition in nuclear matter. The chiral crossing transition determines the baryon fluctuations at much higher $μ_B\approx1.5$ GeV, and at even higher baryon densities $μ_B\approx2.4$ GeV, the behavior of fluctuations is controlled by the deconfinement cross-over. The CMF model also describe well the static properties of high $μ_B$ neutron stars as well as the new neutron star merger observations. The effective EoS presented here describes simultaneously lattice QCD results at $μ_B=0$, as well as observed physical phenomena (nuclear matter and neutron star matter) at $T\cong0$ and high densities, $μ_B>1$ GeV.

preprint2015arXiv

Glueballs amass at RHIC and LHC Colliders! - The early quarkless 1st order phase transition at $T=270$ MeV - from pure Yang-Mills glue plasma to GlueBall-Hagedorn states

The early stage of high multiplicity pp, pA and AA collider is represented by a nearly quarkless, hot, deconfined pure gluon plasma. According to pure Yang-Mills Lattice Gauge Theory, this hot pure glue matter undergoes, at a high temperature, $T_c = 270$ MeV, a first order phase transition into a confined Hagedorn-GlueBall fluid. These new scenario should be characterized by a suppression of high $p_T$ photons and dileptons, baryon suppression and enhanced strange meson production. We propose to observe this newly predicted class of events at LHC and RHIC.