Arxiv:2004.07325V1 [Cond-Mat.Str-El] 15 Apr 2020 3+ Ybalo3 [1, 2]

Experimental observation of magnetic dimers in diluted Yb:YAlO3

S. E. Nikitin,1, 2 Tao Xie,3 A. Podlesnyak,3, ∗ and I. A. Zaliznyak4, † 1Max Planck Institute for Chemical Physics of Solids, D-01187 Dresden, Germany 2Institut fürFestkörper- und Materialphysik, Technische UniversitätDresden, D-01069 Dresden, Germany 3Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 4Condensed Matter Physics and Materials Science Division, Brookhaven National Laboratory, Upton, NY 11973, USA

We present a comprehensive experimental investigation of Yb magnetic dimers in Yb0.04Y0.96AlO3, an Yb-doped Yttrium Aluminum Perovskite (YAP) YAlO3 by means of specific heat, magnetization and high-resolution inelastic neutron scattering (INS) measurements. In our sample, the Yb ions are randomly distributed over the lattice and ∼ 7% of Yb ions form quantum dimers due to nearest- neighbor antiferromagnetic coupling along the c-axis. At zero field, the dimer formation manifests itself in an appearance of an inelastic peak at ∆ ≈ 0.2 meV in the INS spectrum and a Schottky-like anomaly in the specific heat. The structure factor of the INS peak exhibits a cosine modulation along the L direction, in agreement with the c-axis nearest-neighbor intra-dimer coupling. A careful fitting of the low-temperature specific heat shows that the excited state is a degenerate triplet, which indicates a surprisingly small anisotropy of the effective Yb-Yb exchange interaction despite the low crystal symmetry and anisotropic magnetic dipole contribution, in agreement with previous reports for the Yb parent compound, YbAlO3 [1, 2], and in contrast to Yb2Pt2Pb [3, 4]. The obtained results are precisely reproduced by analytical calculations for the Yb dimers.

I. INTRODUCTION information storage schemes with coherent and entangled photons [7–11]. Rare-earth based insulating crystals such Rare Earth based systems have recently been widely as Yb:YAP (Yb-doped Yttrium Aluminum perovskite, considered among the leading candidates for quantum YbxY1−xAlO3), Yb:YAG (Yb-doped Yttrium Aluminum information storage and processing [5–9]. In rare earths, garnet, YbxY3−xAl5O12), or Nd:Y2SiO5 are well known strong relativistic spin-orbit coupling ties the spin and laser materials with high quantum efficiencies, long re- the orbital angular momentum together, so they form a laxation times and narrow optical transition lines in the multiplet of total angular momentum, J , often with large near infrared range overlapping with telecom wavelengths J . In the presence of crystal electric field (CEF) the [12, 13]. Quantum repeaters [6], teleportation between (2J + 1)-fold degeneracy of the orbital multiplet is lifted distant condensed matter qubits [10], optical preparation and the states of the rare-earth ion are split in energy. In of coherent dark states [14], quantum photon teleporta- the case of Kramers ions, such as Ce3+, Nd3+, Dy3+, Er3+, tion and storage in a quantum memory [15], and quantum Yb3+, the CEF level splitting retains Kramers twofold storage of entangled photons [16] and an all-optical re- degeneracy and the lowest doublet can be well separated trieval [17], were all recently demonstrated in rare-earth in energy from the rest of the multiplet, forming quantum laser materials. degree of freedom that can be described by an effective Recently, magnetic quantum doublet degree of free- spin-1/2. These quantum doublets are protected by time- dom in materials with 4f ions with odd number of reversal symmetry and can form entangled spin states. f-electrons, such as Ce3+, Nd3+, Yb3+, etc., which Their peculiar coupling with external electromagnetic filed is described by an effective spin-1/2, came to promi- is governed by the angular momentum, J , conservation, nence in the context of novel quantum magnets. thus offering a route towards protected quantum states in Among these are two-dimensional triangular-lattice sys- rare earth ion magnets, which are suitable for solid state tems [18, 19], three-dimensional pyrochlore-lattice spin- quantum information processing and storage. liquids, e.g. Ce2Zr2O7 [20], quantum dimer mag- In the case of 4f ions with non-zero orbital angular net Yb2Si2O7 [21], one-dimensional spin-chain antifer- momentum, the CEF splitting is typically only a modest romagnets (AFM) Yb2Pt2Pb [3, 4], YbFeO3 [22] and perturbation to the dominant spin-orbit splitting. Con- arXiv:2004.07325v1 [cond-mat.str-el] 15 Apr 2020 3+ YbAlO3 [1, 2]. Traditionally, the physical realization of sequently, in an ion such as Yb , there are a number such systems were found among the 3d-electron (mostly of J → J ± 1 optical transitions at slightly different Cu2+) materials, such as frustrated spin-liquid can- wavelengths, which can be individually addressed by pho- didate herbertsmithites (see Ref. [23] and references tons. The compatibility with optical technology makes therein), two-dimensional Heisenberg AFM PHCC [24] Yb and sister rare earths (Ce, Nd, Er, Dy) attractive for and Cu(DCOO)2·4D2O[25], one-dimensional spin-chain incorporating in optics-based quantum computation and materials SrCuO2 and Sr2CuO3 [26–29], etc. These materials, however, usually have significant disadvantages, which complicate their investigations: (i) a small mag- ∗ Correspondence to: [email protected] netic moment; (ii) a very strong exchange interaction, † [email protected] typically J > 10 meV, which means that one has to ap- 2

(a) (c) ply an exceptionally high magnetic field to influence the Monomer (Yb1 ) |S=-1/2〉 ground state. In rare earth based systems, on the other z +g½ m B|J =-7/2〉 hand, the overlap of the 4f electronic wavefunctions is Y/Yb |S=±1/2〉 B |Jz =±7/2 small and effective spin interactions are often quite weak, 〉 -½ m gBB |S=+1/2〉 c Al z making detailed magnetic field studies possible [1–4]. |J =+7/2〉 YbAlO3, an end member of Yb:YAP family of ma- a b O Dimer (Yb ) terials was shown to be a good realization of S = 1/2 2 |Sz =-1 one-dimensional antiferromagnet, which exhibits spinon (b) (d) tot 〉 |Jz =-7 100 tot 〉 confinement-deconfinement transition and quantum crit- |S =1,Sz =0 m 10-1 tot tot 〉 +gB B icality in different regions of its phase diagram [1]. Its |Jz =0 10-2 tot 〉 physics can be reasonably well described by a combination -3 z 10 |Stot =1,S tot =±1〉 of nearest-neighbor intrachain exchange interaction along 10-4 |Jz =±7 -gm B

