Scott Moreland Ph.D. defense | Duke University | 27 March 2019 Support provided by the Department of Energy National Nuclear Security Administration Stewardship Science Graduate Fellowship (DOE NNSA SSGF) under grant DE-FC52-08NA28752, and by the U.S. Department of Energy (DOE) under grant DE-FG02-05ER41367. Computing resources provided by the Open Science Grid (OSG) and by the National Energy Research Scientific Computing Center (NERSC).
Hadron resonance gas Quark-gluon plasma Increasing temperature Quantum chromodynamics (QCD) predicts that hadrons melt at Tc ∼ 150 MeV to form a deconfined soup of quarks and gluons dubbed the quark-gluon plasma (QGP). J. Scott Moreland (Duke) 1 / 44
what conditions is it formed in nature? How does it recom- bine to form colorless hadrons? Equation of state? Relations between thermal quantities, e.g. P = P( ) Transport properties? Shear/bulk viscosity, probe energy loss, etc J. Scott Moreland (Duke) 2 / 44
a substance depend on the conditions of its environment QCD phase diagram Describes the phases of nuclear matter as a function of the matter’s temperature T and baryon chemical potential µB . Big picture questions • What are the phases of nuclear matter? • Where are they located in the phase diagram? • What are the matter’s intrinsic properties at each point in the phase diagram? When ordinary substances are subjected to variations in tempera- ture or pressure, they will often undergo a phase transition: a physical change from one state to another. At normal atmospheric pressure, for example, water suddenly changes from liquid to vapor as its temperature is raised past 100° C; in a word, it boils. Water also boils if the temperature is held fixed and the pres- sure is lowered—at high altitude, say. The boundary between liquid and vapor for any given substance can be plotted as a curve in its phase diagram, a graph of tem- perature versus pressure. Another curve traces the boundary between solid and liquid. And depending on the substance, still other curves may trace more exotic phase transitions. (Such a phase diagram may also require more exotic variables, as in the figure). One striking fact made apparent by the phase diagram is that the liquid- vapor curve can come to an end. Beyond this “critical point,” the sharp distinction between liquid and vapor is lost, and the transition becomes continuous. The location of this critical point and the phase boundaries represent two of the most fundamental characteristics of any substance. The critical point of water, for a phase explored by the early universe dur- ing the first few micro- seconds after the Big Bang. At low tempera- tures and high baryon density, such as those encountered in the core of neutron stars, the predictions call for color-superconduct- ing phases. The phase transition between a quark-gluon plasma and a gas of ordinary hadrons seems to be continuous for small chemical potential (the dashed line in the figure). However, model studies sug- gest that a critical point appears at higher values of the potential, beyond which the bound- ary between these phases becomes a sharp line (solid line in the figure). Experimentally verifying the location of these fundamental “landmarks” is central to a quantitative understanding ing further experiments in which nuclear matter will be prepared with a broad range of chemical potentials and temperatures, so as to explore the critical point and the Search for the Critical Point: “A Landmark Study” Quark-Gluon Plasma The Phases of QCD Temperature Hadron Gas Early Universe Future FAIR Experiments Future LHC Experiments Nuclear Matter Vacuum Color Superconductor Critical Point Current RHIC Experiments RHIC Energy Scan Crossover Baryon Chemical Potential ~170 MeV 0 MeV 900 MeV 0 MeV Neutron Stars 1st order phase tran sition Schematic QCD phase diagram for nuclear matter. The solid lines show the phase boundaries for the indicated phases. The solid circle depicts the critical point. Possible trajectories for systems created in the QGP phase at different accelerator facilities are also shown. Frontiers of Nuclear Science, A Long Range Plan, arXiv:0809.3137 J. Scott Moreland (Duke) 3 / 44
top energy: compress and heat matter e empera- ndergo nge mal e, water vapor 00° C; if the pres- say. The por for ed as a ph of tem- curve d and bstance, exotic iagram ables, as nt by id- Beyond nction and us. The he a phase explored by the early universe dur- ing the first few micro- seconds after the Big Bang. At low tempera- tures and high baryon density, such as those encountered in the core of neutron stars, the predictions call for color-superconduct- ing phases. The phase transition between a quark-gluon plasma and a gas of ordinary hadrons seems to be continuous for small chemical potential (the dashed line in the figure). However, model studies sug- gest that a critical point appears at higher values of the potential, beyond which the bound- ary between these phases becomes a sharp line (solid line in ing further experiments in which nuclear al Point: “A Landmark Study” Quark-Gluon Plasma The Phases of QCD Temperature Hadron Gas Early Universe Future FAIR Experiments Future LHC Experiments Nuclear Matter Vacuum Color Superconductor Critical Point Current RHIC Experiments RHIC Energy Scan Crossover Baryon Chemical Potential ~170 MeV 0 MeV 900 MeV 0 MeV Neutron Stars 1st order phase tran sition Schematic QCD phase diagram for nuclear matter. The solid lines show the phase boundaries for the indicated phases. The solid circle depicts the critical point. Possible trajectories for systems created in the QGP phase at different accelerator facilities are also shown. Frontiers of Nuclear Science, A Long Range Plan, arXiv:0809.3137 Relativistic Heavy-Ion Collider (RHIC) Large Hadron Collider (LHC) J. Scott Moreland (Duke) 4 / 44
Experiment sees Experiment controls Hidden process projec�les and beam energy Final par�cle proper�es Fig. by Zhi Qiu ALICE Pb-Pb J. Scott Moreland (Duke) 5 / 44
(red glyphs) hydrodynamics (heatmap) initial conditions pre-collision Figure credit Zhi Qiu Initial conditions Starting point of the simulation Minimum requirements: energy density e, flow velocity uµ, initial viscous corrections πµν and Π Each quantity is a three-dimensional field at hydrodynamic starting time τ0 Common simplification: system invariant under limited boosts along the beam direction (boost invariance)—reduces three spatial dimensions to two J. Scott Moreland (Duke) 6 / 44
(red glyphs) hydrodynamics (heatmap) initial conditions pre-collision Figure credit Zhi Qiu Boltzmann transport As system expands and cools, matter begins to break-up and interactions become infrequent Fluid discretized into particles and modeled using non-equilibrium Boltzmann transport: dfi (x, p) dt = Ci (x, p) Kinetic freezeout Last interactions cease, and particles freestream to the “detector” J. Scott Moreland (Duke) 6 / 44
J. S. Moreland, J. E. Bernhard, and S. A. Bass, “Alternative ansatz to wounded nucleon and binary collision scaling in high-energy nuclear collisions”, Phys. Rev. C92, 011901 (2015) arXiv:1412.4708 [nucl-th]
Theory 3 Theory 4 ◎ Experimental data • Color Glass Condensate EFT • pQCD minijet production & gluon saturation • AdS-CFT holography • Color flux tubes and string fragmentation Top-down (data-driven) approach Plausible theory space ◎ Experimental data parameter p • Construct plausible “meta-model” • Parametrize uncertain d.o.f. • Place reasonable prior on parameter • Update prior using global model-to-data comparison (Bayesian approach) J. Scott Moreland (Duke) 10 / 44
Transverse density of participant matter: ˜ TA,B(x⊥) = Npart i=1 γi Tp(x⊥ − x⊥i ± b/2) Random weight γi sampled from Gamma distribution with unit mean and variance 1/k. These ad hoc fluctuations are necessary to describe large minimum-bias proton-proton multiplicity fluctuations. J. Scott Moreland (Duke) 13 / 44
˜ TB Before τ0 dη dx dy ˜ TA ˜ TB After In ultrarelativistic limit, energy and entropy densities are local mappings e0 = fe( ˜ TA, ˜ TB), s0 = fs( ˜ TA, ˜ TB). Scalar functions fe and fs must obey basic physical constraints J. Scott Moreland (Duke) 14 / 44
is participant nucleon scaling e0 s0 ∝ npart A + npart B 2 = ˜ TA + ˜ TB 2 Couldn’t describe multiplicity distribution, so a binary collision term ncoll AB ∝ ˜ TA ˜ TB was added e0 s0 ∝ (1 − α) npart A + npart B 2 + α ncoll AB Parameter α is tuned to fit the data (right) 101 102 103 dNch /dη Pb-Pb, sNN =2.76 TeV dNch /η∝ Npart dNch /η∝ Ncoll Two-component: α=0.08 ALICE, |η| <0.5 0 10 20 30 40 50 60 70 80 Centrality % 0.5 1.0 1.5 Ratio Two-component model fit to Pb-Pb data J. Scott Moreland (Duke) 15 / 44
0.06 ) 2 p(v -1 10 1 10 centrality: 0-1% |<2.5 η >0.5GeV,| T data p 2 ∈ Glauber 0.36 2 ∈ MC-KLN 0.31 ATLAS Pb+Pb =2.76 TeV NN s -1 b µ = 7 int L 0 0.05 0.1 0.15 -1 10 1 10 centrality: 5-10% |<2.5 η >0.5GeV,| T data p 2 ∈ Glauber 0.46 2 ∈ MC-KLN 0.30 ATLAS Pb+Pb =2.76 TeV NN s -1 b µ = 7 int L 0 0.05 0.1 0.15 0.2 0.25 -2 10 -1 10 1 10 centrality: 20-25% |<2.5 η >0.5GeV,| T data p 2 ∈ Glauber 0.40 2 ∈ MC-KLN 0.29 ATLAS Pb+Pb =2.76 TeV NN s -1 b µ = 7 int L 2 v 0 0.1 0.2 ) 2 p(v -2 10 -1 10 1 10 centrality: 30-35% |<2.5 η >0.5GeV,| T data p 2 ∈ Glauber 0.36 2 ∈ MC-KLN 0.28 ATLAS Pb+Pb =2.76 TeV NN s -1 b µ = 7 int L 2 v 0 0.1 0.2 0.3 -2 10 -1 10 1 10 centrality: 40-45% |<2.5 η >0.5GeV,| T data p 2 ∈ Glauber 0.31 2 ∈ MC-KLN 0.26 ATLAS Pb+Pb =2.76 TeV NN s -1 b µ = 7 int L 2 v 0 0.1 0.2 0.3 -3 10 -2 10 -1 10 1 10 centrality: 55-60% |<2.5 η >0.5GeV,| T data p 2 ∈ Glauber 0.24 2 ∈ MC-KLN 0.21 ATLAS Pb+Pb =2.76 TeV NN s -1 b µ = 7 int L ATLAS collaboration JHEP11(2013)183 Two-component ansatz (MC-Glauber) fails to describe the elliptic flow data v2 . Model predicts wrong shape for the event-by-event flow distribution P(v2 ). J. Scott Moreland (Duke) 16 / 44
( ˜ TA, ˜ TB ) Depicted: geometric mean p = 0 Generalized mean ansatz: e0 s0 = Norm × Mp( ˜ TA, ˜ TB), Mp(x, y) ≡ xp + yp 2 1/p Each mean predicts unique geometry: Harmonic p = −1 Geometric p = 0 Arithmetic p = +1 J. Scott Moreland (Duke) 18 / 44
model: Reduced Thickness Event-by-event Nuclear Topology1 Variable model parameters 1. Overall energy (or entropy) normalization 2. Generalized mean parameter p 3. Inverse variance of the nucleon fluctuation gamma weights k 4. Gaussian nucleon width w [fm] 5. Inter-nucleon minimum distance dmin [fm] 1Yes, this is an abuse of terminology J. Scott Moreland (Duke) 19 / 44
0.4 0.6 εn ε2 0 30 60 90 Centrality % 0.0 0.2 0.4 0.6 ε3 Arithmetic: p=1 Geometric: p=0 Harmonic: p= − 1 0 10 20 30 40 Centrality % 0.3 0.6 0.9 1.2 Ratio IP-Glasma • Left and middle: second and third eccentricity harmonics ε2 and ε3 for the Pythagorean means. • Far right: ratio of second and third eccentricities strongly constrained by heavy-ion flow data (gray band) calculated by Retinskaya et al. PRC.89.014902. • TRENTo nicely agrees with flow constraints using geometric mean (p = 0). J. Scott Moreland (Duke) 21 / 44
Bernhard, J. S. Moreland, and S. A. Bass, “Bayesian estimation of the specific shear and bulk viscosity of the quark-gluon plasma”, Nature Physics (submitted), (2019) This study is the capstone of J. Bernhard’s dissertation on Bayesian parameter estimation for heavy-ion collisions. I am not the author of the Bayesian parameter estimation framework.
at 2.76 and 5.02 TeV 3. Generalized mean p for initial energy deposition Harmonic Geometric Arithmetic 4. Std. dev. of nucleon-nucleon multiplicity fluctuations σfluct 5. Gaussian nucleon width w 6. Inter-nucleon minimum distance dmin QGP medium 7. Pre-equilibrium free streaming time τfs 8–10. three parameters (η/s min, slope, and crv) for the temperature dependence of the specific shear viscosity (η/s)(T) 11–13. three parameters (ζ/s max, width, and loc) for the temperature dependence of the specific bulk viscosity (ζ/s)(T) 150 200 250 300 Temperature [MeV] 0.0 0.2 0.4 0.6 η/s 150 200 250 300 Temperature [MeV] 0.00 0.02 0.04 0.06 0.08 0.10 ζ/s 14. Particlization temperature Tswitch J. Scott Moreland (Duke) 24 / 44
2.76 and 5.02 TeV taken from the ALICE experiment Observable Label Particles Reference Integrated Charged particle yield dNch/dη charged [1012.1657] [1512.06104] Identified particle yield dN/dy π, K, p [1303.0737] Transverse energy yield dET /dη [1603.04775] Azimuthal Two-particle flow cumulants vn{k} charged [1105.3865] [1602.01119] Radial Mean transverse momentum pT π, K, p [1303.0737] Transverse momentum fluctuations δpT / p charged [1407.5530] J. Scott Moreland (Duke) 25 / 44
Physics theory initial stages, hydro, and Boltzmann transport Computer model minimum bias event-by- event simulations Gaussian process emulator surrogate model Markov chain Monte Carlo calibrate model to data Posterior distribution quantitative estimates of each parameter Experimental data yields, mean pT , flows J. Scott Moreland (Duke) 27 / 44 Jonah Bernhard, Duke Ph.D. dissertation “Bayesian parameter estimation for heavy-ion collisions” Extensible and adaptable framework
P(E|x) Prior: flat within design range, zero otherwise Likelihood: multivariate Gaussian form P(E|x) = 1 (2π)mdetΣ exp − 1 2 [ym(x) − ye]T Σ−1(x)[ym(x) − ye] Use MCMC to draw samples from unnormalized posterior distribution. Histogram samples to visualize the posterior. J. Scott Moreland (Duke) 28 / 44
Posterior QGP viscosities 1.0 0.5 0.0 0.5 1.0 p KLN EKRT / IP-Glasma Wounded nucleon 150 200 250 300 Temperature [MeV] 0.0 0.1 0.2 0.3 0.4 η/s 1/4π Shear viscosity Posterior median 90% credible region 150 200 250 300 Temperature [MeV] 0.00 0.02 0.04 0.06 0.08 ζ/s Bulk viscosity Evidently, initial energy deposition is well described by geometric mean scaling: e0 ∝ ˜ TA ˜ TB Recall, ˜ TA, ˜ TB are the transverse densities of participant nuclear matter Properly accounting for IC uncertainties (left) allowed for first quantitative estimate of (η/s)(T) and (ζ/s)(T) with meaningful uncertainties (above) These uncertainty bands would be significantly underestimated if we used a bottom-up ab initio initial condition model! J. Scott Moreland (Duke) 31 / 44
Moreland, J. E. Bernhard, and S. A. Bass, “Bayesian calibration of a hybrid nuclear collision model using p-Pb and Pb-Pb data from the LHC”, Phys. Rev. C (submitted), (2019) arXiv:1808.02106 [nucl-th]
-2 0 2 4 (radians) φ ∆ 0 2 4 φ ∆ d η ∆ d pair N 2 d trig N 1 2.4 2.6 2.8 < 260 offline trk N ≤ = 2.76 TeV, 220 NN s (a) CMS PbPb < 3 GeV/c trig T 1 < p < 3 GeV/c assoc T 1 < p η ∆ -4 -2 0 2 4 (radians) φ ∆ 0 2 4 φ ∆ d η ∆ d pair N 2 d trig N 1 3.1 3.2 3.3 3.4 < 260 offline trk N ≤ = 5.02 TeV, 220 NN s (b) CMS pPb < 3 GeV/c trig T 1 < p < 3 GeV/c assoc T 1 < p Pb Pb Pb p • Hydrodynamic description of heavy-ion collisions is now well established. Primary “proof” is remarkable global description of bulk observables, e.g. see previous figure. • The framework surely breaks down at some small length scale. After all, it is an effective theory. A priori, hydro was not expected to describe collisions involving protons and light ions. However, data seems to suggest otherwise! • The exact same collectivity signals used to justify flow in heavy-ion collisions now seen in p-p and p-A collisions. Natural question: Can we apply the same theory machinery, but down to a smaller length scale? J. Scott Moreland (Duke) 33 / 44
light-ion collisions, e.g. p-p and p-Pb, is sensitive to each nucleon’s sampled density profile. • So far I’ve assumed the nucleon is a Gaussian “sphere” for simplicity. • Geometric mean of two Gaussians is also a Gaussian. p=1 p=0 p= − 1 proton-proton collisions → • Flow anisotropy is sensitive to initial state eccentricity, so let’s relax the Gaussian proton assumption. J. Scott Moreland (Duke) 34 / 44
fm .4 fm Lump width Number of lumps Free parameters: • Width of the Gaussian distribution used to sample constituent lump radial positions • Number of Gaussian constituent lumps inside each nucleon • Width of the Gaussian constituent lumps Absent from this work: • Lump spatial correlations
no substructure substructure Transverse proton density 2 fm no substructure substructure Pictured above: transverse nuclear densities T(x⊥ ) = dz ρ(x⊥ , z) J. Scott Moreland (Duke) 36 / 44
at 5.02 TeV 2. Generalized mean p for initial energy deposition Harmonic Geometric Arithmetic 3. Std. dev. of nucleon-nucleon multiplicity fluctuations σfluct 4. Gaussian width of constituent radial position distribution r 5. Gaussian constituent number nc 6. Gaussian width of each constituent v 7. Inter-nucleon minimum distance dmin QGP medium 8. Pre-equilibrium free streaming time τfs 9–11. three parameters (η/s min, slope, and crv) for the temperature dependence of the specific shear viscosity (η/s)(T) 12–14. three parameters (ζ/s max, width, and loc) for the temperature dependence of the specific bulk viscosity (ζ/s)(T) 150 200 250 300 Temperature [MeV] 0.0 0.2 0.4 0.6 η/s 150 200 250 300 Temperature [MeV] 0.00 0.02 0.04 0.06 0.08 0.10 ζ/s 15. Particlization temperature Tswitch J. Scott Moreland (Duke) 37 / 44
QGP viscosities 1.0 0.5 0.0 0.5 1.0 p KLN EKRT / IP-Glasma Wounded nucleon Bayesian calibration on: p-Pb, Pb-Pb 5.02 TeV Pb-Pb 2.76, 5.02 TeV 150 200 250 300 Temperature [MeV] 0.0 0.1 0.2 0.3 0.4 η/s 1/4π 1/4π Shear viscosity 90% credible region Posterior median 150 200 250 300 Temperature [MeV] 0.00 0.02 0.04 0.06 0.08 ζ/s Bulk viscosity Bayesian calibration on: p-Pb, Pb-Pb 5.02 TeV Pb-Pb 2.76, 5.02 TeV Uncertainties are larger, but nevertheless consistent with geometric mean scaling e0 ∝ ( ˜ TA ˜ TB )1/2 as before. Note, this was not guaranteed! If model was wildly wrong, unlikely that similar scaling would hold. Again, viscosity estimates self-consistent with Pb-Pb calibration to 2.76 and 5.02 TeV data without substructure. J. Scott Moreland (Duke) 41 / 44
Constituent sampling radius r [fm] 0.2 0.4 0.6 0.8 1.0 1.2 Constituent width v [fm] 1 3 5 7 9 Constituent number nc Sampling radius r varies constituent dispersion: position distribution → P(x) ∝ exp(−|x|2/2r2) Width v varies each constituent’s size: density → ρ(x) ∝ exp(−|x|2/2v2) Parameter nc varies number of constituents (lumps) inside each nucleon. Evidently nc > 1 favored, but no strong preference for specific number. Measure proton size and shape with ultrarelativistic nuclear collisions? J. Scott Moreland (Duke) 42 / 44
• Applied model to extract QGP transport coefficients from combined p-Pb, Pb-Pb analysis • Demonstrated hydrodynamics is able to simultaneously describe bulk observables in small & large collision systems Harmonic Geometric Arithmetic 150 200 250 300 Temperature [MeV] 0.0 0.1 0.2 0.3 0.4 η/s 1/4π 1/4π Shear viscosity 90% credible region Posterior median 150 200 250 300 Temperature [MeV] 0.00 0.02 0.04 0.06 0.08 ζ/s Bulk viscosity Bayesian calibration on: p-Pb, Pb-Pb 5.02 TeV Pb-Pb 2.76, 5.02 TeV J. Scott Moreland (Duke) 43 / 44
thermal equilibrium using lattice discretization techniques • Method limited to zero baryon chemical potential µB = 0 • Provides QCD equation of state relating energy density, entropy density, pressure, etc. at fixed temperature and baryon chemical potential • Lattice calc. indicate phase transition to QGP is smooth crossover at Tc ∼ 150 MeV for µB = 0 (ε-3p)/T4 p/T4 s/4T4 0 1 2 3 4 130 170 210 250 290 330 370 T [MeV] stout HISQ HotQCD and Wuppertal-Budapest lattice QCD equations of state, 10.1103/PhysRevD.90.094503 J. Scott Moreland (Duke) 1 / 8
the 14-momentum approximation: τΠΠ + ˙ Π = −ζθ − δΠΠΠθ + λΠππµνσµν, (1a) τπ ˙ π µν + πµν = 2ησµν − δπππµνθ + φ7π µ α πν α − τπππ µ α σν α + λπΠΠσµν. (1b) Fluid discretized into particles using the Cooper-Frye formula: E dNi d3p = gi (2π)3 Σ fi (x, p) pµ d3σµ, (2) where i is an index over species, fi the particle species’ distribution function, and d3σµ a volume element (located at spacetime position x) of the isothermal hypersurface Σ defined by Tswitch . J. Scott Moreland (Duke) 2 / 8
on links connecting lattice sites (right). Start by calculating vacuum-to-vacuum transition amplitude Z = DAa µ (x)D ¯ Ψ(x)DΨ(x)ei d4xL[Aa µ ,¯ Ψ,Ψ] Replace time t everywhere with imaginary time, t → iτ Z = DAa µ (x, τ)D ¯ Ψ(x, τ)DΨ(x, τ)e− β 0 dτ d3xLE [Aa µ ,¯ Ψ,Ψ] Associate with the partition function Z = trˆ ρ where ˆ ρ = e−β ˆ H is the density operator for grand canonical ensemble in therm equilibrium J. Scott Moreland (Duke) 3 / 8
layers from inner to outer: first tracking system, then EM and hadronic calorimeters, and finally muon tracking system. • Magnetic field bends tracks of charged particles for momentum and charge determination • Complemented by 18 detectors to specify mass, velocity, and electric sign. • Inner tracking system: (1) silicon pixel, drift, and strip detectors, (2) time projection chamber, and (3) transition radiation detector. • Particle identification: charged hadrons specified by mass & charge. Mass deduced from momentum and velocity. Momentum and charge sign from curvature of track in TPC. Velocity from time-of-flight and radiation characteristics. J. Scott Moreland (Duke) 4 / 8
streaming: • Massless non-interacting parton gas matched to viscous hydrodynamics PRC.91.064906, PRC.80.034902. • Parametrize initial energy density e0 • Initialize hydro with non-zero uµ and πµν $bl; Ə bm= orѴbm] =u;;v|u;-l _7uo =v J. Scott Moreland (Duke) 5 / 8