Single-ion properties of the Seff = 1/2 XY antiferromagnetic pyrochlores, 0 0 2+ 2+ NaA Co2F7 (A = Ca , Sr )

K.A. Ross,1, 2 J.M. Brown,1 R.J. Cava,3 J.W. Krizan,3 S. E. Nagler,4 J.A. Rodriguez-Rivera,5, 6 and M. B. Stone4 1Colorado State University, Fort Collins, Colorado, 80523, USA 2Quantum Materials Program, Canadian Institute for Advanced Research, Toronto, Ontario M5G 1Z8, Canada∗ 3Princeton University, Princeton, New Jersey 08544, USA 4Quantum Condensed Matter Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 5NIST Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA 6Materials Science and Engineering, University of Maryland, College Park, MD 20742, USA (Dated: March 23, 2017)

The antiferromagnetic pyrochlore material NaCaCo2F7 is a thermal spin liquid over a broad 2+ temperature range (≈ 140 K down to TF = 2.4K), in which magnetic correlations between Co dipole moments explore a continuous manifold of antiferromagnetic XY states.1 The thermal spin liquid is interrupted by spin freezing at a temperature that is ∼ 2 % of the mean field interaction strength, leading to short range static XY clusters with distinctive relaxation dynamics. Here we report the low energy inelastic neutron scattering response from the related compound NaSrCo2F7, confirming that it hosts the same static and dynamic spin correlations as NaCaCo2F7. We then present the single-ion levels of Co2+ in these materials as measured by inelastic neutron scattering. An intermediate spin orbit coupling model applied to an ensemble of trigonally distorted octahedral crystal fields accounts for the observed transitions. The single-ion ground state of Co2+ is a Kramers doublet with a strongly XY-like g-tensor (gxy/gz ∼ 3). The local disorder inherent from the mixed pyrochlore A sites (Na+/Ca2+ and Na+/Sr2+) is evident in these measurements as exaggerated broadening of some of the levels. A simple model that reproduces the salient features of the single- ion spectrum produces approximately 8.4% and 4.1% variation in the z and xy components of the g-tensor, respectively. This study confirms that an Seff = 1/2 model with XY antiferromagnetic exchange and weak exchange disorder serves as a basic starting point in understanding the low temperature magnetic behavior of these strongly frustrated magnets.

16,17 I. INTRODUCTION ical disorder, as has been observed in Yb2Ti2O7 and 18 Tb2Ti2O7. The recently discovered 3d transition metal pyrochlore Magnetism on the pyrochlore lattice is a rich field of fluoride NaCaCo F (NCCF) has previously been sug- study that encompasses many unusual magnetic phenom- 2 7 gested as a strongly interacting version of the antifer- ena such as spin ice, spin liquids, and Order by Disorder romagnetic (AFM) XY pyrochlore model,1 potentially (ObD).2 Experimental access to these phenomena has allowing for a detailed investigation of this model and been granted mainly by rare earth oxide pyrochlore mate- its variations at lower effective temperatures. The XY rials. These have the general chemical formula A B O , 2 2 7 anisotropy of effective spin 1/2 (S = 1/2) magnetic mo- with A a magnetic trivalent rare earth cation and B a eff ments in NCCF and the related compound NaSrCo F non-magnetic tetravalent cation such as Ti, Sn, Zr, Ge 2 7 (NSCF) are confirmed here through the modeling of Co2+ or Pt.2–5 Both the A and B sites independently form a single-ion levels by fitting to inelastic neutron scattering pyrochlore-type sublattice, a highly-frustrated three di- measurements. Furthermore, both materials are shown mensional network composed of corner-sharing tetrahe- to have nearly identical spin correlations in their spin dra (Figure 1 a). The anisotropic nature of the rare earth frozen states, consisting of short range ordered XY clus- magnetic moments, arising due to strong spin-orbit cou- ters and associated dynamics, despite differences in the pling (SOC) combined with crystal electric field (CEF) strength of chemical disorder. effects, has offered fascinating variations of magnetic be- havior. For instance, local Ising anisotropy is required for NCCF and NSCF belong to a family of 3d transi- 19–22

arXiv:1703.07468v1 [cond-mat.str-el] 21 Mar 2017 generating spin ice and its emergent magnetic monopoles tion metal fluoride pyrochlores. The chemical for- 6 in Ho2Ti2O7 and Dy2Ti2O7, while XY anisotropy can mula and structure are similar to the rare earth ox- 7–9 lead to ObD, as suggested for the pyrochlore ma- ide pyrochlores, taking the form A2B2F7, where A is 10–15 + 0 0 2+ 2+ terial Er2Ti2O7. However, the 4f electrons which a split site of Na /A (A = Ca , Sr ) and B is are responsible for this magnetism have weak interaction six-fold coordinated Co2+, Ni2+, Fe2+ or Mn2+ (Fig- strengths on the order of 1 K, requiring experiments to be ure 1). Unlike the previously studied fluoride-pyrochlore 23 done at millikelvin temperatures in order to access their CsNiCrF6, the magnetic B-site is chemically uniform, magnetic correlations. This can limit the exploration of hosting just one magnetic species. The average struc- the lowest temperature states in these frustrated systems. tures of these compounds as measured by x-ray diffrac- Furthermore, the ground states selected by these inter- tion show well-ordered pyrochlores (space group F d3¯m), actions can be quite sensitive to small amounts of chem- consistent with random distribution of cations on the A 2

2+ a) Co tetrahedra b) CoF6 octahedra oride pyrochlores. This is due to the quenching of or- bital angular momentum expected for most octahedrally + 2+ Na /A’ a coordinated 3d transition metals. This means that phe- Na+/A’ 2+ c E=[0.2, 5.0] meV, T=14 K nomenology like spinE=[0.2,5] ice meV, orT=[13,15] ObD K will likely not pertain to c) 70 b 4 most members of this series. The60 exceptions to this are the Co2+ compounds2 which we study50 here: NCCF and Co-F-Co bond 2+ 7 40 4 NSCF. The free0 ion Co (3d ) forms a F ground term angle variation 30 [H,0,0]

with S = 3/2[00L](r.l.u.) and L = 3, and when placed in a octahe- c) d) -2 20