% tot 〉 B D D0 z Probability -5 ± the c-axis and weak in-plane interchain dipolar interaction, 10 |S tot =+1〉 -6 z

0.02

92.16 %

7.37 %

0.44 % 10 |Stot =0〉 |J =+7 which stabilizes the magnetic order below TN = 0.88 K. It Monomer DimerTrimerTetramer tot 〉 Yb Yb Yb Yb |Jz =0 is worth noting that this is an unexpected result because 1 2 3 4 tot 〉 the crystal structure of YbAlO3 does not naturally imply any “special” chain direction (this is similar to what was FIG. 1. Crystal structure and energy levels of Yb dimers found in Yb2Pt2Pb [3]). Moreover, while the symmetry of in Yb:YAlO3. (a) Orthorhombic perovskite P bnm crystal structure of Y Yb AlO where Yb form chain segments Yb site is very low (Cs) and the CEF environment induces 1−x x 3 an exceptionally high magnetic anisotropy of the ground along the c-axis. (b) Calculated probability for Yb ions to form n-mers of different length; the ratio n·P (Ybn+1) is constant state doublet (gz >> gxy), the analysis of INS spectra and (n+1)·P (Ybn) critical exponents has shown that the intrachain exchange and equal to the concentration of Yb, x ≈ 0.043. (c) Effective interaction is very isotropic, 0.88J < J < J . spin-1/2 energy levels of the ground state doublet of an isolated xy ∼ z ∼ xy 3+ We note that low-dimensional S = 1/2 magnets re- Yb ion and (d) splitting of the Yb dimer levels in the main at the cutting edge of the solid state physics for presence of XXZ interaction, Eq. (1), and magnetic field, B. the last decades. Depending on dimension of the lattice they can host exotic ground states and excitations. For down to 1.8 K. Specific heat measurements of m = 11.5 mg sample of Yb:YAlO3 were carried out using a commercial instance, two-dimensional triangular AFM is predicted to 3 host entangled quantum spin-liquid ground state. One- PPMS from Quantum Design with He option. Magnetic dimensional spin models show no magnetic ordering down field was applied along the easy a-axis. to zero temperature and exhibit fractionalized spinon Neutron scattering measurements were performed at excitations. Even the simplest toy model – a dimer con- the time-of-flight Cold Neutron Chopper Spectrometer sisting of two spins coupled by an AFM interaction has (CNCS) [31, 32] at the Spallation Neutron Source at an entangled ground state, ψ = √1 |↑↓ − ↓↑i, which Oak Ridge National Laboratory. Data were collected 2 with a single crystal Yb:YAlO sample with mass around has no clear analogy in classical physics and presents a 3 0.6 g, which was aligned in the (0KL) scattering plane. useful two-qubit model for quantum information science Triton - Cryofree dilution refrigeration system from Ox- (QIS). Hence, studies of these systems presents significant ford Instruments was used to cool the sample down to interest. T = 0.1 K. The data were collected using fixed incident In this work, we investigate the magnetic Hamiltonian neutron energies of E = 1.55 meV (λ = 7.27 A)˚ and of Yb-Yb effective spin interactions in Yb:YAlO by mea- i i 3 E = 3.32 meV (λ = 4.96 A)˚ resulting in a full-width-at- suring the diluted sample, Yb Y AlO (x ≈ 0.04). In i i x 1−x 3 half-maximum energy resolution at the elastic position this material, the Yb ions are statistically distributed of 0.04 meV and 0.11 meV, respectively. The measure- over the Y sites. If two Yb ions would occupy nearest ments were performed using the rotating single crystal neighbor positions along the c-axis, they can form a mag- method. All time-of-flight data-sets were combined to netic dimer. By studying the dimer physics one can refine produce a four-dimensional scattering-intensity function the parameters of Yb-Yb exchange interactions with very I(Q, ¯hω), where Q is the momentum transfer and ¯hω is high precision. Thereby, here we aim to answer two main the energy transfer. The software packages [33] questions: (i) Are there magnetic dimers in Yb:YAlO ? Horace 3 and [34] were used for data reduction and (ii) What is the exchange Hamiltonian for the dimers? MantidPlot analysis.

II. EXPERIMENTAL DETAILS III. RESULTS AND ANALYSIS

The YAlO3 single crystal doped with (nominal) 5% 3+ A. Origin of the dimers of Yb (Yb:YAlO3) was procured commercially from Scientific Materials [30]. Magnetic measurements were carried out using MPMS3- The orthorhombically distorted perovskite crystal struc- VSM with magnetic fields up to 7 T and temperatures ture of Yb:YAlO3 (P bnm, a = 5.180 A,˚ b = 5.330 A,˚ 3 c = 7.375 A˚ at room temperature) is illustrated in Fig.1(a) Hamiltonian (1) is straightforwardly diagonalized in terms [1, 2, 35, 36]. In pure YbAlO3, 4f orbitals of Yb ions have of the total effective dimer spin, Stot = S1 + S2, and its z non-negligible overlaps only along the c-axis [1], meaning z-component, Stot, that Yb magnetic moments form chains where each Yb is only coupled with its two nearest neighbors. This is 1 2 ∆ − 1 2 2∆ + 1 H = J (S ) + (Sz ) − . (2) due to the hexagonal, 6-fold symmetry of L = 3 orbital 2 tot 2 tot 4 wave function of a single unpaired electron in 4f 13 shell 3+ of Yb and is similar to the case of Yb2Pt2Pb [3]. In the presence of anisotropy, ∆ 6= 1, the Stot = 1 triplet z z As we show below, in our Yb:YAlO3 sample ≈ 96 % Yb splits into a Stot = ±1 doublet and a Stot = 0 singlet. 1+∆ ions are substituted by nonmagnetic Y ions (the nominal For ∆ < 1, the doublet gap, ∆± = J 2 , is smaller than concentration is 95 %). Assuming that the Yb ions are the gap to a singlet, ∆0 = J (Fig.1). randomly distributed over the sample, we can calculate In magnetic field, B, directed along the quantization z- z the probability to find Yb multimer of length n both ana- axis, the Stot = ±1 doublet splits according to [Fig.1 (d)], lytically and numerically, using Monte-Carlo simulation. 1 + ∆ If we consider a single Yb ion, the probability that it will ∆±(B) = ∆± ∓ gµBB = J ∓ gµBB have only Y neighbors along the chain direction is equal 2 2 ∆0(B) = ∆0 = J. (3) to P1 = (1−x) (throughout the text x represents relative concentration of Yb ions in the sample, YbxY1−xAlO3). 2 Here, Lande g-factor relates magnetic moment along the It is easy to show that for dimer P2 = x(1−x) , for trimer 2 2 n−1 2 field direction with the corresponding effective spin-1/2 P3 = x (1 − x) and finally for n-mer Pn = x (1 − x) . z z P∞ component, M = gµBS . Above the critical field, Bc = The total number of all n-mers in the chain, Pi = i ∆±/(gµB), the gap to ∆+ level closes and it becomes the 1 − x, counts the number of Y ions and the probability of ground state. For B > B , Yb ion to be in an n-mer is nxn−1(1 − x)2. Thereby, in c

