An IsingAnderson model of localisation in hightemperature QCD
Abstract
We discuss a possible mechanism leading to localisation of the lowlying Dirac eigenmodes in hightemperature lattice QCD, based on the spatial fluctuations of the local Polyakov lines in the partially ordered configurations above . This mechanism provides a qualitative explanation of the dependence of localisation on the temperature and on the lattice spacing, and also of the phase diagram of QCD with an imaginary chemical potential. To test the viability of this mechanism we propose a threedimensional effective, Andersonlike model, mimicking the effect of the Polyakov lines on the quarks. The diagonal, onsite disorder is governed by a threedimensional Isinglike spin model with continuous spins. Our numerical results show that localised modes are indeed present in the ordered phase of the Ising model, thus supporting the proposed mechanism for localisation in QCD.
Keywords:
Lattice QCD, Phase Diagram of QCD, Random Systems1 Introduction
The spectrum of the lattice Dirac operator plays a prominent role in current attempts to improve our understanding of the spontaneous breaking of chiral symmetry in QCD. The key relation in this context is the celebrated Banks–Casher formula BC (), which clarifies the relevance of the lowlying Dirac eigenmodes for the development of a nonvanishing chiral condensate in the chiral limit. The Banks–Casher relation also suggests that the (pseudo)critical behaviour of the theory at the chiral transition/crossover, even away from the chiral limit, will be mostly determined by the behaviour of these modes.
It has long been known that in QCD, below the pseudocritical temperature, , the lowlying eigenmodes of the Dirac operator are delocalised, i.e., they extend throughout the whole fourdimensional volume of the system VWrev () (see also Ref. deF ()). In recent years it has been observed that above this is no longer true: the lowestlying modes are localised on the scale of the inverse temperature KP2 (); feri (). Delocalised modes are still present, but only above a (temperaturedependent) critical point in the spectrum. A similar change in the properties of the lowlying modes is found also in other, QCDlike theories, e.g., with gauge group, and/or in the quenched limit GGO (); GGO2 (); KGT (); KP ().
As the temperature is decreased, moves towards the origin, and vanishes at a temperature compatible with , the crossover temperature determined from thermodynamic observables Aoki:2005vt (); Borsanyi:2010cj (). This is certainly not a coincidence. First of all, in QCD GGO2 (); KP2 () and in all the QCDlike theories where localised modes have been observed (puregauge KGT (); KP () and GGO2 () theories), they are present only in the “chirallysymmetric phase”, independently of the specific lattice discretisation of the Dirac operator (staggered GGO2 (); KP (), overlap KGT (), and recently also Möbius domain wall Cossu:2014aua ()). Moreover, further evidence of a very close relationship between the appearence of localised modes and the chiral transition has been found recently GKKP () in a toy model for QCD, namely unimproved staggered fermions unimproved (), which displays a genuine firstorder chiral transition accompanied by the appearence of localised modes at the low end of the spectrum. However, the precise nature of this relationship has not been fully understood yet.
At fixed temperature, the transition in the spectrum from localised to delocalised modes has been shown to be a true, secondorder Andersontype transition crit (). Such transitions have been most extensively studied before using the Anderson model Anderson58 (); LR (); EM (), which is a model for electrons in a crystal with disorder. The Hamiltonian of the Anderson model consists of the usual tightbinding Hamiltonian plus an onsite (diagonal) random potential mimicking the presence of impurities in the crystal. In three dimensions, when the diagonal disorder is switched on, localised modes appear at the band edges, beyond critical energies called “mobility edges”. As the width of the distribution of the diagonal disorder is increased, i.e., as the system is made more disordered, the mobility edges move towards the band center, and beyond a critical value for the disorder all the modes become localised.
The model described above is the socalled orthogonal Anderson model. A variant of this model is the unitary Anderson model, in which the hopping terms are modified with the introduction of random phases to account for the presence of a magnetic field. The specifiers “orthogonal” and “unitary” refer to the symmetry class of the model in the language of random matrix theory (RMT) Mehta (). In this classification, QCD belongs to the unitary class. RMT makes universal predictions for certain statistical properties of the spectrum of a random Hamiltonian with i.i.d. entries, depending on its symmetry class. These predictions correctly describe also the spectral region corresponding to the delocalised modes of the Anderson Hamiltonian, despite this being a sparse matrix. However, this can be understood by noticing that delocalised eigenmodes are easily mixed by fluctuations in the disorder. On the other hand, localised modes are sensitive only to fluctuations taking place within their support, and thus they hardly mix. As a consequence, the corresponding eigenvalues are expected to fluctuate independently, thus obeying Poisson statistics. This behaviour is actually observed in the localised part of the spectrum of the Anderson model. Therefore, a change in the spectral statistics takes place in correspondence to the localisation/delocalisation transition.
The staggered Dirac operator in lattice QCD, being antiHermitian, admits a straightforward interpretation as a random Hamiltonian. As in the Anderson model, the statistical properties of the spectrum change from Poissontype in the spectral region of localised modes, to the appropriate RMTtype in the region of delocalised modes. Besides this, however, at first sight one can find very little in common with the threedimensional Anderson model. For one thing, the dimensionalities apparently do not match. Moreover, the Dirac operator has only offdiagonal non vanishing terms, i.e., the fluctuating hopping terms representing the gauge fields. In this case the location of the mobility edge was found to be controlled by the physical temperature. Finally, unlike in the Anderson model, the disorder in the gauge links is correlated. It is therefore surprising that the critical exponent of the correlation length has been found to be compatible with that of the unitary Anderson model in three dimensions nu_unitary (), suggesting that the two models belong to the same universality class crit (). The unitary class is indeed the symmetry class appropriate for QCD, and correlations in the disorder should not affect the critical behaviour since they are shortranged, as the disorder consists of gaugefield fluctuations. On the other hand, the fact that a fourdimensional model with only offdiagonal disorder behaves like a threedimensional one with diagonal disorder calls for an explanation.^{1}^{1}1Although the unitary Anderson model contains also offdiagonal disorder, this is known to be much less effective in producing localisation offdiag (); offdiag2 ().
An important step towards clarifying these issues would be the determination of the mechanism responsible for the localisation of the lowlying modes. It was initially suggested GGO () that instantons could play the role of localising “defects”, but numerical investigations by one of the authors and collaborators Bruckmann:2011cc () found a large mismatch between the density of localised modes and the density of instantons, thus making this explanation unlikely. In the same paper, an alternative proposal was made to explain localisation at high temperature, which involves the behaviour of local fluctuations of Polyakov loops above . The basic idea is that, in the ordered phase, local fluctuations of the Polyakov loop provide a localising “trap” for eigenmodes. This is because the “islands” where the Polyakov loop fluctuates away from the ordered value provide “energetically” favourable spatial regions for the quark eigenmodes to live in. Numerical support to this proposal was also obtained, in the case of the gauge group.
The purpose of this paper is to refine this proposal, and extend it to the case. This simple generalisation has nevertheless the merit of clarifying some important aspects that went previously unnoticed. To further test the viability of our explanation, we construct a simple threedimensional model which should display localisation precisely through the mechanism proposed for QCD. Our numerical results show that this is indeed the case, and that the qualitative features of the eigenmodes and of the spectrum correspond to those of the eigenmodes and of the spectrum of the Dirac operator. In particular, our toy model displays a transition in the spectrum from localised to delocalised modes, which is likely to become a true phase transition in the thermodynamic limit.
The plan of the paper is the following. In Section 2 we discuss the mechanism for localisation in hightemperature QCD based on fluctuations of Polyakov lines. In Section 3 we construct an effective model (“IsingAnderson” model) which should produce localised states precisely through the mechanism proposed in Section 2. In Section 4 we report on the results of numerical simulations of the IsingAnderson model. Finally, in Section 5 we state our conclusions and discuss possible extensions of the present study.
2 Polyakov lines and localisation
In this Section we argue that hightemperature QCD contains an effectively diagonal and threedimensional disorder, which explains why the critical behaviour at the localisation/delocalisation transition in the spectrum corresponds to that of the threedimensional unitary Anderson model.
An intuitive argument for the effectively reduced dimensionality of hightemperature QCD is the following: as the size of the temporal direction is smaller than the correlation length, the time slices are strongly correlated and the system is effectively threedimensional. This implies that the Dirac eigenmodes should look qualitatively the same on all time slices, in particular for what concerns their localisation properties.
This intuition is strengthened when we consider the role of the antiperiodic boundary conditions in the temporal direction. To this end, it is convenient to work in the temporal gauge, for and , in which equals the local Polyakov line . A further timeindependent gauge transformation allows one to diagonalise each local Polyakov line: we will refer to this as the diagonal temporal gauge. In any temporal gauge, covariant time differences are replaced by ordinary differences for all . One can do the same also at , if at the same time one trades the antiperiodic boundary conditions for effective, dependent boundary conditions, which involve the local Polyakov line,
(1) 
Since the time slices are strongly correlated, these effective, dependent boundary conditions will affect the behaviour at the spatial point for all times . Furthermore, fluctuates in space (and obviously from one configuration to another). From the point of view of a disordered system, QCD above therefore contains effectively a diagonal (onsite), threedimensional source of disorder.
To see how these effective boundary conditions affect the quark wave functions, let us first discuss a simplified setting in which the Dirac equation can be explicitly solved, generalising the argument of Ref. Bruckmann:2011cc () to . We discuss here the naïve lattice Dirac operator for simplicity, but the considerations of the present Section clearly apply to the staggered Dirac operator as well, since the two are essentially connected by a unitary transformation. Consider configurations with constant temporal links and trivial spatial links , . In the diagonal temporal gauge, one has with , so the Dirac operator is diagonal in colour, and the colour components of the quark eigenfunctions decouple. The eigenfunctions of the naïve lattice Dirac operator have the following simple factorised form,
(2)  
Here and throughout the paper the eigenvalues are expressed in lattice units. The spacetime dependence is fully contained in the plane waves , while the spin and colour dependence are encoded in the bispinors , , with , and in the colour vectors , , with , respectively. The spatial momenta satisfy , being the spatial linear size of the lattice (), to fulfill the spatial periodic boundary conditions, while
(3) 
where is the lattice spacing and the temporal size of the lattice, to fulfill the effective temporal boundary conditions. We will refer to as the effective Matsubara frequencies. For a given value of the phase , there are different branches for , corresponding, however, to only different eigenvalues, since the effective Matsubara frequencies obey the relation
(4) 
Notice that the lowest positive branch of eigenvalues is described by the function
(5) 
which decreases as moves away from zero, where it is maximal, and it is minimal when . The bispinors , , satisfy the equation
(6)  
where , are Euclidean Dirac matrices. The spin index labels eigenmodes corresponding to the same eigenvalue . Taking into account the degeneracy with respect to mentioned above, the eigenvalues are therefore fourfold degenerate. Finally, the colour vectors are trivially , , and satisfy .
Let us now consider the qualitative features of the eigenmodes in the general case. Above , a typical gauge configuration consists of a “sea” where the Polyakov line gets ordered around , which percolates through the lattice, with “islands” of “wrong” . If were completely ordered, and spatial links were trivial, there would be a sharp gap in the spectrum at (and a symmetric one at ), corresponding to the lowest positive branch at in Eq. (2) (see also Eq. (5)). This is a “zerothorder” picture, which is modified by the fluctuations of and of the spatial links. The “firstorder” picture is obtained by allowing to fluctuate while still keeping the spatial links trivial. It is clear from Eq. (5) that the quark eigenfunctions can exploit the “islands” of “wrong” , which are “energetically” favourable, to lower their eigenvalues. In particular, localising the wavefunction on the “islands” can achieve a large eigenvalue reduction, if the momentum required to localise the state is not too large, which is the case, for example, if the “islands” are sufficiently big.^{2}^{2}2 This can be seen by considering a configuration with uniform “islands” of “wrong” Polyakov lines, and a test wavefunction of the form given in Eq. (2) but localised and constant on one of such “islands”. More precisely, take , with chosen so that belongs to the lowest positive branch. Here is 1 on the given “island” and 0 elsewhere. Computing the expectation value of , one finds a value smaller than if the “island” is big enough, so that surface effects at the interface with the “sea” do not overbalance the negative difference . Therefore, in the “firstorder” picture both localised and delocalised eigenmodes appear between and , with a few localised low modes well separated from the bulk of delocalised modes. The full, “secondorder” picture is finally obtained by switching on the fluctuations of the spatial links. These fluctuations have most likely a delocalising effect: for example, they allow different colour components to mix (recall that we are working in the diagonal temporal gauge). Nevertheless, the lowest modes can still remain localised, due to the large energy difference with the bulk of delocalised states. The sharp gap of the “zerothorder” picture has turned into an “effective gap” , which we identify with the mobility edge separating localised and delocalised modes.
To better understand this picture, it is useful to set up an analogy between the Dirac operator and the Hamiltonian of an electron in a crystal in the tightbinding approximation (TBA). If one discards the spatial links, setting them to zero, the solutions of the Dirac equation are localised at a given spatial site, and are given by Eq. (2) with , using the local Polyakov line as the gauge background. These “free” solutions correspond to the atomic orbitals in the TBA, which are localised on a crystal site. In turn, the effect of the spatial links corresponds to that of the hopping terms in the TBA, which allow the electron to hop between crystal sites. In this analogy, the “islands” of “wrong” Polyakov lines correspond to “defects” in the crystal (i.e., the diagonal disorder in the Anderson model), and they act in the same way as localising “traps” for the eigenmodes. Given the threedimensional nature of these “islands”, the qualitative picture described above allows one to understand why, in the hightemperature regime, the critical properties of the Dirac spectrum at the localisation/delocalisation transition correspond to those of the threedimensional unitary Anderson model. In order to check the viability of our explanation, in the following Section we study an effective model which should display localisation precisely through the mechanism that we have proposed. In the remaining part of this Section we add a few comments.
The dependence of on the temperature and on the lattice spacing is expected to correspond qualitatively to that of , which should somehow “drag” the mobility edge along. In the “zerothorder” picture, increases when increasing the physical temperature at fixed , and it decreases as the lattice spacing is decreased at fixed physical temperature. This behaviour is qualitatively the same as that of the Polyakovloop expectation value, which is a measure of the ordering of the Polyakov lines in the gauge configurations. Therefore, even though is determined assuming a perfectly ordered system, we can loosely say that it responds in a direct way to the ordering of the Polyakov lines. In the “firstorder” picture, the separation point between localised and delocalised modes lies below . For a delocalised mode with low spatial momentum, which sees all the fluctuations of the Polyakov lines equally, the gain in “energy” will average out, so bringing near the average of the lowest positive branch. In the “secondorder” picture, will be displaced further down to by the mode mixing, but the size of this displacement depends on the effectiveness of the mixing, which is difficult to assess. It is likely, however, that an increase in the ordering of the Polyakov lines will also make more correlated the spatial links at a given spatial site but on different time slices. The fluctuations of the spatial hopping should therefore become less pronounced, which in turn should reduce the effectiveness of the mixing, and so reduce the distance between and . In this case, both the “firstorder” and the “secondorder” effects tend to reduce the distance between and as the system is made more ordered, and to increase it when the ordering is reduced, thus acting in the same direction as the “zerothorder” effect. The bottom line is that the response of to a change in the ordering of the system is expected to be qualitatively the same as that of , i.e., it goes up with at fixed , and it goes down with at fixed . This expectation is matched by the results of numerical simulations KP2 ().
The qualitative correctness of the “zerothorder” picture allows also a qualitative understanding of the dependence of the (pseudo)critical temperature on an imaginary chemical potential. As we have said above, in the hightemperature phase where the Polyakov line gets ordered, the position of the mobility edge, , is affected by the effective Matsubara frequency corresponding to the trivial Polyakov line. At , the “twist” of the quark wavefunction required by the effective boundary conditions is maximal, and so the quantity , defined in Eq. (5), which “drags” , is maximal. Introducing an imaginary chemical potential (small enough in order to avoid the RobergeWeiss critical line RW ()) is equivalent to adding an extra phase to the effective boundary conditions, which reduces the “twist”. As a consequence, , and in particular , i.e., it diminishes (independently of the sign of ), and so one expects to become smaller.^{3}^{3}3It is possible to argue that the expected dependence of on is the same also in the “firstorder” picture. The “firstorder” correction changes the sharp gap into an effective gap , corresponding to the mobility edge in the “firstorder” picture. As we discussed above, should be close to , or, at nonzero , to . The probability distribution function of the Polyakovline phases is a periodic function, and due to chargeconjugation invariance. This implies . In the ordered phase, even at nonzero , is peaked around zero, approximately symmetric, i.e., , and most likely monotonically decreasing in the interval . Exploiting these properties, one can show that is a decreasing function of . If the relation between the chiral transition and the appearence of localised modes holds true, then one expects to increase with . This is actually what has been observed in lattice simulations immu (); immu2 (). The presence of the RobergeWeiss critical lines can also be understood in this framework, noticing that the system prefers to have an effective gap as large as possible due to the fermionic determinant. Comparing for , corresponding to the Polyakov line sectors , one finds that it becomes “energetically” favourable to switch to the sector when exceeds , and similarly to the sector, and again to the trivial sector, when exceeds and , respectively.
As one would expect, our qualitative picture cannot explain all the features of the localised modes, especially when there are competing effects at play. As it has been shown in Ref. KP2 (), the lowestlying Dirac eigenmodes should remain localised also in the continuum limit. While goes to zero as , in agreement with our qualitative explanation, the “renormalised mobility edge”, , with the bare lightquark mass, has a finite continuum limit, and the number of localised modes per unit volume, ,
(7) 
also remains finite in the continuum. Here is the spectral density in lattice units. Within our picture, the fact that localisation survives the continuum limit is the nontrivial outcome of two competing effects. In fact, while the ordering of the Polyakov lines tends to disappear as is reduced, thus lowering , at the same time at the low end of the spectrum is increased, due to the larger amount of fluctuations, and so of “islands” that can support localised modes. In Ref. KP2 () it has also been observed that increases rapidly with . At fixed lattice spacing, increasing the temperature brings up, while at the same time decreasing the spectral density. Both phenomena are in agreement with our explanation, since increasing the ordering of the Polyakov lines, besides pushing up the mobility edge, also reduces the density of fluctuations and thus the spectral density of localised modes.^{4}^{4}4 Notice that since the Polyakovline phases are continuous variables, it is possible to have a finite amount of “islands” also at very high temperature, as their “energy” cost can be contained by making their deviation from the ordered value arbitrarily small. The value of at fixed results from the balance of these two effects. On top of this, to determine the physical value of one has to take the continuum limit, which has an effect on and opposite to that of increasing the temperature, as we have said above. The actual behaviour of as a function of the temperature ultimately depends on the relative magnitude of all these effects, which cannot be determined in our simple picture.
3 IsingAnderson model
The considerations of the previous Section suggest that it should be possible to understand the qualitative features of the Dirac spectrum and eigenfunctions in QCD, for what concerns the localisation properties, by using a genuinely threedimensional model. To construct such a model one has to strip off all the features that are irrelevant to localisation.
The first step is to get rid of the time direction, reducing the lattice to three dimensions, and replacing the time covariant derivative in the Dirac operator with a diagonal noise term, intended to mimic the effective boundary conditions. This is motivated by the strong correlation among timeslices, as already mentioned in the previous Section.
Furthermore, it is known that, in general, offdiagonal disorder is less effective than diagonal disorder in producing localisation offdiag (); offdiag2 (), so it should be safe to replace spatial covariant derivatives with ordinary derivatives. More precisely, in our case the structure of the hopping terms is of the form considered in Ref. GGC (). It is shown there that the width of the disorder distribution has to be rather large to produce localised modes near through offdiagonal disorder only. This is definitely not the case in QCD, where the disorder in the hopping terms involves unitary matrices. As we have already remarked, fluctuations of the spatial links have most likely a delocalising effect on the low modes. While this is certainly important for the detailed features of the spectrum and of the eigenmodes observed in QCD, in our picture it can nevertheless be regarded as a “secondorder” effect, acting on the localised modes which should be produced by the “firstorder” effect, i.e., the fluctuations of the Polyakov lines. To be precise, it is thus the “firstorder” picture, discussed in the previous Section, that we are going to test here. Notice that replacing the spatial covariant derivatives with ordinary derivatives decouples the colour components of the quark wavefunction: this changes the symmetry class (in the sense of random matrix models), as we discuss in detail below.
If our “sea/islands” mechanism is viable, then the main features of localisation should still be captured by a genuine threedimensional model of the following general form,
(8) 
if the diagonal noise term has the same features as the “diagonal noise”, i.e., the effective boundary conditions, appearing in the true “Hamiltonian” . The relevant features are the following:

