Source author record

Sangyong Jeon

Sangyong Jeon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

69works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

69 published item(s)

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

Bayesian quantification of the Quark-Gluon Plasma: Improved design and closure demonstration

We present a demonstration of the design sampling and closure for the first comprehensive Bayesian model-to-data comparison of heavy-ion measurements with IP-Glasma initial conditions, in which we combine with state-of-the-art hydrodynamics (MUSIC), particlization (iS3D), and transport (SMASH). We further introduce a systematically-improvable method of sampling design points with better projection properties to explore the parameter space of the model for the first time in heavy-ion collisions.

preprint2022arXiv

Family of New Exact Solutions for Longitudinally Expanding Ideal Fluids

We report on the discovery of new analytical solutions of the equations of relativistic ideal hydrodynamics. In this solution, the fluid expands in the longitudinal direction and contains a plateau structure that extends over a finite range in rapidity and can be either symmetric or asymmetric in that variable. We further calculate the corresponding pseudo-rapidity distribution of hadron yields, and find decent agreement with experimental measurements in high-energy Pb+Pb, Au+Au, p+Pb, and d+Au collisions.

preprint2022arXiv

Transport Model Comparison Studies of Intermediate-Energy Heavy-Ion Collisions

Transport models are the main method to obtain physics information from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. Calculations under controlled conditions of physical input and set-up were performed with various participating codes. These included both calculations of nuclear matter in a box with periodic boundary conditions, and more realistic calculations of heavy-ion collisions. In this intermediate review, we summarize and discuss the present status of the project. We also provide condensed descriptions of the 26 participating codes, which contributed to some part of the project. These include the major codes in use today. We review the main results of the studies completed so far. They show, that in box calculations the differences between the codes can be well understood and a convergence of the results can be reached. These studies also highlight the systematic differences between the two families of transport codes, known as BUU and QMD type codes. However, when the codes were compared in full heavy-ion collisions using different physical models, as recently for pion production, they still yielded substantially different results. This calls for further comparisons of heavy-ion collisions with controlled models and of box comparisons of important ingredients, like momentum-dependent fields, which are currently underway. We often indicate improved strategies in performing transport simulations and thus provide guidance to code developers. Results of transport simulations of heavy-ion collisions from a given code will have more significance if the code can be validated against benchmark calculations such as the ones summarized in this review.

preprint2021arXiv

Hadronic transport coefficients from the linear sigma model at finite temperature

We investigate general frameworks for calculating transport coefficients for quasiparticle theories at finite temperature. Hadronic transport coefficients are then computed using the linear sigma model (LSM). The bulk viscosity over entropy density ($ζ/s$) is evaluated in the relaxation time approximation (RTA) and the specific shear viscosity ($η/s$) and static electrical conductivity ($σ_{el}/T$) are both obtained in the RTA and using a functional variational approach. Results are shown for different values of the scalar-isoscalar hadron vacuum mass with in-medium masses for the interacting fields. The advantages and limitations of the LSM for studies of strongly interacting matter out of equilibrium are discussed and results are compared with others in the literature.

preprint2021arXiv

The momentum broadening of energetic partons in an anisotropic plasma

The quark-gluon plasma produced in heavy-ion collisions is anisotropic throughout its evolution. This anisotropy changes the physics of jet-medium interaction, making it dependent on the momentum direction of the jet. In this paper we analyze transverse momentum broadening of a jet parton interacting with soft gluons in an anisotropic plasma. Our analysis equally applies to momentum broadening of quasiparticles in kinetic theory. We subtract contribution from instability modes in the deep infrared and discuss how our calculation should be complemented in that regime. The resulting anisotropic collision kernel for momentum broadening is qualitatively different from the equilibrium collision kernel and from the isotropic ansatz used in effective kinetic theory. Because of increased medium screening, there is substantially less transverse broadening at low and intermediate momenta.

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

Electromagnetic radiation from the pre-equilibrium/pre-hydro stage of the quark-gluon plasma

Photons produced in the pre-equilibrium/pre-hydro stage of the quark-gluon plasma produced in relativistic heavy-ion collisions were computed using parton distribution functions obtained from solutions of the Boltzmann equation. The effect of the initial gluon momentum anisotropy $ξ$ and the dependence on the saturation momentum $Q_s$ was investigated. We see that small $Q_s$ results in a photon yield enhancement, whereas a larger $Q_s$ results in a pre-equilibrium photon suppression, owing to the strict constraint of matching to experimental energy density.

preprint2020arXiv

Exploring Longitudinal Observables with 3+1D IP-Glasma

