## Origin of magnetic anisotropy in the spin ladder compound (C H 5 12 N) CuBr 2 4

D. Blosser, 1, ∗ V. K. Bhartiya, 1 D. J. Voneshen, 2 and A. Zheludev 1, †

1

Laboratory for Solid State Physics, ETH Z¨ urich, 8093 Z¨rich, u Switzerland 2 ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxon OX11 0QX, United Kingdom (Dated: October 7, 2019)

The S = 1 2 spin ladder compound (C H / 5 12 N) CuBr 2 4 (BPCB) is studied by means of highresolution inelastic neutron scattering. In agreement with previous studies we find a band of triplet excitations with a spin gap of ∼ 0 8 meV and a bandwidth of . ∼ 0 6 meV. In addition, we observe a . distinct splitting of the triplet band of 50(1) µ eV or 40(2) µ eV at the band minimum or maximum, respectively. By comparison to a strong coupling expansion calculation of the triplet dispersion for a spin ladder with anisotropic exchange, weakly anisotropic leg interactions are identified as the dominant source of magnetic anisotropy in BPCB. Based on these results, we discuss the nature of magnetic exchange anisotropy in BPCB and in related transition-metal insulators.

## I. INTRODUCTION

The Heisenberg model is arguably the most important construct in quantum magnetism. Experimentally, it is usually studied in transition metal oxides and metal-organic salts. Unfortunately, symmetry-breaking magnetic anisotropy is unavoidable in such compounds. Even when weak compared to the dominant Heisenberg exchange, it can have a significant influence on the magnetic and dynamical properties. In understanding these, existing microscopic theories are usually disregarded in favor of minimalistic empirical models. For example, magnetic exchange interactions in insulators are described within Anderson's theory of super exchange 1-3 . In the presence of weak spin-orbit-coupling it predicts both symmetric (Ising) and antisymmetric (Dzyaloshinskii-Moriya, DM) exchange anisotropy always occurring in combination and at a fixed ratio 4-6 . Yet, experimental data is often successfully analyzed in comparison to a Heisenberg model with only an added DM term 7-11 even though such a model lacks microscopic justification. Even more often, signatures of weak magnetic anisotropy are detected in experiments, but their origin remains unknown altogether 12-14 . In short, anisotropy is ubiquitously present, yet pinpointing its microscopic origin is almost never possible.

Here we present a detailed case study: the effect of very weak magnetic anisotropy on spin excitations in the almost ideal Heisenberg quantum spin ladder compound (C H 5 12 N) CuBr 2 4 (BPCB) 15-21 . On this compound we have obtained exceptionally high resolution inelastic neutron scattering data showing a small but distinct anisotropy splitting of the triplet excitations as described in Sec. II. At the same time, BPCB is an excellent realization of a strong rung spin ladder. It is described by only two exchange couplings: the rung exchange and an almost four times smaller leg coupling. For this model, the excitation spectrum in the presence of various possible exchange anisotropies can be precisely calculated as detailed in Sec. III. This combination of high resolution spectroscopic data and precise calculations of the excitation spectrum allows a detailed discussion of the anisotropy's microscopic origins in BPCB

FIG. 1. Schematic of the crystal structure of BPCB. The spin ladders are formed by the magnetic Cu 2+ cations and linking Br -anions. Crystallographic centers of inversion symmetry are marked by green dots. Anisotropic Dzyaloshinskii-Moriya interactions, by symmetry are only allowed on the ladder legs and possible DM vectors D are sketched as blue arrows on these bonds. Note that for this monoclinic structure, in both plots the third axis (projection axis) is not perpendicular to the plane of the figure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00724v2-figure-1c69cc51-5992-4f94-b6df-3db93e81d70a.jpg)

(Sec. IV). These findings we compare to the microscopic theory of anisotropic superexchange and we discuss how the present case study may guide our understanding of magnetic anisotropy in related compounds.

## II. EXPERIMENT AND RESULTS

## A. The compound (C H 5 12 N) CuBr : crystal 2 4 structure and magnetic interaction pathways

The compound (C H 5 12 N) CuBr , 2 4 bis-piperidinium copper bromide or BPCB for short crystallizes in a monoclinic crystal structure with space group P2 /c 1 and room temperature lattice parameters a = 8 49, . b = 17 22, . c = 12 38 ˚ , . A and β = 99 3 . ◦ 15 . The spin ladders are formed by the magnetic Cu 2+ cations car-

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00724v2-figure-e690cf81-3fad-45db-838b-c55bc964136b.jpg)

FIG. 2. False-color maps of the inelastic neutron scattering intensity measured using a) E i = 2 20 meV and b) . E i = 1 35 meV . incident energy neutrons at T = 0 35 K. With the lower incident energy, a superior energy resolution is achieved but only the . low energy parts of the dispersion can be probed.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00724v2-figure-f8b767fd-f75e-4347-9385-b8d87f989ec3.jpg)

rying a spin S = 1 2 linked by super-exchange bridges / via Br -anions 16 . The ladders run along the crystallographic a axis and are well separated by non-magnetic organic piperidinium ions as depicted in Fig. 1 22

