Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Magnetic Excitations of the Classical Spin Liquid MgCr2O4

X. Bai, J. A. M. Paddison, E. Kapit, S. M. Koohpayeh, J.-J. Wen, S. E. Dutton, A. T. Savici, A. I. Kolesnikov, G. E. Granroth, C. L. Broholm, J. T. Chalker, and M. Mourigal
Phys. Rev. Lett. 122, 097201 – Published 5 March 2019
PDFHTMLExport Citation

Abstract

We report a comprehensive inelastic neutron-scattering study of the frustrated pyrochlore antiferromagnet MgCr2O4 in its cooperative paramagnetic regime. Theoretical modeling yields a microscopic Heisenberg model with exchange interactions up to third-nearest neighbors, which quantitatively explains all of the details of the dynamic magnetic response. Our work demonstrates that the magnetic excitations in paramagnetic MgCr2O4 are faithfully represented in the entire Brillouin zone by a theory of magnons propagating in a highly correlated paramagnetic background. Our results also suggest that MgCr2O4 is proximate to a spiral spin-liquid phase distinct from the Coulomb phase, which has implications for the magnetostructural phase transition in MgCr2O4.

  • Figure
  • Figure
  • Figure
  • Received 28 October 2018

DOI:https://doi.org/10.1103/PhysRevLett.122.097201

© 2019 American Physical Society

Physics Subject Headings (PhySH)

Condensed Matter, Materials & Applied Physics

Authors & Affiliations

X. Bai1, J. A. M. Paddison1,2, E. Kapit3,4, S. M. Koohpayeh5, J.-J. Wen5,*, S. E. Dutton6,†, A. T. Savici7, A. I. Kolesnikov7, G. E. Granroth7, C. L. Broholm5,7, J. T. Chalker3, and M. Mourigal1,5,‡

  • 1School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
  • 2Churchill College, University of Cambridge, Storey’s Way, Cambridge CB3 0DS, United Kingdom
  • 3Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Parks Road, Oxford OX1 3NP, United Kingdom
  • 4Department of Physics, Colorado School of Mines, Golden, Colorado 80401, USA
  • 5Institute for Quantum Matter and Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, Maryland 21218, USA
  • 6Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
  • 7Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA

  • *Present address: Stanford Institute for Materials and Energy Sciences, SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA.
  • Present address: Cavendish Laboratory, Department of Physics, University of Cambridge, JJ Thomson Ave., Cambridge CB3 0HE, UK.
  • mourigal@gatech.edu

Article Text (Subscription Required)

Click to Expand

Supplemental Material (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 122, Iss. 9 — 8 March 2019

Reuse & Permissions
Access Options
CHORUS

Article Available via CHORUS

Download Accepted Manuscript
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×

Images

  • Figure 1
    Figure 1

    (a) The pyrochlore lattice of Cr3+ ions (red spheres) in MgCr2O4 and definition of exchange interactions up to third neighbors. Note that J3a and J3b span the same distance but are not equivalent by symmetry. (b) Contour plot of the goodness of fit χ2 between calculations and neutron (blue solid lines) and bulk susceptibility (green dashed lines) measurements. FN exchange interactions J2 and J3a are fixed on a grid with J1 and J3b fitted at each grid point. The choice of J2±J3a as plotting axes highlights the nearly equivalent spin structure factors obtained for J2=J3a. Spin correlations are calculated using the self-consistent Gaussian approximation (SCGA) at T=20K. The red star is the best overall fit. (c) Momentum dependence of I0(Q)=F(|Q|)S(Q) and I1(Q)=F(|Q|)K(Q) along several paths of the Brillouin zone (BZ) at T=20K, and comparisons with SCGA predictions for NN (dashed black line) and FN (solid red line) models. For the NN model, J1=38K. (d)–(e) Selected reciprocal-space planes showing I0(Q) and I1(Q), as labeled in the figure, and comparison between NN and FN models calculated using the SCGA. Throughout, white rings are masked aluminum background lines. In (c)–(e), only Ei=80meV data are shown, but both 40 meV and 80 meV data were included in fits.

    Reuse & Permissions
  • Figure 2
    Figure 2

    Magnetic excitation spectra of MgCr2O4 at T=20K measured with incident neutron energy Ei=40meV, and comparison with linear spin-wave theory (LSWT) calculations for our FN model. (a) Momentum-energy slices through g2S(Q,E) along different paths, comparing data (left column) and FN model (right column). (b) Cuts at constant energies E¯±0.2meV through the data (gray circles) and FN model (red lines), where E¯ is labeled on each plot. The intensity is multiplied by E¯ and offset by 4sr1Cr1 for clarity. (c) Energy dependence of the experimental (colored circles) and modeled (colored lines) dynamical structure factor at selected momenta, normalized to the energy transfer E0=2meV. (d) Slices at constant energies E¯±0.2meV through the data (left column) and the FN model (right column) in the (h,k,k) plane. Throughout, blank space is due to kinematic constraints on the scattering, and the extra intensity at (4,0,0) arises from a strong nuclear Bragg peak and its associated acoustic phonon.

    Reuse & Permissions
  • Figure 3
    Figure 3

    (a) Sets of ordering wave vectors κ with nearly-degenerate energies represented in the Brillouin zone as colored surfaces. The pink surface shows wave vectors with energies within 0.5% of the global energy minimum for the FN interaction parameters of MgCr2O4. The blue lines show wave vectors with energy equal to the global energy minimum for the macroscopically-degenerate phases represented by a blue dashed line in the phase diagram. (b) Mean-field phase diagrams of our FN Heisenberg model as a function of J2 and J3a, showing results for J3b=0 (left) and J3b=0.0085J1 (right). Phases with different ordering wave vectors κ are shown in different colors. Special phases (dashed lines) correspond to a macroscopic number of ordering wave vectors with degenerate energies. They emerge from J3aJ2=0 (yellow line), which corresponds to the NN model, to form two half-planes corresponding to J3aJ2=J3b (green and blue lines).

    Reuse & Permissions
×

Sign up to receive regular email alerts from Physical Review Letters

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×