# Chiral transition, eigenmode localisation and Anderson-like models

###### Abstract

We discuss chiral symmetry restoration and eigenmode localisation in finite-temperature QCD by looking at the lattice Dirac operator as a random Hamiltonian. We argue that the features of QCD relevant to both phenomena are the presence of order in the Polyakov line configuration, and the correlations that this induces between spatial links across time slices. This ties the fate of chiral symmetry and of localisation of the lowest Dirac eigenmodes to the confining properties of the theory. We then show numerical results obtained in a QCD-inspired Anderson-like toy model, derived by radically simplifying the QCD dynamics while keeping the important features mentioned above. The toy model reproduces all the important qualitative aspects of chiral symmetry breaking and localisation in QCD, thus supporting the central role played by the confinement/deconfinement transition in triggering both phenomena.

Chiral transition, eigenmode localisation and Anderson-like models

Ferenc Pittler

HISKP(Theory), University of Bonn,

Nussallee 14-16, D-53115,Bonn, Germany

E-mail: pittler@hiskp.uni-bonn.de

\abstract@cs

## 1 Introduction

As is well known, in QCD both the confining and chiral properties of
the theory undergo a rapid change around a temperature [1]. Moreover, in models where there is
a genuine phase transition, rather than just an analytic crossover
like in QCD, deconfinement and the chiral transition take place at the
same temperature. Such models include SU and SU pure gauge
theory, SU gauge theory with uninmproved staggered fermions on
lattices [2], and SU gauge theory with
adjoint fermions [3].^{1}^{1}1In the latter case the
deconfinement transition is accompanied by a first-order chiral
transition signalled by a jump in the chiral condensate, while a
second transition, fully restoring chiral symmetry, takes place at a
higher temperature.
A convincing explanation of why deconfinement and the chiral transition
take place together is still lacking.

In recent years there has been growing evidence of a third phenomenon taking place in QCD around , namely the change in the localisation properties of the low-lying eigenmodes of the Dirac operator [4, 5, 6, 7, 8, 9, 10, 11]. While below all the modes are delocalised [12], above the lowest modes are spatially localised on the scale of the inverse temperature, up to a critical, temperature-dependent point in the spectrum (“mobility edge”); modes above are again delocalised. The localisation/delocalisation transition at the mobility edge is a second-order phase transition with divergent correlation length [9]. Extrapolating as is decreased, one finds that vanishes at a temperature compatible with . The low-lying modes play a crucial role in physical observables, in particular in the formation of a chiral condensate, and therefore for the spontaneous breaking of chiral symmetry (SSB). It is therefore not so surprising that the localisation properties of these modes change at the chiral transition.

The appearance of localised modes in the deconfined/chirally restored phase is not a unique feature of QCD, but has been observed also in other QCD-like models [4, 6, 13]. In models where the transition is sharp one can ask more meaningfully whether localisation appears right at the critical point. In the above-mentioned model with unimproved staggered fermions [2], which displays a first-order deconfining and chiral transition, the onset of localisation has been shown to take place precisely at the critical temperature [13]. Understanding localisation might therefore help in understanding the relation between deconfinement and the chiral transition.

## 2 QCD, the Anderson Model and the Dirac-Anderson Hamiltonian

The best known model displaying localisation of eigenmodes is the
celebrated Anderson model (AM) [14], which provides an
approximate description of (non-interacting) electrons in “dirty”
conductors. The effect of impurities is modelled by a random on-site
potential, added to the usual hopping terms of the tight-binding
Hamiltonian; the amount of impurities/disorder is controlled by
the width of the range in which the random potential takes its
values. In the AM, eigenmodes at the band edge are localised, for
energies beyond the so-called mobility edge, and the
localisation/delocalisation transition within the spectrum is a
second-order phase transition with divergent correlation
length. Concerning localisation of the eigenmodes, QCD and the AM are
thus completely analogous. Since in the AM the position of the
mobility edge is determined by , while in QCD it is controlled
by the temperature , one is led to identify as the parameter
effectively controlling the amount of disorder in QCD.
There is, however, more than a simple analogy between the two models.
In fact, the correlation-length critical exponent in QCD was found to
be compatible [9] with that of the 3D unitary
AM [15].^{2}^{2}2The unitary AM is obtained by
multiplying the hopping terms in the AM by random phase factors,
mimicking the presence of a magnetic field.
Here “unitary” refers to the symmetry class of the model in the
classification of Random Matrix Theory (RMT) [12]. Moreover,
eigenmodes at the mobility edge are multifractal [10], with
multifractal exponents matching those of the 3D unitary AM [16].
These results have been obtained using the staggered discretisation of
the Dirac operator.