There are two crystallographically equivalent spin ladders in BPCB, related by the glide plane symmetry of the P2 /c 1 space group. Furthermore, the center of each ladder rung corresponds to a crystallographic center of inversion symmetry. For this reason, on the ladder rungs, only symmetric exchange is possible 23 , whilst for the ladder legs there are no symmetry restrictions and both symmetric and antisymmetric exchange are in principle allowed. On the ladder leg, the antisymmetric exchange contribution is parameterized by a Dzyaloshinskii-Moriya vector D . These DM vectors are uniform within every ladder leg and anti-aligned between the two legs of each ladder. The DM vectors of the two crystallographically equivalent ladders are again related by the reflection of the glide plane symmetry as sketched in Fig. 1. Besides these relations, the DM vectors may point in any direction and are not further constrained by symmetry. Similarly, The g -tensors for the magnetic moments residing on the different Cu 2+ ions show different principal axes for the two types of ladders. Thus, in an applied magnetic field, the two types of ladders become inequivalent, except for special orientations of the magnetic field 19,20 . In the present study, however, we exclusively focus on zero-field properties where the two ladders remain equivalent.

## B. Inelastic neutron scattering

Neutron scattering experiments were performed on four fully deuterated single crystals of total mass 2.07 g, co-aligned to better than 1 ◦ effective mosaic spread. The measurements were carried out at the LET cold neutron time-of-flight spectrometer 24 at the ISIS facility, UK.

FIG. 3. Cuts through the data shown in Fig. 2 at the band extrema: a) Cut at the band minimum q ‖ / π 2 = -0 5 . through the 1 35 meV data set. . b,c) Cuts at the band extrema q ‖ / π 2 = -0 5 . or q ‖ / 2 π = 0, respectively, through the 2 20 meV data set. . Solid lines correspond to Gaussian fits. At the band minimum we find a slightly larger splitting of ∆ min = 50(1) µ eV than at the band maximum where ∆ max = 40(2) µ eV.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00724v2-figure-7b638cb5-ec93-4e34-8742-dce3a75a6584.jpg)

The sample was mounted on a 3 He- He dilution refrig4 erator with the crystallographic b axis vertical. Making use of repetition rate multiplication, data was collected simultaneously using neutrons of incident energies E i = 1.35, 2.20, 4.20 and 11.0 meV. For these configurations we find an approximately Gaussian energy resolution of 19, 36, 97 and 410 µ eV at FWHM, respectively, for elastic scattering ( glyph[planckover2pi1] ω = 0) and improving towards higher energy transfer.

## C. Experimental results

An overview of the neutron scattering data 25 collected for E i = 2 2 and 1.35 meV at a temperature of 0 35 K is . . presented as false color intensity plots in Fig. 2. For this quasi-one-dimensional spin ladder system, the measured scattering intensity is plotted versus reduced wave-vector transfer parallel to the ladder q ‖ = Q a · and integrated fully along the non-dispersive b ∗ and c ∗ directions.

We observe a band of excitations with an excitation gap of ∆ ∼ 0 8 meV and a bandwidth of Γ . ∼ 0 6 meV. . The line-width of gapped excitations in 1D systems is known to show activation behavior 26 . With k B T glyph[lessmuch] ∆, in the present experiment these excitations are long lived and their apparent width corresponds to the experimental resolution. Cuts through both data sets at the band extrema are shown in Fig. 3. These correspond to the data of Fig. 2 integrated in thin slices of q ‖ / 2 π = -0 5 . ± 0 02 and . q ‖ / π 2 = 0 ± 0 03, respectively . (the dispersion is less curved at the maximum allowing to integrate a wider slice for improved counting statistics). To these cuts we fit two Gaussian peaks of equal width and a linear background shown as a solid lines in Fig. 3. We find a splitting of ∆ min = 50(1) µ eV at the band minimum and ∆ max = 40(2) µ eV at the band maximum.

## III. TRIPLET DISPERSION FOR A SPIN LADDER WITH EXCHANGE ANISOTROPY

## A. Isotropic Heisenberg ladder

For an isotropic spin S = 1 2 ladder we denote the rung and leg Heisenberg exchange constants as / J ⊥ and J ‖ , respectively. In the limit of λ = J /J ‖ ⊥ glyph[lessmuch] 1, the spin ladder can be understood as a 1D array of weakly interacting spin dimers. For λ = 0 the groundstate corresponds to a product state of singlets and the low energy excitations are triplet excitations. In the presence of weak leg coupling, these excitations become mobile and for small λ their dispersion can be calculated in a strong coupling expansion 27 . For an isotropic ladder with Hamiltonian

$$\mathcal { H } = J _ { \perp } \left ( \mathcal { H } ^ { \perp } + \lambda \mathcal { H } ^ { \| } \right ), \quad \mathcal { H } ^ { \perp } = \sum _ { R = 1 } ^ { L } \mathbf s _ { R, 1 } \cdot \mathbf s _ { R, 2 }, \quad \mathcal { H } ^ { \| } = \sum _ { R = 1 } ^ { L } \sum _ { \alpha = 1 } ^ { 2 } \mathbf s _ { R, \alpha } \cdot \mathbf s _ { R + 1, \alpha },$$

