Showing 1–50 of 116 results for author: Koch, V

Coordinate versus momentum cuts and effects of collective flow on critical fluctuations

Authors: Volodymyr A. Kuznietsov, Mark I. Gorenstein, Volker Koch, Volodymyr Vovchenko

Abstract: We analyze particle number fluctuations in the crossover region near the critical endpoint of a first-order phase transition by utilizing molecular dynamics simulations of the classical Lennard-Jones fluid. We extend our previous study [V.A. Kuznietsov et al., Phys. Rev. C 105, 044903 (2022)] by incorporating longitudinal collective flow. The scaled variance of particle number distribution inside… ▽ More We analyze particle number fluctuations in the crossover region near the critical endpoint of a first-order phase transition by utilizing molecular dynamics simulations of the classical Lennard-Jones fluid. We extend our previous study [V.A. Kuznietsov et al., Phys. Rev. C 105, 044903 (2022)] by incorporating longitudinal collective flow. The scaled variance of particle number distribution inside different coordinate and momentum space acceptances is computed through ensemble averaging and found to agree with earlier results obtained using time averaging, validating the ergodic hypothesis for fluctuation observables. Presence of a sizable collective flow is found to be essential for observing large fluctuations from the critical point in momentum space acceptances. We discuss our findings in the context of heavy-ion collisions. △ Less

Submitted 30 March, 2024; originally announced April 2024.

Comments: 12 pages, 10 figures

Might Normal Nuclear Matter be Quarkyonic?

Authors: Volker Koch, Larry McLerran, Gerald A. Miller, Volodymyr Vovchenko

Abstract: The possibility that nuclear matter might be Quarkyonic is considered. Quarkyonic matter is high baryon density matter that is confined but can be approximately thought of as a filled Fermi sea of quarks surrounded by a shell of nucleons. Here, nuclear matter is described by the IdylliQ sigma model for Quarkyonic matter, generalizing the non-interacting IdylliQ model [Y. Fujimoto et al., Phys. Rev… ▽ More The possibility that nuclear matter might be Quarkyonic is considered. Quarkyonic matter is high baryon density matter that is confined but can be approximately thought of as a filled Fermi sea of quarks surrounded by a shell of nucleons. Here, nuclear matter is described by the IdylliQ sigma model for Quarkyonic matter, generalizing the non-interacting IdylliQ model [Y. Fujimoto et al., Phys. Rev. Lett. 132, 112701 (2024) [arXiv:2306.04304]] to include interactions with a sigma meson and a pion. When such interactions are included, we find that isospin-symmetric nuclear matter binds, with acceptable values of the compressibility and other parameters for nuclear matter at saturation. The energy per nucleon and sound velocity of such matter is computed, and the isospin dependence is determined. Nuclear matter is formed at a density close to but slightly above the density at which Quarkyonic matter forms. Quarkyonic matter predicts a strong depletion of nucleons in normal nuclear matter at low momentum. Such a depletion for nucleon momenta $k \lesssim 120$ MeV is shown to be consistent with electron scattering data. △ Less

Submitted 22 March, 2024; originally announced March 2024.

Comments: 23 pages, 6 figures

Controlling volume fluctuations for studies of critical phenomena in nuclear collisions

Authors: Romain Holzmann, Volker Koch, Anar Rustamov, Joachim Stroth

Abstract: We generalize and extend the recently proposed method to account for contributions of system size (or volume/participant) fluctuations to the experimentally measured moments of particle multiplicity distributions. We find that in the general case there are additional biases which are not directly accessible to experiment. These biases are, however, parametrically suppressed if the multiplicity of… ▽ More We generalize and extend the recently proposed method to account for contributions of system size (or volume/participant) fluctuations to the experimentally measured moments of particle multiplicity distributions. We find that in the general case there are additional biases which are not directly accessible to experiment. These biases are, however, parametrically suppressed if the multiplicity of the particles of interest is small compared to the total charged-particle multiplicity, e.g., in the case of proton number fluctuations at top RHIC and LHC energies. They are also small if the multiplicity distribution of charged particles per wounded nucleon is close to the Poissonian limit, which is the case at low energy nuclear collisions, e.g., at GSI/SIS18. We further find that mixed events are not necessarily needed to extract the correction for volume fluctuations, albeit it can help if event statistics is small, which is typically the case for reconstructing the higher-order cumulants. We provide the formulas to correct pure and mixed cumulants of particle multiplicity distributions up to any order together with their associated biases. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Interplay of baryonic chiral partners in fluctuations of net-baryon number density

Authors: Michał Marczenko, Volker Koch, Krzysztof Redlich, Chihiro Sasaki

Abstract: In this contribution, we use the parity doublet model to investigate the fluctuations of the net-baryon number density. We discuss the systematics of the susceptibilities and their ratios for nucleons of positive and negative parity, as well as their correlator. We demonstrate that the fluctuations of positive-parity nucleon do not reflect the fluctuations of the total net-baryon number at the chi… ▽ More In this contribution, we use the parity doublet model to investigate the fluctuations of the net-baryon number density. We discuss the systematics of the susceptibilities and their ratios for nucleons of positive and negative parity, as well as their correlator. We demonstrate that the fluctuations of positive-parity nucleon do not reflect the fluctuations of the total net-baryon number at the chiral phase transition. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: prepared for the Quark Matter 2023 proceedings

Fluctuations and correlations of baryonic chiral partners

Authors: Volker Koch, Michał Marczenko, Krzysztof Redlich, Chihiro Sasaki

Abstract: Fluctuations and correlations of the net-baryon number play an important role in exploring critical phenomena in phase transitions of strongly interacting matter governed by Quantum chromodynamics (QCD). In this work, we use the parity doublet model to investigate the fluctuations of the net-baryon number density in hot and dense hadronic matter. The model accounts for chiral criticality within th… ▽ More Fluctuations and correlations of the net-baryon number play an important role in exploring critical phenomena in phase transitions of strongly interacting matter governed by Quantum chromodynamics (QCD). In this work, we use the parity doublet model to investigate the fluctuations of the net-baryon number density in hot and dense hadronic matter. The model accounts for chiral criticality within the mean-field approximation. We focus on the qualitative properties and systematics of the first- and second-order susceptibility of the net-baryon number density, and their ratios for nucleons of positive and negative parity, as well as their correlator. We show that the fluctuations of the positive-parity nucleon do not necessarily reflect the fluctuations of the total net-baryon number density at the phase boundary of the chiral phase transition. We also investigate the non-trivial structure of the correlator. Furthermore, we discuss and quantify the differences between the fluctuations of the net-baryon number density in the vicinity of the chiral and liquid-gas phase transition in nuclear matter. We indicate a possible relevance of our results with the interpretation of the experimental data on net-proton number fluctuations in heavy-ion collisions. △ Less

Submitted 29 January, 2024; v1 submitted 30 August, 2023; originally announced August 2023.

Comments: Matches published version

Journal ref: Phys. Rev. D 109, 014033, 2024

Hot QCD White Paper

Authors: M. Arslandok, S. A. Bass, A. A. Baty, I. Bautista, C. Beattie, F. Becattini, R. Bellwied, Y. Berdnikov, A. Berdnikov, J. Bielcik, J. T. Blair, F. Bock, B. Boimska, H. Bossi, H. Caines, Y. Chen, Y. -T. Chien, M. Chiu, M. E. Connors, M. Csanád, C. L. da Silva, A. P. Dash, G. David, K. Dehmelt, V. Dexheimer , et al. (149 additional authors not shown)

Abstract: Hot QCD physics studies the nuclear strong force under extreme temperature and densities. Experimentally these conditions are achieved via high-energy collisions of heavy ions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). In the past decade, a unique and substantial suite of data was collected at RHIC and the LHC, probing hydrodynamics at the nucleon scale, the… ▽ More Hot QCD physics studies the nuclear strong force under extreme temperature and densities. Experimentally these conditions are achieved via high-energy collisions of heavy ions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). In the past decade, a unique and substantial suite of data was collected at RHIC and the LHC, probing hydrodynamics at the nucleon scale, the temperature dependence of the transport properties of quark-gluon plasma, the phase diagram of nuclear matter, the interaction of quarks and gluons at different scales and much more. This document, as part of the 2023 nuclear science long range planning process, was written to review the progress in hot QCD since the 2015 Long Range Plan for Nuclear Science, as well as highlight the realization of previous recommendations, and present opportunities for the next decade, building on the accomplishments and investments made in theoretical developments and the construction of new detectors. Furthermore, this document provides additional context to support the recommendations voted on at the Joint Hot and Cold QCD Town Hall Meeting, which are reported in a separate document. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Comments: 190 pages, 69 figures