The fact that the transitions in the two models share the same universality class calls for an explanation. From the point of view of RMT, the QCD staggered Dirac operator is a random matrix of the unitary class, with off-diagonal noise provided by the fluctuations of the gauge links. Moreover, QCD at finite temperature is a 4D model. The unitary AM, on the other hand, is 3D with mostly diagonal disorder. A qualitative explanation for the matching of universality classes is the following [7, 17, 18]. For time-slices are strongly correlated, QCD is effectively 3D, and the quark eigenfunctions are expected to look qualitatively the same on all the time-slices. Moreover, working in the temporal gauge, one sees that obeys effective boundary conditions in the temporal direction, which involve the local Polyakov line (PL), , namely . Here is the temporal extension of the system in lattice units. Since fluctuates in space, it provides an effective 3D source of diagonal disorder. This explanation leads, however, to another question. The effective boundary conditions apply both above and below , and similarly the PL fluctuates in space in both phases of the theory. One thus has to justify why QCD is effectively 4D at low , and why the effective boundary conditions are ineffective there.

In order to better understand the relation between QCD and the unitary AM, following Ref. [18], we make the connection explicit between the staggered Dirac operator and Anderson-type Hamiltonians. We start from the “Hamiltonian” , split it into a “free” and an “interaction” part, , and then work in the basis of the “unperturbed” eigenvectors of . The physically most sensible choice is to identify with the temporal hoppings, , thus leaving the spatial hoppings as the spatially isotropic interaction part, . Here , are SU gauge links with , and are the usual staggered phases. In the temporal diagonal gauge one has , and , , having diagonalised each PL with a time-independent gauge transformation. We choose , any other choice leading to a unitarily equivalent Hamiltonian. In this gauge it is straightforward to diagonalise , finding for the unperturbed eigenvalues

(2.0) |

where are the effective Matsubara frequencies. The various indices refer to the spatial site where the eigenmode is localised, and the colour component and the temporal momentum (TM) to which it corresponds. In the basis of the unperturbed eigenvectors one finds

(2.0) | ||||

with the spatial gauge links in the temporal diagonal gauge.

The “Dirac-Anderson Hamiltonian” of Eq. (2) is a 3D Anderson-type Hamiltonian with internal degrees of freedom, corresponding to colour and TM, and including both diagonal and off-diagonal noise. Unlike the AM, here the control parameter (the temperature) does not change the overall strength of the disorder (which is always bounded in magnitude), but rather its type and its spatial distribution, which are completely different in the confined and deconfined phases of QCD.

In the deconfined phase, the PL gets ordered along , with
fluctuations away from the ordered value forming “islands” of
“wrong” PLs. This directly reduces the density of small unperturbed
eigenvalues. Moreover, it leads to a reduction of the hopping terms’
entries that are off-diagonal in TM, as a consequence of the
correlations among spatial links on different time-slices induced by
the ordering. Since “wrong” PLs yield smaller unperturbed
eigenvalues, one expects them to provide an “energetically”
favourable place for the quark eigenfunctions to live on, thus
allowing for smaller eigenvalues. Moreover, the mixing of different TM
components is harder outside of these islands. In other words, the
“islands” provide localising “traps” for the low eigenmodes.
The approximate decoupling of different TM components leads to having
in practice weakly coupled, effectively 3D systems. The ordering
of the PLs is also expected to affect the spectral density of the low
modes, and therefore the fate of chiral symmetry. Indeed, mixing of TM
components tends to “push” the modes towards the
origin,^{3}^{3}3This can be seen qualitatively using lowest-order
perturbation theory. and its reduction, combined with a smaller
density of low unperturbed modes, leads us to expect a vanishing
spectral density at zero, .