up to 3 rd order in λ , the following dispersion is obtained in Ref. 27, degenerate for the three σ = { + 0 , , -} triplet branches:

$$\epsilon _ { \sigma } ( k ) / J _ { \perp } = 1 + \gamma \cos ( k ) + \frac { \gamma ^ { 2 } } { 4 } \left [ 3 - \cos ( 2 k ) \right ] + \frac { \gamma ^ { 3 } } { 8 } \left [ 3 - 2 \cos ( k ) - 2 \cos ( 2 k ) + \cos ( 3 k ) \right ].$$

In the following we do the same calculation for a ladder with anisotropic rung or leg exchange, respectively. Some details of these computations are given in the appendix.

## B. Anisotropic rung interactions

In addition to the isotropic rung and leg couplings, here Ising-type anisotropy on the ladder rungs is considered. We choose the z axis parallel the Ising axis and consider the Hamiltonian

$$\mathcal { H } = J _ { \perp } \left ( \mathcal { H } ^ { \perp } + \lambda \mathcal { H } ^ { \prime } \right ), \quad \mathcal { H } ^ { \prime } = \mathcal { H } ^ { \| } + C \, \mathcal { H } ^ { \text{Rung,Using} }, \quad \mathcal { H } ^ { \text{Rung,Using} } = \sum _ { R = 1 } ^ { L } S _ { R, 1 } ^ { z } S _ { R, 2 } ^ { z },$$

where the Ising anisotropy on the ladder rungs is parameterized by C . Using the same strong coupling expansion as applied to the isotropic case in Ref. 27, up to 3 rd order in λ we find

$$\epsilon _ { 0 } ( k ) / J _ { \perp } = 1 + \lambda \cos ( k ) + \frac { \lambda ^ { 2 } } { 4 } \left [ 3 - \cos ( 2 k ) \right ] + \frac { \lambda ^ { 3 } } { 8 } \left [ 3 - 2 C - 2 \cos ( k ) - 2 \cos ( 2 k ) + \cos ( 3 k ) \right ] \quad \ ( 4 )$$

$$\epsilon _ { \pm } ( k ) / J _ { \perp } = 1 + \lambda \left [ \frac { C } { 2 } + \cos ( k ) \right ] + \frac { \lambda ^ { 2 } } { 4 } \left [ 3 - \cos ( 2 k ) \right ] + \frac { \lambda ^ { 3 } } { 8 } \left [ 3 - 2 C - 2 \cos ( k ) + \left ( C - 2 \right ) \cos ( 2 k ) + \cos ( 3 k ) \right ].$$

## C. Anisotropic leg interactions

Finally, we consider anisotropic exchange on the ladder legs with both symmetric and anti-symmetric contributions. Assuming a center of inversion symmetry on the ladder rungs, we consider a uniform DM vector + D on one leg and -D on the other leg. Furthermore we include Ising anisotropy of independent magnitude but pointing in the same direction as D . We choose the z axis parallel to D and consider

$$\mathcal { H } = J _ { \perp } \left ( \mathcal { H } ^ { \perp } + \lambda \mathcal { H } ^ { \prime } \right ), \quad \mathcal { H } ^ { \prime } = \mathcal { H } ^ { \| } + B \, \mathcal { H } ^ { L e g, \text{DM} } + A \, \mathcal { H } ^ { L e g, \text{Ising} }$$

$$\mathcal { H } ^ { L e g, D M } = \sum _ { R = 1 } ^ { L } \sum _ { \alpha = 1 } ^ { 2 } ( - 1 ) ^ { \alpha } \mathbf e _ { z } \cdot ( \mathbf S _ { R, \alpha } \times \mathbf S _ { R + 1, \alpha } ), \quad \mathcal { H } ^ { L e g, \text{ising} } = \sum _ { R = 1 } ^ { L } \sum _ { \alpha = 1 } ^ { 2 } S _ { R, \alpha } ^ { z } S _ { R + 1, \alpha } ^ { z }$$

Up to 3 rd order in λ we find the following triplet dispersion:

$$2$$