The Present and Future of QCD

Authors: P. Achenbach, D. Adhikari, A. Afanasev, F. Afzal, C. A. Aidala, A. Al-bataineh, D. K. Almaalol, M. Amaryan, D. Androić, W. R. Armstrong, M. Arratia, J. Arrington, A. Asaturyan, E. C. Aschenauer, H. Atac, H. Avakian, T. Averett, C. Ayerbe Gayoso, X. Bai, K. N. Barish, N. Barnea, G. Basar, M. Battaglieri, A. A. Baty, I. Bautista , et al. (378 additional authors not shown)

Abstract: This White Paper presents the community inputs and scientific conclusions from the Hot and Cold QCD Town Meeting that took place September 23-25, 2022 at MIT, as part of the Nuclear Science Advisory Committee (NSAC) 2023 Long Range Planning process. A total of 424 physicists registered for the meeting. The meeting highlighted progress in Quantum Chromodynamics (QCD) nuclear physics since the 2015… ▽ More This White Paper presents the community inputs and scientific conclusions from the Hot and Cold QCD Town Meeting that took place September 23-25, 2022 at MIT, as part of the Nuclear Science Advisory Committee (NSAC) 2023 Long Range Planning process. A total of 424 physicists registered for the meeting. The meeting highlighted progress in Quantum Chromodynamics (QCD) nuclear physics since the 2015 LRP (LRP15) and identified key questions and plausible paths to obtaining answers to those questions, defining priorities for our research over the coming decade. In defining the priority of outstanding physics opportunities for the future, both prospects for the short (~ 5 years) and longer term (5-10 years and beyond) are identified together with the facilities, personnel and other resources needed to maximize the discovery potential and maintain United States leadership in QCD physics worldwide. This White Paper is organized as follows: In the Executive Summary, we detail the Recommendations and Initiatives that were presented and discussed at the Town Meeting, and their supporting rationales. Section 2 highlights major progress and accomplishments of the past seven years. It is followed, in Section 3, by an overview of the physics opportunities for the immediate future, and in relation with the next QCD frontier: the EIC. Section 4 provides an overview of the physics motivations and goals associated with the EIC. Section 5 is devoted to the workforce development and support of diversity, equity and inclusion. This is followed by a dedicated section on computing in Section 6. Section 7 describes the national need for nuclear data science and the relevance to QCD research. △ Less

Submitted 4 March, 2023; originally announced March 2023.

Comments: QCD Town Meeting White Paper, as submitted to 2023 NSAC LRP committee on Feb. 28, 2023

Journal ref: Nucl.Phys.A 1047 (2024) 122874

Quarkyonic or baryquark matter? On the dynamical generation of momentum space shell structure

Authors: Volker Koch, Volodymyr Vovchenko