10 Intensity (counts/mon=7.5e4) dral coordination-4 the CEF and SOC conspire to form an 0 26,27 Seff = 1/2 single-ion-3 -2 -1 ground0 1 2 3 state. A distorted oc- !2 !3 tahedral environment,[HH0][0,K,K] such (r.l.u.) as the average environment z = [111] of NCCF and NSCF, will lead to single-ion anisotropy 27–29 in the Seff = 1/2 states (Fig. 1 c). At tempera- y x tures low enough that the ground doublet states are the g-tensor only relevant degrees of freedom, it is then possible model the interactions between effective spin 1/2 operators, in which the single-ion anisotropy encoded by the g-tensor is FIG. 1. a) Depiction of the two cation sublattices in 28,30–32 02+ + projected into the effective exchange interactions. NaCaCo2F7 and NaSrCo2F7. The A = A /Na split site forms a pyrochlore sublattice (connectivity not shown here). Thus, it seems that NCCF and NSCF could serve as The nearest neighbor bonds in the B =Co2+ pyrochlore sub- new “high temperature” examples in which this success- lattice are depicted. b) Relationship between the A and B ful method, well-developed for the rare earth oxide series, sites and the average atomic positions. The A site cations re- could be used.12,31,33 side in the open hexagons of the B site pyrochlore structure, The static and dynamic spin correlations in as viewed from the h111i body diagonals. The Co-F bonds NCCF were previously measured by inelastic neu- are represented as rods. The yellow and blue lines illustrate tron scattering.1 It was found that below the freezing potential Co-F-Co bond angle variation with A cation occu- temperature (T = 2.4 K), short range order (SRO) pation at the central site. c) Illustration of the average local F CoF octahedra, which are trigonally distorted along the lo- of the spins develops with a correlation length of 16 6 ˚ cal h111i directions that point into the center of a pyrochlore A, corresponding to XY AFM configurations from the tetrahedron. The g-tensor is XY-like due to this environment. Γ5 irreducible representation of the tetrahedral point d) The two basis states forming the Γ5 manifold, from which group. This manifold is spanned by two basis states the static and dynamic correlations are selected in NaSrCo2F7 called ψ2 and ψ3, shown in Figure 1 d). In the XY and NaCaCo2F7. AFM pyrochlore model, the q = 0 long range ordered (LRO) magnetic structures based on these two states are accidentally degenerate,7–9,34 and either one can site. Despite the materials being well-ordered on average, be selected by various ObD mechanisms, as has been the local structure near the B cations will be distorted discussed at length in the context of Er2Ti2O7. The away from the average trigonal D3d point group due to observed selection of the non-coplanar ψ2 state in that the A site disorder, leading to variations in exchange pa- material35 has been argued to occur either via quantum rameters and single-ion anisotropy (Figure 1 b). In light and thermal fluctuations (ObD),10–13 or from a non- of this, perhaps it is not surprising that all of the mate- ObD mechanism involving virtual excitations to higher rials in this series eventually display spin freezing tran- crystal field levels.36 Quenched disorder in the form sitions. However, the spin freezing occurs at very low of dilution15 or exchange disorder14 has recently been effective temperatures; their Curie-Weiss temperatures predicted to select the coplanar ψ3 state instead, and 22 19 37 range from θCW ∼ -70 to -140 K, but the spins freeze this has been studied in yttrium-diluted Er2Ti2O7. In only at 2 to 4 K. This gives large frustration indices of NCCF, despite a clear mechanism for exchange disorder f = |θCW | ≈ 19 to 58, suggesting that the strength of and the presence of XY AFM correlations, an LRO state TF the disorder, the presumed cause of the spin freezing, is is not selected. weak compared to the overall interaction strength. An The low energy inelastic neutron scattering (INS) re- interesting comparison in this regard can be made to the sponse of NCCF shows non-dispersive, diffusive excita- B-site disordered Yb-based pyrochlore Yb2GaSbO7; with tions with a distinctive intensity vs. Q~ pattern, extending 1 its θCW = −1.15 K and lack of a freezing transition down to approximately E = 10 meV, as well as a broad dis- to the lowest measured temperature of 20 mK, this ma- tribution of relaxation times in the µeV range, as probed terial demonstrates similarly “ineffective” exchange dis- by NMR.38 These dynamic signatures persist above the order relative to its (much weaker) average interaction freezing temperature up to at least T =14 K (INS1) and strength (f ≥ 58).24,25 20 K (NMR38). These data, combined with the high Unlike the rare earth oxides, nearly ideal Heisenberg Curie-Weiss temperature (θCW = -140 K) relative to the moments with large spins should be expected for most freezing temperature (TF = 2.4 K) giving a frustration of the members in this new series of transition metal flu- index of f = 58, show that NCCF hosts an XY ther- 3 mal spin liquid, i.e., a strongly correlated but disordered ters were used. The resulting energy resolution was δE = state dominated by entropy, over a large temperature 0.17 meV at the elastic line. These data are compared to range. Surprisingly, NCCF resists ordering or freezing those previously published on a 0.87 g crystal of NCCF, to much lower effective temperatures than the canonical which was also taken using the MACS spectrometer in 39 1 XY pyrochlore Er2Ti2O7 (f = 20), suggesting that this a similar configuration. The relative masses cannot be material could be closer to a classical phase boundary in used to compare the intensities directly due to differences the general anisotropic exchange model developed for the in crystal mounting. The NSCF and NCCF data were rare earth oxides.33 To make further progress in modeling scaled relative to each other to match the inelastic inten- NCCF, it is crucial to establish the relevance of the gen- sities near (002) and (111) (Fig. 2). The same scaling eral anisotropic pseudo-spin 1/2 model. This relies upon also produced matching (220) magnetic intensities below the understanding of the single-ion Hamiltonian for Co2+ T = 1.7 K. in these materials. In order to investigate the Co2+ single ion levels in In this article we first present low energy inelastic NCCF and NSCF, higher energy INS experiments were and elastic neutron scattering measurements on the com- performed using the SEQUOIA time-of-flight chopper pound NSCF and compare it to the previously reported spectrometer at the Spallation Neutron Source, Oak measurements on NCCF, demonstrating that despite Ridge National Laboratory.41 The crystals were oriented their different A0 sites the two compounds share the with [HHL] in the horizontal scattering plane. Data were same experimental signatures and could be treated by taken at T = 5 K and 200 K, using incident energies of the same theoretical approach. We then present the ob- Ei = 60, 250, 700, and 2500 meV. This wide range of in- served single-ion levels for both compounds, measured cident energies was employed to probe the large dynamic by INS over the energy range of E = 30 - 1500 meV. range expected for single ion energy levels of Co2+. For The observed INS response as well as the dc magnetic these energies, the T0 Chopper and Fermi Chopper (FC1 susceptibility is well described by a disorder-averaged in- or FC2) speeds were set to the following values: for Ei termediate SOC model. This model confirms that the = 60 meV, 60 Hz (T0) and 420 Hz (FC2); for Ei = 250 single-ion ground state in both materials is a well-isolated meV, 120 Hz (T0) and 480 Hz (FC1); for Ei = 700 meV, Kramers doublet, i.e., Seff = 1/2, and provides the aver- 150 Hz (T0) and 600 Hz (FC1); and for Ei = 2500 meV, age g-tensors, which are observed to be strongly XY-like 180 Hz (T0) and 600 Hz (FC1). The elastic energy resolu- and have approximately 8% variation due to local disor- tions in these configurations were δE = 1.92, 12.71, 44.96 der. and 276.44 meV, respectively. Note that δE decreases The contents of the paper are as follows: in Section as the neutron energy transfer increases for each Ei, and II we give the experimental details of the INS measure- the full energy-dependent resolution for SEQUOIA41 was ments. In Section III we describe single-ion model used employed in our fits. to fit the high energy INS data. Section IV A demon- strates the equivalence of the spin correlations in the two 2+ materials, and presents detailed low energy INS measure- III. SINGLE-ION CALCULATIONS FOR CO ments at temperatures above TF for the first time (i.e., the thermal spin liquid regime). Section IV B details the In this section we present a method for calculating the single-ion results and fits. The results are further dis- single-ion levels in Co2+ within an intermediate spin orbit cussed in Section V, and conclusions from this study are coupling scheme, as well as the resulting INS response presented in Section VI. and dc magnetic susceptibility. This method was used to fit the data presented in Section IV B. The method is similar to that used for the rare earth based magnets,42,43 II. EXPERIMENTAL METHOD and results in similar phenomenology. The single ion levels for Co2+ in a trigonally distorted octahedron are Single crystals of NCCF and NSCF (space group Kramers doublets with anisotropic moments described by F d3¯m, room temperature lattice constants a = 10.4189 a g-tensor. A˚ and 10.545 A,˚ respectively) were grown via the Bridgman-Stockbarger method in an Optical Floating A. Intermediate spin orbit coupling Hamiltonian Zone furnace, as previously reported.19 We studied a 3.527 g single crystal of NSCF using the MACS spectrometer at the NIST Center for Neutron The single-ion Hamiltonian is approximated by two Research.40 The dynamic structure factor, S(Q,~ E), was contributions, measured in the [HHL] reciprocal lattice plane. Neu- trons with a final energy E = 3.7 meV were selected, f Hion = HCEF(L~ ) + HSOC(L,~ S~), (1) and post-sample BeO filters were used to remove higher harmonic contamination and reject neutrons for which where HCEF(L~ ) is the crystal electric field (CEF) Ef > 3.7 meV. For elastic scattering, a Be filter preceded Hamiltonian that acts only on the orbital angular mo- the sample, while for inelastic scattering no incident fil- mentum subspace, and HSOC is the spin orbit coupling 4

(SOC) term. The CEF Hamiltonian can be written in NCCF Elastic 1.7 K NSCF Elastic 1.7 K general as: a) 1500 Intensity Intensity (arb. units)

X 1000 HCEF = Bl,mOˆl,m (2) l,m [00L] r.l.u.

44–46 where Oˆl,m are the Stevens operator equivalents. For transition metal ions, these operators are written in [HH0] (r.l.u.) [HH0] (r.l.u.) terms of the orbital angular momentum matrix opera- b) 60 80 ˆ ˆ ˆ NSCF NSCF tors L+, L− and Lz (as opposed to the total angular mo- NCCF NCCF 50 60 mentum operators Jˆ+, Jˆ− and Jˆz relevant for f-electron systems).45 For the trigonal point group symmetry rele- 40 2+ 40 vant to the average local environments of Co in NCCF 30 and NSCF, (l, m) = (2, 0), (4, 0) and (4, 3) are the only (002) 20 (111) Intensity (arb. units) 45 20 Intensity (arb. units) nonzero terms. The Bl,m values can be extracted by 1.7 K 1.7 K

