## Quantum Order by Disorder: A Key to Understanding the Magnetic Phases of BaCo (AsO ) 2 4 2

Sangyun Lee, 1, ∗ Shengzhi Zhang, 1 S. M. Thomas, 2 L. Pressley, 3 C. A. Bridges, 3 Eun Sang Choi, 4 Vivien S. Zapf, 1, † Stephen M. Winter, 5, ‡ and Minseong Lee 1, §

1 National High Magnetic Field Laboratory, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA.

2 Los Alamos National Laboratory, Los Alamos, NM 87545, USA.

3 Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA.

4 National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32310-3706, USA.

5 Department of Physics and Center for Functional Materials,

Wake Forest University, Winston-Salem, NC 27109, USA.

BaCo (AsO ) 2 4 2 (BCAO), a honeycomb cobaltate, is considered a promising candidate for materials displaying the Kitaev quantum spin liquid state. This assumption is based on the distinctive characteristics of Co 2+ ions (3 d 7 ) within an octahedral crystal environment, resulting in spin-orbit-coupled J eff = 1/2 doublet states. However, recent experimental observations and theoretical analyses have raised questions regarding this hypothesis. Despite these uncertainties, reports of continuum excitations reminiscent of spinon excitations have prompted further investigations. In this study, we explore the magnetic phases of BCAO under both in-plane and out-of-plane magnetic fields, employing dc and ac magnetic susceptibilities, capacitance, and torque magnetometry measurement. Our results affirm the existence of multiple field-induced magnetic phases, with strong anisotropy of the phase boundaries between in-plane and out-of-plane fields. To elucidate the nature of these phases, we develop a minimal anisotropic exchange model. This model, supported by combined first principles calculations and theoretical modeling, quantitatively reproduces our experimental data. In BCAO, the combination of strong bond-independent XXZ anisotropy and geometric frustration leads to significant quantum order by disorder effects that stabilize colinear phases under both zero and finite magnetic fields.

## I. INTRODUCTION

The interplay between crystal field symmetry, spin-orbit coupling, and the electron configuration of cobalt ions (Co 2+ , 3 d 7 ) in a nearly ideal octahedral oxygen environment stabilizes the J eff = 1/2 state, a fundamental building block of numerous intriguing physical phenomena in various configurations. One of the intriguing theoretical proposals is that cobaltates with a J eff = 1/2 state can potentially exhibit bond-dependent exchange interactions, including what is lately called the Kitaev exchange interaction [1-4]. It is well known that when the Kitaev exchange interaction dominates the spin Hamiltonian, the ground state can become a Kitaev quantum spin liquid (KQSL) due to strong magnetic frustration [5-7]. The KQSL is highly sought after among frustrated magnets because it has the potential to host exotic emergent phenomena such as Majorana Fermions and topological order, with significant applications in quantum computing [8-10]. Therefore, several honeycomb cobalt-based magnets, including Na Co TeO 2 2 6 [11-13], and Na 3 Co SbO 2 6 [14] have been investigated intensively.

Recently, another cobalt-based honeycomb-lattice material BaCo (AsO ) 2 4 2 (BCAO) has been proposed as an ideal platform for realizing a KQSL [15, 16] (Figure 1a). BCAO shows a series of magnetic phase transitions under applied magnetic field ( H ). The zero-field order is identified as either an incommensurate spiral or double-stripe (DS), which is suppressed in favor of a collinear up-up-down (UUD) phase under field [17-20]. These spin structures are illustrated in Figure 1b. For magnetic fields in the crystallographic ab -plane, the latter phase is subsequently suppressed by relatively weak fields (0.5 T), above which a KQSL phase has been suggested [15, 21]. For field in the out-of-plane direction, which is theoretically more conducive to a field-induced KQSL phase [22-25], a broad magnetic continuum was observed [16]. However, conflicting reports exist regarding whether BCAO is best described by an extended Kitaev model with strongly bond-dependent magnetic couplings or as an XXZJ 1 -J 3 model with weakly bond-dependent couplings; recent inelastic neutron scattering measurements, as well as ab initio calculations [26], align more with the XXZJ 1 -J 3 model [27]. It has also been pointed out that strong distortions of oxygen octahedra, whose the trigonal term in crystalfield splitting matrix is larger than one fifth of spin-orbit coupling, lead to Ising- and XY-like spin behavior [4, 28, 29]. Additionally, the weak ligand-assisted exchange interactions weaken bond-dependent exchange interactions [4]. Hence, further investigation is still required to characterize the phases of BCAO comprehensively. Investigating the angular dependence of phase boundaries is an effective method to understand the origin of magnetic anisotropy [30]. Recently, Safari et al. explored the phase diagram of BCAO using magnetotropic susceptibility measurements and theoretical models [31]. Their findings indicate that while classical models qualitatively describe the system's behavior, quantum corrections are essential for precise critical field predictions [31].

∗ Present Address: Department of Physics and National High Magnetic Field Laboratory, University of Florida, Florida 32611, USA.

† Corresponding authors:vzapf@lanl.gov

‡ Corresponding authors:winters@wfu.edu

§ Corresponding authors:ml10k@lanl.gov

In this study, we explored the details of magnetic fieldinduced ground states in BCAO using ac and dc magnetic susceptibility, capacitance measurement, and torque magnetometry measurements at various temperatures ( T ). With the aim of establishing the origin of the magnetic anisotropy, we

FIG. 1. Crystal and magnetic structures. a Crystal structure of BaCo 2 (AsO ) 4 2 showing honeycomb magnetic lattice of Co in the ab -plane layers. Blue balls denote Cobalt, green balls Arsenic and red balls oxygen. We omitted Ba as they exist between layers. The a , b , and c are the global coordinate. The a and b point along two zigzag directions of the honeycomb lattice while b ∗ points along the armchair direction. The c points along the out-of-plane direction. The x, y and z are local coordinates used in Eq. (2), where the x -axis points along the bonddirection while z points along the out-of-plane c axis. The y is perpendicular to the x and the z directions. J 1 and J 3 are the nearest and third nearest interactions, respectively. b Energetically competitive magnetic orders at zero and finite field, with in-plane ordering wavevectors ( Q ) indicated.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00622v1-figure-5613686c-3fb5-432b-8ccb-3741b19a8f23.jpg)

study the angular dependent T H -phase boundaries by tuning the magnetic field orientation. With the application of a magnetic field, multiple-phase transitions were identified for H along the in-plane direction, consistent with previous studies [15-20]. In comparison, for H close to the c -axis, the critical field for suppression of ordering is drastically increased, and intermediate phase transitions disappear or are washed out at the lowest temperatures.

To understand the full details of the phase diagram, we employ first-principles-based calculations to provide a starting point for analyzing the magnetic couplings. The resulting model is refined by fitting the experimental phase diagram employing second-order energy corrections [32] to account for quantum fluctuation effects on the stability of various magnetic phases. This combined approach allows for the precise reproduction of the experimental observations. Consistent with recent assertions [26, 31, 33, 34], we conclude that quantum order-by-disorder (QOBD) effects dominate in BCAO, influencing both the critical field values and specific progression of field-induced phases.

## II. RESULTS

## A. Magnetic properties and phase diagrams for in-plane fields

