QCD
in a strong magnetic field
Yoshimasa Hidaka (RIKEN)
Part I: magnetic field and anomaly
QCD in a strong magnetic field Part I: magnetic field and anomaly - - PowerPoint PPT Presentation
QCD in a strong magnetic field Part I: magnetic field and anomaly Yoshimasa Hidaka (RIKEN) Introduction 450,000G Orders of magnitude for magnetic fields wikipedia Typical magnet 50G Neodymium magnet 12,500G (strongest permanent magnet)
Yoshimasa Hidaka (RIKEN)
Part I: magnetic field and anomaly
Orders of magnitude for magnetic fields
Neodymium magnet Typical magnet
50G 12,500G
(strongest permanent magnet)
Strongest continuous magnetic field
produced in a laboratory
450,000G
Magnetars Heavy ion collisions The early Universe
(Electroweak transition)
~104MeV2~1017 G ~1022 G ~1013 G
wikipedia
Orders of magnitude for magnetic fields
Neodymium magnet Typical magnet
50G 12,500G
(strongest permanent magnet)
Strongest continuous magnetic field
produced in a laboratory
450,000G
Magnetars Heavy ion collisions The early Universe
(Electroweak transition)
~104MeV2~1017 G ~1022 G ~1013 G
wikipedia
B
Strong magnetic field
Kharzeev, McLerran, Warringa (2008)
√ eB ∼ 100MeV ∼ 1017 − 1018Gauss
Kharzeev, McLerran, Warringa (’08)
Magnetic field in heavy ion collisions
b = 12 fm b = 8 fm b = 4 fm τ(fm) eB (MeV2) 3 2.5 2 1.5 1 0.5 105 104 103 102 101 100
Strong magnetic filed is the QCD scale.
Synchrotron radiation Real photon decays into dileptons
B
e−
e+
B
e− e−
Synchrotron radiation Real photon decays into dileptons
B
e−
e+
B
e− e−
in heavy ion collisions, Tuchin (’11) (’12)
Synchrotron radiation Real photon decays into dileptons
B
e−
e+
B
e− e−
in heavy ion collisions, Tuchin (’11) (’12)
Vacuum birefringence Hattori, Itakura (’12) →Hattori’s talk
Ji
V =
X
f
qfBiNc 2π2 µA
Ji
A =
X
f
qfBiNc 2π2 µ
Symmetry of QCD in a strong constant magnetic field
Internal symmetry
SU(2)L × SU(2)R × U(1)B × U(1)A
Symmetry of QCD in a strong constant magnetic field
Internal symmetry
SU(2)L × SU(2)R × U(1)B × U(1)A U(1)I3,V × U(1)I3,A × U(1)B × U(1)A B
Symmetry of QCD in a strong constant magnetic field
Internal symmetry
SU(2)L × SU(2)R × U(1)B × U(1)A U(1)I3,V × U(1)I3,A × U(1)B × U(1)A B
SSB, anomaly
U(1)I3,V × U(1)B = U(1)em × U(1)B
Symmetry of QCD in a strong constant magnetic field
Internal symmetry
SU(2)L × SU(2)R × U(1)B × U(1)A SO(3, 1) → SO(1, 1)t,z × SO(2)x,y
Lorentz symmetry
U(1)I3,V × U(1)I3,A × U(1)B × U(1)A B
SSB, anomaly
U(1)I3,V × U(1)B = U(1)em × U(1)B
Symmetry of QCD in a strong constant magnetic field
Internal symmetry
SU(2)L × SU(2)R × U(1)B × U(1)A SO(3, 1) → SO(1, 1)t,z × SO(2)x,y
Lorentz symmetry
U(1)I3,V × U(1)I3,A × U(1)B × U(1)A B
SSB, anomaly
U(1)I3,V × U(1)B = U(1)em × U(1)B
Discrete symmetry C, CP , and T are broken.
Charged scalar particle in a magnetic field
Classical equation of motion
H ¨ x = e(E + ˙ x × B),
H = p p2 + m2
Lorentz force
B
Charged scalar particle in a magnetic field
Classical equation of motion
H ¨ x = e(E + ˙ x × B),
H = p p2 + m2
Lorentz force
B
(Landau) quantization Closed orbital motion in the transverse plane
Charged scalar particle
(−D2
µ − m2)φ(x) = 0
Dµ = ∂µ + ieAµ
B = (0, 0, B)
[Dx, Dy] = −ieB
Klein-Gordon equation in a magnetic field
Charged scalar particle
(−D2
µ − m2)φ(x) = 0
Dµ = ∂µ + ieAµ
B = (0, 0, B)
φ(x) = e−iωt+ipzzϕ(x, y) (−D2
x − D2 y)ϕ(x, y) = λϕ(x, y)
λ = ω2 − p2
z − m2
[Dx, Dy] = −ieB
Klein-Gordon equation in a magnetic field
Charged scalar particle
(−D2
µ − m2)φ(x) = 0
Dµ = ∂µ + ieAµ
B = (0, 0, B)
φ(x) = e−iωt+ipzzϕ(x, y) (−D2
x − D2 y)ϕ(x, y) = λϕ(x, y)
λ = ω2 − p2
z − m2
[X, P] = i
Introducing
X = 1 √ eB iDy, P ≡ −1 √ eB iDx,
[Dx, Dy] = −ieB
Klein-Gordon equation in a magnetic field
Charged scalar particle
(−D2
µ − m2)φ(x) = 0
Dµ = ∂µ + ieAµ
B = (0, 0, B)
φ(x) = e−iωt+ipzzϕ(x, y) (−D2
x − D2 y)ϕ(x, y) = λϕ(x, y)
λ = ω2 − p2
z − m2
[X, P] = i
Introducing
X = 1 √ eB iDy, P ≡ −1 √ eB iDx,
[Dx, Dy] = −ieB
Klein-Gordon equation in a magnetic field
looks like a harmonic oscillator.
(−D2
x − D2 y)ϕ(x, y) = eB(X2 + P 2)ϕ(x, y)
Rx = x − iDy eB Ry = y + iDx eB
[D2
x + D2 y, Rx,y] = 0
[Rx, Ry] = − i eB
Degeneracy
N = eB Z dRxdRy 2π = eBV⊥ 2π
Z dxdp 2π~
[x, p] = i~
(cf. ) Magnetic translation
S = eB 2π
(Rx, Ry) corresponds to the center of motion.
a = 1 √ 2(X + iP) a† = 1 √ 2(X − iP) b = r eB 2 (Rx − iRy) b† = r eB 2 (Rx + iRy)
Continuum
Discrete
B
a = 1 √ 2(X + iP) a† = 1 √ 2(X − iP) b = r eB 2 (Rx − iRy) b† = r eB 2 (Rx + iRy)
Continuum
Discrete
B
−D2
x − D2 y = 2eB
✓ a†a + 1 2 ◆
Energy: Wave function: |n, li = (a†)n
p n! (b†)l p l! |0, 0i En = p eB(2n + 1) + p2
z + m2
a = 1 √ 2(X + iP) a† = 1 √ 2(X − iP) b = r eB 2 (Rx − iRy) b† = r eB 2 (Rx + iRy)
Continuum
Discrete
B
Lz = i(xpy − ypx) = b†b − a†a
Angular momentum (symmetric gauge)
−D2
x − D2 y = 2eB
✓ a†a + 1 2 ◆
Energy: Wave function: |n, li = (a†)n
p n! (b†)l p l! |0, 0i En = p eB(2n + 1) + p2
z + m2
σµ = (1, −σi)
iσµDµψL = 0
i ✓ ∂0 − ∂z −Dx + iDy −Dx − iDy ∂0 + ∂z ◆ ψL = ✓ i∂0 − i∂z i √ 2eBa† −i √ 2eBa i∂0 + i∂z ◆ ψL = 0
Weyl Fermion
σµ = (1, −σi)
iσµDµψL = 0
i ✓ ∂0 − ∂z −Dx + iDy −Dx − iDy ∂0 + ∂z ◆ ψL = ✓ i∂0 − i∂z i √ 2eBa† −i √ 2eBa i∂0 + i∂z ◆ ψL = 0
Weyl Fermion
Positive energy solution: En =
p 2eBn + p2
z
|uL(n, pz, l)i = 1 p2En pEn pz
−i √ 2eB √En−pz a
! |n, li
σµ = (1, −σi)
iσµDµψL = 0
i ✓ ∂0 − ∂z −Dx + iDy −Dx − iDy ∂0 + ∂z ◆ ψL = ✓ i∂0 − i∂z i √ 2eBa† −i √ 2eBa i∂0 + i∂z ◆ ψL = 0
Weyl Fermion
LLL has spin up, and down moving
Higher modes can move up and down directions. |uL(0, pz, l)i = θ(pz) ✓1 ◆ |0, li
related to the chiral magnetic effect.
E0 = |pz|
Positive energy solution: En =
p 2eBn + p2
z
|uL(n, pz, l)i = 1 p2En pEn pz
−i √ 2eB √En−pz a
! |n, li
E pz
B = 0
E
· · ·· · ·
n = 1
n = 3 n = 4 n = 2 pz
B 6= 0
E
· · ·· · ·
n = 1
n = 3 n = 4 n = 2 pz
B 6= 0
LLL
E2
Free Dirac particle in a magnetic field
m2
Orbital quantization
E2
Free Dirac particle in a magnetic field
m2
Orbital quantization Zeeman effect (spin 1/2) up down
E2
Free Dirac particle in a magnetic field
m2
Orbital quantization Zeeman effect (spin 1/2) up down
E2
Free Dirac particle in a magnetic field
m2
Orbital quantization Zeeman effect (spin 1/2) up down Lowest Landau Level (LLL)
E2
Free Dirac particle in a magnetic field
m2
Zeeman effect (spin 1) up zero
E2
Free vector particle in a magnetic field
m2
down Lowest Landau Level (LLL)
Scalar boson
becomes heavier.
has instability if
Weyl and Dirac Fermions Vector boson
eB > m2
LLL: zeromode.
E
· · ·· · ·
LLL
n = 1
n = 3 n = 4 n = 2 pz
Chiral magnetic and separation effects
E
· · ·· · ·
LLL
n = 1
n = 3 n = 4 n = 2 pz µ
At finite density, average of current cannot vanish.
Chiral magnetic and separation effects
Chiral magnetic and separation effects
Jz
L ⌘ 1
V Z d3xhJz
L(x)i
Current average
= Nc V⊥ X
n≥1,l
Z dpz 2π pz En (fq(n, l, pz) − f¯
q(n, l, pz))
+Nc V⊥ X
l
Z dpz 2π (−θ(−pz)fq(0, l, pz) + θ(pz)f¯
q(0, l, pz))
Chiral magnetic and separation effects
Jz
L ⌘ 1
V Z d3xhJz
L(x)i
Current average
Higher orders are cancelled out. Only LLL contributes to JL.
= Nc V⊥ X
n≥1,l
Z dpz 2π pz En (fq(n, l, pz) − f¯
q(n, l, pz))
+Nc V⊥ X
l
Z dpz 2π (−θ(−pz)fq(0, l, pz) + θ(pz)f¯
q(0, l, pz))
Chiral magnetic and separation effects
If the system is thermal equilibrium Jz
L =
fq = 1 e(|pz|−µL)/T + 1
f¯
q =
1 e(|pz|+µL)/T + 1
linear in B, independent of T
Nc V⊥ X
l
Z dpz 2π (−θ(−pz)fq(0, l, pz) + θ(pz)f¯
q(0, l, pz))
= −NceB 2π Z ∞ dpz 2π ✓ 1 e(pz−µL)/T + 1 − 1 e(pz−µL)/T + 1 ◆ = −NceB 2π 1 2π T ⇣ ln(1 + eµL/T ) − ln(1 + e−µL/T ) ⌘ = −NceB 2π 1 2π µL
Fukushima, Kharzeev, Warringa (’08)
Son, Zhitnitsky(’04)
Metlitski, Zhitnitsky(’05)
~ JV = ~ JR + ~ JL = Nc 4⇡2 (µR − µL)e ~ B = Nc 2⇡2 µAe ~ B ~ JA = ~ JR − ~ JL = Nc 4⇡2 (µR + µL)e ~ B = Nc 2⇡2 µV e ~ B
nA = jz
V
nV = jz
A
In 1+1d, Using chemical potential,
Chiral magnetic effect
Chiral separation effect
(Jµ
A = −✏µνJV ν)
LLL projected theory
nV = NceB 2π µV π nA = NceB 2π µA π ~ JV = Nc 2⇡2 µAe ~ B ~ JA = Nc 2⇡2 µV e ~ B
Perturbative vs Holographic
Yee(’09), Rubakov(’10), Rebhan, Schmitt, Stricker(’09), Lifshytz, Lippert(’09), Gorsky, Kopnin, Zayakin(’10)
Holographic models
j = σB
Chiral magnetic conductivity
Yee JHEP11 (2009) 08
10 20 30 40 50 0.5 1.0 5 10 15 20 25 30 Ω T 0.5 0.5 1.0 Σ Σ0Sakai-Sugimoto model
Kharzeev, Warringa, Phys. Rev. D80, 034028(2009)Perturbative
µ/T = 0.1
Yamamoto(’11)
Yannis Burnier, Dmitri E. Kharzeev, Jinfeng Liao, Ho-Ung Yee (’11)
Kharzeev, Yee(’10)
~ JV = 1 2⇡2 µAe ~ B ~ JA = 1 2⇡2 µV e ~ B
Yannis Burnier, Dmitri E. Kharzeev, Jinfeng Liao, Ho-Ung Yee (’11)
Kharzeev, Yee(’10)
~ JV = 1 2⇡2 µAe ~ B ~ JA = 1 2⇡2 µV e ~ B ✓ ∂0 ⌥ NceBα 2π2 ∂z DL∂2
z
◆ j0
L,R = 0
gapless collective mode in chiral symmetric phase
0.0 0.5 1.0 1.5 2.0 eB GeV 2 0.0 0.2 0.4 0.6 0.8 1.0 Kharzeev, Yee, Phys. Rev. D83, 085007 (2011)Sakai-Sugimoto model velocity
T = 150MeV T = 200MeV
T = 250MeV
Chiral magnetic waves in Heavy Ion Collisions
±Observed A
0.05 0.1
(%)
2v
3.6 3.65 3.7 3.75
+
π
STAR Preliminary
±A
0.02 0.04
) (%)
+π (
2) - v
(
2v
0.1 0.2
Au + Au 200 GeV: 30-40% < 0.5 GeV/c
T0.15 < p
STAR Preliminary
v2
± = v2 qe
ρe A± A± = N+ − N− N+ + N−
Gapless mode couples to photon Photon becomes massive
1 2 3 4 5 6 1.0 0.5 0.0 0.5 1.0 1.5 2.0 2.5 t(b) electric field (c) quark number density
ω =
+ m2 γ,
Electric field
Particle number
Ez = E0 cos(ωt − k · x),
1.0
0.5
n/nmax
LLL approximation Schwinger mechanism + Maxwell Eq. From anomaly equation
Son, Yamamoto (’12), 1210.8158, Stephanov, Yin (’12)
Action with Berry phase
S = Z dt( ˙ x · p + ˙ x · A(x) − ˙ p · A(p) − H(x, p))
Son, Yamamoto (’12), 1210.8158, Stephanov, Yin (’12)
Action with Berry phase
˙ x = v + ˙ p × Ω
Eom
˙ p = ˙ x × B + F
v = ∂H/∂p F = −∂H/∂x B = r × A Ω = rp × A S = Z dt( ˙ x · p + ˙ x · A(x) − ˙ p · A(p) − H(x, p))
Son, Yamamoto (’12), 1210.8158, Stephanov, Yin (’12)
Action with Berry phase
˙ x = v + ˙ p × Ω
Eom
˙ p = ˙ x × B + F
v = ∂H/∂p F = −∂H/∂x B = r × A Ω = rp × A S = Z dt( ˙ x · p + ˙ x · A(x) − ˙ p · A(p) − H(x, p))
G = (1 + B · Ω)
√ G ˙ x = v + F × Ω + B(v · Ω) √ G ˙ p = F + v × B + Ω(F · B)
Son, Yamamoto (’12), 1210.8158, Stephanov, Yin (’12)
Action with Berry phase
˙ x = v + ˙ p × Ω
Eom
˙ p = ˙ x × B + F
v = ∂H/∂p F = −∂H/∂x B = r × A Ω = rp × A S = Z dt( ˙ x · p + ˙ x · A(x) − ˙ p · A(p) − H(x, p))
G = (1 + B · Ω)
√ G ˙ x = v + F × Ω + B(v · Ω) √ G ˙ p = F + v × B + Ω(F · B)
Chiral magnetic effect
j = Z
p
√ Gf ˙ x = Z fv + F × Z
p
fΩ + B Z
p
f(v · Ω)
Son, Surówka(’09)
Hydrodynamics
Jµ = nuµ + νµ
∂µT µν = F νλJλ
∂µJµ = CEµBµ
T µν = (✏ + P)uµuν + Pgµν + ⌧ µν
Son, Surówka(’09)
Hydrodynamics
Jµ = nuµ + νµ
∂µT µν = F νλJλ
∂µJµ = CEµBµ
T µν = (✏ + P)uµuν + Pgµν + ⌧ µν
ν = −σTP µν∂ν(µ/T) + σEµ + ξBBµ + ξωµ
sµ = suµ − µ T + Dωµ + DBBµ
!µ = 1 2✏µνρσuν@ρuσ
Son, Surówka(’09)
Hydrodynamics
Jµ = nuµ + νµ
∂µT µν = F νλJλ
∂µJµ = CEµBµ
T µν = (✏ + P)uµuν + Pgµν + ⌧ µν
∂µsµ ≥ 0
⇠ = C ✓ µ2 − 2 3 nµ3 ✏ + P ◆
⇠B = C ✓ µ − 1 2 nµ2 ✏ + P ◆
ν = −σTP µν∂ν(µ/T) + σEµ + ξBBµ + ξωµ
sµ = suµ − µ T + Dωµ + DBBµ
!µ = 1 2✏µνρσuν@ρuσ
Linear response theory Perturbative and Holographic
Landsteiner, Megıas, Pena-Benitez (’11)
⇠ = lim
kn→0
X
ij
✏ijn i 2kn hJi
AT 0jiω=0
Amado, Landsteiner, Pena-Benitez (’11)
ξ = µ2 4π2 + T 2 12
related to gravitational anomaly(?)
Linear response theory Perturbative and Holographic
Landsteiner, Megıas, Pena-Benitez (’11)
⇠ = lim
kn→0
X
ij
✏ijn i 2kn hJi
AT 0jiω=0
Amado, Landsteiner, Pena-Benitez (’11)
ξ = µ2 4π2 + T 2 12
related to gravitational anomaly(?)
Linear response theory Perturbative and Holographic
Landsteiner, Megıas, Pena-Benitez (’11)
⇠ = lim
kn→0
X
ij
✏ijn i 2kn hJi
AT 0jiω=0
Amado, Landsteiner, Pena-Benitez (’11)
ξ = µ2 4π2 + T 2 12
Golkar, Son (’12)
= magnetic (separation) effects
= magnetic (separation) effects
Lorentz force Coriolis force
Generating functional for Fluid
Banerjee, Bhattacharya, Bhattacharyya, Jain, Minwalla, Sharma (’12)
Jensen, Kaminski, Kovtun, Meyer, Ritz, Yarom (’12)
Gravitational and gauge background with a timelike killing vector ds2 = −e2σ(dt + aidxi)2 + gijdxidxj
A = A0dx0 + Aidxi
T = e−σT0 + · · · µ = e−σ(µ0 + A) + · · ·
Local temperature Local chemical potential
Wanom = C 2 ✓Z A0 3T0 AdA + A2 6T0 Ada ◆
Next leading order (First derivative)
Leading order
W0 = Z d3x√g3 eσ T0 P(T0e−σ, e−σA0) T µν = −2T0 δW δgµν
Jµ = T0 δW δAµ
= nuµ = (✏ + P) + Pgµν