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.

ac Wien Effect in Spin Ice, Manifest in Nonlinear, Nonequilibrium Susceptibility.

The Wien effect is a model process for field-induced charge creation. Here it is derived for a nonelectrical system: the spin ice "magnetolyte"-a unique system showing perfect charge symmetry. An entropic reaction field, analogous to the Jaccard field in ice, opposes direct current, but a frequency window exists in which the Wien effect for magnetolyte and electrolyte are indistinguishable. The universal enhancement of monopole density speeds up the magnetization dynamics, which manifests in the nonlinear, nonequilibrium ac susceptibility. This is a rare instance where such effects may be calculated, providing new insights for electrolytes. Experimental predictions are made for Dy2Ti2O7 spin ice.

Onsager's Wien effect on a lattice.

The second Wien effect describes the nonlinear, non-equilibrium response of a weak electrolyte in moderate to high electric fields. Onsager's 1934 electrodiffusion theory, along with various extensions, has been invoked for systems and phenomena as diverse as solar cells, surfactant solutions, water splitting reactions, dielectric liquids, electrohydrodynamic flow, water and ice physics, electrical double layers, non-ohmic conduction in semiconductors and oxide glasses, biochemical nerve response and magnetic monopoles in spin ice. In view of this technological importance and the experimental ubiquity of such phenomena, it is surprising that Onsager's Wien effect has never been studied by numerical simulation. Here we present simulations of a lattice Coulomb gas, treating the widely applicable case of a double equilibrium for free charge generation. We obtain detailed characterization of the Wien effect and confirm the accuracy of the analytical theories as regards the field evolution of the free charge density and correlations. We also demonstrate that simulations can uncover further corrections, such as how the field-dependent conductivity may be influenced by details of microscopic dynamics. We conclude that lattice simulation offers a powerful means by which to model and investigate system-specific corrections to the Onsager theory, and thus constitutes a valuable tool for detailed theoretical studies of the numerous practical applications of the second Wien effect.

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.

Quantum order by disorder and accidental soft mode in Er2Ti2O7.

Motivated by recent neutron scattering experiments, we derive and study an effective "pseudodipolar" spin-1/2 model for the XY pyrochlore antiferromagnet Er(2)Ti(2)O(7). While a bond-dependent in-plane exchange anisotropy removes any continuous symmetry, it does lead to a one-parameter 'accidental' classical degeneracy. This degeneracy is lifted by quantum fluctuations in favor of the noncoplanar spin structure observed experimentally-a rare experimental instance of quantum order by disorder. A non-Goldstone low-energy mode is present in the excitation spectrum in accordance with inelastic neutron scattering data. Our theory also resolves the puzzle of the experimentally observed continuous ordering transition, absent from previous models.

Magnetic monopole dynamics in spin ice.

One of the most remarkable examples of emergent quasi-particles is that of the 'fractionalization' of magnetic dipoles in the low energy configurations of materials known as 'spin ice' into free and unconfined magnetic monopoles interacting via Coulomb's 1/r law (Castelnovo et al 2008 Nature 451 42-5). Recent experiments have shown that a Coulomb gas of magnetic charges really does exist at low temperature in these materials and this discovery provides a new perspective on otherwise largely inaccessible phenomenology. In this paper, after a review of the different spin ice models, we present detailed results describing the diffusive dynamics of monopole particles starting both from the dipolar spin ice model and directly from a Coulomb gas within the grand canonical ensemble. The diffusive quasi-particle dynamics of real spin ice materials within the 'quantum tunnelling' regime is modelled with Metropolis dynamics, with the particles constrained to move along an underlying network of oriented paths, which are classical analogues of the Dirac strings connecting pairs of Dirac monopoles.

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.

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

We examine the statistical mechanics of spin-ice materials with a [100] magnetic field. We show that the approach to saturated magnetization is, in the low-temperature limit, an example of a 3D Kasteleyn transition, which is topological in the sense that magnetization is changed only by excitations that span the entire system. We study the transition analytically and using a Monte Carlo cluster algorithm, and compare our results with recent data from experiments on Dy2Ti2O7.

Universal window for two-dimensional critical exponents.

Two-dimensional condensed matter is realized in increasingly diverse forms that are accessible to experiment and of potential technological value. The properties of these systems are influenced by many length scales and reflect both generic physics and chemical detail. To unify their physical description is therefore a complex and important challenge. Here we investigate the distribution of experimentally estimated critical exponents, ß, that characterize the evolution of the order parameter through the ordering transition. The distribution is found to be bimodal and bounded within a window ∼0.1≤ß≤0.25, facts that are only in partial agreement with the established theory of critical phenomena. In particular, the bounded nature of the distribution is impossible to reconcile with the existing theory for one of the major universality classes of two-dimensional behaviour-the XY model with four-fold crystal field-which predicts a spectrum of non-universal exponents bounded only from below. Through a combination of numerical and renormalization group arguments we resolve the contradiction between theory and experiment and demonstrate how the 'universal window' for critical exponents observed in experiment arises from a competition between marginal operators.

Statistics of extremal intensities for Gaussian interfaces.

The extremal Fourier intensities are studied for stationary Edwards-Wilkinson-type, Gaussian, interfaces with power-law dispersion. We calculate the probability distribution of the maximal intensity and find that, generically, it does not coincide with the distribution of the integrated power spectrum (i.e., roughness of the surface), nor does it obey any of the known extreme statistics limit distributions. The Fisher-Tippett-Gumbel limit distribution is, however, recovered in three cases: (i) in the nondispersive (white noise) limit, (ii) for high dimensions, and (iii) when only short-wavelength modes are kept. In the last two cases the limit distribution emerges in nonconventional scenarios.

Intermittency and non-Gaussian fluctuations of the global energy transfer in fully developed turbulence.

We address the experimentally observed non-Gaussian fluctuations for the energy injected into a closed turbulent flow at fixed Reynolds number. We propose that the power fluctuations mirror the internal kinetic energy fluctuations. Using a stochastic cascade model, we construct the excess kinetic energy as the sum over the energy transfers at different levels of the cascade. We find an asymmetric distribution that strongly resembles the experimental data. The asymmetry is an explicit consequence of intermittency and the global measure is dominated by small scale events correlated over the entire system. Our calculation is consistent with the statistical analogy recently made between a confined turbulent flow and a critical system of finite size.