Figure 2a illustrates the magnetization measured at various temperatures as a function of the field applied along the a -axis. Hysteresis loops emerge between 0.1 T ( H 1 ) and 0.2 T ( H 2 ), which become smaller with increasing temperature and are almost invisible above 4.0 K. Beyond 0.5 T, there is a steep increase in magnetization before reaching saturation around 0.6 T ( H 3 ). The saturation magnetization is about 2.3 µ B/Co 2+ , from which we extracted g ab ∼ 4 6 . . The magnetization between the saturation field and ∼ 0.15 T constitutes a one-third magnetization plateau. These findings are consistent with prior studies [15, 19, 20]. When the temperature exceeds T N ∼ 5.5 K, all features disappear, and a monotonic increase of the magnetization with field is observed, consistent with a thermal paramagnetic state. Figure 2b shows the ac magnetic susceptibility as a function field. The temperature dependent ac magnetic susceptibility is in Supplementary Figure 1. Three peaks are observed: The first two peaks at ∼ H 1 and ∼ H 2 are consistent with the jump in magnetization in the hysteresis loop and the end of the hysteresis loop. The third peak above ∼ H 3 corresponds to the saturation field. ac susceptibility at 6 K does not show any anomalies. The magnetic phase diagram with field along the a -axis deduced from these measurements is presented in Figure 2c. We observed three thermodynamic phases (I, III, and IV) and one intermediate state (II) due to the hysteresis loop [17, 35]. Neutron scattering experiments identified phase I as the DS phase or closely related spiral phase with in-plane wavevector close to Q = [1 / 4 0] , [20, 27]. Phase III is the UUD phase with 1/3 saturation magnetization and Q = [1 / 3 0] , [18-20]. Given the differing wavevectors and the discontinuous change in wavevector through this phase transition, it is not surprising that these phases are separated by a first order phase transition. Phase II corresponds to a coexistence region of Phase I and III. Lastly, phase IV is the asymptotically polarized paramagnetic phase that is smoothly connected to high temperature paramagnetic phase.

Figures 2d-2f show analogous results for field along b ∗ -axis. For this field direction, we observe the same progression of phase transitions. The hysteresis loops are wider but less pronounced, and H 3 is slightly lower in comparison with H a || . Therefore, the width of 1/3 plateau (phase III) is slightly narrower with magnetic field along b ∗ -axis.

FIG. 2. Magnetic properties and phase diagrams for in-plane fields. a Magnetization per Co 2+ at various temperatures with field parallel with a . For the clarity, a 0.5 µ B/Co 2+ offset was applied between the two curves. b Ac susceptibility as a function of field at various temperatures with field parallel with a . c Magnetic field vs. temperature phase diagrams with field along the a -axis. d Magnetization per Co 2+ at various temperatures with field parallel with b ∗ . A 0.5 µ B/Co 2+ offset was applied between the two curves. e Ac susceptibility as a function of field at various temperatures with field parallel with b ∗ direction. f Magnetic field vs. temperature phase diagrams with field along the b ∗ -axis. Data points are taken from dc magnetization M , ac susceptibility χ ′ , and capacitance measurement C (see Supplementary Figure 2) for the phase diagrams.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00622v1-figure-bc425d21-d8c3-4f8d-99fe-8f68af99b1ad.jpg)

## B. Magnetic properties and phase diagrams for out-of-plane field

Figure 3a shows the temperature dependent dc magnetic susceptibility ( M /H c (T)) measured at various constant external magnetic fields. At 1 T, a clear phase transition from the paramagnetic to DS state is observed and marked with an empty circle around 5.3 K. Raising the field to 14 T only suppresses T N to 4 K. The phase transition is therefore very robust against out-of-plane magnetic field, which is in sharp contrast with the behavior for in-plane fields.

Figure 3b shows the magnetization versus out-of-plane field up to 14 T. At 2 K, the magnetization linearly increases without showing any anomaly, demonstrating a gradual spin canting towards the field direction without phase transitions. At 4 K, we observe a weak slope change around 14 T, denoted as a black arrow, which is consistent with the thermal phase transition between the DS and thermal paramagnet observed in Figure 3a. With increasing temperature to 4.5 K, the slope change shifts to slightly lower field of about 11.5 T. Suppression of the phase boundary between AFM and paramagnetic phase with increasing temperature is a typical behavior of antiferromagnets.

Since 14 T is not enough to observe the magnetic saturation, we measured the magnetization with field along c-axis up to

50 T at Pulsed Field Facility in Los Alamos National Laboratory. Results from pulsed field magnetization experiments at various temperature are shown in Figure 3c and 3d. The data are admittedly noisy, but the important features are clearly discernible. In the paramagnetic state at 10 K, the magnetization curve resembles the Brillouin function. Below T N, the magnetization increases linearly and shows a slope change in the vicinity of 25 T, which we assign as the out-of-plane saturation field. Therefore dM /dH c is a constant below ∼ 20 T and above ∼ 35 T. Due to the small crystal field excitation gap of Co 2+ , it is common to observe a linearly increasing magnetization above the saturation field, due to Van Vleck paramagnetism [36, 37]. By extrapolating the Van Vleck contribution to the zero field, we found a saturation moment of about 0.71 µ B /Co 2+ , from which we extract g c ∼ 1 5 . . We note that the saturation field with out-of-plane configuration in this study is higher than the field claimed by Zhang et al. [16] and Safari et al. [31]. As we demonstrated in Supplementary Figure 4 and 5, a tiny misalignment of the crystal can lead to a drastically suppressed saturation field. For example, only 6 ◦ degrees off from the c -axis suppresses the saturation field 27 T to 4.5 T. The 25 T phase transition is smoothly connected to the zero field AFM to paramagnetic phase transition, additionally confirming that it is the saturation field. Thus, the out-of-plane magnetization measurements demon-

FIG. 3. Magnetic properties for out-of-plane fields. a Temperature dependent DC susceptibility at various magnetic field along out-of-plane direction. An offset of 0.001 µ B / Co 2+ /T multiplied by the field is applied to each curve for clarity. b Field dependent magnetization per Co 2+ measured up to 14 T at various different temperatures, with curves offset by 0.05 µ B /Co 2+ . c Pulsed field magnetization up to 50 T at different temperatures. d Differentiated magnetization with respect to the field.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00622v1-figure-86ab095c-34d8-4a25-b178-49c49c0ff8de.jpg)

strate that BCAO is highly anisotropic between in-plane and out-of-plane fields, with a saturation field ratio H c sat /H a sat of about 50.

## C. Angular dependence of magnetic phase transitions

phase transitions to the UUD and polarized phases occur at 1 T and 4 T, respectively. These move towards ∼ H 1 and ∼ H 2 as the field direction approaches b ∗ -axis, consistent with Figure 2d. The magnetic moment also monotonically increases as the field is rotated towards b ∗ direction, consistent with the anisotropy of the g -values.

To investigate how the phase boundaries evolve with the field direction, we measured the magnetization, capacitance, and magnetic torque below T N while varying the magnetic field orientation across the c b -∗ plane. In Figure 4a, we first show the evolution of the magnetization. For an out-of-plane field ( θ = 0 ◦ ), the magnetization does not show any features up to 14 T. However, with a small rotation of the field direction towards the b ∗ direction, the transitions to the intermediate UUD and polarized phases appear. At θ = 6 3 . ◦ , the

