Source author record

Chun Shen

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

nucl-th hep-ph nucl-ex hep-th hep-ex math.AP

Catalog footprint

Research overview

65works

6topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2026arXiv

Signatures of QCD conductivities in heavy-ion collisions

Dissipative processes are pivotal for understanding the hydrodynamic evolution of hot and dense QCD matter created in relativistic nuclear collisions. The interplay of multiple conserved charges -- net baryon, strangeness, and electric charge -- is of particular interest. We simulate the longitudinal hydrodynamic evolution with the three diffusion currents in a hydrodynamic model with a lattice-QCD-based equation of state, NEOS-4D, and estimate rapidity distributions including diffusive corrections to the phase-space distribution in the presence of multiple charges, which ensure charge conservation at particlization. We determine the response of particle yields at midrapidity to changes in the diagonal and off-diagonal conductivities. Inversely, we find that most components of the conductivity matrix can be constrained experimentally using identified particle multiplicities at different collision energies.

preprint2023arXiv

3D Multi-system Bayesian Calibration with Energy Conservation to Study Rapidity-dependent Dynamics of Nuclear Collisions

Considerable information about the early-stage dynamics of heavy-ion collisions is encoded in the rapidity dependence of measurements. To leverage the large amount of experimental data, we perform a systematic analysis using three-dimensional hydrodynamic simulations of multiple collision systems -- large and small, symmetric and asymmetric. Specifically, we perform fully 3D multi-stage hydrodynamic simulations initialized by a parameterized model for rapidity-dependent energy deposition, which we calibrate on the hadron multiplicity and anisotropic flow coefficients. We utilize Bayesian inference to constrain properties of the early- and late- time dynamics of the system, and highlight the impact of enforcing global energy conservation in our 3D model.

preprint2022arXiv

$Λ$ spin polarization in event-by-event relativistic heavy-ion collisions

We present a systematic study of $Λ$ hyperon's polarization observables using event-by-event (3+1)D relativistic hydrodynamics. The effects of initial hot spot size and QGP's specific shear viscosity on the polarization observables are quantified. We examine the effects of the two formulations of the thermal shear tensor on the polarization observables using the same hydrodynamic background. With event-by-event simulations, we make predictions for the Fourier coefficients of $Λ$ 's longitudinal polarization $P^z$ with respect to the event planes of different orders of anisotropic flow. We propose new correlations among the Fourier coefficients of $P^z$ and charged hadron anisotropic flow coefficients to further test the mapping from fluid velocity gradients to hyperon's polarization. Finally, we present a system size scan with Au+Au, Ru+Ru, and O+O collisions at $\sqrt{s_\mathrm{NN}} = 200$ GeV to study the system size dependence of polarization observables at the Relativistic Heavy-ion Collider.

preprint2022arXiv

Bayesian inference of the fluctuating proton shape

Using Bayesian inference, we determine probabilistic constraints on the parameters describing the fluctuating structure of protons at high energy. We employ the color glass condensate framework supplemented with a model for the spatial structure of the proton, along with experimental data from the ZEUS and H1 Collaborations on coherent and incoherent diffractive $\mathrm{J}/ψ$ production in e+p collisions at HERA. This data is found to constrain most model parameters well. This work sets the stage for future global analyses, including experimental data from e+p, p+p, and p+A collisions, to constrain the fluctuating structure of nucleons along with properties of the final state.

preprint2022arXiv

Bayesian inference of the fluctuating proton shape in DIS and hadronic collisions

We determine the likelihood distribution for the model parameters describing the event-by-event fluctuating proton geometry at small $x$ by performing a Bayesian analysis within the Color Glass Condensate framework. The exclusive $\mathrm{J}/ψ$ production data from HERA is found to constrain the model parameters well, and we demonstrate that complementary constraints can be obtained from simulations of Pb+Pb collisions at the LHC.

preprint2022arXiv

Bulk medium evolution has considerable effects on jet observables

We consider the case, in QCD, of a single jet propagating within a strongly interacting fluid, of finite extent. Interactions lead to the appearance of a source of energy-momentum within the fluid. The remnant jet that escapes the container is analyzed along with portions of the medium excited by the jet. We study the effect of a static versus a semi-realistic expanding medium, with jets traveling inward versus outward. We consider the medium response via recoils in partonic scatterings based on a weakly-coupled description and its combination with hydrodynamical medium response based on a strongly-coupled description, followed by incorporation into a jet. The effects of these limits on the reconstructed energy, momentum and mass of the jet, as a function of the angle away from the original parton direction are studied. It is demonstrated that different flow velocity configurations in the medium produce considerable differences in jet observables. This work highlights the importance of accurate dynamical modeling of the soft medium as a foundation on which to calculate jet modification, and casts skepticism on results obtained without such modeling.

preprint2022arXiv

Longitudinal dynamics and particle production in relativistic nuclear collisions

This work presents a three-dimensional dynamical initialization model for relativistic heavy-ion collisions, implementing local energy-momentum conservation and baryon charge fluctuations at string junctions. Constraining parameters using experimental data from p+p collisions at various collision energies, the model provides a very good description of the charged hadron and net proton rapidity distributions in Au+Au collisions from 7.7 to 200 GeV and Pb+Pb collisions at 8.77 and 17.3 GeV. We demonstrate the importance of fluctuations of baryon densities to string junctions for describing net-proton distributions at collision energies of 62.4 and 200 GeV. Including this improved baryon stopping description along with the requirement of strangeness neutrality also yields a good description of identified particle yields as functions of the collision energy above 7.7 GeV. We further study asymmetric p+Al and (p, d, $^3$ He)+Au collisions at the top RHIC energy and p+Pb, Xe+Xe, and Pb+Pb collisions at LHC energies. We identify the produced particle rapidity distributions in asymmetric collision systems as particularly useful for constraining models of the early-time longitudinal dynamics.

preprint2022arXiv

Multi-messenger heavy-ion collision physics

This work studies the production of direct photons in relativistic nuclear collisions, along with the production of hadrons. Radiation from the very first instants to the final moments of the evolution is included. The hybrid model used here describes all stages of relativistic heavy-ion collisions. Chronologically, those are an initial state reflecting the collision of nuclei described within the Color Glass Condensate effective theory; a pre-equilibrium phase based on non-equilibrium linear response; relativistic viscous hydrodynamics, and a hadronic afterburner. The effect of the pre-equilibrium phase on both photonic and hadronic observables is highlighted for the first time. The potential of photon observables -- spectrum, differential elliptic and triangular flow -- to reveal the chemical equilibration time is studied. Finally, we consider "small collision systems", including proton+nucleus collisions and collisions of light nuclei, as probed by hadronic and electromagnetic observables. We demonstrate how photon production can signal the formation of quark-gluon plasma in such small systems.

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

Unveiling the dynamics of nucleosynthesis in relativistic heavy-ion collisions

Like nucleosynthesis during the early universe, light nuclei are also produced in relativistic heavy-ion collisions. Although the deuteron ( $d$ ) yields in these collisions can be well described by the statistical hadronization model (SHM), which assumes that particle yields are fixed at a common chemical freezeout near the phase boundary between the quark-gluon plasma and the hadron gas, the recently measured triton ( $^3\text{H}$ ) yields in Au+Au collisions at $\sqrt{s_{NN}}=7.7-200$ GeV are overestimated systematically by this model. Here, we develop a comprehensive kinetic approach to study the effects of hadronic re-scatterings, such as $πNN\leftrightarrowπd$ and $πNNN\leftrightarrowπ^3\text{H}~(^3\text{He})$ , on $d$ , $^3\text{H}$ , and $^3\text{He}$ production in these collisions. We find that these reactions have little effects on the deuteron yield but reduce the $^3\text{H}$ and $^3\text{He}$ yields by about a factor of 1.8 from their initial values given by the SHM. This finding helps resolve the overestimation of triton production in the SHM and provides the evidence for hadronic re-scattering effects on nucleosynthesis in relativistic heavy-ion collisions.

