Localisation in 2+1 dimensional SU(3) pure gauge theory at finite temperature

I study the localisation properties of low Dirac eigenmodes in 2+1 dimensional SU(3) pure gauge theory, both in the low-temperature, confined and chirally-broken phase and in the high-temperature, deconfined and chirally-restored phase, by means of numerical lattice simulations. While these modes are delocalised at low temperature, they become localised at high temperature, up to a critical point in the Dirac spectrum where a BKT-type Anderson transition takes place. All results point to localisation appearing at the deconfinement temperature, and support previous expectations about the close relation between deconfinement, chiral symmetry breaking, and localisation.


Introduction
The close connection between deconfinement and chiral symmetry restoration is an aspect of the finite-temperature transition in QCD that is still poorly understood. As is well known, both the confining and the chiral properties of QCD at vanishing chemical potential change dramatically in the temperature range T ≃ 145-165 MeV: while in the low-temperature phase quarks and gluons are confined within hadrons and chiral symmetry is spontaneously broken, in the high-temperature phase quarks and gluons are liberated into a plasma, and chiral symmetry is restored. Of course, both the spontaneous breaking and the restoration of chiral symmetry should be understood here as approximate, since chiral symmetry is explicitly (albeit softly) broken by the quark masses. Since the transition is an analytic crossover [1], there is no uniquely defined critical temperature. However, defining a pseudocritical temperature T χ for the chiral transition as the position of the peak of the chiral susceptibility, and a pseudocritical temperature T dec for the deconfinement transition as the position of the peak of the quark entropy, one finds that they agree within errors, T χ = T dec = T c ≈ 155 MeV [2].
This relation between these two seemingly unrelated phenomena is not unique to QCD, but it appears quite generally in a variety of gauge theories. In particular, T χ and T dec coincide in certain theories and models where the finite-temperature transition is a genuine phase transition. Examples are provided by SU(3) pure gauge theory in 3+1 [3] and 2+1 dimensions [4], 1 and N f = 3 QCD with unimproved rooted staggered fermions on N t = 4 lattices [5][6][7]. Recent studies indicate that it is so also for the N = 1 SU(2) super-Yang-Mills theory [8]. Another interesting case is that of 3+1 dimensional SU(3) gauge theory 1 Although, strictly speaking, there is no chiral symmetry here since there are no fermions, one can nevertheless define a valence chiral condensate and study its behaviour in the limit of vanishing valence quark mass. be defined; using SU(3) as gauge group, the symmetry class in the classification of Random Matrix Theory is the unitary one, for which Anderson transitions in two spatial dimensions are known to exist [30].
As there is no obstruction to the sea/islands picture to work, then it should just work, leading to the same situation encountered in 3+1 dimensions: at low temperature, chiral symmetry should be spontaneously broken by the accumulation of delocalised Dirac eigenmodes near the origin; at high temperature, the spectral density should vanish near the origin, leading to chiral symmetry restoration, and modes should be localised up to a mobility edge somewhere in the spectrum. It is already known that chiral symmetry is restored at deconfinement [4]. It is the purpose of this work to verify that the localisation properties of the low modes change there as well.
The plan of this paper is the following. In Section 2 I briefly review localisation in disordered systems, with a special focus on lattice gauge theories in 3+1 dimensions. In Section 3 I describe in some detail the specific model under consideration, and expectations about its behaviour. In Section 4 I show numerical results and their analysis. Conclusions and perspectives on future investigations are discussed in Section 5.