Capacitance and magnetostriction measurements can be sensitive probes of magnetic phase transitions as the change in spin structure accompanies the change in electric susceptibility and lattice parameter of magnetic insulators. We found this to be the case for BaCo 2 (AsO ) . 4 2 Figure 4b shows the change in capacitance as a function of magnetic field at various angles. We applied the electric field at 16 kHz along b ∗ -axis. At 27 , the capacitance change shows clear two phase ◦ transitions corresponding to the phase transition from DS to UUD, and UUD to polarized phase. They evolve to higher

a

FIG. 4. Field-angle dependence of the magnetic properties. a Magnetization per Co 2+ as a function of field at various angles at T = 2 K. 0 ◦ is c -axis and 90 ◦ is b ∗ -axis. b Magnetic field dependent capacitance with electric field along b ∗ -axis at different angles at T = 2.5 K. For clarity, we applied the offset of 1 × 10 -5 pF between two curves. The arrows indicate the phase transition fields. c Field dependent magnetic torque at various angles at T = 0.5 K. d Differentiated magnetic torque with respect to the magnetic field at various angles with an offset. The arrows indicate the peak positions. All these measurements were performed with sweep-up magnetic field.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00622v1-figure-a3467217-ab8e-43f9-92bb-c4ed11195512.jpg)

fields as the field becomes close to c -axis. The anomaly from the phase transition from DS to UUD phase becomes weaker while that from UUD to polarized phase is more or less similar at all angles. The capacitance change and magnetostriction data with in-plane magnetic field can be found in Supplementary Figure 2 and 3.

being only 2 T at 19.7 ◦ .

Figure 4c and 4d show magnetic torque ( τ ) and its differentiation with magnetic field ( dτ/dH ) up to 30 T as a function of magnetic field at different fixed angles. Both τ and dτ/dH display a clear anomaly at phase transitions. When the field is only at 0.36 ◦ from c -axis, a sharp increase in τ and a peak dτ/dH appear at about 23 T. This is reasonably well consistent with the saturation field observed in pulsed field magnetization shown in Figure 3c. The anomaly shifts quickly to lower field as the angle becomes larger than 1 ◦ and ends up

In Figure 5a, we show the magnetic torque measured for continuous field-angle sweeps between the c and b ∗ axes at T = 0.5 K at various fixed magnitudes of H . All curves exhibit 180 ◦ periodicity and a saw-tooth shape due to the strong magnetic anisotropy. The sharp change of sign of the torque near 0 ◦ indicates the c -axis is the hard axis. Additional step-like anomalies are observed near θ = 0 ◦ marking the first order transitions H 1 and H 3 . These two features are only clearly observable up to 10 T and they merge into one above 10 T.

In general, the magnetic torque may be expanded as [30]:

$$\underset { \underset { \underset { n = 1 } { \text{gne-}} } { \text{peak} } } } { \text{peak} } \underset { \underset { n = 1 } { \text{gne-}} } { \text{gne-}} } \frac$$

The experimental field-dependence of | c 2 | and | c 4 | is plotted

FIG. 5. Experimental and theoretical angle-dependence magnetic torque and phase diagrams. a b ∗ c -plane magnetic torque as a function of angle for various magnetic fields measured up to 32 T at T = 0 5 . K. b Theoretical simulation of torque on the basis of second-order corrected state energies employing optimal model discussed in the text. c Coefficients | c 2 | and | c 4 | as a function of magnetic field extracted from 0.5 K experimental data. c The same coefficients extracted from 0.0 K theoretical simulations. e Critical magnetic fields as a function of field angle in the b ∗ c -plane. Data points were obtained through the analysis of the magnetization, capacitance, and torque measurements. The red and black symbols correspond to the H 1 and H 3 , respectively. The dashed lines are guides for the eye. f Theoretical phase diagram for optimized model (see text).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00622v1-figure-d09fe6b1-936d-4ed7-aa30-32f248e5b274.jpg)

in Figure 5c. Both coefficients increase steeply up to 4 T and then slowly decrease above 4 T up to the maximum field. The ratio of c 2 to c 4 remains more or less constant at ∼ 0 4 . across the entire field range due to the dominant saw-tooth curve shapes. In the limit of large field, such that the majority of the field angles correspond to the polarized phase, the c 2 n coefficients have two contributions. The first contribution is field-independent, and depends exclusively on the g -anisotropy. For the implied range of g -values ( g ab ∼ 4 6 . , g c = 1 5 . -2 5 . ), the limiting ratio for this contribution leads to | c 4 /c 2 | = 0 17 . -0 27 . . The next order contribution to the c 2 n scales as 1 /H , and is related to both the XXZ exchange anisotropy and g -anisotropy [30]. The findings that c 2 and c 4 decay with increasing field, and that their ratio exceeds the g -only value both demonstrate significant exchange anisotropy is necessary to explain the observed torque. Below, we develop an exchange model for BCAO, and use these observations to verify the model.

From the phase transitions observed in the angular dependence of magnetic properties in Figure 4 and 5a, an angular dependence of the magnetic phase diagram is constructed, as shown in Figure 5e. The red empty symbols represent the phase transition from the DS to the UUD phase ( H 1 ), and the black empty symbols indicate the saturation field ( H 3 ). As the field direction moves from the b ∗ -axis to the c -axis, both phase boundaries sharply increase below 10 ◦ . The saturation field reaches above 23 T, whereas the phase boundary between the DS and the UUD merges with the saturation field within 5 ◦ above about 10 T.

## D. Microscopic magnetic model

To first define a realistic range of magnetic couplings, we employ the des Cloizeaux effective Hamiltonian (dCEH) approach [38-40], based on exact diagonalization of the electronic d -orbital Hamiltonian for clusters include one and two Co sites [see Methods].

We first estimate the g -tensor by diagonalizing the electronic Hamiltonian on one Co site, and projecting the Zeeman operator into the low-energy space of j 1 / 2 states. This yields g ab = 5 1 . and g c = 2 4 . , which further confirms accurate reproduction of the trigonal crystal field splitting in BCAO. To estimate the magnetic couplings, we exactly diagonalize the electronic Hamiltonian on pairs of sites and project the resulting low-energy Hamiltonian onto j 1 / 2 states. The relatively low symmetry crystalline environment allows for a large number of independent exchange constants. In particular, in the R ¯ 3 space group, each first and third neighbor bond is constrained only by inversion symmetry, allowing six independent exchange parameters per bond:

$$\mathcal { H } & = \sum _ { i j } J _ { i j } ^ { x y } ( S _ { i } ^ { x } S _ { j } ^ { x } + S _ { i } ^ { y } S _ { j } ^ { y } ) + J _ { i j } ^$$

where the bond-dependent local x -axis is oriented parallel to the bond, and the z -axis is parallel to the crystallographic c -axis. Second neighbor couplings tend to be weaker and are ignored in first approximation. By varying the strengths of the Coulomb interactions within a realistic range [see Methods], we establish reasonable ranges for the first and third neighbor couplings: J xy 1 ∼ -8 to -12 meV, J z 1 /J xy 1 ∼ 0 5 . to 0.6, J ±± 1 ∼ -0 6 . to -0 8 . meV, Γ yz 1 ∼ 0 35 . to 0.5 meV, J xy 3 /J xy 1 ∼ -0 25 . to -0 4 . , and J z 3 /J xy 3 ∼ 0 2 . . The remaining couplings are all found to have magnitude smaller than 0.1 meV.

