Multi-particle production in small systems from CGC
Prithwish Tribedy
1
7th International Workshop on Multiple Partonic Interactions at the LHC
The Abdus Salam International Centre for Theoretical Physics (ICTP), Trieste, Italy
Multi-particle production in small systems from CGC Prithwish - - PowerPoint PPT Presentation
Multi-particle production in small systems from CGC Prithwish Tribedy 7th International Workshop on Multiple Partonic Interactions at the LHC The Abdus Salam International Centre for Theoretical Physics (ICTP), Trieste, Italy 1 outline
Prithwish Tribedy
1
7th International Workshop on Multiple Partonic Interactions at the LHC
The Abdus Salam International Centre for Theoretical Physics (ICTP), Trieste, Italy
2
Based on the work done in collaboration with :
Goal : Study correlated production of particles We need :
* dN dy1 d2p⊥1 . . . dyq d2p⊥q + ( ) ⌧ dN dy1 d2p⊥1
* dN dyq d2p⊥q +
3
Focus : Collisions of small systems p+p and p+Pb are interesting as final state effects are minimal
Origin of high multiplicity events
4 arXiv: 1011.5531
Origin of high multiplicity events Systematics of Δη-Δφ correlations Energy dependence of ridge in p+p Similar underlying dynamics must drive these phenomenon
5
p+p p+A
arXiv: 1011.5531 arXiv: 1509.04776, 1210.5482
6
Multi-particle production at high energies in Regge Gribov limit (x→0)
Colliding hadrons/nuclei :
hard saturation scale Qs(x) > ΛQCD
arXiv: 1212.1701
7
Multi-particle production at high energies in Regge Gribov limit (x→0)
Colliding hadrons/nuclei :
semi-hard saturation scale Qs(x) > ΛQCD
0.2 0.4 0.6 0.8 1 1.2 0.1 1 10 Q0
2=0.168 GeV2
Y=0 Y=4 Y=8 Y=12
Qs
N
r e g i m e
linear
un-integrated gluon distribution arXiv: 1212.1701 Dusling, Li, Schenke 1509.07939 parton transverse momentum (kT) GeV
8
Multi-particle production at high energies in Regge Gribov limit (x→0)
Particle production :
like emissions of gluons,
coupling, high occupation of gluonic states f(k) ~ A2~1/g2
Initial configuration JIMWLK evolution dN/d p 3
Color Glass condensate effective field theory → ab-inito framework to this problem
9
Multi-particle production at high energies in Regge Gribov limit (x→0)
Particle production :
like emissions of gluons,
coupling, high occupancy of gluonic states ~1/g2
(classical approximation)
Initial configuration JIMWLK evolution Single gluon emission A (classical field)
McLerran, Venugopalan hep-ph/9309289
Charge density matrices ρa(x⊥,Y) Local Gaussian distribution W[ρ] (MV-Model)
Yang Mills equations [Dμ,Fμν] = Jν for each configuration of source
ρ(x⊥)10
Domains of chromo-electric field
⌦ ρa(x⊥)ρb(y⊥) ↵ = δabδ2(x⊥−y⊥)g 2µ2(x⊥)
Glasma flux tubes —> free streaming gluons
before collisions (τ<0)
hρρi
A
h i hρρi
B
hV †V i
A
hV †V i
B
classical color charge classical color field
D F = J D F = J D F = J
B
A
after collisions (τ>0)
D F = 0 hep-ph/9809433, hep-ph/0303076, arXiv: 1206.6805 arXiv: 1202.6646
Input is constrained by dipole-cross sections in e+p/A collisions Perturbative approach
systems Non-perturbative approach
Monte-Carlo model of initial conditions
11
D
dN dypd2p⊥
E ⇠ ⌦ |M|2↵ ⇠ hρ∗
1ρ1ρ∗ 2ρ2i
M ⇠ ρ1(k⊥) k⊥2 ρ2(p⊥ k⊥) (p⊥ k⊥)2 Lγ(p, k⊥)
Color Averaging
12
Single-Inclusive
p q
C2(p, q) ≡ ⌧ dN2 dypd2p⊥dyqd2q⊥
⌧ dN dypd2p⊥ ⌧ dN dyqd2q⊥
↓ connected disconnected
⊥ ⊥ ⊥ ⊥
⌦ |M|2↵ ! hρ∗
1ρ∗ 1ρ1ρ1ρ∗ 2ρ∗ 2ρ2ρ2i )
⌦ ↵
Double-Inclusive
p q
8 topologies 1 topology It can be shown Dumitru, Gelis, McLerran, Venugopalan 0804.3858
p
1
p
2 ...
p
q
p
q ...
p
2
p
1
p
1 ...
p
q
p
q ...
p
1
Naturally generates Negative Binomial distribution probability distribution
13
CGC framework is extendable to n-particle correlations
High-multiplicity events —> originate from correlated production of n-particles —> Highly non-perturbative
P
NB
n
= Γ(k + n) Γ(k)Γ(n + 1) ¯ nnkk (¯ n + k)n+k
k = κ(Nc
2 − 1)Qs 2S⊥
2π
2n(n-1)! topologies
Gelis, Lappi, McLerran 0905.3234
and impact parameter
configurations
14
IP-Glasma model : combines CGC framework & different sources of initial state fluctuations
Making Nucleus out of proton scattering
SA
dip(r⊥, x, b⊥) = A
Y
i=0
Sp
dip(r⊥, x, b⊥)
Si
p
)
Tp(s⊥) = 1 2πBG exp ✓s⊥2 2BG ◆ Z Tpp(b) = Z d2s⊥ T A
p (s⊥) T B p (s⊥ b⊥).
dP d2b(b) = 1 e−σggN2
gTpp(b)
R d2b ⇣ 1 e−σggN2
g Tpp(b)⌘,
Overlap function Impact parameter distribution
15
color charge distribution in nucleus
Proton profile
Q2
sA ⇠ A1/3Q2 sp
Nuclear saturation scale :
Schenke, Tribedy, Venugopalan 1311.3636
/dy
g
dN
5 10 15 20
Entries
10
2
10
3
10
Entries0<b<0.5 fm 0.5<b<1.0 fm 1.0<b<1.5 fm 1.5<b<2.0 fm
= 0.48 fm τ p+p 7 TeV,
IP-Glasma
For a given geometry fluctuations of color charge —> Negative Binomial distribution at each impact parameter
16
However the distribution is not wide enough to describe data Some sources of fluctuation missing
10-4 10-3 10-2 10-1 100 101 2 4 6 8 10
P(Nch/〈 Nch 〉)
Nch/〈 Nch 〉
IP-Glasma
p+p 7 TeV CMS
Convolution of many NBDs
Input to CGC framework —> dipole cross section e+p/A
r q q z 1-z
* γ
s
With evolution of rapidity each dipole split with probability ~ αs dY —> dipole splitting is however stochastic Color dipole picture : distribution of partons —> dist. of color dipoles Stochastic dipole splitting —> not present in BK/JIMWLK —>beyond CGC
17
2log(r /r )
2α
2α
2r
i 2log(r /r )
2T dipoles T
saturation
r
i1
2log(r /r )
2α
2α
2r
i 2log(r /r )
2T dipoles T
saturation
r
i1
2log(r /r )
2α
2α
2r
i 2log(r /r )
2T dipoles T
saturation
r
i1
2log(r /r )
2α
2α
2r
i 2log(r /r )
2T dipoles T
saturation
r
i1
dipole-probe target
Iancu, Mueller, Munier (hep-ph/0410018) Golec-Biernat, Wusthoff hep-ph/9807513
10-4 10-3 10-2 10-1 100 101 0.5 1 1.5 2 2.5 3 P(QS/〈 QS 〉) QS/〈 QS 〉 σ=0.4 σ=0.5
P(ln(Q2
S/hQ2 Si)) =
1 p 2πσ exp ✓ ln2(Q2
S(s⊥)/hQ2 S(s⊥)i)
2σ2 ◆
σ2(Y ) = σ2
0(Y0) + σ2 1(Y Y0),
0.2 0.4 0.6 0.8 1 0.1 1 10 T(r,Y) log(r0
2/r2)
Y=8 Y=0
Dipole amplitude Saturation scale Stochastic splitting of dipole leads to a distribution of Qs
18
Marquet, Soyez, Xiao hep-ph/0606233
19
10-4 10-3 10-2 10-1 100 101 2 4 6 8 10
P(Nch/〈 Nch 〉)
Nch/〈 Nch 〉
IP-Glasma σ=0.4 p+p 200 GeV UA5
10-4 10-3 10-2 10-1 100 101 2 4 6 8 10
P(Nch/〈 Nch 〉)
Nch/〈 Nch 〉
IP-Glasma σ=0 IP-Glasma σ=0.5
p+p 7 TeV CMS
pp@LHC pp@RHIC
Origin of High multiplicity events (Tail of distributions)
High multiplicity events —> rare configuration of high color charge density (1/g2)
McLerran, Tribedy 1508.03292
∼ Q−1
s ~ E
20
correlation from initial state
color sources/domains Very distinct from Hydrodynamic flow (driven by geometry )
Kovner, Lublinsky 1012.3398 Lappi, Schenke, Schlichting, Venugopalan 1509.03499 Dumitru, Giannini 1406.5781 Dumitru, Dusling, Gelis, Jalilian-Marian, . Lappi, Venugopalan 1009.5295 Dusling, Venugopalan 1201.2658 Kovchegov, Wertepny 1212.1195
Δφ = π trigger
back-to-back in 0 & π
π/2 π
Jet Graph
Y(Δφ)
21
Δφ
Kinematically constrained (back-to-back)
Δφ = π trigger
Symmetric around π/2 Glasma Graph Not kinematically constrained
Dusling, Venugopalan 1201.2658, 1210.3890
Di-Jet Graph
22
π/2 π
Y(Δφ) Δφ
Jet + BFKL emissions
gluon emissions between two triggered hadrons—> broadening of the away side (de-correlation)
η ∆
2 4 φ ∆ 2 4
φ ∆ d η ∆ d
pair
N
2
d
trg
N 1
1.30 1.35 1.40
CMS Preliminary 110 ≥ = 7 TeV, N s pp <3 GeV/c
trig T
2<p <2 GeV/c
assoc T
1<p
φ ∆ 2 4
1.30 1.35 1.40
1 NTrig d2N d∆η d∆φ
∆φ π
1 NTrig d2N d∆φ
q p
Glasma Graph
q p
Jet Graph
23
But why ridge appears in high multiplicity events ?
In CGC, high occupancy ~1/g2 —> effective coupling 1/g2 x g = 1/g
Dusling, Li, Schenke 1509.07939
24
1/g2 g2 g g2 g g2 g2 1/g2 g2 1/g 1/g low density (min-bias events) high density (high multiplicity events) g2 strong color field in CGC g—> ρg~1/g2 g~1/g
0.5 1 1.5 2 2.5 3
1.0 < pT
trig < 2.0 GeV; 1.0 < pT asc < 2.0 GeVATLAS Central ATLAS Peripheral
π π ∆φ
25
Consistent explanation in the CGC picture
Dusling, Venugopalan 1201.2658, 1210.3890, 1211.3701, 1302.7018
0.02 0.04 0.06 0.08 0.1 0.12 0.5 1 1.5 2 2.5 3 d2N/d∆φ ∆φ p+p s1/2 = 7 TeV BFKL + glasma Q2
s0 = 1.008 GeV2
BFKL + glasma Q2
s0 = 0.840 GeV2
BFKL + glasma Q2
s0 = 0.672 GeV2
CMS: Ntrk>110, 1 GeV < pT
a,b < 2 GeV
p+p p+Pb
figure: Dusling, Li, Schenke 1509.07939
0.01 0.02 0.03 0.04 0.05 20 40 60 80 100 120 140
Yint
Nch
rec
1< pasc,trg
T < 2 GeV
Glasma + BFKL
13 TeV (ATLAS acc.) 7 TeV (CMS acc. × 3.6) 13 TeV ATLAS Prelim. 7 TeV CMS × 3.6
Data
0.01 0.02 0.03 0.04 0.05 1 2 3 4 5
Yint
p
T
asc,trg [GeV]
Data: Nch
rec>110
13 TeV ATLAS 7 TeV CMS × 3.6
Glasma + BFKL
13 TeV (ATLAS acc.) 7 TeV (CMS acc. × 3.6) Nch
rec ~ 125
Nch
rec ~ 100
26
ch (√s) = Nrec ch (Qs2) , Yint(√s) = Yint(Qs2)
Energy dependence enters only though the saturation scale (only scale in problem) : Scaling of near side yield is natural in CGC approach
Dusling, Tribedy, Venugopalan 1509.04410
0.02 0.04 0.06 0.08 0.1 1 2 3 4 5 6 7 8 v3 pT [GeV]
Gluons τ=0.4 fm/c ATLAS v3(2PC) 110 < Nch
rec < 140CMS v3(2PC) 120 < Ntrk
Gluons τ=0.0 fm/c v3(EP/2PC) Gluons τ=0.2 fm/c v3(EP) v3(2PC) v3(EP)
0.05 0.1 0.15 0.2 1 2 3 4 5 6 7 8 v2 pT [GeV]
Gluons τ=0.2 fm/c ATLAS v2(2PC) 110 < Nch
rec < 140CMS v2{4} 120 < Ntrk
Gluons τ=0.0 fm/c v2(2PC) v2(2PC) v2(EP)
27
…. + +
π/2 π Y(Δφ) Δφ π/2 π Y(Δφ) Δφ
even harmonics
Schenke, Schlichting, Venugopalan 1502.01331
CGC
space to color neutral strings
final particles
28
rS
z t
1 QS
x x+
A = pure gauge 1
−
A = pure gauge 1 A = 0
T )
fragmentation
gluons quarks anti-quarks strings
px py px py y y
29
Work in progress
Connect the gluons close in phase space to color neutral strings
30 10-4 10-3 10-2 10-1 100 101 1 2 3 4 5 6 7
P(N/〈 N 〉)
N/〈 N 〉
CMS Data Gluons Hadrons
IP-Glasma + PYTHIA 8.2 (p+p 7 TeV)
Preliminary results
31
The ab-initio framework of CGC constrained by HERA DIS data provide successful description of the phenomena seen in high multiplicity events at LHC
32
33
p q
k k k k p − k p − k q − k q − k
Momentum flow in Glasma graph (origin of ridge-like correlation)
Dusling, Li, Schenke 1509.07939
34
Correlation strength Multiplicity for fixed system size
mini jets glasma graphs hydrodynamics parton escape initial state correlations response to initial geometry
8 Soeren Schlichting Quark Matter 2015