$$\text{Up to 3} ^ { r d } & \text{order in } \lambda \, \text{we find the following triplet dispersion:} \\ \epsilon _ { 0 } ( k ) / J _ { \perp } = & 1 + \lambda ( 1 + A ) \cos ( k ) + \frac { \lambda ^ { 2 } } { 4 } \left [ 3 - 4 B ^ { 2 } + 2 A + A ^ { 2 } + 4 B ^ { 2 } \cos ( k ) - \left ( 1 + 2 A + A ^ { 2 } \right ) \cos ( 2 k ) \right ] \\ & \quad + \frac { \lambda ^ { 3 } } { 8 } \left [ 3 + 3 A - \left ( 2 + 4 B ^ { 2 } + 4 A + 4 A B ^ { 2 } + 3 A ^ { 2 } + A ^ { 3 } \right ) \cos ( k ) \\ & \quad - \left ( 2 + 2 A - 4 B ^ { 2 } - 4 A B ^ { 2 } \right ) \cos ( 2 k ) + \left ( 1 + 3 A + 3 A ^ { 2 } + A ^ { 3 } \right ) \cos ( 3 k ) \right ] \\ \epsilon _ { \pm } ( k ) / J _ { \perp } = & 1 + \lambda \cos ( k ) + \frac { \lambda ^ { 2 } } { 4 } \left [ 3 - 2 B ^ { 2 } + 2 A + A ^ { 2 } - \cos ( 2 k ) \right ] \\ & \quad + \frac { \lambda ^ { 3 } } { 8 } \left [ 3 + 3 A - \left ( 2 + 1 2 B ^ { 2 } + A + \frac { A ^ { 2 } } { 2 } \right ) \cos ( k ) - \left ( 2 + 2 B ^ { 2 } + 2 A + 2 A B ^ { 2 } \right ) \cos ( 2 k ) + \frac { 1 } { 8 } \cos ( 3 k ) \right ] \\ & \quad \coprod$$

## IV. DISCUSSION

## A. Discussion of experimental results

The experimental data shown in Fig. 2 closely resembles previously published results on BPCB with dispersive triplet excitations 16,19 . In these experiments the resolution was not sufficient to detect the small band splitting and those data were very well described using a model of an isotropic spin ladder with J ⊥ /k B = 12 7(1) K and . J /k ‖ B = 3 54(3) K . 21 . Only with the superior resolution of the present experiment, the small splitting of the triplet band is observed.

teractions present in all materials 28 . For classical magnetic moments at the shortest Cu-Cu distance found in BPCB, dipolar interactions are estimated as J dipolar /k B ≈ 6 mK ( &lt; 1 µ eV), much smaller than the observed band splitting. Dipolar interactions clearly can be ruled out as the primary cause of magnetic anisotropy in BPCB. Since for the S = 1 2 Cu / 2+ ions there cannot be any single-ion anisotropy, exchange anisotropy must be the dominant cause of the observed small splitting of the triplet band. Before further analyzing our data in this regard, we first discuss the nature of anisotropic super-exchange interactions in the following section.

In our data, perpendicular to the one-dimensional ladders, no dispersion is observed within the experimental resolution. This validates BPCB as an exceptionally onedimensional ladder compound and it corroborates, that the observed band splitting is inherent to the spin ladders and is not related to residual 3D interactions. In addition, inter-ladder interactions were previously estimated on the order of a few 10s of mK 18 (a few µ eV) much smaller than the observed band splitting.

A previous ESR study on BPCB 20 found clear signatures of magnetic anisotropy, of exactly the same magnitude as the band splitting observed in the present experiments. By a careful study of the angular dependence of the ESR signal, an anisotropy axis tilted approx. 50 ◦ from the b axis in the ( b, c glyph[star] ) plane could be identified as a special direction. However, the origin of this anisotropy remained unclear.

One source of magnetic anisotropy are dipolar in-

## B. Anisotropic super-exchange

Most generally, the effective exchange interaction between two spins takes the form H 1 2 , = S 1 Γ S 2 where Γ is a tensor. Decomposing Γ into its symmetric and anti-symmetric parts and choosing a basis such that the symmetric part is diagonal, this can be written as

$$\underset { \underset { \underset { \text{min}} { \sigma } } { \text{min} } } { \text{max} } \quad \mathcal { H } _ { 1, 2 } = J \, \mathbf S _ { 1 } \cdot \mathbf S _ { 2 } + \mathbf D \cdot ( \mathbf S _ { 1 } \times \mathbf S _ { 2 } ) + \sum _ { \alpha = x, y, z } G ^ { \alpha } \, S _ { 1 } ^ { \alpha } S _ { 2 } ^ { \alpha }. \ ( 1 0 )$$

Here, J is the isotropic (Heisenberg) exchange. The so-called Dzyaloshinskii-Moriya vector D quantifies the anti-symmetric exchange and G α ( α = x, y, z ) determine the symmetric (Ising) exchange anisotropy where ∑ α G α = 0.

Including spin-orbit coupling (SOC) into Anderson's theory of super-exchange 1 indeed all these contributions

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00724v2-figure-c907c8b9-12e2-4240-84f1-62d09a48b48d.jpg)

FIG. 4. Triplet dispersion for a spin ladder with exchange anisotropy: a) Ising anisotropy on the ladder rungs, b) Ising anisotropy on the legs, c) DM anisotropy on the legs, and d) both Ising and DM anisotropy on the ladder legs, in the ratio as predicted for a simple super-exchange mechanism (Eq. 11). For these plots, isotropic J ⊥ /k B = 12 7 K and . J /k ‖ B = 3 54 K . were used, approximately describing the ladders in BPCB. Further, in all four cases the anisotropy parameter was chosen such that the triplet splitting at the band minimum amounts to 50 µ eV. In all plots, the blue line shows the doubly degenerate glyph[epsilon1] ± ( k ) dispersion and the orange line shows the glyph[epsilon1] 0 ( k ) triplet.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00724v2-figure-313b61ad-2534-4666-a1ee-f6f5c89544e2.jpg)

