Source author record

Bjoern Schenke

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

52works
5topics
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

52 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

Rapidity dependence of initial state geometry and momentum correlations in p+Pb collisions

Event geometry and initial state correlations have been invoked as possible explanations of long range azimuthal correlations observed in high multiplicity p+p and p+Pb collisions. We study the rapidity dependence of initial state momentum correlations and event-by-event geometry in $\sqrt{s}=5.02~\rm{TeV}$ p+Pb collisions within the 3+1D IP-Glasma model~\cite{Schenke:2016ksl}, where the longitudinal structure is governed by JIMWLK rapidity evolution of the incoming nuclear gluon distributions. We find that the event geometry is correlated across large rapidity intervals whereas initial state momentum correlations are relatively short range in rapidity. Based on our results, we discuss implications for the relevance of both effects in explaining the origin of collective phenomena in small systems.

preprint2021arXiv

Energy dependence of the chiral magnetic effect in expanding holographic plasma

Based on a holographic far-from-equilibrium calculation of the chiral magnetic effect~(CME) in an expanding quark gluon plasma, we study collisions at various energies. We compute the time evolution of the CME current in the presence of a time-dependent axial charge density and subject to a time-dependent magnetic field. The plasma expansion leads to a dilution of the CME current. We study distinct combinations of how the initial magnetic field and initial axial charge behave with changing initial energy as proposed in previous literature. Most scenarios we consider lead to an increasing time-integrated CME current, when increasing the initial energy. This would make it more likely to observe the CME at higher collision energies.

preprint2020arXiv

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

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

preprint2020arXiv

Diffractive dijet production in impact parameter dependent saturation models

We study coherent diffractive dijet production in electron-hadron and electron-nucleus collisions within the dipole picture. We provide semi-analytic results for the differential cross section and elliptic anisotropy in the angle between dijet transverse momentum and hadron recoil momentum. We demonstrate the direct relation between angular moments of the dipole amplitude in coordinate space and angular moments of the diffractive dijet cross sections. To perform explicit calculations we employ two different saturation models, extended to include the target geometry. In the limit of large photon virtuality or quark masses, we find fully analytic results that allow direct insight into how the differential cross section and elliptic anisotropy depend on the saturation scale, target geometry, and kinematic variables. We further provide numerical results for more general kinematics in collisions at a future electron-ion collider, and study the effects of approaching the saturated regime on diffractive dijet observables.

preprint2020arXiv

Inclusive prompt photon-jet correlations as a probe of gluon saturation in electron-nucleus scattering at small $x$

We compute the differential cross-section for inclusive prompt photon$+$quark production in deeply inelastic scattering of electrons off nuclei at small $x$ ($e+A$ DIS) in the framework of the Color Glass Condensate effective field theory. The result is expressed as a convolution of the leading order (in the strong coupling $α_{\mathrm{s}}$) impact factor for the process and universal dipole matrix elements, in the limit of hard photon transverse momentum relative to the nuclear saturation scale $Q_{s,A}(x)$. We perform a numerical study of this process for the kinematics of the Electron-Ion Collider (EIC), exploring in particular the azimuthal angle correlations between the final state photon and quark. We observe a systematic suppression and broadening pattern of the back-to-back peak in the relative azimuthal angle distribution, as the saturation scale is increased by replacing proton targets with gold nuclei. Our results suggest that photon+jet final states in inclusive $e+A$ DIS at high energies are in general a promising channel for exploring gluon saturation that is complementary to inclusive and diffractive dijet production. They also provide a sensitive empirical test of the universality of dipole matrix elements when compared to identical measurements in proton-nucleus collisions. However because photon+jet correlations at small $x$ in EIC kinematics require jet reconstruction at small $k_\perp$, it will be important to study their feasibility relative to photon-hadron correlations.

preprint2020arXiv

Initial State Summary of Hard Probes 2020

The description of the initial state of heavy ion collisions, which covers the description of the incoming nuclei, the initial hard and soft interactions, the resulting spatial geometry of the produced matter, as well as the dynamic approach to a medium well described by hydrodynamics, has important consequences for the study of hard and electromagnetic probes. I will review new developments presented at Hard Probes 2020 that have an impact on these aspects of our understanding of the initial state of heavy ion and smaller system collisions.

