Three-dimensional Kasteleyn transition: spin ice in a [100] field

Fragmented Spin Ice and Multi-k Ordering in Rare-Earth Antiperovskites

We study near-neighbor and dipolar Ising models on a lattice of corner-sharing octahedra. In an extended parameter range of both models, frustration between antiferromagnetism and a spin-ice-like three-in-three-out rule stabilizes a Coulomb phase with correlated dipolar and quadrupolar spin textures, both yielding distinctive neutron-scattering signatures. Strong further-neighbor perturbations cause the two components to order independently, resulting in unusual multi-k orders. We propose experimental realizations of our model in rare-earth antiperovskites.

Proton strings and rings in atypical nucleation of ferroelectricity in ice

Ordinary ice has a proton-disordered phase which is kinetically metastable, unable to reach, spontaneously, the ferroelectric (FE) ground state at low temperature where a residual Pauling entropy persists. Upon light doping with KOH at low temperature, the transition to FE ice takes place, but its microscopic mechanism still needs clarification. We introduce a lattice model based on dipolar interactions plus a competing, frustrating term that enforces the ice rule (IR). In the absence of IR-breaking defects, standard Monte Carlo (MC) simulation leaves this ice model stuck in a state of disordered proton ring configurations with the correct Pauling entropy. A replica exchange accelerated MC sampling strategy succeeds, without open path moves, interfaces, or off-lattice configurations, in equilibrating this defect-free ice, reaching its low-temperature FE order through a well-defined first-order phase transition. When proton vacancies mimicking the KOH impurities are planted into the IR-conserving lattice, they enable standard MC simulation to work, revealing the kinetics of evolution of ice from proton disorder to partial FE order below the transition temperature. Replacing ordinary nucleation, each impurity opens up a proton ring generating a linear string, an actual FE hydrogen bond wire that expands with time. Reminiscent of those described for spin ice, these impurity-induced strings are proposed to exist in doped water ice too, where IRs are even stronger. The emerging mechanism yields a dependence of the long-time FE order fraction upon dopant concentration, and upon quenching temperature, that compares favorably with that known in real-life KOH doped ice.

Anomalous out-of-equilibrium dynamics in the spin-ice material Dy2Ti2O7under moderate magnetic fields

We study experimentally and numerically the dynamics of the spin ice material Dy2Ti2O7in the low temperature (T) and moderate magnetic field (B) regime (T∈ [0.1, 1.7]  K,B∈ [0, 0.3]  T). Our objective is to understand the main physics shaping the out-of-equilibrium magnetisation vs temperature curves in two different regimes. Very far from equilibrium, turning on the magnetic field after having cooled the system in zero field (ZFC) can increase the concentration of magnetic monopoles (localised thermal excitations present in these systems); this accelerates the dynamics. Similarly to electrolytes, this occurs through dissociation of bound monopole pairs. However, for spin ices the polarisation of the vacuum out of which the monopole pairs are created is a key factor shaping the magnetisation curves, with no analog. We observe a threshold field near 0.2 T for this fast dynamics to take place, linked to the maximum magnetic force between the attracting pairs. Surprisingly, within a regime of low temperatures and moderate fields, an extended Ohm's law can be used to describe the ZFC magnetisation curve obtained with the dipolar spin-ice model. However, in real samples the acceleration of the dynamics appears even sharper than in simulations, possibly due to the presence of avalanches. On the other hand, the effect of the field nearer equilibrium can be just the opposite to that at very low temperatures. Single crystals, as noted before for powders, abandon equilibrium at a blocking temperatureTBwhich increases with field. Curiously, this behaviour is present in numerical simulations even within the nearest-neighbours interactions model. Simulations and experiments show that the increasing trend inTBis stronger forB‖[100]. This suggests that the field plays a part in the dynamical arrest through monopole suppression, which is quite manifest for this field orientation.

The history of spin ice

This review is a study of how the idea of spin ice has evolved over the years, with a focus on the scientific questions that have come to define the subject. Since our initial discovery of spin ice in 1997, there have been well over five thousand papers that discuss it, and in the face of such detail, it must be difficult for the curious observer to 'see the wood for the trees'. To help in this task, we go in search of the biggest insight to have emerged from the study of spin ice. On the way, we identify highlights and outstanding puzzles, and celebrate the inspirational role that Roger Cowley played in the early years.

Husimi-cactus approximation study on the diluted spin ice