preprint2021arXiv

Constraining the nucleon size with relativistic nuclear collisions

The notion of the "size" of nucleons and their constituents plays a pivotal role in the current paradigm of the formation and the fluctuations of the quark-gluon plasma produced in high-energy nuclear collision experiments. We report on state-of-the-art hydrodynamic results showing that the correlation between anisotropic flow, $v_n^2$ , and the mean transverse momentum of hadrons, $[p_t]$ , possesses a unique sensitivity to the nucleon size in off-central heavy-ion collisions. We argue that existing experimental measurements of this observable support a picture where the relevant length scale characterizing the colliding nucleons is of order 0.5 fm or smaller, and we discuss the broad implications of this finding for future global Bayesian analyses aimed at extracting initial-state and medium properties from nucleus-nucleus collision data, including $v_n^2$ - $[p_t]$ correlations. Determinations of the nucleon size in heavy-ion collisions will provide a solid independent constraint on the initial state of small system collisions, and will establish a deep connection between collective flow data in nucleus-nucleus experiments and data on deep inelastic scattering on protons and nuclei.

preprint2021arXiv

Dynamic modeling for heavy-ion collisions

Recent theory progresses in (3+1)D dynamical descriptions of relativistic nuclear collisions at finite baryon density are reviewed. Heavy-ion collisions at different collision energies produce strongly coupled nuclear matter to probe the phase structure of Quantum Chromodynamics (QCD). Dynamical frameworks serve as a quantitative tool to study properties of hot QCD matter and map collisions to the QCD phase diagram. Outstanding challenges are highlighted when confronting theoretical models with the current and forthcoming experimental measurements from the RHIC beam energy scan program.

preprint2020arXiv

A collision geometry-based 3D initial condition for relativistic heavy-ion collisions

We present a simple way to construct 3D initial conditions for relativistic heavy-ion collisions based on the Glauber collision geometry. Local energy and momentum conservation conditions are imposed to set non-trivial constraints on our parameterizations of longitudinal profiles for the system's initial energy density and flow velocity. After calibrating parameters with charged hadron rapidity distributions in central Au+Au collisions, we test model predictions for particle rapidity distributions in d+Au and peripheral Au+Au collisions in the Beam Energy Scan (BES) program at Relativistic Heavy-Ion Collider (RHIC). Simulations and comparisons with measurements are also made for Pb+Pb collisions at Super Proton Synchrotron (SPS) energies. We demonstrate that elliptic flow measurements in heavy-ion collisions at $\sqrt{s} \sim 10$ GeV can set strong constraints on the dependence of Quark-Gluon Plasma shear viscosity on temperature and net baryon chemical potential.

preprint2020arXiv

Bulk properties and multi-particle correlations in large and small systems

Charged particle production is calculated in a hybrid framework consisting of the IP-Glasma initial state, MUSIC viscous relativistic fluid dynamics, and the UrQMD microscopic hadronic cascade. Using one set of parameters, we compute observables for a large variety of collision systems. We compare to experimental data in p+p, p+Pb, Xe+Xe, Au+Au and Pb+Pb collisions at various energies, and make predictions for potential O+O runs at RHIC and LHC.

preprint2020arXiv

Effects of the mean field on fluid dynamics in the relaxation time approximation

In this paper the nonequilibrium correction to the distribution function containing a time and space dependent mass is obtained. Given that, fully consistent fluid dynamic equations are formulated. Then, the physics of the bulk viscosity is elaborated for Boltzmann and Bose-Einstein gases within the relaxation time approximation. It is found that the parametric form of the ratio $ζ/τ_R$ for the quantum gas is affected by the infrared cut-off. This may be an indication that the relaxation time approximation is too crude to obtain a reliable form of bulk viscosity.

preprint2020arXiv

Event-plane decorrelation of photons produced in the early stage of heavy-ion collisions

We study photon production in the early stage of heavy-ion collisions using a multistage model combining IP-Glasma, KøMPøST and relativistic hydrodynamics. We discuss the small momentum anisotropy of these photons, highlighting the role of the photon-hadron event-plane decorrelation. We comment that this singular characteristic of early photons could be used to provide dynamical information on the complex pre-hydrodynamics phase of heavy-ion collisions.

preprint2020arXiv

Exploring the influence of bulk viscosity of QCD on dilepton tomography

The collective behavior of hadrons and of electromagnetic radiation in heavy-ion collisions has been widely used to study the properties of the the Quark-Gluon Plasma (QGP). Indeed this collectivity, as measured by anisotropic flow coefficients, can be used to constrain the transport properties of QGP. The goal of this contribution is to investigate the influence of the specific bulk viscosity ( $ζ/s$ ) on dilepton production, both at Relativistic Heavy-Ion Collider (RHIC) and Large Hadron Collider (LHC) energies. We explore the sensitivity of dileptons to dynamical features that bulk viscosity induces on the evolution of a strongly-interacting medium, and highlight what makes them a valuable probe in the pursuit to also constrain $ζ/s$ .

preprint2020arXiv

Probing Early-Time Dynamics and Quark-Gluon Plasma Transport Properties with Photons and Hadrons

We use a hybrid model which relies on QCD effective kinetic theory -- K\oMP\oST -- to dynamically bridge the gap between IP-Glasma initial states and viscous hydrodynamics, for the theoretical interpretation of heavy-ion collision results. The hydrodynamic phase is then followed by dynamical freeze-out handled by UrQMD. It is observed that the new pre-equilibrium/pre-hydro phase does influence the extraction of transport coefficients. We also show how the early dynamics is reflected in photon observables.

preprint2020arXiv

QCD equation of state at finite densities for nuclear collisions

We construct the QCD equation of state at finite chemical potentials including net baryon, electric charge, and strangeness based on the results of lattice QCD simulations and the hadron resonance gas model. The situation of strangeness neutrality and a fixed charge-to-baryon ratio, which resembles that of heavy nuclei, is considered for the application to relativistic heavy-ion collisions. This increases the values of baryon chemical potential compared to the case of vanishing strangeness and electric charge chemical potentials, modifying the fireball trajectory in the phase diagram. We perform viscous hydrodynamic simulations and demonstrate the importance of multiple chemical potentials for identified particle production in heavy-ion collisions at the RHIC and SPS beam energy scan energies.

preprint2020arXiv

Studying QGP with flow: A theory overview

I review recent developments in the phenomenological study of the quark-gluon plasma (QGP) transport properties based on a personal selection of results that were presented at Quark Matter 2019. The constraints on the temperature dependence of QGP shear and bulk viscosity are summarized. I discuss new theory advancements towards more realistic 3D dynamical simulations of heavy-ion collisions at finite baryon density. The challenges and opportunities of applying hydrodynamics to small collision systems are highlighted.

preprint2020arXiv

Transverse momentum fluctuations and their correlation with elliptic flow in nuclear collision