are obtained 2,3 . However, in this setting, symmetric and anti-symmetric interactions are not independent and for a single bond, to a good approximation the superexchange interaction Hamiltonian reads 4-6

## C. Exchange anisotropy in BPCB

For the compound BPCB with centers of inversion symmetry on the ladder rungs, we consider the following cases compatible with the crystallographic symmetry:

- a) Ising anisotropy on the ladder rungs
- b) Ising anisotropy on the ladder legs
- c) DM anisotropy on the ladder legs
- d) Both DM and Ising anisotropy on the ladder legs as predicted for simple super-exchange (Eq. 11).

Using the results of Sec. III, for these four cases, the triplet dispersion is plotted in Fig. 4. For these plots, we have used the Heisenberg exchange constants J ⊥ /k B = 12 7 . K and J /k ‖ B = 3 54 K approximately describing . BPCB 21 . In all cases the anisotropy parameter was chosen such that the triplet splitting at the band minimum amounts to 50 µ eV.

Comparing the dispersion calculated for Ising-type exchange anisotropy on the ladder rungs (Fig. 4a) to the data of Fig. 2, we can clearly rule out this case as the dominant source of anisotropy in BPCB. The data qualitatively disagree with the calculated dispersion. For the case of anisotropic leg exchange, all three types of anisotropy considered give qualitatively the same dispersion (Fig. 4b,c,d). There is a doubly degenerate band with a smaller bandwidth and a non-degenerate band with larger bandwidth. In all cases the calculated anisotropy splitting is slightly larger at the band minimum than at the band maximum, just as in our data 31 .

In Fig. 5 we show the triplet dispersion extracted from our data. These points were obtained from Gaussian fits to constantq ‖ -cuts, similar to the ones shown in Fig. 3. The vertical bars denote the width of the observed peaks. Fitting the calculated dispersions with anisotropic leg exchange (cases b,c,d) to this data, we find excellent agreement for all three cases. Indeed, the three calculated dispersions are so similar, that from our data, it is impossible to determine which type of exchange anisotropy on the ladder legs actually causes the observed band splitting in BPCB. As an example, in