the diagonal noise involves continuous variables, namely the phases of the local Polyakov lines;

the diagonal noise terms are correlated, and are governed by the dynamics of the local Polyakov lines;

as the system is made more ordered, the diagonal noise tends to introduce a gap in the spectrum.
The simplest way to incorporate these features in is to base the diagonal noise on a spin model with continuous spins: this obviously satisfies the first requirement, and also the second one as Polyakov lines in the hightemperature phase display indeed a spinmodel type dynamics Yaffe:1982qf (); DeGrand:1983fk (). Finally, also the third requirement is easily implemented if we take
(9) 
where is the spin variable at point , and is a constant determining the strength of the coupling of the fermions to the spins. The distribution of the spins is determined by the dynamics of an Isinglike model with continuous spins, which we take of the simplest possible form with nearestneighbour interactions only,
(10) 
Since is 1 for “aligned” spins (i.e., for ), this choice provides indeed an effective spectral gap when the spins are ordered.^{5}^{5}5It is understood that we work with magnetic field . In the ordered phase of the Ising model there is a “sea” of spins with “islands” of spins, so the underlying configurations have the same features as the Polyakov line configurations in QCD. In QCD we have a single parameter governing the ordering of the configuration and the size of the effective gap. In the effective model the ordering of the spin configuration is governed by , while the size of the gap is mainly determined by the spinfermion coupling , although one expects it to be affected also by the magnetisation of the system, and thus by .
We notice in passing a curious feature of this model: it belongs to different symmetry classes for lattices of even or odd size , namely to the orthogonal class for even, and to the symplectic class for odd (see Appendix A). For our purposes it is convenient to work with evensized lattices, which allow to get rid of the Dirac degree of freedom through a spin diagonalisation.
Let us now describe the spectrum of our effective model in more detail. The spectrum is symmetric with respect to the origin on a configurationbyconfiguration basis, due to the fact that . In order to find the eigenvalues, it is convenient to first perform a spin diagonalisation through the following unitary transformation,
(11) 
which yields
(12) 
where is here the identity in Dirac space, and
(13) 
The Hamiltonian obviously commutes with , so one can diagonalise them simultaneously. Denoting with the common eigenvectors, , , one has for the bispinors
(14) 
and the identical singlecomponent eigenvalue equations
(15) 
where
(16) 
It is immediate to notice that there is an exact twofold degeneracy of the eigenvalues. It is also straightforward to prove that , so setting one finds
(17) 
i.e., is an eigenvector of with eigenvalue . This shows again that the spectrum is symmetric with respect to the origin.
It is less straightforward to show that for a given noise configuration and a given eigenvector with eigenvalue , there is another configuration with noise and with the same Boltzmann weight, and an eigenvector corresponding to the eigenvalue . Indeed, taking and as follows,
(18) 
one can explicitly verify that this is the case. This implies that the average spectrum of (and so that of ), obtained by integrating over the disorder (i.e., over the spin configurations), has an exact symmetry with respect to the origin.
4 Numerical results
We have performed numerical simulations of the effective IsingAnderson model discussed in the previous Section, Eqs. (8), (9) and (10). For our purposes it was sufficient to focus on , Eq. (16), which contains all the relevant information on the localisation properties of the eigenmodes. To form an overall picture of the properties of the effective model, we have performed full diagonalisation of the Hamiltonian on mediumsize lattices ( and ), for inverse temperatures , in the ordered phase of the continuousspin Ising system, and spinfermion coupling , on 300 spin configurations. We have furthermore fully diagonalised on lattices of the same size for , well above the critical inverse temperature, setting .
In Fig. 1 we show the spectral density per unit volume of the system. Low modes have small spectral density, which rapidly increases as one goes up in the spectrum, as expected. Decreasing the temperature, thus making the system more ordered, decreases the spectral density of low modes, as one expects if these modes are localised on fluctuations of the spin variables. Increasing the spinfermion coupling enlargens the region where the spectral density is small, again as expected, since it should push the effective gap up in the spectrum. The same effect is obtained by increasing , i.e., by making the system more ordered. Notice the symmetry under of the average spectrum, which has been discussed in the previous Section. There is apparently also a sharp gap in the spectrum, which is however not very sensitive to , , or to the size of the system. This gap is probably a feature of the effective model, due to the drastic simplifications of the spatial hopping terms; in any case, it is irrelevant for our purposes.
In Figs. 2 and 3 we show the average participation ratio , which gives a measure of the fraction of space occupied by a given mode, as a function of its location in the spectrum.^{6}^{6}6Normalisation to 1 of the eigenfunction is understood. The changes by two orders of magnitude as one moves up in the spectrum starting from the lowdensity region; more importantly, when increasing the size of the system it remains almost constant in the bulk of the spectrum, while it visibly decreases near the origin. This signals that modes in the bulk are delocalised, while modes near the origin are localised. Notice in Fig. 3 the upturn of the at the low end of the spectrum: a similar phenomenon has been observed in QCD with staggered fermions KP2 ().
Finally, in Fig. 4 we show a suitably defined local spectral statistics across the spectrum. Spectral statistics can be used to detect a localisation/delocalisation transition, since the eigenvalues corresponding to localised or delocalised eigenmodes are expected to obey different statistics, namely Poisson or WignerDyson statistics, respectively. More precisely, is defined as
(19) 
where the unfolded level spacing distribution (ULSD) is the probability distribution, computed locally in the spectrum, of the socalled unfolded level spacing , i.e., the level spacing divided by the local average level spacing . The quantity is simply the integrated ULSD up to the crossing point of the exponential distribution, corresponding to Poisson statistics, and the orthogonal Wigner surmise, which accurately describes the ULSD for (orthogonal) WignerDyson statistics. The results confirm that the eigenmodes change from localised to delocalised when one moves up in the spectrum. Moreover, the mobility edge separating localised and delocalised modes goes up in the spectrum when the Ising system is made more ordered by decreasing the temperature, or when the spinfermion coupling is increased. Our results give also a first indication that the slope of the curve increases as the volume is increased, thus hinting at the existence of a true phase transition in the spectrum. Notice that the transition takes place near the point in the spectrum, thus indicating that the position of the mobility edge is connected to the magnetisation of the system, as expected (see Table 1). This is also in agreement with our expectation for the position of the mobility edge in the “firstorder” picture, discussed previously in Section 2.
To further investigate the onset of critical behaviour, we have determined eigenvalues and eigenvectors on larger volumes () in the vicinity of the point in the spectrum where the localisation properties show a dramatic change, for and . We used 20k configurations for , scaling down the statistics for larger volumes by keeping approximately constant the total number of eigenvalues computed in the relevant spectral window. In Fig. 5 we show as obtained from these simulations: one can easily appreciate the tendency of the crossover from Poisson to WignerDyson statistics to become steeper as the volume is increased.
\backslashbox[0pt][lr]  1.5  2.0 