We propose observables $v_0$ and $v_0(p_T)$ which quantify the relative fluctuations in the total transverse momentum at fixed multiplicity. We first study the factorization of the fixed multiplicity momentum dependent two particle correlation function into a product of $v_0(p_T^a)$ and $v_0(p_T^b)$ within realistic hydrodynamic simulations. Then we present computations of $v_0(p_T)$ for different particle types. We determine the relation between the integrated $v_0$ and previously measured observables, and compare results from a hybrid hydrodynamics based model to experimental data. The effects of bulk viscosity and an initial pre-equilibrium stage on the results are quantified. We find that $v_0$ is strongly correlated with the initial state entropy per elliptic area, $S/A$ . Using this result, we explain how the observed correlations between the elliptic flow and the transverse momentum (both in simulations and experiment) reflect the initial state correlations between $1/A$ and ellipticity $\varepsilon_2$ at fixed multiplicity. We argue that the systematic experimental study of $v_0$ , with the same sophistication as used for the other $v_n$ , can contribute significantly to our understanding of quark gluon plasma properties.

preprint2019arXiv

Hybrid Color Glass Condensate and hydrodynamic description of the Relativistic Heavy Ion Collider small system scan

Multi-particle correlation observables in the Relativistic Heavy Ion Collider small system scan are computed in a framework that contains both initial state momentum anisotropies from the Color Glass Condensate effective theory and final state hydrodynamic evolution. The initial state is computed using the IP-Glasma model and coupled to viscous relativistic hydrodynamic simulations, which are followed by microscopic hadronic transport. All parameters of the calculation were previously constrained using experimental data on Au+Au collisions at the same center of mass energy. We find that the qualitative features of the experimental data, such as the system and centrality dependence of the charged hadron momentum anisotropy, can only be reproduced when final state interactions are present. On the other hand, we also demonstrate that the details of the initial state are crucially important for the quantitative description of observables in the studied small systems, as neglecting the initial transverse flow profile or the initial shear stress tensor, which contain information on the momentum anisotropy from the Color Glass Condensate, has dramatic effects on the produced final state anisotropy. We further show that the initial state momentum anisotropy is correlated with the observed elliptic flow in all small systems, with the effect increasing with decreasing multiplicity. We identify the precise measurement of $v_2$ in d+Au and Au+Au collisions at RHIC energy at the same multiplicity as a means to reveal effects of the initial state momentum anisotropy.

preprint2019arXiv

The energy and scale dependence of $\hat{q}$ and the JET puzzle

We present an attempt to probe the underlying structure of the quark-gluon plasma (QGP) at high resolution, based on the extracted jet transport coefficient $\hat{q}$ . We argue that the exchanged momentum $k$ between the hard parton and the medium varies over a range of scales, and for $k\geq$ 1 GeV, $\hat{q}$ can be expressed in terms of a parton distribution function (PDF). Because the mass of a QGP constituent is unknown, we define a scaling variable $x_N$ to represent the ratio of the parton momentum to the momentum of a self-contained section of the plasma which has a mass of 1 GeV. This scaling variable is used to parametrize the QGP-PDF. Calculations, based on this reconstructed $\hat{q}$ are compared to data sensitive to the hardcore of jets $i.e.,$ the single hadron suppression in terms of the nuclear modification factor $R_{AA}$ and the azimuthal anisotropy parameter $v_{2}$ , as a function of transverse momentum $p_{\mathrm{T}}$ , centrality and energy of the collision. It is demonstrated that the scale evolution of the QGP-PDF is responsible for the reduction in the normalization of $\hat{q}$ between fits to Relativistic Heavy-Ion Collider (RHIC) and Large Hadron Collider (LHC) data; a puzzle, first discovered by the JET collaboration.

preprint2016arXiv

Bulk viscous corrections to photon production in the quark-gluon plasma

Photons radiated in heavy-ion collisions are a penetrating probe, and as such can play an important role in the determination of the quark-gluon plasma (QGP) transport coefficients. In this work we calculate the bulk viscous correction to photon production in two-to-two scattering reactions. Phase-space integrals describing the bulk viscous correction are evaluated explicitly in order to avoid the forward scattering approximation which is shown to be poor for photons at lower energies. We furthermore present hydrodynamical simulations of AA collisions focusing on the effect of this calculation on photonic observables. Bulk corrections are shown to reduce the elliptic flow of photons at higher $p_T$ .

preprint2016arXiv

Direct photon production and jet energy-loss in small systems

Two types of penetrating probes, direct photon and QCD jets, are investigated in the background of a small and rapidly expanding droplet of quark-gluon plasma. The additional thermal electromagnetic radiation results in a $\sim$ 50\% enhancement of the direct photons. In high multiplicity p+Pb collisions, jets can lose a sizeable fraction of their initial energy, leading to a charged hadron $R_\mathrm{pA}$ of $\sim$ 0.8 at a transverse momentum around 10\,GeV. Those two proposed measurements can help understand the apparent collective behaviour observed in small collision systems.

preprint2016arXiv

Electromagnetic radiation and collectivity in small quark-gluon droplets

We study the multiplicity and rapidity dependence of thermal and prompt photon production in p+Pb collisions at 5.02 TeV, using a (3+1)D viscous hydrodynamic framework. Direct photon anisotropic flow coefficients $v^γ_{2,3}$ and nuclear modification factor $R^γ_\mathrm{pPb}(p_T)$ are presented in both the p-going (backward) and the Pb-going (forward) directions. The interplay between initial state cold nuclear effect and final state thermal enhancement at different rapidity regions is discussed. The proposed rapidity dependent thermal photon enhancement and direct photon anisotropic flow observables can elucidate non-trivial longitudinal dynamics of hot quark-gluon plasma droplets created in small collision systems.

preprint2016arXiv

Electromagnetic Radiation from QCD Matter: Theory Overview

Recent theory developments in electromagnetic radiation from relativistic heavy-ion collisions are reviewed. Electromagnetic observables can serve as a thermometer, a viscometer, and tomographic probes to the collision system. The current status of the "direct photon flow puzzle" is highlighted.

preprint2016arXiv

Holographic Photon Production and Anisotropic Flow

The thermal-photon emission from strongly coupled gauge theories at finite temperature via the bottom-up models in holographic QCD in the deconfined phase is studied. The models are constructed to approximately reproduce the electric conductivity obtained from lattice simulations for the quark gluon plasma (QGP). The emission rates are then embedded in hydrodynamic simulations combined with prompt photons and hadronic contributions to analyze the spectra and anisotropic flow of direct photons in RHIC and LHC. In general, the holographic models enhance the yield and improve the agreement in spectra, while they reduce the flow in low $p_T$ and increase it in high $p_T$ .

preprint2016arXiv

Rapidity-dependent jet energy loss in small systems with finite-size effects and running coupling

Longitudinal dynamics of particle production and rapidity-dependent jet energy loss are investigated in small and asymmetric colliding systems. We utilize an improved version of \textsc{martini} in which two improvements are implemented to calculate the effect of the strongly coupled QGP droplet on jet energy loss. We show that those realistic prescriptions improve the results of nuclear modification factor calculations. We also observe visible energy loss of jets in a thermal background of high-multiplicity p-Pb collisions, and a clear correlation between the energy loss and elliptic flow coefficients for energetic particles. We conclude that systematic measurements of jet quenching in central collisions of small systems can support the formation of the QGP droplet.

preprint2016arXiv

The production of photons in relativistic heavy-ion collisions