Yb:YAlO3 92 % of all Yb ions occupy single sites, 7.4 % ∆−(B) = ∆± + 2gµB(B − Bc) form magnetic dimers and ∼ 0.5 % form longer structures. ∆ (B) = ∆ + gµ (B − B ) (4) To crosscheck these conclusions classical Monte-Carlo 0 0 B c calculations were used. We generated a chain with length ∆0,0(B) = gµB(B − Bc), 8 6 of 10 ions with 4 · 10 randomly distributed impurity where ∆ is the gap between the new ∆ ground state atoms and calculated number of single Yb, dimers, etc. 0,0 + and the Stot = 0 singlet. Figure1 (b) shows the resulting histogram, which is in While all three excited states contribute to entropy and perfect agreement with the analytical results. the heat capacity both below and above Bc, the dipole selection rules, ∆J z = 0, ±1 (angular momentum conservation), only allow magnetic transitions between the B. Effective spin Hamiltonian and level splitting z singlet, |Stot = 0i, and the |Stot = 1,Stot = 0i component of the triplet. Hence, only this, ∆J z = 0 transition con- The crystal electric field (CEF) splits the J = 7/2 Yb tributes to neutron scattering cross-section. The energy spin-orbit octet in YbAlO3 into four Kramers doublets. of this transition is ∆0 = J, independent of the anisotropy The lowest doublet corresponds to nearly pure J z = ±7/2 and magnetic field, but its intensity fades above Bc where and is separated by a substantial energy gap, ∆E/kB ∼ the transition becomes thermally activated. At zero mag- 300 K (kB is the Boltzmann constant), from the first netic field, the intensity of this singlet-triplet transition excited doublet [1, 3]. This leads to an effective spin- is determined by the level populations [38, 39] and also 1/2 description of the low-energy and low-temperature decreases with the increasing temperature, according to behavior for T ∆E/kB where the higher doublets are the temperature prefactor, not populated. The effective Hamiltonian is that of an S = 1/2 XXZ chain [3], which for a dimer reduces to, 1 p(T ) = , (5) ∆0 y y − x x z z 1 + 3e kB T H = J (S1 S2 + S1 S2 + ∆S1 S2 ) . (1) Here, J is the effective exchange coupling and ∆ is its which changes from 1 at T = 0 to ' 0.25 in the equipar- anisotropy. In Yb2Pt2Pb, it was found that the effective tition regime at high temperature. Hamiltonian has Ising anisotropy, ∆ ∼ 2−3 and J ∼ 0.1− 0.2 meV [3]. In the case of YbAlO3, Ising gap is absent in C. Specific heat the magnetically disordered state above TN and isotropic XXX Hamiltonian with J = 0.2 meV adequately describes < magnetic excitations, albeit an XY-type anisotropy, ∆ ∼ We start exposition of our experimental results with 1, could not be definitively ruled out [1,2]. the low-temperature magnetic specific heat measured at Measuring the excitation spectra of a diluted system different magnetic fields, which is presented in Fig.2. dominated by dimers provides a simple and direct way to The phonon contribution was fitted with the standard probe the Hamiltonian, such as Eq. (1) [37]. The dimer Debye power-law dependence, C ∝ T 3, and subtracted 4

0.02 The analysis revealed ∆± = 0.2146(9) meV, ∆0 = 0.2 (a) B = 0.2 T (b) 0.2142(9) meV, and concentration x = 0.0439(3). Surpris- B = 0.5 T 0.015 0.18 B = 1.0 T ingly, within rather high precision of our ﬁt we did not B = 3.0 T 0.16 B = 6.0 T 0.01 observe any appreciable splitting of the triplet. Hence,

0.14 0.005 we find ∆± ≈ ∆0, implying effective isotropic Heisenberg Specific heat (J/moleK) 0.12 exchange interaction, ∆ ≈ 1 in Eq. (1), in agreement with 0 100 0.1 what was inferred in Ref. 1 from the analysis of magnetic Temperature (K) excitations and free fermion behavior near the saturation 0.08 (c) field in YbAlO3. This is surprising because it requires 0.06 100 Specific heat (J/moleK) fine tuning of the anisotropic superexchange [3] and mag- 0.04 netic dipole interaction between the dimer moments. In

