We consider a system of spins on the sites of a three-dimensional pyrochlore lattice of corner-sharing tetrahedra interacting with a predominant effective $xy$ exchange. In particular, we investigate the selection of a long-range ordered state with b
roken discrete symmetry induced by thermal fluctuations near the critical region. At the standard mean-field theory (s-MFT) level, in a region of the parameter space of this Hamiltonian that we refer to as $Gamma_5$ region, the ordered state possesses an accidental $U(1)$ degeneracy. In this paper, we show that fluctuations beyond s-MFT lift this degeneracy by selecting one of two states (so-called $psi_2$ and $psi_3$) from the degenerate manifold, thus exposing a certain form of order-by-disorder (ObD). We analytically explore this selection at the microscopic level and close to criticality by elaborating upon and using an extension of the so-called TAP method, originally developed by Thouless, Anderson and Palmer to study the effect of fluctuations in spin glasses. We also use a single-tetrahedron cluster-mean-field theory (c-MFT) to explore over what minimal length scale fluctuations can lift the degeneracy. We find the phase diagrams obtained by these two methods to be somewhat different since c-MFT only includes the shortest-range fluctuations. General symmetry arguments used to construct a Ginzburg-Landau theory to lowest order in the order parameters predict that a weak magnetic moment, $m_z$, along the local $langle 111 rangle$ (${hat z}$) direction is generically induced for a system ordering into a $psi_2$ state, but not so for $psi_3$ ordering. Both E-TAP and c-MFT calculations confirm this weak fluctuation-induced $m_z$ moment. Using a Ginzburg-Landau theory, we discuss the phenomenology of multiple phase transitions below the paramagnetic phase transition and within the $Gamma_5$ long-range ordered phase.
Motivated by recent experimental and theoretical progress on the Er2Ti2O7 pyrochlore XY antiferromagnet, we study the problem of quantum order-by-disorder in pyrochlore XY systems. We consider the most general nearest-neighbor pseudo spin-1/2 Hamilto
nian for such a system characterized by anisotropic spin-spin couplings J_e = [J_pm, J_{pmpm}, J_{zpm}, J_{zz}] and construct zero-temperature phase diagrams. Combining symmetry arguments and spin-wave calculations, we show that the ground state phase boundaries between the two candidate ground states of the Gamma_5 irreducible representation, the psi_2 and psi_3 (basis) states, are rather accurately determined by a cubic equation in J_{pm}J_{pmpm})/J_{zpm}^2. Depending on the value of J_{zz}, there can be one or three phase boundaries that separate alternating regions of psi_2 and psi_3 states. In particular, we find for sufficiently small J_{zz}/J_{pm} a narrow psi_2 sliver sandwiched between two psi_3 regions in the J_{pmpm}/J_pm vs J_{zpm}/J_pm phase diagram. Our results further illustrate the very large potential sensitivity of the ground state of XY pyrochlore systems to minute changes in their spin Hamiltonian. Using the experimentally determined J_3 and g-tensor values for Er2Ti2O7, we show that the heretofore neglected long-range 1/r^3 magnetostatic dipole-dipole interactions do not change the conclusion that Er2Ti2O7 has a psi_2 ground state induced via a quantum order-by-disorder mechanism. We propose that the CdDy2Se4 chalcogenide spinel, in which the Dy^{3+} ions form a pyrochlore lattice and may be XY-like, could prove interesting to investigate.
The spin ice materials, including Ho2Ti2O7 and Dy2Ti2O7, are rare earth pyrochlore magnets which, at low temperatures, enter a constrained paramagnetic state with an emergent gauge freedom. Remarkably, the spin ices provide one of very few experiment
ally realised examples of fractionalization because their elementary excitations can be regarded as magnetic monopoles and, over some temperature range, the spin ice materials are best described as liquids of these emergent charges. In the presence of quantum fluctuations, one can obtain, in principle, a quantum spin liquid descended from the classical spin ice state characterised by emergent photon-like excitations. Whereas in classical spin ices the excitations are akin to electrostatic charges, in the quantum spin liquid these charges interact through a dynamic and emergent electromagnetic field. In this review, we describe the latest developments in the study of such a quantum spin ice, focussing on the spin liquid phenomenology and the kinds of materials where such a phase might be found.
We study the problem of partially ordered phases with periodically arranged disordered (paramagnetic) sites on the pyrochlore lattice, a network of corner-sharing tetrahedra. The periodicity of these phases is characterized by one or more wave vector
s k=(1/2 1/2 1/2). Starting from a general microscopic Hamiltonian including anisotropic nearest-neighbor exchange, long-range dipolar interactions and second- and third-nearest neighbor exchange, we identify using standard mean-field theory (s-MFT) an extended range of interaction parameters that support partially ordered phases. We demonstrate that thermal fluctuations ignored in s-MFT are responsible for the selection of one particular partially ordered phase, e.g. the 4-k phase over the 1-k phase. We suggest that the transition into the 4-k phase is continuous with its critical properties controlled by the cubic fixed point of a Ginzburg-Landau theory with a 4-component vector order-parameter. By combining an extension of the Thouless-Anderson-Palmer method originally used to study fluctuations in spin glasses with parallel-tempering Monte-Carlo simulations, we establish the phase diagram for different types of partially ordered phases. Our results elucidate the long-standing puzzle concerning the origin of the 4-k partially ordered phase observed in the Gd2Ti2O7 dipolar pyrochlore antiferromagnet below its paramagnetic phase transition temperature.
Several rare earth magnetic pyrochlore materials are well modeled by a spin-1/2 quantum Hamiltonian with anisotropic exchange parameters Js. For the Er2Ti2O7 material, the Js were recently determined from high-field inelastic neutron scattering measu
rements. Here, we perform high-temperature (T) series expansions to compute the thermodynamic properties of this material using these Js. Comparison with experimental data show that the model describes the material very well including the finite temperature phase transition to an ordered phase at Tc~1.2 K. We show that high temperature expansions give identical results for different q=0 xy order parameter susceptibilities up to 8th order in beta=1/T (presumably to all orders in beta). Conversely, a non-linear susceptibility related to the 6th power of the order parameter reveals a thermal order-by-disorder selection of the same non-colinear psi_2 state as found in Er2Ti2O7.
Spin ice materials, such as Dy2Ti2O7 and Ho2Ti2O7, have been the subject of much interest for over the past fifteen years. Their low temperature strongly correlated state can be mapped onto the proton disordered state of common water ice and, consequ
ently, spin ices display the same low temperature residual Pauling entropy as water ice. Interestingly, it was found in a previous study [X. Ke {it et. al.} Phys. Rev. Lett. {bf 99}, 137203 (2007)] that, upon dilution of the magnetic rare-earth ions (Dy^{3+} and Ho^{3+}) by non-magnetic Yttrium (Y^{3+}) ions, the residual entropy depends {it non-monotonically} on the concentration of Y^{3+} ions. In the present work, we report results from Monte Carlo simulations of site-diluted microscopic dipolar spin ice models (DSIM) that account quantitatively for the experimental specific heat measurements, and thus also for the residual entropy, as a function of dilution, for both Dy2Ti2O7 and Ho2Ti2O7. The main features of the dilution physics displayed by the magnetic specific heat data are quantitatively captured by the diluted DSIM up to, and including, 85% of the magnetic ions diluted (x=1.7). The previously reported departures in the residual entropy between Dy2Ti2O7 versus Ho2Ti2O7, as well as with a site-dilution variant of Paulings approximation, are thus rationalized through the site-diluted DSIM. For 90% (x=1.8) and 95% (x=1.9) of the magnetic ions diluted, we find a significant discrepancy between the experimental and Monte Carlo specific heat results. We discuss some possible reasons for this disagreement.
The thermodynamic properties of the pyrochlore Yb2Ti2O7 material are calculated using the numericallinked-cluster (NLC) calculation method for an effective anisotropic-exchange spin-1/2 Hamiltonian with parameters recently determined by fitting the n
eutron scattering spin wave data obtained at high magnetic field h. Magnetization, M(T,h), as a function of temperature T and for different magnetic fields h applied along the three high symmetry directions [100], [110] and [111], are compared with experimental measurements on the material for temperature T>1.8K. The excellent agreement between experimentally measured and calculated M(T,h) over the entire temperature and magnetic field range considered provides strong quantitative validation of the effective Hamiltonian. It also confirms that fitting the high-field neutron spin wave spectra in the polarized paramagnetic state is an excellent method for determining the microscopic exchange constants of rare-earth insulating magnets that are described by an effective spin-1/2 Hamiltonian. Finally, we present results which demonstrate that a recent analysis of the polarized neutron scattering intensity of Yb2Ti2O7 using a random phase approximation (RPA) method [Chang et al., Nature Communications {3}, 992 (2012)] does not provide a good description of M(T,h) for $Tlesssim 10$ K, that is in the entire temperature regime where correlations become non-negligible.
We use numerical linked cluster (NLC) expansions to compute the specific heat, C(T), and entropy, S(T), of a quantum spin ice model of Yb2Ti2O7 using anisotropic exchange interactions recently determined from inelastic neutron scattering measurements
and find good agreement with experimental calorimetric data. In the perturbative weak quantum regime, this model has a ferrimagnetic ordered ground state, with two peaks in C(T): a Schottky anomaly signalling the paramagnetic to spin ice crossover followed at lower temperature by a sharp peak accompanying a first order phase transition to the ferrimagnetic state. We suggest that the two C(T) features observed in Yb2Ti2O7 are associated with the same physics. Spin excitations in this regime consist of weakly confined spinon-antispinon pairs. We suggest that conventional ground state with exotic quantum dynamics will prove a prevalent characteristic of many real quantum spin ice materials.
We present new high resolution inelastic neutron scattering data on the candidate spin liquid Tb2Ti2O7. We find that there is no evidence for a zero field splitting of the ground state doublet within the 0.2 K resolution of the instrument. This resul
t contrasts with a pair of recent works on Tb2Ti2O7 claiming that the spin liquid behavior can be attributed to a 2 K split singlet-singlet single-ion spectrum at low energies. We also reconsider the entropy argument presented in Chapuis {it et al.} as further evidence of a singlet-singlet crystal field spectrum. We arrive at the conclusion that estimates of the low temperature residual entropy drawn from heat capacity measurements are a poor guide to the single ion spectrum without understanding the nature of the correlations.
The rare earth pyrochlore magnet Yb2Ti2O7 is among a handful of materials that apparently exhibit no long range order down to the lowest explored temperatures and well below the Curie-Weiss temperature. Paramagnetic neutron scattering on a single cry
stal sample has revealed the presence of anisotropic correlations and recent work has led to the proposal of a detailed microscopic Hamiltonian for this material involving significantly anisotropic exchange. In this article, we compute the local sublattice susceptibility of Yb2Ti2O7 from the proposed model and compare with the measurements of Cao and coworkers [Physical Review Letters, {103}, 056402 (2009)], finding quite good agreement. In contrast, a model with only isotropic exchange and long range magnetostatic dipoles gives rise to a local susceptiblity that is inconsistent with the data.