In this work it is shown that the use of a hydrodynamical model of heavy ion collisions which incorporates recent developments, together with updated photon emission rates greatly improves agreement with both ALICE and PHENIX measurements of direct photons, supporting the idea that thermal photons are the dominant source of direct photon momentum anisotropy. The event-by-event hydrodynamical model uses IP-Glasma initial states and includes, for the first time, both shear and bulk viscosities, along with second order couplings between the two viscosities. The effect of both shear and bulk viscosities on the photon rates is studied, and those transport coefficients are shown to have measurable consequences on the photon momentum anisotropy.

preprint2016arXiv

The Vanishing Pressure Limits of Riemann Solutions to the Chaplygin Euler Equations

The Riemann solutions to Chaplygin Euler equations with a scaled pressure are considered. When the pressure vanishes, there are three cases. The Riemann solution containing two shock waves converges to the delta shock wave solution of the transport equations. During this process, both the strength and propagation speed of the delta shock are investigated in detail. We find that there is something different from that for polytropic or isothermal gas. The Riemann solution containing two rarefaction waves tends to the two contact discontinuity solution to the transport equations as the pressure goes to zero. The intermediate state between the two contact discontinuities is a vacuum state. The Riemann solution containing one rarefaction wave and one shock wave tends to the contact discontinuity solution to transport equations as the pressure vanishes.

preprint2015arXiv

Collision Geometry and Flow in Uranium+Uranium Collisions

Using event-by-event viscous fluid dynamics to evolve fluctuating initial density profiles from the Monte-Carlo Glauber model for U+U collisions, we report a "knee"-like structure in the elliptic flow as a function of collision centrality, located around the 0.5% most central collisions as measured by the final charged multiplicity. This knee is due to the preferential selection of tip-on-tip collision geometries by a high-multiplicity trigger. Such a knee structure is not seen in the STAR data. This rules out the two-component MC-Glauber model for initial energy and entropy production. Hence an enrichment of tip-tip configurations by triggering solely on high-multiplicity in the U+U collisions does not work. On the other hand, by using the Zero Degree Calorimeters (ZDCs) coupled with event-shape engineering such a selection is possible. We identify the selection purity of body-body and tip-tip events in full-overlap U+U collisions. By additionally constraining the asymmetry of the ZDC signals we can further increase the probability of selecting tip-tip events in U+U collisions.

preprint2015arXiv

Collision Geometry and Flow in Uranium+Uranium Collisions

Using event-by-event viscous fluid dynamics to evolve fluctuating initial density profiles from the Monte-Carlo Glauber model for U+U collisions, we report a "knee"-like structure in the elliptic flow as a function of collision centrality, located near 0.5% centrality as measured by the final charged multiplicity. This knee is due to the preferential selection of tip-on-tip collision geometries by a high-multiplicity trigger. Such a knee structure is not seen in the STAR data. This rules out the two-component MC-Glauber model for initial energy and entropy production. An enrichment of tip-tip configurations by triggering solely on high-multiplicity in the U+U collisions thus does not work. On the other hand, using the Zero Degree Calorimeters (ZDCs) coupled with event-shape engineering, we identify the selection purity of body-body and tip-tip events in the full-overlap U+U collisions. With additional constraints on the asymmetry of the ZDC signals one can further increases the probability of selecting tip-tip events in U+U collisions.

preprint2015arXiv

Dilepton emission in high-energy heavy-ion collisions with dissipative hydrodynamics

In this contribution we study the effects of three transport coefficients of dissipative hydrodynamics on thermal dilepton anisotropic flow observables. The first two transport coefficients investigated influence the overall size and growth rate of shear viscous pressure, while the last transport coefficient governs the magnitude of net baryon number diffusion in relativistic dissipative fluid dynamics. All calculations are done using state-of-the-art 3+1D hydrodynamical simulations. We show that thermal dileptons are sensitive probes of the transport coefficients of dissipative hydrodynamics.

preprint2015arXiv

Physics perspectives of heavy-ion collisions at very high energy

Heavy-ion collisions at very high colliding energies are expected to produce a quark-gluon plasma (QGP) at the highest temperature obtainable in a laboratory setting. Experimental studies of these reactions can provide an unprecedented range of information on properties of the QGP at high temperatures. We report theoretical investigations of the physics perspectives of heavy-ion collisions at a future high-energy collider. These include initial parton production, collective expansion of the dense medium, jet quenching, heavy-quark transport, dissociation and regeneration of quarkonia, photon and dilepton production. We illustrate the potential of future experimental studies of the initial particle production and formation of QGP at the highest temperature to provide constraints on properties of strongly interaction matter.

preprint2015arXiv

Pre-equilibrium evolution effects on heavy-ion collision observables

In order to investigate the importance of pre-equilibrium dynamics on relativistic heavy-ion collision observables, we match a highly non-equilibrium early evolution stage, modeled by free-streaming partons generated from the Monte Carlo Kharzeev-Levin-Nardi (MC-KLN) and Monte Carlo Glauber (MC-Glb) models, to a locally approximately thermalized later evolution stage described by viscous hydrodynamics, and study the dependence of final hadronic transverse momentum distributions, in particular their underlying radial and anisotropic flows, on the switching time between these stages. Performing a 3-parameter fit of the measured values for the average transverse momenta $\langle p_\perp \rangle$ for pions, kaons and protons as well as the elliptic and triangular flows of charged hadrons $v_{2,3}^\mathrm{ch}$ , with the switching time $τ_s$ , the specific shear viscosity $η/s$ during the hydrodynamic stage, and the kinetic decoupling temperature $T_\mathrm{dec}$ as free parameters, we find that the preferred "thermalization" times $τ_s$ depend strongly on the model of the initial conditions. MC-KLN initial conditions require an earlier transition to hydrodynamic behavior (at $τ_s \approx$ 0.13 fm/ $c$ ) , followed by hydrodynamic evolution with a larger specific shear viscosity $η/s\approx$ 0.2, than MC-Glb initial conditions which prefer switching at a later time ( $τ_s\approx$ 0.6 fm/ $c$ ) followed by a less viscous hydrodynamic evolution with $η/s\approx$ 0.16. These new results including pre-equilibrium evolution are compared to fits without a pre-equilbrium stage where all dynamic evolution before the onset of hydrodynamic behavior is ignored. In each case, the quality of the dynamical descriptions for the optimized parameter sets, as well as the observables which show the strongest constraining power for the thermalization time, are discussed.

preprint2015arXiv

Probing the dissipative properties of a strongly interacting medium with dileptons

We investigate the effects of the presence of a non-vanishing net baryon number density and its diffusion on dilepton production, within a hydrodynamical description of the medium created at $\sqrt{s_{NN}}=7.7$ GeV collision energy. This energy value is explored within the Beam Energy Scan (BES) program at Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory. Particular attention is devoted to a new dissipative degree of freedom: the net baryon number diffusion ( $V^μ$ ), and to the net baryon number conductivity ( $κ$ ) | a transport coefficients governing the overall magnitude of $V^μ$ . The effects of $κ$ on dilepton production are assessed, with an outlook on how future experimental dilepton data can be used to learn more about $κ$ .

preprint2015arXiv

Recent developments in the theory of electromagnetic probes in relativistic heavy-ion collisions

The theoretical developments in the study of electromagnetic radiation in relativistic heavy-ion collisions are reviewed. The recent progress in the rates for photon and lepton pair production is discussed. Together with the improvements in the hydrodynamic descriptions of the bulk medium, the combined effort is discussed to resolve the "direct photon flow puzzle" in the RHIC and the LHC experiments. Further prediction of the direct photon production in high multiplicity proton-nucleus collisions at the LHC energy can serve as a signature of the quark gluon plasma formation in these small systems. Phenomenological study of dilepton production at finite net baryon density is highlighted at the collision energies available for the RHIC beam energy scan program.