We present a formulation of the initial state of heavy ion collisions that generalizes the 2+1D boost invariant IP-Glasma \cite{Schenke:2012wb} to 3+1D through JIMWLK rapidity evolution of the pre-collision Wilson lines. The rapidity dependence introduced by the JIMWLK evolution leads us to modify the initial condition for the gauge fields, and to solve Gauss' law iteratively in order to allow for temporal evolution on a 3-dimensional lattice. While the transverse physics of QGP has been studied nearly exhaustively, the effect of longitudinal fluctuations introduced by the JIMWLK evolution has yet to be studied in detail phenomenologically. Hence, we couple our 3+1D IP-Glasma model to MUSIC+UrQMD, for completely 3+1D simulations of heavy ion collisions. Specifically, we consider Pb-Pb collisions at $\sqrt{s} = 2.76\, {\rm TeV}$ and study the rapidity dependence of the charged hadron $v_n(η)$ via the $η$-dependent flow factorization ratios $r_n(η_a,η_b)$ as measured by CMS \cite{Khachatryan:2015oea}, as well as the charged hadron multiplicity $dN_{ch}/dη$.

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

From Chiral Kinetic Theory To Spin Hydrodynamics

The total angular momentum is conserved in the evolution of the Quark-Gluon Plasma (QGP) created in heavy-ion collision, and consists of two sectors: the orbital angular momentum (OAM) caused by kinetic motion, and the spin, an intrinsic degree of freedom of quarks and gluons. Microscopic scattering processes allow the conversion between these two components, hence the spin density could eventually have non-trivial influence on the QGP evolution. A hydrodynamic theory, with the spin polarization effect properly taken into account, is required to quantitatively study the polarization rate for observed hadrons, e.g. the $Λ$-hyperon. In this work, we start with chiral kinetic theory and construct the spin hydrodynamic framework for a chiral spinor system. We obtain the equations of motion of second-order dissipative relativistic fluid dynamics with non-trivial spin polarization density.

preprint2020arXiv

Hard probes of non-equilibrium quark-gluon plasma

Jets and photons could play an important role in finding the transport coefficients of the quark-gluon plasma. To this end we analyze their interaction with a non-equilibrium quark-gluon plasma. Using new field-theoretical tools we derive two-point correlators for the plasma which show how instabilities evolve in time. This allows us, for the first time, to derive finite rates of interaction with the medium. We furthermore show that coherent, long-wavelength instability fields in the Abelian limit do not modify the rate of photon emission or jet-medium interaction.

preprint2020arXiv

The emission of electromagnetic radiation from the early stages of relativistic heavy-ion collisions

We estimate the production of electromagnetic radiation (real and virtual photons) from the early, pre-equilibrium, stage of relativistic heavy-ion collisions. The parton dynamics are obtained as a solution of the Boltzmann equation in the Fokker-Planck diffusion limit. The photon and dilepton rates are integrated and the obtained yields are compared with those from standard sources and with available experimental data. Non-equilibrium photon spectra are predicted for Pb+Pb at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

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 as a probe of the initial state and of viscous dynamics in relativistic nuclear collisions

The penetrating nature of electromagnetic signals makes them suitable probes to explore the properties of the strongly-interacting medium created in relativistic nuclear collisions. We examine the effects of the initial conditions and shear relaxation time on the spectra and flow coefficients of electromagnetic probes, using an event-by-event 3+1D viscous hydrodynamic simulation (MUSIC).

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.

preprint2015arXiv

Introduction to Hydrodynamics

We give a pedagogical review of relativistic hydrodynamics relevant to relativistic heavy ion collisions. Topics discussed include linear response theory derivation of 2nd order viscous hydrodynamics including the Kubo formulas, kinetic theory derivation of 2nd order viscous hydrodynamics, anisotropic hydrodynamics and a brief review of numerical algorithms. Emphasis is given to the theory of hydrodynamics rather than phenomenology.

preprint2015arXiv

Kinetic theory of a longitudinally expanding system of scalar particles

A simple kinematical argument suggests that the classical approximation may be inadequate to describe the evolution of a system with an anisotropic particle distribution. In order to verify this quantitatively, we study the Boltzmann equation for a longitudinally expanding system of scalar particles interacting with a $ϕ^4$ coupling, that mimics the kinematics of a heavy ion collision at very high energy. We consider only elastic $2\to 2$ scatterings, and we allow the formation of a Bose-Einstein condensate in overpopulated situations by solving the coupled equations for the particle distribution and the particle density in the zero mode. For generic CGC-like initial conditions with a large occupation number and a moderate coupling, the solutions of the full Boltzmann equation do not follow a classical attractor behavior.

preprint2015arXiv

Perturbative Vacuum Wavefunctional for Gauge Theories in the Milne Space

The spectrum of vacuum fluctuations in the Milne space (i.e. the tau-eta coordinate system) is an important ingredient in the thermalization studies in relativistic heavy ion collisions. In this paper, the Schrodinger functional for the gauge theory perturbative vacuum is derived for the Milne space. The Wigner-transform of the corresponding vacuum density functional is also found together with the propagators. We finally identify the fluctuation spectrum in vacuum, and show the equivalence between the present approach and the symplectic product based method.

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

The domain of validity of fluid dynamics and the onset of cavitation in ultrarelativistic heavy ion collisions

If the bulk viscosity of QCD matter is large, the effective pressure of the hot and dense matter created in ultrarelativistic heavy ion collisions can become negative, leading to instabilities in the evolution of the plasma. In the context of heavy ion collisions, this effect is sometimes referred to as cavitation. In this contribution we discuss the onset of cavitation in event-by-event hydrodynamic simulations of ultrarelativistic heavy ion collisions at LHC energies. We estimate how large the bulk viscosity of the QGP has to be in the QCD (pseudo) phase transition region in order for the effective pressure of the system to actually become negative.