There are several aspects of the computed couplings that bear mentioning. In BCAO, the local trigonal distortion of the CoO 6 octahedra has a significant impact on the magnetic Hamiltonian, inducing both the large g -anisotropy, and bondindependent XXZ anisotropy of the intersite exchange. Microscopically, the third neighbor couplings are dominated by a single exchange path involving the e g electrons, which leads to antiferromagnetic interactions that are rendered anisotropic only by the crystal field effects. As such, the off-diagonal anisotropies should be nearly vanishing and the third neighbor interactions should be parameterized by a pure XXZ anisotropy, with J xy 3 being the same sign as J z 3 . For the first neighbor couplings, the dominant anisotropy is also of the bond-independent XXZ type. While Co 2+ materials can, in principle, host strongly bond-dependent couplings [1-4], in BCAO we find that these are suppressed by a combination of factors. First, the large trigonal distortion from a pure octahedral environment significantly modifies the composition of the J eff = 1/2 moments, and imposes significant XXZ anisotropy on all couplings instead of bond-dependent exchange interactions. Second, the largest contribution to the first neighbor exchange comes from ferromagnetic 'chargetransfer' contributions arising from hybridization of the Co e g orbitals with the oxygen p -orbitals. Unlike the couplings between t 2 g electrons, these interactions do not display large bond-dependent anisotropies. As such, the bond-dependent anisotropies ( J ±± 1 and Γ yz 1 ) are subleading. Although some of these features are found in previously proposed models for BCAO [26, 27, 31, 33], not all are reproduced. This raises the question of whether the phase diagram can be understood by employing a model that follows more closely these microscopic considerations.

To address this question, we made a comprehensive search of the parameter space around the predicted range, with a focus on fitting the experimental phase diagram. BCAO necessarily lies in a parameter region with multiple competing phases, making it necessary to consider quantum corrections to the classical state energies. We employ the method outlined in Ref. 32 [see Methods], which considers states defined by classical spin vectors, but accounts for local fluctuation

M

a

g

n

e

t

i

c

M

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00622v1-figure-5c222b89-81af-487f-8dee-7bb560b50de7.jpg)

a

g

n

e

t

i

c

M

a

g

n

e

t

i

c

FIG. 6. Energy per site calculations for the energetically competing phases. a , b Comparison of 2nd order corrected energies per site for different magnetic states for fields H a || and H b || ∗ , respectively. c d , Theoretical magnetization per site for different field directions showing 1/3 magnetization plateau.

effects on the state energies at second order in perturbation theory. This naturally includes the leading quantum orderby-disorder (QOBD) effects. From this search, we propose a simplified, but extremely finely tuned, six-parameter model: J xy 1 = -11 75 . , J z 1 = -5 89 . , J ±± 1 = -0 33 . , Γ yz 1 = -0 15 . , J xy 3 = +3 37 . , and J z 3 = +0 79 . meV. The parameters are more consistent with microscopic considerations but otherwise broadly compatible with previously proposed models for BCAO with dominant XXZ anisotropy [26, 31].

In Figure 6, we show the computed evolution of the state energies and magnetization per site of the above model for inplane fields. We evaluate m = | ∑ i g · S i | from the classical spin vectors for which the second order energies correspond. The zero-field ground state is found to be a colinear DS phase, with spins oriented in the plane, as shown in Figure 1b. Under applied field, there are two first order transitions. The first is to the magnetization plateau UUD state, and the second to the polarized state. We estimate for H || a : H 1 = 0 16 . T, H 3 = 0 55 . T, and for H || b ∗ : H 1 = 0 18 . T, H 3 = 0 52 . T. Thus, the model reproduces the slight anisotropy in the in-field critical fields observed in experiment. The predicted magnetization as a function of field is shown in Figure 6c and 6d, and reproduces well the experimental magnetization curves in Figure 2.

In Figure 5e and 5f, we show a comparison of the experimental and theoretical phase diagram for fields oriented in

Q

Q

Q

Q

M

a

g

n

e

t

i

c

Q

Q

Q

M

a

g

n

e

t

i

z

o

Energy per site (meV)

the b ∗ c -plane. For strictly out-of-plane fields H c || , the model reproduces the direct transition from the DS state to the polarized phase. The critical field depends very strongly on the field angle; for θ = 0 ◦ , we compute H 3 = 36 T, while for θ = 1 ◦ , we compute H 3 = 18 T. This range is comparable to the broadened saturation observed in the pulsed-field measurements in Figure 3. The sharp increase of the critical fields in the vicinity of the c -axis is mainly due to the large bond-independent XXZ anisotropy, which leads to an easy ab -plane. Overall, the model reproduces the phase boundaries well, confirming that the specific progression of field-induced phases can be understood within a microscopically compatible model.

Finally, to further check that the model reproduces the characteristic magnetic anisotropy of BCAO, we simulated the magnetic torque data from the second-order corrected energies, where τ = ∂E/∂θ . Results are shown in Figure 5b alongside the experimental data. The simulations reproduce very well both the positions and magnitudes of the first order jumps in τ corresponding to the various phase boundaries, as well as the rounding of τ ( θ ) at high field. A precise comparison can be made on the basis of the coefficients c 2 n . In Figure 5d, we show the theoretical evolution of these coefficients with field. The finding that c 2 and c 4 are well reproduced by the model implies the relative magnitudes of both the g -anisotropy and XXZ exchange anisotropy are well captured.

Overall, we conclude that the microscopically motivated model provides a very good reproduction of the experimental phase diagram of BCAO and therefore may serve as a starting point for future refinement. The present results confirm the validity of recently proposed models with primarily bondindependent XXZ anisotropy [26, 27, 31], but show that the complex phase diagram of BCAO illuminated by various recent studies is compatible with the microscopic expectations for this material.

## III. DISCUSSION

As a consequence of the exchange anisotropy and frustration between first and third neighbor couplings, the phase diagram of BCAO is uniquely dominated by QOBD effects. Over the region of couplings we considered for BCAO, we observe several trends when comparing classical state energies with the second order corrected energies. First, primarily due to the strong XXZ anisotropy, both antiferromagnetic and polarized (ferromagnetic) states have finite energy corrections. Similar to the DMRG studies in Ref. [31, 34], for the physically relevant magnitudes of the XXZ anisotropy, we find that fluctuations tend to stabilize the polarized phase over the other phases. Overall, we find that quantum fluctuations tend to reduce the critical fields by favoring the polarized state, consistent with the conclusions of [31]. Second, we emphasize that the collinear DS and UUD phases only appear as ground states due to QOBD effects. We did not find any realistic parameters for which these are classical ground states for any magnetic field. Conversely, the spiral phase is never found to be stable once fluctuation effects are included. These findings are consistent with [33, 34]. Finally, we emphasize that the present model for BCAO is extremely finely tuned. As can be anticipated from the (necessarily) small differences of the state energies in Figure 6a and 6b deviations of any coupling on the order of ± 0 01 . meV may be sufficient to ruin the agreement with experiment. This sensitivity likely also leads to variations in the results of different numerical approaches. In this work, we opted to use the inexpensive approximation of second-order energy corrections to facilitate a fine parameter search. Future refinement of the model could employ more sophisticated numerical methods, provided those approaches adequately treat the QOBD effects.

