Frustrated magnetism: Thesis motivation

From water to spin ice

Frustration is ubiquitous in nature. So much so that the properties of water, including its unusually high melting and boiling points, its large heat capacity, latent heat, and anomalous density curve—all crucial for the development of life as we know it—are all consequences of frustration.

It is well known that the hydrogen atoms in \(\mathrm{H_2O}\) span an angle of about \(104^\circ\), slightly less than the \(109^\circ\) corresponding to the unperturbed tetrahedral coordination of hybridized \(sp^3\) orbitals of the central oxygen atom. This leads to a nonvanishing dipolar moment due to the partial charges induced by the larger electronegativity of oxygen as compared to hydrogen. In water ice (specifically the ice \(\mathrm{I}_h\) phase formed under terrestrial conditions), the molecules are arranged in a lattice with close to tetrahedral bond angles as well. Due to the partial charges, it is energetically favorable for exactly one hydrogen to sit between any neighboring oxygen atoms, being covalently bound to one of them and forming a weak “hydrogen” bond with the other. The orientation of the water molecules in the lattice is then determined by the Bernal-Fowler “ice rules” [1] which state that each molecule should be faced by the hydrogen atoms of two of its neighbors and the oxygen sides of the other two. An example of a configuration respecting these ice rules is shown in Fig. 1.

The ice rules do not uniquely determine the ground-state configuration of water ice. Rather, the number of compatible configurations, i.e. the degeneracy of the ground state, grows exponentially with system size. This then gives rise to an extensive contribution to the residual entropy which was famously estimated by Pauling as \(S_0/n = \ln(3/2) R\approx 3.4\,\mathrm{J}/(\mathrm{mol\,K})\) [2]. The extensive residual entropy consequently affects the thermodynamics of bulk properties and leads to the above anomalies associated with water. An exponential ground-state degeneracy is one of the hallmarks of frustrated systems. Somewhat sloppily, it is more commonly referred to as an extensive ground state degeneracy (exGSD) because it gives rise to an extensive residual entropy.

An analogous situation is realized in a class of magnetic systems which—for this reason—are known as spin ices. Experimental realizations of spin ices have been discovered in the rare-earth oxide materials \(\mathrm{Ho_2Ti_2O_7}\) [3] and \(\mathrm{Dy_2Ti_2O_7}\) [4,5]. In these compounds, only the rare-earth ions (\(\mathrm{Ho}^{3+}\) and \(\mathrm{Dy}^{3+}\)) are magnetic. These form a pyrochlore lattice, a network of corner-sharing tetrahedra; its geometry is depicted in Fig. 2. Due to a strong crystal-field anisotropy along the axis connecting the centers of tetrahedra, the magnetic moments behave as Ising spins. The energy of such a configuration is then minimized if within each of the tedrahedra, two of the spins point inwards, while the other two point outwards. This “two-in-two-out” rule is the spin-ice analogue of the ice rules in water. Indeed, Fig. 2 depicts a configuration obeying the two-in-two-out rule which is equivalent to the water ice configuration in Fig. 1 by having the spins in the former point towards the location of each of the hydrogen atoms in the latter. The residual entropy of \(\mathrm{Dy_2Ti_2O_7}\) has been found from its heat capacity to be consistent with that of water ice \(\mathrm{I}_h\) and Pauling’s estimate [6].

Gauge theoretical description and magnetic monopoles

The fluctuations permitted by the spin-ice rule may also be understood by analogy to magnetostatics. By identifying the spins \(\mathbf S_i\) at positions \(\mathbf r_i\) in the lattice with the local strength of an artificial magnetic field \(\mathbf B(\mathbf r_i)=\mathbf S_i\), the spin-ice rule can be summarized in a generalize Gauss’ law, \(\mathbf\nabla\cdot\mathbf B=0\) [7]. Hence, two magnetic “field lines” enter into each lattice tetrahedron and likewise two also exit it. After coarse graining, this immediately implies that the correlations between magnetic fields a distance \(\mathbf r\) apart take on a dipolar form [8,9,10],

and thus decays as \(1/\|\mathbf r\|^3\) at large distances. Power-law correlations like this are unusual away from critical regimes. They manifest themselves as sharp features called “pinch points” in the spin structure factor and have been observed in neutron scattering experiments on \(\mathrm{Ho_2Ti_2O_7}\) [11].