preprint2015arXiv

Shape and flow fluctuations in ultra-central Pb+Pb collisions at the LHC

In ultra-central heavy-ion collisions, anisotropic hydrodynamic flow is generated by density fluctuations in the initial state rather than by geometric overlap effects. For a given centrality class, the initial fluctuation spectrum is sensitive to the method chosen for binning the events into centrality classes. We show that sorting events by total initial entropy or by total final multiplicity yields event classes with equivalent statistical fluctuation properties, in spite of viscous entropy production during the fireball evolution. With this initial entropy-based centrality definition we generate several classes of ultra-central Pb+Pb collisions at LHC energies and evolve the events using viscous hydrodynamics with non-zero shear but vanishing bulk viscosity. Comparing the predicted anisotropic flow coefficients for charged hadrons with CMS data we find that both the Monte Carlo Glauber (MC-Glb) and Monte Carlo Kharzeev-Levin-Nardi (MC-KLN) models produce initial fluctuation spectra that are incompatible with the measured final anisotropic flow power spectrum, for any choice of the specific shear viscosity. In spite of this failure, we show that the hydrodynamic model can qualitatively explain, in terms of event-by-event fluctuations of the anisotropic flow coefficients and flow angles, the breaking of flow factorization for elliptic, triangular and quadrangular flow measured by the CMS experiment. For elliptic flow, this factorization breaking is large in ultra-central collisions. We conclude that the bulk of the experimentally observed flow factorization breaking effects are qualitatively explained by hydrodynamic evolution of initial-state fluctuations, but that their quantitative description requires a better understanding of the initial fluctuation spectrum.

preprint2015arXiv

The iEBE-VISHNU code package for relativistic heavy-ion collisions

The iEBE-VISHNU code package performs event-by-event simulations for relativistic heavy-ion collisions using a hybrid approach based on (2+1)-dimensional viscous hydrodynamics coupled to a hadronic cascade model. We present the detailed model implementation, accompanied by some numerical code tests for the package. iEBE-VISHNU forms the core of a general theoretical framework for model-data comparisons through large scale Monte-Carlo simulations. A numerical interface between the hydrodynamically evolving medium and thermal photon radiation is also discussed. This interface is more generally designed for calculations of all kinds of rare probes that are coupled to the temperature and flow velocity evolution of the bulk medium, such as jet energy loss and heavy quark diffusion.

preprint2015arXiv

The road to precision: Extraction of the specific shear viscosity of the quark-gluon plasma

Recent progress in extracting the specific shear viscosity of the Quark-Gluon Plasma from phenomenological studies is highlighted.

preprint2014arXiv

Anisotropic flow of thermal photons as a quark-gluon plasma viscometer

We present state-of-the-art calculations of viscous photon emission from nuclear collisions at RHIC and LHC. Fluctuating initial density profiles are evolved with event-by-event viscous hydrodynamics. Momentum spectra of thermal photons radiated by these explosively expanding fireballs and their $p_T$ -differential anisotropic flow coefficients $v_n(p_T)$ are computed, both with and without accounting for viscous corrections to the standard thermal emission rates. Viscous corrections to the rates are found to have a larger effect on the $v_n$ coefficients than the viscous suppression of hydrodynamic flow anisotropies. The benefits of taking the ratio of elliptic to triangular flow, $v_2/v_3$ , are discussed, and the spacetime regions which contribute dominantly to the photon flow harmonics are identified. The directed flow $v_1$ of thermal photons is predicted for RHIC and LHC energies.

preprint2014arXiv

Electromagnetic fingerprints of the Little Bang

Measurements of thermal photons emitted from the rapidly expanding hot and dense medium ("Little Bang") formed in ultra relativistic heavy-ion collisions, and their current theoretical interpretation, are reviewed.

preprint2014arXiv

Event-by-event direct photon anisotropic flow in relativistic heavy-ion collisions

We consider directly emitted and hadronic decay photons from event-by-event hydrodynamic simulations. We compute the direct photon anisotropic flow coefficients and compare with recent experimental measurements. We find that it is crucial to include the photon multiplicity as a weighting factor in the definition of $v^γ_n$ . We also investigate the sensitivity of the direct photon spectrum and elliptic flow to the theoretical uncertainty of the photon emission rate in the quark-hadron transition region and to the pre-equilibrium dynamics of relativistic heavy-ion collisions.

preprint2014arXiv

Photon Emission from a Momentum Anisotropic Quark-Gluon Plasma

We compute the photon emission rate from a quark-gluon plasma with an anisotropic particle momentum distribution induced by a non-vanishing local shear pressure tensor. Our calculation includes photon production through Compton scattering and quark-antiquark annihilation at leading order in $α_s$ , with all off-equilibrium corrections to leading order in the momentum anisotropy. For fermions we prove that the Kubo-Martin-Schwinger (KMS) relation holds in the hard loop regime for any particle momentum distribution function that is reflection-symmetric. This supports the equivalence, for 2 to 2 scattering processes, of the diagrammatic and kinetic approaches to calculating the photon emission rate. We compare the viscous rates from these two approaches at weak and realistic coupling strengths and provide parameterizations of the equilibrium and viscous photon emission rates for phenomenological studies in relativistic heavy-ion collisions.

preprint2014arXiv

Thermal photon anisotropic flow serves as a quark-gluon plasma viscometer

Photons are a penetrating probe of the hot and dense medium created in heavy-ion collisions. We present state-of-the-art calculations of viscous photon emission from nuclear collisions at RHIC and LHC. Thermal photons anisotropic flow coefficients v_n are computed, both with and without accounting for viscous corrections to the standard thermal emission rates. These corrections are found to have a larger effect on the v_n coefficients than the viscous suppression of hydrodynamic flow anisotropies. For thermal photons, the ratio v_2{SP}/v_3{SP} shows stronger sensitivity to the quark-gluon plasma (QGP) shear viscosity than for hadrons, and it can thus serve as a sensitive QGP viscometer.

preprint2014arXiv

Thermal photons as a quark-gluon plasma thermometer revisited

Photons are a penetrating probe of the hot medium formed in heavy-ion collisions, but they are emitted from all collision stages. At photon energies below 2-3 GeV, the measured photon spectra are approximately exponential and can be characterized by their inverse logarithmic slope, often called "effective temperature" $T_\mathrm{eff}$ . Modelling the evolution of the radiating medium hydrodynamically, we analyze the factors controlling the value of $T_\mathrm{eff}$ and how it is related to the evolving true temperature $T$ of the fireball. We find that at RHIC and LHC energies most photons are emitted from fireball regions with $T{\,\sim\,}T_\mathrm{c}$ near the quark-hadron phase transition, but that their effective temperature is significantly enhanced by strong radial flow. Although a very hot, high pressure early collision stage is required for generating this radial flow, we demonstrate that the experimentally measured large effective photon temperatures $T_\mathrm{eff}{\,>\,}T_\mathrm{c}$ , taken alone, do not prove that any electromagnetic radiation was actually emitted from regions with true temperatures well above $T_\mathrm{c}$ . We explore tools that can help to provide additional evidence for the relative weight of photon emission from the early quark-gluon and late hadronic phases. We find that the recently measured centrality dependence of the total thermal photon yield requires a larger contribution from late emission than presently encoded in our hydrodynamic model.

preprint2013arXiv