fitting to INS data. 10 0 0 2 4 6 8 10 0 2 4 6 8 10 The spin orbit coupling term is given by, E (meV) E (meV)

H = pλ(S~ · L~ ) = pλ(Sˆ Lˆ + Sˆ Lˆ + Sˆ Lˆ ), (3) FIG. 2. Comparison of neutron scattering intensities from SOC x x y y z z NSCF and NCCF. Data taken on the MACS spectrometer with energy resolution δE = 0.17 meV. a) Magnetic elastic where for the free ion Co2+ λ = −22.32 meV,47 and scattering in NaCaCo2F7 (left) and NaSrCo2F7 (right) at T p is the “orbital reduction parameter” that can be used = 1.7 K. Intensities from symmetrically equivalent quadrants to account for changes in effective SOC strength due of the scattering plane were averaged. The sharp feature near to covalency.26 The x and y spin operators are linear (002) in NSCF is due to multiple scattering and is not repre- ˆ ˆ ˆ 1 ˆ ˆ ~ combinations of S+ and S− as usual: Sx = 2 (S+ + S−) sentative of spin correlations. b) Constant-Q scans near the 1 (002) and (111) positions. Errorbars represent one standard and Sˆy = (Sˆ+ − Sˆ−) and similarly for Lˆx and Lˆy. 2i deviation.

The full single-ion Hamiltonian (Eqn. 1) can be diago- B. Comparison to inelastic neutron scattering nalized within the 28×28 manifold of states formed by the 4 |Lz,Szi basis of the S = 3/2 and L = 3 free ion term ( F ) ~ for Co2+. A trigonally distorted octahedral coordination, The dynamic structure factor, S(Q, E), is related to as found in NCCF and NSCF, produces 14 Kramers dou- the observed intensity of inelastic neutron scattering via kf 2 I(Q,~ E) = f(|Q|) S(Q,~ E), where ki and kf are the blets. The presence of a doublet ground state permits the ki description of the magnetic moments as pseudo-spin 1/2 initial and final wavenumbers of neutrons scattered with energy transfer E = ~ (k2 − k2), and f(|Q|) is the mag- (Seff = 1/2) at sufficiently low temperatures. This dou- 2m i f blet can be described as an anisotropic magnetic moment netic form factor.48 The single-ion dynamic structure fac- with strength given by the g-tensor.27 In the trigonal tor at constant |Q| associated with Eqns. 1 to 3 can be symmetry appropriate to NCCF and NSCF, the g-tensor calculated as follows:49 has two independent components, gz and gxy, which re- fer, respectively, to the local [111] direction pointing into −βE 2 n Γ 0 |hν |Lˆ + 2Sˆ |ν 0 i| the center of the tetrahedron, and the plane perpendicu- X e hn,n i n α α n S(E) = C 2 2 , (4) Z ([En0 − En] − E) + Γ 0 lar to that (Fig. 1c). The g-tensor components are given n,n0,α hn,n i by the matrix elements of the magnetic moment operator 0 in the subspace of the two ground state wavevectors (|ν1i where α = x, y, z, while n and n label eigenstates and |ν2i), of Eqn. 1. C is a scale factor applied to match the measured intensity of constant |Q| data (note that at constant |Q| the form factor f(|Q|)2 can be absorbed ˆ ˆ P gxy = −2hν1|(Lx + 2Sx)|ν2i, into C). The partition function is Z = n exp (−βEn) with β = 1/kBT . |νni is the wavevector from the sin- and gle ion calculation corresponding to energy eigenvalue En. A Lorentzian half width at half maximum (HWHM), ˆ ˆ gz = 2hν2|(Lz + 2Sz)|ν2i. Γhn,n0i, accounts for line broadening due to instrumental resolution, finite excitation lifetimes, or averaged disper- The resulting saturated magnetic moment sizes in the sion of the transitions between state n and n0. To fit gz gxy principal directions are µz = 2 µB and µxy = 2 µB. the data presented in Section IV B, we set Γhn,n0i to the 5 energy transfer-dependent instrument resolution for all IV. RESULTS transitions except those involving the ground doublet, n = 1, 2, to the first excited doublet n0 = 3, 4. In that case, Γhn,n0i = 13.5 meV accounts for the increased width due to the finite dispersion of the first excited level (Fig. A. Low energy dynamic structure factor in 4). An elastic peak with a full width at half maximum NaSrCo2F7 (FWHM) fixed at the instrumental resolution was also added to account for all sources of nuclear elastic scat- tering, both coherent and incoherent. To account for disorder in the local environment of The low temperature (T = 1.7 K) and low energy Co2+ brought about by the randomly mixed A site in (E < 5 meV) magnetic neutron scattering response from the crystal structures of NCCF and NSCF, the average NSCF is compared to that from NCCF in Figure 2. The S(E) can be calculated over an ensemble of ions with NCCF data was previously published in Ref. 1. For the slightly varying Bl,m parameters in Eqn. 2. If the lo- elastic scattering presented for both materials, the equiv- cal site symmetry remains trigonal, as expected based alent high temperature (T = 14 K > TF ) elastic maps on results of vibrational spectroscopy from the related were subtracted in order to isolate the magnetic scatter- 50 compound NaCaMg2F7, there remain only three pa- ing. For the inelastic scattering, an empty can subtrac- rameters that enter into HCEF. For simplicity we chose tion was performed. The neutron scattering intensity a linear distribution of the B4,0 parameter, allowing the of NCCF has been scaled such that the (220) magnetic width of the distribution (∆B4,0) to be fit. The reasoning elastic peak intensity matches NSCF, which also pro- for this choice is described in Section V. duces equivalent intensities for the inelastic scattering. The neutron scattering patterns are nearly identical over the whole energy range, indicating that the static and dynamic spin correlations are equivalent in both com- pounds. An interpretation of this scattering in terms of C. Magnetic Susceptibility static and dynamic XY correlations has been presented in Ref. 1; this interpretation is expected to qualitatively ap- The eigenvectors and eigenvalues of Eqn. 1 also per- ply to both materials. The sharp features near the (002) ¯ mit the calculation of the (powder averaged) single ion and (002) positions in NSCF seen in Fig. 2 a) appear magnetic susceptibility as a function of temperature,49,51 in the elastic channel using several choices of incident energies, but can be made to vanish for other choices, signifying that their origin is multiple scattering. Such (002) scattering has been observed before in rare earth

2 P ˆ ˆ 2 −En/T titanates and was previously thought to be caused by NAµ X |hνn|Lα + 2Sα|νni| e χ (T ) = B n symmetry lowering, but has recently been identified as ion 3k Z T B α multiple scattering in those materials as well.52 −En/T −Em/T  X X 2 e − e + |hνm|Lˆα + 2Sˆα|νni| . Em − En n m6=n The data shown in Fig. 2 was collected below the (5) freezing temperatures of NCCF and NSCF (2.4 K and 3.0 K, respectively). Above the freezing transition, it was shown previously in NCCF that the elastic magnetic This expression includes the van Vleck susceptibility, scattering vanishes, while inelastic magnetic scattering at but excludes the diamagnetic susceptibility χdia, and E = 0.5 meV persists to at least T = 14 K and retains the does not account for the mean field interaction between same distinctive pattern in reciprocal space.1 In NSCF we magnetic moments. To compare to the measured dc sus- have confirmed the persistence of this inelastic scattering ceptibility, the following mean field approximation can 30 above TF at T = 14 K. The inelastic scattering at E = be used, 0.25 meV is broader at 14 K compared to 1.7 K, implying shorter ranged dynamic correlations (Fig. 3 a). The full energy dependence of the scattering up to 5 meV is χion compared for the two temperatures in Fig. 3 b), which χMF = χdia + . (6) 1 + λW χion shows that the spectrum has similar features above and below TF , but with enhanced low energy spectral weight at T =14 K. The energy integrated scattering at T = 14 Here, λW is the Weiss molecular field constant, which K in NSCF (E = 0.2 to 5.0 meV) over the [HHL] plane is accounts for the mean exchange interactions. With this shown in Fig. 3 c). This scattering could be compared to sign convention, a positive value of λW indicates AFM equal time correlation functions, such as those calculated interactions. in Ref. 33. update panel c