In conclusion, this work is a comprehensive study of the magnetic field-induced ground states of BCAO, employing various techniques such as dc magnetization, capacitance, ac magnetic susceptibility at low temperatures. We mapped out the full T H -phase diagram for different magnetic field orientations, including H along a , b ∗ , c , and c ′ . To understand the specific progression of field-induced phases, we employed first-principles-based calculations to derive a model that reproduces all essential features of the low-temperature phase diagram. The stabilization of two colinear phases (I and III) in both zero and finite fields is attributed to the significant influence of QOBD effects. BCAO presents an intriguing material where quantum fluctuation effects dominate the phase diagram at both zero and finite field. This investigation highlights new avenues for understanding the complex behavior of BCAO under external magnetic fields.

## IV. METHODS

## A. Experimental details

High-quality single crystals of BaCo 2 (AsO ) 4 2 synthesized through a flux method at Oak Ridge National Laboratory in Oak Ridge, for which the detailed methods are described elsewhere [27]. At low temperature, dc magnetization ( M ) is measured up to 60 T through millisecond timescales in pulsed magnet using custom-built 3 He pumping (400 mK-3 K), and ac magnetic susceptibility ( χ ac ) and capacitance ( C ) measurements are shown here performed in PPMS for 2 K-30 K under the magnetic fields up to 14 T at Pulsed Field Facility in Los Alamos. For dc magnetization and ac magnetic susceptibility measurements, the sample with 1 mm 2 was aligned along a , b ∗ , c , and c ′ and mounted on the sample holder inside the compensation and pick-up coil, where c ′ -axis is 6 ◦ from c -axis on the b ∗ -c plane. The electrical contacts were made by painting silver epoxy on the two parallel opposite surfaces for the capacitance measurement and Andeen-Hagerling AH-2700A capacitance bridge was employed to measure the capacitance.

Magnetic torque was measured with the 31 T resistive magnet at the National High Magnetic Field Laboratory in Tallahassee, Florida using a home-made torque magnetometer. The torque magnetometer was formed from two parallel BeCu plates. A tiny aligned crystal with the mass of 0.13 mg was attached using epoxy on the cantilever plates. The capac-

itance change was monitored while the field or angle was swept using an AH-2700A capacitance bridge. More details can be found in https://nationalmaglab.org/user-facilities/dcfield/measurement-techniques/torque-magnetometry-dc/.

We then converted the measured capacitance to the deflected angle ( α ) using the following equation

$$C ( \alpha ) = \frac { \epsilon _ { 0 } \epsilon _ { r } w } { \tan \alpha } \ln \frac { l \tan \alpha - h _ { 0 } } { h _ { 0 } }, \quad \quad ( 3 ) \quad \quad \quad$$

where α is the defected angle that is proportional to the magnetic torque, ϵ 0 and ϵ r are vacuum and relative permittivity, respectively, and w l , , and h 0 are width, length of the BeCu plate and the gap distance between two BeCu plates at zero magnetic field. We extracted the coefficients c 2 n using the following equation

$$c _ { 2 n } = \frac { 2 } { \pi } \int _ { \theta _ { 1 } } ^ { \theta _ { i } + \pi } \frac { \tau } { H } \sin ( 2 n \theta ) d \theta, \quad \quad ( 4 ) \$$

where θ i is the initial angle, τ is the measured magnetic torque signal.

## B. Theoretical phase diagrams

In order to account for quantum fluctuation effects, we employ the method of second-order energy corrections outlined in Ref. 32. At each site i , the local quantization axis (denoted ˆ z ′ i ) is along the local spin vector, such that ⟨ S z ′ i ⟩ = + 1 2 for all sites. Fluctuations are then driven by terms acting like J --ij S -′ i S -′ j . The second-order energy correction due to such fluctuations is:

$$\Delta E ^ { ( 2 ) } = \sum _ { \text{bonds } i j } \frac { J _ { i j } ^ { + + } J _ { i j } ^ { - - } | \langle m | S _ { i } ^ { - \prime } S _ { j } ^ { -$$

where E gs, classical is the classical energy of the spin configuration, and E m, classical is the classical energy of the state with the two spins along the bond ij flipped. For the example of a Heisenberg honeycomb antiferromagnet, the classical ground state energy is E gs, classical = -0 375 . J per site, whereas the present approach provides E gs, classical +∆ E (2) = -0 563 . J per site. The latter estimate (which includes local singlet fluctuations) is in much better agreement with E gs = -0 544 . J per site from series expansion [41] and quantum Monte Carlo [42]. In order to capture the energetically competitive magnetic orders depicted in Figure 1b, we employ a combination of 4-, 6-, and 8-site clusters of spins with appropriate boundary conditions, and numerically minimize E classical +∆ E (2) with respect to the spin orientations. We consider the commensurate spiral state with Q = [1 4 0] / , as a proxy for the proposed incommensurate spiral with Q = [0 27 0] . , [27]. We find the Q = [1 / , 4 0] spiral to be the classical ground state for a range of parameters, but it is completely replaced by the colinear ferromagnetic and Q = [1 / 4 0] , DS phases once ∆ E (2) is included.

## C. First principles calculations

To estimate the magnetic couplings, we consider an electronic model for the Co 3 d -orbitals, which is a sum of, respectively, one- and two-particle terms: H el = H 1 p + H 2 p . The one-particle terms include intersite hopping, intrasite crystal field, and spin-orbit coupling, H 1 p = H hop + H CF + H SO :

$$\underset { i \omega } { \underset { i \omega } { \underset { j \omega } { \leftarrow } } } \quad \mathcal { H } _ { h o p } = \sum _ { i j \alpha \beta \sigma } t _ { i j } ^ { \alpha \beta } c _ { i, \alpha, \sigma } ^ { \dagger } c _ { j, \beta, \sigma }$$

$$\overset { \text{ity} } { \text{Cu} } _ { \text{ero} } \quad \mathcal { H } _ { \text{CF} } = \sum _ { i \alpha \beta \sigma } d _ { i } ^ { \alpha \beta } c _ { i, \alpha, \sigma } ^ { \dagger } c _ { i, \beta, \sigma }$$

$$\text{the} \quad \mathcal { H } _ { \text{SO} } = \sum _ { i \alpha \beta \sigma \sigma ^ { \prime } } \lambda _ { \text{Mn} } \langle \phi _ { i } ^ { \alpha } ( \sigma ) | \mathbf L \cdot \mathbf S | \phi _ { i } ^ { \beta } ( \sigma ^ { \prime } ) \rangle c _ { i, \alpha, \sigma } ^ { \dagger } c _ { i, \beta, \sigma ^ { \prime } } \ ( 8 )$$

where c † i,α,σ creates an electron at Co site i , in d -orbital α , with spin σ . To parameterize the one-particle terms, we performed fully relativistic density functional theory (DFT) calculations at the GGA (PBE [43]) level using FPLO [44], a 12 × 12 × 12 k -grid, and the structure from Ref. 45. The resulting Kohn-Sham Bloch functions were projected onto Co d -orbitals to yield single-particle terms from the Kohn-Sham Fock matrix written in the Wannier basis [46]. For BCAO, this approach yields accurate crystal-field excitation energies [47], confirming the validity of the computed single-particle parameters.

For the two-particle terms, we consider a combination of on-site and nearest neighbor intersite contributions: H 2 p = H U + H nn J , respectively. The on-site contributions are given by:

$$\mathcal { H } _ { U } = \sum _ { i \alpha \beta \delta \gamma } \sum _ { \sigma \sigma ^ { \prime } } U _ { \alpha \beta \gamma \delta } \, c _ { i, \alpha, \sigma } ^ { \dagger } c _ { i, \beta, \sigma ^ { \prime } } ^ { \dagger } c _ { i, \gamma, \sigma ^ { \prime } } c _ { i, \delta, \sigma } \quad ( 9 )$$

In the spherically symmetric approximation [48], U αβγδ are parameterized by the Slater parameters F dd 0 , F dd 2 , F dd 4 . For the (screened) on-site metal Coulomb interactions, we use F dd 2 = 11 55 . eV, and F dd 4 = 7 68 . eV, following the fitted parameters for LiCoO 2 from Ref. 49. We consider U t 2 g = F dd 0 +(4 49)( / F dd 2 + F dd 4 ) within the realistic range of 5.0 to 5.5 eV.

The nearest neighbor intersite Hund's coupling is given by:

$$\mathcal { H } _ { J } ^ { \mathrm n n } = \sum _ { i j \alpha \beta \sigma \sigma ^ { \prime } } J _ { H, i j } ^ { \alpha \beta } \, c _ { i, \alpha, \sigma } ^ { \dagger } c _ { j, \beta, \sigma ^ { \prime } } ^ { \dagger } c _ { i, \alpha, \sigma ^ { \prime } } c _ { j, \beta, \sigma } \quad ( 1 0 )$$

A significant contribution to the intersite magnetic couplings comes from ligand 'charge transfer' exchange terms [1, 4], which account for perturbative processes where two holes meet on a single oxygen ligand. When downfolded into the d -orbital Wannier basis, these appear as intersite Hund's coupling terms. They give an overall ferromagnetic contribution to the intersite magnetic couplings between the low-energy j 1 / 2 doublets, which may be rendered anisotropic (bond-independent XXZ anisotropy) by the crystal-field splitting of the t 2 g orbitals. In a previous study of BCAO [26],

some of the authors accounted for this term via an estimated correction added to the computed couplings. Here, we employ the improved approach utilized in [40, 50, 51] that captures these effects at the level of the electronic Hamiltonian. The d -orbital Wannier functions can be written:

$$| \tilde { d } _ { i, \alpha, \sigma } \rangle = \phi _ { i, \alpha } ^ { 0 } | d _ { i, \alpha, \sigma } \rangle + \sum _ { n, \beta } \phi _ { i, \alpha } ^ { n, a } | p _ { n, a, \sigma } \rangle \quad ( 1 1 ) \quad \text{and} \quad \text{(11)}$$

where i labels a given Co metal site, α is a d -orbital, n labels an oxygen ligand site, and a ∈ { x, y, z } is a p -orbital. The tilde on the left side of the equation indicates a Wannier function, whereas the absence of a tilde indicates a bare atomic orbital. Rotating the p -orbital Coulomb terms into the ˜ d Wannier basis results in residual intersite Hund's coupling given by:

$$J ^ { \alpha \beta } _ { H, i j } = \sum _ { n } \sum _ { \alpha \beta } & \left ( U ^ { p } - J ^ { p } _ { H } \right ) \phi ^ { n, a } _ { i, \alpha } \phi ^ { n, b } _ { i, \alpha } \phi ^ { n, a } _ { j, \beta } \phi ^ { n, b } _ { j, \beta } \\ & + J ^ { p } _ { H } \, \phi ^ { n, a } _ { i, \alpha } \phi ^ { n, a } _ { i, \alpha } \phi ^ { n, b } _ { j, \beta } \phi ^ { n, b } _ { j, \beta }$$

where U p and J p H are the atomic ligand Hubbard and Hund's couplings, respectively. The computed magnetic couplings depend most strongly on the overall magnitude of J αβ ij , rather than the specific orbital dependence, which provides for some inconsequential flexibility in the choice of U p and J p H . We take U p = 2 J p H , and consider J p H within the range of 0.4 to 0.6 eV. This choice corresponds to roughly 60% to 80% of the atomic value for oxygen [52]. The sum over oxygen sites in Eq. (12) is then evaluated from the ab-initio Wannier function coefficients for the bridging oxygen atoms for each bond.

## ACKNOWLEDGEMENTS

The authors would like to thank A. L. Chernyshev, P. A. Maksimov, Y. B. Kim, F. Desrochers, E. Z. Zhang, K. A. Modic, and S. A. Crooker for insightful discussions. The experimental work, analysis and manuscript writing (S.L., S.Z., S.M.T., L.P., C.A.B., V.S.Z., and M.L.) was funded by the US Department of Energy, Office of Science, National Quantum Information Science Research Centers. The theoretical work, analysis and manuscript writing of S. M. W. was funded by the Oak Ridge Associated Universities (ORAU) through the Ralph E. Powe Junior Faculty Enhancement Award. Computations were performed using the Wake Forest University (WFU) High-Performance Computing Facility [53], a centrally managed computational resource available to WFU researchers including faculty, staff, students, and collaborators. A portion of this work was performed at the National High Magnetic Field Laboratory, which is supported by National Science Foundation Cooperative Agreement No. DMR2128556, the State of Florida, and the U.S. Department of Energy.

## AUTHOR CONTRIBUTIONS

V.S.Z. and M.L. conceived and designed the project. S.L., S.Z., S.M.T., E.S.C. and M.L. performed measurements. S.M.W. performed the theoretical calculations. L.P. and C.A.B. prepared high quality samples. S.L., V.S.Z., S.M.W., and M.L. wrote the manuscript, with contributions from all authors.

## COMPETING INTERESTS

The authors declare no competing interests.

Correspondence and requests for materials should be addressed to Vivien S. Zapf, Stephen M. Winter or Minseong Lee.

## V. SUPPLEMENTARY MATERIAL 'QUANTUM ORDER BY DISORDER: A KEY TO UNDERSTANDING THE MAGNETIC PHASES OF BACO (ASO ) ' 2 4 2

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00622v1-figure-3861c358-f4e3-4754-a751-e7329a3f23e3.jpg)

Supplemental Figure 1. In-plane magnetic susceptibility of BaCo (AsO ) 2 4 2 Temperature-dependent ac magnetic susceptibility at various magnetic field along in-plane direction with offsets for magnetic field along a ( a ) and b ∗ -axis ( b ), respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00622v1-figure-e77a96f3-b743-4a59-a422-a662ae33e012.jpg)

Supplemental Figure 2. In-plane capacitance and magnetization of BaCo (AsO ) 2 4 2 Temperature-dependent capacitance ( a ) and magnetization ( b ) at various magnetic field with offsets for magnetic field along b ∗ -axis. The ac electrical field excitation is applied along c -axis.

Supplemental Figure 3. Thermal expansion and magnetostriction for H b || ∗ (a) Normalized thermal expansion along c -axis at various magnetic field for H b || ∗ . (b) Magnetic field-dependent capacitance (open circles) and normalized magnetostriction along c -axis (solid lines) are plotted together at 2 K for H b || ∗ . For capacitance, electrical field is applied along c -axis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00622v1-figure-f4b0f0fe-8f37-47d0-839e-5762343ef7ce.jpg)

Supplemental Figure 4. Magnetization per Co 2+ for H || c ′ , 6 ◦ away from c -axis on b ∗ c plane. (a) Magnetization at various temperatures with field (b) Magnetization with field along H c || ′ at various magnetic field with temperature.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00622v1-figure-8201096f-999b-4001-88aa-6575da4b8a55.jpg)