preprint2014arXiv

Probing the early-time dynamics of relativistic heavy-ion collisions with electromagnetic radiation

Using 3+1D viscous relativistic fluid dynamics, we show that electromagnetic probes are sensitive to the initial conditions and to the out-of-equilibrium features of relativistic heavy-ion collisions. Within the same approach, we find that hadronic observables show a much lesser sensitivity to these aspects. We conclude that electromagnetic observables allow access to dynamical regions that are beyond the reach of soft hadronic probes.

preprint2014arXiv

Probing the non-equilibrium dynamics of hot and dense QCD with dileptons

It is argued that, in heavy ion collisions, thermal dileptons are good probes of the transport properties of the medium created in such events, and also of its early-time dynamics, usually inaccessible to hadronic observables. In this work we show that electromagnetic azimuthal momentum anisotropies do not only display a sensitivity to the shear relaxation time and to the initial shear-stress tensor profile, but also to the temperature dependence of the shear viscosity coefficient.

preprint2014arXiv

Signatures of collective behavior in small systems

We perform 3+1D viscous hydrodynamics calculations of proton-nucleus (pA) and nucleus-nucleus (AA) collisions. Our goal is to understand the apparent collective behavior recently observed in pA collisions and to verify whether the highest multiplicity collision systems can be accurately described as a relativistic fluid. We compare our calculations of flow variables to existing measurements, and demonstrate that hydrodynamics correctly captures the measured trends. We show that our predictions for the pair correlation observable $r_n$ are validated by recent experimental pA measurements, and that our results are sensitive to the granularity of the initial state. We also compare our results with measurements done for nucleus-nucleus collisions.

preprint2014arXiv

Transverse momentum structure of pair correlations as a signature of collective behavior in small collision systems

We perform 3+1D viscous hydrodynamic calculations of proton-lead and lead-lead collisions at top LHC energy. We show that existing data from high-multiplicity p-Pb events can be well described in hydrodynamics, suggesting that collective flow is plausible as a correct description of these collisions. However, a more stringent test of the presence of hydrodynamic behavior can be made by studying the detailed momentum dependence of two-particle correlations. We define a relevant observable, $r_n$, and make predictions for its value and centrality dependence if hydrodynamics is a valid description. This will provide a non-trivial confirmation of the nature of the correlations seen in small collision systems, and potentially to determine where the hydrodynamic description, if valid anywhere, stops being valid. Lastly, we probe what can be learned from this observable, finding that it is insensitive to viscosity, but sensitive to aspects of the initial state of the system that other observables are insensitive to, such as the transverse length scale of the fluctuations in the initial stages of the collision.

preprint2013arXiv

Color Glass Condensate in Schwinger-Keldysh QCD

Within the Schwinger-Keldysh representation of many-body QCD, it is shown that the leading quantum corrections to the strong classical color field is classical in the sense that the fluctuation field still obeys the classical Jacobi-field equation, while the quantum effects solely resides in the fluctuations of the initial field configurations. Within this context, a systematic derivation of the JIMWLK renormalization group equation is presented. A clear identification of the correct form of gauge transformation rules and the correct form of the matter-field Lagrangian in the Schwinger-Keldysh QCD is also presented.

preprint2013arXiv

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

The invariant mass spectrum and the elliptic flow of lepton pairs produced in relativistic heavy-ion collisions at RHIC are studied with viscous hydrodynamics. The effects of viscous corrections on dilepton observables are explored. The lepton pairs originating from charm quarks evolving in the viscous background are seen to be a good probe of quark energy loss and gain, as quantified by the dilepton spectrum and by the dilepton elliptic flow.

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

Hydrodynamic Modeling of Heavy-Ion Collisions

We review progress in the hydrodynamic description of heavy-ion collisions, focusing on recent developments in modeling the fluctuating initial state and event-by-event viscous hydrodynamic simulations. We discuss how hydrodynamics can be used to extract information on fundamental properties of quantum-chromo-dynamics from experimental data, and review successes and challenges of the hydrodynamic framework.

preprint2013arXiv

On the fluid behavior of a baryon rich hadron resonance gas

We investigate the effects of finite baryon chemical potential on the transport properties of a hadron resonance gas. We find that a hadron resonance gas with large baryon number density is closer to the ideal fluid limit than the corresponding gas with zero baryon number. This suggests that the system created at the Relativistic Heavy Ion Collider (RHIC) at lower collision energies may behave as a fluid, with an effective fluidity close to the one found at RHIC's highest energy near phase transition. This might explain why the differential elliptic flow coefficient measured at lower collisional energies at RHIC is similar to the one observed at high energies.

preprint2013arXiv

Solutions of Conformal Israel-Stewart Relativistic Viscous Fluid Dynamics