preprint2019arXiv

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

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

preprint2016arXiv

3-D Glasma initial state for relativistic heavy ion collisions

We extend the impact parameter dependent Glasma model (IP-Glasma) to three dimensions using explicit small x evolution of the two incoming nuclear gluon distributions. We compute rapidity distributions of produced gluons and the early time energy momentum tensor as a function of space-time rapidity and transverse coordinates. We study rapidity correlations and fluctuations of the initial geometry and multiplicity distributions and compare to existing models for the three dimensional initial state.

preprint2016arXiv

Initial State Quantum Fluctuations in the Little Bang

We review recent developments in the ab-initio theoretical description of the initial state in heavy-ion collisions. We emphasize the importance of fluctuations, both for the phenomenological description of experimental data from the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC), and the theoretical understanding of the non-equilibrium early time dynamics and thermalization of the medium.

preprint2016arXiv

Mass ordering of spectra from fragmentation of saturated gluon states in high multiplicity proton-proton collisions

The mass ordering of mean transverse momentum $\left<p_T\right>$ and of the Fourier harmonic coefficient $v_2 (p_T)$ of azimuthally anisotropic particle distributions in high energy hadron collisions is often interpreted as evidence for the hydrodynamic flow of the matter produced. We investigate an alternative initial state interpretation of this pattern in high multiplicity proton-proton collisions at the LHC. The QCD Yang-Mills equations describing the dynamics of saturated gluons are solved numerically with initial conditions obtained from the Color Glass Condensate based IP-Glasma model. The gluons are subsequently fragmented into various hadron species employing the well established Lund string fragmentation algorithm of the PYTHIA event generator. We find that this ab initio initial state approach reproduces characteristic features of bulk spectra, in particular the particle mass dependence of $\left<p_T\right>$ and $v_2 (p_T)$.

preprint2016arXiv

Novel Collective Phenomena in High-Energy Proton-Proton and Proton-Nucleus Collisions

The observation of long-range collective correlations for particles emitted in high-multiplicity pp and pPb collisions has opened up new opportunities of investigating novel high-density QCD phenomena in small colliding systems. We review experimental results related to the studies of collective phenomena in small systems from RHIC and the LHC over the past several years. Latest development in theoretical interpretations motivated by different frameworks are also reviewed, and confronted with the experimental data. Perspectives on possible future directions are discussed, with the aim of further exploring the rich emergent QCD phenomena.

preprint2015arXiv

A Tale of Tails: Photon Rates and Flow in Ultra-Relativistic Heavy Ion Collisions

We consider the possibility that quark and gluon distributions in the medium created in high energy heavy ion collisions may be modified by a power law tail at energies much higher than the temperature. We parametrize such a tail by Tsallis distributions with an exponent motivated by phenomenology. These distributions are characterized by an effective temperature scale that we assume to evolve in time like the temperature for thermal distributions. We find that including such a tail increases the rates for photon production and significantly delays the emission times for photons of a fixed energy. We argue that these effects are sufficiently large that they should be able to account for photon yields and flow patterns seen in LHC and RHIC experiments.

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

Moving forward to constrain the shear viscosity of QCD matter

We demonstrate that measurements of rapidity differential anisotropic flow in heavy ion collisions can constrain the temperature dependence of the shear viscosity to entropy density ratio η/s of QCD matter. Comparing results from hydrodynamic calculations with experimental data from RHIC, we find evidence for a small η/s $\approx$ 0.04 in the QCD cross-over region and a strong temperature dependence in the hadronic phase. A temperature independent η/s is disfavored by the data. We further show that measurements of the event-by-event flow as a function of rapidity can be used to independently constrain the initial state fluctuations in three dimensions and the temperature dependent transport properties of QCD matter.

preprint2015arXiv

New insights from 3D simulations of heavy ion collisions

Viscous relativistic hydrodynamics in 3+1 dimensions is applied to describe heavy ion collisions at RHIC and LHC. We present calculations of observables that are sensitive to the longitudinal structure of the created system. In particular we present pseudo-rapidity correlations and demonstrate their dependence on both the initial state and short range correlations introduced via a microscopic transport description. We further demonstrate the effect of a varying temperature dependence of the shear viscosity to entropy density ratio on rapidity dependent flow harmonics.

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

