Fluctuation of Conserved Quantities to look for Critical Point in - - PowerPoint PPT Presentation

fluctuation of conserved quantities to look for critical
SMART_READER_LITE
LIVE PREVIEW

Fluctuation of Conserved Quantities to look for Critical Point in - - PowerPoint PPT Presentation

Fluctuation of Conserved Quantities to look for Critical Point in Phase Diagram TGSW2016 Toshihiro Nonaka University of Tsukuba Outline RHIC Beam Energy Scan Phase I Search for Critical Point with Higher order Fluctuations STAR


slide-1
SLIDE 1

Fluctuation of Conserved Quantities to look for Critical Point in Phase Diagram

TGSW2016 Toshihiro Nonaka University of Tsukuba

slide-2
SLIDE 2
  • T. Nonaka, TGSW2016, Sep.18

Outline

2

✓ RHIC Beam Energy Scan Phase I ✓ Search for Critical Point with Higher

  • rder Fluctuations

✓ STAR Detector and Particle Identification ✓ Published results and Recent studies ✓ Beam Energy Scan Phase II

slide-3
SLIDE 3
  • T. Nonaka, TGSW2016, Sep.18

QCD phase diagram

3

✓ Crossover at μB=0

✓ 1st order phase transition at large μB?

✓ Critical point? ✓ Beam Energy Scan Phase I

at RHIC, √sNN=7.7, 11.5, 14.5, 19.6, 27, 39, 62.4 and 200 GeV. μB T

slide-4
SLIDE 4
  • T. Nonaka, TGSW2016, Sep.18

QCD phase diagram

4

μB T

baseline

✦ We measure the higher order fluctuation of conserved quantities as a function of beam energy, and see “non-monotonic” behaviour with respect to the baseline.

  • X. Luo (STAR collaboration) arXiv:1503.02558v2
slide-5
SLIDE 5
  • T. Nonaka, TGSW2016, Sep.18

RHIC Beam Energy Scan Phase I

5

✓ BES I was performed in 2010, 2011 and 2014.

µB , T : J. Cleymans et al., Phys. Rev. C 73, 034905 (2006)

✓ √sNN=14.5 GeV in 2014 in order to fill in the large μB gap between 11.5 and 19.6 GeV.

√sNN (GeV)

Year

Statistics(Million) 0-80%

μB (MeV)

7.7

2010 ~3 422

11.5

2010 ~6.6 316

14.5

2014 ~13 266

19.6

2011 ~15 206

27

2011 ~32 156

39

2010 ~86 112

62.4

2010 ~45 73

200

2010 ~238 24

slide-6
SLIDE 6
  • T. Nonaka, TGSW2016, Sep.18

Higher order fluctuations

6

✓ Moments : Mean(M), sigma(σ), skewness(S) and kurtosis(κ). ✓ S and κ are non-gaussian fluctuations.

κ > 0 κ < 0

skewness→asymmetry kurtosis→shapness

from wikipedia

✓ Cumulant ⇄ Moment ✓ Cumulant : additivity

Volume dependence

✦ Moments and Cumulants are mathematical measures of “shape”

  • f a histogram which probe the fluctuation of observables.
slide-7
SLIDE 7
  • T. Nonaka, TGSW2016, Sep.18

Fluctuations of conserved quantities

7

(1) Sensitive to correlation length (2) Direct comparison with susceptibilities.

Volume dependence can be canceled by taking ratio.

  • M. Cheng et al, PRD 79, 074505 (2009)

✦ Net-baryon, net-charge and net-strangeness

  • No. of “positively charged”

particles in one collision

  • No. of “negatively charged”

particles in one collision

“Net” : positive - negative →neutrons cannot be measured Fill in histograms

  • ver many collisions
slide-8
SLIDE 8
  • T. Nonaka, TGSW2016, Sep.18

Statistical baselines

8