Localisation in lattice gauge theories
Localisation is a well known phenomenon in condensed matter physics. It has long been known, since the seminal work of Anderson [31], that the addition of a random on-site potential to the usual tight-binding Hamiltonian causes the localisation of the energy states at the band edge, beyond a critical energy E c called mobility edge, while states in the band centre remain extended (see Ref. [20] for a review). Such "disordered Hamiltonians" aim at describing metals with impurities, and in this context the width W of the distribution of the on-site potential is a measure of the amount of impurities in the system. As the disorder parameter W is increased, the mobility edge moves towards the band centre, and all modes become localised beyond some critical disorder, turning the metal into an insulator. In three (and higher) dimensions the transition between localised and delocalised modes is a secondorder quantum phase transition known as Anderson transition (see Ref. [21] for a review), with a divergent correlation length ξ ∼ |E − E c | −ν characterised by the critical exponent ν.
Localisation in gauge theories was initially studied at zero temperature in investigations of the topological structure of the QCD vacuum (see the review Ref. [32] and references therein). The idea that the finite-temperature transition of QCD could be related to localisation of the low Dirac modes dates back to Refs. [33][34][35]. The first numerical results supporting this idea appeared in Ref. [12], coming from the effective description of QCD via an Instanton Liquid Model, and in Ref. [13], coming from numerical simulations of QCD on a lattice, both in the quenched approximation and with 2+1 flavours of staggered fermions. Further evidence of localisation of the low Dirac modes in the high-temperature phase of QCD was provided by the absence of correlations in the low-lying spectrum of the overlap operator [14], typical of localised modes.
A detailed study of localisation in lattice QCD was undertaken in Ref. [16], using 2+1 flavours of 2-stout improved rooted staggered fermions. There, it was shown that above T c the eigenmodes of the staggered Dirac operator are localised, for (imaginary part of the) eigenvalue up to a temperature-dependent mobility edge, λ c (T ), beyond which modes become extended. The extrapolation of λ c (T ) vanishes at a temperature compatible with T c , in agreement with the absence of localised modes in the low-temperature phase. Localisation was shown to survive the continuum limit, indicating that it is not a lattice artefact. This is also supported by the fact that localised modes have been found with other fermion discretisations, namely with domain-wall fermions [17] and with overlap fermions on twisted-mass Wilson fermion backgrounds [18].
Localisation in the high-temperature, deconfined phase has been observed also in other gauge theories (see Ref. [36] for a recent review), namely SU(2) [15] and SU(3) pure gauge theory [22], and SU(3) with N f = 3 flavours of unimproved rooted staggered fermions on N t = 4 lattices [23]. In these models the finite-temperature transition is a genuine phase transition, which provides a clean-cut setting for investigating the possible coincidence of deconfinement, chiral-symmetry restoration and localisation of the low Dirac modes. While in the first case a detailed study of the temperature dependence of the mobility edge is missing, in the other two cases it was indeed found that the low Dirac modes start localising precisely at the critical temperature where deconfinement and chiral-symmetry restoration take place. A similar result was found in a toy model for QCD [24], where the Polyakov-line dynamics is mimicked by a spin model: as the spins get ordered, the low Dirac modes localise and their density near the origin drops to zero.
The simplest way to identify localised modes is by studying their so-called participation ratio (PR). Given a normalised eigenmode of a lattice Dirac operator, Ψ(n), one defines the inverse participation ratio (IPR) as where Ψ(n) † Ψ(n) denotes the scalar product in colour and (possibly) Dirac space, and the sum is over the lattice sites n. 5 The PR is just the inverse of the IPR divided by the lattice size, where we have assumed that the lattice is a (d+1)-dimensional hypercube of spatial extension N s and temporal extension N t in lattice units. In all the models discussed above d = 3, while in this work I will consider d = 2. For a fully delocalised mode, Ψ(n) † Ψ(n) ∼ 1/V , so that the PR remains constant in the large volume limit. 6 For a mode localised in a spatial region of fixed size v one finds instead PR ∼ v/V → 0 in the large-volume limit.
A useful observation is that in a random matrix model the localisation properties of the eigenmodes and the statistical properties of the corresponding eigenvalues are closely related [37]. For localised modes the eigenvalues fluctuate independently, following Poisson statistics, while for extended modes the eigenvalue statistics are those of the appropriate Gaussian ensemble of Random Matrix Theory (RMT; see Ref. [38] for a general introduction, and Refs. [39,40] for the application of RMT to QCD). This is made evident after unfolding the spectrum, i.e., after mapping is the spectral density and . . . denotes averaging over the random matrix ensemble. This mapping makes the spectral density identically 1 throughout the spectrum. It is well known that for dense matrices the (bulk) statistical properties of the unfolded spectrum are universal and uniform throughout the spectrum, i.e., do not depend on the details of the random matrix model under study, but only on the symmetry class of the ensemble [38]. The main classes are the orthogonal, unitary, and symplectic classes. 7 Universal analytic results can therefore be obtained using the so-called Gaussian ensembles. In particular, the probability distribution of consecutive unfolded level spacings s i ≡ x i+1 −x i is known, and very accurately approximated by the so-called Wigner surmise, which for the unitary class reads In contrast, for eigenvalues that fluctuate independently the unfolded level spacing distribution is the exponential function, appropriate for Poisson statistics, The lattice Dirac operator in a gauge-field background is in practice a sparse random matrix, with fluctuations provided by the gauge links and ensemble averaging corresponding to integration over gauge fields with the appropriate measure. For ensembles of sparse matrices the statistical properties can depend on the spectral region under consideration, and so it is customary to compute the spectral statistics locally in the spectrum, i.e., restricting to small spectral intervals around the chosen point. Looking at Dirac spectra in the high-temperature phase of QCD one sees indeed a transition from Poisson to RMT statistics. By a finite-size scaling study of statistical spectral observables it is possible to determine precisely the location of the mobility edge and the critical exponent ν [42]. This has been used to show that the localisation/delocalisation transition in the Dirac spectrum in high-temperature lattice QCD with staggered fermions is a genuine second-order phase transition in the same universality class [19] as the three-dimensional unitary Anderson model [43]. Further support to this conclusion came from a study of the multifractal properties of eigenmodes at the mobility edge [44]. This matching can be easily understood in the light of the sea/islands picture [26,27] and of the Dirac-Anderson approach [25], since the spatial fluctuations of the Polyakov lines precisely provide the kind of three-dimensional on-site disorder present in the Anderson model, while the unitary symmetry class is the one to which the staggered Dirac operator belongs.

