Hidenori Fukaya (Osaka U.)
for JLQCD Collaboration
JLQCD collaboration, arXiv:1111.0417
S. Aoki (Tsukuba) & HF, PRD 84, 014501 (2011) [arXiv:1105.1606][P.H. Damgaard (NBI,NBIA) & HF, JHEP0901:052, 2009]
新学術領域「素核宇宙融合による計算科学に基づいた重層的物質構造の解明」
公募研究 23105710「カイラル対称性の破れにともなう物理の定量評価」

1

1. Introduction
Q.
What is the “epsilon regime” ?
A.
Quantum Chromo Dynamics (QCD)in (vicinity of) thelimit
in a finite volume.
2

1. Introduction [Nambu, 1961 ] (2008 Nobel Prize)
Chiral symmetry is important near
• Chiral symmetry breaking and constituent mass
→ “effective” quark action
acquires constituent mass
→ hadron masses are ~ O(1) GeV.

• Pion effective theory（pseudo) Nambu-Goldstone boson = pion
described by Chiral perturbation theory （ChPT)[Weinberg 1979]3

1. Introduction	
Origin of mass = chiral symmetry breaking.
 Hadrons consist of quarks.

uud

 But

proton mass >> quark mass ×3(1GeV)(3-6MeV) Chiral symmetry breaking generate~90% of mass. 	
4

1. Introduction	JLQCD (+ TWQCD) collaboration
have been simulating lattice QCD
with exact chiral symmetry
using overlap fermion action.

5

1. Introduction	JLQCD (+ TWQCD) collaboration
have been simulating lattice QCD
with exact chiral symmetry
using overlap fermion action.

6

1. Introduction	
Chiral symmetry is expensive.
differenct action = different errors and cost.
Fermion action	

ChiralDiscretizat
symmetry	 ion error	

Numerical cost	

Lattice size	Overlap

Exact	O(a2)	

Very expensive	< 2.4 fm	Domain-wall	

Weaklybroken	O(a2)	

Expensive	~4fm	Wilson 	Broken	O(a)	

Marginal	~5fm	Staggered	Broken

O(a2)	

Cheap	~6fm	

(U(1) remains)	

JLQCD = Small volume QCD…	

7

1. Introduction
Finite volume = Pion physics.
Correlation length (1/M) of QCD particles
Pions(~140MeV) ~ 1.4fm
Kaons(~500MeV)~ 0.4fm
Rho (~800MeV)~0.26fm
Proton (~1GeV) ~0.2fm
At
QCD = pion (+ kaon) theory.
Finite volume correction in QCD =chiral perturbation theory (ChPT)
weakly coupled = analytically calculable.

81. Introduction
Our strategy
Numerical calculation

Analytic calculation

Small lattice QCD
with exact chiral symmetry	

Finite volumecorrection by pion
theory	Chiral symmetry isimportant but expensive.

QCD at V=∞	9Contents✓ 1. Introduction

2. Analytic calculation of finite V effects(on the pions)
3. Numerical lattice QCD results
4. Summary102. Analytic calculation of finite V
effects
Pion correlator at V=(In Euclidean space-time,)

11

2. Analytic calculation of finite Veffects
Pion correlator at finite V ( in the p-expansion )(periodic boundary for t-direction)

12

2. Analytic calculation of finite VeffectsBUT…

In the limit,Infra-red divergence due to finite V ???
despite we have IR cut-off 1/V1/4?

Something wrong ! Exp->Cosh is not enough !
13

2. Analytic calculation of finite Veffects
Many vacua contributes at finite V
This fake IR divergence is due to a fixed vacuum:

but at finite V, the vacuum is not uniquelydetermined: vacuum= moduli = dynamical variable

U0 should be non-perturbatively treated.
→ ε-expansion (is needed for)14

2. Analytic calculation of finite Veffects
ε and p expansions are really useful ?
ε expansion

p expansion

Group	Nf	Action	

a(fm)	L	Mπ (MeV)	

ETMC

2	

Twisted mass	

0.05-0.100	

~3fm	

280	

MILC	

2+1	

Staggered	

0.045-0.12	

3~6 fm	250	

RBC/UKQCD	

2+1	

Domain wall	

0.085-0.11	

3~4fm	290	JLQCD	

2+1	Overlap	

0.11	

1.8fm	

310	PACS-CS	

2+1	

Wilson	0.09	~3fm	

140	

BMW	

2+1	

Wilson	0.065-0.125	

3~5fm	

190	

ALV	

2+1	

DW on MILC	

0.06-0.12	

3~4fm	

250	HPQCD	

2+1	

HISQ	

0.045-0.15	

3~4fm	360	But on the lattice,	

No requirement in original theory ( if):both expansions are bad. → Better way of expansion ?
15