✓ Poisson - Poisson = Skellam ✓ Odd(even) order cumulant of Skellam distribution is difference(sum) between means of two Poissons. N1-N2 N1 N2

μ1 = 10 μ2 = 8

  • No. of negatively

charged particles

μ1, μ2 : mean parameter of Poisson

  • No. of positively

charged particles

slide-9
SLIDE 9
  • T. Nonaka, TGSW2016, Sep.18

Statistical baselines

9

✓ Poisson - Poisson = Skellam ✓ Odd(even) order cumulant of Skellam distribution is difference(sum) between means of two Poissons. N1-N2 N1 N2

μ1 = 10 μ2 = 8

μ1, μ2 : mean parameter of Poisson

  • No. of negatively

charged particles

  • No. of positively

charged particles

slide-10
SLIDE 10
  • T. Nonaka, TGSW2016, Sep.18

Solenoidal Tracker At RHIC

10

3-dimensional trajectory reconstruction of charged particles

Particle identification from energy loss.

Particle identification from mass squared.

slide-11
SLIDE 11
  • T. Nonaka, JPS 2016 fall meeting, Sep.22

Particle Identification

11

✓ Charged particles are counted using the reconstructed tracks by TPC. ✓ Protons can be identified by using dE/dx from TPC. STAR TPC dE/dx Proton Phase Space

slide-12
SLIDE 12
  • T. Nonaka, TGSW2016, Sep.18

Event by event distribution

12

PRL 112, 032302 (2014)

✓ Event by event net-proton distribution. ✓ Low collision energy, small number of antiproton.

slide-13
SLIDE 13
  • T. Nonaka, TGSW2016, Sep.18

Published results in 2014

13

PRL 112, 032302 (2014) PRL 113, 092301 (2014)

✓ It seems to be interesting around 20 GeV for net-proton results. ✓ Net-charge results are consistent with the baseline due to large

  • errors. → A wide distribution gives large statistical errors.

✦ Finite tracking efficiency is corrected.

slide-14
SLIDE 14
  • T. Nonaka, TGSW2016, Sep.18

Published results in 2014

14

PRL 112, 032302 (2014) PRL 113, 092301 (2014)

✓ It seems to be interesting around 20 GeV for net-proton results. ✓ Net-charge results are consistent with the baseline due to large

  • errors. → A wide distribution gives large statistical errors.

✦ Finite tracking efficiency is corrected.

slide-15
SLIDE 15
  • T. Nonaka, TGSW2016, Sep.18

Extending pT coverage

15

✓ pT region can be extended up to 2.0 GeV by using m2 cut from Time Of Flight detector. ✓ We gain factor two (anti)protons with respect to the published results.

slide-16
SLIDE 16
  • T. Nonaka, TGSW2016, Sep.18

Recent results

16

✓ We can obtain larger signals with larger acceptance. ✓ Acceptance is crucial.

STAR Preliminary Net-proton

  • X. Luo, CPOD2014, Bielefeld, Germany

✦ Finite tracking efficiency is corrected.

slide-17
SLIDE 17
  • T. Nonaka, TGSW2016, Sep.18

Recent results

17

σ model, M.A. Stephanov, PRL107, 052301 (2011)

Signal from the critical point?

  • X. Luo (STAR collaboration) arXiv:1503.02558v2

✓ κσ2 (C4/C2) shows a non-monotonic

  • behaviour. The trend is consistent

with the theoretical calculation. ✓ Measurement at the lower energy is important.

✦ Finite tracking efficiency is corrected.

slide-18
SLIDE 18
  • T. Nonaka, TGSW2016, Sep.18

Beam Energy Scan Phase||

18

  • D. Cebra, 34th Reimei workshop, J-PARC, Tokai, Japan

✓ BES II is planned in 2019 and 2020. ✓ Luminosity will be improved with electron cooling system. ✓ Some detector upgrades will be done by BESII. Pseudo-rapidity coverage will be extended from 1.0 to 1.5. ✓ Higher order fluctuation measurement with small errors and large acceptance.