0.02 Energy (meV) YbAlO3 the effective g-factor, which relates magnetic 10-1 α α 0 moment to the effective spin-1/2, M = µBgαS , is very 0 0 10 10 Magnetic Field (T) anisotropic, gz gx,y ' 0 [1, 2]. Hence, magnetic dipole Temperature (K) z z interaction, ∼ M 1M 2, provides anisotropic, ∼ S1 S2 contribution to the effective spin coupling. FIG. 2. Low-temperature specific heat of Yb:YAlO3. (a) Tem- perature dependences of specific heat measured at different To account for the contribution of the longer structures magnetic field (B k [100]) along with the results of fitting we also included the trimer subsystem, and now the total using Eq. (9). (b) Specific heat taken at 0 T. Note that the specific heat is described as: y-scale is an order of magnitude smaller compare to panel 2 2 (a). Blue and red lines show calculated specific heat of Yb CMag(T ) = x (1 − x) (CDimer(T ) + xCTrimer(T )). (8) subsystem taking into account only the dimer and the dimer The energy splitting for the trimer subsystem at zero plus trimer contributions using Eq. (7) and Eq. (8), respec- field was determined by solving Schrodinger equation for tively. (c) Energy splitting for Yb1 subsystem as a function of magnetic field extracted from the specific heat measurements the linear spin trimer with isotropic exchange interaction: and its linear fitting E(B) = gµBBS. J(S1S2 + S2S3)ψ = Eψ [the trimer Hamiltonian can 1 2 2 be diagonalized as 2 J (S1 + S2 + S3) − (S1 + S3) ]. from the data. The curve taken at 0 T is shown in Fig.2 The splitting scheme is as follows: a doublet ground state, 3 (b) and one can see a clear Schottky-like anomaly with a doublet at E1 = J, and a quartet at E2 = 2 J (note that a maximum at T ≈ 0.85 K. Figure2 (a) presents the due to the half-integer total spin of the trimer all states data measured at different magnetic fields up to 6 T. have Kramers degeneracy). The fitting of specific heat One can see that: (i) a new, strong, Schottky-like peak with Eq. (8) revealed slightly lower Yb concentration, x = appears, which arises from the contribution of individual 0.0430(2), whereas the energy splitting remains essentially Yb moments, Yb1, to the specific heat and significantly unchanged, ∆± = 0.212(5) meV, ∆0 = 0.211(5) meV. exceeds the zero-field signal in intensity; (ii) magnetic Specific heat of Yb:YAlO3 at B = 0 calculated for field gradually shifts the peak due to the Zeeman effect. x = 0.043 is shown in Fig.2 (b). The blue and red To describe the measured specific heat of Yb:YAlO3 lines show only the dimer and the dimer plus trimer we use a standard equation for specific heat of a discrete contributions, respectively. Good agreement between the p-level system, observed and calculated curves proves the predominantly dimer origin of the observed signal and gives an estimate p=2n of the contribution from the longer structures. d 1 Ei X − k T Cn−mer(T ) = R Eie B , (6) Having determined the real concentration of Yb in the dT Z i=1 sample, let us turn to the analysis of the specific heat measured in magnetic field. Application of magnetic field where R is universal gas constant, E are the energy levels, i modifies the dimer spectrum according to Eqs. (3),(5) Z is the partition function and summation runs over all and splits the degenerate doublet state of Yb monomers states of the system (p = 2n for an n-mer subsystem). according to E(B) = g µ B [Fig.1 (c,d)]. Thereby, we At low magnetic field, B < B , the splitting of the Yb a B c fit the measured specific heat as a sum of Yb monomer doublet is described by Eqs. (3). Note that at zero field and dimer contributions [40]: Yb monomers have a degenerate doublet ground state and 2 do not contribute to the specific heat. Thus, the signal CMag(T ) = x(1 − x) (CMono(T ) + xCDimer(T )). (9) can be considered as a specific heat of dimers plus a weak contribution of longer n-mers with n > 2. At B = 0 the In the fit, we fixed x = 0.043, while the doublet splitting, dimer level splitting scheme is a doublet at ∆± and a Ei, which determines the g factor for Yb1 subsystem were singlet at ∆0 [see Fig.1(d)]. In order to evaluate these varied as free parameters. The experimental data along parameters, as well as the actual concentration of Yb in with the fitted curves are shown in Fig.2 (a) by symbols our sample we used Eq. (6) to fit the B = 0 data, taking and solid lines, respectively, and one can see that they are into account the dimer concentration [Fig.2 (b)], in good agreement. The extracted g-factor that describes the Zeeman splitting of Yb1 [Fig.2 (c)] is ga = 6.21(3), 2 2 CMag(T ) = x (1 − x) CDimer(T ). (7) somewhat below ga = 6.9 observed in pure YbAlO3 [1]. 5

4 9 0 ( a ) T = 2 K measured on diﬀerent small single crystals with m ∼ 8 mg; 3 B = 7 T 4 5 ) a slight discrepancy of M and g obtained from these

b S a 3 ) ( b ) Y b

/ 2 B Y /