Supplemental Figure 5. Phase diagram for H || c and c ′ . c ′ -axis is 6 ◦ away from c -axis on b ∗ -c plane. Magnetic field vs. temperature phase diagrams with field along c -axis (a) and c ′ -axis (b) , respectively. Data points are taken from dc magnetization M (blue symbols), ac susceptibility χ (red symbols), and capacitance measurement C (black symbols).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00622v1-figure-0190d95e-89c2-479f-a09f-7b684c399844.jpg)

- [1] H. Liu and G. Khaliullin, Pseudospin exchange interactions in d 7 cobalt compounds: Possible realization of the Kitaev model, Phys. Rev. B 97 , 014407 (2018).
- [2] R. Sano, Y. Kato, and Y. Motome, Kitaev-Heisenberg Hamiltonian for high-spin d 7 Mott insulators, Phys. Rev. B 97 , 014408 (2018).
- [3] H. Liu, J. c. v. Chaloupka, and G. Khaliullin, Kitaev Spin Liquid in 3 d Transition Metal Compounds, Phys. Rev. Lett. 125 , 047201 (2020).
- [4] S. M. Winter, Magnetic couplings in edge-sharing high-spin d 7 compounds, J. Phys. Mater. 5 , 045003 (2022).
- [5] G. Jackeli and G. Khaliullin, Mott Insulators in the Strong SpinOrbit Coupling Limit: From Heisenberg to a Quantum Compass and Kitaev Models, Phys. Rev. Lett. 102 , 017205 (2009).
- [6] J. G. Rau, E. K.-H. Lee, and H.-Y. Kee, Generic Spin Model for the Honeycomb Iridates beyond the Kitaev Limit, Phys. Rev. Lett. 112 , 077204 (2014).
- [7] Chaloupka, Jiˇ r´ ı and Jackeli, George and Khaliullin, Giniyat, Kitaev-Heisenberg Model on a Honeycomb Lattice: Possible Exotic Phases in Iridium Oxides A 2 IrO 3 , Phys. Rev. Lett. 105 , 027204 (2010).
- [8] A. Kitaev, Anyons in an exactly solved model and beyond, Ann. Phys. 321 , 2 (2006).
- [9] C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das Sarma, Non-Abelian anyons and topological quantum computation, Rev. Mod. Phys. 80 , 1083 (2008).
- [10] A. Y. Kitaev, Fault-tolerant quantum computation by anyons, Ann. Phys. 303 , 2 (2003).
- [11] G. Lin, J. Jeong, C. Kim, Y. Wang, Q. Huang, T. Masuda, S. Asai, S. Itoh, G. G¨nther, M. Russina, et al., Field-induced u quantum spin disordered state in spin-1/2 honeycomb magnet Na Co TeO , Nat. Commun. 2 2 6 12 , 5559 (2021).
- [12] S. Zhang, S. Lee, A. J. Woods, W. K. Peria, S. M. Thomas, R. Movshovich, E. Brosha, Q. Huang, H. Zhou, V. S. Zapf, and M. Lee, Electronic and magnetic phase diagrams of the Kitaev quantum spin liquid candidate Na Co TeO 2 2 6 , Phys. Rev. B 108 , 064421 (2023).
- [13] S. Zhang, S. Lee, E. Brosha, Q. Huang, H. Zhou, V. S. Zapf, and M. Lee, Out-of-plane magnetic phase diagram of Kitaev quantum spin liquid candidate Na Co TeO , arXiv preprint 2 2 6 arXiv:2405.13935 (2024).
- [14] J.-Q. Yan, S. Okamoto, Y. Wu, Q. Zheng, H. D. Zhou, H. B. Cao, and M. A. McGuire, Magnetic order in single crystals of Na Co SbO 3 2 6 with a honeycomb arrangement of 3d 7 Co 2+ ions, Phys. Rev. Mater. 3 , 074405 (2019).
- [15] R. Zhong, T. Gao, N. Ong, and R. Cava, Weak-field induced nonmagnetic state in a Co-based honeycomb, Sci. Adv. 6 , 1 (2020).
- [16] X. Zhang, Y. Xu, T. Halloran, R. Zhong, C. Broholm, R. Cava, N. Drichko, and N. Armitage, A magnetic continuum in the cobalt-based honeycomb magnet BaCo (AsO ) , Nat. Mater. 2 4 2 22 , 58 (2023).
- [17] P. Mukharjee, B. Shen, S. Erdmann, A. Jesche, J. Kaiser, B. PR, O. Zaharko, P. Gegenwart, and A. Tsirlin, Intermediate fieldinduced phase of the honeycomb magnet BaCo2(AsO4)2, arXiv preprint arXiv:2403.04466v1 (2024).
- [18] L. Regnault, P. Burlet, and J. Rossat-Mignod, Magnetic ordering in a planar XY model: BaCo2(AsO4)2, Physica B+C 86 , 660 (1977).
- [19] L. Regnault and J. Rossat-Mignod, Effect of a magnetic field on the magnetic ordering of BaCo2(AsO4)2, J. Magn. Magn.

Mater. 14 , 194 (1979).

- [20] L.-P. Regnault, C. Boullier, and J. Lorenzo, Polarized-neutron investigation of magnetic ordering and spin dynamics in BaCo2(AsO4)2 frustrated honeycomb-lattice magnet, Heliyon 4 (2018).
- [21] C. Tu, D. Dai, X. Zhang, C. Zhao, X. Jin, B. Gao, T. Chen, P. Dai, and S. Li, Evidence for gapless quantum spin liquid in a honeycomb lattice, arXiv preprint arXiv:2212.07322 (2022).
- [22] C. Hickey and S. Trebst, Emergence of a field-driven U(1) spin liquid in the Kitaev honeycomb model, Nat. Commun. 10 , 530 (2019).
- [23] Y.-F. Jiang, T. P. Devereaux, and H.-C. Jiang, Field-induced quantum spin liquid in the Kitaev-Heisenberg model and its relation to α -RuCl 3 , Phys. Rev. B 100 , 165123 (2019).
- [24] N. Patel and N. Trivedi, Magnetic field-induced intermediate quantum spin liquid with a spinon Fermi surface, Proc. Natl. Acad. Sci. 116 , 12199 (2019).
- [25] Gordon, Jacob S and Catuneanu, Andrei and Sørensen, Erik S and Kee, Hae-Young, Theory of the field-revealed kitaev spin liquid, Nat. Commun. 10 , 2470 (2019).
- [26] P. A. Maksimov, A. V. Ushakov, Z. V. Pchelkina, Y. Li, S. M. Winter, and S. V. Streltsov, Ab initio guided minimal model for the 'Kitaev' material BaCo 2 ( AsO ) 4 2 : Importance of direct hopping, third-neighbor exchange, and quantum fluctuations, Phys. Rev. B 106 , 165131 (2022).
- [27] T. Halloran, F. Desrochers, E. Z. Zhang, T. Chen, L. E. Chern, Z. Xu, B. Winn, M. Graves-Brook, M. Stone, A. I. Kolesnikov, et al., Geometrical frustration versus Kitaev interactions in BaCo (AsO ) , 2 4 2 Proc. Natl. Acad. Sci. 120 , e2215509119 (2023).
- [28] M. E. Lines, Magnetic properties of co cl 2 and ni cl 2 , Phys. Rev. 131 , 546 (1963).
- [29] T. Oguchi, Theory of magnetism in cocl2· 2h2o, J. Phys. Soc. Jpn. 20 , 2236 (1965).
- [30] K. Riedl, Y. Li, S. M. Winter, and R. Valent´ ı, Sawtooth Torque in Anisotropic j eff = 1 2 / Magnets: Application to α -RuCl 3 , Phys. Rev. Lett. 122 , 197202 (2019).
- [31] S. Safari, W. Bateman-Hemphill, A. Mitra, F. Desrochers, E. Z. Zhang, L. Shafeek, A. Ferrenti, T. M. McQueen, A. Shekhter, Z. K¨ oll¨, o et al., Quantum Fluctuations Suppress the Critical Fields in BaCo 2 (AsO ) , arXiv preprint arXiv:2403.15315 4 2 (2024).
- [32] G. Jackeli and A. Avella, Quantum order by disorder in the Kitaev model on a triangular lattice, Phys. Rev. B 92 , 184416 (2015).
- [33] P. A. Maksimov, Proximity-induced sequence of field transitions in the Kitaev candidate BaCo 2 (AsO ) , Phys. Rev. B 4 2 108 , L180405 (2023).
- [34] S. Jiang, S. R. White, and A. L. Chernyshev, Quantum phases in the honeycomb-lattice J 1 -J 3 ferro-antiferromagnetic model, Phys. Rev. B 108 , L180406 (2023).
- [35] X. Wang, R. Sharma, P. Becker, L. Bohat´ y, and T. Lorenz, Single-crystal study of the honeycomb XXZ magnet BaCo 2 ( PO 4 2 ) in magnetic fields, Phys. Rev. Mater. 7 , 024402 (2023).
- [36] M. Lee, J. Hwang, E. S. Choi, J. Ma, C. R. Dela Cruz, M. Zhu, X. Ke, Z. L. Dun, and H. D. Zhou, Series of phase transitions and multiferroicity in the quasi-two-dimensional spin1 2 triangular-lattice antiferromagnet Ba CoNb O 3 2 9 , Phys. Rev. B 89 , 104420 (2014).