slide-19
SLIDE 19
  • T. Nonaka, TGSW2016, Sep.18

Beam Energy Scan Phase||

19

✓ BES II is planned in 2019 and 2020. ✓ Luminosity will be improved with electron cooling system. ✓ Some detector upgrades will be done by BESII. Pseudo-rapidity coverage will be extended from 1.0 to 1.5. ✓ Higher order fluctuation measurement with small errors and large acceptance.

  • J. Thader, NPA 00(2016)
slide-20
SLIDE 20
  • T. Nonaka, TGSW2016, Sep.18

Summary

20

✓ Beam Energy Scan Phase | was carried out at √sNN = 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4 and 200 GeV in 2010, 2011 and 2014. ✓ STAR experiment has measured up to 4th order fluctuation of net-charge and net-proton multiplicity distributions for searching the critical point. ✓ Net-proton results with extended pT region show the non-monotonic behaviour. However there is still large errors at low beam energies. ✓ Beam Energy Scan Phase II is planned in 2019 and 2020 focusing on low energy region.

slide-21
SLIDE 21

Back up

21

slide-22
SLIDE 22
  • T. Nonaka, TGSW2016, Sep.18

Acceptance dependence (pT)

22

From the talk by X. Luo at CPOD2014, Bielefeld, Germany

slide-23
SLIDE 23
  • T. Nonaka, TGSW2016, Sep.18

Acceptance dependence (y)

23

From the talk by X. Luo at CPOD2014, Bielefeld, Germany

slide-24
SLIDE 24
  • T. Nonaka, TGSW2016, Sep.18

Centrality dependence

24

  • X. Luo (STAR collaboration) arXiv:1503.02558v2
slide-25
SLIDE 25
  • T. Nonaka, TGSW2016, Sep.18

Beam Energy Scan Phase||

25

  • J. Thader, NPA 00(2016)

✓ Some detector upgrades will be done by BESII. Pseudo- rapidity coverage will be extended from 1.0 to 1.5. ✓ Higher order fluctuation measurement with small errors and large signals.

slide-26
SLIDE 26
  • T. Nonaka, TGSW2016, Sep.18

iTPC upgrade

26

From the talk by W. J. Llope at AGS/RHIC Annual User’s Meeting, BNL, June 7, 2016

slide-27
SLIDE 27
  • T. Nonaka, TGSW2016, Sep.18

Analysis technique

27

  • 1. Centrality determination
  • 2. Centrality Bin Width Correction
  • 3. Statistical error calculation

Use charged particles except protons in

  • rder to avoid the auto correlation.

Calculate cumulants at each multiplicity bin in order to suppress the volume fluctuation.

Analysis : |y|<0.5, p and pbar Centrality : |η|<1.0, exclude p and pbar

Data Glauber χ2/NDF=1.9

Refmult3

X.Luo et al. J. Phys.G40,105104(2013)

  • B. Efron,R. Tibshirani, An introduction to the bootstrap,

Chapman & Hall (1993).

✓ Bootstrap ✓ Delta theorem

STAR preliminary

slide-28
SLIDE 28
  • T. Nonaka, TGSW2016, Sep.18

Efficiency correction

28

✓ Based on the assumption of binomial efficiency. ✓ Simple relationship between measured and true factorial moments. ✓ It can be extended to the case of multi-number of phase spaces.

A.Bzdak and V. Koch PRC.86.044904 M.Kitazawa PRC.86.024904

  • X. Luo PRC.91.034907

A.Bzdak and V. Koch PRC.91.027901

slide-29
SLIDE 29
  • T. Nonaka, TGSW2016, Sep.18

Baryon # vs Proton #

29

  • M. Kitazawa, M. Asakawa
  • PRC. 86. 024904

✦ Net-proton cumulants can be corrected to the net-baryon cumulants by assuming the binomial distribution function.