A picture of me.

Tom Hodson

Physicist, Programmer, Maker and Baker

This thesis has focussed on two strongly correlated systems. In these systems the many-body ground state can be complex and often cannot be reduced to or even adiabatically connected to a product state. I looked at the Falicov-Kimball (FK) model and the Kitaev Honeycomb (KH) model and defined extensions to them: the Long-Range Falicov-Kimball (LRFK) model and the Amorphous Kitaev (AK) model.

These models are all exactly solvable. They contain extensively many conserved charges which allow their Hamiltonians, and crucially, the interaction terms within them, to be written in quadratic form. This allows them to be solved using the theoretical machinery of non-interacting systems. In the case of the FK and LRFK models, this solvability arises from what is essentially a separation of timescales. The heavy particles move so slowly that they can be treated as stationary. In the KH and AK models, on the other hand, the origin of the conserved degrees of freedom is more complex. Here, the algebra of the Pauli matrices interacts with the trivalent lattices on which the models are defined, to give rise to an emergent \(\mathbb{Z}_2\) gauge field whose fluxes are conserved. This latter case is a beautiful example of emergence at play in condensed matter. The gauge and Majorana physics of the KH and AK models seems to arise spontaneously from nothing. Of course, this physics was hidden within the structure and local symmetries of the spin Hamiltonian all along.

At first glance, exactly solvable models can seem a little too fine tuned to be particularly relevant to the real world. Surely, these models don’t spontaneously arise in nature? The models studied here provide two different ways to answer this. As we saw in chapter 2, the FK model arises quite naturally as a limit of the Hubbard model. The Hubbard model is not exactly solvable. In fact, the FK model has been used as a way to understand more about the behaviour of the Hubbard model itself and of the Mott insulating state. We have also seen that it can provide insight into other phenomena such as disorder-free localisation. The KH model was not originally proposed as a model of any particular physical system. It was nevertheless a plausible microscopic Hamiltonian and, given its remarkable properties, it is little wonder that material candidates for Kitaev physics were quickly found. In neither case is the model expected to be a perfect description of any material. Indeed, more realistic corrections to each model are likely to break their integrability. Despite this, exactly solvable models, by virtue of being solvable, can provide important insights into the diverse physics of strongly correlated materials.

In chapter 3, we looked at a generalised FK model in 1D: the LRFK model. Metal-insulator transitions are a key theme of work on the FK model and our 1D extension to it was no exception. With the addition of long-range interactions, the model showed a similarly rich phase diagram as its higher dimensional cousins, allowing us to look at transitions between metallic states and insulators of the band, Anderson and Mott variety. We also looked at thermodynamics in 1D and how thermal fluctuations of the conserved charges can lead to disorder-free localisation in the FK and LRFK models. The initial surprising results that suggested the presence of a mobility edge in 1D turned out to be a weak localisation effect present in finite sized systems.

That the Mott insulating state is a key part of the work on the LRFK model is fitting because Mott insulators are the main route to the formation of Quantum Spin Liquid (QSL) states. QSLs are themselves a primary driver of interest in the KH model and our extension, the AK model. In chapter 4, we addressed the question of whether frustrated magnetic interactions on amorphous lattices can give rise to quantum phases such as the QSL state and found that indeed they can. The AK model, a generalisation of the KH model to random lattices with fixed coordination number three, supports a kind of symmetry broken QSL state called a chiral spin liquid. We showed numerically that the ground state of the model follows a simple generalisation of Lieb’s theorem  [13]. As with other extensions of the KH model  [4], we found that removing the chiral symmetry of the lattice allows the model to support a gapped phase with non-Abelian anyon excitations. The broken chiral symmetry plays the role of the external magnetic field in the original KH model. Finally, like the KH model, finite temperature causes vortex defects to proliferate. This triggers a transition to a thermal metal state.

Figure 1: An example of a Metal Organic Framework (MOF). MOFs are large synthetic organic ligands coordinated with metal ions that can cross link to form crystals. It has been demonstrated that these materials can sustain Kitaev-like interactions  [5]. They can also be made into amorphous materials  [6]. This image from  [7] is in the public domain.

