580 likes | 816 Views
Space-time Evolution of Bulk QCD Matter at RHIC. Chiho NONAKA Nagoya University. Contents Introduction Hydrodynamic models at RHIC Freezeout process, finite state interactions 3D hydro+UrQMD Results Single particle spectra, elliptic flow Summary.
E N D
Space-time Evolution of Bulk QCD Matter at RHIC Chiho NONAKA Nagoya University • Contents • Introduction • Hydrodynamic models at RHIC • Freezeout process, finite state interactions • 3D hydro+UrQMD • Results • Single particle spectra, elliptic flow • Summary In collaboration with Steffen A. Bass (Duke University & RIKEN BNL) HQ2006
Success of Perfect Hydrodynamic Models at RHIC Single particle spectra Huovinen, Kolb, Heinz, Hirano, Teaney, Shuryak, Hama, Morita, ……. Strong elliptic flow Strong coupled (correlated) QGP • However…. • Elliptic flow vs. Hirano and Tsuda, PRC66 Discrepancy at large : • Insufficient thermalization? • Viscosity effect? • Simple freezeout process? • freezeout • viscosity Hydrodynamic Models at RHIC Huovinen et.al, PLB503 at mid rapidity
Freezeout process in Hydro • Single freezeout temperature? • Difference between chemical andthermal freezeout Heinz,NPA • chemical freezeout • statistical model • Tch ~170 MeV • hadron ratio • thermal freezeout • hydro • Tf ~ 110~140 MeV • Possible solutions: • Partial chemical equilibrium (PCE) • Hirano, Kolb, Rapp • Hydro + Micro Model • Bass, Dumitru, Teaney, Shuryak
Final State Interactions Thermal model T = 177 MeV = 29 MeV UrQMD 3D-Hydro +UrQMD 2. In UrQMD final state interactions are included correctly. Markert @QM2004
UrQMD • Full 3-d Hydrodynamics • EoS :1st order phase transition • QGP + excluded volume model Hadronization Cooper-Frye formula final state interactions Monte Carlo TC TSW t fm/c 3D-Hydro + UrQMD Model Bass and Dumitru, PRC61,064909(2000) Teaney et al, nucl-th/0110037 • Hadron phase: viscosity effect • Freezeout process: • Chemical freezeout & thermal freezeout • Final state interactions 3D-Hydro +UrQMD Key: • Treatment of freezeout is determined by mean free path. • Brake up thermalization: viscosity effect TC:critical temperature TSW:Hydro UrQMD
nucleus nucleus 3-D Hydrodynamic Model • Relativistic hydrodynamic equation • Baryon number conservation • Coordinates • Lagrangian hydrodynamics • Tracing the adiabatic path of each volume element • Effects of phase transition on observables • Computational time • Easy application to LHC • Algorithm • Focusing on the conservation law energy momentum tensor Lagrangian hydrodynamics Flux of fluid
0=0.6 fm/c max=40 GeV/fm3, nBmax=0.15 fm-3 0=0.5 =1.5 Parameters • Initial Conditions • Energy density • Baryon number density • Parameters • Flow • Switching temperature • longitudinal direction • transverse plane Different from Pure Hydro ! vL=Bjorken’s solution); vT=0 TSW=150 [MeV]
PT spectra • PT spectra at central collisions
PT Spectra (Pure Hydro) Tf=110 MeV normalization of K and p: ratio at Tchem Heinz and Kolb, hep-ph/0204061
Rapidity Distribution • Impact parameter dependence of rapidity distributions
Centrality Dependence • Impact parameter dependence of PT
decay Reaction Dynamics • At mid rapidity • Slope becomes steeper
Final State Interactions p, : flatter, pion wind : small cross section
Hydrodynamic expansion final state interaction <PT> vs mass • At mid rapidity
v2 in hadron phase ? • Quark number scaling of v2 v2: early stage of expansion
Reaction Dynamics in v2 II • Hydro+decay • ~ Hydro@TSW • v2 grows in hadron • phase a bit. • v2 builds up in QGP • phase.
Summary • We present the 3-D hydrodynamic +cascade model • 3-D Hydrodynamic Model • Single particle distribution • Centrality dependence • Elliptic flow • Low Tf : single spectra, elliptic flow, hadron ratios Necessity of improvement of freezeout process in Hydro • 3-D Hydro + UrQMD • Single particle distribution • Centrality dependence • Elliptic flow • Different initial conditions from pure Hydro • Hadron ratios Switching temperature from Hydro to UrQMD • Reaction dynamics in UrQMD • ,small cross section, v2:improvement at forward/backward • Work in progress • EoS dependence: ex. lattice QCD • Switching temperature dependence
0=0.6 fm/c max=55 GeV/fm3, nBmax=0.15 fm-3 0=0.5 =1.5 EOS(entropy density) =0 Parameters • Initial Conditions • Energy density • Baryon number density • Parameters (pure hydro) • Flow • Equation of State • 1st order phase transition Bag Model + Exclude volume model • Freezeout Temperature • longitudinal direction • transverse plane Different from hydo + UrQMD ! vL=Bjorken’s solution); vT=0 Tf=110 [MeV]
PT Spectra Tf=110 MeV normalization of K and p: ratio at Tchem Heinz and Kolb, hep-ph/0204061
V2 vs PT • Elliptic Flow b=6.3 fm b=4.5 fm • b=4.5 fm: consistent with experimental data • b=6.3 fm: proton overestimate
v2 vs 3-15 % 15-25 % • Forward/backward rapidity: overestimate • Pure hydro is valid around mid rapidity.
Hydrodynamic Model at RHIC • Relativistic Heavy Ion Collision & Hydrodynamic models • Schematic sketch t collision thermalization hydrodynamical expansion hadronization freeze-out QGP production hadron phase phase transition Hydro Model • initial conditions • parametrization • color glass • condensate… • equation of states • bag model • lattice QCD… • freezeout process • chemical equilibrium • partial chemical • equilibrium • cascade model…
Introduction 1 • Success of Ideal Hydrodynamic Models at RHIC • Single particle spectra • PT spectra up to ~ 2GeV • Huovinen, Kolb, Heinz, Hirano, Teaney, • Shuryak, Hama, Morita, ……. • Strong elliptic flow • strong coupled QGP at mid rapidity Huovinen et.al, PLB503 Nonaka and Bass
Introduction 2 • Success of Ideal Hydrodynamic Models at RHIC However… • Elliptic flow as a function of Hirano and Tsuda, PRC66 Discrepancy at large : • Insufficient thermalization? • Mean free path • Viscosity effect? • freezeout • viscosity
RHIC sQGP V2 vs multiplicity SPS NA49:PRC68,034903(2003)
Hydrodynamic Models Au + Au GeV PHENIX:Nucl.Phys. A757 (2005) 184 PT spectra p • hadron ratio • CE • PCE • RQMD
Hydrodynamic Models Au + Au GeV PHENIX:Nucl.Phys. A757 (2005) 184 Elliptic Flow p CE PCE RQMD
EoS Initial Conditions Freeze-out Hydrodynamic Models Au + Au GeV PHENIX:Nucl.Phys. A757 (2005) 184 p p PT spectra Elliptic Flow [1]Huovinen et.al., PLB(130) [2]Kolb et.al. PRC, [3]Hirano et.al.PRC(130), [4]Teaney et al.
x y Trajectories on the Phase Diagram • Lagrangian hydrodynamics transverse plane (ix,iy,iz) effect of phase transition C.N et al., Eur. Phys.J C17,663(2000)
Numerical Calculation Step 1. Step 2. Step 3. Coordinates move in parallel with baryon number current and entropy density current. local velocity: from hydro eq. temperature and chemical potential CPU time is almost proportional of # of lattice points.
Improvement of freezeout process Summary 1 • Pure hydro with single freezeout temperature Hadron ratio Elliptic flow at forward/backward rapidity
f distribution =0 =2.2 =3.2
f distribution b=2.4 fm b=4.5 fm b=6.3 fm at mid rapidity
Distribution of # of Collisions • finite b: • dN/ncoll becomes small • steep drop