Pseudorapidity correlations in heavy ion collisions from viscous fluid dynamics

We demonstrate by explicit calculations in 3+1 dimensional viscous relativistic fluid dynamics how two-particle pseudorapidity correlation functions in heavy ion collisions at the LHC and RHIC depend on the number of particle producing sources and the transport properties of the produced medium. In particular, we present results for the Legendre coefficients of the two-particle pseudorapidity correlation function in Pb+Pb collisions at 2760 GeV and Au+Au collisions at 200 GeV from viscous hydrodynamics with three dimensionally fluctuating initial conditions. Our results suggest that these coefficients provide important constraints on initial state fluctuations and the transport properties of the quark gluon plasma.

preprint2015arXiv

Rapidity dependence in holographic heavy ion collisions

We present an attempt to closely mimic the initial stage of heavy ion collisions within holography, assuming a decoupling of longitudinal and transverse dynamics in the very early stage. We subsequently evolve the obtained initial state using state-of-the-art hydrodynamic simulations, and compare results to experimental data. We present results for charged hadron pseudo-rapidity spectra and directed and elliptic flow as functions of pseudo-rapidity for $\sqrt{s_{NN}}=200\,{\rm GeV}$ Au-Au and $2.76\,{\rm TeV}$ Pb-Pb collisions. The directed flow interestingly turns out to be quite sensitive to the viscosity. The results can explain qualitative features of the collisions, but the rapidity spectra in our current model is narrower than the experimental data.

preprint2014arXiv

Collective effects in light-heavy ion collisions

We present results for the azimuthal anisotropy of charged hadron distributions in A+A, p+A, d+A, and He3+A collisions within the IP-Glasma+MUSIC model. Obtained anisotropies are due to the fluid dynamic response of the system to the fluctuating initial geometry of the interaction region. While the elliptic and triangular anisotropies in peripheral Pb+Pb collisions at root-s=2.76 TeV are well described by the model, the same quantities in root-s=5.02 TeV p+Pb collisions underestimate the experimental data. This disagreement can be due to neglected initial state correlations or the lack of a detailed description of the fluctuating spatial structure of the proton, or both. We further present predictions for azimuthal anisotropies in p+Au, d+Au, and He3+Au collisions at root-s=200 GeV. For d+Au and 3He+Au collisions we expect the detailed substructure of the nucleon to become less important.

preprint2014arXiv

Eccentric protons? Sensitivity of flow to system size and shape in p+p, p+Pb and Pb+Pb collisions

We determine the transverse system size of the initial non-equilibrium Glasma state and of the hydrodynamically evolving fireball as a function of produced charged particles in p+p, p+Pb and Pb+Pb collisions at the Large Hadron Collider. Our results are consistent with recent measurements of Hanbury-Brown-Twiss (HBT) radii by the ALICE collaboration. Azimuthal anisotropy coefficients v_n generated by combining the early time Glasma dynamics with viscous fluid dynamics in Pb+Pb collisions are in excellent agreement with experimental data for a wide range of centralities. In particular, event-by-event distributions of the v_n agree with the experimental data out to fairly peripheral centrality bins. In striking contrast, our results for p+Pb collisions significantly underestimate the magnitude and do not reproduce the centrality dependence of data for v_2 and v_3 coefficients. We argue that the measured v_n data and HBT radii strongly constrain the shapes of initial parton distributions across system sizes that would be compatible with a flow interpretation in p+Pb collisions. Alternately, additional sources of correlations may be required to describe the systematics of long range rapidity correlations in p+p and p+Pb collisions.

preprint2014arXiv

Initial state geometry and fluctuations in Au+Au, Cu+Au and U+U collisions at RHIC

We study within the IP-Glasma and two-component MC-Glauber models the effects of initial state geometry and fluctuations on multiplicities and eccentricities for several collision species at the Relativistic Heavy Ion Collider (RHIC). These include copper-gold (Cu+Au), gold-gold (Au+Au), and uranium-uranium (U+U) collisions. The multiplicity densities per participant pair are very similar in all systems studied. Ellipticities vary strongly between collision systems, most significantly for central collisions, while fluctuation driven odd moments vary little between systems. Event-by-event distributions of eccentricities in mid-central collisions are wider in Cu+Au relative to Au+Au and U+U systems. An anti-correlation between multiplicity and eccentricity is observed in ultra central U+U collisions which is weaker in the IP-Glasma model than the two-component MC-Glauber model. In ultra central Au+Au collisions the two models predict opposite signs for the slope of this correlation. Measurements of elliptic flow as a function of multiplicity in such central events can therefore be used to discriminate between models with qualitatively different particle production mechanisms.

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