We use symmetry arguments developed by Gubser to construct the first radially-expanding explicit solutions of the Israel-Stewart formulation of hydrodynamics. Along with a general semi-analytical solution, an exact analytical solution is given which is valid in the cold plasma limit where viscous effects from shear viscosity and the relaxation time coefficient are important. The radially expanding solutions presented in this paper can be used as nontrivial checks of numerical algorithms employed in hydrodynamic simulations of the quark-gluon plasma formed in ultra-relativistic heavy ion collisions. We show this explicitly by comparing such analytic and semi-analytic solutions with the corresponding numerical solutions obtained using the MUSIC viscous hydrodynamics simulation code.

preprint2012arXiv

Comparison of Jet Quenching Formalisms for a Quark-Gluon Plasma "Brick"

We review the currently available formalisms for radiative energy loss of a high-momentum parton in a dense strongly interacting medium. The underlying theoretical framework of the four commonly used formalisms is discussed and the differences and commonalities between the formalisms are highlighted. A quantitative comparison of the single gluon emission spectra as well as the energy loss distributions is given for a model system consisting of a uniform medium with a fixed length of L=2 fm and L=5 fm (the `Brick'). Sizable quantitative differences are found. The largest differences can be attributed to specific approximations that are made in the calculation of the radiation spectrum.

preprint2012arXiv

Event-by-event anisotropic flow in heavy-ion collisions from combined Yang-Mills and viscous fluid dynamics

Anisotropic flow coefficients v_1-v_5 in heavy ion collisions are computed by combining a classical Yang-Mills description of the early time glasma flow with the subsequent relativistic viscous hydrodynamic evolution of matter through the quark-gluon plasma and hadron gas phases. The glasma dynamics, as realized in the IP-Glasma model, takes into account event-by-event geometric fluctuations in nucleon positions and intrinsic sub-nucleon scale color charge fluctuations; the pre-equilibrium flow of matter is then matched to the MUSIC algorithm describing viscous hydrodynamic flow and particle production at freeze-out. The IP-Glasma+MUSIC model describes well both transverse momentum dependent and integrated v_n data measured at the Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC). The model also reproduces the event-by-event distributions of v_2, v_3 and v_4 measured by the ATLAS collaboration. The implications of our results for better understanding of the dynamics of the glasma as well as for the extraction of transport properties of the quark-gluon plasma are outlined.

preprint2012arXiv

Initial state fluctuations and higher harmonic flow in heavy-ion collisions

A framework combining Yang-Mills dynamics of the pre-equilibrium glasma with relativistic viscous hydrodynamic evolution of the quark-gluon plasma and hadron gas phases is presented. Event-by-event fluctuations of nucleon positions and color charges are taken into account, leading to negative binomial fluctuations of gluon multiplicities. Experimental anisotropic flow coefficients v2-v5 of charged hadron distributions in heavy-ion collisions at the Large Hadron Collider are well described. Furthermore, event-by-event distributions of v2, v3 and v4 measured by the ATLAS collaboration are reproduced.

preprint2012arXiv

MUSIC with the UrQMD Afterburner

As RHIC is entering the precision measurement era and the LHC is producing a copious amount of new data, the role of 3+1D event-by-event viscous hydrodynamics is more important than ever to understand the bulk data as well as providing the background for hard probes. For more meaningful comparison with the experimental data, it is also important that hydrodynamics be coupled to the hadronic afterburner. In this proceeding we report on preliminary results of coupling MUSIC with UrQMD.

preprint2012arXiv

Photon Production from a Quark-Gluon-Plasma at Finite Baryon Chemical Potential

We compute the photon production of a QCD plasma at leading order in the strong coupling with a finite baryon chemical potential. Our approach starts from the real time formalism of finite temperature field theory. We identify the class of diagrams contributing at leading order when a finite chemical potential is added and resum them to perform a full treatment of the LPM effect similar to the one performed by Arnold, Moore, and Yaffe at zero chemical potential. Our results show that the contribution of $2\mapsto 3$ and $3\mapsto 2$ processes grows as the chemical potential grows.

preprint2012arXiv

Realistic modelling of jets in heavy-ion collisions

The reconstruction of jets in heavy-ion collisions provides insight into the dynamics of hard partons in media. Unlike the spectrum of single hadrons, the spectrum of jets is highly sensitive to $\hat{q}_{\perp}$, as well as being sensitive to partonic energy loss and radiative processes. We use MARTINI, an event generator, to study how finite-temperature processes at leading order affect dijets.

preprint2011arXiv

A holographic model for large N thermal QCD

We summarize the dual gravity description for a thermal gauge theory, reviewing the key features of our holographic model of large N QCD and elaborating on some new results.The theory has matter in the fundamental representation and the gauge coupling runs logarithmically with energy scale at low energies. At the highest energies the theory becomes approximately scale invariant, much like what we would expect for large N QCD although not with asymptotic freedom. In this limit the theory has a gravity dual captured by an almost classical supergravity description with a controlled quantum behavior, such that by renormalizing the supergravity action, we can compute the stress tensor of the dual gauge theory. From the stress tensor we obtain the shear viscosity and the entropy of the medium at a temperature T, and the violation of the bound for the viscosity to the entropy ratio is then investigated. By considering dynamics of open strings in curved spacetime described by the supergravity limit, we compute the drag and diffusion coefficients for a heavy parton traversing the thermal medium. It is shown that both coefficients have a logarithmic dependence on momentum, consistent with pQCD expectations. Finally, we study the confinement/deconfinement mechanism for quarks by analyzing open strings in the presence of the flavor seven branes. We find linear confinement of quarks at low temperatures, while at high temperatures the quarkonium states melt, a behavior consistent with the existence of a deconfined phase.

preprint2011arXiv

A shear spectral sum rule in a non-conformal gravity dual

A sum rule which relates a stress-energy tensor correlator to thermodynamic functions is examined within the context of a simple non-conformal gravity dual. Such a sum rule was previously derived using AdS/CFT for conformal $\mathcal{N} = 4$ Supersymmetric Yang-Mills theory, but we show that it does not generalize to the non-conformal theory under consideration. We provide a generalized sum rule and numerically verify its validity. A useful byproduct of the calculation is the computation of the spectral density in a strongly coupled non-conformal theory. Qualitative features of the spectral densities and implications for lattice measurements of transport coefficients are discussed.

preprint2011arXiv

Anisotropic flow in sqrt(s)=2.76 TeV Pb+Pb collisions at the LHC

The results on elliptic flow in sqrt(s)=2.76 TeV Pb+Pb collisions at the Large Hadron Collider (LHC) reported by the ALICE collaboration are remarkably similar to those for sqrt(s)=200 GeV gold-gold collisions at the Relativistic Heavy Ion Collider (RHIC). This result is surprising, given the expected longer lifetime of the system at the higher collision energies. We show that it is nevertheless consistent with 3+1 dimensional viscous event-by-event hydrodynamic calculations, and demonstrate that elliptic flow at both RHIC and LHC is built up mostly within the first 5 fm/c of the evolution. We conclude that an "almost perfect liquid" is produced in heavy-ion collisions at the LHC. Furthermore, we present predictions for triangular flow as a function of transverse momentum for different centralities.

preprint2011arXiv

Dijet asymmetry at the Large Hadron Collider

The MARTINI numerical simulation allows for direct comparison of theoretical model calculations and the latest results for dijet asymmetry from the ATLAS and CMS collaborations. In this paper, partons are simulated as undergoing radiative and collisional processes throughout the evolution of central lead-lead collisions at the Large Hadron Collider. Using hydrodynamical background evolution determined by a simulation which fits well with the data on charged particle multiplicities from ALICE and a value of $α_s\approx 0.25-0.3$, the dijet asymmetry is found to be consistent with partonic energy loss in a hot, strongly-interacting medium.

preprint2011arXiv

Higher flow harmonics from (3+1)D event-by-event viscous hydrodynamics

We present event-by-event viscous hydrodynamic calculations of the anisotropic flow coefficients v_2 to v_5 for heavy-ion collisions at the Relativistic Heavy-Ion Collider (RHIC). We study the dependence of different flow harmonics on shear viscosity and the morphology of the initial state. v_3 and higher flow harmonics exhibit a particularly strong dependence on both the initial granularity and shear viscosity. We argue that a combined analysis of all available flow harmonics will allow to determine eta/s of the quark gluon plasma precisely. Presented results strongly hint at a value (eta/s)_QGP < 2/4pi at RHIC. Furthermore, we demonstrate the effect of shear viscosity on pseudo-rapidity spectra and the mean transverse momentum as a function of rapidity.

preprint2011arXiv

Hydrodynamic evolution and jet energy loss in Cu+Cu collisions

We present results from a hybrid description of Cu+Cu collisions using 3+1 dimensional hydrodynamics (MUSIC) for the bulk evolution and a Monte-Carlo simulation (MARTINI) for the evolution of high momentum partons in the hydrodynamical background. We explore the limits of this description by going to small system sizes and determine the dependence on different fractions of wounded nucleon and binary collisions scaling of the initial energy density. We find that Cu+Cu collisions are well described by the hybrid description at least up to 20% central collisions.

preprint2011arXiv

MARTINI for heavy quarks: initialization, parton evolution, and hadronization with event generation

We present additions to the MARTINI event generator for examining heavy quarks and quarkonia in heavy-ion collisions. All stages of a heavy-ion collision affect the observables associated with heavy quarks: the initial phase space of the heavy quarks are sampled with PYTHIA8.1, the heavy quarks are evolved using Langevin dynamics and a 3+1-dimensional hydrodynamical description of the heavy-ion collision, and are fragmented and hadronized using a modified version of the color evaporation model that takes into account non-trivial evolution in position space, as well as the possibility of recombinant quarkonium production in heavy-ion collisions. We use this to re-examine the production of quarkonium at RHIC, and anticipating vertex detection we predict yields of B_c mesons at RHIC and the LHC.

preprint2011arXiv

Monte-Carlo simulation of jets in heavy-ion collisions

We present Monte-Carlo simulations of jet evolution in lead-lead collisions at the Large Hadron Collider (LHC) at CERN focusing on the dijet asymmetry measured by the ATLAS and CMS collaborations. In the simulation, hard partons are interacting with the hydrodynamical background medium, undergoing radiative and collisional processes. The measured dijet asymmetry is well described by the simulation and is hence consistent with partonic energy loss in a hot, strongly-interacting medium.

preprint2011arXiv

Photons at RHIC: The role of viscosity and of initial state fluctuations

We study real photons produced in heavy ion collisions at RHIC, and we calculate their spectrum and its azimuthal momentum anisotropy. The photons from a variety of sources are included, and the interplay and the time-evolution of those sources are modelled in a full 3D hydrodynamic simulation. We quantify the v_2 of thermal photons produced in ideal and viscous fluids, and the consequences of using different initial conditions are explored.

preprint2011arXiv

Viscous photons in relativistic heavy ion collisions

Theoretical studies of the production of real thermal photons in relativistic heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) are performed. The space-time evolution of the colliding system is modelled using MUSIC, a 3+1D relativistic hydrodynamic simulation, using both its ideal and viscous versions. The inclusive spectrum and its azimuthal angular anisotropy are studied separately, and the relative contributions of the different photon sources are highlighted. It is shown that the photon v2 coefficient is especially sensitive to the details of the microscopic dynamics like the equation of state, the ratio of shear viscosity over entropy density, eta/s, and to the morphology of the initial state.

preprint2010arXiv

3+1D hydrodynamic simulation of relativistic heavy-ion collisions

We present MUSIC, an implementation of the Kurganov-Tadmor algorithm for relativistic 3+1 dimensional fluid dynamics in heavy-ion collision scenarios. This Riemann-solver-free, second-order, high-resolution scheme is characterized by a very small numerical viscosity and its ability to treat shocks and discontinuities very well. We also incorporate a sophisticated algorithm for the determination of the freeze-out surface using a three dimensional triangulation of the hyper-surface. Implementing a recent lattice based equation of state, we compute p_T-spectra and pseudorapidity distributions for Au+Au collisions at root s = 200 GeV and present results for the anisotropic flow coefficients v_2 and v_4 as a function of both p_T and pseudorapidity. We were able to determine v_4 with high numerical precision, finding that it does not strongly depend on the choice of initial condition or equation of state.

preprint2010arXiv

Bulk spectral functions in single and multi-scalar gravity duals

We examine two point correlation functions involving the trace of the energy momentum tensor in five dimensional gravity dual theories supported by one or more scalar fields. A prescription for determining bulk channel spectral functions is developed. This prescription generalizes previous work which centered on one scalar field. As an application of these techniques, we investigate the bulk spectral function and corresponding sum rule in the Chamblin-Reall background. We show that, when expressed in terms of the beta function, the sum rule for the Chamblin-Reall background can be written in a form similar to the sum rule in Yang-Mills theory.

preprint2010arXiv

Elliptic and triangular flow in event-by-event (3+1)D viscous hydrodynamics

We present results for the elliptic and triangular flow coefficients in Au+Au collisions at root-s=200 AGeV using event-by-event (3+1)D viscous hydrodynamic simulations. We study the effect of initial state fluctuations and finite viscosities on the flow coefficients v_2 and v_3 as functions of transverse momentum and pseudo-rapidity. Fluctuations are essential to reproduce the measured centrality dependence of elliptic flow. We argue that simultaneous measurements of v_2 and v_3 can determine eta/s more precisely.

preprint2010arXiv

Five Easy Pieces: The Dynamics of Quarks in Strongly Coupled Plasmas

We revisit the analysis of the drag a massive quark experiences and the wake it creates at a temperature T while moving through a plasma using a gravity dual that captures the renormalisation group runnings in the dual gauge theory. Our gravity dual has a black hole and seven branes embedded via Ouyang embedding, but the geometry is a deformation of the usual conifold metric. In particular the gravity dual has squashed two spheres, and a small resolution at the IR. Using this background we show that the drag of a massive quark receives corrections that are proportional to powers of log T when compared with the drag computed using AdS/QCD correspondence. We use the perturbation produced by the quark strings to compute the wake and compare with the results obtained using AdS/QCD correspondence. We also study the shear viscosity with running couplings, analyze the viscosity to entropy ratio and compare the result with the known bound. In the presence of higher order curvature square corrections from the back-reactions of the embedded D7 branes, we argue the possibility of the entropy to viscosity bound being violated. Finally, we show that our set-up could in-principle allow us to study a family of gauge theories at the boundary by cutting off the dual geometry respectively at various points in the radial direction. All these gauge theories can have well defined UV completions, and more interestingly, we demonstrate that any thermodynamical quantities derived from these theories would be completely independent of the cut-off scale and only depend on the temperature at which we define these theories. Such a result would justify the holographic renormalisabilities of these theories which we, in turn, also demonstrate. We give physical interpretations of these results and compare them with more realistic scenarios.

preprint2010arXiv

Heavy Quarkonium Melting in Large N Thermal QCD

Large N QCD is mostly governed by planar diagrams and should show linear confinement when these diagrams are suitably summed. The linear confinement of quarks in a class of these theories using gravity duals that capture the logarithmic runnings of the coupling constants in the IR and strongly coupled asymptotic conformal behavior in the UV was studied in our previous work. We also extended the theories to high temperatures and argued the possibilities of meltings and suppressions of heavy quarkonium states. In this paper we give a formal proof of melting using very generic choices of UV completions, and point out some subtleties associated with meltings in generic large N theories. Our proof requires only the existence of well defined UV behaviors that are devoid of Landau poles and UV divergences of the Wilson loops, allowing degrees of freedom to increase monotonously with energy scale. We determine the melting temperatures of heavy quarkonium states, which could suggest the presence of deconfinement phase transitions in these theories.

preprint2010arXiv

Toward Large N Thermal QCD from Dual Gravity: The Heavy Quarkonium Potential

We continue our study on the gravity duals for strongly coupled large N QCD with fundamental flavors both at zero and non-zero temperatures. The gravity dual at zero temperature captures the logarithmic runnings of the coupling constants at far IR and the almost conformal, albeit strongly coupled, behavior at the UV. The full UV completion of gauge theory is accomplished in the gravity side by attaching an AdS cap to the IR geometry described in our previous work. Attaching such an AdS cap is highly non-trivial because it amounts to finding the right interpolating geometry and sources that take us from a gravity solution with non-zero three-form fluxes to another one that has almost vanishing three-form fluxes. In this paper we give a concrete realisation of such a scenario, completing the program advocated in our earlier paper. One of the main advantage of having such a background, in addition to providing a dual description of the required gauge theory, is the absence of Landau poles and consequently the UV divergences of the Wilson loops. The potential for the heaviest fundamental quark anti-quark pairs, which are like the heavy quarkonium states in realistic QCD, can be computed and their linear behavior at large separations and zero temperature could be demonstrated. At small separations the expected Coulombic behavior appears to dominate. On the other hand, at non-zero temperatures interesting properties like heavy quarkonium type suppressions and melting are shown to emerge from our gravity dual. We provide some discussions of the melting temperature and compare our results with the Charmonium spectrum and lattice simulations. We argue that, in spite of the large N nature of our construction, certain model-independent predictions can be made.

preprint2009arXiv

Jet energy loss, photon production, and photon-hadron correlations at RHIC

Jet energy loss, photon production and photon-hadron correlations are studied together at high transverse momentum in relativistic heavy-ion collisions at RHIC energies. The modification of hard jets traversing a hot and dense nuclear medium is evaluated by consistently taking into account induced gluon radiation and elastic collisions. The production of high transverse momentum photons in Au+Au collisions at RHIC is calculated by incorporating a complete set of photon-production channels. Comparison with experimental photon production and photon-hadron correlation data is performed, using a (3+1)-dimensional relativistic hydrodynamic description of the thermalized medium created in these collisions. Our results demonstrate that the interaction between the hard jets and the soft medium is important for the study of photon production and of photon-hadron correlation at RHIC.

preprint2009arXiv

MARTINI - Monte Carlo simulation of jet evolution

We present the Modular Algorithm for Relativistic Treatment of heavy IoN Interactions (MARTINI), an event generator for the hard and penetrating probes in high energy nucleus-nucleus collisions. The simulation consists of a time evolution model for the soft background, such as hydrodynamics, PYTHIA 8.1 to generate and hadronize the hard partons after the medium evolution, which is based on the McGill-AMY formalism and includes both radiative and elastic processes. MARTINI allows for the generation of full event configurations in the high transverse momentum region. We present results for the neutral pion and photon nuclear modification factor in Au+Au collisions at RHIC.

preprint2005arXiv

A statistical model analysis of yields and fluctuations in 200 GeV Au-Au collisions

We show that the simultaneous measurement of yields and fluctuations is capable of falsifying and constraining the statistical hadronization model. We show how such a measurement can test for chemical non-equilibrium, and distinguish between a high temperature chemically equilibrated freeze-out from a supercooled freeze-out with an over-saturated phase space. We perform a fit, and show that both yields and fluctuations measured at RHIC 200 GeV can be accounted for within the second scenario, with both the light and strange quark phase space saturated significantly above detailed balance. We point to the simultaneous fit of the $K/π$ fluctuation and the $K^*/K^-$ ratio as evidence that the effect of hadronic re-interactions after freeze-out is small.

preprint2005arXiv

Energy Loss of Leading Hadrons and Direct Photon production in Evolving Quark-Gluon Plasma

We calculate the nuclear modification factor of neutral pions and the photon yield at high p_T in central Au-Au collisions at RHIC (\sqrt{s}=200 GeV) and Pb-Pb collisions at the LHC (\sqrt{s}=5500 GeV). A leading-order accurate treatment of jet energy loss in the medium has been convolved with a physical description of the initial spatial distribution of jets and a (1+1) dimensional expansion. We reproduce the nuclear modification factor of pion R_{AA} at RHIC, assuming an initial temperature T_i=370 MeV and a formation time τ_i=0.26 fm/c, corresponding to dN/dy=1260. The resulting suppression depends on the particle rapidity density dN/dy but weakly on the initial temperature. The jet energy loss treatment is also included in the calculation of high p_T photons. Photons coming from primordial hard N-N scattering are the dominant contribution at RHIC for p_T > 5 GeV, while at the LHC, the range 8<p_T<14 GeV is dominated by jet-photon conversion in the plasma.

preprint2005arXiv

Particle multiplicities and fluctuations in 200 GeV Au-Au collisions

We use the statistical hadronization model (SHM) to describe hadron multiplicity yields and fluctuations. We consider 200 GeV Au-Au collisions, and show that both event averaged yields of stable particles and resonances, and event-by-event fluctuation of the $K/π$ ratio can be described within the SHM using the same set of thermal parameters, provided that the phase space occupancy parameter value is significantly above chemical equilibrium, and the freeze-out temperature is $\sim 140$ MeV. We present predictions that allow to test the consistency of our results.

preprint2004arXiv

Universal Transition Curve in Pseudo-Rapidity Distribution

We show that an unambiguous way of determining the universal limiting fragmentation region is to consider the derivative ($d^2n/dη^2$) of the pseudo-rapidity distribution per participant pair. In addition, we find that the transition region between the fragmentation and the central plateau regions exhibits a second kind of universal behavior that is only apparent in $d^2n/dη^2$. The $\sqrts$ dependence of the height of the central plateau $(dn/dη)_{η=0}$ and the total charged particle multiplicity $n_{\rm total}$ critically depend on the behavior of this universal transition curve. Analyzing available RHIC data, we show that $(dn/dη)_{η=0}$ can be bounded by $\ln^2 s$ and $n_{\rm total}$ can be bounded by $\ln^3 s$. We also show that the deuteron-gold data from RHIC has the exactly same features as the gold-gold data indicating that these universal behaviors are a feature of the initial state parton-nucleus interactions and not a consequence of final state interactions. Predictions for LHC energy are also given.

preprint2001arXiv

Balance Functions, Correlations, Charge Fluctuations and Interferometry

Connections between charge balance functions, charge fluctuations and correlations are presented. It is shown that charge fluctuations can be directly expressed in terms of a balance functions under certain assumptions. The distortion of charge balance functions due to experimental acceptance is discussed and the effects of identical boson interference is illustrated with a simple model.

preprint2000arXiv

Wong's equations and the small x effective action in QCD

We propose a new form for the small x effective action in QCD. This form of the effective action is motivated by Wong's equations for classical, colored particles in non-Abelian background fields. We show that the BFKL equation, which sums leading logarithms in x, is efficiently reproduced with this form of the action. We argue that this form of the action may be particularly useful in computing next-to- leading-order results in QCD at small x.

preprint1999arXiv

Coulomb distortion of pi+/pi- as a tool to determine the fireball radius in central high energy heavy ion collisions

We compute the Coulomb distortion produced by an expanding and highly charged fireball on the spectra of low transverse momenta and mid rapidity pions produced in central high energy heavy ion collisions. We compare to data on Au+Au at 11.6A GeV from E866 at the BNL AGS and of Pb+Pb at 158A GeV from NA44 at the CERN SPS. We match the fireball expanison velocity with the average transverse momentum of protons and find a best fit to the charged pion ratio when the fireball radius is about 10 fm at freeze-out. This value is common to both AGS and SPS data.

preprint1999arXiv

Drell-Yan and J/psi Production in High Energy Proton-Nucleus and Nucleus-Nucleus Collisions

The distributions of outgoing protons and charged hadrons in high energy proton-nucleus collisions are described rather well by a linear extrapolation from proton-proton collisions. This linear extrapolation is applied to precisely measured Drell-Yan cross sections for 800 GeV protons incident on a variety of nuclear targets. The deviation from linear scaling in the atomic number A can be accounted for by energy degradation of the proton as it passes through the nucleus if account is taken of the time delay of particle production due to quantum coherence. We infer an average proper coherence time of 0.4 +/- 0.1 fm/c. Then we apply the linear extrapolation to measured J/psi production cross sections for 200 and 450 GeV/c protons incident on a variety of nuclear targets. Our analysis takes into account energy loss of the beam proton, the time delay of particle production due to quantum coherence, and absorption of the J/psi on nucleons. The best representation is obtained for a coherence time of 0.5 fm/c, which is consistent with Drell-Yan production, and an absorption cross section of 3.6 mb, which is consistent with the value deduced from photoproduction of the J/psi on nuclear targets. Finally, we compare to recent J/psi data from S+U and Pb+Pb collisions at the SPS. The former are reproduced reasonably well with no new parameters, but not the latter.

preprint1997arXiv

Hot Neutron Stars as a Source for Gamma Ray Bursts at Cosmological Distance Scales

We discuss the possibility that the sources for gamma ray bursts are hot neutron stars at cosmological distance scales. The temperature of such stars would be $T \sim 1 \MeV$. Such hot stars can produce an electromagnetic blast wave provided that the ratio of baryon and photon numbers $N_{B}/N_γ \le 10^{-6}$. The typical time scale for such blasts, the total luminosity, and correlation of gamma ray energy with time of arrival are shown to be roughly consistent with observation. The spectrum of photons also appears to be consistent with known data.