Abstract: We study the equation of state of a mixture of (quasi-)free constituent quarks and nucleons with hard-core repulsion at zero temperature. Two opposite scenarios for the realization of the Pauli exclusion principle are considered: (i) a Fermi sea of quarks surrounded by a shell of baryons -- the quarkyonic matter, and (ii) a Fermi sea of nucleons surrounded by a shell of quarks which we call \emph{… ▽ More We study the equation of state of a mixture of (quasi-)free constituent quarks and nucleons with hard-core repulsion at zero temperature. Two opposite scenarios for the realization of the Pauli exclusion principle are considered: (i) a Fermi sea of quarks surrounded by a shell of baryons -- the quarkyonic matter, and (ii) a Fermi sea of nucleons surrounded by a shell of quarks which we call \emph{baryquark matter}. In both scenarios, the sizes of the Fermi sea and shell are fixed through energy minimization at fixed baryon number density. While both cases yield a qualitatively similar transition from hadronic to quark matter, we find that baryquark matter is energetically favored in this setup and yields a physically acceptable behavior of the speed of sound without the need to introduce an infrared regulator. In order to retain the theoretically more appealing quarkyonic matter as the preferred form of dense QCD matter will thus require modifications to the existing dynamical generation mechanisms, such as, for example, the introduction of momentum-dependent nuclear interactions. △ Less

Submitted 8 May, 2023; v1 submitted 26 November, 2022; originally announced November 2022.

Comments: 6 pages, 3 figures, to be published in Physics Letters B

Journal ref: Phys. Lett. B 841, 137942 (2023)

Constraints on hadron resonance gas interactions via first-principles Lattice QCD susceptibilities

Authors: Jamie M. Karthein, Volker Koch, Claudia Ratti, Volodymyr Vovchenko

Abstract: We investigate extensions of the Hadron Resonance Gas (HRG) Model beyond the ideal case by incorporating both attractive and repulsive interactions into the model. When considering additional states exceeding those measured with high confidence by the Particle Data Group, attractive corrections to the overall pressure in the HRG model are imposed. On the other hand, we also apply excluded-volume c… ▽ More We investigate extensions of the Hadron Resonance Gas (HRG) Model beyond the ideal case by incorporating both attractive and repulsive interactions into the model. When considering additional states exceeding those measured with high confidence by the Particle Data Group, attractive corrections to the overall pressure in the HRG model are imposed. On the other hand, we also apply excluded-volume corrections, which ensure there is no overlap of baryons by turning on repulsive (anti)baryon-(anti)baryon interactions. We emphasize the complementary nature of these two extensions and identify combinations of conserved charge susceptibilities that allow us to constrain them separately. In particular, we find interesting ratios of susceptibilities that are sensitive to one correction and not the other. This allows us to constrain the excluded volume and particle spectrum effects separately. Analysis of the available lattice results suggests the presence of both the extra states in the baryon-strangeness sector and the repulsive baryonic interaction, with indications that hyperons have a smaller repulsive core than non-strange baryons. We note that these results are interesting for heavy-ion-collision systems at both the LHC and RHIC. △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: Contribution to the 20th International Conference on Strangeness in Quark Matter (SQM2022)

Centrality dependence of proton and light nuclei yields as a consequence of baryon annihilation in the hadronic phase

Authors: Volodymyr Vovchenko, Volker Koch

Abstract: The centrality dependence of the $p/πぱい$ ratio measured by the ALICE Collaboration in 5.02 TeV Pb-Pb collisions indicates a statistically significant suppression with the increase of the charged particle multiplicity once the centrality-correlated part of the systematic uncertainty is eliminated from the data. We argue that this behavior can be attributed to baryon annihilation in the hadronic phase… ▽ More The centrality dependence of the $p/πぱい$ ratio measured by the ALICE Collaboration in 5.02 TeV Pb-Pb collisions indicates a statistically significant suppression with the increase of the charged particle multiplicity once the centrality-correlated part of the systematic uncertainty is eliminated from the data. We argue that this behavior can be attributed to baryon annihilation in the hadronic phase. By implementing the $B\bar{B} \leftrightarrow 5πぱい$ reaction within a generalized partial chemical equilibrium framework, we estimate the annihilation freeze-out temperature at different centralities, which decreases with increasing charged multiplicity and yields $T_{\rm ann} = 132 \pm 5$ MeV in 0-5% most central collisions. This value is considerably below the hadronization temperature of $T_{\rm had} \sim 160$ MeV but above the thermal (kinetic) freeze-out temperature of $T_{\rm kin} \sim 100$ MeV. Baryon annihilation reactions thus remain relevant in the initial stage of the hadronic phase but freeze out before (pseudo-)elastic hadronic scatterings. One experimentally testable consequence of this picture is a suppression of various light nuclei to proton ratios in central collisions of heavy ions. △ Less

Submitted 21 November, 2022; v1 submitted 27 October, 2022; originally announced October 2022.

Comments: 7 pages, 5 figures, to be published in Physics Letters B

Journal ref: Physics Letters B 835, 137577 (2022)

QCD Phase Structure and Interactions at High Baryon Density: Continuation of BES Physics Program with CBM at FAIR

Authors: D. Almaalol, M. Hippert, J. Noronha-Hostler, J. Noronha, E. Speranza, G. Basar, S. Bass, D. Cebra, V. Dexheimer, D. Keane, S. Radhakrishnan, A. I. Sheikh, M. Strickland, C. Y. Tsang, . X. Dong, V. Koch, G. Odyniec, N. Xu, F. Geurts, D. Hofman, M. Stephanov, G. Wilks, Z. Y. Ye, H. Z. Huang, G. Wang , et al. (19 additional authors not shown)

Abstract: We advocate for an active US participation in the international collaboration of the CBM experiment that will allow the US nuclear physics program to build on its successful exploration of the QCD phase diagram, use the expertise gained at RHIC to make complementary measurements at FAIR, and contribute to achieving the scientific goals of the beam energy scan (BES) program. We advocate for an active US participation in the international collaboration of the CBM experiment that will allow the US nuclear physics program to build on its successful exploration of the QCD phase diagram, use the expertise gained at RHIC to make complementary measurements at FAIR, and contribute to achieving the scientific goals of the beam energy scan (BES) program. △ Less

Submitted 21 December, 2022; v1 submitted 11 September, 2022; originally announced September 2022.

Comments: 30 pages, 7 figures

Measuring the speed of sound using cumulants of baryon number

Authors: Agnieszka Sorensen, Dmytro Oliinychenko, Larry McLerran, Volker Koch

Abstract: We show that the values of the first three cumulants of the baryon number distribution can be used to calculate the isothermal speed of sound and its logarithmic derivative with respect to the baryon number density. We discuss applications of this result to heavy-ion collision experiments and address possible challenges, including effects due to baryon number conservation, differences between prot… ▽ More We show that the values of the first three cumulants of the baryon number distribution can be used to calculate the isothermal speed of sound and its logarithmic derivative with respect to the baryon number density. We discuss applications of this result to heavy-ion collision experiments and address possible challenges, including effects due to baryon number conservation, differences between proton and baryon cumulants, and the influence of finite number statistics on fluctuation observables in both experiment and hadronic transport simulations. In particular, we investigate the relation between quantities calculated in infinite, continuous matter and observables obtained in simulations using a finite number of particles. △ Less

Submitted 11 September, 2022; originally announced September 2022.

Comments: 6 pages, 1 figure, presented at the XXIXth International Conference on Ultra-relativistic Nucleus-Nucleus Collisions (Quark Matter 2022)

Report number: INT-PUB-22-030

Journal ref: Acta Phys. Pol. B Proc. Suppl. 16, 1-A48 (2023)

The high-density equation of state in heavy-ion collisions: Constraints from proton flow

Authors: Jan Steinheimer, Anton Motornenko, Agnieszka Sorensen, Yasushi Nara, Volker Koch, Marcus Bleicher

Abstract: A set of different equations of state is implemented in the molecular dynamics part of a non-equilibrium transport simulation (UrQMD) of heavy-ion collisions. It is shown how different flow observables are affected by the density dependence of the equation of state. In particular, the effects of a phase transition at high density are explored, including an expected reduction in mean $m_T$. We also… ▽ More A set of different equations of state is implemented in the molecular dynamics part of a non-equilibrium transport simulation (UrQMD) of heavy-ion collisions. It is shown how different flow observables are affected by the density dependence of the equation of state. In particular, the effects of a phase transition at high density are explored, including an expected reduction in mean $m_T$. We also show that an increase in $v_2$ is characteristic for a strong softening of the equation of state. The phase transitions with a low coexistence density, $n_{\mathrm{CE}}<4 n_0$, show a distinct minimum in the slope of the directed flow as a function of the beam energy, which would be a clear experimental signal. By comparing our results with experimental data, we can exclude any strong phase transition at densities below $4n_0$. △ Less

Submitted 14 October, 2022; v1 submitted 25 August, 2022; originally announced August 2022.

Comments: 12 pages, 5 figures, version accepted for publication

Journal ref: Eur. Phys. J. C 82, 911 (2022)

Sensitivity of Au+Au collisions to the symmetric nuclear matter equation of state at 2 -- 5 nuclear saturation densities

Authors: Dmytro Oliinychenko, Agnieszka Sorensen, Volker Koch, Larry McLerran

Abstract: We demonstrate that proton and pion flow measurements in heavy-ion collisions at incident energies ranging from 1 to 20 GeV per nucleon in the fixed target frame can be used for an accurate determination of the symmetric nuclear matter equation of state at baryon densities equal 2--4 times nuclear saturation density $n_0$. We simulate Au+Au collisions at these energies using a hadronic transport m… ▽ More We demonstrate that proton and pion flow measurements in heavy-ion collisions at incident energies ranging from 1 to 20 GeV per nucleon in the fixed target frame can be used for an accurate determination of the symmetric nuclear matter equation of state at baryon densities equal 2--4 times nuclear saturation density $n_0$. We simulate Au+Au collisions at these energies using a hadronic transport model with an adjustable vector mean-field potential dependent on baryon density $n_B$. We show that the mean field can be parametrized to reproduce a given density-dependence of the speed of sound at zero temperature $c_s^2(n_B, T = 0)$, which we vary independently in multiple density intervals to probe the differential sensitivity of heavy-ion observables to the equation of state at these specific densities. Recent flow data from the STAR experiment at the center-of-mass energies $\sqrt{s_{NN}} = \{3.0, 4.5 \}\ $ GeV can be described by our model, and a Bayesian analysis of these data indicates a hard equation of state at $n_B \in (2,3) n_0$ and a possible phase transition at $n_B \in (3,4) n_0$. More data at $\sqrt{s_{NN}} = 2-5$ GeV, as well as a more thorough analysis of the model systematic uncertainties will be necessary for a more precise conclusion. △ Less

Submitted 25 September, 2023; v1 submitted 25 August, 2022; originally announced August 2022.

Comments: 26 pages, 18 figures, accepted to publication in Phys Rev C, minor fixes

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

Authors: Volodymyr Vovchenko, Volker Koch, Chun Shen

Abstract: 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 colli… ▽ More 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. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: 6 pages, 2 figures, contribution to the proceedings of Quark Matter 2022

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

Authors: Volodymyr Vovchenko, Volker Koch

Abstract: 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 me… ▽ More 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. △ Less

Submitted 29 August, 2022; v1 submitted 31 March, 2022; originally announced April 2022.

Comments: 7 pages, 1 figure. v3: typos corrected, version published in Physics Letters B

Journal ref: Phys. Lett. B 833, 137368 (2022)

Critical point particle number fluctuations from molecular dynamics

Authors: Volodymyr A. Kuznietsov, Oleh Savchuk, Mark I. Gorenstein, Volker Koch, Volodymyr Vovchenko

Abstract: 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 se… ▽ More 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. △ Less

Submitted 23 March, 2022; v1 submitted 20 January, 2022; originally announced January 2022.

Comments: 16 pages, 9 figures, to be published in Physical Review C

Cumulants: It's More Than You Think

Authors: Agnieszka Sorensen, Dmytro Oliinychenko, Volker Koch, Larry McLerran

Abstract: Cumulants of baryon number are given considerable attention in analyses of heavy-ion collision experiments as possible signatures of the QCD critical point. In this work, we show that the values of the lowest three cumulants can also be utilized to recover information about the isothermal speed of sound and its logarithmic derivative with respect to the baryon number density. This result provides… ▽ More Cumulants of baryon number are given considerable attention in analyses of heavy-ion collision experiments as possible signatures of the QCD critical point. In this work, we show that the values of the lowest three cumulants can also be utilized to recover information about the isothermal speed of sound and its logarithmic derivative with respect to the baryon number density. This result provides a new method for obtaining information about fundamental properties of nuclear matter studied in heavy-ion collisions, with consequences for both the search for the QCD critical point and neutron star studies. While the approximations and the model comparison we considered apply to experiments at low energies, the approach itself can be used at any collision energy provided that measurements of cumulants of baryon number distribution as well as their temperature dependence are available. △ Less

Submitted 21 October, 2021; originally announced October 2021.

Comments: Submission to Proceedings of CPOD 2021. New discussion on leading contributions to the isothermal and isentropic speeds of sound. Contributed talk on arXiv:2103.07365

The BEST framework for the search for the QCD critical point and the chiral magnetic effect

Authors: Xin An, Marcus Bluhm, Lipei Du, Gerald V. Dunne, Hannah Elfner, Charles Gale, Joaquin Grefa, Ulrich Heinz, Anping Huang, Jamie M. Karthein, Dmitri E. Kharzeev, Volker Koch, Jinfeng Liao, Shiyong Li, Mauricio Martinez, Michael McNelis, Debora Mroczek, Swagato Mukherjee, Marlene Nahrgang, Angel R. Nava Acuna, Jacquelyn Noronha-Hostler, Dmytro Oliinychenko, Paolo Parotto, Israel Portillo, Maneesha Sushama Pradeep , et al. (18 additional authors not shown)

Abstract: The Beam Energy Scan Theory (BEST) Collaboration was formed with the goal of providing a theoretical framework for analyzing data from the Beam Energy Scan (BES) program at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory. The physics goal of the BES program is the search for a conjectured QCD critical point as well as for manifestations of the chiral magnetic effect. W… ▽ More The Beam Energy Scan Theory (BEST) Collaboration was formed with the goal of providing a theoretical framework for analyzing data from the Beam Energy Scan (BES) program at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory. The physics goal of the BES program is the search for a conjectured QCD critical point as well as for manifestations of the chiral magnetic effect. We describe progress that has been made over the previous five years. This includes studies of the equation of state and equilibrium susceptibilities, the development of suitable initial state models, progress in constructing a hydrodynamic framework that includes fluctuations and anomalous transport effects, as well as the development of freezeout prescriptions and hadronic transport models. Finally, we address the challenge of integrating these components into a complete analysis framework. This document describes the collective effort of the BEST Collaboration and its collaborators around the world. △ Less

Submitted 22 November, 2021; v1 submitted 31 August, 2021; originally announced August 2021.

Comments: 103 pages

Journal ref: Nuclear Physics A, Volume 1017, January 2022, 122343

Constraining the hadronic spectrum and repulsive interactions in a hadron resonance gas via fluctuations of conserved charges

Authors: Jamie M. Karthein, Volker Koch, Claudia Ratti, Volodymyr Vovchenko

Abstract: We simultaneously incorporate two common extensions of the hadron resonance gas model, namely the addition of extra, unconfirmed resonances to the particle list and the excluded volume repulsive interactions. We emphasize the complementary nature of these two extensions and identify combinations of conserved charge susceptibilities that allow to constrain them separately. In particular, ratios of… ▽ More We simultaneously incorporate two common extensions of the hadron resonance gas model, namely the addition of extra, unconfirmed resonances to the particle list and the excluded volume repulsive interactions. We emphasize the complementary nature of these two extensions and identify combinations of conserved charge susceptibilities that allow to constrain them separately. In particular, ratios of second-order susceptibilities like $χかい_{11}^{BQ}/χかい_2^B$ and $χかい_{11}^{BS}/χかい_2^B$ are sensitive only to the baryon spectrum, while fourth-to-second order ratios like $χかい_4^B/χかい_2^B$, $χかい_{31}^{BS}/χかい_{11}^{BS}$, or $χかい_{31}^{BQ}/χかい_{11}^{BQ}$ are mainly determined by repulsive interactions. Analysis of the available lattice results suggests the presence of both the extra states in the baryon-strangeness sector and the repulsive baryonic interaction, with indications that hyperons have a smaller repulsive core than non-strange baryons. The modified hadron resonance gas model presented here significantly improves the description of lattice QCD susceptibilities at chemical freeze-out and can be used for the analysis of event-by-event fluctuations in heavy-ion collisions. △ Less

Submitted 8 November, 2021; v1 submitted 1 July, 2021; originally announced July 2021.

Comments: 10 pages, 5 figures, version published in Phys. Rev. D

Journal ref: J. M. Karthein, V. Koch, C. Ratti, and V. Vovchenko, "Constraining the hadronic spectrum and repulsive interactions in a hadron resonance gas via fluctuations of conserved charges," Phys. Rev. D, vol. 104, p. 094009, Nov 2021

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

Authors: Volodymyr Vovchenko, Volker Koch, Chun Shen

Abstract: 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… ▽ More 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. △ Less

Submitted 6 January, 2022; v1 submitted 30 June, 2021; originally announced July 2021.

Comments: 18 pages, 12 figures, published in Physical Review C

Journal ref: Phys. Rev. C 105, 014904 (2022)

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

Authors: Oleh Savchuk, Volodymyr Vovchenko, Volker Koch, Jan Steinheimer, Horst Stoecker

Abstract: 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 reac… ▽ More 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. △ Less

Submitted 20 March, 2022; v1 submitted 15 June, 2021; originally announced June 2021.

Comments: 6 pages, 4 figures

Journal ref: Phys. Lett. B 827, 136983 (2022)

Speed of sound and baryon cumulants in heavy-ion collisions

Authors: Agnieszka Sorensen, Dmytro Oliinychenko, Volker Koch, Larry McLerran

Abstract: We present a method that may allow an estimate of the value of the speed of sound as well as its logarithmic derivative with respect to the baryon number density in matter created in heavy-ion collisions. To this end, we utilize well-known observables: cumulants of the baryon number distribution. In analyses aimed at uncovering the phase diagram of strongly interacting matter, cumulants gather con… ▽ More We present a method that may allow an estimate of the value of the speed of sound as well as its logarithmic derivative with respect to the baryon number density in matter created in heavy-ion collisions. To this end, we utilize well-known observables: cumulants of the baryon number distribution. In analyses aimed at uncovering the phase diagram of strongly interacting matter, cumulants gather considerable attention as their qualitative behaviour along the explored range of collision energies is expected to aid in detecting the QCD critical point. We show that the cumulants may also reveal the behavior of the speed of sound in the temperature and baryon chemical potential plane. We demonstrate the applicability of such estimates within two models of nuclear matter, and explore what might be understood from known experimental data. △ Less

Submitted 26 July, 2021; v1 submitted 12 March, 2021; originally announced March 2021.

Comments: 6 pages, 3 figures; a few minor changes; published version

Journal ref: Phys. Rev. Lett. 127, 042303 (2021)

Exploring the QCD phase diagram with fluctuations

Authors: Volker Koch, Volodymyr Vovchenko

Abstract: We report on recent progress concerning theoretical description of event-by-event fluctuations in heavy-ion collisions. Specifically we discuss a new Cooper-Frye particlization routine -- the subensemble sampler -- which is designed to incorporate effects of global conservation laws, thermal smearing and resonance decays on fluctuation measurements in various rapidity acceptances. First applicatio… ▽ More We report on recent progress concerning theoretical description of event-by-event fluctuations in heavy-ion collisions. Specifically we discuss a new Cooper-Frye particlization routine -- the subensemble sampler -- which is designed to incorporate effects of global conservation laws, thermal smearing and resonance decays on fluctuation measurements in various rapidity acceptances. First applications of the method to heavy-ion collisions at the LHC energies are presented and further necessary steps to analyze fluctuations from RHIC beam energy scan are outlined. △ Less

Submitted 16 February, 2021; originally announced February 2021.

Comments: 13 pages, 3 figures. Prepared for the 60th anniversary of the Cracow School of Theoretical Physics. Contains overview of the material presented in arXiv:2012.09954

Efficiency corrections for factorial moments and cumulants of overlapping sets of particles

Authors: Volodymyr Vovchenko, Volker Koch

Abstract: In this note we discuss subtleties associated with the efficiency corrections for measurements of off-diagonal cumulants and factorial moments for a situation when one deals with overlapping sets of particles, such as correlations between numbers of protons and positively charged particles. In particular, we discuss the situation commonly encountered in heavy-ion experiments, where first all charg… ▽ More In this note we discuss subtleties associated with the efficiency corrections for measurements of off-diagonal cumulants and factorial moments for a situation when one deals with overlapping sets of particles, such as correlations between numbers of protons and positively charged particles. In particular, we discuss the situation commonly encountered in heavy-ion experiments, where first all charges are reconstructed and then protons are selected from these charges by an additional particle identification procedure. We present the efficiency correction formulas for the case when the detection efficiencies follow a binomial distribution. △ Less

Submitted 11 March, 2021; v1 submitted 6 January, 2021; originally announced January 2021.

Comments: 13 pages, 1 figure, to appear in Nucl. Phys. A

Journal ref: Nucl. Phys. A 1010, 122179 (2021)

Particlization of an interacting hadron resonance gas with global conservation laws for event-by-event fluctuations in heavy-ion collisions

Authors: Volodymyr Vovchenko, Volker Koch

Abstract: We revisit the problem of particlization of a QCD fluid into hadrons and resonances at the end of the fluid dynamical stage in relativistic heavy-ion collisions in a context of fluctuation measurements. The existing methods sample an ideal hadron resonance gas, therefore, they do not capture the non-Poissonian nature of the grand-canonical fluctuations, expected due to QCD dynamics such as the chi… ▽ More We revisit the problem of particlization of a QCD fluid into hadrons and resonances at the end of the fluid dynamical stage in relativistic heavy-ion collisions in a context of fluctuation measurements. The existing methods sample an ideal hadron resonance gas, therefore, they do not capture the non-Poissonian nature of the grand-canonical fluctuations, expected due to QCD dynamics such as the chiral transition or QCD critical point. We address the issue by partitioning the particlization hypersurface into locally grand-canonical fireballs populating the space-time rapidity axis that are constrained by global conservation laws. The procedure allows to quantify the effect of global conservation laws, volume fluctuations, thermal smearing and resonance decays on fluctuation measurements in various rapidity acceptances, and can be used in fluid dynamical simulations of heavy-ion collisions. As a first application, we study event-by-event fluctuations in heavy-ion collisions at the LHC using an excluded volume hadron resonance gas model matched to lattice QCD susceptibilities, with a focus on (pseudo)rapidity acceptance dependence of net baryon, net proton, and net charge cumulants. We point out large differences between net proton and net baryon cumulant ratios that make direct comparisons between the two unjustified. We observe that the existing experimental data on net-charge fluctuations at the LHC shows a strong suppression relative to a hadronic description. △ Less

Submitted 17 December, 2020; originally announced December 2020.

Comments: 25 pages, 9 figures

Journal ref: Phys. Rev. C 103, 044903 (2021)

Phase transitions and critical behavior in hadronic transport with a relativistic density functional equation of state

Authors: Agnieszka Sorensen, Volker Koch

Abstract: We develop a flexible, relativistically covariant parameterization of dense nuclear matter equation of state suited for inclusion in computationally demanding hadronic transport simulations. Within an implementation in the hadronic transport code SMASH, we show that effects due to bulk thermodynamic behavior are reproduced in dynamic hadronic systems, demonstrating that hadronic transport can be u… ▽ More We develop a flexible, relativistically covariant parameterization of dense nuclear matter equation of state suited for inclusion in computationally demanding hadronic transport simulations. Within an implementation in the hadronic transport code SMASH, we show that effects due to bulk thermodynamic behavior are reproduced in dynamic hadronic systems, demonstrating that hadronic transport can be used to study critical behavior in dense nuclear matter, both at and away from equilibrium. We also show that two-particle correlations calculated from hadronic transport simulation data follow theoretical expectations based on the second order cumulant ratio, and constitute a clear signature of crossing the phase diagram above the critical point. △ Less

Submitted 3 September, 2021; v1 submitted 12 November, 2020; originally announced November 2020.

Comments: 31 pages, 14 figures. Revised version to be published in PRC (in press). We added a discussion of correlations in hadronic transport and improved the proof of the relationship between the second order cumulant and the pair distribution function at close distances

Journal ref: Phys. Rev. C 104, 034904 (2021)

Cumulants of multiple conserved charges and global conservation laws

Authors: Volodymyr Vovchenko, Roman V. Poberezhnyuk, Volker Koch

Abstract: We analyze the behavior of cumulants of conserved charges in a subvolume of a thermal system with exact global conservation laws by extending a recently developed subensemble acceptance method (SAM) [V. Vovchenko et al., arXiv:2003.13905] to multiple conserved charges. Explicit expressions for all diagonal and off-diagonal cumulants up to sixth order that relate them to the grand canonical suscept… ▽ More We analyze the behavior of cumulants of conserved charges in a subvolume of a thermal system with exact global conservation laws by extending a recently developed subensemble acceptance method (SAM) [V. Vovchenko et al., arXiv:2003.13905] to multiple conserved charges. Explicit expressions for all diagonal and off-diagonal cumulants up to sixth order that relate them to the grand canonical susceptibilities are obtained. The derivation is presented for an arbitrary equation of state with an arbitrary number of different conserved charges. The global conservation effects cancel out in any ratio of two second order cumulants, in any ratio of two third order cumulants, as well as in a ratio of strongly intensive measures $Σしぐま$ and $Δでるた$ involving any two conserved charges, making all these quantities particularly suitable for theory-to-experiment comparisons in heavy-ion collisions. We also show that the same cancellation occurs in correlators of a conserved charge, like the electric charge, with any non-conserved quantity such as net proton or net kaon number. The main results of the SAM are illustrated in the framework of the hadron resonance gas model. We also elucidate how net-proton and net-$Λらむだ$ fluctuations are affected by conservation of electric charge and strangeness in addition to baryon number. △ Less

Submitted 7 September, 2020; v1 submitted 7 July, 2020; originally announced July 2020.

Comments: 39 pages, 5 figures, to be published in JHEP

Journal ref: JHEP 10 (2020) 089

Connecting fluctuation measurements in heavy-ion collisions with the grand-canonical susceptibilities

Authors: Volodymyr Vovchenko, Oleh Savchuk, Roman V. Poberezhnyuk, Mark I. Gorenstein, Volker Koch

Abstract: We derive the relation between cumulants of a conserved charge measured in a subvolume of a thermal system and the corresponding grand-canonical susceptibilities, taking into account exact global conservation of that charge. The derivation is presented for an arbitrary equation of state, with the assumption that the subvolume is sufficiently large to be close to the thermodynamic limit. Our framew… ▽ More We derive the relation between cumulants of a conserved charge measured in a subvolume of a thermal system and the corresponding grand-canonical susceptibilities, taking into account exact global conservation of that charge. The derivation is presented for an arbitrary equation of state, with the assumption that the subvolume is sufficiently large to be close to the thermodynamic limit. Our framework -- the subensemble acceptance method (SAM) -- quantifies the effect of global conservation laws and is an important step toward a direct comparison between cumulants of conserved charges measured in central heavy ion collisions and theoretical calculations of grand-canonical susceptibilities, such as lattice QCD. As an example, we apply our formalism to net-baryon fluctuations at vanishing baryon chemical potentials as encountered in collisions at the LHC and RHIC. △ Less

Submitted 2 November, 2020; v1 submitted 30 March, 2020; originally announced March 2020.

Comments: 7 pages, 2 figures, published version

Journal ref: Phys. Lett. B 811, 135868 (2020)

The QCD phase diagram and statistics friendly distributions

Authors: Volker Koch, Adam Bzdak, Dmytro Oliinychenko, Jan Steinheimer

Abstract: The preliminary STAR data for proton cumulants for central collisions at \sqrt{s}=7.7 GeV are consistent with a two-component proton multiplicity distribution. We show that this two-component distribution is statistics friendly in that factorial cumulants of surprisingly high orders may be extracted with a relatively small number of events. As a consequence the two-component model can be tested an… ▽ More The preliminary STAR data for proton cumulants for central collisions at \sqrt{s}=7.7 GeV are consistent with a two-component proton multiplicity distribution. We show that this two-component distribution is statistics friendly in that factorial cumulants of surprisingly high orders may be extracted with a relatively small number of events. As a consequence the two-component model can be tested and verified right now with the presently available STAR data from the first phase of the RHIC beam energy scan. △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: Proceedings, Quark Matter 2019, Wuhan

Dynamics of critical fluctuations: Theory -- phenomenology -- heavy-ion collisions

Authors: M. Bluhm, M. Nahrgang, A. Kalweit, M. Arslandok, P. Braun-Munzinger, S. Floerchinger, E. S. Fraga, M. Gazdzicki, C. Hartnack, C. Herold, R. Holzmann, Iu. Karpenko, M. Kitazawa, V. Koch, S. Leupold, A. Mazeliauskas, B. Mohanty, A. Ohlson, D. Oliinychenko, J. M. Pawlowski, C. Plumberg, G. W. Ridgway, T. Schäfer, I. Selyuzhenkov, J. Stachel , et al. (5 additional authors not shown)

Abstract: This report summarizes the presentations and discussions during the Rapid Reaction Task Force "Dynamics of critical fluctuations: Theory -- phenomenology -- heavy-ion collisions", which was organized by the ExtreMe Matter Institute EMMI and held at GSI, Darmstadt, Germany in April 2019. We address the current understanding of the dynamics of critical fluctuations in QCD and their measurement in he… ▽ More This report summarizes the presentations and discussions during the Rapid Reaction Task Force "Dynamics of critical fluctuations: Theory -- phenomenology -- heavy-ion collisions", which was organized by the ExtreMe Matter Institute EMMI and held at GSI, Darmstadt, Germany in April 2019. We address the current understanding of the dynamics of critical fluctuations in QCD and their measurement in heavy-ion collision experiments. In addition, we outline what might be learned from studying correlations in other physical systems, such as cold atomic gases. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Comments: 89 pages, 11 figures, report on an ExtreMe Matter Institute EMMI Rapid Reaction Task Force

Effects of local event-by-event conservation laws in ultra-relativistic heavy ion collisions at the particlization

Authors: Dmytro Oliinychenko, Shuzhe Shi, Volker Koch

Abstract: Many simulations of relativistic heavy ion collisions involve the switching from relativistic hydrodynamics to kinetic particle transport. This switching entails the sampling of particles from the distribution of energy, momentum and conserved currents provided by hydrodynamics. Usually this sampling ensures the conservation of these quantities only on the average, i.e. the conserved quantities ma… ▽ More Many simulations of relativistic heavy ion collisions involve the switching from relativistic hydrodynamics to kinetic particle transport. This switching entails the sampling of particles from the distribution of energy, momentum and conserved currents provided by hydrodynamics. Usually this sampling ensures the conservation of these quantities only on the average, i.e. the conserved quantities may actually fluctuate among the sampled particle configurations and only their averages over many such configurations agree with their values from hydrodynamics. Here we apply a recently invented method [Oliinychenko, Koch; PRL 123, 18, 182302 (2019)] to ensure conservation laws for each sampled configuration in spatially compact regions (patches) and study their effects: from the well-known (micro-)canonical suppression of means and variances to little studied (micro-)canonical correlations and higher order fluctuations. Most of these effects are sensitive to the patch size. Many of them do not disappear even in the thermodynamic limit, when the patch size goes to infinity. The developed method is essential for particlization of stochastic hydrodynamics. It is useful for studying the chiral magnetic effect, small systems, and in general for fluctuation and correlation observables. △ Less

Submitted 25 August, 2020; v1 submitted 22 January, 2020; originally announced January 2020.

Comments: 40 pages, 9 figures, version accepted to Phys. Rev. C

Journal ref: Phys. Rev. C 102, 034904 (2020)

Critical point signatures in the cluster expansion in fugacities

Authors: Volodymyr Vovchenko, Carsten Greiner, Volker Koch, Horst Stoecker

Abstract: The QCD baryon number density can formally be expanded into a Laurent series in fugacity, which is a relativistic generalization of Mayer'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 as… ▽ More The QCD baryon number density can formally be expanded into a Laurent series in fugacity, which is a relativistic generalization of Mayer'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. △ Less

Submitted 23 December, 2019; v1 submitted 5 September, 2019; originally announced September 2019.

Comments: 14 pages, 5 figures. Title adjusted, appendix with NJL model analysis added, to be published in Physical Review D

Journal ref: Phys. Rev. D 101, 014015 (2020)

A machine learning study to identify spinodal clumping in high energy nuclear collisions

Authors: Jan Steinheimer, LongGang Pang, Kai Zhou, Volker Koch, Jørgen Randrup, Horst Stoecker

Abstract: The coordinate and momentum space configurations of the net baryon number in heavy ion collisions that undergo spinodal decomposition, due to a first-order phase transition, are investigated using state-of-the-art machine-learning methods. Coordinate space clumping, which appears in the spinodal decomposition, leaves strong characteristic imprints on the spatial net density distribution in nearly… ▽ More The coordinate and momentum space configurations of the net baryon number in heavy ion collisions that undergo spinodal decomposition, due to a first-order phase transition, are investigated using state-of-the-art machine-learning methods. Coordinate space clumping, which appears in the spinodal decomposition, leaves strong characteristic imprints on the spatial net density distribution in nearly every event which can be detected by modern machine learning techniques. On the other hand, the corresponding features in the momentum distributions cannot clearly be detected, by the same machine learning methods, in individual events. Only a small subset of events can be systematically differentiated if only the momentum space information is available. This is due to the strong similarity of the two event classes, with and without spinodal decomposition. In such scenarios, conventional event-averaged observables like the baryon number cumulants signal a spinodal non-equilibrium phase transition. Indeed the third-order cumulant, the skewness, does exhibit a peak at the beam energy ($\mathrm{E}_{\mathrm{lab}}= 3-4$ A GeV), where the transient hot and dense system created in the heavy ion collision reaches the first-order phase transition. △ Less

Submitted 20 December, 2019; v1 submitted 15 June, 2019; originally announced June 2019.

Comments: 24 pages, 14 figures, accepted for publication by JHEP

Journal ref: JHEP12(2019)122

Mapping the Phases of Quantum Chromodynamics with Beam Energy Scan

Authors: Adam Bzdak, ShinIchi Esumi, Volker Koch, Jinfeng Liao, Mikhail Stephanov, Nu Xu

Abstract: We review the present status of the search for a phase transition and critical point as well as anomalous transport phenomena in Quantum Chromodynamics (QCD), with an emphasis on the Beam Energy Scan program at the Relativistic Heavy Ion Collider at Brookhaven National Laboratory. We present the conceptual framework and discuss the observables deemed most sensitive to a phase transition, QCD criti… ▽ More We review the present status of the search for a phase transition and critical point as well as anomalous transport phenomena in Quantum Chromodynamics (QCD), with an emphasis on the Beam Energy Scan program at the Relativistic Heavy Ion Collider at Brookhaven National Laboratory. We present the conceptual framework and discuss the observables deemed most sensitive to a phase transition, QCD critical point, and anomalous transport, focusing on fluctuation and correlation measurements. Selected experimental results for these observables together with those characterizing the global properties of the systems created in heavy ion collisions are presented. We then discuss what can be already learned from the currently available data about the QCD critical point and anomalous transport as well as what additional measurements and theoretical developments are needed in order to discover these phenomena. △ Less

Submitted 16 April, 2020; v1 submitted 3 June, 2019; originally announced June 2019.

Comments: Review article on the physics of beam energy scan in relativistic nuclear collisions, updated version-2 matching the published paper [Physics Reports 853 (2020) pp. 1-87]

Journal ref: Physics Reports 853 (2020) pp. 1-87

All the Fun of the FAIR: Fundamental physics at the Facility for Antiproton and Ion Research

Authors: M. Durante, P. Indelicato, B. Jonson, V. Koch, K. Langanke, Ulf-G. Meißner, E. Nappi, T. Nilsson, Th. Stöhlker, E. Widmann, M. Wiescher

Abstract: The Facility for Antiproton and Ion Research (FAIR) will be the accelerator-based flagship research facility in many basic sciences and their applications in Europe for the coming decades. FAIR will open up unprecedented research opportunities in hadron and nuclear physics, in atomic physics and nuclear astrophysics as well as in applied sciences like materials research, plasma physics and radiati… ▽ More The Facility for Antiproton and Ion Research (FAIR) will be the accelerator-based flagship research facility in many basic sciences and their applications in Europe for the coming decades. FAIR will open up unprecedented research opportunities in hadron and nuclear physics, in atomic physics and nuclear astrophysics as well as in applied sciences like materials research, plasma physics and radiation biophysics with applications towards novel medical treatments and space science. FAIR is currently under construction as an international facility at the campus of the GSI Helmholtzzentrum for Heavy-Ion Research in Darmstadt, Germany. While the full science potential of FAIR can only be harvested once the new suite of accelerators and storage rings is completed and operational, some of the experimental detectors and instrumentation are already available and will be used starting in summer 2018 in a dedicated research program at GSI, exploiting also the significantly upgraded GSI accelerator chain. The current manuscript summarizes how FAIR will advance our knowledge in various research fields ranging from a deeper understanding of the fundamental interactions and symmetries in Nature to a better understanding of the evolution of the Universe and the objects within. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: invited comment, published in Physica Scripta 94 (2019) 033001

Particlization with local event-by-event conservation laws

Authors: Dmytro Oliinychenko, Volker Koch

Abstract: We present a sampling method for the transition from relativistic hydrodynamics to particle transport, commonly referred to as particlization, which preserves the local event-by-event conservation of energy, momentum, baryon number, strangeness, and electric charge. The proposed method is essential for studying fluctuations and correlations by means of stochastic hydrodynamics. It is also useful f… ▽ More We present a sampling method for the transition from relativistic hydrodynamics to particle transport, commonly referred to as particlization, which preserves the local event-by-event conservation of energy, momentum, baryon number, strangeness, and electric charge. The proposed method is essential for studying fluctuations and correlations by means of stochastic hydrodynamics. It is also useful for studying small systems. The method is based on Metropolis sampling applied to particles within distinct patches of the hypersurface, where hydrodynamics and kinetic evolutions are matched. △ Less

Submitted 30 October, 2019; v1 submitted 26 February, 2019; originally announced February 2019.

Comments: 4 pages, 1 figure

Journal ref: Phys. Rev. Lett. 123, 182302 (2019)

Centrality dependence of deuteron production in Pb+Pb collisions at 2.76 TeV via hydrodynamics and hadronic afterburner

Authors: Dmytro Oliinychenko, Long-Gang Pang, Hannah Elfner, Volker Koch

Abstract: The deuteron binding energy is only 2.2 MeV. At the same time, its yield in Pb+Pb collisions at $\sqrt{s_{NN}} = $2.76 TeV corresponds to a thermal yield at the temperature around 155 MeV, which is too hot to keep deuterons bound. This puzzle is not completely resolved yet. In general, the mechanism of light nuclei production in ultra-high energy heavy ion collisions remains under debate. In a pre… ▽ More The deuteron binding energy is only 2.2 MeV. At the same time, its yield in Pb+Pb collisions at $\sqrt{s_{NN}} = $2.76 TeV corresponds to a thermal yield at the temperature around 155 MeV, which is too hot to keep deuterons bound. This puzzle is not completely resolved yet. In general, the mechanism of light nuclei production in ultra-high energy heavy ion collisions remains under debate. In a previous work [1] we suggest a microscopic explanation of the deuteron production in central ultra-relativistic Pb+Pb collisions, the main mechanism being $πぱいpn \leftrightarrow πぱいd$ reactions in the hadronic phase of the collision. We use a state-of-the-art hybrid approach, combining relativistic hydrodynamics for the hot and dense stage and hadronic transport for a later, more dilute stage. Deuteron rescattering in the hadronic stage is implemented explicitly, using its experimentally measured vacuum cross-sections. In these proceedings we extend our previous work to non-central collisions, keeping exactly the same methodology and parameters. We find that our approach leads to a good description of the measured deuteron transverse momentum spectra at centralities up to 40%, and underestimates the amount of deuterons at low transverse momentum at higher centralities. Nevertheless, the coalescence parameter $B_2$, measured by ALICE collaboration, is reproduced well in our approach even for peripheral collisions. △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: Proceedings of Hot Quarks 2018, 4 pages, 2 figures

Mapping the QCD phase diagram with statistics friendly distributions

Authors: Adam Bzdak, Volker Koch

Abstract: We demonstrate that the multiplicity distribution of a system located in the vicinity of a first-order phase transition can be successfully measured in terms of its factorial cumulants with a surprisingly small number of events. This finding has direct implications for the experimental search of a QCD phase transition conjectured to be located in the high baryon density region of the QCD phase dia… ▽ More We demonstrate that the multiplicity distribution of a system located in the vicinity of a first-order phase transition can be successfully measured in terms of its factorial cumulants with a surprisingly small number of events. This finding has direct implications for the experimental search of a QCD phase transition conjectured to be located in the high baryon density region of the QCD phase diagram. △ Less

Submitted 14 November, 2019; v1 submitted 11 November, 2018; originally announced November 2018.

Comments: 6 pages, 3 figures; comments added

Journal ref: Phys. Rev. C 100, 051902 (2019)

Net-baryon multiplicity distribution consistent with lattice QCD

Authors: Adam Bzdak, Volker Koch

Abstract: We determine the net-baryon multiplicity distribution which reproduces all cumulants measured so far by lattice QCD. We present the dependence on the volume and temperature of this distribution. We find that for temperatures and volumes encountered in heavy ion reactions, the multiplicity distribution is very close to the Skellam distribution, making the experimental determination of it rather cha… ▽ More We determine the net-baryon multiplicity distribution which reproduces all cumulants measured so far by lattice QCD. We present the dependence on the volume and temperature of this distribution. We find that for temperatures and volumes encountered in heavy ion reactions, the multiplicity distribution is very close to the Skellam distribution, making the experimental determination of it rather challenging. We further provide estimates for the statistics required to measure cumulants of the net-baryon and net-proton distributions. △ Less

Submitted 19 December, 2018; v1 submitted 3 October, 2018; originally announced October 2018.

Comments: 13 pages, 6 figures; Extended version. Now include statistics estimate for RHIC and LHC based on delta method

Journal ref: Phys. Rev. C 99, 024913 (2019)

Evaluation of particle--anti-particle scaled correlation within effective models

Authors: Andre Garcia, Volker Koch, Marcus Benghi Pinto

Abstract: Correlations and fluctuations of physical quantities are known to play an important role in phase transitions and critical phenomena. In recent years some experimental attempts were made in the scope of the Beam Energy Scan program to locate a possible critical point in the QCD phase diagram. In this work we use the Nambu--Jona-Lasinio model to investigate the off-diagonal quark susceptibility, wh… ▽ More Correlations and fluctuations of physical quantities are known to play an important role in phase transitions and critical phenomena. In recent years some experimental attempts were made in the scope of the Beam Energy Scan program to locate a possible critical point in the QCD phase diagram. In this work we use the Nambu--Jona-Lasinio model to investigate the off-diagonal quark susceptibility, which is related to the quark--anti-quark scaled correlation at the mean field level. We show that this correlation has a significant peak near the critical point and, therefore, may be a useful quantity to measure in experiment. We further study the effects of a repulsive vector coupling, which reduces the strength of the scaled correlation near the critical point. △ Less

Submitted 7 September, 2018; originally announced September 2018.

Comments: 17 pages, 10 figures

Large proton cumulants from the superposition of ordinary multiplicity distributions

Authors: Adam Bzdak, Volker Koch, Dmytro Oliinychenko, Jan Steinheimer

Abstract: We construct a multiplicity distribution characterized by large factorial cumulants (integrated correlation functions) from a simple combination of two ordinary multiplicity distributions characterized by small factorial cumulants. We find that such a model, which could be interpreted as representing two event classes, reproduces the preliminary data for the proton cumulants measured by the STAR c… ▽ More We construct a multiplicity distribution characterized by large factorial cumulants (integrated correlation functions) from a simple combination of two ordinary multiplicity distributions characterized by small factorial cumulants. We find that such a model, which could be interpreted as representing two event classes, reproduces the preliminary data for the proton cumulants measured by the STAR collaboration at $7.7$ GeV very well. This model then predicts very large values for the fifth and sixth order factorial cumulants, which can be tested in experiment. △ Less

Submitted 9 October, 2018; v1 submitted 12 April, 2018; originally announced April 2018.

Comments: 15 pages, 4 figures; comments and addendum (2 pages) added

Journal ref: Phys. Rev. C 98, 054901 (2018)

Femtoscopy of stopped protons

Authors: Andrzej Bialas, Adam Bzdak, Volker Koch

Abstract: The longitudinal proton-proton femtoscopy (HBT) correlation function, based on the idea that in a heavy ion collision at $\sqrt{s} \lesssim 20$ GeV stopped protons are likely to be separated in configuration space, is evaluated. It shows a characteristic oscillation which appears sufficiently pronounced to be accessible in experiment. The proposed measurement is essential for estimating the baryon… ▽ More The longitudinal proton-proton femtoscopy (HBT) correlation function, based on the idea that in a heavy ion collision at $\sqrt{s} \lesssim 20$ GeV stopped protons are likely to be separated in configuration space, is evaluated. It shows a characteristic oscillation which appears sufficiently pronounced to be accessible in experiment. The proposed measurement is essential for estimating the baryon density in the central rapidity region, and can be also viewed as an (almost) direct verification of the Lorentz contraction of the fast-moving nucleus. △ Less

Submitted 26 February, 2019; v1 submitted 26 November, 2017; originally announced November 2017.

Comments: 8 pages, 3 figures; discussion on strong and Coulomb interactions included

Journal ref: Phys. Rev. C 99, 034906 (2019)

Rapidity dependence of proton cumulants and correlation functions

Authors: Adam Bzdak, Volker Koch

Abstract: The dependence of multi-proton correlation functions and cumulants on the acceptance in rapidity and transverse momentum is studied. We find that the preliminary data of various cumulant ratios are consistent, within errors, with rapidity and transverse momentum independent correlation functions. However, rapidity correlations which moderately increase with rapidity separation between protons are… ▽ More The dependence of multi-proton correlation functions and cumulants on the acceptance in rapidity and transverse momentum is studied. We find that the preliminary data of various cumulant ratios are consistent, within errors, with rapidity and transverse momentum independent correlation functions. However, rapidity correlations which moderately increase with rapidity separation between protons are slightly favored. We propose to further explore the rapidity dependence of multi-particle correlation functions by measuring the dependence of the integrated reduced correlation functions as a function of the size of the rapidity window. △ Less

Submitted 10 November, 2017; v1 submitted 9 July, 2017; originally announced July 2017.

Comments: 13 pages, 5 figures

Journal ref: Phys. Rev. C 96, 054905 (2017)

Effect of finite particle number sampling on baryon number fluctuations

Authors: Jan Steinheimer, Volker Koch

Abstract: The effects of finite particle number sampling on the net baryon number cumulants, extracted from fluid dynamical simulations, are studied. The commonly used finite particle number sampling procedure introduces an additional Poissonian (or multinomial if global baryon number conservation is enforced) contribution which increases the extracted moments of the baryon number distribution. If this proc… ▽ More The effects of finite particle number sampling on the net baryon number cumulants, extracted from fluid dynamical simulations, are studied. The commonly used finite particle number sampling procedure introduces an additional Poissonian (or multinomial if global baryon number conservation is enforced) contribution which increases the extracted moments of the baryon number distribution. If this procedure is applied to a fluctuating fluid dynamics framework one severely overestimates the actual cumulants. We show that the sampling of so called test-particles suppresses the additional contribution to the moments by at least one power of the number of test-particles. We demonstrate this method in a numerical fluid dynamics simulation that includes the effects of spinodal decomposition due to a first order phase transition. Furthermore, in the limit where anti-baryons can be ignored, we derive analytic formulas which capture exactly the effect of particle sampling on the baryon number cumulants. These formulas may be used to test the various numerical particle sampling algorithms. △ Less

Submitted 29 September, 2017; v1 submitted 23 May, 2017; originally announced May 2017.

Comments: 9 pages 3 figures

Report number: INT-PUB-17-018

Journal ref: Phys. Rev. C 96, 034907 (2017)

Correlated stopping, proton clusters and higher order proton cumulants

Authors: Adam Bzdak, Volker Koch, Vladimir Skokov

Abstract: We investigate possible effects of correlations between stopped nucleons on higher order proton cumulants at low energy heavy-ion collisions. We find that fluctuations of the number of wounded nucleons $N_{\mathrm{part}}$ lead to rather nontrivial dependence of the correlations on the centrality; however, this effect is too small to explain the large and positive four-proton correlations found in… ▽ More We investigate possible effects of correlations between stopped nucleons on higher order proton cumulants at low energy heavy-ion collisions. We find that fluctuations of the number of wounded nucleons $N_{\mathrm{part}}$ lead to rather nontrivial dependence of the correlations on the centrality; however, this effect is too small to explain the large and positive four-proton correlations found in the preliminary data collected by the STAR collaboration at $\sqrt{s}=7.7$ GeV. We further demonstrate that, by taking into account additional proton clustering, we are able to qualitatively reproduce the preliminary experimental data. We speculate that this clustering may originate either from collective/multi-collision stopping which is expected to be effective at lower energies or from a possible first-order phase transition, or from (attractive) final state interactions. To test these ideas we propose to measure a mixed multi-particle correlation between stopped protons and a produced particle (e.g. pion, antiproton). △ Less

Submitted 15 December, 2016; originally announced December 2016.

Comments: 12 pages, 5 figures

Report number: INT-PUB-16-053, RBRC 1221

Stopped nucleons in configuration space

Authors: Andrzej Bialas, Adam Bzdak, Volker Koch

Abstract: In this note, using the colour string model, we study the configuration space distribution of stopped nucleons in heavy-ion collisions. We find that the stopped nucleons from the target and the projectile end up separated from each other by the distance increasing with the collision energy. In consequence, for the center of mass energies larger than 6 or 10 GeV (depending on the details of the mod… ▽ More In this note, using the colour string model, we study the configuration space distribution of stopped nucleons in heavy-ion collisions. We find that the stopped nucleons from the target and the projectile end up separated from each other by the distance increasing with the collision energy. In consequence, for the center of mass energies larger than 6 or 10 GeV (depending on the details of the model) it appears that the system created is not in thermal and chemical equilibrium, and the net baryon density reached is likely not much higher than that already present in the colliding nuclei. △ Less

Submitted 25 August, 2016; originally announced August 2016.

Comments: 8 pages, 2 figures

Journal ref: Acta Phys. Polon. B49 (2018) 103

Chiral Magnetic Effect Task Force Report

Authors: Vladimir Skokov, Paul Sorensen, Volker Koch, Soeren Schlichting, Jim Thomas, Sergei Voloshin, Gang Wang, Ho-Ung Yee

Abstract: In this report, we briefly examine the current status of the study of the chiral magnetic effect including theory and experimental progress. We recommend future strategies for resolving uncertainties in interpretation including recommendations for theoretical work, recommendations for measurements based on data collected in the past five years, and recommendations for beam use in the coming years… ▽ More In this report, we briefly examine the current status of the study of the chiral magnetic effect including theory and experimental progress. We recommend future strategies for resolving uncertainties in interpretation including recommendations for theoretical work, recommendations for measurements based on data collected in the past five years, and recommendations for beam use in the coming years of RHIC. We have specifically investigated the case for colliding nuclear isobars (nuclei with the same mass but different charge) and find the case compelling. We recommend that a program of nuclear isobar collisions to isolate the chiral magnetic effect from background sources be placed as a high priority item in the strategy for completing the RHIC mission. △ Less

Submitted 12 August, 2016; v1 submitted 2 August, 2016; originally announced August 2016.

Comments: 28 pages, 7 figures; Task Force Report on Chiral Magnetic Effect; v2: references added, minor corrections

Journal ref: Chin.Phys. C41 (2017) no.7, 072001

Cumulants and Correlation Functions vs the QCD phase diagram

Authors: Adam Bzdak, Volker Koch, Nils Strodthoff

Abstract: In this paper we discuss the relation of particle number cumulants and correlation functions. It is argued that measuring couplings of the genuine multi-particle correlation functions could provide cleaner information on possible non-trivial dynamics in heavy-ion collisions. We extract integrated multi-proton correlation functions from the presently available experimental data on proton cumulants.… ▽ More In this paper we discuss the relation of particle number cumulants and correlation functions. It is argued that measuring couplings of the genuine multi-particle correlation functions could provide cleaner information on possible non-trivial dynamics in heavy-ion collisions. We extract integrated multi-proton correlation functions from the presently available experimental data on proton cumulants. We find that the STAR data contain significant four-proton correlations, at least at the lower energies, with indication of changing dynamics in central collisions. We also find that these correlations are rather long-ranged in rapidity. Finally, using the Ising model, we demonstrate how the signs of the multi-proton correlation functions may be used to exclude certain regions of the phase diagram close to the critical point. △ Less

Submitted 5 May, 2017; v1 submitted 25 July, 2016; originally announced July 2016.

Comments: 13 pages, 8 figures; PRC version, improved discussion

Journal ref: Phys. Rev. C 95, 054906 (2017)

Multiplicity dependent and non-binomial efficiency corrections for particle number cumulants

Authors: Adam Bzdak, Romain Holzmann, Volker Koch

Abstract: In this note we extend previous work on efficiency corrections for cumulant measurements [1,2]. We will discuss the limitations of the methods presented in these papers. Specifically we will consider multiplicity dependent efficiencies as well as a non-binomial efficiency distributions. We will discuss the most simple and straightforward methods to implement those corrections. In this note we extend previous work on efficiency corrections for cumulant measurements [1,2]. We will discuss the limitations of the methods presented in these papers. Specifically we will consider multiplicity dependent efficiencies as well as a non-binomial efficiency distributions. We will discuss the most simple and straightforward methods to implement those corrections. △ Less

Submitted 30 March, 2016; originally announced March 2016.

Comments: 11 pages, 7 figures

Journal ref: Phys. Rev. C 94, 064907 (2016)