Lattice SU(3) pure gauge theory in 2+1 dimensions
Finite-temperature SU(3) pure gauge theory in 2+1 dimensions has been studied on the lattice in several papers [45][46][47][48][49][50][51][52]. This theory shows a deconfining second-order phase transition at finite temperature, in the same universality class as the two-dimensional 3-colour Potts model. These works used a hypercubic lattice and the Wilson action, which up to an irrelevant additive constant reads

2)
µ, ν = 1, 2, 3 are the lattice directions, U µ (n) are SU(3) gauge links living on the lattice edges (n, n +μ), andμ denotes the unit vector in direction µ. Periodic boundary conditions both in the temporal and in the spatial directions are imposed. In 2+1 dimensions the gauge coupling g has dimensions of [mass] 1 2 , and the lattice coupling β is related to g and the lattice spacing as β = 6/(g 2 a). 8 The partition function is with dU µ (n) the SU(3) Haar measure. The critical temperature was determined precisely in Ref. [49], and for lattices of temporal extension N t = 4 it corresponds to the lattice gauge couplingβ c = 4.9057 (57). Chiral symmetry breaking, as discussed in the Introduction, was first studied in Ref. [4] with staggered fermions in the quenched approximation, so studying the Dirac operator on pure-gauge theory backgrounds. The staggered Dirac operator reads with periodic boundary conditions in the spatial directions and antiperiodic boundary conditions in the temporal direction, and in the continuum limit it describes N f = 2 degenerate species of fermions. The staggered operator D stag has purely imaginary eigenvalues, iλ, and since its spectrum is symmetric about zero it is enough to consider only λ ≥ 0. Fermions break the centre symmetry of the pure-gauge theory, selecting the vacuum with trivial Polyakov lines, but when working in the quenched approximation this has to be done by hand, for example by multiplying all the temporal links in the last time slice by the appropriate centre element. This was the approach adopted in Ref. [4], which we will also use when needed. While the full chiral subgroup of the U(N f ) flavour symmetry discussed in the Introduction is explicitly broken at finite lattice spacing, a remnant U(1) × U(1) chiral symmetry still survives [29], and can be spontaneously broken by the formation of a quark-antiquark condensate. The authors of Ref. [4] observed that while in the low β (low temperature), confined phase a finite chiral condensate was found in the limit of vanishing valence quark mass, this vanished in the high β, deconfined phase, and that the chiral transition coincided with the deconfinement transition. The purpose of this work is to study the localisation properties of the low Dirac modes in the two phases of the theory, both by measuring the PR of the eigenmodes, Eq. (2.2), and by studying the statistical properties of the unfolded spectrum. After obtaining the Dirac spectra for an ensemble of gauge configurations, unfolding is done by sorting all the eigenvalues (for a given coupling and lattice size) by increasing size, and replacing them by their rank divided by the number of configurations. This automatically makes the spectral density equal to 1 throughout the spectrum. I use two quantities extracted from the unfolded level spacing distribution measured locally in the spectrum, P λ (s). One of them is the integrated distribution I λ , which takes the values s 2 RMT ≃ 3π 8 and s 2 Poisson = 2 for the two types of statistics. In the thermodynamic limit, the spectral density and the local average level spacing ∆λ λ ≡ λ i+1 − λ i λ are related as ρ(λ) ∆λ λ = 1. The same relation holds for the unfolded spectrum, and so, since the unfolded spectral density equals 1, one has s λ ≡ ∆x λ = 1 in infinite volume. This might however fail in a finite volume, where the local averaging is necessarily done in a small but finite spectral interval, 9 in regions where the spectral density is small and rapidly varying. Indeed, it is easy to see that the average level spacing in an interval of size ∆ is given by the difference between the smallest eigenvalue λ ′ right above the end of the interval and the smallest eigenvalue λ inside the interval, divided by the number N of eigenvalues inside the interval. The spectral density associated to the given interval is simply N/∆. If the spectral density is increasing, then one typically finds λ ′ − λ < ∆, since eigenvalues get closer going up in the spectrum, and so (N/∆) · [(λ ′ − λ)/N ] < 1, from which we are led to expect ρ(λ) ∆λ λ < 1, and so s λ < 1. This is what happens near the origin in the high temperature phase, where ρ(λ) ∼ λ α for some positive power α: the small spectral density requires the use of relatively large intervals over which the spectral density increases non-negligibly. The value of s λ can then be used to assess the reliability of estimates of spectral statistics based on unfolding, and only those spectral regions where s λ ≃ 1 should be considered in further analyses.
From the point of view of the Dirac-Anderson approach, the model at hand should behave like a two-dimensional unitary Anderson model with on-site (diagonal) disorder. 10 This model has been studied in Ref. [30], where it was shown that it displays an Anderson transition of Berezinskiȋ-Kosterlitz-Thouless (BKT) type [55][56][57], with an exponentially divergent localisation length. We expect to find the same behaviour in our model in the high-β phase, with a mobility edge λ c separating localised (low) and delocalised (high) Dirac modes, where the localisation length should diverge in the infinite-volume limit as for some constant α, as one approaches λ c from the localised side, λ < λ c . Points on the delocalised side, λ > λ c , are all critical in a disorder-driven BKT transition, and so there the statistical properties of the spectrum should be independent of the volume (see Ref. [58]).
To verify if this is the case, I did a finite-size scaling analysis of the spectral observables O λ = I λ , s 2 λ , making the scaling hypothesis with F some analytic function and L the linear size of the system. The one-parameter scaling hypothesis for localisation [59] is usually motivated by assuming that the behaviour of the system near the transition in a finite volume is determined only by the ratio ξ(λ)/L. However, since we are working on a lattice, ξ and L can be made dimensionless dividing them by the lattice spacing, and so one can use instead a function F (r) of ratios of the more general form r = h(ξ(λ))/h(L). Indeed, since near the transition and for large volumes only the leading behaviour of h matters, r is invariant under the rescaling ξ → bξ, L → bL, as it should. On the other hand, analyticity in a finite volume constrains the form of h, which in the case at hand must be of the form

