3. Quantum Chromodynamis on the Lattie 15
3.3. Twisted Mass F ermions
The modiation of Wilson type fermions by a hirally twisted mass term has been
introdued about a deade ago [81, 82℄. Some years later Frezzotti and Rossi realised
that thetheoryis automatially
O (a)
-improved iftheoriginal untwisted massisset toits ritialvalue[83 ℄. Ourpresentation oftwisted masslattiefermions isbasedon the
reviewbyShindler[84 ℄.
3.3.1. Twisted Mass Formulation
The entral keyto the introdution of twisted massfermions isthetransformation
ψ =
eiωγ 5 τ 3 /2 χ
andψ = χ
eiωγ 5 τ 3 /2
(3.20)whihisasymmetryoftheontinuumationifthemassterm isreinterpreted
appropri-ately.
ψ
,inthe so-alledphysialbasis,andχ
,inthetwistedbasis,areavour doubletsand thethirdPauli matrix
τ 3
ats intheaording avour spae. It isthenpossibletoreformulate the QCDation as
S F = Z
d
4 x χ(x) γ µ D µ + m 0 +
iµ 0 γ 5 τ 3
χ(x) .
(3.21)The barequarkmassisrelated to the untwisted andtwisted omponents,
µ 0 = m q sin ω ,
(3.22a)m 0 = m q cos ω ,
(3.22b)m q = q
m 2 0 + µ 2 0 ,
(3.22)withthe orrespondingtwistangle
ω
,tan ω = µ 0
m 0 .
(3.23)Thetwist angle obviously hanges underrenormalisation astheuntwisted and twisted
massomponentsdonotsharethesamerenormalisationfator. Therenormalisedquark
massis determined as
m R = q
Z m 2 (m 0 − m c ) 2 + Z µ 2 µ 2 0 .
(3.24)Only if the untwisted quark mass is set to its ritial value, the angle,
ω = π/2
, isindependent ofrenormalisation. This aseisalledmaximal twist.
Applyingthe Wilson disretisation to theontinuumation withtwisted massterm,
weobtain
S F [U, ψ, ψ] = X
x
χ(x) 1 − κD W [U] + 2
iκaµ 0 γ 5 τ 3
χ(x) .
(3.25)Apartiular feature ofthetwisted massfermionmatrix,
M
tm= 1 − κD W [U] + 2
iκaµ 0 γ 5 τ 3 ,
(3.26)thatwasimportant fortheoriginalmotivationtointroduethis modiationofWilson
fermionsis that the twisted massparameter ats as aregulator proteting simulations
fromzeromodes. Thesezeromodeshavebeenasoure ofonernintheuntwisted
the-ory[82 ℄. Thispropertyan beseenwhen evaluatingtheavour determinant expliitly,
Det
avour
M
tm= M W M W † + (aµ 0 ) 2 ,
(3.27)where
M W
denotes the standard Wilson fermionmatrix for one avour.Inversionof equation(3.26) inmomentumspae for afreeeldleads to the
determi-nationof the propagator,
S
tm(p) = a −
iP
ν γ ν p ν + 1 2 p ˆ 2 + am 0 −
iaµ 0 γ 5 τ 3 p 2 + 1 2 p ˆ 2 + am 0 2
+ (aµ 0 ) 2 ,
(3.28)whih isimportant for the later analytial alulations in thefree and weakly oupled
limits.
3.3.2. Automati Improvement
Sinethe twistrotations(3.20) arenomoreasymmetryofthelattieationdue tothe
Wilsonterm, thefermiondisretisation providesdierent regularisations dependingon
thetwistangle. Thatmeansinpartiularthattheutoeetsarefuntionsofthetwist.
Ithasbeen shownthattheleading
O (a)
eetsvanishinmostrelevantquantities iftheationis tunedto maximaltwist. Morepreisely,this statement appliesto parityeven
quantities as an be derived from the symmetries of the Symanzikexpansion [84, 85℄.
Thisfeature is known as automati
O (a)
-improvement and is the ruial property for urrent appliationsof twisted massfermions.We now briey sketh the main idea for the proof of automati
O (a)
improvement following the reviewbyShindler[84℄. ThestartingpointistheSymanzikexpansionforthelattie ation,
S
e= S 0 + aS 1 + . . . ,
(3.29)where
S 0
isthe ontinuum ation. Thehigher order ontributions,S n = Z
d
4 y L n ,
(3.30)are onstruted from all possible ounterterms. Espeially, we have innext-to-leading
order
L 1 =
5
X
k=1
c k O k .
(3.31)We do not give the expliit form of the operators
O k
here but refer to [84 ℄ where theomplete set is given. The oeients
c k
have then to be tuned suh that ontinuumalulationswith
S
ereprodue thelattietheoryuptothehosenorderofa
. Asimilarexpansion holds foranyeldobservable,
φ
e= φ 0 + aφ 1 + . . . .
(3.32)With theabove expansions, one an identify the possible disretisation eetsto
O (a)
for an arbitraryexpetation value,
h φ(x 1 ) . . . φ(x n ) i = 1 Z
Z
D [ψ, ψ, A µ ] φ(x 1 ) . . . φ(x n )
e−S 0 (1 − aS 1 + . . .)
= h φ 0 (x 1 ) . . . φ 0 (x n ) i 0
− a Z
d
4 y h φ 0 (x 1 ) . . . φ 0 (x n ) L 1 (y) i 0
(3.33)+ a
n
X
k=1
h φ 0 (x 1 ) . . . φ 1 (x k ) . . . φ 0 (x n ) i 0 + O (a 2 ) ,
where ontat termsan be absorbed into a redenition of theelds. For theproof of
automati improvement one anutilise thefollowing symmetryof thelattieation,
R 5 1 × D × [µ → − µ] ,
(3.34)where
D
multiplies all terms by( − 1) d m
aording to their mass dimensiond m
.R 1 5
transforms theelds as
χ →
iγ 5 τ 1 χ ,
(3.35a)χ → χ
iγ 5 τ 1 .
(3.35b)This symmetryholds for the lattie ationeven inthe innitevolume limit. Therefore
all termson the right handside of equation(3.33) have to respetthe symmetry. This
neessitates the absene of
O (a)
termsfor anyeldφ
thatrespetsR 1 5
sine thenbothontributions to thatorder in equation(3.33) areodd undertheabove symmetry. For
the rst part,
φ 0 (x 1 ) . . . φ 0 (x n ) L 1 (y)
, this is true beause of the operator insertions inL 1
. Theseond part,φ 0 (x 1 ) . . . φ 1 (x k ) . . . φ 0 (x n )
is odd beause of themassdimensionof
φ 1
thathas to be dierent by 1 fromthat ofφ 0
so thatall terms inequation (3.32)possess the same symmetries. Introduing an untwisted quark massto the lattie
a-tion expliitly breaksthe ruial symmetry (3.34) and therefore there is no automati
improvement exept for maximal twist. Note that aording to the above arguments
alsoquantities withvanishingexpetationvaluean have disretisationeets in
O (a)
.Apossible aveatisposedbypionpoles,i.e.byutoeetsthatareof
O (a 2k /m 2h π )
with
2k ≥ h ≥ 1
. Thesepolesan beavoidedbya suitabledenitionof maximaltwist.For this reason ETMC have adopted the ondition of a vanishing PCAC (partially
onserved axialurrent)quarkmass[86℄,
m
PCAC= P
x h ∂ 4 A a 4 (x, t)P a (0) i 2 P
x h P a (x, t)P a (0) i ,
(3.36)where
P a = ψγ 5 τ a /2ψ
,A a µ = ψγ µ γ 5 τ a /2ψ
. From thenumerial sidea quenhed test ofa 2
-saling withthis onditionhasbeen presented in[87 ,88, 89℄.On the level of free lattie fermions the improvement an be observed easily, for
instane, for the dispersion relation that is obtained from the pole of the propagator,
S F (p, p 4 =
iE)
,E(p) = q
p 2 + m 2 q − a m 3 q 2
cos(ω) q
p 2 + m 2 q
+ O (a 2 ) .
(3.37)The improvement for maximal twist,
ω = π/2
is apparent. Note that for a theorywithoutspontaneousbreakingofhiral symmetry,like inthefreelimit,thelinearuto
eets have to be proportional to the quark mass. This is not neessarily the ase in
thepresene ofspontaneous hiral symmetrybreaking [85℄.
Automati improvement also holds at nite temperature. Thisis expeted sine the
relevant symmetries do not mixspatial and temporaldiretions sothat thenite time
extent whihharaterises simulations atnite temperature doesnotinterferewiththe
improvement. A report of observed
a 2
-saling at nite temperature in the quenhed ase has been given by our tmfT ollaborators [90℄. Our perturbative analysis of thepressure,disussed insetion 4.2,nds improvement at maximaltwist aswell.
Another possible drawbak of twisted mass fermions is thebreaking of avour
sym-metry. Flavour singlet and doublet quantities suh as
m 0 π
andm ± π
may show largesplittingsfornitelattiespaing. Indeedinaseofthepiontheselargesplittingshave
been observed [91 ℄. A theoretial analysisintheSymanzikexpansion shemeallows to
identify the leading orders ofuto eets[92 ℄,
m 0 π 2
= m 2 π + a 2 ζ π + O (a 2 m π , a 4 ) ,
(3.38a)m ± π 2
= m 2 π + O (a 2 m 2 π , a 4 ) .
(3.38b)Theleading order of thesplitting,
ζ π
,originates solely in theneutral setorand turnsout to be large. Important to note is that the large size of splitting uto eets is
restritedto thepionmassand relatedquantities.
3.3.3. Phase Diagram
Asfor all Wilson type fermions, twisted mass fermions an exhibit unphysial phases
in their bare parameter spae that is spanned by the hopping parameter, the
lat-tie oupling and the twisted mass. Its vauum struture has been studied
exten-sively [93 , 94 ,95 ,96 , 97 ,98 ℄. From ordinary Wilsonfermions it isknownthat there is
Figure 3.3.:Phasediagramin
(κ, β, µ 0 )
parameterspae, from[96℄. The strongouplingre-gion, i.e. theregion of small lattieoupling
β
, exhibits the Aoki phasein theplaneofvanishingtwistedmassparameter
µ 0
. Forweakerouplingsasurfaeofbulktransitionsisenounteredthat ontainsthelineofritialhopping
parame-ter,
κ c (β)
andextentstoµ 0 6 = 0
. Thewidthof thissurfaevanishestowardstheontinuumlimit.
the so-alled Aoki phase of broken parity-avour symmetry in the strong oupling
re-gion[99 ,100℄. Fortwistedmassfermions,thatadditionallyhavethe
µ 0
-axisintheirbareparametersphasediagram,anotherunphysialphaseappearsforintermediateouplings
that shrinks towards the ontinuum limit. This phaseis a surfae of rst order
transi-tions thatinorporatesthe lineofritialhoppingparameters
κ c (β)
. Thesendingsaresummarised inthe skethingure 3.3.
Thetheoretialexplanation forthesetwo phasesanbeworked outintheframework
ofhiral perturbationtheory[101 ,102℄whihallowstoinvestigatethevauumstruture
of its eetive Lagrangian by means of minimising the orresponding potential. For
this purpose, the ontinuum Lagrangian given in equation (2.30) has, of ourse, to be
supplementedbythe leading utoeets,
L χ = f 2 4
Tr∂ µ Σ∂ µ Σ †
− f 2 4
TrχΣ † + Σχ †
− f 2 4
TrρΣ † + Σρ †
,
(3.39)with
χ = 2B 0 (m1 −
iµτ 3 )
andρ = 2aW 0 1
, slightly adjusting thenotation in [102 ℄.ρ
parametrisestheutoeets. NotethatitenterstheLagrangianinthesamewayasthe
quarkmasssothatthe leadinguto eetsessentiallyleadtoashiftinthebarequark
mass,f.alsothedisussionin[101 ℄. Thepotentialthatanbeextratedfromtheabove
Lagrangian anbeexpressed interms ofa unitfour vetor,
u = (u 0 , u 1 , u 2 , u 3 )
[102 ℄,V = − c 1 u 0 + c 2 u 2 0 + c 3 u 3 .
(3.40)This vetorparametrises
Σ
,Σ = u 0 1 +
iu a τ a ,
(3.41)sothat itsomponents determine the vauum struture,inpartiular,
u 0 ∼ ψψ
,
(3.42a)u 3 ∼
ψγ 5 τ 3 ψ
.
(3.42b)m 0 µ
B
m 0 µ
µ c
Figure3.4.:Possiblephasetransition senariosfor Wilson fermionswith twistedmass term.
Left: Aokiphase. Right: Normalsenario,alsoalledSharpe-Singletonsenario.
(from[102℄)
The oeients
c 1 = 2f 2 (B 0 m q + W 0 a)
andc 3 = 2f 2 B 0 µ 0
are diretly related to thequark mass and uto eets. However, the ruial observation is that the sign of
c 2
is not determined from the Lagrangian. Therefore, the vauum struture an be qualitatively dierent depending on that sign. The two possible senarios have beenmapped out in [102 ℄. The rst possibility,
c 2 > 0
, leads to thepreviously known Aokiphase of broken avour symmetry,
ψγ 5 τ 3 ψ
6
= 0
, if the twisted mass vanishes. Theseondsenario,
c 2 < 0
, orrespondsto the Sharpe-Singleton plane of bulk transitions that has been found for weaker ouplings numerially. This plane inorporates theritial untwisted quark mass and extends to non-vanishing twisted mass. It ends for
some ritial
µ c ∼ a 2
and thus vanishes in the ontinuum limit. A sketh of the twosenariosis showningure3.4.