Source author record

Charles Gale

Charles Gale 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

81works
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

81 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

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.

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

A Non-Equilibrium Approach to Photon Emission from the Late Stages of Relativistic Heavy-Ion Collisions

Cross sections for photon production in hadronic scattering processes have been calculated according to an effective chiral field theory. For $π+ ρ\to π+ γ$ and $π+ π\to ρ+ γ$ processes, these cross sections have been implemented into a novel hadronic transport approach (SMASH), which is suitable for collisions at low and intermediate energies. The implementation is verified by systematically comparing the thermal photon rate to theoretical expectations. The photon rates we obtain are compared to previous works, where scattering processes mediated by $ω$ mesons are found to contribute significantly to the total photon production. Finally, the impact of considering the finite width of the $ρ$ meson is investigated, and a significant enhancement of photon production in the low-energy region is observed. This work is the first step towards a consistent treatment of photon emission in hybrid hydrodynamics+transport approaches. The quantification of the importance of the hadronic stage for the resolution of the direct photon flow puzzle is a next step and can be applied to identify equilibrium and non-equilibrium effects in the hadronic afterburner.

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

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 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

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

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.

preprint2020arXiv

The QCD Trace Anomaly at Strong Coupling from M-Theory

In this short note we address the issue of obtaining the temperature dependence of the QCD conformal anomaly from a top-down holographic dual consistent with very recent lattice results for both, T<Tc and T>Tc. As the holographic dual, we use the M-theory uplift of the SYZ type IIA mirror obtained in arXiv:1306.4339 at finite gauge/string coupling (as part of the 'MQGP' limit) of the type IIB holographic dual of large-N thermal QCD of arXiv:0902.1540. We also show that after a tuning of the (small) Ouyang embedding parameter and radius of a blown-up S^2 when expressed in terms of the horizon radius, a QCD deconfinement temperature Tc = 150 MeV from a Hawking-Page phase transition at vanishing baryon chemical potential, consistent with lattice QCD in the heavy-quark limit, can be obtained.

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

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

Infrared Dynamics of a Large N QCD Model, the Massless String Sector and Mesonic Spectra

A consistency check for any UV complete model for large N QCD should be, among other things, the existence of a well-defined vector and scalar mesonic spectra. In this paper, we use our UV complete model in type IIB string theory to study the IR dynamics and use this to predict the mesonic spectra in the dual type IIA side. The advantage of this approach is two-fold: not only will this justify the consistency of the supergravity approach, but it will also give us a way to compare the IR spectra and the model with the ones proposed earlier by Sakai and Sugimoto. Interestingly, the spectra coming from the massless stringy sector are independent of the UV physics, although the massive string sector may pose certain subtleties regarding the UV contributions as well as the mappings to actual QCD. Additionally, we find that a component of the string landscape enters the picture: there are points in the landscape where the spectra can be considerably improved over the existing results in the literature. These points in the landscape in-turn also determine certain background supergravity components and fix various pathologies that eventually lead to a consistent low energy description of the theory.

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

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

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

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

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.

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

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

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

Finite-size effects on the radiative energy loss of a fast parton in hot and dense strongly interacting matter

We consider finite-size effects on the radiative energy loss of a fast parton moving in a finite temperature strongly interacting medium, using the light cone path integral formalism put forward by Zakharov. We present a convenient reformulation of the problem which makes possible its exact numerical analysis. This is done by introducing the concept of a radiation rate in the presence of finite-size effects. This effectively extends the finite-temperature approach of AMY (Arnold, Moore, and Yaffe) to include interference between vacuum and medium radiation. We compare results with those obtained in the regime considered by AMY, with those obtained at leading order in an opacity expansion, and with those obtained deep in the LPM (Landau-Pomeranchuk-Migdal) regime.

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.

preprint2007arXiv

Electromagnetic Radiation from Broken Symmetries in Relativistic Nuclear Collisions

A new channel of direct photon production from a quark gluon plasma (QGP) is explored in the framework of high-temperature QCD. This process appears at next-to-leading order, in the presence of a charge asymmetry in the excited matter. The photon production rate from this new mechanism is suppressed compared to the QCD annihilation and Compton scattering at low baryon density but assumes importance in baryon-rich matter.