Numerical results
In this Section I report on the results of a numerical investigation of Dirac spectra in finitetemperature SU(3) pure-gauge theory in 2+1 dimensions on the lattice. The gauge back- of configurations  32  100  40000  40  140  40000  48  216  20000  56  280  8000  64  380  4000  72 480 2000 Right above the deconfinement transition, and for the lattice sizes used here, a system initialised in the trivial Polyakov-loop sector is still able to switch to one of the complex sectors, and to temporarily tunnel to the confined phase, even though this becomes less and less likely as the volume is increased. Forβ = 5.05, 5.10, 5.15, 5.20 I then always analysed the gauge configuration obtained by rotating the centre sector to the trivial one, as explained above in Section 3. Forβ = 5.25, 5.50, 5.75, 6.00, 6.25, 6.50 this was not necessary, as a system initialised in the trivial Polyakov-loop sector never left it. Of course, no sector switching was used forβ = 2.50, 3.00, 4.00 in the confined phase. The first few eigenvalues of the staggered Dirac operator were obtained by means of the ARPACK routine [60]. In the following I denote with λ the (imaginary part of the) eigenvalues of the staggered Dirac operator in lattice units. Monte Carlo simulations were performed with standard heatbath [61,62] and overrelaxation [63,64] algorithms. Details about statistics and the number of computed eigenvalues for each volume can be found in Table 1. All statistical errors were obtained through a jackknife analysis with 100 samples. All fits in this work were done using the MINUIT routine [65], establishing the accuracy of the errors with the MINOS subroutine. Whenever the difference is negligible, the symmetric parabolic error is used instead.

