Transcript QCD Vacuum

Role of Anderson-Mott localization in the
QCD phase transitions
Antonio M. García-García
[email protected]
Princeton University
ICTP, Trieste
We investigate in what situations Anderson localization may be relevant in the
context of QCD. At the chiral phase transition we provide compelling evidence
from lattice and phenomenological instanton liquid models that the QCD Dirac
operator undergoes a metal - insulator transition similar to the one observed in
a disordered conductor. This suggests that Anderson localization plays a
fundamental role in the chiral phase transition. Based on a recent relation
between the Polyakov loop and the spectral properties of the Dirac operator we
discuss how the confinement-deconfinement transition may be related to a
metal-insulator transition in the bulk of the spectrum of the Dirac operator.
James Osborn
In collaboration with
PRD,75 (2007) 034503 ,NPA, 770, 141 (2006) PRL 93 (2004) 132002
Outline
1. A few words about localization and spectral analysis.
2. Disorder in QCD, random matrix theory.
Dyakonov - Petrov ideas.
3. A few words about QCD phase transitions.
4. Role of localization in the QCD phase transitions. Results
from instanton liquid models and lattice.
4.1 The chiral phase transition.
4.2 The deconfinement transition. In progress.
5. What’s next. Relation confinement and chiral symmetry
breaking. Quark diffusion in LHC.
A few words on localization
Quantum particle in a random potential
Anderson-Mott localization
Quantum destructive interference, tunneling
or interactions can induce a transition
to an insulating state.
Insulator
For d < 3 or, at strong disorder,
in d > 3 all eigenstates are
localized in space.
Classical diffusion is
eventually stopped
Metal
d > 2, Weak disorder
Eigenstates delocalized.
Quantum effects do not alter
significantly the classical diffusion.
Insulator
Metal
How do we know that a metal is a metal?
Texbook answer: Look at the conductivity or other transport properties
Other options: Look at eigenvalue and eigenvectors

 D
QCD
 n  in  n
H n  En n
1. Eigenfunctions moments:
1 Insulator
IPR    ( r ) d r ~ 
Metal
V
4
d
1
n
2. Decay of the eigenfunctions:

e


 ( r ) ~ 1 / V
 1/ r 

r/
n
Insulator
 d
Metal

  d
?
 d

Insulator
Critical
Metal
3. Spectral characterization:
Spectral correlations in a metal are given by random
matrix theory up to the Thouless energy. Matrix
elements are only constrained by symmetry
2
2
Metal 
n ( E ) = n  n ~ log E

β  As 2
( RMT ) 


P
s
~
s
e

2
In units of the mean level spacing,
the Thouless energy,
d 2
Ec  L
n ( E )  E
2
d /2
Eigenvalues in an insulator are not correlated.
Insulator n ( E )  E

( Poisson )  P( s)  e s
2
In the context of QCD the metallic
region corresponds with the infrared
limit (constant fields) of the Dirac
operator" (Verbaarschot,Shuryak)
Random Matrix Theory: Basic Facts
Idea: Spectrum of a RM is highly correlated... not random


Hij Random (real/complex Gaussian number)
PH  = A exp  TrH 2 dH
ρ x1  x N    xi  x j 
i< j
β
N

2
exp   β  xi 
i


Spectral correlators
β=2
Rn x1, .., xn  ~ detK xi , x j 
N 
Unfolding
sin π  x  y 
K x, y  =
π x  y 
Signatures of RM spectrum
Spectral rigidity
Σ 2 n  = n
Level Repulsion
n n
2
 n
2
~   R2 x, y dxdy ~ log n
o 0
n  1 Ps  ~ s e
β
 As 2
QCD vacuum and instantons
Diakonov, Petrov, later Verbaarschot, Shuryak, Osborn, Zahed, Osborn & AGG
Dirac operator has a zero mode in the field of an instanton

Dψ0 r   0 ψ0 r  1/ r
ins 
μ
D =  μ + gA 
QCD vacuum saturated by weakly interacting (anti) instantons
3
(Shuryak)
Density and size of instantons are fixed phenomenologically
(


)
T   d x ( x  z )iD  ( x  z ) ~ i(u  Rˆ )
R
4
IA