Unlike the 2D FK model and 1D LRFK models, the KH and AK models don’t have a finite-temperature phase transition (FTPT). They immediately disorder at any finite temperature  [3]. However, generalisations of the KH model to 3D do in general have an FTPT. Indeed, the role of dimensionality has been a key theme in this work. Both localisation and thermodynamic phenomena depend crucially on dimensionality, with thermodynamic order generally suppressed and localisation effects strengthened in low dimensions. The graph theory that underpins the KH and AK models itself also changes strongly with dimension. Voronisation in 2D produces trivalent lattices, on which the spin-\(1/2\) AK model is exactly solvable. Meanwhile in 3D, voronisation gives us \(z=4\) lattices upon which a spin-\(3/2\) generalisation to the KH model is exactly solvable  [810]. Similarly, planar and toroidal graphs are a uniquely 2D construct. Satisfying planarity imposes constraints on the connectivity of planar graphs. This leads amorphous planar graphs to have strong anti-correlations which can violate otherwise robust bounds like the Harris criterion  [11]. Contrast this with Anderson localisation in 1D where only longer range correlations in the disorder can produce surprising effects  [1217].

Looking towards future work, the LRFK model provides multiple possible routes. One interesting idea is to park the model at a thermal critical point. This would generate a scale-free disorder potential, which could potentially lead to complex localisation physics not often seen in 1D  [12,18]. Like other solvable models of disorder-free localisation, the LRFK model should also exhibit intriguing out-of-equilibrium physics, such as slow entanglement dynamics. This could be used to help understand these phenomena in more generic interacting systems  [19]. There is also the rich ground state phenomenology of the FK model as a function of filling  [20], such as the devil’s staircase  [21] as well as superconductor like states  [22]. Could the LRFK model stabilise these at finite temperature? Finally, a topological variant of the LRFK model akin to the Su-Schrieffer-Heeger (SSH) model could be an interesting way to probe the interplay of topological bound states and thermal domain wall defects.

Looking at the AK model, we discussed whether its physics might be realisable in amorphous versions of known KH candidate materials  [23]. Alternatively, we might be able to engineer them in synthetic materials, such as Metal Organic Frameworks (MOFs). Work on MOFs has already explored the possibility of both Kitaev-like interactions and amorphous lattices  [5,6]. It is an open question whether the superexchange couplings that generate Kitaev interactions could survive the transition to an amorphous lattice. If the interactions do survive, there will likely be many defects of different kinds present in the resulting material. These might take the form of dangling bonds, vertex degree \(> 3\) or violations of the colouring conditions. I therefore hope that future work examines how robust the QSL and CSL states of the AK model are to these kinds of disorder. This is a difficult task, as many of these classes of defects would break the integrability of the AK model that we relied on to make the work computationally feasible  [2428]. While we are considering models with defects, we might consider alternate lattice construction techniques  [29,30]. Lattices constructed using these methods would have defects but may have other desirable properties when compared to Voronoi lattices.

In terms of experimental signatures, we discussed the quantised thermal Hall effect  [3134], local probes such as spin-polarised scanning tunnelling microscopy  [3537], and longitudinal heat transport signatures  [38]. One possible difficulty is that the introduction of topological disorder may dilute some of these signatures. On the brighter side, topological disorder may also suppress competing interactions that would otherwise induce magnetic ordering. This could potentially widen the class of materials that could host a QSL or CSL ground state.

Both the LRFK and AK models could be interesting targets for quantum simulation. Cold atom experiments can naturally generate long-range interactions which would be well suited to simulations of the LRFK model  [39]. Optical lattice experiments are capable of supporting both Kitaev interactions  [40,41] and amorphous lattices  [42] so it is natural to ask whether they could simulate the AK model too.

Xiao-Gang Wen makes the case that one of the primary reasons to study QSLs is as a stepping stone to understanding the high-\(T_c\) superconductors  [43]. His logic is that since the high-\(T_c\) superconductors are believed to arise from doped Mott insulators, the QSLs, which arise from undoped Mott insulators, make a good jumping off point. This is where exactly solvable models like the FK and KH model shine. The FK model provides a tractable means to study superconductor like states in doped Mott insulators  [22], while the KH model gives us a tangible QSL state to play with. The extensions introduced in this work serve to explore how far we can push these models. To what extent they are a useful theoretical tool for understanding the behaviour of generic strongly correlated materials remains to be seen. Either way, I hope that this works spurs more interest in amorphous quantum many-body phases and strongly correlated materials.