Participation ratio
I discuss first the PR of the low modes. In Fig. 1 I show how the average PR of the Dirac modes computed locally in the spectrum, PR λ , changes along the spectrum forβ = 5.50 (in the deconfined phase) and the various lattice volumes. All locally averaged quantities, like PR λ , are obtained by averaging the relevant observable over the modes in spectral intervals of size w = 0.005 and over gauge configurations; the result is assigned to the average eigenvalue (computed similarly) in that interval. In the top panel of Fig. 1 I show PR λ : it is clear that this quantity keeps decreasing with the volume for the lowest modes, while it remains almost constant for the higher modes. In the bottom panel I show instead PR λ · V , i.e., the spatial "size" of the mode in lattice units: this remains constant for the lowest modes, and blows up for the higher modes. This shows that the lowest modes are localised, and the higher ones are delocalised. Furthermore, the size of the modes is approximately constant in the localised part of the spectrum. A similar behaviour is observed at allβ in the deconfined phase. For the lowest two values ofβ there, i.e.,β = 5.05, 5.10, the volume dependence appears to be non-monotonic. This is probably caused by large finite-size effects near the critical temperature (see below).
To give a well-defined estimate of the size of a localised mode in the deconfined phase, I have measured the PR of the lowest mode averaged over configurations times the spatial size of the lattice, PR 1 · V , and extrapolated it to infinite volume. Data points lie approximately β = 5.05 β = 5.10 β = 5.15 β = 5.20 β = 5.25 β = 5.50 β = 5.75 β = 6.00 β = 6.25 β = 6.50   Fig. 1 (bottom). Not surprisingly, this is less so for the lowest two values ofβ. In Fig. 2 (4.1) The linear extrapolation inβ/N s is fully satisfactory for all values ofβ, except for the lowest two, for which the linear fit to the data is of poorer quality. Nonetheless, the result of the extrapolation is in agreement with the general trend of the data, which shows the size of the lowest mode increasing as one gets closer to the critical coupling. If, as expected, all the modes become delocalised at the deconfinement point, then this quantity should blow up there. Near the critical point, the size of the lowest mode in lattice units, ℓ ∞β , can then become comparable to N s for the available lattice sizes, causing sizeable finite-size effects.
The results for ℓ ∞β vary from around 7 for the highest values ofβ, to around 8 for the ones closest toβ c , so there is probably nothing to worry about for the lattice sizes used in this work. On the other hand, if localisation is driven by the behaviour of the Polyakov line, then near the critical point one expects further large finite-size effects due to the large fluctuations of this quantity (see also the discussion below in Section 4.2). In particular, during the numerical simulations the system can still tunnel between the confined and deconfined phases for not large enough volumes. Although these effects are harder to quantify, they are probably responsible for the larger error bars of ℓ Ns nearβ c . The results of the infinite-volume extrapolation of ℓ Ns are shown in Fig. 2 (bottom). For the lowest values ofβ, even though one may not fully trust quantitatively the results due to the above-mentioned finite-size effects, there is a clear sign of a divergence as the critical coupling is approached. At largeβ instead this quantity decreases as 1/β, reflecting the fact that there the extrapolated size of the lowest mode in lattice units is almost independent of the coupling. A constant fit using onlyβ ≥ 5.50 yields ℓ ∞β = 7.155(10) with χ 2 /d.o.f. = 1.4. This leads to a linear dependence of the physical linear size of the mode on the inverse temperature at high temperature, as one would expect given that this is the most important scale there. The same behaviour has been observed in QCD [16]. Plugging in the fit result, one gets ℓ ∞ T = 1.7887 (25), where T ≡β/N t . Despite the limited number of points close tō β c , I have tried to fit the divergent behaviour with the following functional forms, The first one is a simple power-law, while the second one is inspired by the behaviour of the localisation length 12 at criticality. 13 Only points up toβ = 5.75 were included. Both fits converge, with suspiciously small χ 2 /d.o.f. ≃ 0.2 and χ 2 /d.o.f. ≃ 0.5, respectively, indicating 12 It must be remarked that in general the localisation length, ξ, and the quantities ℓ (q) = (IPRq) built out of the generalised IPRs (see footnote 5), do not have the same critical behaviour. In fact, even though all these quantities provide an estimate of the size of the localised modes and should blow up at the critical point, they are sensitive to different features of the modes, which usually results in different critical exponents. 13 The critical behaviour of ξ should not depend on whether the critical line is crossed in the direction of the eigenvalues at fixedβ, or in the direction ofβ at fixed eigenvalue. Having extrapolated the lowest mode to infinite volume, here we are looking at around λ = 0. nontrivial. In Fig. 3 I show PR λ and PR λ · V forβ = 3.00 for the various lattice sizes.
None of these two quantities shows a clear sign of converging to a constant anywhere in the available spectrum, indicating that they show a somewhat intermediate behaviour between localised and fully delocalised. To investigate this issue quantitatively, one can compute the fractal dimension D 2 , related to the PR as [21] PR where ∼ denotes the asymptotic large-volume behaviour and L is the linear size of the system. For fully delocalised modes D 2 = 2, while for localised modes D 2 = 0. The fractal dimension D 2 can be estimated by comparing the PR computed on pairs of lattices of linear spatial sizes N s1 and N s2 , using the formula From this estimator one obtains D 2 in the limit of large N s1,2 . The results forβ = 2.50, 3.00, 4.00 are shown in Fig. 4 for three different pairs of volumes. The volume dependence is mild to non-existent. Low modes have a nontrivial fractal dimension between 0 and 2, while for higher modes D 2 approaches 2. Moreover, low modes are closer and closer to being localised as one approaches the critical temperature from below.