Violations of the spin-ice rule lead to defects for which \(\mathbf\nabla\cdot\mathbf B\) takes on a nonvanishing value which may be interpreted as a magnetic “charge” density. Excitations above the ground state manifold thus correspond to the creation of a pair of magnetic monopoles of opposite charge at which a magnetic field line—in this context referred to as a Dirac string—originates and terminates, respectively [7,12,13]. Once created, the monopoles can move away from each other at no additional energy cost. They thus constitute fractionalized excitations and the energy associated with each monopole is half that of a single spin flip [14]. Since the spins carry a magnetic moment, these monopoles exert a magnetostatic Coulomb interaction onto each other. This is weak compared to the electrostatic force but its signature was nonetheless first measured in diffuse neutron scattering experiments on \(\mathrm{Dy_2Ti_2O_7}\) in 2009 [15]. The ensuing onslaught of experiments has focussed primarily on the dynamics of monopoles both in classical [16] and quantum [17] spin ices. Monopoles have also been observed in artificial spin ice realized in lithographically fabricated magnetic nanodevices [18].

Classical and quantum spin liquids

Spin ices are the posterchild of a class of correlated spin states, nowadays collectively known as classical spin liquids (CSLs). More traditionally, these have also been called correlated or cooperative paramagnets as they crucially do not exhibit long-range order, but rather strong algebraically decaying correlations. As in the case of spin ice, these originate from fluctuations within an extensive ground-state manifold which is characterized by local constraints such as the spin-ice rule. The constraints in turn signal the emergence of a gauge theoretical structure, complete with monopole-like fractionalized excitations.

At low enough temperatures, even fluctuations within the ground-state manifold eventually slow down and the CSL may freeze out of equilibrium. For example, in spin ice, at the very least hexagonal rings of six spins can cooperatively flip at once while remaining in the ground state. The ever-smaller tunneling probabilities for such zero modes hence cause the spins to freeze into an amorphous glassy state; in \(\mathrm{Dy_2Ti_2O_7}\) this reportedly happens below \(T_f=0.5\,\mathrm{K}\) [19]. Other types of CSLs may instead enter into an ordered phase below a transition temperature \(T_c\) (see “Hidden order” below for a possible mechanism). Either way, the relevant temperatures are expected to lie significantly below the scale set by the Curie-Weiss temperature \(|\Theta_\text{CW}|\) which provides a rough estimate for the ordering temperature in unfrustrated systems. \(\Theta_\text{CW}\) can be determined from the asymptotic high-temperature behavior of the magnetic susceptibility, \(\chi^{-1}\sim (T - \Theta_\text{CW})\). The frustration-induced CSL behavior thus takes hold in the regime \(T_c\)<\(T\)<\(|\Theta_\text{CW}|\) and its extent is commonly quantified by the frustration parameter \(f=|\Theta_\text{CW}|/T_c\) (or \(f=|\Theta_\text{CW}|/T_f\) in case of freezing). While \(f\sim 5\) may otherwise still occur, frustration can easily result in \(f\gg 100\) [20].

Materials with small atomic spins might, however, display an entirely different behavior at the lowest temperatures: quantum fluctuations due to a nonvanishing zero-point energy may allow for the formation of a quantum spin liquid (QSL) which persists down to \(T=0\), corresponding to \(f=\infty\).

The energy between any two antiferromagnetically coupled spins is locally minimized by a singlet. One may naively assume that the ground state thus consists of a product of such valence bond singlets, periodically tiled over the bonds of the lattice. Such a valence bond solid (VBS) has a vanishing total spin and is hence nonmagnetic, yet by choosing a specific dimer covering, it breaks lattice symmetries. Anderson proposed that a system might rather form a superposition of all possible dimer coverings of spins into valence bonds on the lattice [21–23], thereby retaining the lattice symmetry. In the language of quantum chemistry, such a state would exhibit valence bonds which constantly fluctuate amongst each other, but do so in a phase-coherent (highly entangled) way, thus coining the name resonating valence bond (RVB) state.

The RVB state is the archetype of QSLs, but a superposition of dimer coverings is not the only possible microscopic realization of a QSL state and a consensus on the defining characteristics of a QSL was not reached until recently (and even that is debateable). Quite generally, however, a “true” QSL would be expected to host fractionalized excitations [20,24]. This makes them attractive for quantum computation as it allows for the development of topological error correction codes, as exemplified by Kitaev’s toric code [25,26]. Unfortunately, as of yet, not a single QSL has been experimentally identified in an unambiguous way, despite tremendous effort over the past decade. Extensive reviews of the state of the art of experiments into both CSLs and QSL candidates are provided in Refs. 24 and 27.