$$\mathcal { H } _ { 1, 2 } & = \left ( J - \frac { | \mathbf D | ^ { 2 } } { 4 J } \right ) \mathbf S _ { 1 } \mathbf S _ { 2 } + \mathbf D \cdot ( \mathbf S _ { 1 } \times \mathbf S _ { 2 } ) + \frac { 1 } { 2 J } ( \mathbf D \cdot \mathbf S _ { 1 } ) ( \mathbf D \cdot \mathbf S _ { 2 } ). \quad \ \ \ c ) \ I \\ \mathbf U \dots \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d } { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n D } { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { c e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { [ \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { a n d e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ a n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ o n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ ele n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ s e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { } e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d {  e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { $ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { n e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { i e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e n d { \ e }$$

Here a single vector D defines both the symmetric and antisymmetric exchange. This expression, although not obvious, is fully invariant under spin rotations 5,6 : While the term D · ( S 1 × S 2 ) acts to confine the spins to a plane perpendicular to the vector D , the Ising contribution acts to align the spins with D . Thus, for a single bond, anisotropic super-exchange will never single out a particular direction. Nonetheless, for multiple connected bonds, frustration of the D vectors may still beak spin rotation symmetry and lead to the appearance of anisotropy 5,6 .

Besides this mechanism, considering multiple orbitals of the ligands transmitting the super-exchange interactions, in the presence of SOC may result in additional sources of exchange anisotropy, beyond Eq. 11. However, these additional effects have been argued to be much smaller in magnitude 29,30 .

This leaves us with the following situation: When considering super-exchange interactions, to first approximation symmetric and antisymmetric anisotropy contributions always come together pointing in the same direction and at a fixed ratio. In this case, if the center of a bond corresponds to a center of inversion symmetry, not only is antisymmetric exchange anisotropy forbidden by symmetry as is commonly known 23 , but we also expect symmetric exchange anisotropy to vanish. Only additional mechanisms may then lead to Ising-type anisotropy which is in principal allowed on a bond with inversion symmetry.

FIG. 5. Position of the triplet bands extracted from the neutron scattering data of Fig. 2. The dispersion drawn as a solid line represents a model with anisotropic interactions on the ladder legs (case d) as described in the text.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00724v2-figure-cb893f84-5785-47b1-9813-0945f4e40158.jpg)

Fig. 5 the solid line shows the dispersion calculated for anisotropic exchange on the ladder legs given by Eq. 11 (case d) with fitted parameters

$$J _ { \perp } / k _ { B } = & 1 2. 7 7 ( 1 ) \, K, \quad ( 1 2 ) \quad \text{$\text{$type$} \\ J _ { \| } / k _ { B } = 3. 5 5 ( 1 ) \, K, \ \ D _ { \| } / k _ { B } = 1. 4 4 ( 2 ) \, K. \quad \text{$\text{$of $\substack{$us$} \\ \text{We $a$} }$$

In the parametrization of the Hamiltonian of Sec. III C, the coefficients are ˜ J ‖ = J ‖ -D / 2 ‖ (4 J ‖ ), λ = ˜ J /J ‖ ⊥ , A = D / 2 ‖ (2 J J ‖ ˜ ) and ‖ B = D /J ‖ ˜ . ‖ Here, we stress that the numerical value of D /k ‖ B = 1 44 K cannot serve as . an estimate for the overall magnitude of the magnetic anisotropy which should be estimated from the band splitting of 50 µ eV, i.e. ≈ 0 6 K in units of Kelvin. .

This leaves the question, how one might experimentally determine the type of leg exchange anisotropy dominant in BPCB. The most promising route to answering this question would be to measure the triplet dispersion as a function of orientation and strength of a small applied magnetic field. In combination with calculations also including an applied magnetic field, such data might indeed allow to determine the precise type of exchange anisotropy on the ladder legs in BPCB. However, this is clearly beyond scope of the present work.

## D. Final remarks

In BPCB the interactions on the ladder rungs are the strongest ones by far ( J ⊥ ≈ 3 6 . J ‖ ). On these bonds, antisymmetric exchange anisotropy is prohibited by symmetry and only Ising-type exchange anisotropy is allowed. It may thus seem plausible to expect this to be the dominant source of anisotropy. Yet, experimentally we find that the anisotropy is predominantly due to the much weaker leg interactions of the ladder. Indeed, this is fully consistent with the theory of superexchange interactions, where symmetric and antisymmetric anisotropy contributions are always coupled. If either is prohibited by symmetry, also the other will vanish. The negligible Ising anisotropy on the ladder rungs in BPCB we thus see as supporting evidence, for the theory of anisotropic super-exchange 5,6 .

Since the theory of super-exchange predicts DM anisotropy B ∝ | D | and Ising anisotropy A ∝ | D | 2 , initially the latter was considered subdominant 2,3 . However, both contributions are now understood to be equally important 4-6 . Whilst not immediately intuitive, the triplet dispersions calculated in Sec. III C illustrate this effect: In the obtained expressions A appears linearly, while B only appears to second order as B 2 .

Finally, we mention an apparent contradiction: Superexchange can never give rise to DM anisotropy exclusively. Yet, oftentimes pure DM interactions are used very successfully to explain anisotropy effects observed experimentally. Here, the present case study offers some illustration. For the spin ladder, different types of exchange anisotropy on the ladder legs lead to almost identical triplet dispersions that all explain our data equally well. We speculate that also in other systems one might encounter similar situations and it is for this reason that pure DM anisotropy is used so successfully to explain experimental findings, even though (within the theory of super-exchange) it lacks any microscopic justification. We are only aware of one experimental study on the helimagnet Ba CuGe O , where a model employing the full 2 2 7 Hamiltonian for anisotropic super-exchange interactions (Eq. 11) was compared to experimental data and indeed better agreement was found than for a model employing DM anisotropy only 32 .

Further case studies of the nature of magnetic exchange anisotropy would certainly be interesting. We also suggest, that whenever a model with pure DM interactions is considered, it would be enlightening to also consider a model with both Ising and DM anisotropy, compatible with super-exchange.

## V. CONCLUSION

The prototypical spin ladder compound (C H 5 12 N) CuBr 2 4 (BPCB) has been studied by means of high resolution inelastic neutron scattering. We find a small splitting of the triplet band of 50(1) µ eV at the band minimum and 40(2) µ eV at the band maximum. Further, for a spin ladder with exchange anisotropy, the triplet dispersion is calculated in a strong coupling expansion.

The Ising-type anisotropy allowed by crystallographic symmetry on the ladder rungs (by far the strongest bonds) we find to be negligible, in line with the theory of anisotropic super-exchange. Three models with exchange anisotropy on the ladder legs, all describe the data equally well. Whilst we cannot distinguish these by

comparison to our data, we note that only one is compatible with super-exchange interactions.

Acknowledgement. This work is partially supported by the Swiss National Science Foundation under Division II. Experiments at the ISIS Neutron and Muon Source were supported by a beamtime allocation RB1720009 from the Science and Technology Facilities Council, UK.

## Appendix: Strong coupling expansion

In this appendix we give some details on the strong coupling expansion calculation of the triplet dispersion in a strong rung spin ladder following Ref. 27: First, consider a single Heisenberg S = 1 2 spin dimer with / Hamiltonian H 0 = J ⊥ S 1 · S 2 . The eigenstates are given by a singlet state | s 〉 with energy E s = -3 4 J ⊥ and three triplet states | t σ 〉 with energy E t = 1 4 J ⊥ and S z spin component σ = -1 0 1. , , In all cases our calculations start from the strong rung limit, where the spin ladder is nothing but an array of L independent Heisenberg spin dimers

$$\mathcal { H } ^ { \perp } = J _ { \perp } \sum _ { R = 1 } ^ { L } \mathbf S _ { R, 1 } \cdot \mathbf S _ { R, 2 }.$$

The ground state is a product of singlets | 0 〉 = | s . . . s 〉 with energy E 0 = LE s . The first excited states are the 3 L degenerate single triplet states | R,σ 〉 = | s . . . t σ R . . . s 〉 where the R th rung is excited to an S z = σ triplet.

Now we introduce the much weaker interaction H ′ as

∗ dblosser@phys.ethz.ch

- † zhelud@ethz.ch; http://www.neutron.ethz.ch/
- 1 P. W. Anderson, Phys. Rev. 115 , 2 (1959).
- 2 T. Moriya, Phys. Rev. Lett. 4 , 228 (1960).
- 3 T. Moriya, Phys. Rev. 120 , 91 (1960).
- 4 T. A. Kaplan, Zeitschrift f¨r Physik B Condensed Matter u 49 , 313 (1983).
- 5 L. Shekhtman, O. Entin-Wohlman, and A. Aharony, Phys. Rev. Lett. 69 , 836 (1992).
- 6 L. Shekhtman, A. Aharony, and O. Entin-Wohlman, Phys. Rev. B 47 , 174 (1993).
- 7 M. Oshikawa and I. Affleck, Phys. Rev. Lett. 79 , 2883 (1997).
- 8 I. Affleck and M. Oshikawa, Phys. Rev. B 60 , 1038 (1999).
- 9 S. A. Zvyagin, A. K. Kolezhuk, J. Krzystek, and R. Feyerherm, Phys. Rev. Lett. 93 , 027201 (2004).
- 10 R. Coldea, D. A. Tennant, K. Habicht, P. Smeibidl, C. Wolters, and Z. Tylczynski, Phys. Rev. Lett. 88 , 137203 (2002).
- 11 M. H¨ alg, W. E. A. Lorenz, K. Y. Povarov, M. M˚nsson, a Y. Skourski, and A. Zheludev, Phys. Rev. B 90 , 174413 (2014).
- 12 M. Kenzelmann, A. Zheludev, S. Raymond, E. Ressouche, T. Masuda, P. B¨ oni, K. Kakurai, I. Tsukada, K. Uchinokura, and R. Coldea, Phys. Rev. B 64 , 054422 (2001).
- 13 H. Jacobsen, S. M. Gaw, A. J. Princep, E. Hamilton,

a perturbation. It connects the individual dimers giving the spin ladder Hamiltonian H = J ⊥ ( H ⊥ + λ H ′ ) with λ glyph[lessmuch] 1. In section III A, H ′ describes Heisenberg exchange connecting the dimers via the ladder legs. In addition to this, in Secs. III B and III C, H ′ also contains anisotropic exchange contributions on the ladder rungs or legs, respectively.

The ground state energy in the presence of H ′ we denote as ˜ . E 0 It is computed using standard nondegenerate Rayleigh-Schr¨dinger o perturbation theory. The triplet states however are 3 L -fold degenerate and some care is required. The L fold degeneracy is lifted to first order in λ and the appropriate eigenstates are Bloch waves

$$| k, \sigma \rangle = \frac { 1 } { \sqrt { L } } \sum _ { R = 1 } ^ { L } e ^ { i k R } \, | R, \sigma \rangle$$

for k = n 2 π , n = 0 , . . . , L -1. In the absence of H ′ they have energy E k,σ = E 0 + J ⊥ . Further we note that all H ′ commute with S z = ∑ L R =1 S z R and thus with the unperturbed H ⊥ . Therefore, they will never mix triplets with different quantum numbers σ . Starting with the | k, σ 〉 triplet Bloch waves we can therefore again use nondegenerate perturbation theory to obtain their energy ˜ E k,σ in the presence of the perturbation H ′ . Calculating the ground state energy ˜ E 0 and the triplet energy ˜ E k,σ to third order perturbation theory in λ , we obtained the triplet dispersions glyph[epsilon1] σ ( k ) = ˜ E k,σ -˜ E 0 . For the different H ′ , these are given in Sec. III.

- S. T´ oth, R. A. Ewings, M. Enderle, E. M. H. Wheeler, D. Prabhakaran, and A. T. Boothroyd, Phys. Rev. B 97 , 144401 (2018).
- 14 S. Bettler, F. Landolt, O. M. Aksoy, Z. Yan, S. Gvasaliya, Y. Qiu, E. Ressouche, K. Beauvois, S. Raymond, A. N. Ponomaryov, S. A. Zvyagin, and A. Zheludev, Phys. Rev. B 99 , 184437 (2019).
- 15 B. R. Patyal, B. L. Scott, and R. D. Willett, Phys. Rev. B 41 , 1657 (1990).
- 16 A. T. Savici, G. E. Granroth, C. L. Broholm, D. M. Pajerowski, C. M. Brown, D. R. Talham, M. W. Meisel, K. P. Schmidt, G. S. Uhrig, and S. E. Nagler, Phys. Rev. B 80 , 094411 (2009).
- 17 C. R¨ uegg, K. Kiefer, B. Thielemann, D. F. McMorrow, V. Zapf, B. Normand, M. B. Zvonarev, P. Bouillot, C. Kollath, T. Giamarchi, S. Capponi, D. Poilblanc, D. Biner, and K. W. Kr¨ amer, Phys. Rev. Lett. 101 , 247202 (2008).
- 18 M. Klanjˇ sek, H. Mayaffre, C. Berthier, M. Horvati´, c B. Chiari, O. Piovesana, P. Bouillot, C. Kollath, E. Orignac, R. Citro, and T. Giamarchi, Phys. Rev. Lett. 101 , 137207 (2008).
- 19 B. Thielemann, C. R¨egg, H. M. Rønnow, A. M. L¨uchli, u a J.-S. Caux, B. Normand, D. Biner, K. W. Kr¨mer, H.-U. a G¨ udel, J. Stahn, K. Habicht, K. Kiefer, M. Boehm, D. F. McMorrow, and J. Mesot, Phys. Rev. Lett. 102 , 107204 (2009).

- 20 E. ˇ iˇm´r, M. Ozerov, J. Wosnitza, B. Thielemann, K. W. C z a Kr¨ amer, C. R¨ uegg, O. Piovesana, M. Klanjˇ sek, M. Horvati´, C. Berthier, c and S. A. Zvyagin, Phys. Rev. B 82 , 054431 (2010).
- 21 D. Blosser, V. K. Bhartiya, D. J. Voneshen, and A. Zheludev, Phys. Rev. Lett. 121 , 247201 (2018).
- 22 The crystal structure plots shown in Fig. 1 are based on renderings generated by the VESTA software package 33 .
- 23 I. Dzyaloshinsky, Journal of Physics and Chemistry of Solids 4 , 241 (1958).
- 24 R. Bewley, J. Taylor, and S. Bennington., Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 637 , 128 (2011).
- 25 The raw data is available in the STFC Research Data Repository 34 .
- 26 A. Zheludev, V. O. Garlea, L.-P. Regnault, H. Manaka, A. Tsvelik, and J.-H. Chung, Phys. Rev. Lett. 100 , 157204 (2008).
- 27 M. Reigrotzki, H. Tsunetsugu, and T. M. Rice, Journal of Physics: Condensed Matter 6 , 9235 (1994).
- 28 P. Fazekas, Lecture Notes on Electron Correlation and Magnetism (WORLD SCIENTIFIC, 1999)
- https://www.worldscientific.com/doi/pdf/10.1142/2945.
- 29 W. Koshibae, Y. Ohta, and S. Maekawa, Phys. Rev. Lett. 71 , 467 (1993).
- 30 L. Shekhtman, O. Entin-Wohlman, and A. Aharony, Phys. Rev. Lett. 71 , 468 (1993).
- 31 In Fig. 4 the anisotropy parameter was chosen so as to ensure ∆ min = 50 µ eV. For the different models of leg anisotropy (Ising only, DM only, both), this gives ∆ max = 37 32 and 35 , µ eV. However, since from 2 nd to 3 rd order perturbation expansion the dispersion still slightly changed, these small differences are probably not reliable enough to use as a criterion to select one of these three cases as the most suitable model describing the compound BPCB.
- 32 A. Zheludev, S. Maslov, I. Tsukada, I. Zaliznyak, L. P. Regnault, T. Masuda, K. Uchinokura, R. Erwin, and G. Shirane, Phys. Rev. Lett. 81 , 5410 (1998).
- 33 K. Momma and F. Izumi, Journal of Applied Crystallography 44 , 1272 (2011).
- 34 D. Blosser, V. K. Bhartiya, and A. Zheludev, RB1720009, STFC ISIS Facility (2018), 10.5286/ISIS.E.87813753.
D. Blosser V. K. Bhartiya D. J. Voneshen A. Zheludev
Origin of magnetic anisotropy in the spin ladder compound (C$_5$H$_{12}$N)$_2$CuBr$_4$
The $S=1/2$ spin ladder compound (C$_5$H$_{12}$N)$_2$CuBr$_4$ (BPCB) is
studied by means of high-resolution inelastic neutron scattering. In agreement
with previous studies we find a band of triplet excitations with a spin gap of
$\sim0.8$~meV and a bandwidth of $\sim0.6$~meV. In addition, we observe a
distinct splitting of the triplet band of $50(1)$~$\mu$eV or $40(2)$~$\mu$eV at
the band minimum or maximum, respectively. By comparison to a strong coupling
expansion calculation of the triplet dispersion for a spin ladder with
anisotropic exchange, weakly anisotropic leg interactions are identified as the
dominant source of magnetic anisotropy in BPCB. Based on these results, we
discuss the nature of magnetic exchange anisotropy in BPCB and in related
transition-metal insulators.