- [37] T. Susuki, N. Kurita, T. Tanaka, H. Nojiri, A. Matsuo, K. Kindo, and H. Tanaka, Magnetization Process and Collective Excitations in the S =1 2 / Triangular-Lattice Heisenberg Antiferromagnet Ba CoSb O 3 2 9 , Phys. Rev. Lett. 110 , 267201 (2013).
- [38] J. Des Cloizeaux, Extension d'une formule de Lagrange ` a des probl` emes de valeurs propres, Nucl. Phys. 20 , 321 (1960).
- [39] K. Riedl, Y. Li, R. Valent´ ı, and S. M. Winter, Ab Initio approaches for low-energy spin Hamiltonians, Phys. Status Solidi B 256 , 1800684 (2019).
- [40] R. Dhakal, S. Griffith, K. Choi, and S. M. Winter, Spin-Phonon Coupling in Transition Metal Insulators: General Computational Approach and Application to MnPSe , arXiv preprint 3 arXiv:2407.00659 (2024).
- [41] J. Oitmaa, C. J. Hamer, and Z. Weihong, Quantum magnets on the honeycomb and triangular lattices at T =0, Phys. Rev. B 45 , 9834 (1992).
- [42] U. L¨ ow, Properties of the two-dimensional spin-1/2 Heisenberg model on a honeycomb lattice with interlayer coupling, Condensed Matter Physics (2009).
- [43] J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized Gradient Approximation Made Simple, Phys. Rev. Lett. 77 , 3865 (1996).
- [44] I. Opahle, K. Koepernik, and H. Eschrig, Full-potential bandstructure calculation of iron pyrite, Phys. Rev. B 60 , 14035 (1999).
- [45] T. Dordevi´ c, BaCo (AsO ) , 2 4 2 Acta Crystallogr. E 64 , i58 (2008).
- [46] K. Koepernik, O. Janson, Y. Sun, and J. van den Brink, Symmetry-conserving maximally projected Wannier functions, Phys. Rev. B 107 , 235135 (2023).
- [47] B. S. Mou, X. Zhang, L. Xiang, Y. Xu, R. Zhong, R. J. Cava, H. Zhou, Z. Jiang, D. Smirnov, N. Drichko, et al., Comparative Raman Scattering Study of Crystal Field Excitations in Co-based Quantum Magnets, arXiv preprint arXiv:2403.11980 (2024).
- [48] S. Sugano, Y. Tanabe, and H. Kamimura, eds., Multiplets of Transition-Metal Ions in Crystals, Pure and Applied Physics, Vol. 33 (Elsevier, 1970) pp. 294-301.
- [49] J. van Elp, J. L. Wieland, H. Eskes, P. Kuiper, G. A. Sawatzky, F. M. F. de Groot, and T. S. Turner, Electronic structure of CoO, Li-doped CoO, and LiCoO 2 , Phys. Rev. B 44 , 6090 (1991).
- [50] R. Dhakal, S. Griffith, and S. M. Winter, Hybrid SpinOrbit Exciton-Magnon Excitations in FePS 3 , arXiv preprint arXiv:2403.05958 (2024).
- [51] A. A. Konieczna, D. A. Kaib, S. M. Winter, and R. Valent´ ı, Understanding the microscopic origin of the magnetic interactions in CoNb O , arXiv preprint arXiv:2406.17854 2 6 (2024).
- [52] T. Anno and H. Teruya, Systematic determination of the SlaterCondon parameters of atoms and ions with 1s 2s 2 m 2p n configurations, Theor. Chim. Acta 21 , 127 (1971).
- [53] I. Systems and W. F. University, WFU High Performance Computing Facility (2021).
Sangyun Lee Shengzhi Zhang S. M. Thomas L. Pressley C. A. Bridges Eun Sang Choi Vivien S. Zapf Stephen M. Winter Minseong Lee
Quantum Order by Disorder: A Key to Understanding the Magnetic Phases of BaCo$_2$(AsO$_4$)$_2$
BaCo$_2$(AsO$_4$)$_2$ (BCAO), a honeycomb cobaltate, is considered a
promising candidate for materials displaying the Kitaev quantum spin liquid
state. This assumption is based on the distinctive characteristics of Co$^{2+}$
ions (3$d^7$) within an octahedral crystal environment, resulting in
spin-orbit-coupled $J_{\rm eff}$~=~1/2 doublet states. However, recent
experimental observations and theoretical analyses have raised questions
regarding this hypothesis. Despite these uncertainties, reports of continuum
excitations reminiscent of spinon excitations have prompted further
investigations. In this study, we explore the magnetic phases of BCAO under
both in-plane and out-of-plane magnetic fields, employing dc and ac magnetic
susceptibilities, capacitance, and torque magnetometry measurement. Our results
affirm the existence of multiple field-induced magnetic phases, with strong
anisotropy of the phase boundaries between in-plane and out-of-plane fields. To
elucidate the nature of these phases, we develop a minimal anisotropic exchange
model. This model, supported by combined first principles calculations and
theoretical modeling, quantitatively reproduces our experimental data. In BCAO,
the combination of strong bond-independent XXZ anisotropy and geometric
frustration leads to significant quantum order by disorder effects that
stabilize collinear phases under both zero and finite magnetic fields.