Abundance of Hard-Hexagon Crystals in the Quantum Pyrochlore Antiferromagnet.
Physical review letters 131:9 (2023) 096702
Abstract:
We propose a simple family of valence-bond crystals as potential ground states of the S=1/2 and S=1 Heisenberg antiferromagnet on the pyrochlore lattice. Exponentially numerous in the linear size of the system, these can be visualized as hard-hexagon coverings, with each hexagon representing a resonating valence-bond ring. This ensemble spontaneously breaks rotation, inversion, and translation symmetries. A simple, yet accurate, variational wave function allows a precise determination of the energy, confirmed by the density matrix renormalization group and numerical linked cluster expansion, and extended by an analysis of excited states. The identification of the origin of the stability indicates applicability to a broad class of frustrated lattices, which we demonstrate for the checkerboard and ruby lattices. Our work suggests a perspective on such quantum magnets, in which unfrustrated motifs are effectively uncoupled by the frustration of their interactions.Quantum spin ice response to a magnetic field in the dipole-octupole pyrochlore Ce2Zr2 O7
Physical Review B 108:5 (2023)
Abstract:
The pyrochlore magnet Ce2Zr2O7 has attracted much attention as a quantum spin ice candidate whose novelty derives in part from the dipolar-octupolar nature of the Ce3+ pseudospin-1/2 degrees of freedom it possesses. We report heat capacity measurements on single crystal samples of Ce2Zr2O7 down to T∼0.1K in a magnetic field along the [1,1¯,0] direction. These measurements show that the broad hump in the zero-field heat capacity moves higher in temperature with increasing field strength and is split into two separate humps by the [1,1¯,0] magnetic field at ∼2T. These separate features are due to the decomposition of the pyrochlore lattice into effectively decoupled chains for fields in this direction: One set of chains (α chains) is polarized by the field while the other (β chains) remains free. This situation is similar to that observed in the classical spin ices Ho2Ti2O7 and Dy2Ti2O7, but with the twist that here the strong transverse exchange interactions produce substantial quantum effects. Our theoretical modeling suggests that the β chains are close to a critical state, with nearly-gapless excitations. We also report elastic and inelastic neutron scattering measurements on single crystal Ce2Zr2O7 in [1,1¯,0] and [0,0,1] magnetic fields at temperatures down to T=0.03K. The elastic scattering behaves consistently with the formation of independent chains for a [1,1¯,0] field, while the [0,0,1] field produces a single field-induced elastic magnetic Bragg peak at (0,2,0) and equivalent wavevectors, indicating a polarized spin ice state for fields above ∼3T. For both [1,1¯,0] and [0,0,1] magnetic fields, our inelastic neutron scattering results show an approximately-dispersionless continuum of scattering that increases in both energy and intensity with increasing field strength. By modeling the complete set of experimental data using numerical linked cluster and semiclassical molecular dynamics calculations, we demonstrate the dominantly multipolar nature of the exchange interactions in Ce2Zr2O7 and the smallness of the parameter θ, which controls the mixing between dipolar and octupolar degrees of freedom. These results 91̽»¨ previous estimates of the microscopic exchange parameters and place strong constraints on the theoretical description of this prominent spin ice candidate.Arresting dynamics in hardcore spin models
Physical Review B American Physical Society (APS) 107:18 (2023) l180302
Random-bond Ising model and its dual in hyperbolic spaces.
Physical review. E 107:2-1 (2023) 024125
Abstract:
We analyze the thermodynamic properties of the random-bond Ising model (RBIM) on closed hyperbolic surfaces using Monte Carlo and high-temperature series expansion techniques. We also analyze the dual-RBIM, that is, the model that in the absence of disorder is related to the RBIM via the Kramers-Wannier duality. Even on self-dual lattices this model is different from the RBIM, unlike in the Euclidean case. We explain this anomaly by a careful rederivation of the Kramers-Wannier duality. For the (dual-)RBIM, we compute the paramagnet-to-ferromagnet phase transition as a function of both temperature T and the fraction of antiferromagnetic bonds p. We find that as temperature is decreased in the RBIM, the paramagnet gives way to either a ferromagnet or a spin-glass phase via a second-order transition compatible with mean-field behavior. In contrast, the dual-RBIM undergoes a strongly first-order transition from the paramagnet to the ferromagnet both in the absence of disorder and along the Nishimori line. We study both transitions for a variety of hyperbolic tessellations and comment on the role of coordination number and curvature. The extent of the ferromagnetic phase in the dual-RBIM corresponds to the correctable phase of hyperbolic surface codes under independent bit- and phase-flip noise.Critical properties of the Ising model in hyperbolic space.
Physical review. E 101:2-1 (2020) 022124