6

a) E=0.25 meV, [HH0]=[1.9, 2.1] b) c) E=[0.2,6]E=[0.2, meV, 5.0] T=[13,15]meV, T=14 K K 100 70 1.7 K 1.7 K (002) 4 150 14 K 1.7 K (111) 60

80 14 K (002) Intensity (arb. units) 14 K (111) ξ(1.7 K) = 11.2 Å 2 50 ξ(14 K) = 6.96 Å 60 100 40 0

40 [H,0,0] 30 [00L](r.l.u.) 50 -2 20 Intensity (arb. units) Intensity (arb. units) 20 10 -4 Intensity (counts/mon=7.5e4) 0 0 0 -2 -1 0 1 2 0 2 4 6 8 10 -3 -2 -1 0 1 2 3 [22L] (r.l.u.) E (meV) [HH0][0,K,K] (r.l.u.)

FIG. 3. Temperature dependence of inelastic scattering in NSCF. a) Constant energy cuts at 0.25 meV along [22L] reveal how the correlation length of the low energy dynamics changes with increasing temperature. The dynamic correlation length along this direction decreases from 11.2 A˚ to 6.96 A˚ upon heating (dotted lines are fitted Lorentzians). b) Constant Q cuts near the (002) and (111) Bragg positions at T=14K and 1.7K. c) Diffuse inelastic scatteringSet in to NSCF 0.0014 above for NCCF the freezing temperature, at T =14 K (data folded and symmetrized, integrated from E = 0.2 to 5.0 meV).set to 0.0014/1.19 Errorbars represent= 0.00117 one for NSCF standard deviation. old

Note: the maximum on the color scale for NSCF was divided by [K,-K,0]=[-0.1,0.1],[K,-K,0]=[-0.1,0.1], [H,H,0]=[-2.1,-1.9] [H,H,0]=[-2.1,-1.9] [K,-K,0]=[-0.1,0.1], [H,H,0]=[-2.1,-1.9] B. Single ion levels and average g-tensors of 0.0014 0.0014 1.19 to match the scaling in the line cuts of the CEF figure (panel c). This 5050 5K 50 200K NaSrCo2F7 and NaCaCo2F7 0.00120.0012 0.00120.8 was done to match elastic lines. NCCF NCCF 4040 40 0.00100.0010 0.0010 I copied and pasted the colorbar from NCCF to NSCF. Also 0.00080.0008 2+ 3030 30 0.00080.6 redefined max # since it is clearer this way. Transitions between single ion levels of Co were ob- Intensity (arb. units) 0.0006 E (meV) 0.0006 E (meV) 0.0006 2020 20 served in NCCF and NSCF using the SEQUOIA time E (meV) 0.0004 [K,-K,0]=[-0.1,0.1], [H,H,0]=[-2.1,-1.9] Intensity (arb. unit) 0.00040.4 [K,-K,0]=[-0.1,0.1], [H,H,0]=[-2.1,-1.9] 0.0004 Intensity (arb. unit) [K,-K,0]=[-0.1,0.1], [H,H,0]=[-2.1,-1.9][K,-K,0]=[-0.1,0.1],Intensity (arb. unit) [H,H,0]=[-2.1,-1.9] of flight neutron spectrometer. The energies of the ob- 1010 10 0.00020.0002 0.0002 50 5K 50 0.0014200K 1.40.0014 [K,-K,0]=[-0.1,0.1],[K,-K,0]=[-0.1,0.1], [H,H,0]=[-2.1,-1.9] [H,H,0]=[-2.1,-1.9] [K,-K,0]=[-0.1,0.1], [H,H,0]=[-2.1,-1.9] served transitions were found range from 28 to 908 meV, 500 50 0.00000.0012 0.00000.00.0012 NCCF 0.0012NCCF 0.0012 -6 -4 -2 0 2 4 6 0.0000 0.0014 40 40 -6 -4 -2 0 2 4 6 1.0 requiring the use of multiple incident energies. The first 5050 -6 -4 -2 [0,0,L]05K2 4 6 50 200K[0,0,L] 0.0010 0.0010 0.00120.00100.0010 0.00100.00120.80.0010 40 [0,0,L] 40 30 30 0.0008 0.0008 NSCF NSCF Intensity (arb. units) excited level measured using Ei = 60 meV in both NCCF 4040 40 0.0010 0.0010 0.00080.0008 0.00080.0008 20 20 0.0006 0.60.0006 E (meV) 30 30 E (meV) and NSCF is shown in Figure 4, at both T = 5 and 200 0.0008 0.0008 Intensity (arb. unit) 30 30 0.6 0.0004 0.0004 Intensity (arb. unit) 30 0.00060.0006 0.00060.0006 10 10 0.0006 0.2 K. In both materials this level is located near 30 meV, 20E (meV) 20E (meV) 0.0006 0.0002 0.0002 E (meV) 2020 E (meV) 20 0.0004 0.0004 E (meV) 0.0004 0.0004 0[K,-K,0]=[-0.1,0.1], [H,H,0]=[-2.1,-1.9]0 [K,-K,0]=[-0.1,0.1], [H,H,0]=[-2.1,-1.9] Intensity (arb. unit) 0.0004 Intensity (arb. unit) 0.0000 0.0000 0.0004 Intensity (arb. unit) 0.4 Intensity (arb. unit) persists to temperatures above |θCW | (demonstrating its 10 10 Intensity (arb. unit) -6 -4 -2 0 2 4 6 -6 -4 -2 0 2 4 6 1010 10 5K 200K 1.40.0014 0.00020.00020.0002 0.00020.0002 50 50 0.0012 0.0012 single-ion origin), and displays the characteristic decreas- [0,0,L]NSCF [0,0,L]NSCF 0.0012 00 0 0.00000.0000 0.00000.00.0000 40 40 0.0010 0.0010 ing intensity vs. |Q| dependence of magnetic excitations -6 -4 -2 0 2 4 6 -60.0000-4 -2 0 2 44 66 1.00.0010 -6-6 -4 -2-2 00 224 46 6 -6 -4 -2 0 2 4 6 0.0008 0.0008 (Fig. 5a)). At T = 5 K the first excited level displays [0,0,L] [0,0,L] 30 30 0.0008 [22L][0,0,L][0,0,L] (r.l.u.) [22L][0,0,L] (r.l.u.) 0.0006 0.0006 20 20 0.60.0006 distinct dispersion, suggesting exchange-induced mixing E (meV) E (meV) 0.0004 0.0004 Intensity (arb. unit) 0.0004 Intensity (arb. unit) 26 10 10 of the ground state and first excited level. Note that the 0.0002 0.20.0002 low energy spin excitations presented in Section IV A are FIG. 4. The first excited single ion levels in NCCF (top row) 0 0 0.0000 0.0000 also identifiable in Fig. 4 below ∼ 10 meV; these are due and NSCF (bottom row), observed by inelastic neutron scat- -6 -4 -2 0 2 4 6 -6 -4 -2 0 2 4 6 [22L][0,0,L] (r.l.u.) [22L][0,0,L] (r.l.u.) to the magnetic correlations discussed in Ref. 1, rather tering (Ei = 60 meV). Data are averaged over the following than single-ion levels. ranges in perpendicular reciprocal lattice directions (r.l.u.): [HH0], H= [1.9, 2.1] and [K-K0], K= [-0.1, 0.1]. The full set of observed single-ion levels in NCCF at T = 5 K are shown in Figure 5 along with results of fit- ting the model described in Section III. Some of these features are significantly broader than the energy reso- excited level, and this dispersion is well-accounted for in lution as expected, since the local variation in Co envi- the model by an increased FWHM, as discussed below. ronment should produce finite energy ranges over which For a uniform Co2+ environment, 14 levels (including transitions are observed. The measured energies of the the ground state) are expected in total, each of them single-ion levels and corresponding energy widths are tab- doubly degenerate according to Kramers theorem. Note ulated in Table I for both NCCF and NSCF.53 All neu- that at T = 5 K, a low temperature compared to the tron scattering data in Fig. 5 have a corresponding back- energy of the first excited level (∼ 300 K), only transi- ground subtracted and are corrected for the kf /ki factor tions from the ground state to the excited states will be described in Section III B, thus representing the quan- observable. Additionally, some of these transitions may tity f(|Q|)2S(Q, E) in arbitrary units. The orientational have weak intensities, depending on the strength of the average of these single crystal data are presented in or- transition matrix elements between eigenstates of Eqn. 1. der to provide a concise overview of the excitations; no The intensities are weighted by the square of the tran- significant dispersion was detectable in any but the first sition matrix elements of the magnetic moment opera- P = -25.7809288710721 -0.455513285837226 -14.0096981480374 1.02146454540241 0.126644904813703 5001.33565514524 23.4850001834351 errP = 0.523850994443288 0.00530897134718740 0.0705415649952567 0.00413396271890303 0.00358858313503159 71.6050101069163 0.2133307311537947