Extracting jet transport coefficient from jet quenching at RHIC and LHC

Within five different approaches to parton propagation and energy loss in dense matter, a phenomenological study of experimental data on suppression of large $p_T$ single inclusive hadrons in heavy-ion collisions at both RHIC and LHC was carried out. The evolution of bulk medium used in the study for parton propagation was given by 2+1D or 3+1D hydrodynamic models which are also constrained by experimental data on bulk hadron spectra. Values for the jet transport parameter $\hat q$ at the center of the most central heavy-ion collisions are extracted or calculated within each model, with parameters for the medium properties that are constrained by experimental data on the hadron suppression factor $R_{AA}$ . For a quark with initial energy of 10 GeV we find that $\hat q\approx 1.2 \pm 0.3$ GeV $^2$ /fm at an initial time $τ_0=0.6$ fm/ $c$ in Au+Au collisions at $\sqrt{s}=200$ GeV/n and $\hat q\approx 1.9 \pm 0.7$ GeV $^2$ /fm in Pb+Pb collisions at $\sqrt{s}=2.76$ TeV/n. Compared to earlier studies, these represent significant convergence on values of the extracted jet transport parameter, reflecting recent advances in theory and the availability of new experiment data from the LHC.

preprint2013arXiv

Fluctuating flow angles and anisotropic flow measurements

Event-by-event fluctuations in the initial density distributions of the fireballs created in relativistic heavy-ion collisions lead to event-by-event fluctuations of the final anisotropic flow angles, and density inhomogeneities in the initial state cause these flow angles to vary with the transverse momentum of the emitted particles. It is shown that these effects lead to characteristically different transverse momentum dependencies for anisotropic flow coefficients extracted from different experimental methods. These differences can be used to experimentally constrain flow angle fluctuations in the final state of heavy-ion collisions which, in turn, are sensitive to the initial state density fluctuations and the shear viscosity of the expanding fireball medium.

preprint2013arXiv

HBT interferometry relative to the triangular flow plane in heavy-ion collisions

The PHENIX Collaboration has reported third-order harmonic oscillations of the source radius parameters when measuring the Hanbury Brown-Twiss correlation function for charged hadrons relative to the triangular flow angle. We explore possible origins of such third-order oscillations with a simple Gaussian source featuring both a triangular geometric deformation and triangular flow. Third-order oscillations of the HBT radii can arise from a purely geometric triangular deformation superimposed on an azimuthally symmetric radial flow, or from a radially symmetry spatial distribution which expands anisotropically with a triangular component in the flow velocity profile. In both cases the final particle momentum distribution features triangular flow. We show that the two alternatives can be distinguished experimentally through the phase of the azimuthal oscillations of the HBT radii relative to the triangular flow plane.

preprint2012arXiv

Collision Energy Dependence of Viscous Hydrodynamic Flow in Relativistic Heavy-Ion Collisions

Using a (2+1)-d viscous hydrodynamical model, we study the dependence of flow observables on the collision energy ranging from sqrt(s)=7.7 A GeV at the Relativistic Heavy Ion Collider (RHIC) to sqrt(s)=2760 A GeV at the Large Hadron Collider (LHC). With a realistic equation of state, Glauber model initial conditions and a small specific shear viscosity eta/s = 0.08, the differential charged hadron elliptic flow v_2^{ch}(p_T,sqrt(s)) is found to exhibit a very broad maximum as a function of sqrt(s) around top RHIC energy, rendering it almost independent of collision energy for 39 < sqrt(s) < 2760 A GeV. Compared to ideal fluid dynamical simulations, this "saturation" of elliptic flow is shifted to higher collision energies by shear viscous effects. For color-glass motivated MC-KLN initial conditions, which require a larger shear viscosity eta/s = 0.2 to reproduce the measured elliptic flow, a similar "saturation" is not observed up to LHC energies, except for very low p_T. We emphasize that this "saturation" of the elliptic flow is not associated with the QCD phase transition, but arises from the interplay between radial and elliptic flow which shifts with sqrt(s) depending on the fluid's viscosity and leads to a subtle cancellation between increasing contributions from light and decreasing contributions from heavy particles to v_2 in the sqrt(s) range where v_2^{ch}(p_T,sqrt(s)) at fixed p_T is maximal. By generalizing the definition of spatial eccentricity epsilon_x to isothermal hyper-surfaces, we calculate epsilon_x on the kinetic freeze-out surface at different collision energies. Up to top RHIC energy, sqrt(s)=200 A GeV, the fireball is still out-of-plane deformed at freeze out, while at LHC energy the final spatial eccentricity is predicted to approach zero.

preprint2012arXiv

Resonance Decay Contributions to Higher-Order Anisotropic Flow Coefficients

We show that in hydrodynamic simulations for relativistic heavy-ion collisions, strong resonance decay calculations can be performed with fewer species of particle resonances while preserving good accuracy in single particle spectra and flow anisotropies. Such partial resonance calculations boost computation efficiency by a factor of 10 which is essential for large scale event-by-event simulations.

preprint2012arXiv

Viscous Flow in Heavy-Ion Collisions from RHIC to LHC

We present a systematic hydrodynamic study of the evolution of hadron spectra and their azimuthal anisotropy from the lowest collision energy studied at the Relativistic Heavy Ion Collider (RHIC), sqrt(s) = 7.7 A GeV, to the highest energy reachable at the Large Hadron Collider (LHC), sqrt(s) = 5500 A GeV. The energy dependence of the flow observables are quantitatively studied for both the Monte-Carlo Glauber and Monte-Carlo Kharzeev-Levin-Nardi (MC-KLN) models. For MC-Glauber model initial conditions with η/s = 0.08, the differential charged hadron elliptic flow v_2^{ch}(p_T, sqrt(s)) is found to exhibit a very broad maximum in the region 39 < sqrt(s) < 2760 A GeV. For MC-KLN initial conditions with η/s = 0.2, a similar "saturation" is not observed up to LHC energies. We emphasize that this "saturation" of elliptic flow arises from the interplay between radial flow and elliptic flow which shifts with sqrt(s) depending on the fluid's viscosity. By generalizing the definition of spatial eccentricity to isothermal hyper-surface, we also calculate ε_x on the kinetic freeze-out surface at different collision energies.

preprint2011arXiv

200 A GeV Au+Au collisions serve a nearly perfect quark-gluon liquid

The specific shear viscosity (eta/s)_QGP of a Quark-Gluon-Plasma (QGP) at temperatures T_c < T < 2T_c is extracted from the centrality dependence of the eccentricity-scaled elliptic flow measured in ultra-relativistic heavy-ion collisions. Coupling viscous fluid dynamics for the QGP with a microscopic transport model for hadronic freeze-out we find that the eccentricity-scaled elliptic flow is a universal function of charged multiplicity per unit overlap area, (1/S)(dN_ch/dy), that depends only on the viscosity but not on the model used for computing the initial fireball eccentricity. Comparing with measurements we find 1 < (4pi)(eta/s)_QGP < 2.5 where the uncertainty range is dominated by model uncertainties for the eccentricity values used to normalize the measured elliptic flow.

preprint2011arXiv

Hadron spectra and elliptic flow for 200 A GeV Au+Au collisions from viscous hydrodynamics coupled to a Boltzmann cascade