preprint2005arXiv

A Transport Model for Nuclear Reactions Induced by Radioactive Beams

Major ingredients of an isospin and momentum dependent transport model for nuclear reactions induced by radioactive beams are outlined. Within the IBUU04 version of this model we study several experimental probes of the equation of state of neutron-rich matter, especially the density dependence of the nuclear symmetry energy. Comparing with the recent experimental data from NSCL/MSU on isospin diffusion, we found a nuclear symmetry energy of $% E_{sym}(ρ)\approx 31.6(ρ/ρ_{0})^{1.05}$ at subnormal densities. Predictions on several observables sensitive to the density dependence of the symmetry energy at supranormal densities accessible at GSI and the planned Rare Isotope Accelerator (RIA) are also made.

preprint2005arXiv

Bremsstrahlung neutrinos from electron-electron scattering in a relativistic degenerate electron plasma

We present a calculation of neutrino pair bremsstrahlung due to electron-electron scattering in a relativistic degenerate plasma of electrons. Proper treatment of the in-medium photon propagator, i.e., inclusion of Debye screening of the longitudinal part and Landau damping of the transverse part, leads to a neutrino emissivity which is several orders of magnitude larger than when Debye screening is imposed for the tranverse part. Our results show that this in-medium process can compete with other sources of neutrino radiation and can, in some cases, even be the dominant neutrino emission mechanism. We also discuss the natural extension to quark-quark bremsstrahlung in gapped and ungapped quark matter.

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.

preprint2003arXiv

Hadronic Production of Thermal Photons

We study the thermal emission of photons from hot and dense strongly interacting hadronic matter at temperatures close to the expected phase transition to the Quark-Gluon Plasma (QGP). Earlier calculations of photon radiation from ensembles of interacting mesons are re-examined with additional constraints, including new production channels as well as an assessment of hadronic form factor effects. Whereas strangeness-induced photon yields turn out to be moderate, the hitherto not considered t-channel exchange of omega mesons is found to contribute appreciably for photon energies above ~1.5 GeV. The role of baryonic effects is assessed using existing many-body calculations of lepton pair production. We argue that our combined results constitute a rather realistic emission rate, appropriate for applications in relativistic heavy-ion collisions. Supplemented with recent evaluations of QGP emission, and an estimate for primordial (hard) production, we compute photon spectra at SPS, RHIC and LHC energies.

preprint1999arXiv

$ρ$ meson broadening and dilepton production in heavy ion collisions

The modification of the width of the rho meson due to in-medium decays and collisions is evaluated. In high temperature and/or high density hadronic matter, the collision width is much larger than the one-loop decay width. The large width of the $ρ$ meson in matter seems to be consistent with some current interpretations of the $e^+ e^-$ mass spectra measured at the CERN/SPS.

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.

preprint1995arXiv

Another Look at ``soft'' lepton pair production in nucleon-nucleon bremsstrahlung

We compare different formalisms for the calculation of lepton pair emission in hadron-hadron collisions and discuss the consequences of the approximations inherent to each of them. Using a Lorentz-covariant and gauge-invariant formalism proposed by Lichard, we calculate lepton pair emission via bremsstrahlung in proton-proton and neutron-proton reactions at energies between 1 and 5 GeV to leading order in the virtual photon four-momentum. These new results suggest that some previous bremsstrahlung calculations based on varieties of the soft-photon approximation might have somewhat overestimated the DLS low-mass dilepton cross sections. We find that the respective intensities of dilepton production through $p p$ and $n p$ bremsstrahlung are energy-dependent and become comparable at $E_{\rm kin} \agt $ 2 GeV. We also calculate the polar angle anisotropy of the lepton spectrum.

preprint1993arXiv

Measuring velocity ratios with correlation functions

We show how to determine the ratio of the transverse velocity of a source to the velocity of emitted particles, using split-bin correlation functions. The technique is to measure $S_2$ and $S_2^ϕ$, subtract the contributions from the single-particle distribution, and take the ratio as the bin size goes to zero. We demonstrate the technique for two cases: each source decays into two particles, and each source emits a large number of particles.