Hidden order

The experimental identification of spin liquids, both classical and quantum, typically has to resort to a negative definition, i.e. ruling out the presence of any symmetry-breaking order, in lieu of experimentally accessible genuine spin-liquid characteristics. This is further complicated by the fact that such orders in frustrated lattice geometries may too elude common experimental probes as they do not give rise to a macroscopically observable magnetization. For this reason, they are referred to as “hidden” orders [28–33].

Indeed, even in Monte Carlo simulations of classical spin models where one has access to the full microscopic spin configuration, identifying the presence of hidden orders remains a highly nontrivial task. This is due to their multipolar nature: in contrast to dipolar orders whose symmetry-breaking is characterized by a nonvanishing vector-like (rank-1) order parameter such as the (staggered) magnetization in the case of (anti)ferromagnetism, frustrated lattice geometries may favor orders which only break part of the full spin-rotational symmetry while staying disordered with respect to the remaining symmetry, giving rise to order parameter tensors of higher rank such as matrices (quadrupolar order) or third-rank tensors (octupolar order). Examples can be found in spin nematic orders [34–45]. Conventional numeric probes then often suffer the same fate as their experimental counterparts and one has to resort to (weak) signatures in thermodynamic properties such as the heat capacity to pin down the presence of a phase transition. Even then, this does not yield any insight into the nature of the tentative phase.

It may seem paradoxical at first that a spin system whose ground state manifold is defined by constraints, which in and of themselves do not break any symmetries, would develop an ordered phase at low temperatures. One mechanism to facilitate symmetry-breaking in classical systems was first described by Villain et al. [46]. While by definition all ground states have the same energy, some of them may be entropically selected for their lower free energy. This is the case when a certain subset of ground state configurations admits thermal fluctuations into a vast region of phase space surrounding these ground states at small energy expenditure by means of soft modes¹. The phenomenon was coined order-by-disorder [8] and famously accounts–among other cases [47]—for the selection of coplanar spin order in the Heisenberg antiferromagnet on the kagome lattice which was postulated by Chalker et al. in 1992 [48]. However, his initial description in terms of a quadrupolar order parameter was later recognized as incomplete and rather called for an additional octupolar order parameter [49–51].

Premise of my PhD thesis

The study of frustrated magnets, particularly the quest for spin liquids—both classical and quantum—remains an active field of research in both theory and experiment. While QSLs have proved more elusive, the occurrence and stability of CSLs, especially as compared to spin nematic phases, is far from understood. Advances in experimental techniques over the past decade have accelerated the pace of the synthesis and charaterization of new materials to the point where it has become a matter of months, not years. Meanwhile the process of piecing together the phase diagram of each new model remains a laborious and manual process which heavily depends on the ingenuity and intuition of individual researchers, rendering it a serious bottleneck.

Devising a machine learning scheme to aid this process, inspired by the recent sprawl of activity, thus seemed like a worthwhile endeavor. However, merely recognizing different phases in frustrated models—while useful—would not have sufficed; rather, a highly interpretable machine is in demand which allows for the reconstruction of the physical order parameter. For this reason, we turned to kernel methods and proposed a tensorial kernel (TK) function. As the name suggests, its functional form was chosen to be capable of expressing tensorial order parameters like those characteristic of multipolar order, while not imposing any further restrictions other than locality. Its design thus does not require the problem to be already “solved”, but covers the realm of multipolar tensor order parameters exhaustively. When combined with a support vector machine, it is jointly referred to as TK-SVM, but the use in conjunction with kernel principal component analysis (kPCA) has also been demonstrated.

In my thesis, the TK is applied to the infamous case of the aforementioned Heisenberg model on the kagome lattice as well as the XXZ model on the pyrochlore lattice which features a diverse phase diagram, including a spin nematic phase and three types of CSLs [52,53]. Indeed, the order parameter tensors of both the spin nematic phase in the case of the XXZ pyrochlore model and the complete characterization of the hidden order in the kagome model have been identified by TK-SVM.