I
I
I
A
A
A
4
Long range hopping in the instanton liquid model (ILM)
Diakonov - Petrov
As a consequence of the long range hopping the QCD vacuum is a metal:
Zero modes initially bounded to an instanton get delocalized due to the
overlapping with the rest of zero modes. By increasing temperature (or
other parameters) the QCD vacuum will eventually undergo a metal
insulator transition.
3
QCD vacuum as a disordered medium
METAL
Diakonov, Petrov, later Verbaarschot
Osborn, Shuryak, Zahed
An electron initially bounded to a single atom gets delocalized due to
the overlapping with nearest neighbors.
QCD Vacuum
Zero modes initially bounded to an instanton get delocalized due to
the overlapping with the rest of zero modes.
Impurities
Instantons
Electron
Quarks
Differences
Dis.Sys: Exponential decay
QCD vacuum Power law decay
Nearest neighbors
Long range hopping!
Conductivity versus chiral symmetry breaking
Spectral properties of the smallest eigenvalues of the Dirac operator are
controled by instantons" Is that important? Yes.
    lim 
m  m0
 (m)
Banks-Casher (Kubo)
V
1 3N N 
   
  (240MeV )
  2 V 
1/ 2
Metallic behavior means chiSB in the ILM
c
3
Recent developments:
1. RMT and Thouless energy in QCD. Metallic limit of QCD.
Verbaarschot,Osborn, PRL 81 (1998) 268, Zahed, Janik et.al., PRL. 81 (1998) 264, Verbaarschot, Shuryak
Nuclear Physics A 560 306 ,1993
2.The QCD Dirac operator can be described by a random matrix with
long range hopping even beyond the Thouless energy.
AGG and Osborn, PRL, 94 (2005) 244102
3. Localization in the QCD Dirac operator:
C. Gattringer, M. Gockeler, et.al. Nucl. Phys.
B618, 205 (2001),R.V. Gavai, S. Gupta and R. Lacaze, Phys. Rev. D 65, 094504 (2002),M. Golterman and Y.
Shamir, Phys. Rev. D 68, 074501 (2003);V. Weinberg, E.-M. Ilgenfritz, et.al, PoS { LAT2005}, 171 (2005);F.
Farchioni, Ph. de Forcrand, I. Hip, C. B. Lang and K. Splittorff,Phys. Rev. D 62, 014503 (2000);P.H. Damgaard,
U.M. Heller, et.al.,Nucl. Phys. B583, 347 (2000);M. Gockeler, P.E.L. Rakow,et.al,Phys. Rev. Lett. 87, 042001
(2001);J. Greensite, A.V. Kovalenko, et.al Phys. Rev. D 74,094507 (2006),
1. QCD, random matrix theory,
Thouless energy:
Spectral correlations of the QCD Dirac operator in the infrared limit are
universal (Verbaarschot, Shuryak Nuclear Physics A 560 306 ,1993). They can be
obtained from a RMT with the symmetries of QCD.
1. The microscopic spectral density is universal, it depends only on the
global symmetries of QCD, and can be computed from random matrix
theory.
 (s) ~ sJ
2
0
( s)  J ( s)
2
1

2. RMT describes the eigenvalue correlations of the full QCD Dirac
operator up to Ec. This is a finite size effect. In the thermodynamic
limit the spectral window in which RMT applies vanishes but at the
same time the number of eigenvalues, g, described by RMT diverges.
Fπ2
Ec ~ 2  0
L
L
Ec
g  ~ Fπ2 L2  

L
2. Random matrix approach
beyond the Thouless energy:
D
QCD
 0 W

  
W 0 
WIA ~ aIA
The spectral correlations of the QCD Dirac operator agree with
those of a (chiral) random matrix with the appropiate symmetry.
(Verbaarschot,Shuryak)
W ~a /R
IA
IA
3
IA
We describe the spectral correlations of the QCD Dirac operator
analytically well beyond this limit by combining the random
matrix approach with the phenomenological power-law decay of
matrix elements from ILM. Agreement with chiral perturbation
theory, Phys.Rev.Lett. 93 (2004) 132002
 
2
2
Phase transitions in QCD
At which temperature does the transition occur ? What is the nature of transition ?
Péter Petreczky
J. Phys. G30 (2004) S1259
Quark- Gluon Plasma
perturbation theory only for
T>>Tc
Deconfinement and chiral restoration
They must be related but nobody* knows exactly how
Deconfinement
•Linear confining
potential vanishes.
L
 0
Chiral Restoration
•Matter becomes light
•QCD still non
 ~ 0
