2 Background
Disorder and Localisation
Disorder is a fact of life for the condensed matter physicist. No sample will ever be completely free of contamination or of structural defects. The classical Drude theory of electron conductivity envisages electrons as scattering off impurities. In this model, one would expect the electrical conductivity to be proportional to the mean free path [1], decreasing smoothly as the number of defects increases. However, Anderson in 1958 [2] showed that in a simple model, there is some critical level of disorder at which all single particle eigenstates localise.
What would later be known as Anderson localisation is characterised by exponentially localised eigenfunctions \(\psi(x) \sim e^{-x/\lambda}\) which cannot contribute to transport processes. The localisation length \(\lambda\) is the typical scale of localised states and can be extracted with transmission matrix methods [3]. Anderson localisation provided a different kind of insulator to that of the band insulator.
The Anderson model is about the simplest model of disorder one could imagine \[ H = -t\sum_{\langle jk \rangle} c^\dagger_j c_k + \sum_j V_j c_j^\dagger c_j. \qquad{(1)}\]
It is one of non-interacting fermions subject to a disorder potential \(V_j\) drawn uniformly from the interval \([-W,W]\). The discovery of localisation in quantum systems was surprising at the time given the seeming ubiquity of extended Bloch states. Within the Anderson model, all the states localise at the same disorder strength \(W\). Later, Mott showed that in other contexts extended Bloch states and localised states can coexist at the same disorder strength but at different energies. The transition in energy between localised and extended states is known as a mobility edge [4].
Localisation phenomena are strongly dimension dependent. In 3D the scaling theory of localisation [5,6] shows that Anderson localisation is a critical phenomenon with critical exponents both for how the conductivity vanishes with energy when approaching the mobility edge and for how the localisation length increases below it. By contrast, in 1D disorder generally dominates. Even the weakest disorder exponentially localises all single particle eigenstates in the 1D Anderson model. Only long-range spatial correlations of the disorder potential can induce delocalisation [7–12].
Later localisation was found in disordered interacting many-body systems:
\[ H = -t\sum_{\langle jk \rangle} c^\dagger_j c_k + \sum_j V_j c_j^\dagger c_j + U\sum_{jk} n_j n_k. \] Here, in contrast to the Anderson model, localisation phenomena are robust to weak perturbations of the Hamiltonian. This is called many-body localisation [13,14].
Both many-body localisation and Anderson localisation depend crucially on the presence of quenched disorder. Quenched disorder takes the form a static background field drawn from an arbitrary probability distribution to which the model is coupled. Disorder may also be introduced into the initial state of the system rather than the Hamiltonian. This has led to ongoing interest in the possibility of disorder-free localisation where the disorder is instead annealed. In this scenario, the disorder necessary to generate localisation is generated entirely from the thermal fluctuations of the model.
The concept of disorder-free localisation was first proposed in the context of Helium mixtures [15] and then extended to heavy-light mixtures in which multiple species with large mass ratios interact. The idea is that the heavier particles act as an effective disorder potential for the lighter ones, inducing localisation. Two such models [16,17] instead find that the models thermalise exponentially slowly in system size, which Ref. [16] dubs Quasi-MBL.
True disorder-free localisation does occur in exactly solvable models with extensively many conserved quantities [18]. As conserved quantities have no time dynamics, this can be thought of as taking the separation of timescales to the infinite limit. The localisation phenomena present in the Falicov-Kimball model are instead the result of annealed disorder. A strong separation of timescales means that the heavy species is approximated as immobile with respect to the lighter itinerant species. At finite temperature the heavy species acts as a disorder potential for the lighter one. However, in contrast to quenched disorder, the probability distribution of annealed disorder is entirely determined by the thermodynamics of the Hamiltonian. In the 2D FK model this leads to multiple phases where localisation effects are relevant. At low temperatures, the heavy species orders to a symmetry broken CDW phase, leading to a traditional band gap insulator. At higher temperatures, however, thermal disorder causes the light species to localise. At weak coupling, the localisation length can be very large, so finite sized systems may still conduct, an effect known as weak localisation [19].
In Chapter 3 we will consider a generalised FK model in 1D and study how the disorder generated near a 1D thermodynamic phase transition interacts with localisation physics.
Topological Disorder
So far we have considered disorder as a static or dynamic field coupled to a model defined on a translation invariant lattice. Another kind of disordered system that is worthy of study are amorphous systems. Amorphous systems have disordered bond connectivity, so called topological disorder. As discussed in the introduction these include amorphous semiconductors such as amorphous germanium and silicon [20–23]. While materials do not have long-range lattice structure they can enforce local constraints such as the approximate coordination number \(z = 4\) of silicon.
Topological disorder can be qualitatively different from other disordered systems. Disordered graphs are constrained by fixed coordination number and the Euler equation. A standard method for generating such graphs with coordination number \(d+1\) is Voronoi tessellation [24,25]. The Harris [26] and the Imry-Mar [27] criteria are key results on the effect of disorder on thermodynamic phase transitions. The Harris criterion signals when disorder will affect the universality of a thermodynamic critical point while the Imry-Ma criterion simply forbids the formation of long-range ordered states in \(d \leq 2\) dimensions in the presence of disorder. Both these criteria are modified for the case of topological disorder. This is because the Euler equation and vertex degree constraints lead to strong anti-correlations which mean that topological disorder is effectively weaker than standard disorder in 2D [28,29]. This does not apply to 3D Voronoi lattices where the Euler equation contains an extra volume term and so is effectively a weaker constraint. It is worth exploring how QSLs and disorder interact. The KH model has been studied subject to both flux [30] and bond [31] disorder. In some instances it seems that disorder can even promote the formation of a QSL ground state [32]. It has also been shown that the KH model exhibits disorder-free localisation after a quantum quench [33].
In chapter 4 we will put the Kitaev model onto 2D Voronoi lattices and show that much of the rich character of the model is preserved despite the lack of long-range order.
Diagnosing Localisation in practice
Looking at practical tools for diagnosing localisation, there are a few standard methods [6].
The most direct method would be to fit a function of the form \(\psi(x) = f(x) e^{-|x-x_0|/\lambda}\) to each single particle wavefunction to extract the localisation length \(\lambda\). This method is little used in practice since it requires storing and processing full wavefunctions which quickly becomes computationally expensive for large systems.
For low dimensional systems with quenched disorder, transfer matrix methods can be used to directly extract the localisation length. These work by turning the time independent Schrödinger equation \(\hat{H}|\psi\rangle = E|\psi\rangle\) into a matrix equation linking the amplitude of \(\psi\) on each \(d-1\) dimensional slice of the system to the next and looking at average properties of this transmission matrix. This method is less useful for systems like the FK model where the disorder as a whole must be sampled from the thermodynamic ensemble.
A more versatile method is based on the Inverse Participation Ratio (IPR). The IPR is defined for a normalised wave function \(\psi_i = \psi(x_i), \sum_i |\psi_i|^2 = 1\) as its fourth moment [6]:
\[ P^{-1} = \sum_i |\psi_i|^4. \]
The name derives from the fact that this operator acts as a measure of the volume where the wavefunction is significantly different from zero. They can alternatively be thought of as providing a measure of the average diameter \(R\) from \(R = P^{1/d}\). See fig. 1 for the distinction between \(R\) and \(\lambda\).
For localised states, the inverse participation ratio \(P^{-1}\) is independent of system size while for plane wave states in \(d\) dimensions \(P^{-1} = L^{-d}\). States may also be intermediate between localised and extended, described by their fractal dimensionality \(d > d* > 0\):
\[ P(L)^{-1} \sim L^{-d*}. \]
Such intermediate states tend to appear as critical phenomena near mobility edges [34]. For finite size systems, these relations only hold once the system size \(L\) is much greater than the localisation length. When the localisation length is comparable to the system size, the states still contribute to transport, this is the aforementioned weak localisation effect [35,36].
In the following two chapters I will use an energy resolved IPR \[ \begin{aligned} DOS(\omega) &= \sum_n \delta(\omega - \epsilon_n)\\ IPR(\omega) &= DOS(\omega)^{-1} \sum_{n,i} \delta(\omega - \epsilon_n) |\psi_{n,i}|^4, \end{aligned} \] where \(\psi_{n,i}\) is the wavefunction corresponding to the energy \(\epsilon_n\) at the ith site. In practice I bin the IPRs into finely spaced bins in energy space and use the mean IPR within each bin.
Chapter Summary
In this chapter we have covered the Falicov-Kimball model, the Kitaev Honeycomb model and the theory of disorder and localisation. We saw that the FK model is one of immobile species (spins) interacting with an itinerant quantum species (electrons). While the KH model is specified in terms of spins on a honeycomb lattice interacting via a highly anisotropic Ising coupling, it can be transformed into one of Majorana fermions interacting with a classical gauge field that supports immobile flux excitations. In each case, it is the immobile species that make the model exactly solvable. Both models have rich ground states and thermodynamic phase diagrams. The last part of this chapter dealt with disorder and how it almost inevitably leads to localisation. Both the FK and KH models are effectively disordered at finite temperatures by their immobile species. In the next chapter we will look at a version of the FK model in 1D augmented with long-range interactions in order to retain its ordered phase. The model is translation invariant but we will see that it exhibits disorder-free localisation. After that we will look at the KH model defined on an amorphous lattice with vertex degree \(z=3\).
Next Chapter: 3 The Long Range Falicov-Kimball Model