It is shown that the recently developed hybrid code VISHNU, which couples a relativistic viscous fluid dynamical description of the quark-gluon plasma (QGP) with a microscopic Boltzmann cascade for the late hadronic rescattering stage, yields an excellent description of charged and identified hadron spectra and elliptic flow measured in 200 A GeV Au+Au collisions at the Relativistic Heavy-Ion Collider (RHIC). Using initial conditions that incorporate event-by-event fluctuations in the initial shape and orientation of the collision fireball and values eta/s for the specific shear viscosity of the quark-gluon plasma that were recently extracted from the measured centrality dependence of the eccentricity-scaled, p_T-integrated charged hadron elliptic flow, we obtain universally good agreement between theory and experiment for the p_T-spectra and differential elliptic flow v_2(p_T) for both pions and protons at all collision centralities.

preprint2011arXiv

Hydrodynamic elliptic and triangular flow in Pb-Pb collisions at sqrt(s)=2.76ATeV

It is shown that a simultaneous comparison of both elliptic and triangular flow from (2+1)-dimensional viscous fluid dynamics with recent measurements in Pb+Pb collisions at the Large Hadron Collider (LHC) favors a small specific shear viscosity (eta/s)_QGP approximately equal to 1/(4pi) for the quark-gluon plasma. Using this viscosity value, the relative magnitude of the elliptic and triangular flow is well described with Monte-Carlo Glauber (MC-Glauber) initial conditions while Monte-Carlo Kharzeev-Levin-Nardi (MC-KLN) initial conditions require twice as large viscosity to reproduce the elliptic flow and then underpredict triangular flow by about 30%.

preprint2011arXiv

Hydrodynamic flow in heavy-ion collisions with large hadronic viscosity

Using the (2+1)-dimensional viscous hydrodynamic code VISH2+1 with a temperature dependent specific shear viscosity (eta/s)(T), we present a detailed study of the influence of a large hadronic shear viscosity and its corresponding relaxation time tau_pi on the transverse momentum spectra and elliptic flow of hadrons produced in 200 A GeV Au+Au collisions. Although theory, in principle, predicts a well-defined relation T*tau_pi = kappa(T)*(eta/s)(T), the precise form of kappa(T) for the matter created in relativistic heavy-ion collisions is not known. For the popular choice kappa=3 the hadron spectra are found to be insensitive to a significant rise of eta/s in the hadronic stage, whereas their differential elliptic flow v_2(p_T) is strongly suppressed by large hadronic viscosity. The large viscous effects on v_2 are strongly reduced if (as theoretically expected) kappa(T) is allowed to grow with decreasing temperature in the hadronic stage. This implies that, until reliable calculations of kappa(T) become available, an extraction of the hadronic shear viscosity from a comparison between VISH2+1 and a microscopic hadron cascade or experimental data requires a simultaneous fit of (eta/s)(T) and kappa(T).

preprint2011arXiv

Radial and elliptic flow in Pb+Pb collisions at the Large Hadron Collider from viscous hydrodynamic

A comprehensive viscous hydrodynamic fit of spectra and elliptic flow for charged hadrons and identified pions and protons from Au+Au collisions of all centralities measured at the Relativistic Heavy Ion Collider is performed and used as the basis for predicting the analogous observables for Pb+Pb collisions at the Large Hadron Collider at sqrt(s)=2.76 and 5.5 A TeV. Comparison with recent measurements of the elliptic flow of charged hadrons by the ALICE experiment shows that the model slightly over-predicts the data if the same (constant) specific shear viscosity eta/s is assumed at both collision energies. In spite of differences in our assumptions for the equation of state, the freeze-out temperature, the chemical composition at freeze-out, and the starting time for the hydrodynamic evolution, our results agree remarkably well with those of Luzum [M. Luzum, Phys. Rev. C 83, 044911 (2011)], indicating robustness of the hydrodynamic model extrapolations. Future measurements of the centrality and transverse momentum dependence of spectra and elliptic flow for identified hadrons predicted here will further test the model and shed light on possible variations of the quark-gluon transport coefficients between RHIC and LHC energies.

preprint2011arXiv

Systematics of parton-medium interaction from RHIC to LHC

Despite a wealth of experimental data for high-P_T processes in heavy-ion collisions, discriminating between different models of hard parton-medium interactions has been difficult. A key reason is that the pQCD parton spectrum at RHIC is falling so steeply that distinguishing even a moderate shift in parton energy from complete parton absorption is essentially impossible. In essence, energy loss models are effectively only probed in the vicinity of zero energy loss and, as a result, at RHIC energies only the pathlength dependence of energy loss offers some discriminating power. At LHC however, this is no longer the case: Due to the much flatter shape of the parton p_T spectra originating from 2.76 AGeV collisions, the available data probe much deeper into the model dynamics. A simultaneous fit of the nuclear suppression at both RHIC and LHC energies thus has great potential for discriminating between various models that yield equally good descriptions of RHIC data alone.

preprint2011arXiv

The QGP shear viscosity -- elusive goal or just around the corner?

With the new viscous hydrodynamic + hadron cascade hybrid code VISHNU, a rather precise (O(25%)) extraction of the QGP shear viscosity (eta/s)_QGP from heavy-ion elliptic flow data is possible if the initial eccentricity of the collision fireball is known with <5% accuracy. At this point, eccentricities from initial state models differ by up to 20%, leading to an O(100%) uncertainty for (eta/s)_QGP. It is shown that a simultaneous comparison of elliptic and triangular flow, v_2 and v_3, puts strong constraints on initial state models and can largely eliminate the present uncertainty in (eta/s)_QGP. The variation of the differential elliptic flow v_2(p_T) for identified hadrons between RHIC and LHC energies provides additional tests of the evolution model.

preprint2011arXiv

The viscosity of quark-gluon plasma at RHIC and the LHC

The specific shear viscosity (eta/s)_QGP of quark-gluon plasma (QGP) can be extracted from elliptic flow data in heavy-ion collisions by comparing them with the dynamical model VISHNU which couples a viscous fluid dynamic description of the QGP with a microscopic kinetic description of the late hadronic rescattering and freeze-out stage. A robust method for fixing (eta/s)_QGP from the collision centrality dependence of the eccentricity-scaled charged hadron elliptic flow is presented. The systematic uncertainties associated with this extraction method are discussed, with specific attention to our presently restricted knowledge of initial conditions. With the (eta/s)_QGP extracted in this way, VISHNU yields an excellent description of all soft-hadron data from Au+Au collisions at top RHIC energy. Extrapolations to Pb+Pb collisions at the LHC, using both a purely hydrodynamic approach and VISHNU, are presented and compared with recent experimental results from the ALICE Collaboration. The LHC data are again well described by VISHNU, with the same (eta/s)_QGP value as at RHIC energies.

preprint2010arXiv

A systematic comparison of jet quenching in different fluid-dynamical models

Comparing four different (ideal and viscous) hydrodynamic models for the evolution of the medium created in 200 AGeV Au-Au collisions, combined with two different models for the path length dependence of parton energy loss, we study the effects of jet quenching on the emission-angle dependence of the nuclear suppression factor R_AA(phi) and the away-side per trigger yield I_AA(phi). Each hydrodynamic model was tuned to provide a reasonable description of the single-particle transverse momentum spectra for all collision centralities, and the energy loss models were adjusted to yield the same pion nuclear suppression factor in central Au-Au collisions. We find that the experimentally measured in-plane vs. out-of-plane spread in R_AA(phi) is better reproduced by models that shift the weight of the parton energy loss to later times along its path. Among the models studied here, this is best achieved by energy loss models that suppress energy loss at early times, combined with hydrodynamic models that delay the dilution of the medium density due to hydrodynamic expansion by viscous heating. We were unable to identify a clear tomographic benefit of a measurement of I_AA(phi) over that of R_AA(phi).