We investigate dilution effects on the classical spin-ice materials such as Ho_{2}Ti_{2}O_{7} and Dy_{2}Ti_{2}O_{7}. In particular, we derive a formula of the thermodynamic quantities as functions of the temperature and a nonmagnetic ion concentration based on a Husimi-cactus approximation. We find that the formula predicts a dilution-induced crossover from the cooperative to the conventional paramagnets in a ground state, and that it also reproduces the "generalized Pauling's entropy" given by Ke et al. To verify the formula from a numerical viewpoint, we compare these results with Monte Carlo simulation calculation data, and then find good agreement for all parameter values.

Adiabatic cooling processes in frustrated magnetic systems with pyrochlore structure

We investigate in detail the process of adiabatic cooling in the framework of the exactly solvable antiferromagnetic spin-1/2 Ising model in the presence of the external magnetic field on an approximate lattice with pyrochlore structure. The behavior of the entropy of the model is studied and exact values of the residual entropies of all ground states are found. The temperature variation of the system under adiabatic (de)magnetization is investigated and the central role of the macroscopically degenerated ground states in cooling processes is explicitly demonstrated. It is shown that the model parameter space of the studied geometrically frustrated system is divided into five disjunct regions with qualitatively different processes of the adiabatic cooling. The effectiveness of the adiabatic (de)magnetization cooling in the studied model is compared to the corresponding processes in paramagnetic salts. It is shown that the processes of the adiabatic cooling in the antiferromagnetic frustrated systems are much more effective especially in nonzero external magnetic fields. It means that the frustrated magnetic materials with pyrochlore structure can be considered as very promising refrigerants mainly in the situations with nonzero final values of the magnetic field.

Fragmentation in spin ice from magnetic charge injection

The complexity embedded in condensed matter fertilizes the discovery of new states of matter, enriched by ingredients like frustration. Illustrating examples in magnetic systems are Kitaev spin liquids, skyrmions phases, or spin ices. These unconventional ground states support exotic excitations, for example the magnetic charges in spin ices, also called monopoles. Here, we propose a mechanism to inject monopoles in a spin ice at equilibrium through a staggered magnetic field. We show theoretically, and demonstrate experimentally in the Ho₂Ir₂O₇ pyrochlore iridate, that it results in the stabilization of a monopole crystal, which exhibits magnetic fragmentation. In this new state of matter, the magnetic moment fragments into an ordered part and a persistently fluctuating one. Compared to conventional spin ices, the different nature of the excitations in this fragmented state opens the way to tunable field-induced and dynamical behaviors.

Exploring unconventional magnetism facilities both fundamental understanding of materials and their real applications. Here the authors demonstrate that a magnetic monopole crystal is stabilized by a staggered magnetic field in the pyrochlore iridate Ho₂Ir₂O₇, leading to a fragmented magnetization.

Spin ice Thin Film: Surface Ordering, Emergent Square ice, and Strain Effects

Motivated by recent realizations of Dy_{2}Ti_{2}O_{7} and Ho_{2}Ti_{2}O_{7} spin ice thin films, and more generally by the physics of confined gauge fields, we study a model spin ice thin film with surfaces perpendicular to the [001] cubic axis. The resulting open boundaries make half of the bonds on the interfaces inequivalent. By tuning the strength of these inequivalent "orphan" bonds, dipolar interactions induce a surface ordering equivalent to a two-dimensional crystallization of magnetic surface charges. This surface ordering may also be expected on the surfaces of bulk crystals. For ultrathin films made of one cubic unit cell, once the surfaces have ordered, a square ice phase is stabilized over a finite temperature window. The square ice degeneracy is lifted at lower temperature and the system orders in analogy with the well-known F transition of the 6-vertex model. To conclude, we consider the addition of strain effects, a possible consequence of interface mismatches at the film-substrate interface. Our simulations qualitatively confirm that strain can lead to a smooth loss of Pauling entropy upon cooling, as observed in recent experiments on Dy_{2}Ti_{2}O_{7} films.

The 3D Kasteleyn transition in dipolar spin ice: a numerical study with the conserved monopoles algorithm