The Glasma, Photons and the Implications of Anisotropy

We introduce distribution functions for quarks and gluons in the Glasma and discuss how they satisfy various relationships of statistical physics. We use these distributions to compute photon production in the early stages of heavy ion collisions. Photon rates satisfy geometric scaling, that is, the emission rate per unit area scales as a function of the saturation momenta divided by the transverse momentum of the photon. Photon distributions from the Glasma are steeper than those computed in the Thermalized Quark Gluon Plasma (TQGP). Both the delayed equilibration of the Glasma and a possible anisotropy in the pressure lead to slower expansion and mean times of photon emission of fixed energy are increased. This delayed emission might allow for larger photon elliptic flow.

preprint2014arXiv

The shape of the proton at high energies

We present first calculations of the fluctuating gluon distribution in a proton as a function of impact parameter and rapidity employing the functional Langevin form of the JIMWLK renormalization group equation. We demonstrate that when including effects of confinement by screening the long range Coulomb field of the color charges, the evolution is unitary. The large-x structure of the proton, characterized by the position of three valence quarks, retains an effect on the proton shape down to very small values of x. We determine the dipole scattering amplitude as a function of impact parameter and dipole size and extract the rapidity evolution of the saturation scale and the proton radius.

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

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

Initial state geometry and the role of hydrodynamics in proton-proton, proton-nucleus and deuteron-nucleus collisions

We apply the successful Monte Carlo Glauber and IP-Glasma initial state models of heavy ion collisions to the much smaller size systems produced in proton-proton, proton-nucleus and deuteron- nucleus collisions. We observe a significantly greater sensitivity of the initial state geometry to details of multi-particle production in these models compared to nucleus-nucleus collisions. In particular, we find that the size of the system produced in p+A collisions is very similar to the one produced in p+p collisions, and predict comparable Hanbury-Brown-Twiss radii in the absence of flow in both systems. Differences in the eccentricities computed in the models are large, while differences amongst the generated flow coefficients v_2 and v_3 are smaller. For a large number of participants in proton-lead collisions, the v_2 generated in the IP-Glasma model is comparable to the value obtained in proton-proton collisions. Viscous corrections to flow are large over characteristic lifetimes in the smaller size systems. In contrast, viscous contributions are significantly diminished over the longer space-time evolution of a heavy ion collision.

preprint2013arXiv

Multiplicity distributions in p+p, p+A and A+A collisions from Yang-Mills dynamics

We compute transverse momentum and momentum integrated multiplicity distributions consistently in the IP-Glasma model for proton-proton and proton-lead collisions at the LHC, in deuteron-gold collisions at RHIC, and in heavy ion collisions at both RHIC and LHC energies. Several sources of sub-nucleon scale contributions to the multiplicity distributions are identified. Our results, which are constrained by inclusive and diffractive deeply inelastic scattering data from HERA, are compared to measured distributions for a range of collision energies. These results are an essential first step in quantifying the relative role of initial and final state effects on multiparticle correlations in light and heavy ion collisions.

preprint2013arXiv

Transverse Momentum of Protons, Pions and Kaons in High Multiplicity pp and pA Collisions: Evidence for the Color Glass Condensate?

The CMS experiment at the Large Hadron Collider has recently presented data on the average transverse momentum of protons, kaons and pions as a function of particle multiplicity. We relate the particle multiplicity to saturation momentum using recently published computations of the interaction radius determined from the theory of the Color Glass Condensate. We show that the pp and the pA experimental data scale in terms of these saturation momenta. Computing transverse momentum spectra for identified particles using Boltzmann-type distributions and relating different associated multiplicities using geometric scaling, these simple distributions reproduce the observed dependence of the mean transverse momentum on particle multiplicities seen in both pp and pA interactions for pions to good accuracy, and to fair agreement for protons and kaons.

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