In the low- phase, on the other hand, there is no spatial structure in the diagonal noise, and no strong correlations among time-slices, so no localisation is expected. The various TM components of the wave function are effectively mixed by the hopping terms, and so TM acts effectively as a fourth dimension, i.e., one deals with a single effectively 4D system. The larger density of low unperturbed modes and the fact that they can easily mix is expected to lead to a nonvanishing .

Summarising, the relevant feature both for chiral symmetry restoration and localisation of the lowest modes is expected to be the ordering of the PLs. In the absence of ordering we expect instead to have SSB and delocalisation. This ultimately links both phenomena to the confining properties of the theory, thus providing an explanation for the coincidence of the chiral transition, the onset of localisation, and deconfinement, the latter being the “fundamental” phenomenon on which the other two depend.

## 3 Toy model

In order to check the arguments of the previous Section, we have constructed a QCD-based toy model, keeping only the features that are expected to be relevant for chiral symmetry breaking/restoration and delocalisation/localisation of the lowest modes [18]. As we discussed above, these should be the ordering of the PLs, and the correlation between spatial links on different time slices that this induces. We thus considered a Hamiltonian with the same structure as the Dirac-Anderson Hamiltonian, Eq. (2), and replaced the PL phases and the spatial gauge links with appropriate toy-model counterparts and . In particular, are taken to be the phases of a set of complex-spin variables, constrained by , and are SU matrices. Writing the Dirac-Anderson Hamiltonian as a functional , the toy model Hamiltonian is then defined by .

The dynamics obeyed by the toy-model variables is constructed mimicking that of their counterparts in gauge theory. The important feature of the PL dynamics is the existence of an ordered and of a disordered phase. This is easily achieved for our complex spins by taking a spin model of the following form:

(3.0) |

For this model has a symmetry, which is spontaneously broken at large . The system can then select one of the vacua (corresponding to the center elements of SU along which the PL can align). The important feature of the dynamics of spatial gauge links is instead the appearance of correlations among time-slices when the PLs get ordered. This is achieved by using the Wilson action without spatial plaquettes, and dropping the fermion determinant. The effect of the PLs on the links is reproduced by coupling the toy-model links to the complex spin phases, treated as an external field. Denoting with , we thus define the toy-link action

(3.0) |

where plays the role of gauge coupling, and the expectation value of an observable as

(3.0) |

For our numerical investigation we have employed the minimal version of this toy model with only two colours and two time-slices, . In this case there is a single relevant phase, , and only one relevant Matsubara frequency , and the “unperturbed” eigenvalues are simply . Since we are interested mainly in the dependence on , i.e., on the ordering of the spin system, we fixed , which leads to a critical , separating the disordered () and ordered () phases, corresponding to the confined and deconfined phases of a gauge theory. The “gauge coupling” was also fixed to . We then studied the localisation properties of the low modes, and “chiral symmetry breaking”, defined here as the presence of a nonzero spectral density at the origin.

In Fig. 1, left panel, we show the spectral density at the low end of the spectrum, both in the disordered and in the ordered phases. While in the disordered phase one has a nonvanishing and so “chiral symmetry breaking”, in the ordered phase and “chiral symmetry” is restored. In order to detect localisation of the eigenmodes we have exploited the fact that delocalised modes obey RMT statistics, while localised modes obey Poisson statistics. We have then computed the integrated probability density of the unfolded level spacings , , with . Here is the probability density of computed locally in the spectrum in a small interval around . Our results are shown in Fig. 1, right panel, together with the predictions of Poisson statistics and of the appropriate ensemble of RMT, which in the case at hand is the Gaussian Symplectic Ensemble. The transition from Poisson to RMT statistics is evident, and becomes sharper as the volume of the system is increased, thus signaling the presence of a phase transition in the thermodynamic limit. The mobility edge is determined by the crossing point of the curves corresponding to different volumes, and increases as the ordering of the spins increases.