50 a) NCCF Ei = 250 meV b) c) 50 n=3 Ei = 250 meV n=5 Ei = 250 meV

100 (arb. unit) Log(Intensity) 200 40 40 n=7 meang = 150 n=9 30 30 n=11 1.8702 10 6.0780 100 6.0780 NCCF

E (meV) 20 NSCF 20 50 Intensity (arb. unit) 8.4506% variation in gz S(E) (arb. units) -1 4.1285% variation in gxy Intensity (arb. units) 10 10 Q=[4-5] Å 0 1 calculation 2 4 6 8 10 12 -1 0 0 |Q| (Å ) 0 50 100 150 0 50 100 150 200 E (meV) E (meV) d) 3 e) 2 f) 100 Ei = 700 meV n=21 Ei = 2500 meV n=23 measured 2.5 364 ± 1.3 meV 908 ± 27 meV calculated 1.5 80 0.2 2 0.1 Oe/ emu) 1.5 1 " 60

I (arb. units) 0 n=13 1 Q=[5-14]Å-1 500 1000 1500 (molCo

S(E) (arb. units) n=15 S(E) (arb. units) 0.5 E (meV) -1 @ 40 n=17 n=25 Q=[5-15]Å 1/ 0.5 n=27 -1 n=19 Q=[15-25]Å NCCF 0 0 20 300 350 400 450 500 600 800 1000 1200 1400 1600 1800 0 100 200 300 E (meV) E (meV) T (K)

2+ FIG. 5. Measurements and fits of single ion levels of Co in NaCaCo2F7. Panels a) through e) show the observed single-ion levels from orientationally averaged single crystal inelastic neutron scattering measurements (SEQUOIA spectrometer), at T = 5 K. Shaded regions in panels c) to e) indicate the energy ranges spanned by the doubly-degenerate eigenvalues of Eqn. 1; the eigenvalues vary for different local environments, producing a range of values (see main text). The first eigenstate label, n, in the corresponding doublet is indicated for each range. a) Energy vs. |Q| pseudo-color plot, with Ei = 250 meV. Intensity is displayed on a logarithmic scale. Arrows indicate the positions of relatively sharp magnetic modes, and the rectangle indicates a broad band of magnetic scattering. b) Cut through data in panel a), and equivalent data from NSCF, averaging from |Q| −1 = 4 - 5 A˚ (empty can background subtracted). c) Fit to the NCCF data shown in panel b). d) Ei = 700 meV, averaging −1 from |Q| = 5 − 14 A˚ . e) Ei = 2500 meV, cuts taken at both high and low |Q|. On top of a steep non-magnetic background, additional intensity is seen for low |Q| near 900 meV. The inset shows the subtraction of the high-|Q| background, which gives a peak at 908 ± 27 meV. f) Single-ion susceptibility with mean field interactions (Eqns. 5 and 6) compared to the measured dc susceptibility of NCCF from Ref. 19.

2 P ˆ ˆ 2 tor, Mhn,n0i = α |hνn|Lα + 2Sα|νn0 i| , which enter into at 364 ± 1 meV and 908 ± 30 meV. In these high energy Eqn. 4. ranges, multi-phonon scattering creates a strong, sloping background and obscures these magnetic features. Nev- Fig. 5 a) shows measurements at Ei = 250 meV, with ertheless, their magnetic nature is confirmed by their |Q|- the intensity presented on a logarithmic color scale. This dependence. For the data shown in panel d), the single- shows several of the magnetic excitations, including the ion level energy was extracted using a fit that included first excited level near 30 meV (as in Fig. 4), a broad a sloping background plus a Gaussian. For the 908 meV band of magnetic features spanning approximately 46 to feature, intensity vs. E for two |Q| ranges were com- 100 meV (indicated by a rectangle), and two higher en- pared, as shown in panel e). The precision of the energy ergy features near 120 and 140 meV. Note that phonon determination for these higher modes is relatively poor scattering is also visible below 30 meV, but its intensity due to decreased energy resolution of the higher incident increases as function of |Q|, while the magnetic features energy instrument configurations. show the opposite trend in intensity vs. |Q| due to the magnetic form factor (see also Appendix A). Intensity vs. Fits to these measured single-ion levels were carried energy cuts from Q = [4, 5] A˚−1 are compared for NCCF out using a least squares minimization routine. The dy- and NSCF in Fig. 5 b), showing the overall similarity namic structure factor at constant |Q|, i.e, S(E), was between the single ion levels in the two compounds. The calculated from the single-ion model with local trigonal first excited level is at a slightly higher energy for NSCF disorder described in Section III, and was compared to −1 (29.31(2) meV) compared to NCCF (28.05(2) meV). Two the Ei = 250 meV data averaged from |Q| = 4 − 5 A˚ higher energy levels are shown in Fig. 5 d) and e), located (Fig. 5 c)). The measured energies and uncertainties of a) b) 1200 200 n = 27,28 n = 25,26

1000 n = 11,12 n = 23,24 150 n = 21,22 n = 9,10 800

n = 7,8 600 100 n = 5,6 E (meV) n = 19,20 E (meV) n = 17,18 400 n = 15,16 n = 13,14 50

200 n = 11,12 n = 9,10 n = 3,4 n = 7,8 n = 5,6 n = 3,4 0 0 0.35 0.4 0.45 0.5 0.55 0.35 0.4 0.45 0.5 0.55 B (meV) B (meV) 4,0 4,0 8

NCCF NSCF a) b) 1200 200 Energy (meV) FWHM (meV) Energy (meV) FWHM (meV) n = 25-28 1000 n = 11,12 n = 21-24 150 a a n = 9,10 28.05(2) 27.50(7) 29.31(2) 29.53(8) 800

a a n = 7,8 ∼ 46 - 101 broad ∼ 41 - 100 broad 600 100 n = 5,6 E (meV) a a ∗ E (meV) 121.9(1) 9.1(4) 120.7(2) 7.2(9) 400 n = 13-20 50 a a 200 ∼ 139 broad ∼ 134 broad n = 3,4 n = 3-12 364(1)b 24(3) – – 0 0 0.35 0.4 0.45 0.5 0.55 0.35 0.4 0.45 0.5 0.55 B (meV) B (meV) 9.1(3)×102 c 3.2(6)×102 9.0(3)×102 c 3.7(6)×102 4,0 4,0

TABLE I. Measured single-ion levels in NaCaCo2F7 (NCCF) FIG. 6. a) Energy eigenvalues as a function of B4,0. The and NaSrCo2F7 (NSCF), with central peak positions and full eigenstates are labelled by n. b) Expanded view of the low width at half maxima (FWHM). Due to time constraints, the energy part of panel a). peak expected near 360 meV in NSCF was not measured. The quoted error on the peak centers is the standard deviation of the fitted peak center. The FWHM are marked with an as- Parameter Fitted value terisk (∗) when they are resolution-limited. Instrument con- B2,0 25.8(5) (meV) figurations are indicated by a : Ei = 250 meV, b : Ei = 700 meV, c : Ei = 2500 meV. B4,0 0.455(5) (meV)

