1 / 63

本審査・公開発表 (10/31/2007)

Proton and Antiproton Production in High Energy Heavy Ion Collisions at RHIC (RHIC での高エネルギー重イオン衝突における陽子反陽子生成). 金野正裕 (数理物質科学研究科). 本審査・公開発表 (10/31/2007). Outline. 1. Introduction 2. Motivation - Baryon/Meson difference - Hadron production in heavy ion collisions 3. Methods

berke
Télécharger la présentation

本審査・公開発表 (10/31/2007)

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Proton and Antiproton Production in High Energy Heavy Ion Collisions at RHIC (RHICでの高エネルギー重イオン衝突における陽子反陽子生成) 金野正裕 (数理物質科学研究科) 本審査・公開発表 (10/31/2007)

  2. Outline 1. Introduction 2. Motivation - Baryon/Meson difference - Hadron production in heavy ion collisions 3. Methods - PHENIX detector - Charged hadron measurement & PID 4. Results & Discussions - Freeze-out properties - (Anti-)Proton production at intermediate pT 5. Conclusion

  3. Introduction - QGP • Quarks can exist as an apparent degree of freedom? • => Quark Gluon Plasma (QGP) • - Matter under high temperature and high energy density. • - Quarks and gluons are freely moving in a large volume. • Relativistic heavy ion collisions is a method to approach the QGP. • QCD transition : Hadron gas <=> QGP Tc ~ 175 MeV εc ~ 0.6 GeV/fm3 QCD phase diagram JHEP 04 (2004) 050 (predicted by lattice QCD calculations) PLB 478 (2000) 447

  4. Relativistic Heavy Ion Collider (RHIC) • RHIC located at Brookhaven National Lab (USA) • A first collider for heavy ion beam • 2 circulating rings (circumference: 3.83 km) • Colliding nuclei: Au+Au, Cu+Cu, d+Au, p+p • Top energy (Au+Au): sNN = 200 GeV • Peak luminosity: ~ 3x1027 cm-2s-1 • Experiments: PHENIX, STAR, BRAHMS, PHOBOS

  5. Hadron gas QGP Pre equilibrium Incoming nuclei Picture of Relativistic Heavy Ion Collisions Freeze-out Hadronic scatterings Hadronization Expansion & Cooling Thermalization Hard scattering Initial collision Space-time evolution of a heavy-ion collision (time scale: ~10 fm/c) There are some stages and dynamic changes.

  6. RHIC Findings (1) Rapidity distributions of charged hadron multiplicity(PRL 91, 052303 (2003)) • Longitudinal expansion (parallel to beam axis) is dominant. • Longitudinal boost invariance partly holds at mid-rapidity. • Hydrodynamic calculations reproduce • elliptic flow behavior at low pT. • Small viscosity (/s) estimated • => Nearly perfect fluid A bulk system is expanding longitudinally & transversely. Phys. Rev. C 67 (03) 044903, Phys.Rev.Lett.91 182301 (2003)

  7. RHIC Findings (2) Jet Quenching • In central Au+Au collisions, hadrons are suppressed at high pT. • - The suppression is a final state effect (parton energy loss). • Away-side jet peak disappeared in central Au+Au collisions.

  8. 1.0 0.5 0 RHIC Findings (3) Particles & Medium Effects Baryon enhanced B/M Splitting of v2 - Suppression/Enhancement has particle-type dependence. => Baryon/Meson difference in yields and emission patterns at intermediate pT(2-5 GeV/c).

  9. Hadron Production in RHI Collisions Current understanding: - There are multiple hadronization mechanisms at intermediate pT. - The relative contributions and particle-type dependence are not yet fully understood.

  10. Motivation Outstanding questions: - Understanding Baryon/Meson difference at intermediate pT. => What is the origin? - What pT does hydrodynamic contribution exist up to? - Quark recombination process is really necessary? - Can we separate hadron radial flow and quark radial flow ? What we should do: - Measurement of Proton and Antiproton pT Spectra - Sensitive to collective flow due to its relatively large mass. - Indicator of baryon number transport at lower energies. => Enhance the high-pT PID capability with new detector. - Systematic Study - Au+Au/Cu+Cu/p+p collision systems at √sNN = 62.4/200 GeV (system size, energy dependence).

  11. My Activities • Construction and Installation of • Aerogel Cherenkov detector in PHENIX • for high-pT PID upgrade (2002~2004) • - Participation in PHENIX experiment • during data taking periods (Run3-Run7) • - Staying at BNL for ~2 years • Data analysis (2005~current): • + Calibrations, Software developments • + Proton spectra using Aerogel detector • Presentations • (QM05, HQ06, QM06, etc.) • Papers, proceedings • (NPA 774 (2006) 461, EPJC 49, 29 (2007)) • (Preparing full papers)

  12. PHENIX EM Calorimeter (PID) TOF (PID) • Global detectors • (event characterization) • - Central Arm Detectors (||<0.35) • (magnetic spectrometer) Aerogel Cherenkov (PID) Pad Chambers (tracking) Drift Chamber (momentum meas.) Global detectors: Beam Beam Counter (trigger, centrality, t0, z-vertex, RP) (efficiency: A+A: ~90%, p+p: ~50%) Zero Degree Calorimeter (centrality)

  13. Data Analysis Data sets: Au+Au at sNN = 62.4, 200 GeV (Run-4) Cu+Cu at sNN = 62.4, 200 GeV (Run-5) p+p at sNN = 200 GeV (Run-5) p+p at sNN = 62.4 GeV (Run-6) Analysis methods: (1) Event selection (z-vertex, centrality) (2) Tracking, Momentum determination (3) Track selection (4) Particle Identification (TOF, ACC) (5) MC corrections (acceptance, efficiency) => Invariant yield pT distributions (/K/p) at mid rapidity ||<0.35

  14. Spectator Participant Spectator Event Selection - Minimum Bias Trigger (BBC coin.) - Centrality determination (BBC, ZDC) Participant-Spectator model ZDC BBC

  15. y R x r  z Track Reconstruction • Drift chamber provides 12 hits in (x,y) plane • Giving the bending angle  (R=220cm) • after passing in magnetic field • Giving pT with field-integral value • PC1 hits and collision z-vertex fix the polar angle  • Momentum resolution: • Find intersection points between the trajectory and • outer detectors. Projected points are then matched • to measured points.

  16. Track Selection DC tracks at PC3 • Residual distribution between hit point • and projection point. • Centroid and width are parameterized • as a function of pT(position ~8 mm at r = 5m). • Require tracks to be within 2. Background subtraction for charged hadron measurement • Matching residual distribution has a tail. • Asymmetric shape comes from residual bend. • Background is subtracted with shape • of the distribution. MC study was done. Background sources (dominant at high pT): - e+, e- from  conversion in materials - Weak decays mostly K+, K-

  17. pT [GeV/c] 0 2 4 6 8 TOF ACC p K+ Clear proton line + proton & antiproton ID Particle Identification Time of Flight (~120 ps), p(p) ID up to 4 GeV/c m2 distributions (3.5-4.0 GeV/c) Aerogel Cherenkov (n=1.011), p(p) ID up to 7 GeV/c Veto for proton ID

  18. Monte Carlo Corrections Acceptance (TOF) Corrections for: - Geometrical acceptance - Decays in flight - Momentum resolution - Detector efficiency - Occupancy effect (tracking efficiency is reduced in high multiplicity environment.) Acceptance (ACC) * MC simulation based on Geant-3. Real data / MC matching: - Dead areas are removed - Detector stability is checked - Same cuts are applied to obtain efficiency Occupancy (TOF, ACC)

  19. Feed-down from weak decays Proton and Lambda pT spectra Feed-down: Weak decays take place very close to the vertex. Since heavier particles take most of the decay momentum, these tracks are inseparable from tracks coming from the vertex of a collision. p+p 200 GeV TOF Evaluation of the fraction: - Effective lambda spectra measured including higher resonances (~33%) - Decays in PHENIX acceptance (MC) - Fraction in measured p(p) : ~15% Fraction of Feed-down ( from ’s)

  20. Systematic Errors Evaluation: Systematic errors are evaluated by varying cut conditions in data analysis. Some parts can be canceled when taking particle ratios etc. Systematic errors (TOF) Systematic errors (ACC)

  21. Proton and Antiproton pT spectra Au+Au sNN = 200 GeV Cu+Cu sNN = 200 GeV p+p sNN = 200 GeV NOTE: No weak decay feed-down correction applied. pT reach extended up to 6 GeV/c for p(p) with fine centrality bins. (1) Aerogel Cherenkov (2) High statistics

  22. Freeze-out Properties Characterizing bulk properties: - Chemical Freeze-out - Kinetic Freeze-out

  23. Particle Yield dN/dy at mid rapidity  K p Au+Au/Cu+Cu/p+p (sNN = 200 GeV) Au+Au/Cu+Cu/p+p (sNN = 62.4 GeV) • Particle yields are (roughly) scaled with Npart btw. Au+Au and Cu+Cu. • dN/dy(Cu+Cu) >~ dN/dy(Au+Au) at smaller Npart. • Statistical model describes their ratios with few parameters (T,).

  24. Statistical Model Fit Phys. Rev. C71 054901, 2005 nucl-th/0405068 • Extracting chemical freeze-out properties with statistical model fit. • Fitting particle ratios of dN/dy (/K/p) at y~0. • Assuming local chemical equilibrium of light quarks (u,d,s), s=1. • Partial feed-down correction taken into account. • - Tch, q : relatively stable • s, s : not determined with this set of ratios (/K/p). Strangeness info is short. Au+Au 62.4 GeV (0-5%) -/+: 0.84 +/- 0.04 K+/K-: 1.19 +/- 0.06 p/p: 0.48 +/- 0.03 K-/-: 0.17 +/- 0.02 p/-: 0.08 +/- 0.01 Au+Au 200 GeV (0-5%) -/+: 1.02 +/- 0.05 K+/K-: 1.09 +/- 0.05 p/p: 0.74 +/- 0.05 K-/-: 0.16 +/- 0.02 p/-: 0.08 +/- 0.01 data model data model 1.01 +/- 0.01 1.20 +/- 0.13 0.48 +/- 0.09 0.17 +/- 0.02 0.08 +/- 0.02 1.00 +/- 0.01 1.09 +/- 0.08 0.74 +/- 0.08 0.16 +/- 0.02 0.08 +/- 0.02 Tch: 167 +/- 10 MeV q: 24 +/- 3 MeV 2/ndf: 9.2/2 Tch: 157 +/- 8 MeV q: 9 +/- 1 MeV 2/ndf: 1.1/2

  25. pp annihilation rate limiting case thermal motion (+ collective flow) T = 120 MeV p ~ 500 MeV,  ~ 0.4 ann ~ 100 mb small for high-p particles cross section dNch/d ~ 700 at mid rapidity dNp/d ~ 18 (2.5 %) • V = 8.4 x 103 fm3 (freeze-out, 10 fm/c) •  ~ 8 • n = dNp/d/(V/) ~ 1.7 x 10-2 fm-3 • = 1/n ~ 5.9 fm particle density survival rate (hadronic phase) exp(- (t)/) = exp(- (0.4x4)/5.9) ~ 0.76

  26. Tch ~160 MeV Chemical Freeze-out Temperature Au+Au/Cu+Cu/p+p (sNN = 200 GeV) Au+Au/Cu+Cu/p+p (sNN = 62.4 GeV) • Tch ~160 MeV, flat centrality dependence. • Npart scaling of Tch between Au+Au and Cu+Cu. • Almost same Tch at √sNN = 62.4, 200 GeV.

  27. Mean Transverse Momentum p K  Au+Au/Cu+Cu/p+p (sNN = 200 GeV) Au+Au/Cu+Cu/p+p (sNN = 62.4 GeV) • Clear hadron mass dependence: larger <pT> for heavier particles. • => Consistent with radial flow picture. • <pT> increases with Npart. it is clearly seen for (anti)proton.

  28. Blast-wave Model Fit Blast-wave model is a parameterization within a simple boost-invariant model with transverse collective flow. pT spectra reflecting thermal freeze-out temperature and transverse flow at final state. 2 map * Ref: PRC48(1993)2462 Tfo ~120 MeV, T ~0.7 (* Resonance decay feed-down correction not applied. Instead, tighter pT fitting range used. ; 0.6-1.2 GeV/c K; 0.4-1.4 GeV/c, p/pbar; 0.6-1.7 GeV/c) Spectra for heavier particles has a convex shape due to radial flow.

  29. <T> ~0.5 Transverse Flow Velocity Au+Au/Cu+Cu/p+p (sNN = 200 GeV) Au+Au/Cu+Cu/p+p (sNN = 62.4 GeV) • <T>: increasing with Npart. • - Npart scaling of <T> between Au+Au and Cu+Cu. • Almost same <T> at √sNN = 62.4, 200 GeV.

  30. Tfo ~120 MeV Kinetic Freeze-out Temperature Au+Au/Cu+Cu/p+p (sNN = 200 GeV) Au+Au/Cu+Cu/p+p (sNN = 62.4 GeV) • Tfo: decreasing with Npart. • Npart scaling of Tfo between Au+Au and Cu+Cu. • Almost same Tfo at √sNN = 62.4, 200 GeV.

  31. Summary - Freeze-out properties Characterizing bulk properties: - Chemical freeze-out - Kinetic freeze-out => Hadron production at low pT : “Thermal emission + Radial flow” Scaling properties between different systems: - Chemical/kinetic freeze-out properties show similarities between different collision systems. - Npart scaling of freeze-out properties (Au+Au, Cu+Cu), * even though the overlapped region has a different shape. => System volume Npart is a control parameter. * Particle yield: (Cu+Cu) > (Au+Au) at smaller Npart - Similarity at sNN = 200 and 62.4 GeV.

  32. Baryon Enhancement

  33. Baryon enhancement at sNN = 200 GeV p/ p/ • (Anti-)proton enhancement observed/confirmed in 200 GeV Au+Au/Cu+Cu. • Larger than expected from jet fragmentation (measured in pp, e+e-). • p/ (p/) ratios turn over at 2~3 GeV/c ,and fall towards the ratio in p+p.

  34. Baryon enhancement at sNN = 62.4 GeV p/ p/ • (Anti-)proton enhancement observed/confirmed in 62.4 GeV Au+Au/Cu+Cu. • Similar pT dependence as at 200 GeV.

  35. p/ ratio vs. Npart1/3 Cu+Cu vs. Au+Au (62.4 GeV) Cu+Cu vs. Au+Au (200 GeV) • Npart scaling of p/ (p/) at same √sNN. • The ratios are controlled by the initial overlap size of colliding nuclei, • even though overlap region has a different geometrical shape.

  36. Nuclear Modification Factor RAA • Comparison with p+p spectra (reference) in binary collision scaling. • Proton, antiproton are enhanced at 1.5 - 4 GeV/c for all centralities. • - Suppression is seen for , K.

  37. Comparison of RAA in Au+Au/Cu+Cu Pion RAA (pT=2.25 GeV/c) Proton RAA (pT=2.25 GeV/c) RAA (Cu+Cu) > RAA (Au+Au) - Proton is enhanced for all centralities, while /K are suppressed.

  38. RAA (Cu+Cu) > RAA (Au+Au) • - Geometrical shape : Au+Au more deformed • - No. of N-N scatterings per N : narrow peak in Cu+Cu Comparison of Au+Au and Cu+Cu Glauber model calculation (Cu+Cu: b=0.0 fm, Au+Au: b=8.6 fm) <Npart> ~117 Cu+Cu: good resolution at smaller Npart Even though Ncoll-Npart relation is almost same between Au+Au and Cu+Cu, the geometrical overlap shape is different. <Npart> ~100

  39. Beam energy dependence

  40. Beam energy dependence of enhancement p/ p/ * No weak decay feed-down correction applied. - p/+ ratio:decreasing as a function of sNN. - p/-ratio: increasing as a function sNN. • Antiproton is a good probe to study the baryon enhancement.

  41. p/ ratio vs. (dET/d)1/3 • No Npart scaling of p/ (p/) in Au+Au between 62.4 and 200 GeV. • Transverse energy density dET/d scaling of p/ is favored. • - dET/d is a connection key between different √sNN. Proton production at 62.4 GeV is partly from baryon number transport, not only proton-antiproton pair production.

  42. Net proton distribution BRAHMS, PRL 93 (2004) 102301 it drastically changes with beam energy. Energy loss per nucleon: 73±6 GeV

  43. Chemical Potential Au+Au/Cu+Cu/p+p (sNN = 200 GeV) Au+Au/Cu+Cu/p+p (sNN = 62.4 GeV) • q (200 GeV) : ~8 MeV, independent of Npart • q (62.4 GeV) : increasing with Npart • => more baryon stopping at central

  44. Summary - Baryon enhancement Baryon enhancement: - Proton and antiproton enhancement confirmed at intermediate pT (2-5 GeV/c) in Au+Au/Cu+Cu. A turnover of p/ ratio seen at pT = 2-3 GeV/c. - In terms of binary collision scaling, (anti)protons are enhanced while pions/kaons are suppressed. Low energy 62.4 GeV data: - At lower energy 62.4 GeV, proton production seems to be more affected by baryon number transport process. => Antiproton is a good indicator of the baryon enhancement. Scaling properties between different systems: - Npart scaling of p/ (p/) - dET/d scaling of p/

  45. Two-component model(Soft+Hard)

  46. Particle production in expanding matter time high-pT particles low-pT particles z-axis time x-axis

  47. Two-component Model (Soft+Hard) Soft component :Thermal emission + Radial flow - Described by Blast-wave model - Npart scaling seen - Thermal distribution extrapolated up to high pT Hard component : Jet fragmentation + Jet suppression - Measured p+p spectra - Ncoll scaling - Constant suppression factor (power-law distribution & fractional energy loss)

  48. 200 GeV 62.4 GeV Hard component in p+p and Au+Au p+p sNN = 200 GeV Au+Au 200 GeV (pi0: diamond, h+h-: circle) • Hard component (in p+p) at high pT depends on s. • In Au+Au, suppression effect should be taken into account.

  49. Pion pT spectra Au+Au 200 GeV + - Soft Line Hard Line Reproduce the measured pion pT spectra.

  50. Pion fraction vs. pT Au+Au 200 GeV + - Hard Soft Residual

More Related