Spectral statistics
As mentioned above in Section 2, spectral statistics can be used to study the localisation properties of the Dirac modes and to determine efficiently the mobility edge and the related critical properties. The typical behaviour of I λ and s 2 λ in the deconfined phase is shown in Fig. 5 forβ = 6.25. Low modes are close to having Poisson statistics, and more and more so as the size of the lattice increases. For every lattice size the statistics change from Poisson to RMT-type, more precisely to GUE-type, as one moves up along the spectrum. At some point the curves corresponding to the various volumes merge, and show little to no dependence on the volume. This is precisely the behaviour expected for a disorder-driven BKT transition in the spectrum. To make this statement quantitative, I verified the scaling hypothesis Eq. (3.8) by fitting the data for I λ for various volumes with the rational function Fits to the data with Eq. (4.5) are quite sensitive to the initial values of the fitting parameters and the choice of fitting interval, and sometimes the error estimate of the parameters provided by MINOS is not accurate, especially close to the critical coupling (the missing point atβ = 5.15 for N s min = 48 in Fig. 6 is precisely a case of inaccurate error estimate). Despite these numerical shortcomings, the results are all quite consistent with the hypothesis ν = 1 2 . In fact, I performed these fits to obtain an unbiased (or at least as little biased as possible) estimate of the critical exponent, as a preliminary step for a more accurate and more systematic finite-size scaling study using constrained fits [66]. These fits turn out to be unstable if ν is not constrained strongly, for which reason I preferred to first check the viability of the theoretical value with unconstrained fits, and then perform the constrained fits fixing ν = 1 2 . Constrained fits were done by fitting the data with polynomial functions, where y(λ, N s ) is given in Eq. (4.5), which is nothing but a truncated Taylor expansion around λ c . The coefficients F n and the mobility edge λ c are the fitting parameters. In general, in a constrained fit one minimises an augmented χ 2 , χ 2 aug = χ 2 + χ 2 prior , where n and σ λc , σ ν , σ n , reflecting one's knowledge of the parameters. The order n max of the polynomial is then increased until the error on the parameters stabilises. This allows to estimate accurately the systematic error due to truncation [66]. In my analysis I used no prior for the first four coefficients and for the mobility edge (so formally setting σ λc = σ 1 = σ 2 = σ 3 = σ 4 = ∞), I fixed ν = 1 2 (so formally setting σ ν = 0), and set F  14 The symmetric parabolic error was used, since it never differed from the MINOS errors by more than 10%. and (0.06, 48,8), respectively. Among the various fitting parameters, the most important ones are the mobility edge, λ c , and the values of the statistics at the critical point, I c ≡ I λc and s 2 c ≡ s 2 λc . These are shown in Figs. 9 and 10. The mobility edge in physical units, λ c phys (β) =βλ c (β), obtained multiplying λ c by the lattice coupling (assuming perfect scaling) is shown in Fig. 11. The values of λ c (β) obtained from the two spectral statistics agree with each other within errors. The values of the spectral statistics at the mobility edge differ from those corresponding to Poisson or RMT statistics, and depend on the lattice coupling (see Figs. 9 and 10). As functions ofβ they are rather flat at largeβ, while they seem to grow as one approaches the critical temperature. This is different from what has been observed in QCD and QCD-like models in 3+1 dimensions: there the statistics at the mobility edge have always been found to be compatible with the critical statistics of the three-dimensional unitary Anderson model [23,67].
Finally, I fitted the mobility edge in lattice units, λ c (β), as obtained from the two spectral statistics, as a function ofβ, trying to establish if it extrapolates to zero at a coupling compatible with the critical one. In order to estimate accurately the errors, I performed another constrained fit with a function of the form increasing m max up to m max = 5. While no prior was assumed on u, v andβ 0 , the parameters h m were constrained to be small, centring their distributions at zero and using widths σ 1 = σ 2 = 0.05, σ 3 = 0.005, σ 4 = 0.001, and σ 5 = 0.0005. These values were chosen so that the fit would converge and errors would be estimated accurately by MINOS. While the resulting errors on h m are bigger than the central values, so that these coefficients are compatible with zero, nevertheless their presence in the fit has a visible impact on the errors of the important parameters. This can be seen in Fig. 12, where I show the dependence on m max ofβ 0 and its error. The final results for u, v andβ 0 , determined from both the I λ and the s 2 λ analyses, are reported in Table 3. The χ 2 aug /n data is around 0.2 in both cases, a suspiciously small value that can however be explained as follows. 15 In Figs. 9 and 10 one sees that the data points atβ = 5.15, 5.20 have a considerably larger error than the other points. This is due to a rather large finite-volume effect: values of λ c determined with N s min = 48 show a visible jump upwards with respect to those determined with N s min = 40, something that does not happen for the other values ofβ. This suggests that it is probably due to a fit instability that would be cured by increasing the statistics, and thus that the finite-volume error is actually overestimated. The results obtained atβ = 5.15, 5.20 with N s min = 40 are probably more reliable than those obtained with N s min = 48, and so data with N s min = 40 have been used as central values. Repeating the fit with Eq. (4.8) using these data and only the statistical errors leads to similar results with a χ 2 aug /n data of about 0.7-0.8.  Table 3. Results forβ 0 , u, and v obtained with a constrained fit of λ c (β), obtained from I λ (left) and s 2 λ (right), with the function Eq. (4.8), using m max = 5.
In Figs. 9 and 10 I also show the fit to the λ c data with Eq. (4.8), marking the critical pointβ 0 at which the mobility edge vanishes, and the corresponding error band. The critical couplingβ c for deconfinement, as determined in Ref. [49], is also shown. The two values agree within one standard deviation, although they just do so. We could just be happy with that and conclude that the mobility edge vanishes at the deconfinement transition, but it is worth trying to explain why the difference is sizeable. According to the sea/islands picture and the Dirac-Anderson approach, the source of disorder leading to localisation in the deconfined phase are the local fluctuations of the Polyakov lines around the ordered value. One then expects the mobility edge to be sensitive to the average Polyakov loop, since this provides a measure of how ordered the system is. Finite-size effects affecting this quantity should then reflect themselves both on the mobility edge and on the pseudocritical coupling in a finite volume. To try to quantify finite-size effects on the latter quantity, I used the results of Ref. [49]. There, the critical coupling in the infinite-volume limit for the 2+1 dimensional SU(3) pure gauge theory is obtained via extrapolation from finite volumes using the formula  Table 3, it should not be taken too literally. The important point is that the deviation from the infinitevolume limit is still about 4% for N s = 72, which is just about the same deviation we find betweenβ 0 and the infinite-volume result forβ c . If anything, this is reinforcing rather than weakening the claim that deconfinement and localisation of the low modes happen together. In this respect, it might be worth mentioning that the exponent v in Table 3, governing the approach to zero of the mobility edge, agrees (within the rather large errors) with the magnetisation critical exponent governing the vanishing of the Polyakov loop expectation value, which in turn has been found to agree with the magnetisation critical exponent 1 9 of the two-dimensional q = 3 Potts model [46]. It would be interesting to investigate this relation further, although this requires numerical simulations close to the phase transition, where they are known to be difficult.
Once the mobility edge is known, the density of localised modes in physical units can be computed as where both λ and the spectral density ρ(λ) = n δ(λ − λ n ) are in lattice units. In  This quantity is independent of temperature within errors, and a simple constant fit gives 16 R = 0.642(85). This means that the density of localised modes behaves like N loc /V phys ∝ λ c phys /T ≃ (T − T c ) v /T , with T c =β c /N t , and since v ∼ 0.1 ÷ 0.2 one has that it rises steeply to a maximum, and then decreases quite fast with temperature. This is different from what was found in QCD in 3+1 dimensions [16]: there, the density of localised modes was seen to keep increasing, up to the highest available temperatures of about 5 T c , while here the decrease begins already at 1.1 ÷ 1.2 T c .
To conclude the discussion of spectral statistics, I present an independent test, usually referred to as shape analysis [68], for the one-parameter scaling hypothesis, Eq. (3.8), exploited above to determine the mobility edge by means of a finite-size scaling analysis. If indeed a single quantity, namely the ratio log ξ/ log L, determines the statistical properties of the spectrum, then plotting one spectral statistic against another should yield universal curves, on which the data points coming from different volumes and lattice couplings should all (approximately) lie, at least for sufficiently large volumes. In Fig. 14 I show plots of s 2 λ against I λ in the deconfined phase for N s = 64, 72. Only points for which | s λ − 1| < 0.05 are plotted. Data from different lattice couplings lie on a common curve connecting the RMT point and the Poisson point. The same kind of behaviour has been observed in 3+1-dimensional QCD [67,69]. The dashed line corresponds to the statistics determined by the so-called "sinh-kernel", 17 a one-parameter family of two-level connected correlators that describes the statistical behaviour in the bulk of the spectrum for several different random matrix models [70][71][72]. This curve runs close to the numerical data obtained in 3+1 dimensional QCD, intersecting them at the critical point [67]. In the present case, this curve describes well the numerical data close to the RMT point, up to (and possibly even beyond) the critical points found for the different values of the lattice coupling. This suggests that the spectral statistics on the line of critical points above the mobility edge belong to the family parameterised by the sinh-kernel. Numerical errors are however still too large to make conclusive statements.
After this long discussion of the high-temperature phase, a few words about the lowtemperature phase are in order. In Fig. 15 I show the spectral statistics I λ and s 2 λ in the confined phase. Analogously to what was found for the fractal dimension of the eigenmodes, the spectral statistics do not show a uniform RMT-type behaviour, but they get closer to Poisson behaviour in the low end of the spectrum, the more so as one gets closer to the deconfinement transition. The transition in the spectrum however does not seem to be a genuine phase transition: the distance from Poisson behaviour remains large even for the largest volumes (as can be seen in the shape-analysis plot in Fig. 16), and attempts at an unconstrained finite-size scaling analysis fail. It is likely that in the infinite-volume limit the spectral statistics converge to a non-trivial behaviour interpolating between Poisson and RMT: this would reflect the non-trivial fractal dimension of the eigenmodes. 18 β = 5.05 β = 5.10 β = 5.15 β = 5.20 β = 5.25 β = 5.50 β = 5.75 β = 6.00 β = 6.25 β = 6.50 critical  Combining the behaviour of the spectral statistics observed in the confined phase with that observed at the mobility edge in the deconfined phase, one is led to conjecture the following development with temperature. The spectral statistics at the low end of the spectrum the RMT predictions for the microscopic statistics and lattice data (at zero temperature) cf. Ref. [73]. keeps approaching the Poisson statistics as the deconfinement transition is approached from below, and reaches it at the critical lattice coupling. At the same time, the crossover along the spectrum turns into a true BKT-type, localisation/delocalisation phase transition, and a mobility edge appears at the origin. Moving further into the deconfined phase, the mobility edge moves up in the spectrum and the statistics at the mobility edge moves towards RMT, apparently stabilising at an intermediate point between Poisson and RMT for sufficiently large temperature.