B4,3 14.01(7) (meV) the higher excited levels (364±1 meV and 908±30 meV) ∆B4,0 0.127(4) (meV) were included as constraints on the minimization, but no p 1.021(4) % variation attempt was made to compare their relative intensities, g 1.87 8.4 % which are difficult to accurately determine due to instru- z mental effects. The HWHM parameters Γhn,n0i in Eqn. gxy 6.08 4.1 % 4 were held fixed to the calculated instrumental resolu- tion for most values of n and n0, but were increased to TABLE II. Parameters for the single-ion Hamiltonian, fit to 13.5 meV for n = 1, 2 and n0 = 3, 4 to account for the inelastic neutron scattering data in NCCF. The values of gz dispersion of the first excited level. The coherent and and gxy are those obtained from the mean value of the CEF parameters (i.e., when ∆B = 0). incoherent nuclear elastic scattering was modeled by a 40 Gaussian with FWHM fixed to the instrumental resolu- tion (δE = 12.7 meV), and the area of that elastic signal NCCF from Ref. 19, using the average B parameters was used as an additional fitting parameter (C ). In to- l,m E shown in Table II. Fitting the inverse susceptibility from tal there were seven fitting parameters; the three average T =100-300 K using Eqns. 5 and 6, the two parame- crystal field parameters (B , B , and B ), a width 2,0 4,0 4,3 ters χ = −3.1(1) × 10−4 emu/mol and λ = 23.99(5) of the linear distribution of B (∆B ), the orbital re- dia W 4,0 4,0 mol/emu were determined, and excellent agreement is duction parameter p, the elastic line area C , and an E found with the high temperature data (Fig. 5 f). This overall scale factor C. These fitted parameters and the agreement further confirms the validity of the single-ion associated average g-tensors with their percent variation scheme determined from the neutron scattering results. are listed in Table II. The energy eigenvalues for the av- erage Bl,m, their variation over the full ∆B4,0 range, as well as transition matrix elements, are listed in Table III. V. DISCUSSION The energy levels as a function of B4,0 are visualized in Figure 6. Note that due to the overall similarity between the single-ion levels of NCCF and NSCF, as well as the 1. Single-ion model limited accuracy of the simple model of local disorder (as discussed in Section V), only one set of parameters are The single-ion fits suggest that the reduction of spin presented. Although this fit is not likely to be unique in orbit coupling (quantified by p) is not significant, since it the sense that the broad observed transitions and mini- refines to 1.021(4), which is nearly 1.0. This is in contrast 26 mal modeling of disorder prevent a precise determination to KCoF3, for which p =0.93 ), indicating orbital mo- of the parameters, the obtained parameters adequately ment reduction due to covalency. Thus, naively, NCCF capture the main features of the single-ion levels in both and NSCF seem to have strongly ionic bonds. materials and provide the key physical insights that can The model used to fit the single-ion levels, presented in be gained from these data. Section III, accounts for the local disorder in a minimal Finally, the high temperature dc magnetic susceptibil- way. The local symmetry at the Co2+ site is assumed to ity can be compared to the measured susceptibility of remain trigonal, based on vibrational spectroscopy of the 9

2 2 State (νn0 ) En0 for ∆B4,0 = 0 Range of En0 Mh1,n0i+Mh2,n0i ∼ 40 and 100 meV in INS, without creating excessive widening (beyond the measured widths) of the higher ν ,ν 0.0 0.0 19.48 1 2 energy levels. It may be noted that this model does not ν3,ν4 27.7 3.5 17.06 perfectly account for the intensities and positions of all features, particularly near the modes near 121 meV. This ν5,ν6 69.4 36.1 3.55 is not unexpected, as the model we are using for the dis- ν7,ν8 90.0 25.4 0.37 order is simplified. The numerical values we report for ν9,ν10 135.4 3.03 0.12 the g-tensor and its variation (gz = 1.87 ± 8.4 % and ν ,ν 154.3 13.0 0.04 gxy = 6.08 ± 4.1%) can nevertheless be used as good es- 11 12 timates. ν13,ν14 355.6 22.6 0.51 A more accurate model of the disorder could perhaps ν ,ν 367.6 28.6 0.67 15 16 be developed based on additional experimental informa- ν17,ν18 389.3 35.6 0.83 tion about the local structure of these materials. To this 23 ν19,ν20 411.1 41.1 0.17 end, a Na NMR study of NCCF has revealed an unusual high temperature response suggesting that there are two ν21,ν22 863.9 13.5 0.11 main crystallographic environments for Na.38 So far no ν23,ν24 896.7 12.2 0.71 model to explain this observation has been proposed. A measurement of the neutron or x-ray Pair Distribution ν25,ν26 1103.9 16.5 0.02 Function (PDF) could provide additional information. ν27,ν28 1133.3 15.8 0.01 Although there remains some uncertainty in the dis- TABLE III. Calculated energy levels (in meV) from the fit tribution of g-tensor values, the results of our single- shown in Fig. 5, and the corresponding energy ranges arising ion study unequivocally indicate that NSCF and NCCF from the variation in a crystal field parameter (∆B4,0) which are XY pyrochlores with well-isolated Kramers doublets. 2+ was used to model the local disorder at the Co site. The This doublet nature of the low energy degrees of freedom sum of the square of the transition matrix elements from the is consistent with estimates of the spin entropy from spe- ground doublet to each excited doublet are listed, for the av- 19,20 2 P 2 cific heat which reach nearly Rln2 per mole of Co, as erage CEF parameters (M 0 = |hν |Lˆ + 2Sˆ |ν 0 i| ). hn,n i α n α α n well as with the sum rule analysis of low energy neutron Note that not all of the transitions listed in this table have 1 been observed; refer to Fig. 5 for the correspondence between scattering intensities. Our results show that interactions the measured and calculated peak positions. in these materials may be treated within a pseudo-spin 1/2 model with XY effective exchange. If the domi- nant interactions occur between nearest-neighbors, this non-magnetic structural analog NaCaMg F from Ref. would produce a model similar to the one used to de- 2 7 11,12,32 50. In that study it was found that the number and se- scribe Er2Ti2O7 and Er2Sn2O7. The variation in lection rules of vibrational modes was inconsistent with the g-tensor implies that exchange disorder is needed in the full F d3¯m symmetry of the average structure, but the effective model. could be accounted for using the F 43¯m subgroup. Using the online software package isodistort,54 and assum- ing atomic displacements consistent with this symmetry lowering, we find that the point group of the Co site re- mains trigonal, but is reduced from D3d to C3v (note that on average, these local distortions would cancel out and 2. Susceptibility produce the F d3¯m space group that is observed for the average structure). This implies that only three CEF parameters are required, but that there may be some The dc susceptibility in both NSCF and NCCF can variation in their values from site to site, as the Co-F be well-accounted for by the single-ion model using with bond angles and distances (and therefore the CEF) are the average Bl,m parameters reported in Table II. The modified by the local disorder. In the absence of a more model naturally explains the high effective moments of realistic model for the distribution of local disorder, we these two materials (∼ 6 µB) as determined from Curie- initially allowed a linear variation in each of the crystal Weiss fits of the inverse susceptibility over the range 100 field parameters in Eqn. 2 and averaged the calculated K to 300 K.19 Another feature of the model is a downturn S(E) for each set of parameters. We found that only in 1/χ near 50 K; this experimentally observed feature a variation in B4,0 (which we call ∆B4,0) was necessary was previously assumed to be due to ferromagnetic cor- to reproduce the observed widths and positions of the relations, but it can now be understood as changes in the peaks observed via INS (Fig. 5). Specifically, the ener- thermal population of excited single-ion levels. The fitted gies of the ν5, ν6 and ν7, ν8 doublets are more widely Weiss molecular field constant, λW = 23.99(4) mol/emu, distributed compared to the other levels, consistent with confirms that AFM interactions are dominant, in agree- the broad band of magnetic intensity observed between ment with the low energy INS results. 10

3. Comparison of NaCaCo2F7 and NaSrCo2F7 VI. CONCLUSIONS