E. H. Lieb, Flux Phase of the Half-Filled Band, Physical Review Letters 73, 2158 (1994).
K. O’Brien, M. Hermanns, and S. Trebst, Classification of Gapless Z₂ Spin Liquids in Three-Dimensional Kitaev Models, Phys. Rev. B 93, 085101 (2016).
T. Eschmann, P. A. Mishchenko, K. O’Brien, T. A. Bojesen, Y. Kato, M. Hermanns, Y. Motome, and S. Trebst, Thermodynamic Classification of Three-Dimensional Kitaev Spin Liquids, Phys. Rev. B 102, 075125 (2020).
H. Yao and S. A. Kivelson, An Exact Chiral Spin Liquid with Non-Abelian Anyons, Phys. Rev. Lett. 99, 247203 (2007).
M. G. Yamada, H. Fujita, and M. Oshikawa, Designing Kitaev Spin Liquids in Metal-Organic Frameworks, Phys. Rev. Lett. 119, 057202 (2017).
T. D. Bennett and A. K. Cheetham, Amorphous Metal–Organic Frameworks, Acc. Chem. Res. 47, 1555 (2014).
N. L. Rosi, J. Eckert, M. Eddaoudi, D. T. Vodak, J. Kim, M. O’Keeffe, and O. M. Yaghi, Hydrogen Storage in Microporous Metal-Organic Frameworks, Science 300, 1127 (2003).
H. Yao, S.-C. Zhang, and S. A. Kivelson, Algebraic Spin Liquid in an Exactly Solvable Spin Model, Phys. Rev. Lett. 102, 217202 (2009).
S. Ryu, Three-Dimensional Topological Phase on the Diamond Lattice, Phys. Rev. B 79, 075124 (2009).
A. B. Harris, Effect of Random Defects on the Critical Behaviour of Ising Models, J. Phys. C: Solid State Phys. 7, 1671 (1974).
S. Aubry and G. André, Analyticity Breaking and Anderson Localization in Incommensurate Lattices, Proceedings, VIII International Colloquium on Group-Theoretical Methods in Physics 3, 18 (1980).
D. H. Dunlap, H.-L. Wu, and P. W. Phillips, Absence of Localization in a Random-Dimer Model, Phys. Rev. Lett. 65, 88 (1990).
F. M. Izrailev and A. A. Krokhin, Localization and the Mobility Edge in One-Dimensional Potentials with Correlated Disorder, Phys. Rev. Lett. 82, 4062 (1999).
A. Croy, P. Cain, and M. Schreiber, Anderson Localization in 1d Systems with Correlated Disorder, Eur. Phys. J. B 82, 107 (2011).
F. M. Izrailev, A. A. Krokhin, and N. M. Makarov, Anomalous Localization in Low-Dimensional Systems with Correlated Disorder, Physics Reports 512, 125 (2012).
T. Shimasaki, M. Prichard, H. E. Kondakci, J. Pagett, Y. Bai, P. Dotti, A. Cao, T.-C. Lu, T. Grover, and D. M. Weld, Anomalous Localization and Multifractality in a Kicked Quasicrystal, arXiv:2203.09442.
O. Hart, S. Gopalakrishnan, and C. Castelnovo, Logarithmic Entanglement Growth from Disorder-Free Localization in the Two-Leg Compass Ladder, Phys. Rev. Lett. 126, 227202 (2021).
C. Gruber, J. Iwanski, J. Jedrzejewski, and P. Lemberger, Ground States of the Spinless Falicov-Kimball Model, Phys. Rev. B 41, 2198 (1990).
C. Micheletti, A. B. Harris, and J. M. Yeomans, A Complete Devil’s Staircase in the Falicov-Kimball Model, J. Phys. A: Math. Gen. 30, L711 (1997).
P. Cai, W. Ruan, Y. Peng, C. Ye, X. Li, Z. Hao, X. Zhou, D.-H. Lee, and Y. Wang, Visualizing the Evolution from the Mott Insulator to a Charge-Ordered Insulator in Lightly Doped Cuprates, Nature Phys 12, 11 (2016).
S. Trebst and C. Hickey, Kitaev Materials, Physics Reports 950, 1 (2022).
J. G. Rau, E. K.-H. Lee, and H.-Y. Kee, Generic Spin Model for the Honeycomb Iridates Beyond the Kitaev Limit, Phys. Rev. Lett. 112, 077204 (2014).
J. Chaloupka, G. Jackeli, and G. Khaliullin, Kitaev-Heisenberg Model on a Honeycomb Lattice: Possible Exotic Phases in Iridium Oxides A₂IrO₃, Phys. Rev. Lett. 105, 027204 (2010).
J. Chaloupka, G. Jackeli, and G. Khaliullin, Zigzag Magnetic Order in the Iridium Oxide Na₂IrO₃, Phys. Rev. Lett. 110, 097204 (2013).
S. M. Winter, Y. Li, H. O. Jeschke, and R. Valentí, Challenges in Design of Kitaev Materials: Magnetic Interactions from Competing Energy Scales, Phys. Rev. B 93, 214431 (2016).
A. Agarwala, Topological Insulators in Amorphous Systems, in Excursions in Ill-Condensed Quantum Matter: From Amorphous Topological Insulators to Fractional Spins, edited by A. Agarwala (Springer International Publishing, Cham, 2019), pp. 61–79.
N. H. Christ, R. Friedberg, and T. D. Lee, Random Lattice Field Theory: General Formulation, Nuclear Physics B 202, 89 (1982).
Y. Kasahara, T. Ohnishi, Y. Mizukami, et al., Majorana Quantization and Half-Integer Thermal Quantum Hall Effect in a Kitaev Spin Liquid, Nature 559, 227 (2018).
T. Yokoi, S. Ma, Y. Kasahara, et al., Majorana Quantization and Half-Integer Thermal Quantum Hall Effect in a Kitaev Spin Liquid, Science 373, 568 (2021).
M. Yamashita, J. Gouchi, Y. Uwatoko, N. Kurita, and H. Tanaka, Sample Dependence of Half-Integer Quantized Thermal Hall Effect in the Kitaev Spin-Liquid Candidate α-RuCl₃, Phys. Rev. B 102, 220404 (2020).
J. A. N. Bruin, R. R. Claus, Y. Matsumoto, N. Kurita, H. Tanaka, and H. Takagi, Majorana Quantization and Half-Integer Thermal Quantum Hall Effect in a Kitaev Spin Liquid, Nature Physics 18, 401 (2022).
J. Feldmeier, W. Natori, M. Knap, and J. Knolle, Local Probes for Charge-Neutral Edge States in Two-Dimensional Quantum Magnets, Phys. Rev. B 102, 134423 (2020).
E. J. König, M. T. Randeria, and B. Jäck, Tunneling Spectroscopy of Quantum Spin Liquids, Phys. Rev. Lett. 125, 267206 (2020).
M. Udagawa, S. Takayoshi, and T. Oka, Scanning Tunneling Microscopy as a Single Majorana Detector of Kitaev’s Chiral Spin Liquid, Phys. Rev. Lett. 126, 127201 (2021).
C. W. J. Beenakker, Search for Majorana Fermions in Superconductors, Annual Review of Condensed Matter Physics 4, 113 (2013).
M. Lepers and O. Dulieu, Long-Range Interactions Between Ultracold Atoms and Molecules, arXiv:1703.02833.
L.-M. Duan, E. Demler, and M. D. Lukin, Controlling Spin Exchange Interactions of Ultracold Atoms in Optical Lattices, Phys. Rev. Lett. 91, 090402 (2003).
A. Micheli, G. K. Brennen, and P. Zoller, A Toolbox for Lattice-Spin Models with Polar Molecules, Nature Phys 2, 5 (2006).
X.-G. Wen, Quantum Orders and Symmetric Spin Liquids, Phys. Rev. B 65, 165113 (2002).