0.65  1.1153(5)  1.4870(7) 
0.70  1.1903(3)  1.5871(4) 
1.00  1.3386(1)  1.7849(2) 
Although the data are not of sufficiently good quality to allow a proper finitesizescaling analysis, it is still possible to make a rough qualitative attempt at showing compatibility of the critical properties with those of the threedimensional orthogonal Anderson model. According to the oneparameter scaling hypothesis, data from different volumes should collapse on a single curve if plotted against the scaling variable , where is the critical point and the critical exponent of the correlation length. The pretty large distance between the crossing point of the data and of the data indicates the presence of sizable finitesize corrections to the position of the critical point, so we limit ourselves to the two largest volumes. In Fig. 5 we show the data collapse for the data. We determined as the crossing point of secondorder polynomial fits to the data, and used the value corresponding to the threedimensional orthogonal Anderson model nu_orth (). Although far from conclusive, this plot shows at least that agreement between the critical properties of the two models is possible.
5 Conclusions
In this paper we have discussed a possible mechanism leading to the localisation of lowlying modes of the Dirac operator in the hightemperature phase of QCD. Our proposal is that the localising “traps” are provided by spatial fluctuations (“islands”) of the local Polyakov lines away from the ordered value (i.e., the identity in colour space). In this paper we have generalised the original proposal of Ref. Bruckmann:2011cc () to , which has allowed us to identify the relevant variables for localisation, namely the phases of the eigenvalues of the local Polyakov lines, which govern the “depth” of the localising “traps”. These phases enter the Dirac equation in the form of effective boundary conditions for the quark eigenmodes, and provide a threedimensional source of disorder. The threedimensional nature of the “islands”, which constitute disconnected regions in a “sea” of (almost) ordered Polyakov lines, sheds light on the compatibility of the critical properties of the Dirac spectrum at the localisation/delocalisation transition with those of the threedimensional unitary Anderson model crit ().
To further substantiate our proposal, and check the viability of the “sea/islands” explanation, we have reproduced the qualitative features of the localised modes in a threedimensional effective model, designed to produce localisation precisely through the proposed mechanism. This “IsingAnderson” model is built replacing the fourdimensional (anisotropic) lattice, appropriate for finitetemperature QCD, with a threedimensional one, and replacing the temporal covariant derivative in the Dirac operator with a diagonal (onsite) noise. The dynamics of the diagonal noise is governed by an Ising model with continuous spins in the ordered phase, in order to reproduce the qualitative features of the Polyakov line configurations. Localised modes are indeed present, and respond to changes in the parameters (i.e., the temperature of the Ising system and the strength of the spinfermion coupling) as expected theoretically, thus reinforcing our confidence in the validity of the “sea/islands” explanation in QCD.
There are several possible extensions of the present work. Besides a thorough check of the critical properties of the model at the localisation/delocalisation transition, to be compared with those of the Anderson model in the appropriate symmetry class, the most interesting issue to study is how to properly model the effect of the spatial hopping. As a matter of fact, the diagonal disorder which we believe is responsible for localisation in the hightemperature phase, i.e., the effective boundary conditions, is present also in the lowtemperature phase of QCD, but in that case it is ineffective in producing localisation. Indeed, in the lowtemperature phase the system looks “fully fourdimensional”, and is basically insensitive to the temporal boundary conditions. In our opinion, in order to understand why this happens, and how the transition between the two phases is realised, adopting the point of view of “QCD as a randommatrix model” held in the present paper, it is necessary to understand how the spatial hopping acts in the two phases, especially concerning its capability of producing delocalised states around the origin. This could lead to useful insights in the subject of chiral symmetry breaking in QCD.
Acknowledgements
MG and TGK are supported by the Hungarian Academy of Sciences under “Lendület” grant No. LP2011011. FP is supported by OTKA under the grant OTKANF104034.
Appendix A Symmetry class of the effective model
The Hamiltonian of our toy model is of the form
(20) 
with , , and . In this Appendix Latin indices run from 1 to 3, Greek indices from 1 to 4. We want to find an antiunitary “timereversal” symmetry of this system. We denote with the corresponding operator, which we parameterise as , with the usual complex conjugation and unitary to be determined, and such that , with , and . We have to impose
(21) 
Since the last relation has to hold for any choice of , we see that
(22) 
by setting , and also
(23) 
by setting . must be of the general form
(24) 
where . Imposing Eq. (23) restricts the form of to
(25) 
Plugging this into the first equation in Eq. (22) we obtain
(26) 
which implies
(27) 
Using now the second equation in Eq. (22) we obtain, for the nonvanishing entries with ,
(28) 
This imposes the following equations,
(29)  
which are solved by
(30)  
However, since must respect the periodic boundary conditions, for lattices of odd linear size the only possibility is . Computing the square of the timereversal operator,
(31) 
so for odd lattices the symmetry class is the symplectic one. For even lattices, it is possible to choose for any , which yields
(32) 
so that the symmetry class is the orthogonal one (the existence of at least one antiunitary symmetry operator with square one is sufficient to see this). Even more simply, after spin diagonalisation the Hamiltonian reads
(33) 
and one can easily identify the timereversal operator , which has .
References
 (1) T. Banks and A. Casher, Nucl. Phys. B 169 (1980) 103.
 (2) J. J. M. Verbaarschot and T. Wettig, Ann. Rev. Nucl. Part. Sci. 50 (2000) 343 [hepph/0003017].
 (3) P. de Forcrand, AIP Conf. Proc. 892 (2007) 29 [heplat/0611034].
 (4) T. G. Kovács and F. Pittler, Phys. Rev. D 86 (2012) 114515 [arXiv:1208.3475 [heplat]].
 (5) M. Giordano, T. G. Kovács and F. Pittler, PoS LATTICE 2013 (2013) 212 [arXiv:1311.1770 [heplat]].
 (6) A. M. GarcíaGarcía and J. C. Osborn, Nucl. Phys. A 770 (2006) 141 [heplat/0512025].
 (7) A. M. GarcíaGarcía and J. C. Osborn, Phys. Rev. D 75 (2007) 034503 [heplat/0611019].
 (8) T. G. Kovács, Phys. Rev. Lett. 104 (2010) 031601 [arXiv:0906.5373 [heplat]].
 (9) T. G. Kovács and F. Pittler, Phys. Rev. Lett. 105 (2010) 192001 [arXiv:1006.1205 [heplat]].
 (10) Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabó, Jour. High Energy Phys. 01 (2006) 089 [heplat/0510084].
 (11) S. Borsányi, G. Endrődi, Z. Fodor, A. Jakovác, S. D. Katz, S. Krieg, C. Ratti and K. K. Szabó, Jour. High Energy Phys. 11 (2010) 077 [arXiv:1007.2580 [heplat]].
 (12) G. Cossu et al. [JLQCD Collaboration], arXiv:1412.5703 [heplat].
 (13) M. Giordano, T. G. Kovács, S. D. Katz and F. Pittler, arXiv:1410.8392 [heplat].
 (14) P. de Forcrand and O. Philipsen, Jour. High Energy Phys. 11 (2008) 012 [arXiv:0808.1096 [heplat]].
 (15) M. Giordano, T. G. Kovács and F. Pittler, Phys. Rev. Lett. 112 (2014) 102002 [arXiv:1312.1179 [heplat]].
 (16) P. W. Anderson, Phys. Rev. 109 (1958) 1492.
 (17) P. A. Lee and T. V. Ramakrishnan, Rev. Mod. Phys. 57 (1985) 287.
 (18) F. Evers and A. D. Mirlin, Rev. Mod. Phys. 80 (2008) 1355 [arXiv:0707.4378 [condmat.meshall]].
 (19) M. Mehta, Random Matrices (Academic Press, San Diego, 1991).
 (20) K. Slevin and T. Ohtsuki, Phys. Rev. Lett. 78 (1997) 4083 [condmat/9704192 [condmat.disnn]].
 (21) E. N. Economou, P. D. Antoniou, Solid State Commun. 21 (1977) 285.
 (22) D. Weaire and V. Srivastava, Solid State Commun. 23 (1977) 863.
 (23) F. Bruckmann, T. G. Kovács and S. Schierenberg, Phys. Rev. D 84 (2011) 034505 [arXiv:1105.5336 [heplat]].
 (24) A. Roberge and N. Weiss, Nucl. Phys. B 275 (1986) 734.
 (25) P. de Forcrand and O. Philipsen, Nucl. Phys. B 642 (2002) 290 [heplat/0205016].
 (26) M. D’Elia and M. P. Lombardo, Phys. Rev. D 67 (2003) 014505 [heplat/0209146].
 (27) A. M. GarcíaGarcía and E. Cuevas, Phys. Rev. B 74 (2006) 113101 [condmat/0602331 [condmat.disnn]].
 (28) L. G. Yaffe and B. Svetitsky, Phys. Rev. D 26 (1982) 963.
 (29) T. A. DeGrand and C. E. DeTar, Nucl. Phys. B 225 (1983) 590.
 (30) K. Slevin and T. Ohtsuki, Phys. Rev. Lett. 82 (1999) 382 [condmat/9812065 [condmat.disnn]].