The neutron scattering results presented in Section We have presented inelastic neutron scattering results IV indicate that the magnetic correlations and single on the pyrochlore materials NaCaCo2F7 and NaSrCo2F7. ion properties of NSCF and NCCF are nearly identical. The low energy response of NaSrCo2F7 (below 10 meV) NCCF and NSCF have freezing temperatures of 2.4 K confirms that it hosts the same type of spin correlations and 3.0 K, respectively. The increased freezing temper- as NaCaCo2F7, which are fully dynamic above the freez- ature of NSCF seems likely to be related to the larger ing temperatures (2.4 - 3.0 K) but static and short range differential ionic radius on the A site: for 8-fold coordi- correlated below. These correlations are well-described in nated Na+, r = 1.18 A,˚ for Ca2+ r = 1.12 A˚ and for terms of static and dynamic antiferromagnetic XY spin 2+ 55 Sr r = 1.26. The bigger this size difference, the more clusters, as discussed in Ref. 1. Below TF , the XY config- severe the local disorder is expected to be. The measured urations are frozen, short range correlated versions of the single-ion levels are very similar, with the only notable ψ2 and ψ3 states which are known to be selected through difference being the position of the first excited level, various types of order-by-disorder in the XY AFM mod- which is slightly higher in energy for NSCF compared to els of the pyrochlore lattice. NCCF. However, within the accuracy of our model, there Our main result is the measurement and analysis of is no discernible difference in the single-ion ground state the single-ion levels of Co2+ in these materials. These wavefunctions. Thus, it seems likely that both NSCF and are well-modeled by an intermediate spin orbit coupling NCCF could be treated within the same low energy effec- Hamiltonian. We incorporated a distribution of crystal tive model, but the strength of exchange disorder should field parameters to account for local disorder brought be enhanced for NSCF. about by the split non-magnetic A site. The ground state of Co2+ in these environments is always a Kramers doublet with an XY g-tensor. The average crystal field parameters produce gz = 1.87 and gxy = 6.08, with a variation of 8.4% and 4.1%, respectively. The single- 4. Co2+ pyrochlores ion ground state doublets are separated from the first excited states by ∼300 K. Thus, at low temperatures (T < θCW ∼ 150 K), where spin correlations develop, a Given the promise of extending the phenomenology low energy effective theory built from S = 1/2 operators of anisotropic rare earth pyrochlores to materials with can be used to describe these materials. Due to the XY higher interaction strengths, searching for other Co2+- g-tensor and its variation, the effective exchange interac- based pyrochlores seems desirable. To the best of our tions are likely to be XY-like with bond (exchange) disor- knowledge, only one other material has been reported der. A basic model of the AFM XY pyrochlore with ex- 2+ with a Co pyrochlore sublattice: the spinel GeCo2O4. change disorder has been predicted to lead to a long range This material has been discussed in terms of its Seff = 1/2 ordered state. However, NaCaCo2F7 and NaSrCo2F7 do single-ion ground state, and some unusual details of the not seem to conform to this prediction. higher energy excitation spectrum have been modeled In summary, our results show that these materials in terms of molecular magnetism involving neighboring can be thought of as “high temperature” versions of the tetrahedra of these S = 1/2 moments.56 In GeCo O eff 2 4 S = 1/2 AFM XY pyrochlore, which has been studied the CoO octahedra are not distorted, and the moment eff 6 in relation to the rare earth oxide series of pyrochlores. is therefore isotropic. Their strong interactions allow the exploration of a wider An interesting question for future exploration is temperature range of this model. The inherent local dis- whether there exists a Co2+-based spin ice material (ei- order arising from the split non-magnetic A-site is likely ther quantum or classical). Such a material would have to produce exchange disorder, in part through the local the advantage that the emergent monopole dynamics variations of the g-tensor that we have estimated here. would occur at temperatures much easier to access, pos- The role of this exchange disorder in a model appropri- sibly even high enough to be useful for applications. If ate to these materials requires more investigation, as the the U(1) quantum spin liquid phase of quantum spin ice ground states appear to be contrary to the prediction could be produced in a Co2+ pyrochlore, its emergent that quenched order by disorder selects long range order. photon modes would have a more accessible bandwidth Finally, NaCaCo2F7 and NaSrCo2F7 are more frustrated compared to the rare earth-based candidates, making (f ∼ 50) versions of the AFM XY pyrochlore than the them more easily identifiable with low energy inelastic canonical example, Er2Ti2O7 (f = 20). This may in- neutron scattering and thermodynamic measurements. dicate that they are closer to a phase boundary in the To achieve the necessary Ising anisotropy within a Co2+ general anisotropic exchange model developed for the py- pyrochlore, the octahedral environment would need to rochlore lattice; this proximity would make the materials be elongated along the h111i directions rather than com- more sensitive to quantum fluctuations, and may even pressed as it is for NCCF and NSCF. put them in proximity to a quantum disordered phase. 11

ACKNOWLEDGMENTS Appendix A: Q-dependence of modes observed with Ei = 700 meV and 250 meV

The inelastic modes shown in Figure 5 a)-d) can be identified as being magnetic in origin due to the |Q|- dependence of their intensities. The intensity of mag- The authors gratefully acknowledges discussions with netic features (such as single-ion levels) should decrease K. Kimura, the use of software adapted from work by with increasing |Q| due to the magnetic form factor.59 G. Hester and J. O’Brien, and the encouragement of C. Meanwhile the intensity of features arising from phonons L. Broholm. The neutron scattering data were reduced (whether single- or multi-phonon processes) should show using the Mantid57 and DAVE58 software packages. The the opposite trend. For the mode identified near 360 crystal growth activities at Princeton were supported by meV in NCCF (using Ei = 700 meV), the |Q-dependence the US Department of Energy, office of Basic Energy of intensity reveals a combination of both magnetic and Sciences, Division of Material Sciences and Engineering multi-phonon contributions. This is demonstrated in Fig- under grant DE-FG02-08ER46544. This work utilized ure 7 a), which compares intensity near 360 meV to that neutron scattering facilities supported in part by the Na- in energy ranges just above and just below the observed tional Science Foundation under Agreement No. DMR- mode. This demonstrates that the peak observed near 1508249. Research using ORNL’s Spallation Neutron 360 meV is magnetic in origin. The intensity of the fea- Source was sponsored by the Scientific User Facilities Di- tures observed using Ei = 250 meV are also shown for vision, Office of Basic Energy Sciences, U.S. Department three different Q ranges in panel b), demonstrating their of Energy. magnetic nature.

∗ author to whom correspondence should be addressed: NCCF Ei = 700 meV NCCF Ei = 250 meV a) b) [email protected] 2 50 6 1 K. A. Ross, J. W. Krizan, J. A. Rodriguez-Rivera, R. J. 40 4 Cava, and C. L. Broholm, Physical Review B 93, 014433 2 1.5 30 (2016). 0 2 20 I (arb. units) J. S. Gardner, M. J. P. Gingras, and J. E. Greedan, Re- 100 150 views of Modern Physics 82, 53 (2010). 1 10 E (meV) I (arb. units) 3 300 - 400 meV I (arb. units) -1 R. D. Shannon and A. W. Sleight, Inorganic Chemistry 7, 0 2-4 Å 250 - 300 meV -1 4-6 Å 1649 (1968). 0.5 400 - 500 meV -10 6-8 Å-1 4 Z. L. Dun, X. Li, R. S. Freitas, E. Arrighi, C. R. D. Cruz, 5 10 15 0 50 100 150 200 M. Lee, E. S. Choi, H. B. Cao, H. J. Silverstein, C. R. -1 Q (A ) E (meV) Wiebe, et al., Physical Review B 92, 140407 (2015). 5 Y. Q. Cai, Q. Cui, X. Li, Z. Dun, J. Ma, C. dela Cruz, Y. Jiao, J. Liao, P. Sun, Y. Q. Li, et al., Physical Review FIG. 7. a) Q dependence of intensity for energy ranges span- B 93, 014443 (2016). ning the identified single-ion mode at ∼ 360 meV (300 to 400 6 C. Castelnovo, R. Moessner, and S. L. Sondhi, Annual meV), just above the single-ion mode (400 to 500 meV, inten- Review of Condensed Matter Physics 3, 35 (2012). sity shifted up for comparison), and just below the single-ion 7 S. T. Bramwell, M. J. P. Gingras, and J. N. Reimers, mode (250-300 meV, intensity shifted down for comparison). Journal of Applied Physics 75, 5523 (1994). b) Energy dependence of neutron scattering intensity in dif- 8 A. W. C. Wong, Z. Hao, and M. J. P. Gingras, Physical ferent Q ranges for Ei =250 meV. Review B 88, 144402 (2013). 9 P. A. McClarty, P. Stasiak, and M. J. P. Gingras, Physical Review B 89, 024425 (2014). 10 J. D. M. Champion, M. J. Harris, P. C. W. Holdsworth, 14 V. S. Maryasin and M. E. Zhitomirsky, Physical Review B A. S. Wills, G. Balakrishnan, S. T. Bramwell, E. Ciˇzm´ar,ˇ 90, 094412 (2014). T. Fennell, J. S. Gardner, J. Lago, et al., Physical Review 15 A. Andreanov and P. A. McClarty, Physical Review B 91, B 68, 020401 (2003). 064401 (2015). 11 M. E. Zhitomirsky, M. V. Gvozdikova, P. C. W. 16 K. A. Ross, T. Proffen, H. A. Dabkowska, J. A. Quilliam, Holdsworth, and R. Moessner, Physical Review Letters L. R. Yaraskavitch, J. B. Kycia, and B. D. Gaulin, Phys- 109, 077204 (2012). ical Review B 86, 174424 (2012). 12 L. Savary, K. A. Ross, B. D. Gaulin, J. P. Ruff, and L. Ba- 17 A. Yaouanc, P. D. de R´eotier, C. Marin, and V. Glazkov, lents, Physical Review Letters 109, 167201 (2012). Physical Review B 84, 172408 (2011). 13 J. Oitmaa, R. R. P. Singh, B. Javanparast, A. G. R. Day, 18 T. Taniguchi, H. Kadowaki, H. Takatsu, B. F˚ak,J. Ol- B. V. Bagheri, and M. J. P. Gingras, Physical Review B livier, T. Yamazaki, T. J. Sato, H. Yoshizawa, Y. Shimura, 88, 220404 (2013). T. Sakakibara, et al., Physical Review B 87, 060408 (2013). 19 J. W. Krizan and R. J. Cava, Physical Review B 89, 214401 (2014). 12