We study the three-dimensional Kasteleyn transition in both nearest neighbours and dipolar spin ice models using an algorithm that conserves the number of excitations. We first limit the interactions range to nearest neighbours to test the method in the presence of a field applied along [Formula: see text], and then focus on the dipolar spin ice model. The effect of dipolar interactions, which is known to be greatly self screened at zero field, is particularly strong near full polarization. It shifts the Kasteleyn transition to lower temperatures, which decreases  ≈0.4 K for the parameters corresponding to the best known spin ice materials, [Formula: see text] and [Formula: see text]. This shift implies effective dipolar fields as big as 0.05 T opposing the applied field, and thus favouring the creation of 'strings' of reversed spins. We compare the reduction in the transition temperature with results in previous experiments, and study the phenomenon quantitatively using a simple molecular field approach. Finally, we relate the presence of the effective residual field to the appearance of string-ordered phases at low fields and temperatures, and we check numerically that for fields applied along [Formula: see text] there are only three different stable phases at zero temperature.

Maxwell electromagnetism as an emergent phenomenon in condensed matter

The formulation of a complete theory of classical electromagnetism by Maxwell is one of the milestones of science. The capacity of many-body systems to provide emergent mini-universes with vacua quite distinct from the one we inhabit was only recognized much later. Here, we provide an account of how simple systems of localized spins manage to emulate Maxwell electromagnetism in their low-energy behaviour. They are much less constrained by symmetry considerations than the relativistically invariant electromagnetic vacuum, as their substrate provides a non-relativistic background with even translational invariance broken. They can exhibit rich behaviour not encountered in conventional electromagnetism. This includes the existence of magnetic monopole excitations arising from fractionalization of magnetic dipoles; as well as the capacity of disorder, by generating defects on the lattice scale, to produce novel physics, as exemplified by topological spin glassiness or random Coulomb magnetism.This article is part of the themed issue 'Unifying physics and technology in light of Maxwell's equations'.

Intermediate magnetization state and competing orders in Dy2Ti2O7 and Ho2Ti2O7

Among the frustrated magnetic materials, spin-ice stands out as a particularly interesting system. Residual entropy, freezing and glassiness, Kasteleyn transitions and fractionalization of excitations in three dimensions all stem from a simple classical Hamiltonian. But is the usual spin-ice Hamiltonian a correct description of the experimental systems? Here we address this issue by measuring magnetic susceptibility in the two most studied spin-ice compounds, Dy₂Ti₂O₇ and Ho₂Ti₂O₇, using a vector magnet. Using these results, and guided by a theoretical analysis of possible distortions to the pyrochlore lattice, we construct an effective Hamiltonian and explore it using Monte Carlo simulations. We show how this Hamiltonian reproduces the experimental results, including the formation of a phase of intermediate polarization, and gives important information about the possible ground state of real spin-ice systems. Our work suggests an unusual situation in which distortions might contribute to the preservation rather than relief of the effects of frustration.

A classical Hamiltonian captures key properties of spin ice materials such as residual entropy and fractionalized excitations. Here, the authors present experimental results of the polarization transition that motivate a Hamiltonian with lattice distortions, which predicts an intermediate magnetization state and competing ground state orders.

Topological Aspects of Symmetry Breaking in Triangular-Lattice Ising Antiferromagnets

Using a specially designed Monte Carlo algorithm with directed loops, we investigate the triangular lattice Ising antiferromagnet with coupling beyond the nearest neighbors. We show that the first-order transition from the stripe state to the paramagnet can be split, giving rise to an intermediate nematic phase in which algebraic correlations coexist with a broken symmetry. Furthermore, we demonstrate the emergence of several properties of a more topological nature such as fractional edge excitations in the stripe state, the proliferation of double domain walls in the nematic phase, and the Kasteleyn transition between them. Experimental implications are briefly discussed.

Emergence of cooperative dynamics in fully packed classical dimers

The classical dimer model on the square lattice is a paradigmatic example of a system subject to strong local constraints. We study its behavior under local stochastic dynamics, by means of Monte Carlo simulations and theoretical arguments. We observe clear signatures of correlated dynamics in both global and local observables and over a broad range of time scales, indicating a breakdown of the simple continuum description that approximates well the statics. We show that this collective dynamics can be understood in terms of one-dimensional "strings" of high mobility, which govern both local and long-wavelength dynamical properties. We introduce a coarse-grained description of the strings, based on the Edwards-Wilkinson model, which leads to exact results in the limit of low string density and provides a detailed qualitative understanding of the dynamics in all flux sectors. We discuss the implications of our results for the dynamics of constrained systems more generally.

Crystal shape-dependent magnetic susceptibility and Curie law crossover in the spin ices Dy2Ti2O7 and Ho2Ti2O7