preprint2010arXiv

Systematic parameter study of hadron spectra and elliptic flow from viscous hydrodynamic simulations of Au+Au collisions at sqrt(s_NN) = 200 GeV

Using the (2+1)-dimensional viscous hydrodynamic code VISH2+1, we present systematic studies of the dependence of pion and proton transverse momentum spectra and their elliptic flow in 200 A GeV Au+Au collisions on the parameters of the hydrodynamic model (thermalization time, initial entropy density distribution, decoupling temperature, equation of state, and specific shear viscosity eta/s). We identify a tension between the slope of the proton spectra, which (within hydrodynamic simulations that assume a constant shear viscosity to entropy density ratio) prefer larger eta/s values, and the slope of the p_T-dependence of charged hadron elliptic flow, which prefers smaller values of eta/s. Changing other model parameters does not appear to permit dissolution of this tension

preprint2009arXiv

Incompressibility of asymmetric nuclear matter

The incompressibility $K_sat(δ)$ of isospin asymmetric nuclear matter at its saturation density. Our results show that in the expansion of $K_sat(δ)$ in powers of isospin asymmetry $δ$ , i.e., $K_sat(δ)$ =K_{0}+K_{sat,2}δ^{2}+K_{sat,4}δ^{4}+O(δ^{6})$, the magnitude of the 4th-order K_{sat,4} parameter is generally small. The 2nd-order K_{sat,2} parameter thus essentially characterizes the isospin dependence of the incompressibility of asymmetric nuclear matter at saturation density. Furthermore, the K_{sat,2} can be expressed as K_{sat,2}=K_{sym}-6L-J_{0}/{K_{0}L in terms of the slope parameter $L $and the curvature parameter$ K_{\mathrm{sym}} $of the symmetry energy and the third-order derivative parameter$ J_0 $of the energy of symmetric nuclear matter at saturation density, and we find the higher order$ J_0 $contribution to K_{sat,2} generally cannot be neglected. Also, we have found a linear correlation between K_{sym} and$ L $as well as between$ J_{0}/K_{0} $and$ K_{0} $. Using these correlations together with the empirical constraints on$ K_{0} $and$ L $, the nuclear symmetry energy$ E_sym(ρ_{0})$ at normal nuclear density, and the nucleon effective mass, we have obtained an estimated value of K_{sat,2}=-370 +- 120 MeV for the 2nd-order parameter in the isospin asymmetry expansion of the incompressibility of asymmetric nuclear matter at its saturation density.

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Chun Shen

Research overview

Is this your work?

65 publications

Signatures of QCD conductivities in heavy-ion collisions

3D Multi-system Bayesian Calibration with Energy Conservation to Study Rapidity-dependent Dynamics of Nuclear Collisions

ΛΛΛ spin polarization in event-by-event relativistic heavy-ion collisions

Bayesian inference of the fluctuating proton shape

Bayesian inference of the fluctuating proton shape in DIS and hadronic collisions

Bulk medium evolution has considerable effects on jet observables

Longitudinal dynamics and particle production in relativistic nuclear collisions

Multi-messenger heavy-ion collision physics

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

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

Unveiling the dynamics of nucleosynthesis in relativistic heavy-ion collisions

Constraining the nucleon size with relativistic nuclear collisions

Dynamic modeling for heavy-ion collisions

A collision geometry-based 3D initial condition for relativistic heavy-ion collisions

Bulk properties and multi-particle correlations in large and small systems

Effects of the mean field on fluid dynamics in the relaxation time approximation

Event-plane decorrelation of photons produced in the early stage of heavy-ion collisions

Exploring the influence of bulk viscosity of QCD on dilepton tomography

Probing Early-Time Dynamics and Quark-Gluon Plasma Transport Properties with Photons and Hadrons

QCD equation of state at finite densities for nuclear collisions

Studying QGP with flow: A theory overview

Transverse momentum fluctuations and their correlation with elliptic flow in nuclear collision

Hybrid Color Glass Condensate and hydrodynamic description of the Relativistic Heavy Ion Collider small system scan

The energy and scale dependence of q^\hat{q}q^​ and the JET puzzle

Bulk viscous corrections to photon production in the quark-gluon plasma

Direct photon production and jet energy-loss in small systems

Electromagnetic radiation and collectivity in small quark-gluon droplets

Electromagnetic Radiation from QCD Matter: Theory Overview

Holographic Photon Production and Anisotropic Flow

Rapidity-dependent jet energy loss in small systems with finite-size effects and running coupling

The production of photons in relativistic heavy-ion collisions

The Vanishing Pressure Limits of Riemann Solutions to the Chaplygin Euler Equations

Collision Geometry and Flow in Uranium+Uranium Collisions

Collision Geometry and Flow in Uranium+Uranium Collisions

Dilepton emission in high-energy heavy-ion collisions with dissipative hydrodynamics

Physics perspectives of heavy-ion collisions at very high energy

Pre-equilibrium evolution effects on heavy-ion collision observables

Probing the dissipative properties of a strongly interacting medium with dileptons

Recent developments in the theory of electromagnetic probes in relativistic heavy-ion collisions

Shape and flow fluctuations in ultra-central Pb+Pb collisions at the LHC

The iEBE-VISHNU code package for relativistic heavy-ion collisions

The road to precision: Extraction of the specific shear viscosity of the quark-gluon plasma

Anisotropic flow of thermal photons as a quark-gluon plasma viscometer

Electromagnetic fingerprints of the Little Bang

Event-by-event direct photon anisotropic flow in relativistic heavy-ion collisions

Photon Emission from a Momentum Anisotropic Quark-Gluon Plasma

Thermal photon anisotropic flow serves as a quark-gluon plasma viscometer

Thermal photons as a quark-gluon plasma thermometer revisited

Extracting jet transport coefficient from jet quenching at RHIC and LHC

Fluctuating flow angles and anisotropic flow measurements

HBT interferometry relative to the triangular flow plane in heavy-ion collisions

Collision Energy Dependence of Viscous Hydrodynamic Flow in Relativistic Heavy-Ion Collisions

Resonance Decay Contributions to Higher-Order Anisotropic Flow Coefficients

Viscous Flow in Heavy-Ion Collisions from RHIC to LHC

200 A GeV Au+Au collisions serve a nearly perfect quark-gluon liquid

Hadron spectra and elliptic flow for 200 A GeV Au+Au collisions from viscous hydrodynamics coupled to a Boltzmann cascade

Hydrodynamic elliptic and triangular flow in Pb-Pb collisions at sqrt(s)=2.76ATeV

Hydrodynamic flow in heavy-ion collisions with large hadronic viscosity

Radial and elliptic flow in Pb+Pb collisions at the Large Hadron Collider from viscous hydrodynamic

Systematics of parton-medium interaction from RHIC to LHC

The QGP shear viscosity -- elusive goal or just around the corner?

The viscosity of quark-gluon plasma at RHIC and the LHC

A systematic comparison of jet quenching in different fluid-dynamical models

Systematic parameter study of hadron spectra and elliptic flow from viscous hydrodynamic simulations of Au+Au collisions at sqrt(s_NN) = 200 GeV

Incompressibility of asymmetric nuclear matter

See the researcher in context

Building this map preview

$Λ$ spin polarization in event-by-event relativistic heavy-ion collisions

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

The energy and scale dependence of $\hat{q}$ and the JET puzzle