20 J. W. Krizan and R. J. Cava, Journal of Physics: Con- 42 J. Gaudet, D. D. Maharaj, G. Sala, E. Kermarrec, K. A. densed Matter 27, 296002 (2015). Ross, H. A. Dabkowska, A. I. Kolesnikov, G. E. Granroth, 21 J. W. Krizan and R. J. Cava, Physical Review B 92, 014406 and B. D. Gaulin, Physical Review B 92, 134420 (2015). (2015). 43 P. Babkevich, A. Finco, M. Jeong, B. Dalla Piazza, I. Ko- 22 M. B. Sanders, J. W. Krizan, K. W. Plumb, T. M. Mc- vacevic, G. Klughertz, K. Kr¨amer, C. Kraemer, D. T. Queen, and R. J. Cava, Journal of Physics: Condensed Adroja, E. Goremychkin, et al., Physical Review B 92, Matter 29, 045801 (2016). 144422 (2015). 23 M. P. Zinkin, M. J. Harris, and T. Zeiske, Physical Review 44 K. W. H. Stevens, Proceedings of the Physical Society. B 56, 11786 (1997). Section A 65, 209 (1952). 24 H. Bl¨ote,R. Wielinga, and W. Huiskamp, Physica 43, 549 45 M. T. Hutchings, Solid state physics 16, 227 (1964). (1969). 46 M. Rotter, D. Manh Le, J. Keller, L. G. Pascut, T. Hoff- 25 J. A. Hodges, P. D. de R´eotier,A. Yaouanc, P. C. M. mann, M. Doerr, R. Schedler, P. Fabi, S. Rotter, M. Banks, Gubbens, P. J. C. King, and C. Baines, Journal of Physics: and N. Kluver, “Mcphase users manual: Symmetry consid- Condensed Matter 23, 164217 (2011). erations for crystal field parameters,” . 26 W. J. L. Buyers, T. M. Holden, E. C. Svensson, R. A. 47 J. Kanamori, Progress of Theoretical Physics 17, 177 Cowley, and M. T. Hutchings, Journal of Physics C: Solid (1957). State Physics 4, 2139 (1971). 48 G. L. Squires, Introduction to the theory of thermal neutron 27 A. Abragam and B. Bleaney, Electron paramagnetic reso- scattering (Cambridge university press, 2012). nance of transition ions (OUP Oxford, 2012). 49 J. Jensen and A. R. Mackintosh, Rare earth magnetism 28 M. E. Lines, Physical Review 131, 546 (1963). (Clarendon Oxford, 1991). 29 J. P. Goff, D. A. Tennant, and S. E. Nagler, Physical 50 E. A. Oliveira, I. Guedes, A. P. Ayala, J. Y. Gesland, J. El- Review B 52, 15992 (1995). lena, R. L. Moreira, and M. Grimsditch, Journal of Solid 30 R. M. White, Quantum theory of magnetism, 3rd ed., State Chemistry 177, 2943 (2004). edited by Springer (Springer, 2006). 51 K. Kimura, S. Nakatsuji, and T. Kimura, Physical Review 31 K. A. Ross, L. Savary, B. D. Gaulin, and L. Balents, B 90, 060414 (2014). Physical Review X 1, 021002 (2011). 52 K. Baroudi, B. D. Gaulin, S. H. Lapidus, J. Gaudet, and 32 S. Guitteny, S. Petit, E. Lhotel, J. Robert, P. Bonville, R. J. Cava, Physical Review B 92, 024110 (2015). A. Forget, and I. Mirebeau, Physical Review B 88, 134408 53 The E ∼ 360 meV mode was measured only in NCCF due (2013). to time limitations. Given the overall similarity if the other 33 H. Yan, O. Benton, L. D. C. Jaubert, and N. Shannon, single-ion levels in the two materials, this mode is assumed arXiv preprint arXiv:1603.09466 (2016). to lie in same energy range in NSCF. 34 J. D. M. Champion and P. C. W. Holdsworth, Journal of 54 B. J. Campbell, H. T. Stokes, D. E. Tanner, and D. M. Physics: Condensed Matter 16, S665 (2004). Hatch, Journal of Applied Crystallography 39, 607 (2006). 35 A. Poole, A. S. Wills, and E. Lelievre-Berna, Journal of 55 R. D. Shannon, Acta Crystallographica Section A: Crystal Physics: Condensed Matter 19, 452201 (2007). Physics, Diffraction, Theoretical and General Crystallog- 36 S. Petit, J. Robert, S. Guitteny, P. Bonville, C. Decorse, raphy 32, 751 (1976). J. Ollivier, H. Mutka, M. J. Gingras, and I. Mirebeau, 56 K. Tomiyasu, M. K. Crawford, D. T. Adroja, P. Manuel, Physical Review B 90, 060410 (2014). A. Tominaga, S. Hara, H. Sato, T. Watanabe, S. I. Ikeda, 37 J. Gaudet, A. Hallas, D. Maharaj, C. Buhariwalla, E. Ker- J. W. Lynn, et al., Physical Review B 84, 054405 (2011). marrec, N. Butch, T. Munsie, H. Dabkowska, G. Luke, 57 O. Arnold, J.-C. Bilheux, J. Borreguero, A. Buts, S. I. and B. Gaulin, Physical Review B 94, 060407 (2016). Campbell, L. Chapon, M. Doucet, N. Draper, R. F. Leal, 38 R. Sarkar, J. W. Krizan, F. Br¨uckner, E. C. Andrade, M. A. Gigg, et al., Nuclear Instruments and Methods in S. Rachel, M. Vojta, R. J. Cava, and H.-H. Klauss, arXiv Physics Research Section A: Accelerators, Spectrometers, preprint arXiv:1604.00814 (2016). Detectors and Associated Equipment 764, 156 (2014). 39 J. P. C. Ruff, J. P. Clancy, A. Bourque, M. A. White, 58 R. T. Azuah, L. R. Kneller, Y. Qiu, P. L. Tregenna- M. Ramazanoglu, J. S. Gardner, Y. Qiu, J. R. D. Copley, Piggott, C. M. Brown, J. R. Copley, and R. M. Dimeo, M. B. Johnson, H. A. Dabkowska, et al., Physical review Journal of Research of the National Institute of Standards letters 101, 147205 (2008). and Technology 114, 341 (2009). 40 J. A. Rodriguez, D. M. Adler, P. C. Brand, C. Broholm, 59 S. W. Lovesey, Theory of neutron scattering from con- J. C. Cook, C. Brocker, R. Hammond, Z. Huang, P. Hun- densed matter, Vol. 2 (United Kingdom, 1984). dertmark, J. W. Lynn, et al., Measurement Science and Technology 19, 034023 (2008). 41 G. E. Granroth, D. H. Vandergriff, and S. E. Nagler, Phys- ica B: Condensed Matter 385, 1104 (2006).