We present an experimental determination of the isothermal magnetic susceptibility of the spin ice materials Dy2Ti2O7 and Ho2Ti2O7 in the temperature range 1.8-300 K. The use of spherical crystals has allowed accurate correction for demagnetizing fields and allowed the true bulk isothermal susceptibility χT(T) to be estimated. This has been compared against a theoretical expression based on a Husimi tree approximation to the spin ice model. Agreement between experiment and theory is excellent at T > 10 K, but systematic deviations occur below that temperature. Our results largely resolve an apparent disagreement between neutron scattering and bulk measurements that has been previously noted. They also show that the use of non-spherical crystals in magnetization studies of spin ice may introduce very significant systematic errors, although we note some interesting--and possibly new--systematics concerning the demagnetizing factor in cuboidal samples. Finally, our results show how experimental susceptibility measurements on spin ices may be used to extract the characteristic energy scale of the system and the corresponding chemical potential for emergent magnetic monopoles.

Thermal phase transitions in artificial spin ice

We use the sixteen-vertex model to describe bidimensional artificial spin ice. We find excellent agreement between vertex densities in 15 differently grown samples and the predictions of the model. Our results demonstrate that the samples are in usual thermal equilibrium away from a critical point separating a disordered and an antiferromagnetic phase in the model. The second-order phase transition that we predict suggests that the spatial arrangement of vertices in near-critical artificial spin ice should be studied in more detail in order to verify whether they show the expected space and time long-range correlations.

Universal monopole scaling near transitions from the Coulomb phase

Certain frustrated systems, including spin ice and dimer models, exhibit a Coulomb phase at low temperatures, with power-law correlations and fractionalized monopole excitations. Transitions out of this phase, at which the effective gauge theory becomes confining, provide examples of unconventional criticality. This Letter studies the behavior at nonzero monopole density near such transitions, using scaling theory to arrive at universal expressions for the crossover phenomena. For a particular transition in spin ice, quantitative predictions are made by mapping to the XY model and confirmed using Monte Carlo simulations.

Quantum strings in quantum spin ice

We study quantum spin ice in an external magnetic field applied along a <100> direction. When quantum spin fluctuations are weak, elementary excitations are quantum strings with monopoles at their ends manifested as multiple spin wave branches in the dynamical structure factor. Strong quantum fluctuations make the string tension negative and give rise to the deconfinement of monopoles. We discuss our results in the light of recent neutron scattering experiments in Yb2Ti2O7.

Generalized Monte Carlo loop algorithm for two-dimensional frustrated Ising models

We introduce a generalized loop move (GLM) update for Monte Carlo simulations of frustrated Ising models on two-dimensional lattices with bond-sharing plaquettes. The GLM updates are designed to enhance Monte Carlo sampling efficiency when the system's low-energy states consist of an extensive number of degenerate or near-degenerate spin configurations, separated by large energy barriers to single spin flips. Through implementation on several frustrated Ising models, we demonstrate the effectiveness of the GLM updates in cases where both degenerate and near-degenerate sets of configurations are favored at low temperatures. The GLM update's potential to be straightforwardly extended to different lattices and spin interactions allows it to be readily adopted on many other frustrated Ising models of physical relevance.

Spin ice under pressure: symmetry enhancement and infinite order multicriticality

We study the low-temperature behavior of spin ice when uniaxial pressure induces a tetragonal distortion. There is a phase transition between a Coulomb liquid and a fully magnetized phase. Unusually, it combines features of discontinuous and continuous transitions: the order parameter exhibits a jump, but this is accompanied by a divergent susceptibility and vanishing domain wall tension. All these aspects can be understood as a consequence of an emergent SU(2) symmetry at the critical point. We map out a possible experimental realization.

Field-induced spin-ice-like orders in spin liquid Tb2Ti2O7

We have studied the field-induced magnetic structures in Tb2Ti2O7, in a wide temperature (0.3 < T < 270 K) and field (0 < H < 7 T) range, by single crystal polarized and unpolarized neutron diffraction, with H parallel[110] axis. A ferromagneticlike structure with k = 0 propagation vector is induced, whose local order at low field and low temperature is akin to spin ice. The four Tb ions separate in alpha and beta chains having different values of the magnetic moments, which is quantitatively explained by taking the crystal field anisotropy into account. Above 2 T and below 2 K, an antiferromagneticlike structure with k = (0,0,1) is induced besides the k = 0 structure. It shows a reentrant behavior and extends over a finite length scale. It occurs together with a broadening of the nuclear peaks, which suggests a field-induced distortion and magnetostriction effect.