Serendipitously, while the original design of the TK was aimed at expressing multipolar order, it turned out to be also capable of expressing the emergent constraints which characterize the ground-state manifold of CSLs and underpin their gauge-theoretical description. In that regard, the method went above and beyond expectations and reconstructed the phase diagram of the XXZ pyrochlore model, including the location of the crossover regimes between different types of CSLs, with little human guidance.

References

• [1] J. D. Bernal and R. H. Fowler. A Theory of Water and Ionic Solution, with Particular Reference to Hydrogen and Hydroxyl Ions. J. Chem. Phys. 1, 515–548 (1933).
• [2] Linus Pauling. The Structure and Entropy of Ice and of Other Crystals with Some Randomness of Atomic Arrangement. J. Am. Chem. Soc. 57, 2680–2684 (1935).
• [3] M. J. Harris, S. T. Bramwell, D. F. McMorrow, T. Zeiske, and K. W. Godfrey. Geometrical Frustration in the Ferromagnetic Pyrochlore Ho₂Ti₂O₇. Phys. Rev. Lett. 79, 2554–2557 (1997).
• [4] Steven T. Bramwell and Michel J. P. Gingras. Spin Ice State in Frustrated Magnetic Pyrochlore Materials. Science 294, 1495–1501 (2001).
• [5] Jason S. Gardner, Michel J. P. Gingras, and John E. Greedan. Magnetic pyrochlore oxides. Rev. Mod. Phys. 82, 53–107 (2010).
• [6] A. P. Ramirez, A. Hayashi, R. J. Cava, R. Siddharthan, and B. S. Shastry. Zero-point entropy in ‘spin ice’. Nature 399, 333–335 ((1999).).
• [7] Claudio Castelnovo, Roderich Moessner, and Shivaji L. Sondhi. Magnetic monopoles in spin ice. Nature 451, 42–45 (2008).
• [8] John T. Chalker. Geometrically Frustrated Antiferromagnets: Statistical Mechanics and Dynamics. In Claudine Lacroix, Philippe Mendels, and Frédéric Mila, editors, Introduction to Frustrated Magnetism, chapter 1. Springer Berlin Heidelberg (2011). ISBN 978-3-642-10588-3.
• [9] Christopher L. Henley. Power-law spin correlations in pyrochlore antiferromagnets. Phys. Rev. B 71, 014424 (2005).
• [10] Christopher L. Henley. The “Coulomb Phase” in Frustrated Systems. Annu. Rev. Condens. Matter Phys. 1, 179–210 (2010).
• [11] T. Fennell, P. P. Deen, A. R. Wildes, K. Schmalzl, D. Prabhakaran, A. T. Boothroyd, R. J. Aldus, D. F. McMorrow, and S. T. Bramwell. Magnetic Coulomb Phase in the Spin Ice Ho₂Ti₂O₇. Science 326, 415–417 (2009).
• [12] Ivan A. Ryzhkin. Magnetic relaxation in rare-earth oxide pyrochlores. J. Exp. Theor. Phys. 101, 481–486 (2005).
• [13] Claudio Castelnovo, Roderich Moessner, and Shivaji L. Sondhi. Spin Ice, Fractionalization, and Topological Order. Annu. Rev. Condens. Matter Phys. 3, 35–55 (2012).
• [14] Ludovic D. C. Jaubert and Peter C. W. Holdsworth. Signature of magnetic monopole and Dirac string dynamics in spin ice. Nat. Phys. 5, 258–261 (2009).
• [15] D. J. P. Morris, D. A. Tennant, S. A. Grigera, B. Klemke, C. Castelnovo, R. Moessner, C. Czternasty, M. Meissner, K. C. Rule, J.-U. Hoffmann, K. Kiefer, S. Gerischer, D. Slobinsky, and R. S. Perry. Dirac Strings and Magnetic Monopoles in the Spin Ice Dy₂Ti₂O₇. Science 326, 411–414 (2009).
• [16] Sean R. Giblin, Steven T. Bramwell, Peter C. W. Holdsworth, Dharmalingam Prabhakaran, and Ian Terry. Creation and measurement of long-lived magnetic monopole currents in spin ice. Nat. Phys. 7, 252–258 (2011).
• [17] Bruno Tomasello, Claudio Castelnovo, Roderich Moessner, and Jorge Quintanilla. Correlated Quantum Tunneling of Monopoles in Spin Ice. Phys. Rev. Lett. 123, 067204 (2019).
• [18] Yann Perrin, Benjamin Canals, and Nicolas Rougemaille. Extensive degeneracy, Coulomb phase and magnetic monopoles in artificial square ice. Nature 540, 410–413 (2016).
• [19] H. Fukazawa, R. G. Melko, R. Higashinaka, Y. Maeno, and M. J. P. Gingras. Magnetic anisotropy of the spin-ice compound Dy₂Ti₂O₇. Phys. Rev. B 65, 054410 (2002).
• [20] Leon Balents. Spin liquids in frustrated magnets. Nature 464, 199–208 (2010).
• [21] G. Baskaran, Z. Zou, and P. W. Anderson. The resonating valence bond state and high-T_c superconductivity — A mean field theory. Solid State Commun. 63, 973–976 (1987).
• [22] P. W. Anderson, G. Baskaran, Z. Zou, and T. Hsu. Resonating-valence-bond theory of phase transitions and superconductivity in La₂CuO₄-based compounds. Phys. Rev. Lett. 58, 2790–2793 (1987).
• [23] G. Baskaran and P. W. Anderson. Gauge theory of high-temperature superconductors and strongly correlated Fermi systems. Phys. Rev. B 37, 580–583 (1988).
• [24] Lucile Savary and Leon Balents. Quantum spin liquids: a review. Rep. Prog. Phys. 80, 016502 (2017).
• [25] Alexei Yu. Kitaev. Fault-tolerant quantum computation by anyons. Ann. Phys. (N. Y.) 303, 2–30 (2003).
• [26] Alexei Kitaev and Chris Laumann. Topological quantum phases and quantum computation. In Jesper Jacobsen, Stephane Ouvry, Vincent Pasquier, Didina Serban, and Leticia Cugliandolo, editors, Exact Methods in Low-dimensional Statistical Physics and Quantum Computing, volume 89 of Lecture Notes of the Les Houches Summer School, chapter 4. Oxford University Press, Oxford, U.K. (2008), arXiv:0904.2771.
• [27] Yi Zhou, Kazushi Kanoda, and Tai-Kai Ng. Quantum spin liquid states. Rev. Mod. Phys. 89, 025003 (2017).
• [28] Nic Shannon, Tsutomu Momoi, and Philippe Sindzingre. Nematic Order in Square Lattice Frustrated Ferromagnets. Phys. Rev. Lett. 96, 027213 (2006).
• [29] John A. Mydosh and Peter M. Oppeneer. Colloquium: Hidden order, superconductivity, and magnetism: The unsolved case of URu₂Si₂. Rev. Mod. Phys. 83, 1301–1322 (2011).
• [30] Joseph A. M. Paddison, Henrik Jacobsen, Oleg A. Petrenko, Maria Teresa Fernández-Díaz, Pascale P. Deen, and Andrew L. Goodwin. Hidden order in spin-liquid Gd₃Ga₅O₁₂. Science 350, 179–181 (2015).
• [31] Yao-Dong Li, Xiaoqun Wang, and Gang Chen. Hidden multipolar orders of dipole-octupole doublets on a triangular lattice. Phys. Rev. B 94, 201114 (2016).
• [32] H. Takatsu, S. Onoda, S. Kittaka, A. Kasahara, Y. Kono, T. Sakakibara, Y. Kato, B. Fåk, J. Ollivier, J. W. Lynn, T. Taniguchi, M. Wakita, and H. Kadowaki. Quadrupole Order in the Frustrated Pyrochlore Tb_2+xTi_2-xO_7+y. Phys. Rev. Lett. 116, 217201 (2016).
• [33] Qiang Luo, Shijie Hu, Bin Xi, Jize Zhao, and Xiaoqun Wang. Ground-state phase diagram of an anisotropic spin-½ model on the triangular lattice. Phys. Rev. B 95, 165110 (2017).
• [34] A. F. Andreev and I. A. Grishchuk. Spin Nematics. J. Exp. Theor. Phys. 87, 467–475 (1984).
• [35] Andrey V. Chubukov. Fluctuations in spin nematics. J. Phys.: Condens. Matter 2, 1593 (1990).
• [36] Tsutomu Momoi, Philippe Sindzingre, and Nic Shannon. Octupolar Order in the Multiple Spin Exchange Model on a Triangular Lattice. Phys. Rev. Lett. 97, 257204 (2006).
• [37] Tsutomu Momoi, Philippe Sindzingre, and Kenn Kubo. Spin Nematic Order in Multiple-Spin Exchange Models on the Triangular Lattice. Phys. Rev. Lett. 108, 057206 (2012).
• [38] Roderich Moessner and John T. Chalker. Low-temperature properties of classical geometrically frustrated antiferromagnets. Phys. Rev. B 58, 12049–12062 (1998).
• [39] Andreas Läuchli, Frédéric Mila, and Karlo Penc. Quadrupolar Phases of the S=1 Bilinear-Biquadratic Heisenberg Model on the Triangular Lattice. Phys. Rev. Lett. 97, 087205 (2006).
• [40] Congjun Wu. Orbital Ordering and Frustration of p-Band Mott Insulators. Phys. Rev. Lett. 100, 200406 (2008).
• [41] J. Yamaura, K. Ohgushi, H. Ohsumi, T. Hasegawa, I. Yamauchi, K. Sugimoto, S. Takeshita, A. Tokuda, M. Takata, M. Udagawa, M. Takigawa, H. Harima, T. Arima, and Z. Hiroi. Tetrahedral Magnetic Order and the Metal-Insulator Transition in the Pyrochlore Lattice of Cd₂Os₂O₇. Phys. Rev. Lett. 108, 247205 (2012).
• [42] M. Mourigal, M. Enderle, B. Fåk, R. K. Kremer, J. M. Law, A. Schneidewind, A. Hiess, and A. Prokofiev. Evidence of a Bond-Nematic Phase in LiCuVO₄. Phys. Rev. Lett. 109, 027203 (2012).
• [43] O. Janson, J. Richter, P. Sindzingre, and H. Rosner. Coupled frustrated quantum spin-½ chains with orbital order in volborthite Cu₃V₂O₇(OH)₂ · 2H₂O. Phys. Rev. B 82, 104434 (2010).
• [44] R. Wawrzyńczak, Y. Tanaka, M. Yoshida, Y. Okamoto, P. Manuel, N. Casati, Z. Hiroi, M. Takigawa, and G. J. Nilsen. Classical Spin Nematic Transition in LiGa_0.95In_0.05Cr₄O₈. Phys. Rev. Lett. 119, 087201 (2017).
• [45] A. Orlova, E. L. Green, J. M. Law, D. I. Gorbunov, G. Chanda, S. Krämer, M. Horvatić, R. K. Kremer, J. Wosnitza, and G. L. J. A. Rikken. Nuclear Magnetic Resonance Signature of the Spin-Nematic Phase in LiCuVO₄ at High Magnetic Fields. Phys. Rev. Lett. 118, 247201 (2017).
• [46] Jacques Villain, R. Bidaux, J.-P. Carton, and R. Conte. Order as an effect of disorder. J. Phys. France 41, 1263–1272 (1980).
• [47] Doron Bergman, Jason Alicea, Emanuel Gull, Simon Trebst, and Leon Balents. Order-by-disorder and spiral spin-liquid in frustrated diamond-lattice antiferromagnets. Nat. Phys. 3, 487–491 (2007).
• [48] John T. Chalker, Peter C. W. Holdsworth, and E. F. Shender. Hidden order in a frustrated system: Properties of the Heisenberg Kagomé antiferromagnet. Phys. Rev. Lett. 68, 855–858 (1992).
• [49] I. Ritchey, P. Chandra, and P. Coleman. Spin folding in the two-dimensional Heisenberg kagomé antiferromagnet. Phys. Rev. B 47, 15342–15345 (1993).
• [50] Michael E. Zhitomirsky. Field-Induced Transitions in a Kagomé Antiferromagnet. Phys. Rev. Lett. 88, 057204 (2002).
• [51] Michael E. Zhitomirsky. Octupolar ordering of classical kagome antiferromagnets in two and three dimensions. Phys. Rev. B 78, 094423 (2008).
• [52] Mathieu Taillefumier, Owen Benton, Han Yan, L. D. C. Jaubert, and Nic Shannon. Competing Spin Liquids and Hidden Spin-Nematic Order in Spin Ice with Frustrated Transverse Exchange. Phys. Rev. X 7, 041057 (2017).
• [53] Owen Benton, Ludovic D. C. Jaubert, Rajiv R. P. Singh, Jaan Oitmaa, and Nic Shannon. Quantum Spin Ice with Frustrated Transverse Exchange: From a π-Flux Phase to a Nematic Quantum Spin Liquid. Phys. Rev. Lett. 121, 067201 (2018).