2. Analytic calculation of finite Veffects
New expansion
p-expansion
ε-expansion

16

2. Analytic calculation of finite Veffects
New expansion
p-expansion
ε-expansion

New i (interpolating)- expansion

17[Damgaard & HF, 2008]	

2. Analytic calculation of finite Veffects
The pion (chiral) Lagrangian at finite V
Zero-mode = SU(N) matrix model	
Non-zero-mode = massive bosons	(perturbative) interactions	= a hybrid system of matrix model and bosonic fields
18

2. Analytic calculation of finite Veffects
Pion correlator at 1-loop[Aoki & HF, 2011]	

Because of the mixing of zero and non-zero modes,
the calculation is fairly tedious :19

2. Analytic calculation of finite Veffects
where[Aoki & HF, 2011]	

20

2. Analytic calculation of finite Veffects
Pion correlator at 1-loop[Aoki & HF, 2011]	

cancels IR divergence:
disappears in the p-regime:(functions of)	21

3. Numerical lattice QCD results(preliminary)
QCD simulation with exact chiral symmetry[JLQCD & TWQCD collaborations, 2006-2011]	

2+1-flavor overlap Dirac fermions [Neuberger 98]
Iwasaki gauge action, β=2.3, 1/a ~ 1.759 GeV.
Lattice size : L=16 [1.8 fm], T=48.
Topology fixed: Q=0 (or 1)
Quark masses : ms=0.08, 0.100,
mud= 0.002, 0.015,0.025,0.035,0.050,0.080,0.100(~3 MeV)ε-regime,(30MeV <)p-regime	

22

3. Numerical lattice QCD results(preliminary)Really inside the ε-regime?→ Scalar channel knows.

mud=0.002	

mud=0.015	

23

3. Numerical lattice QCD results(preliminary)Really inside the ε-regime? → Yes.
ε expansion analysis :mud=0.002	

[Damgaard et al, 2002]	

PP →SS →consistent with our previous value
from Dirac spectrum.

Next, let us extract

with the “i” expansion formula.
( Naïve GMOR relation suggests)	 24

3. Numerical lattice QCD results(preliminary)
Pion mass and decay constant in the ε regime
Fitting withafter 1-loop (of non-zeromodes) volume
correction, we obtain25

3. Numerical lattice QCD results(preliminary)
Chiral “interpolation” forbigger than ourprevious analysis using
p-regime data only :Linear fit looks better
than NLO ChPT fit, though…
26

3. Numerical lattice QCD results(preliminary)
Chiral “interpolation” for
ChPT fit looks bad :deviation.
2-loop (1/V) effects fromnon-zero modes undercontrol ?

~7.5%	

27

3. Numerical lattice QCD results(preliminary)
Bad convergence of ChPT for
In the limit?

,Bad.
(big correction/LO)
Good.

28

4. Summary
Small lattice QCD with exact chiral symmetry+
Finite V correction from pion effective theory(ChPT)

can extract physics at V=∞.
29

A. A short-cut prescriptionAny interpolation is a good interpolation.
[Aoki & HF, 2011]	

1. The epsilon-expansion and p-expansion aredifferent expansions of the same theory.
2. They should converge in the limit V→∞.
3. Any interpolation gives the same results in thatlimit.
4. No additional information is needed : all thenecessary ingredients can be read off from the εand p expansions.

30

A. A short-cut prescription
A quick recipe[Aoki & HF, 2011]	

Starting from the p-regime result,
1. Remove IR divergent part from the propagator:

2. Multiply a factor from zero-mode integrals, which
can be read off from the ε-regime results.
3. Add a constant term if it exists in the ε-expansion.
31

A. A short-cut prescription
Results[Aoki & HF, 2011]	

(The same result as our full calculation !)	

Axial-Ward-Takahashi identities also confirmed.

32

B. Comparison with Colangelo etal. [2005]
Differences from p-regime calculation:“Finite volume effects for meson masses and decay
constants,” G. Colangelo, S.Durr, C. Haefeli NPB 721[2005]
Loop
corrections	

Chiral limit	

Temporaldirection	

Colangelo et al.
2005	2-loop	

fake IR
divergence	neglected	

Our work	

1-loop	finite	finite	33

B. Comparison with Colangelo etal. [2005]
Pion mass:

34

B. Comparison with Colangelo etal. [2005]
Pion decay constant:35

3. Numerical lattice QCD results(preliminary)
How to control convergence of ChPT
p-expansion needs
keeping
New i-expansion:
We can separately take2-loop corrections ->Physical point simulation on L~3 fm lattice isinteresting on next generation machines.	

36

2. New chiral expansion
Mass term plays a key role.
Mass term decomposition into 3 pieces :Third term is NLO in both of p and ε regimes.
It is natural to assume it is perturbative everywhere:
37[Damgaard & HF, 2008]