perturbative
How to explain these transitions?
1. Effective model of QCD close to the chiral restoration (Wilczek,Pisarski):
Universality, epsilon expansion.... too simple?
2. QCD but only consider certain classical solutions (t'Hooft):
Instantons (chiral), Monopoles (confinement)
No monopoles found, instantons only after lattice cooling, no from QCD
We propose that quantum interference and tunneling, namely, Anderson
localization plays an important role. Nuclear Physics A, 770, 141 (2006)
Localization and chiral transition: Why do we expect
a metal insulator transition close to the origin at finite
temperature?
1.The effective QCD coupling constant g(T) decreases as
temperature increases. The density of instantons also
decreases  (T )   (0) / 2
c
2. Zero modes are exponentially localized  ( R)  exp(TR )
in space but oscillatory in time.
3. Amplitude TIA ~ exp( ATR) hopping restricted to
neighboring instantons.
4. Localization will depend strongly on the temperature. There
must exist a T = TLsuch that a MIT takes place.
Dyakonov,
5. There must exist a T = Tc such that
   0 Petrov
6. This general picture is valid beyond the instanton liquid
approximation (KvBLL solutions) provided that the hopping
induced by topological objects is short range.
Is TL = Tc ?...Yes
Does the MIT occur at the origin? Yes
Main Result
D =  + gA
QCD
μ
 D

QCD

μ
  i 
n
At Tc , Chiral phase transition
A  A , A
lat
n
ins

  0
but also the low lying,



n
undergo a metal-insulator transition.
n
"A metal-insulator transition in the Dirac operator induces the
chiral phase transition "
Signatures of a metal-insulator transition
1. Scale invariance of the spectral correlations.
A finite size scaling analysis is then carried
out to determine the transition point.
Skolovski, Shapiro, Altshuler
2.
P( s) ~ s 
P(s) ~ e  As
s  1
s  1
 (n) ~ n
2
  (r ) d r ~ L
2q
3. Eigenstates are multifractals.
d
 Dq ( q 1)
n
Mobility edge
Anderson transition
1. Finite size
scaling analysis:
Dynamical 2+1
Unquenched
var  s  s
2
2
s   s P(s)ds
n
Quenched
n
Quenched Lattice
IPR versus eigenvalue
2. Level repulsion s << 1
Exponential decay s > 1
Unquenched ILM, 2 m = 0
3. Transiton from
metal to insulator
Unquenched lattice, close to the origin,
2+1 flavors, 163X6
INSULATOR
METAL
Unquenched ILM, close to the origin,
2+1 flavors, N = 200
4. Localization versus
chiral transition
Instanton liquid model Nf=2, masless
Lattice: and inverse participation ratio versus T
Unquenched, 2+1
Quenched
For zero mass, transition sharper with the volume
First order
For finite mass, the condensate is volume independent
Crossover
Localization and order of the chiral phase transition
    lim 
m  m0
 (m)
V
1. Metal insulator transition always occur close to the origin.
2. Systems with chiral symmetry the spectral density is sensitive to
localization.
3. For zero mass localization predicts a (first) phase transition not
crossover.
4. For a non zero mass m, eigenvalues up to m contribute to the
condensate but the metal insulator transition occurs close to the
origin only. Larger eigenvalue are delocalized se we expect a
crossover.
5. Multifractal dimension m=0 should modify susceptibility
exponents.
Confinement and spectral properties
Idea: Polyakov loop is expressed as the response of the Dirac operator to a
change in time boundary conditions Gattringer,PRL 97 (2006) 032003, hep-lat/0612020


U 4 ( x, N )  zU 4 ( x, N )


 2 N  ( x )  (1  z1 ) Nz1  ( x ) z1 

 z1

1  

L( x ) 

N
8 N   (1  z 2 )  z  ( x ) z

2
2



z
2


Politely Challenged in:
heplat/0703018,
Synatschke, Wipf,
Wozar


1 
N
N
N 
P  L( x ) 
2   (1  z1 ) z1  (1  z2 ) z2

8V  
z1
 z2






 ( x )   v  ( x, t )  v  ( x, t )
N
t 1
L,
R,
…. but sensitivity to spatial boundary conditions
is a criterium (Thouless) for localization!
Localization and confinement
The dimensionless conductance, g, a measure of localization, is related to
the sensitivity of eigenstates to a change in boundary conditions.
Metal
g  Ld 2
Insulator
g 0 L 
MI transition
L
g  g c (d ) L  
1.What part of the spectrum contributes the most to the
Polyakov loop?.Does it scale with volume?
2. Does it depend on temperature?
3. Is this region related to a metal-insulator transition at Tc?
4. What is the estimation of the P from localization theory?
Accumulated Polyakov loop versus eigenvalue
Confinement is controlled by the ultraviolet part of the spectrum
P

Localization and Confinement
IPR (red), Accumulated Polyakov loop (blue) for T>Tc as a
function of the eigenvalue.
Metal
prediction
MI
transition?
Quenched ILM, IPR, N = 2000
Metal
IPR X N= 1
Insulator
IPR X N = N
Multifractal
Similar to overlap prediction
Origin
Morozov,Ilgenfritz,Weinberg, et.al.
Bulk
IPR X N =
D2~2.3(origin)
N
D2
Quenched ILM, T =200, bulk
Mobility edge in the Dirac
operator. For T =200 the
transition occurs around the
center of the spectrum
D2~1.5 similar to the
3D Anderson model.
Not related to chiral
symmetry
Unquenched ILM, 2+1 flavors
We have observed a metal-insulator transition at T ~ 125 Mev
Conclusions
●
●
●
Eigenvectors of the QCD Dirac operator becomes more
localized as the temperature is increased.
For a specific temperature we have observed a metalinsulator transition in the QCD Dirac operator.
For lattice and ILM, and for quenched and unquenched we
have found two transitions close to the origin and in the UV
part of the spectrum and.
MAIN
"The Anderson transition occurs at the same T than the chiral
phase transition and in the same spectral region"
“ Confinement-Deconfinemente transition has to do with
localization-delocalization in time direction”
What's next?
1. How critical exponents are affected by
localization?
2. Confinement and localization, analytical
result?
3. How are transport coefficients in the quark
gluon plasma affected by localization?
4. Localization in finite density. Color
superconductivity.
QCD : The Theory of the strong interactions
L   (i   gA )  m    G


QCD
1
4
q
2
High Energy g << 1 Perturbative
1. Asymptotic freedom
Quark+gluons, Well understood
Low Energy
g ~ 1 Lattice simulations
The world around us
2. Chiral symmetry breaking
 ~ (240MeV )
3
Massive constituent quark
3. Confinement
Colorless hadrons
V (r )  a / r  r
Analytical information? Instantons , Monopoles, Vortices
Quenched ILM, Origin, N = 2000
For T < 100 MeV we expect (finite size scaling) to
see a (slow) convergence to RMT results.
T = 100-140, the metal insulator transition occurs
IPR, two massless flavors
D2 ~ 1.5 (bulk) D2~2.3(origin)
A A
 W  =
A A
RMT
P
RMT
A =  s Ps ds

2
0
How to get information from a bunch of levels
Spectrum
Unfolding
Spectral Correlators
Quenched ILM, Bulk, T=200
Nuclear (quark) matter at finite temperature
1. Cosmology 10-6 sec after Bing Bang, neutron stars (astro)
3 Analytical, 4N=4 super YM ?
2. 1
Lattice QCD finite2size effects.
3. High energy Heavy Ion Collisions. RHIC, LHC
Colliding Nuclei
Hard
Collisions
QG Plasma ?
Hadron Gas &
Freeze-out
sNN = 130, 200 GeV
(center-of-mass energy per nucleon-nucleon collision)
Multifractality
Intuitive: Points in which the modulus of the
wave function is bigger than a (small) cutoff M. If
the fractal dimension depends on the cutoff M,
Kravtsov, Chalker,Aoki,
the wave function is multifractal.
Schreiber,Castellani
IPR  I =  ψ r  d r  L
4
2
Ld
n
d
 D2
Instanton liquid models T = 0
"QCD vacuum saturated by interacting (anti) instantons"
Density and size of (a)instantons are fixed phenomenologically
The Dirac operator D, in a basis of single I,A:
0
iD  
T
T 

0

IA
AI
  200MeV ,
N
 1 fm
V
4
(


)
T   d x ( x  z )iD  ( x  z ) ~ i(u  Rˆ )
R
4
IA


I
I
I
A
A
3
A
4
1. ILM explains the chiSB
2. Describe non perturbative effects in hadronic correlation functions
(Shuryak,Schaefer,dyakonov,petrov,verbaarchot)
QCD Chiral Symmetries
L,R  (1   5 )
Classical
SU A (3)  SUV (3)  UV (1)  U A (1)

Quantum
SUV (3)  UV (1)
U(1)A explicitly broken by the anomaly.
SU(3)A spontaneously broken by the QCD vacuum
qq  (250 MeV )
3
Dynamical mass
Eight light Bosons (,K,), no parity doublets.
Quenched lattice QCD simulations
Symanzik 1-loop glue with asqtad valence