Event-by-event gluon multiplicity, energy density and eccentricities at RHIC and LHC

The event-by-event multiplicity distribution, the energy densities and energy density weighted eccentricity moments epsilon_n (up to n=6) at early times in heavy-ion collisions at both RHIC (root-s=200 GeV) and LHC (root-s=2.76 TeV) are computed in the IP-Glasma model. This framework combines the impact parameter dependent saturation model (IP-Sat) for nucleon parton distributions (constrained by HERA deeply inelastic scattering data) with an event-by-event classical Yang-Mills description of early-time gluon fields in heavy-ion collisions. The model produces multiplicity distributions that are convolutions of negative binomial distributions without further assumptions or parameters. The eccentricity moments are compared to the MC-KLN model; a noteworthy feature is that fluctuation dominated odd moments are consistently larger than in the MC-KLN model.

preprint2012arXiv

Fluctuating Glasma initial conditions and flow in heavy ion collisions

We compute initial conditions in heavy-ion collisions within the Color Glass Condensate (CGC) framework by combining the impact parameter dependent saturation model (IP-Sat) with the classical Yang-Mills description of initial Glasma fields. In addition to fluctuations of nucleon positions, this IP-Glasma description includes quantum fluctuations of color charges on the length-scale determined by the inverse nuclear saturation scale Q_s. The model naturally produces initial energy fluctuations that are described by a negative binomial distribution. The ratio of triangularity to eccentricity is close to that in a model tuned to reproduce experimental flow data. We compare transverse momentum spectra and v_(2,3,4)(p_T) of pions from different models of initial conditions using relativistic viscous hydrodynamic evolution.

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.

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

Flow in heavy-ion collisions - Theory Perspective

I review recent developments in the field of relativistic hydrodynamics and its application to the bulk dynamics in heavy-ion collisions at the Relativistic Heavy- Ion Collider (RHIC) and the Large Hadron Collider (LHC). In particular, I report on progress in going beyond second order relativistic viscous hydrodynamics for conformal fluids, including temperature dependent shear viscosity to entropy density ratios, as well as coupling hydrodynamic calculations to microscopic hadronic rescattering models. I describe event-by-event hydrodynamic simulations and their ability to compute higher harmonic flow coefficients. Combined comparisons of all harmonics to recent experimental data from both RHIC and LHC will potentially allow to determine the desired details of the initial state and the medium properties of the quark-gluon plasma produced in heavy-ion collisions.

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

Renormalization group evolution of multi-gluon correlators in high energy QCD

Many-body QCD in leading high energy Regge asymptotics is described by the Balitsky-JIMWLK hierarchy of renormalization group equations for the x evolution of multi-point Wilson line correlators. These correlators are universal and ubiquitous in final states in deeply inelastic scattering and hadronic collisions. For instance, recently measured di-hadron correlations at forward rapidity in deuteron-gold collisions at the Relativistic Heavy Ion Collider (RHIC) are sensitive to four and six point correlators of Wilson lines in the small x color fields of the dense nuclear target. We evaluate these correlators numerically by solving the functional Langevin equation that describes the Balitsky-JIMWLK hierarchy. We compare the results to mean-field Gaussian and large N_c approximations used in previous phenomenological studies. We comment on the implications of our results for quantitative studies of multi-gluon final states in high energy QCD.

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

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.

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.

preprint2009arXiv

Transverse momentum diffusion and jet energy loss in non-Abelian plasmas

We consider momentum broadening and energy loss of high momentum partons in a hot non-Abelian plasma due to collisions. We solve the coupled system of Wong-Yang-Mills equations on a lattice in real time, including binary hard elastic collisions among the partons. The collision kernel is constructed such that the total collisional energy loss and momentum broadening are lattice spacing independent. We find that the transport coefficient $\hat{q}$ corresponding to transverse momentum broadening receives sizable contributions from a power-law tail in the $p_\perp$-distribution of high-momentum partons. We establish the scaling of $\hat{q}$ and of $dE/dx$ with density, temperature and energy in the weak-coupling regime. We also estimate the nuclear modification factor $R_{AA}$ due to elastic energy loss of a jet in a classical Yang-Mills field.