Conclusions and outlook
The connection between deconfinement, chiral symmetry restoration and localisation of the low Dirac modes has become increasingly evident in recent years, both in QCD [12-14, 16-19, 44] and in QCD-like theories [15,[22][23][24][25]. The theoretical arguments of Refs. [25][26][27] suggest that the driving force behind both chiral symmetry restoration and localisation is the ordering of the Polyakov lines causing deconfinement. This is further supported by numerical studies of dedicated toy models [25]. Given the generality of these arguments, the connection between these three phenomena is expected to be of rather general nature. A non-trivial test of this idea has been provided in the present paper, devoted to the study of localisation of Dirac modes in 2+1 dimensional SU(3) pure gauge theory at finite temperature on the lattice. This model differs from QCD and from the other QCD-like models mentioned above in several aspects: different dimensionality (2+1 instead of 3+1), different critical behaviour at deconfinement (second-order phase transition instead of crossover or first order), 19 different expected type of localisation/delocalisation transition (BKT instead of second order). While the simultaneity of deconfinement and chiral restoration has long been known [4], no previous studies existed about localisation of the Dirac modes. The numerical results presented here, obtained with the staggered discretisation, indicate that the lowest Dirac modes are delocalised (although with non-trivial fractal dimension) in the confined phase, and localised in the deconfined phase. A BKT-type Anderson transition is seen to take place at a critical point (mobility edge) in the spectrum in the high-temperature phase, in agreement with expectations based on universality arguments and on known results for the two-dimensional unitary Anderson model [30]. Although this is perfectly natural in the framework of the Dirac-Anderson approach [25], it is by itself a rather nontrivial finding, which provides nontrivial support to the related sea/islands mechanism [26,27] for localisation in high-temperature gauge theories. Both the inverse of the typical size of localised modes and the mobility edge extrapolate to zero at temperatures compatible with the deconfinement temperature, as determined in Ref. [49], indicating that the onset of localisation coincides with deconfinement and chiral symmetry restoration. This work thus provides further support to the idea of deconfinement driving the system to a chirally restored phase with localised low Dirac modes.
There are other theories where it would be worth studying the relation between deconfinement, chiral symmetry restoration and localisation. The case of gauge theories in the presence of an imaginary chemical potential is interesting both for its theoretical aspects, and for the possible relevance to the study of hadronic matter at finite density. Since the imaginary chemical potential changes the effective boundary conditions affecting the Dirac eigenmodes, according to the sea/islands picture and the Dirac-Anderson approach it should control the density and the localisation properties of the low modes. Preliminary results [36] show that in SU(3) pure gauge theory in 2+1 dimensions an imaginary chemical potential leads to an increase both in the spectral density near the origin and in the size of the low modes, eventually leading to chiral-symmetry breaking and delocalisation of the low modes for sufficiently large imaginary chemical potential, in agreement with expectations.
Another interesting case is that of U(1) pure gauge theory in 2+1 dimensions, for various reasons. The Abelian nature of the gauge group and the different nature (BKT) of the finitetemperature deconfining transition provide a setup qualitatively different from any other investigated so far. Moreover, in contrast with the case of SU(N c ) gauge groups, in the high-temperature phase there is no infinite barrier in the thermodynamic limit separating Polyakov loop sectors that differ by the phase of the spatially-averaged Polyakov loop. Since it is the phases of the Polyakov lines that provide the effective boundary conditions for the Dirac modes, in turn determining their localisation properties, one expects the coexistence of localised and delocalised modes at the low end of the spectrum. Also in this case the available preliminary results [36] confirm the expectations.
An important consequence of localisation of the low modes is that it prevents a spontaneously broken continuous symmetry from generating Goldstone bosons [74,75]. If the low Dirac modes become localised at deconfinement then the Goldstone mechanism does not apply anymore, independently of the spectral density near the origin, and Goldstone bosons disappear from the spectrum. From this point of view, the connection between deconfinement and localisation is possibly even more important than that between localisation and chiral symmetry breaking. A setting in which these issues can be investigated is the SU(3) gauge theory with adjoint fermions in 3+1 dimensions, which has long been known to display separate deconfining and chirally-restoring phase transitions [9]. Preliminary (unpublished) results show that this happens also in 2+1 dimensions with gauge group SU(2) in the quenched approximation, and that while the spectral density remains finite above the deconfinement transition, the near-zero modes become localised.
The results of these paper increase the confusion about the role played by topology in making the low Dirac modes localised. In Refs. [12,13] the localised modes were understood as coming from the localised zero modes supported by instantons at finite temperature. However, the authors of Ref. [15] claimed that the estimated density of topological objects was an order of magnitude smaller than that of localised modes (in quenched SU(2) configurations). More recently, in Ref. [22] (see also Ref. [76]) the density of topological near-zero modes in the pure gauge SU(3) theory was shown to account only up to 60% of the localised modes at the deconfinement transition, a fraction rapidly falling with temperature, and which overestimates the corresponding result in the presence of dynamical fermions. The results of Ref. [17] (for QCD with domain-wall fermions) show that localised modes prefer locations with larger action density and topological charge density, which may be related to the positions of L-L monopoleinstanton pairs. On the other hand, in 2+1 dimensions there is no topological charge for SU(3) fields, since the homotopy group π 2 (SU(3)) is trivial, and nonetheless localised modes appear in the deconfined phase. Further studies are needed to clarify this issue, which, given the close relation between localisation and deconfinement, might even help in understanding the role played by topology in the deconfinement transition.
In conclusion, further studies of localisation in gauge theories may help in better understanding the relation between confinement and chiral symmetry breaking.