µ measurements reﬂects an uncertainty of sample weight. B (

µ (

n T = 2 K 1

o n i

t 2 B | | [ 1 0 0 ] o i a t 0 1 8 0 0 z B | | [ 0 0 1 ] i a t z i e t n

e 1 g n E. Neutron scattering a

1 g a M 2 M

3 2 2 5 3 1 5 0 0 2 4 6 Inelastic neutron scattering (INS) is a powerful tool for 2 7 0 M a g n e t i c F i e l d ( T ) investigation of magnetic excitations (see Ref. [42–45] and FIG. 3. (a) Field dependence of magnetization measured at references therein) and has proven extremely successful T = 2 K along [100] and [001] directions. Green solid line in characterizing magnetic clusters [37, 46]. We write the represents the ﬁt of Ma(B) with Brillouin function. (b) An- INS cross-section in the following form [42–45], gular dependence of magnetization in ab-plane measured at ◦ ◦ 2 B = 7 T, T = 2 K. 0 and 90 corresponds to B k b and B k a, d σ g Q Q ∝|F (Q)|2 P gα β δ − α β P eiQ·rij respectively. Red line shows the ﬁt with Eq. (10). dΩω α,β 2 2 αβ Q2 i,j X ˆα 0 0 ˆβ D. Magnetization × pλhλ|Si |λ ihλ |Sj |λiδ(¯hω + Eλ − Eλ0 ), (11) λ,λ0

The local environment of Yb ions in Yb:YAlO3 is very corresponding to unpolarized neutrons used in our mea- 3+ similar to the parent compound YbAlO3, thus one can surements. Here, F (Q) is magnetic form factor of Yb ; expect that the crystalline electric field produces a similar α, β = x, y, z and gα, gβ are the corresponding g-factors; type of magnetic anisotropy in both compounds. To |λi and |λ0i are the eigenvectors of the initial and final verify this, we measured magnetization of Yb:YAlO3 at state of the system; pλ is the probability to find the T = 2 K using MPMS3-VSM with horizontal rotator and ˆα ˆβ system in state λ [cf Eq. (5)]; Si and Sj are the effec- the results are summarized in Fig.3. tive spin-1/2 angular momentum operators at sites i and In orthorhombic perovskite structure, Yb moments can j, respectively. The term P = P δ − QαQβ is either point along the c-axis, or lie in ab-plane. As the first α,β α,β Q2 a so-called polarization factor of neutron scattering. It step, we measured the field dependence of magnetization takes into account the property of magnetic dipole inter- with B k [001] and B k [100] and the results are shown action that only magnetic moment components that are in Fig.3 (a). One can see that the signal for B along perpendicular to Q contribute to magnetic neutron scat- [001] direction is very small (∼ 0.2 µ at 7 T) compared B tering cross-section. However, the polarization factor is to the B k [100] case, indicating that the moments are irrelevant for the present study because our spectra were indeed confined in the ab-plane in a similar fashion to that collected within (0KL) scattering plane, whereas the Yb observed in other Yb-based orthorhombic perovskites [1, moments have a preferred orientation along the [100]-axis. 2, 22]. By fitting M (B) and M (B) curves with the a c Therefore, the dominant fluctuations that are favorably Brillouin function for S = 1/2, we found g = 6.23, a weighted by the g-factor in Eq. (11) have P ' 1 for our in good agreement with the heat capacity results, and scattering plane. eiQ·rij gives rise to the structure factor, g = g = 0.6. c ⊥ which causes the modulation of INS intensity in reciprocal To further investigate magnetic anisotropy within the space according to the geometry of the cluster. For Yb ab plane, we measured the angular dependence of magneti- dimers in Yb:YAlO , r ≈ (0, 0, 1 )[47] and the structure zation using the horizontal rotator, as shown in Fig.3 (b). 3 12 2 factor can be reduced to S(Q) = 1 − cos(Q π)[48]. Note that at high magnetic field, above the saturation, L At low temperature, k T J, the spectrum of a Yb moments are confined in ab-plane by CEF and the B dimer measured by INS should consist of a single peak at angular dependence of magnetization can be described z E = ∆0 = J, which corresponds to ∆S = 0 transition. as [1, 41], 1+∆ z The peak at ∆± = J 2 corresponding to ∆S = ±1 is 2 > M relatively suppressed by the g-factor ratio, (gc/ga,b) ∼ M(θ) ≈ S (|cos (θ − ϕ)| + |cos (θ + ϕ)|) , (10) 2 100 [2], and is undetectable. In pure YbAlO3, exchange interaction was found to be J = 0.21 meV [1, 2]. To where MS is the saturation moment of Yb, θ is the angle explore the dynamic response of Yb:YAlO3, we performed between the applied field and the a-axis, and ϕ is the high-resolution INS measurements. The INS spectrum angle between the a-axis and the Yb3+ moment direction. shows a sharp, resolution-limited inelastic peak, which can Result of the fit is shown in Fig.3 by red line. Our be fitted by a Gaussian function with the center energy ◦ analysis revealed MS = 3.57 µB/Yb and ϕ = 26 , which E = 0.2034(4) meV [Fig.4 (a)], very close to the energy ◦ are close to to MS = 3.8 µB/Yb and ϕ = 23.5 observed in of exchange interaction found in pure YbAlO3. pure YbAlO3. In agreement with the above heat capacity To further elucidate the dimer origin of the observed results, the g-factor along the quantization axis is gz = peak, we studied the Q-dependence of the excitation 7.14 (ga = 6.4), somewhat lower than gz = 7.6 found in using the same instrument with orange cryostat (base YbAlO3 [2]. M(B) and M(θ) curves shown in Fig.3 were temperature ∼2 K) and higher neutron incident energies 3 0.5 (a) (b) max 2 0.4 1 .l.u.) 0 0.3 ), (r INS Intensity L (0 0

-1 Energy (meV) min 0.2 -2

-3 0.1 -2 -1 0 1 2 -2 -1 0 1 2 (0 K 0), (r.l.u.) (0 0 L), (r.l.u.)

3 6 (b) (a) Ei = 1.55 meV Ei = 3.32 meV 10 T = 0.1 K T = 2 K (c) 2 8 4 6

1 4 2 Intensity (arb. u.)

INS intensity (arb. u.) 2 INS intensity (arb. u.) 0 0.1 0.2 0.3 -1 0 1 2 3 0 10 20 30 40 Energy (meV) (0 0 L), (r.l.u.) Temperature (K)

FIG. 4. Inelastic neutron scattering spectra of Yb:YAlO3. (a) Energy spectrum measured at T ≈ 0.1 K, Ei = 1.55 meV. The data were integrated within H = 0 ± 0.1, K = 1 ± 0.5, L = 1 ± 0.5. (b) Constant-energy cut along the (00L) direction taken at T ≈ 2 K, Ei = 3.32 meV. The data were integrated in energy window ¯hω = [0.17–0.25] meV; H = 0 ± 0.1, K = 0 ± 1. (c) Intensity of the inelastic peak at Q = (001) as a function of temperature and its fit with Boltzmann function. Error bars represent one standard deviation. in order to cover large area of Q-space. The Q-dependence dimer structure factor, experimental energy resolution, of the observed excitation is presented in Fig.5 (a), which and the observed transition energy E = 0.2 meV. The shows the constant-energy map in (0KL) plane integrated results of the modeling show a nearly perfect agreement within the energy window of the excitation [49]. One can between the experimental and calculated patterns, which see that the spectral intensity is modulated along the L- provides another evidence of the dimer origin of the ob- direction because of the structure factor and shows stripe- served excitation [Figs.5 (a-b)]. like features elongated along the K-direction. Figure5 (b) presents the INS spectrum along (00L) direction where one can see the gapped excitation at E ≈ 0.2 meV whose IV. DISCUSSION AND CONCLUSION intensity is concentrated close to (001). The INS intensity along the (00L) direction integrated within the ¯hω = [0.17, 0.25] meV energy range is shown To summarize, we performed detailed experimental in Fig.4 (b). We fitted the obtained curve using the investigation of Yb-Yb magnetic interactions in diluted equation for the dimer’s structure factor plus a constant Yb:YAlO3, which sheds the light on magnetic Hamiltonian of Yb dimers, n-mers, and pure YbAlO3. For Yb-Yb background, I(QL) = a0 + b(1 − cos(QLπ)), and found a very good agreement between the experimental spectrum dimers, we determined the level splitting and the effective and the fitted curve. The temperature dependence of the spin-1/2 coupling, J, from the heat capacity and high- intensity of the inelastic peak at Q = (001) fits well with resolution neutron scattering measurements. From the magnetic field dependent heat capacity and magnetization, the thermal population factor of the ground state, Eq. (5) 3+ [Fig.4 (c)]. The gap, ∆ = 0 .21(5) meV, obtained from we determined the concentration of the isolated Yb this fit agrees well with that determined directly. moments and their effective g-factor. To corroborate that the measurement of substituted3 6 (a) E = 1.55 meV Ei = 3.32 meV (b) 10 We modeled the experimental spectrum and constant- i T = 0.1 K T = 2 K (c) 3+ Yb:YAlO3 indeed reflects physics of YbAlO3,2 we want to 8 energy map taking into account Yb form factor and 4 stress several key points of our analysis. First of all, the 6 1 4 saturation magnetization per Yb ion and the magnetic 2 Intensity (arb. u.)

INS intensity (arb. u.) 2 INS intensity (arb. u.) 0 3 0.5 anisotropy are very similar for both samples, which0.1 means0.2 0.3 -1 0 1 2 3 0 10 20 30 40 (a) (b) Energy (meV) (0 0 L), (r.l.u.) Temperature (K) max that the CEF and the ground state doublet wavefunc- 2 0.4 tions are also virtually identical. Second, the inelastic 1 peak in the INS spectrum of Yb:YAlO3 is located at .l.u.) 0 0.3 ), (r E = 0.2 meV, very close to the value of exchange inter- INS Intensity L action determined by ﬁtting the spectrum of magnetic (0 0

-1 Energy (meV) min 0.2 excitations in YbAlO3 [1]. The structure factor of the -2 observed excitation exhibits a simple cosine modulation,

-3 0.1 -2 -1 0 1 2 -2 -1 0 1 2 as expected for magnetic dimer. Third, we measured (0 K 0), (r.l.u.) (0 0 L), (r.l.u.) magnetic speciﬁc heat with and without mangetic ﬁeld.

FIG. 5. Inelastic neutron scattering spectra of Yb:YAlO3. In zero magnetic field the dominant contribution is from (a) Observed (top) and modeled (bottom) constant energy Yb-Yb dimers and determines the dimer level splitting. slice in (0KL) scattering plane at T ≈ 2 K, Ei = 3.32 meV. In magnetic field, it is dominated by monomers (single The data are integrated within the energy window ¯hω = Yb moments) and can be quantitatively described using [0.17, 0.25] meV; H = 0 ± 0.1. (b) Observed (right) and a single concentration parameter, x, which corroborates modeled (left) spectra along the (00L) direction at T ≈ 2 K, random statistical distribution of the Yb ions in the sam- Ei = 3.32 meV; K = 0 ± 1; H = 0 ± 0.1. ple. In addition to x, we have also refined the g-factor 7 describing Zeeman splitting of a single Yb3+ doublet in ions, in agreement with what was inferred in the pre- Yb:YAlO3 in magnetic field and obtained the value that vious report on parent compound YbAlO3 [1]. As dis- is within 10% of that determined in YbAlO3. cussed above, this is surprising because a non-negligible We would like to point out that in-field low-temperature magnetic dipole interaction of the nearest Yb moments, µ0 specific heat measurements, such as described in Sec.IIIC, 4π(c/2)3 M 1M 2 (µ0 is the magnetic constant), yields a provide an extremely sensitive quantitative experimental measurable anisotropic contribution to the effective spin method for determining the concentration of magnetic −3 P 2 α α interaction, ∼ (10 meV) · α gαS1 S2 , which appears impurities in a bulk sample. Its sensitivity and capability to be balanced by the anisotropic exchange coupling, sug- of averaging over the large sample volume are well be- gesting serendipitous fine-tuning. Our results provide yond the standard non-destructive analytical techniques, important basis for understanding quantum states of Yb such as energy dispersive x-ray spectroscopy analysis. An n-mers and chains in Yb:YAlO3 system and for their pos- alternative bulk-sensitive thermodynamic probe is the sible future applications in quantum information science. magnetization measurement. Magnetic impurities exhibit paramagnetic behavior, which is well understood and can be described by Brillouin function, allowing to extract the total number of magnetic ions in the sample. How- ACKNOWLEDGMENTS ever, such an analysis requires a-priori knowledge of the magnetic moment of the impurity ion in the given CEF environment. This might not be known, especially for We acknowledge A. S. Sukhanov for stimulating dis- the 4f ions with strong spin-orbit coupling, where CEF cussions. This research used resources at the Spallation determines the magnetic moment. The specific heat mea- Neutron Source, a DOE Office of Science User Facility surements have no such limitation and we suggest this operated by Oak Ridge National Laboratory. S.E.N. ac- method for measurements of impurity concentration for knowledges support from the International Max Planck samples with 4f magnetic ions, especially for the cases Research School for Chemistry and Physics of Quantum where concentration of magnetic ions is relatively low. Materials (IMPRS-CPQM). Laue x-ray diffraction mea- Finally and perhaps most surprisingly, the analysis surements were conducted at the Center for Nanophase of the zero-field magnetic specific heat have shown that Materials Sciences (CNMS) (CNMS2019-R18) at Oak there is no noticeable splitting of the excited triplet state, Ridge National Laboratory (ORNL), which is a DOE which is further corroborated by comparing the obtained Office of Science User Facility. The work at Brookhaven energy of the triplet level with the high-resolution INS National Laboratory was supported by Office of Basic En- result. Hence, we conclude that the dominant, c-axis ergy Sciences (BES), Division of Materials Sciences and exchange interaction of the effective spins-1/2 is sur- Engineering, U.S. Department of Energy (DOE), under prisingly isotropic despite the low site symmetry of Yb contract DE-SC0012704.

[1] L. S. Wu, S. E. Nikitin, Z. Wang, W. Zhu, C. D. Batista, [5] S. Bertaina, S. Gambarelli, A. Tkachuk, I. N. Kurkin, A. M. Tsvelik, A. M. Samarakoon, D. A. Tennant, B. Malkin, A. Stepanov, and B. Barbara, “Rare-earth M. Brando, L. Vasylechko, M. Frontzek, A. T. Savici, solid-state qubits,” Nat. Nanotechnol. 2, 39–42 (2007). G. Sala, G. Ehlers, A. D. Christianson, M. D. Lumsden, [6] C. Simon, M. Afzelius, J. Appel, A. Boyer de la Giro- and A. Podlesnyak, “Tomonaga-Luttinger liquid behavior day, S. J. Dewhurst, N. Gisin, C. Y. Hu, F. Jelezko, and spinon confinement in YbAlO3,” Nat. Commun. 10, S. Kröll,J. H. Müller,J. Nunn, E. S. Polzik, J. G. Rarity, 698 (2019). H. De Riedmatten, W. Rosenfeld, A. J. Shields, N. Sköld, [2] L. S. Wu, S. E. Nikitin, M. Brando, L. Vasylechko, R. M. Stevenson, R. Thew, I. A. Walmsley, M. C. Weber, G. Ehlers, M. Frontzek, A. T. Savici, G. Sala, A. D. Chris- H. Weinfurter, J. Wrachtrup, and R. J. Young, “Quantum tianson, M. D. Lumsden, and A. Podlesnyak, “Antiferro- memories,” Eur. Phys. J. D 58, 1–22 (2010). magnetic ordering and dipolar interactions of YbAlO3,” [7] C.W. Thiel, Thomas Böttger, and R.L. Cone, “Rare- Phys. Rev. B 99, 195117 (2019). earth-doped materials for applications in quantum infor- [3] L. S. Wu, W. J. Gannon, I. A. Zaliznyak, A. M. Tsvelik, mation storage and signal processing,” J. Lumin. 131, M. Brockmann, J.-S. Caux, M. S. Kim, Y. Qiu, J. R. D. 353 – 361 (2011). Copley, G. Ehlers, A. Podlesnyak, and M. C. Aronson, [8] D. D. Awschalom, R. Hanson, J. Wrachtrup, and B. B. “Orbital-exchange and fractional quantum number exci- Zhou, “Quantum technologies with optically interfaced tations in an f-electron metal, Yb2Pt2Pb,” Science 352, solid-state spins,” Nat. Photonics 12, 516–527 (2018). 1206 (2016). [9] N. Kunkel and P. Goldner, “Recent Advances in Rare [4] W. J. Gannon, I. A. Zaliznyak, L. S. Wu, A. E. Feiguin, Earth Doped Inorganic Crystalline Materials for Quantum A. M. Tsvelik, F. Demmel, Y. Qiu, J. R. D. Copley, M. S. Information Processing,” Z. Anorg. Allg. Chem 644, 66– Kim, and M. C. Aronson, “Spinon confinement and a 76 (2018). sharp longitudinal mode in Yb2Pt2Pb in magnetic fields,” [10] S. Olmschenk, D. N. Matsukevich, P. Maunz, D. Hayes, Nat. Commun. 10, 1123 (2019). L.-M. Duan, and C. Monroe, “Quantum Teleportation 8

Between Distant Matter Qubits,” Science 323, 486–489 YbFeO3: Evolution of magnetic excitations through the (2009). spin-reorientation transition,” Phys. Rev. B 98, 064424 [11] H.-J. Lim, S. Welinski, A. Ferrier, P. Goldner, and J. J. L. (2018). Morton, “Coherent spin dynamics of ytterbium ions in [23] M. R. Norman, “Colloquium: Herbertsmithite and the yttrium orthosilicate,” Phys. Rev. B 97, 064409 (2018). search for the quantum spin liquid,” Rev. Mod. Phys. 88, [12] D. Fagundes-Peters, N. Martynyuk, K. Lanstedt, V. Pe- 041002 (2016). ters, K. Petermann, G. Huber, S. Basun, V. Laguta, and [24] M. B. Stone, I. A. Zaliznyak, T. Hong, C. L. Broholm, A. Hofstaetter, “High quantum efficiency YbAG-crystals,” and D. H. Reich, “Quasiparticle breakdown in a quantum J. Lumin. 125, 238 – 247 (2007). spin liquid,” Nature 440, 187–190 (2006). [13] G. Boulon, Y. Guyot, H. Canibano, S. Hraiech, and [25] D. P. Bastien, M. Mourigal, N. B. Christensen, G. J. A. Yoshikawa, “Characterization and comparison of Yb3+- Nilsen, P. Tregenna-Piggott, T. G. Perring, M. Enderle, 3+ doped YAlO3 perovskite crystals (Yb:YAP) with Yb - D. F. McMorrow, D. A. Ivanov, and H. M. Rønnow, doped Y3Al5O12 garnet crystals (Yb:YAG) for laser ap- “Fractional excitations in the square-lattice quantum anti- plication,” J. Opt. Soc. Am. B 25, 884–896 (2008). ferromagnet,” Nat. Phys. 11, 62–68 (2015). [14] K. Xia, R. Kolesov, Ya Wang, P. Siyushev, R. Reuter, [26] C. Kim, A. Y. Matsuura, Z.-X. Shen, N. Motoyama, T. Kornher, N. Kukharchyk, A. D. Wieck, B. Villa, H. Eisaki, S. Uchida, T. Tohyama, and S. Maekawa, S. Yang, and J. Wrachtrup, “All-Optical Preparation “Observation of spin-charge separation in one-dimensional of Coherent Dark States of a Single Rare Earth Ion Spin SrCuO2,” Phys. Rev. Lett. 77, 4054 (1996). in a Crystal,” Phys. Rev. Lett. 115, 093602 (2015). [27] N. Motoyama, H. Eisaki, and S. Uchida, “Magnetic sus- [15] F. Bussières,C. Clausen, A. Tiranov, B. Korzh, V. B. ceptibility of ideal spin-1/2 Heisenberg antiferromagnetic Verma, S. W. Nam, F. Marsili, A. Ferrier, P. Goldner, chain systems, Sr2CuO3 and SrCuO2,” Phys. Rev. Lett. H. Herrmann, C. Silberhorn, W. Sohler, M. Afzelius, 76, 3212 (1996). and N. Gisin, “Quantum teleportation from a telecom- [28] I. A. Zaliznyak, H. Woo, T. G. Perring, C. L. Broholm, wavelength photon to a solid-state quantum memory,” C. D. Frost, and H. Takagi, “Spinons in the Strongly Nat. Photonics 8, 775–778 (2014). Correlated Copper Oxide Chains in SrCuO2,” Phys. Rev. [16] E. Saglamyurek, J. Jin, V. B. Verma, M. D. Shaw, F. Mar- Lett. 93, 087202 (2004). sili, S. W. Nam, D. Oblak, and W. Tittel, “Quantum [29] A. C. Walters, T. G. Perring, J.-S. Caux, A. T. Savici, storage of entangled telecom-wavelength photons in an G. D. Gu, C.-C. Lee, W. Ku, and I. A. Zaliznyak, “Effect erbium-doped optical fibre,” Nat. Photonics 9, 83–87 of covalent bonding on magnetism and the missing neutron (2015). intensity in copper oxide compounds,” Nature Physics 5, [17] T. Zhong, J. M. Kindem, J. G. Bartholomew, J. Rochman, 867 – 872 (2009). I. Craiciu, E. Miyazono, M. Bettinelli, E. Cavalli, [30] https://www.scientificmaterials.com. V. Verma, S. W. Nam, F. Marsili, M. D. Shaw, A. D. [31] G. Ehlers, A. Podlesnyak, J. L. Niedziela, E. B. Iverson, Beyer, and A. Faraon, “Nanophotonic rare-earth quan- and P. E. Sokol, “The new cold neutron chopper spec- tum memory with optically controlled retrieval,” Science trometer at the spallation neutron source: design and 357, 1392–1395 (2017). performance,” Rev. Sci. Instrum. 82, 085108 (2011). [18] J. A. M. Paddison, M. Daum, Z. Dun, G. Ehlers, Y. Liu, [32] G. Ehlers, A. Podlesnyak, and A. I. Kolesnikov, “The cold M. B. Stone, H. Zhou, and M. Mourigal, “Continuous neutron chopper spectrometer at the Spallation Neutron excitations of the triangular-lattice quantum spin liquid Source - A review of the first 8 years of operation,” Rev. YbMgGaO4,” Nat. Phys. 13, 117–122 (2017). Sci. Instrum. 87, 093902 (2016). [19] M. Baenitz, P. Schlender, J. Sichelschmidt, Y. A. Onyki- [33] R. A. Ewings, A. Buts, M. D. Le, J. van Duijn, I. Bustin- ienko, Z. Zangeneh, K. M. Ranjith, R. Sarkar, L. Hozoi, duy, and T. G. Perring, “HORACE: software for the anal- H. C. Walker, J.-C. Orain, H. Yasuoka, J. van den Brink, ysis of data from single crystal spectroscopy experiments H. H. Klauss, D. S. Inosov, and Th. Doert, “NaYbS2:A at time-of-flight neutron instruments,” Nucl. Instrum. planar spin-1/2 triangular-lattice magnet and putative Methods Phys. Res. Sect. A 834, 3132–142 (2016). spin liquid,” Phys. Rev. B 98, 220409(R) (2018). [34] O. Arnold, J. C. Bilheux, J. M. Borreguero, A. Buts, S. I. [20] B. Gao, T. Chen, D. W. Tam, C.-L. Huang, K. Sasmal, Campbell, L. Chapon, M. Doucet, N. Draper, R. Fer- D. T. Adroja, F. Ye, H. Cao, G. Sala, M. B. Stone, raz Leal, M. A. Gigg, V. E. Lynch, A. Markvardsen, C. Baines, J. A. T. Verezhak, H. Hu, J.-H. Chung, X. Xu, D. J. Mikkelson, R. L. Mikkelson, R. Miller, K. Palmen, S.-W. Cheong, M. Nallaiyan, S. Spagna, M. B. Maple, P. Parker, G. Passos, T. G. Perring, P. F. Peterson, S. Ren, A. H. Nevidomskyy, E. Morosan, G. Chen, and P. Dai, M. A. Reuter, A. T. Savici, J. W. Taylor, R. J. Taylor, “Experimental signatures of a three-dimensional quantum R. Tolchenov, W. Zhou, and J. Zikovsky, “Mantid – Data spin liquid in effective spin-1/2 Ce2Zr2O7 pyrochlore,” analysis and visualization package for neutron scattering Nature Phys. 15, 1052–1057 (2019). and µSR experiments,” Nucl. Instrum. Methods Phys. [21] G. Hester, H.S. Nair, T. Reeder, D. R. Yahne, T. N. Res. Sect. A 764, 156 (2014). DeLazzer, L. Berges, D. Ziat, J. R. Neilson, A. A. Aczel, [35] O. Buryy, Ya. Zhydachevskii, L. Vasylechko, D. Sugak, G. Sala, J. A. Quilliam, and K. A. Ross, “Novel Strongly N. Martynyuk, S. Ubizskii, and K. D. Becker, “Ther- Spin-Orbit Coupled Quantum Dimer Magnet: Yb2Si2O7,” mal changes of the crystal structure and the influence of Phys. Rev. Lett. 123, 027201 (2019). thermo-chemical annealing on the optical properties of [22] S. E. Nikitin, L. S. Wu, A. S. Sefat, K. A. Shaykhutdi- YbAlO3 crystals,” J. Phys. Condens. Matter 22, 055902 nov, Z. Lu, S. Meng, E. V. Pomjakushina, K. Conder, (2010). G. Ehlers, M. D. Lumsden, A. I. Kolesnikov, S. Bar- [36] M. A. Noginov, G. B. Loutts, K. Ross, T. Grandy, N. Nogi- ilo, S. A. Guretskii, D. S. Inosov, and A. Podlesnyak, nova, B. D. Lucas, and T. Mapp, “Role of traps in pho- “Decoupled spin dynamics in the rare-earth orthoferrite tocoloration of Mn: YAlO3,” J. Opt. Soc. Am. B 18, 9

931–941 (2001). [43] Stephen W Lovesey, “Theory of neutron scattering from [37] A. Furrer and O. Waldmann, “Magnetic cluster excita- condensed matter,” (Oxford Science Publishers, 1984). tions,” Rev. Mod. Phys. 85, 367 (2013). [44] J. Jensen and A.R. Mackintosh, Rare earth magnetism: [38] B. Leuenberger, A. Stebler, H. U. Güdel,A. Furrer, structures and excitations, International series of mono- R. Feile, and J. K. Kjems, “Spin dynamics of an graphs on physics (Clarendon Press, 1991). isotropic singlet-ground-state antiferromagnet with alter- [45] I. A. Zaliznyak and S.-H. Lee, “Magnetic neutron scatter- nating strong and weak interactions: An inelastic-neutron- ing,” in Modern Techniques for Characterizing Magnetic scattering study of the dimer compound cs3cr2br9,” Phys. Materials, edited by Yimei Zhu (Springer US, 2005) pp. Rev. B 30, 6300–6307 (1984). 3–64. [39] A. Zheludev, G. Shirane, Y. Sasago, M. Hase, and [46] M. Kenzelmann, G. Xu, I. A. Zaliznyak, C. Broholm, K. Uchinokura, “Dimerized ground state and magnetic ex- J. F. DiTusa, G. Aeppli, T. Ito, K. Oka, and H. Takagi, citations in CaCuGe2O6,” Phys. Rev. B 53, 11642–11646 “Structure of end states for a haldane spin chain,” Phys. (1996). Rev. Lett. 90, 087202 (2003). 1 [40] Trimer contribution is below 1 % of the in-filed specific [47] In real structure r12 ≈ (δa, δb, 2 ), with δa ≈ ±0.03 and heat, which is beyond the precision of our measurements. δb ≈ ±0.13. This causes an additional intensity modula- [41] L. S. Wu, S. E. Nikitin, M. Frontzek, A. I. Kolesnikov, tion I(Q) ∝ cos(2πδa) · cos(2πδb). E.g. it is responsible G. Ehlers, M. D. Lumsden, K. A. Shaykhutdinov, E-J. for the second “shadow mode” in YbFeO3 spectrum [22]. Guo, A. T. Savici, Z. Gai, A. S. Sefat, and A. Podlesnyak, However, it has only a minor influence at our range of Q. “Magnetic ground state of the Ising-like antiferromagnet [48] J. T. Haraldsen, T. Barnes, and J. L. Musfeldt, “Neutron DyScO3,” Phys. Rev. B 96, 144407 (2017). scattering and magnetic observables for s = 1/2 spin [42] G. L. Squires, Introduction to the Theory of Thermal clusters and molecular magnets,” Phys. Rev. B 71, 064403 Neutron Scattering, 2012th ed. (Cambridge University (2005). Press, England, 1978). [49] To increase the signal-to-noise ratio we symmetrized the data according to the crystal symmetry and then unfolded it back, which means that data at negative and positive K in Fig.5 (a) are equivalent.