Our minimal toy model is thus able to reproduce the qualitative features of QCD regarding chiral symmetry and localisation. By tweaking the toy model one can further check the relative importance of the two effects caused by ordering, i.e., the decreased density of low unperturbed modes and the decreased mixing of the wave function’s TM components. The same consequences are expected to be observed in QCD. In Fig. 2 we show the spectral density obtained in several modifications of the toy model. In the top left panel, correlations across time-slices are removed by setting , in the ordered phase of the spin model: as a result, chiral symmetry gets broken. In the top right panel, the mixing between TM components is removed by making the hopping terms diagonal in TM, in the disordered phase: this restores chiral symmetry. In the bottom left panel, boundary conditions in the temporal direction are switched from antiperiodic to periodic, in the ordered phase, so increasing the density of low unperturbed modes: chiral symmetry is broken. Finally, in the bottom right panel we artificially increase the unperturbed eigenvalues, in the disordered phase: chiral symmetry is restored. In all cases, chiral symmetry restoration/breaking is accompanied by localisation/delocalisation of the low modes. In conclusion, both mixing of different TM components of the quark wave function (or lack thereof) and the density of low unperturbed modes are crucial for the fate of chiral symmetry and localisation.

## 4 Conclusions

Recasting the staggered Dirac operator in the form of an Anderson-like Hamiltonian clarifies the relation between the localisation/delocalisation transitions in QCD and in the unitary Anderson model. Moreover, it allows a better understanding of the localisation mechanism at high temperature, and sheds some light on mode delocalisation and on the formation of a chiral condensate via accumulation of small eigenmodes at low temperature. The fate of localisation and chiral symmetry seems to be determined by the ordering of the Polyakov lines, and thus ultimately by the confining properties of the theory. Further work is certainly needed along these lines.

## References

- [1] S. Borsányi, G. Endrődi, Z. Fodor, A. Jakovác, S. D. Katz, S. Krieg, C. Ratti and K. K. Szabó, JHEP 11 (2010) 077 [arXiv:1007.2580 [hep-lat]].
- [2] P. de Forcrand and O. Philipsen, JHEP 11 (2008) 012 [arXiv:0808.1096 [hep-lat]].
- [3] F. Karsch and M. Lütgemeier, Nucl. Phys. B 550 (1999) 449 [hep-lat/9812023].
- [4] A. M. García-García and J. C. Osborn, Phys. Rev. D 75 (2007) 034503 [hep-lat/0611019].
- [5] T. G. Kovács, Phys. Rev. Lett. 104 (2010) 031601 [arXiv:0906.5373 [hep-lat]].
- [6] T. G. Kovács and F. Pittler, Phys. Rev. Lett. 105 (2010) 192001 [arXiv:1006.1205 [hep-lat]].
- [7] F. Bruckmann, T. G. Kovács and S. Schierenberg, Phys. Rev. D 84 (2011) 034505 [arXiv:1105.5336 [hep-lat]].
- [8] T. G. Kovács and F. Pittler, Phys. Rev. D 86 (2012) 114515 [arXiv:1208.3475 [hep-lat]].
- [9] M. Giordano, T. G. Kovács and F. Pittler, Phys. Rev. Lett. 112 (2014) 102002 [arXiv:1312.1179 [hep-lat]].
- [10] L. Ujfalusi, M. Giordano, F. Pittler, T. G. Kovács and I. Varga, Phys. Rev. D 92 (2015) 094513 [arXiv:1507.02162 [cond-mat.dis-nn]].
- [11] G. Cossu and S. Hashimoto, JHEP 06 (2016) 056 [arXiv:1604.00768 [hep-lat]].
- [12] J. J. M. Verbaarschot and T. Wettig, Ann. Rev. Nucl. Part. Sci. 50 (2000) 343 [hep-ph/0003017].
- [13] M. Giordano, T. G. Kovács, S. D. Katz and F. Pittler, PoS LATTICE 2014 (2014) 214 [arXiv:1410.8392 [hep-lat]].
- [14] P. W. Anderson, Phys. Rev. 109 (1958) 1492.
- [15] K. Slevin and T. Ohtsuki, Phys. Rev. Lett. 78 (1997) 4083 [cond-mat/9704192 [cond-mat.dis-nn]].
- [16] L. Ujfalusi and I. Varga, Phys. Rev. B 91 (2015) 184206 [arXiv:1501.02147 [cond-mat.dis-nn]].
- [17] M. Giordano, T. G. Kovács and F. Pittler, JHEP 04 (2015) 112 [arXiv:1502.02532 [hep-lat]].
- [18] M. Giordano, T. G. Kovács and F. Pittler, JHEP 06 (2016) 007 [arXiv:1603.09548 [hep-lat]].