## Resonant inelastic x-ray scattering study of α -RuCl : a progress report 3

## Blair W Lebert , Subin Kim , Valentina Bisogni , Ignace 1 1 2 Jarrige , Andi M Barbour , Young-June Kim 2 2 1

- 1 Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 1A7
- 2 National Synchrotron Light Source II, Brookhaven National Laboratory, Upton, NY, 11973, USA

E-mail:

yjkim@physics.utoronto.ca

July 11th, 2017

Abstract. Ru M -edge resonant 3 inelastic x-ray scattering (RIXS) measurements of α -RuCl 3 with 27 meV resolution reveals a spin-orbit exciton without noticeable splitting. We extract values for the spin-orbit coupling constant ( λ = 154 ± 2 meV) and trigonal distortion field energy ( ∆ | | &lt; 65 meV) which support the j eff = 1 2 nature of / α -RuCl . 3 We demonstrate the feasibility of Medge RIXS for 4 d systems, which allows ultra high-resolution RIXS of 4 d systems until instrumentation for L-edge RIXS improves.

Keywords : Kitaev quantum spin liquid, resonant inelastic x-ray scattering, honeycomb lattice, spin-orbit coupling, crystal field

Submitted to: J. Phys.: Condens. Matter

## 1. Introduction

Study of Kitaev materials has blossomed in the last decade as a promising direction to find quantum spin liquids [1]. The Kitaev honeycomb model is an exactly solvable model with quantum spin liquid (QSL) ground states [1-4]. One of the key developments in this field was the proposal by Jackeli and Khaliullin that the crucial ingredient of the Kitaev model, the bonddependent Kitaev interaction, can be realized in a magnetic insulator with a strong spin-orbit coupling (SOC), such as the iridates where Ir 4+ ions take on a low-spin d 5 state due to large octahedral crystal field splitting [5]. In particular, Chaloupka, Jackeli, and Khaliullin pointed out that a honeycomb lattice formed by edge-shared IrO 6 octahedra could be a realization of Kitaev's honeycomb model [6]. The successful synthesis of a honeycomb iridate compound Na IrO 2 3 [7, 8] and other iridate materials based on a honeycomb lattice or its three-dimensional variants opened a new genre of Kitaev materials research [9, 10], which has become an important branch of QSL research. For a comprehensive review on Kitaev materials, see Refs. [3, 11].

One of the most promising Kitaev materials is α -RuCl , which has been drawing much attention ever 3 since it was first suggested as a candidate material for Kitaev QSL by Plumb et al [12]. Subsequent observation of a broad excitation continuum by inelastic neutron scattering [13, 14] and Raman scattering [15] suggests fractionalization of spin excitations expected for a QSL. Recently, the observation of quantized thermal Hall conductivity in the field-induced state of α -RuCl 3 has caused much excitement in the field of QSLs [16].

The fact that α -RuCl 3 has emerged as a prime candidate for Kitaev QSL is somewhat surprising since its SOC is smaller than that of the iridates. However, one should realize that what is important is not necessarily the bare value of SOC, but the size of SOC in comparison to a non-cubic (i.e. trigonal or tetragonal) crystal field energy scale caused by the small distortions of RuCl 6 or IrO 6 octahedra. A large enough trigonal and/or tetragonal distortion could quench orbital angular momentum and spoil the Kitaev interaction between the j eff = 1 2 pseudospins. / Earlier structural report for α -RuCl 3 indicated that the RuCl 6 octahedron is free of trigonal distortion [17], which was one of the reasons behind the original proposal by Plumb et al [12]. However, recent powder neutron diffraction [18] and single crystal xray diffraction [19] studies revealed that the crystal structure is quite complicated due to stacking disorder, and most importantly, there exists small trigonal distortion. Therefore, it is in fact very important to re-examine whether the j eff = 1 2 description survives / this trigonal crystal field effect, which is a necessary condition for realizing Kitaev physics in α -RuCl . 3

Resonant inelastic x-ray scattering (RIXS) has emerged as an invaluable technique for measuring electronic and magnetic excitations. In particular, RIXS is an exquisite tool for studying transitions between crystal-field-split d orbitals. The electronic transitions between the d orbitals, or dd transitions, are dipole-forbidden and therefore typically difficult to study with optical spectroscopy. However, RIXS is a second-order process in which two dipole-allowed transitions are combined to provide a highly sensitive probe of dd transitions. For example, Gretarsson et al measured the splitting of dd transitions in Na 2 IrO 3 and showed that the splitting due to the trigonal crystal field (0.1 eV) is much smaller than the overall splitting between the j eff = 1 2 and / j eff = 3 2 states, confirming / the j eff = 1 2 nature of its magnetism [20]. /

In this paper, we report our RIXS investigation of the dd transitions in α -RuCl . Similar to the Na IrO 3 2 3 case, our RIXS results show that the trigonal splitting is too small to observe with the current instrumental resolution. Our results thus suggest that the trigonal splitting is much smaller than the SOC energy scale, confirming the j eff = 1 2 nature of this / compound. Our result also indicates that the g factor in α -RuCl 3 should be close to isotropic limit g ⊥ /g ‖ ≈ 1, in agreement with a recent x-ray absorption study [21].

We would like to also note that our measurements were carried out at the Ru M 3 edge, not at the L 3 edge which is more widely used for RIXS studies of cuprates and iridates. This study illustrates that using the M 3 edge instead of the L 3 edge is an option to be explored for the study of 4 d elements, which allows researchers to use the existing soft x-ray RIXS beamlines instead of tender/intermediate energy range RIXS beamlines. We will discuss the advantages and disadvantages of using the M 3 edge instead of the L 3 edge.

## 2. Experimental

Single crystals of α -RuCl 3 were prepared by vacuum sublimation of commerical RuCl 3 powder (SigmaAldrich, Ru content 45% - 55%) in sealed quartz tubes, as in Ref. [12]. The three-zone tube furnace was cooled from 700 C to 400 C at 0.8 ◦ ◦ ◦ C/hour while maintaining a 10 ◦ C difference between the tube ends. The resulting platelike crystals have their crystallographic c axis perpendicular to the surface and display hexagonal facets. Unless otherwise noted, we use the hexagonal notation corresponding to the the trigonal P 3 12 space 1 group from Ref. [17] in this paper ( a = b = 5 96 ˚ , . A c = 17 2 ˚ , . A α = β = 90 , ◦ γ = 120 ). ◦ The magnetic and thermodynamic properties of these single crystals have been well-characterized [22].

Single crystals were aligned and characterized using a four-circle x-ray diffractometer with a Motube source. The sample we report on in this paper was relatively large, with a ≈ × 5 5 mm 2 area. Rocking curves taken at the (0,0,9) and (3,0,12) Bragg reflections showed a 1.5 ◦ mosaic. The crystal facets were confirmed to be high-symmetry hexagonal axes, which were used to align the sample during mounting. The sample was mounted on a copper holder using silver epoxy and cured for 3 hours at 80 C. ◦ X-ray diffraction of the sample on the holder confirmed that the sample was aligned with 〈 1 0 0 〉 and 〈 0 0 1 〉 in the horizontal scattering plane.

We performed RIXS at the Ru M 3 edge which is at ≈ 462 eV, i.e. in the soft x-ray regime. Therefore, we used the SIX beamline [23, 24] at the NSLSII. The VLS plane grating on the monochromator and spectrometer had line spacing of 500 mm -1 and 1250 mm -1 respectively, which coupled with a 20 µ m exit slit gave a combined resolution of 27 meV FWHM as measured on carbon tape. The x-ray beam size at the sample was 6 (H) × 13 (V) µ m . 2 A 2 θ = 90 ◦ scattering angle and π -polarized incident x-rays were used to reduce the elastic line. The scattering angle gave a consant momentum transfer of q = 2 k sin θ = 0 33 . ˚ A -1 . All measurements reported here were performed at 20 ◦ grazing incidence giving q ‖ = 0 14 ˚ . A -1 in-plane. In terms of reciprocal lattice units (r.l.u) we measured in-plane at q = (0 130), . shown as a black diamond in the Brillouin zone inset of figure 3. All measurements were performed at ambient temperature. The samples were cleaved with tape before putting the holder in the vacuum load-lock.

## 3. Results

In α -RuCl , the Ru 3 3+ ions are surrounded by a nearly ideal octahedral arrangement of Cl -ions (figure 1). These RuCl 6 octahedra share an edge to form a honeycomb net. The octahedral cubic crystal field,

CFS (∆ &gt; 0)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01190v1-figure-5ba243c8-f923-41f4-b715-bb3920d31715.jpg)

glyph[negationslash]

Figure 1. Schematic energy level diagram of Ru 4 d orbitals in α -RuCl 3 under the influence of crystal field splitting (CFS) and spin-orbit coupling (SOC). On the left we show an ideal RuCl 6 octahedron, i.e. with no trigonal distortion (∆ = 0). Under this octahedral crystal field the degeneracy of the 4 d energy levels is lifted and they split into t 2g and e g (not shown) manifolds. Turning on SOC splits the t 2g manifold in α -RuCl 3 into a lower j eff = 3 2 level and higher / j eff = 1 2 level separated / by 3 λ/ 2, where λ is the SOC constant. A trigonal distortion field (∆ = 0) will lead to further splitting. Here we show the case of trigonal compression (∆ &gt; 0) which is illustrated by the arrows compressing the RuCl 6 octahedron on the right. Under small trigonal compression the j eff = 3 2 levels split giving in the end / three non-degenerate doublets. We have labeled these levels as X Y , , and Z corresponding to A B , , and C respectively in the notation of Chaloupka et al [25] to avoid confusion with our labeling of the RIXS peaks. The splitting of the lower two levels defines D trig = E Y -E Z , which is therefore positive (negative) for trigonal compression (elongation).

10 Dq , splits the Ru 4 d levels into lower t 2g and upper e g manifolds. It is believed that SOC further splits the t 2g manifold into upper j eff = 1 2 and / lower j eff = 3 2 states. / The 4 d 5 electronic configurations of the Ru ions should take on a low-spin state, filling up the j eff = 3 2 level completely leaving a half-filled / j eff = 1 2 level / and leading to the so-called j eff = 1 / 2 pseudospin. The RuCl 6 octahedra however have a slight trigonal distortion which might quench the oribital angular momentum and destroy the j eff = 1 2 / nature of α -RuCl . 3

X-ray absorption spectroscopy (XAS) is a useful tool to shed light on the influence of SOC on electronic structure. XAS measurements by Plumb et al [12] on α -RuCl 3 at the Ru L 2 (2 p 1 / 2 ) and L 3 (2 p 3 / 2 ) edges observed a difference in their lineshape and found an anomalously large L /L 3 2 intensity ratio, the so-called branching ratio (BR). Assuming nonnegligible SOC, atomic dipole transitions must follow the J selection rules. Therefore at the L 2 edge, 2 p 1 / 2 → 4 d 3 / 2 is allowed and 2 p 1 / 2 → 4 d 5 / 2 is forbidden, while at the L 3 edge, 2 p 3 / 2 → 4 d 3 / 2 and 2 p 3 / 2 → 4 d 5 / 2 are both allowed. This is manifested in the XAS measurements of Plumb et al , where they find the L 2 absorption edge has a single peak, compared to double peaks at the L 3 edge. Of these two peaks, the lower (higher) energy one is due to t 2g ( e g ) states and correspondingly has a lower (higher) intensity because of its one (four) hole(s). The absence of a t 2g peak at the L 2 edge is because it

acquires a J = 5 2 character due to SOC. The BR is / one way of expressing this effect of SOC, where a BR &gt; 2 is expected with SOC effects due to a reduction of dipole-allowed transitions for the L 2 edge as explained above. BR = 2 is typical without SOC since the L 3 edge has twice as many electrons available as the L 2 edge: 2 p 3 / 2 is a quartet and 2 p 1 / 2 is a doublet. Plumb et al reported a BR = 3 ± 0 5 for . α -RuCl . 3

We performed XAS at the Ru M 2 (3 p 1 / 2 ) and M 3 (3 p 3 / 2 ) edge in total electron yield (TEY) mode, shown as a black line in figure 2(b). The first noticeable difference is that the t 2g -e g double peak structure is not visible at the M 3 edge. The M edges correspond to shallower core holes than the L edges, i.e. lower binding energy, and therefore have shorter lifetimes. The decreased lifetime leads to broadening in the energy domain which make it appear as a single peak. Nonetheless, the peak appears asymmetric and has a slight shoulder near 460.5 eV, which with the main peak at 462.5 eV would give an estimate of the splitting between t 2g and e g of ≈ 2 0 eV with a large error bar. . In comparison, a splitting of ≈ 2 3 eV was observed in . the L-edge XAS data [12]. We also roughly estimated the BR by subtracting a linear background at each peak and integrating the intensity. We found BR = 3 6 . ± 0 8 . which is consisent with the BR from Ledge XAS [12].

We first studied the resonant behavior of RIXS by scanning the incident energy E i as shown in figure 2(a). We measured RIXS spectra at five different energies, which are also shown relative to the TEY-XAS signal as vertical lines with the same color in figure 2(b). The Ru M -edge RIXS directly 3 probes the crystalfield split Ru 4 d levels, the so-called dd excitations. The two dipole transitions, 3 p 3 / 2 → 4 d followed by 4 d → 3 p 3 / 2 , allow dipole-forbidden transitions from occupied t 2g states to empty t 2g and e g states. The resonant enhancement is indeed strong, over an order of magnitude, and we find two different resonance regimes. The excitations to unoccupied t 2g states correspond to the spectral weight below &lt; 1 eV while excitations to e g states occur above this energy. We quantitatively studied this resonant behavior by integrating the spectral weight in each region. The integration ranges are shown as solid bars at the bottom of figure 2(a): the t 2g region in green from 0.1 to 0.4 eV and the e g region in magenta from 1.9 to 2.7 eV. The integrated intensities are plotted in arbitrary units in figure 2(b) as green squares for t 2g and magenta diamonds for e g (the five points are fit with a Lorentzian function and constant background as shown by the corresponding solid lines). The t 2g and e g states resonate at approximately 460.8 eV and 462.4 eV respectively. The ≈ 1 6 . eV difference between these energies corresponds roughly to the 10 Dq octahedral crystal-field energy, estimated to be ≈ 2 0 eV earlier. .

To search for splitting due to trigonal distortions we focused on the t 2g region. Therefore, we measured a high statistics spectrum with E i = 461 eV, i.e. at the resonance of the t 2g states. In figure 3, we show a high-resolution Ru M -edge RIXS spectrum on 3 α -RuCl 3 which has many features visible in the energy loss region: three low-energy features A , A , and A 1 2 3 (zoom in figure 4) at 231 ± 3 meV, 524 ± 10 meV, and 745 ± 10 meV; a charge gap around 1 eV; and, highenergy features α , β , B, γ , and δ above 1 eV. The features have been labeled with the same notation as Sandilands et al [26] with the addition of the B peak.

The energies of the dd excitations are consistent with those found with optical spectroscopy [26] and quantum chemistry calculations [27]. The quantum chemistry calculation results are indicated in figure 3. Yadav and coworkers calculated dd transition energies using multireference configuration-interaction (MRCI) calculations with SOC for two different structures: P 3 12 1 [17] shown in red and C /m 2 [19] shown in blue. Besides the B peak, the other high-energy RIXS features, α β , , and γ , are broad and difficult to compare with calculations. Nonetheless, the α β , , and B peaks correspond well with predicted energies for the C /m 2 structure, although the RIXS measurements are at a slightly higher energy. On the other hand, the P 3 12 structure predicts energies which are far too low 1 compared to the RIXS data. Recent structural studies seem to agree on the fact that the room temperature structure is described by the C /m 2 symmetry [19, 28, 29]. Our measurements performed at 300 K clearly agrees with the C /m 2 structure prediction.

The low-energy features of the RIXS spectrum are shown in figure 4(a). The A 1 peak is a dd excitation from j eff = 3 2 to / j eff = 1 2, and is therefore / also known as a spin-orbit exciton [30]. Ignoring trigonal distortions, the A 1 peak position corresponds to the energy difference between these levels, 231 ± 3 meV. The splitting between these levels is equal to 3 λ/ 2 (figure 1) from which we find a SOC constant of λ = 154 ± 2 meV. There is no apparent energy splitting due to trigonal distortion - the A 1 excitation is best fit with a single Lorentzian peak. The low energy region was fit with four peaks as shown in figure 4(a). The elastic line (orange), A 2 (red), and A 3 (green) peaks were fit with Gaussian functions, while the A 1 peak (blue) was fit with a Lorentzian function. The fit also included a background function, shown as a dashed black line, which included a constant offset with a linear slope only in the energy loss region.

We also performed a series of fits of the A 1 feature using two peaks with fixed splitting, shown in figure 4(b-i), to estimate an upper bound on the trigonal splitting | D trig | (figure 1). Each panel

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01190v1-figure-e4f7386a-434f-49de-8c50-a4c3181b7ce2.jpg)

Figure 2. (a) Ru M -edge RIXS spectra on 3 α -RuCl 3 as a function incident energy, E i . The lower (higher) energy features resonate at a lower (higher) energy and correspond to t 2g ( e g ) excitations. The green and magenta bars represent the energy range integrated over to extract the intensity for t 2g and e g respectively. (b) Total electron yield x-ray absorption spectroscopy (TEYXAS) measurement on α -RuCl . The integrated intensities of the 3 t 2g and e g RIXS regions are plotted as a function of E i in green and magenta respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01190v1-figure-e3973c64-e855-44c3-a5b4-d3ab81c9c043.jpg)

data at | D trig | &lt; 40 meV which we define as our upper bound.

Figure 3. Ru M -edge 3 RIXS ( E i = 461 eV) spectrum of α -RuCl 3 single crystals measured at 300 K. A 2 θ = 90 ◦ scattering angle was used and the incident beam was in 20 ◦ grazing incidence with the sample, giving a measurement at q = (0 130) . r.l.u. (black diamond in Brillouin zone inset). The features are labeled using the notation of Sandilands et al [26]. Predicted dd excitation energies from quantum chemistry calculations [27] are shown for the P 3 12 (red) [17] and 1 C /m 2 (blue) [19] structures. The inset shows the Brillouin zone of α -RuCl 3 in hexagonal notation. The dark (light) blue represent areas which can be probed by Ru M 3 -edge RIXS with 2 θ = 90 ◦ (2 θ = 150 ). ◦

corresponds to a fit using a fixed | D trig | ranging from 10 to 80 meV. The relative energy position of the peaks is fixed by this splitting, however they are free to move in energy together during the fit. Both peaks are described by a pseudo-Voigt function with the same intensity, width, and Lorentzian/Gaussian ratio, however these values are fit for each of the individual | D trig | values. We notice the fit diverging from our

The energy difference between the trigonally split j eff = 3 2 levels / is given by | D trig | , however this is not equivalent to the trigonal field energy ∆. Furthermore, the sign of both of these values depends on whether it is trigonal compression (∆ &gt; 0, D trig &gt; 0) or elongation (∆ &lt; 0, D trig &lt; 0). In figure 1 we show the case of trigonal compression where the energy level of Y is above Z , while in the case of trigonal elongation the opposite would be true. Our current measurements cannot differentiate between compression or elongation, therefore we have to estimate ∆ separately for each case. The splitting between the levels is given by D trig /λ = 1 4 [ √ 8 + (1 + δ ) 2 -3 + δ ], where δ = 2∆ /λ [25]. From this we find for compression | ∆ | &lt; 55 meV and for elongation | ∆ | &lt; 65 meV

## 4. Discussion

We have confirmed that the trigonal splitting in α -RuCl 3 is small enough, such that the energy hierarchy needed for a j eff = 1 2 ground / state is satisfied: | ∆ ( | &lt; 0.065 eV) glyph[lessmuch] λ (0.154 eV) glyph[lessmuch] 10 Dq ( ≈ 2 eV). Our results are consistent with the MRCI calculations which predict a splitting of | D trig | = 39 meV for the C /m 2 structure. Ru L-edge XAS reports a vanishingly small linear dichroism effect corresponding to D trig = -12 ± 10 meV [21]. It is interesting to note this report of trigonal elongation is at odds with structural studies [19, 28, 29] which show trigonal compression. Regardless of the sign of D trig , the small magnitude found in this study and by

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01190v1-figure-9e5e62e6-c9f8-48b8-a81d-89b8d1ec7827.jpg)

-

Figure 4. (a) Low-energy region of the RIXS spectrum from figure 3. The elastic (orange), A 2 (red), and A 3 (green) peaks are fit with Gaussian functions while the A 1 peak is fit with a Lorentzian function. A constant background and linear background only in the energy loss region are also included in the fit (dashed black line). The energies extracted from the fit are 233 meV, 529 meV, and 743 meV for A , A , and A . 1 2 3 (b-i) Fit of the A 1 peak using two peaks separated by a fixed ∣ ∣ D trig ∣ ∣ value. The width and intensity of the peaks varies for different ∣ ∣ D trig ∣ ∣ values, but the two peaks are constrained to the same values for each individual fit. The peaks are pseudo-Voigt functions, where both peaks are fixed with the same Lorentzian/Gaussian ratio. The peaks are best fit with 100% Lorentzian for 10-60 meV, 93% for 70 meV, and 78% for 80 meV.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01190v1-figure-221ed526-8416-4cdf-8913-f845a3f69f92.jpg)

Agrestini et al indicates the g factor should be nearly isotropic.

λ = 154 ± 2 meV, agrees almost exactly with the tabulated Ru 3+ free ion value ( λ = 155 meV) [34].

The limit | D trig | &lt; 40 meV is likely an overestimation since our energy resolution was 27 meV and we observed no splitting of the spin-orbit exciton. The use of equal intensity peaks for our fitting in figure 4(b-i) is an approximation, the peaks should actually have different relative intensities depending on the RIXS matrix elements. Therefore, future high-resolution Ru M -edge RIXS studies on 3 α -RuCl 3 can use different polarizations and incident angles to determine more accurately the magnitude and perhaps even the sign of D trig [25, 31].

Optical spectroscopy observes low-energy features at similar energies as our RIXS results: A 1 ≈ 270 meV, A 2 ≈ 530 meV, and A 3 ≈ 740 meV at room temperature [26, 32]. Additionally, an extremely sharp peak was observed with Raman spectroscopy at A 0 ≈ 145 meV ( λ ≈ 96 meV) [26]. The A 0 peak was initially attributed as the spin-orbit exciton, however it is at lower energy than our RIXS results and MRCI calculations [27]. The 2 meV width of the A 0 Raman mode was a surprise since coupling with phonons is expected to broaden the spin-orbit exciton peak [33]. Furthermore, a recent Raman spectroscopy study on α -RuCl 3 finds that the A 0 peak vanishes with increasing temperatures, which is unexpected behavior for a spin-orbit exciton [32]. This study also reports on a higher energy Raman mode at ≈ 235 meV which agrees very well with our RIXS results. Inelastic neutron scattering has observed hints of an excitation around this energy region, albeit with poor statistics [18]. Curiously, the SOC constant we report here,

Previous studies incorrectly assigned the A -A 1 3 peaks as transitions to SOC-split e g states [12, 26, 35]. The resonant behavior we observe, coupled with MRCI calculations [27], demonstrates that their energy is too low to correspond to e g excitations. The A 1 peak observed with optical spectroscopy is the spin-orbit exciton but since it is dipole-forbidden it is a phononassited transition. Indeed, the optical peak is shifted 30-40 meV higher in energy, which is consistent with the energy of optical phonons measured in α -RuCl 3 [15]. The temperature dependence of all three peaks was studied by Sandilands et al [26], integrating the spectral weight from 0.1-0.87 eV, and was consistent with a phonon-assisted mechanism. However, the A 1 peak is dominant in this region and Borgwardt et al [32] found that the A 2 and A 3 peak do not have a strong temperature dependence, i.e. A 2 and A 3 are not phonon-assisted transitions. This is shown clearly by comparing their energies ( ≈ 530 meV and ≈ 740 meV) to our RIXS energies (523 meV and 745 meV), which are nearly identical.

Borgwardt et al [32] interpreted these peaks as multiparticle excitations, i.e. double and triple spinorbit excitons, but this seems unlikely. For example, the energy difference between multiples of a single excitation and a multiple excitation should correspond to a phonon energy, but A 2 -2 · A 1 ≈ 60 meV which is higher than any observed phonon in α -RuCl . 3 As well, we observe the same peaks with RIXS as optical spectroscopy strongly indicating a common origin. Multiparticle excitations are a nonlinear effect and

thus their observation usually requires high-intensity photon beams. This could be possible in RIXS or Raman spectroscopy, however the infrared absorption experiment used a low-intensity lamp [26].

We are not completely certain of the true nature of these two peaks, however we believe that Ru 4 d -Cl 3 p hybridization plays a key role. One picture could be Cl 3 p → Ru 4 d charge-transfer type excitations. For example, recent density functional theory calculations found two sharp peaks in the DOS with significant RuCl hybridization at 590 meV and 730 meV above the Fermi level [21].

To our knowledge, this is the first report of Medge RIXS on a 4 d transition metal system. Actually, there are very few reports of M-edge RIXS at all. The first M-edge RIXS experiment cleverly leveraged the low energy of the Cu M 2 3 , edge to achieve improved resolution on existing instrumentation [36] before the development of next-generation soft x-ray RIXS beamlines [23, 24, 37-40]. Studying 3 d transition metal systems with M-edge RIXS is rare since the lifetime of the shallow core hole has a complex dependence on incident photon energy, the elastic line due to offspecular reflectivity is extremely strong and obscures low-energy features, and the inelastic cross-section is lower due to increased Auger emission [41].

These first two disadvantages are related to the &lt; 120 eV incident energy used for M-edge RIXS of 3 d systems and not so important for 4 d systems which have higher energy M edges. In our experiment, the elastic line intensity was three times the A 1 intensity, but we measured at 300 K and we found considerable quasi-elastic weight (50 meV FWHM vs. 27 meV FWHM resolution). These results are promising but further studies at higher resolutions and lower temperatures, as well as varying the scattering angle away from 2 θ = 90 , ◦ are important to determine if elastic intensity will be a limiting factor of M-edge RIXS in 4 d systems. The inelastic cross-section is generally lower for M vs. L edges due to the decreased fluorescence yield. For example, in Ru the fluorescence yield is 2 3 . · 10 -4 and 4 5 . · 10 -2 for the M and L 3 3 edge respectively [42], however our results show that the resonant enhancement is still sufficient to perform experiments.

One disadvantage is that the lower energy of the M edge limits the area of the Brillouin zone which can be probed. The inset of figure 3 shows the area available with Ru M -edge RIXS in 3 α -RuCl 3 as dark (light) blue for 2 θ = 90 ◦ (2 θ = 150 ), while Ru L -◦ 3 edge RIXS can probe a few Brillouin zones. We note that the accessible Brillouin zone in α -RuCl 3 at the M edge is still enough to search for the gapless Majorana fermions predicted by theory [43]. As well, in general the larger supercells of magnetically ordered materials will have a correspondingly smaller magnetic Brillouin zone which could possibly be probed completely at the M edge. However, the decreased SOC and lifetime of the 3 p core hole in M-edge RIXS does make it less effective at measuring magnons [44].

Nonetheless, M-edge RIXS for 4 d systems has the enormous advantage that many soft x-ray RIXS beamlines already exist or are under development which can provide high flux and sub-30 meV resolution [23, 24, 37-40]. For Ru L -edge RIXS (2840 eV) there 3 is currently only one instrument in the world, IRIXS at P1/DESY, which is currently unavailable to general users. The latest results from IRIXS [45, 46] have a resolution approximately a factor of 5 worse than the resolution routinely available for Ru M 3 -edge RIXS.

## 5. Conclusion

We have performed Ru M -edge 3 resonant inelastic x-ray scattering (RIXS) on α -RuCl . 3 We observe dd excitations in agreement with optical spectroscopy and quantum chemistry calculations. Our observation of a spin-orbit exciton allows us to extract a very accurate value for the spin-orbit coupling constant λ = 154 ± 2 meV. The spin-orbit exciton shows no splitting due trigonal distortions and overall we find the energy hierarchy necessary for j eff = 1 2 physics / is satisfied in α -RuCl . 3 Our results resolve some previous misconceptions about the electronic structure of α -RuCl 3 and provide a springboard for future calculations and experiments to further elucidate its nature. Measurement of 4 d systems with M-edge RIXS is a novel technique which we believe will be an important part of the x-ray spectroscopist's toolbox since it allows ultra high-resolution RIXS measurements here and now.

## Acknowledgments

Work at the University of Toronto was supported by the Natural Science and Engineering Research Council (NSERC) of Canada, Canadian Foundation for Innovation, and Ontario Innovation Trust. This research used the SIX beamline of the National Synchrotron Light Source II, a U.S. Department of Energy (DOE) Office of Science User Facility operated for the DOE Office of Science by Brookhaven National Laboratory under Contract No. DE-SC0012704. This work was performed in part at Aspen Center for Physics, which is supported by National Science Foundation grant PHY-1607611. BWL acknowledges the support from the University of Toronto Faculty of Arts and Sciences Postdoctoral Fellowship.

## References

- [1] Kitaev A 2006 Annals of Physics 321 2 - 111 ISSN 0003-4916 january Special Issue
- [2] Trebst S 2017 arXiv preprint arXiv:1701.07056
- [3] Winter S M, Tsirlin A A, Daghofer M, van den Brink J, Singh Y, Gegenwart P and Valenti R 2017 Journal of Physics - Condensed Matter 29 ISSN 0953-8984
- [4] Hermanns M, Kimchi I and Knolle J 2018 Annual Review of Condensed Matter Physics 9 17-33
- [5] Jackeli G and Khaliullin G 2009 Phys. Rev. Lett. 102 (1) 017205
- [6] Chaloupka J, Jackeli G and Khaliullin G 2010 Phys. Rev. Lett. 105 (2) 027204
- [7] Singh Y and Gegenwart P 2010 Phys. Rev. B 82 064412
- [8] Singh Y, Manni S, Reuther J, Berlijn T, Thomale R, Ku W, Trebst S and Gegenwart P 2012 Phys. Rev. Lett. 108 (12) 127203
- [9] Modic K A, Smidt T E, Kimchi I, Breznay N P, Biffin A, Choi S, Johnson R D, Coldea R, Watkins-Curry P, McCandless G T, Chan J Y, Gandara F, Islam Z, Vishwanath A, Shekhter A, McDonald R D and Analytis J G 2014 Nature Communications 5 4203
- [10] Takayama T, Kato A, Dinnebier R, Nuss J, Kono H, Veiga L S I, Fabbris G, Haskel D and Takagi H 2015 Phys. Rev. Lett. 114 (7) 077202
- [11] Takagi H, Takayama T, Jackeli G, Khaliullin G and Nagler S E 2019 Nature Reviews Physics 1 264-280
- [12] Plumb K W, Clancy J P, Sandilands L J, Shankar V V, Hu Y F, Burch K S, Kee H Y and Kim Y J 2014 Phys. Rev. B 90 (4) 041112
- [13] Banerjee A, Yan J, Knolle J, Bridges C A, Stone MB, Lumsden M D, Mandrus D G, Tennant D A, Moessner R and Nagler S E 2017 Science 356 1055-1059 ISSN 0036-8075
- [14] Do S H, Park S Y, Yoshitake J, Nasu J, Motome Y, Kwon Y S, Adroja D T, Voneshen D J, Kim K, Jang T H, Park J H, Choi K Y and Ji S 2017 Nat. Phys. 13 1079-1084
- [15] Sandilands L J, Tian Y, Plumb K W, Kim Y J and Burch K S 2015 Phys. Rev. Lett. 114 (14) 147201
- [16] Kasahara Y, Ohnishi T, Mizukami Y, Tanaka O, Ma S, Sugii K, Kurita N, Tanaka H, Nasu J, Motome Y, Shibauchi T and Matsuda Y 2018 Nature 559 227-231 ISSN 1476-4687
- [17] Stroganov E V and Ovchinnikov K V 1957 Ser. Fiz. i Khim. 12 152
- [18] Banerjee A, Bridges C A, Yan J Q, Aczel A A, Li L, Stone M B, Granroth G E, Lumsden M D, Yiu Y, Knolle J, Bhattacharjee S, Kovrizhin D L, Moessner R, Tennant D A, Mandrus D G and Nagler S E 2016 Nature Materials 15 733-740
- [19] Cao H B, Banerjee A, Yan J Q, Bridges C A, Lumsden M D, Mandrus D G, Tennant D A, Chakoumakos B C and Nagler S E 2016 Phys. Rev. B 93 (13) 134423
- [20] Gretarsson H, Clancy J P, Liu X, Hill J P, Bozin E, Singh Y, Manni S, Gegenwart P, Kim J, Said A H, Casa D, Gog T, Upton M H, Kim H S, Yu J, Katukuri V M, Hozoi L, van den Brink J and Kim Y J 2013 Phys. Rev. Lett. 110 (7) 076402
- [21] Agrestini S, Kuo C Y, Ko K T, Hu Z, Kasinathan D, Vasili H B, Herrero-Martin J, Valvidares S M, Pellegrin E, Jang L Y, Henschel A, Schmidt M, Tanaka A and Tjeng L H 2017 Phys. Rev. B 96 (16) 161107
- [22] Sears J A, Songvilay M, Plumb K W, Clancy J P, Qiu Y, Zhao Y, Parshall D and Kim Y J 2015 Phys. Rev. B 91 (14) 144420
- [23] Dvorak J, Jarrige I, Bisogni V, Coburn S and Leonhardt W 2016 Review of Scientific Instruments 87 115109
- [24] Jarrige I, Bisogni V, Zhu Y, Leonhardt W and Dvorak J 2018 Synchrotron Radiation News 31 7-13
- [25] Chaloupka J and Khaliullin G 2016 Phys. Rev. B 94 (6) 064435
- [26] Sandilands L J, Tian Y, Reijnders A A, Kim H S, Plumb K W, Kim Y J, Kee H Y and Burch K S 2016 Phys. Rev. B 93 (7) 075144
- [27] Yadav R, Bogdanov N A, Katukuri V M, Nishimoto S, Van Den Brink J and Hozoi L 2016 Scientific reports 6 37925
- [28] Johnson R D, Williams S C, Haghighirad A A, Singleton J, Zapf V, Manuel P, Mazin I I, Li Y, Jeschke H O, Valent´ ı R and Coldea R 2015 Phys. Rev. B 92 (23) 235119
- [29] Sears J A 2017 Neutron and X-ray Diffraction Studies of Structure and Magnetism in α -RuCl 3 Ph.D. thesis University of Toronto
- [30] Kim J, Casa D, Upton M H, Gog T, Kim Y J, Mitchell J F, van Veenendaal M, Daghofer M, van den Brink J, Khaliullin G and Kim B J 2012 Phys. Rev. Lett. 108 (17) 177003
- [31] Kim J, Daghofer M, Said A H, Gog T, van den Brink J, Khaliullin G and Kim B J 2014 Nature Communications 5 4453
- [32] Borgwardt N 2019 Optics on materials with strong spin-orbit coupling: topological insulators Bi2xSbxTe3-ySey and the j=1/2 compounds Na2IrO3

and alpha-RuCl3 Ph.D. thesis Universit¨t zu K¨ln a o URL https://kups.ub.uni-koeln.de/9187/

- [33] Plotnikova E M, Daghofer M, van den Brink J and Wohlfeld K 2016 Phys. Rev. Lett. 116 (10) 106401
- [34] Porterfield W 2013 Inorganic Chemistry (Elsevier Science) ISBN 9780323138949
- [35] Reschke S, Mayr F, Wang Z, Do S H, Choi K Y and Loidl A 2017 Phys. Rev. B 96 (16) 165120
- [36] Kuiper P, Guo J H, S˚the C, Duda L C, Nordgren a J, Pothuizen J J M, de Groot F M F and Sawatzky G A 1998 Phys. Rev. Lett. 80 (23) 5204-5207
- [37] Lai C H, Fung H S, Wu W B, Huang H Y, Fu H W, Lin S W, Huang S W, Chiu C C, Wang D J, Huang L J, Tseng T C, Chung S C, Chen C T and Huang D J 2014 Journal of Synchrotron Radiation 21 325-332
- [38] Brookes N, Yakhou-Harris F, Kummer K, Fondacaro A, Cezar J, Betto D, Velez-Fort E, Amorese A, Ghiringhelli G, Braicovich L, Barrett R, Berruyer G, Cianciosi F, Eybert L, Marion P, van der Linden P and Zhang L 2018 Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 903 175 - 192 ISSN 0168-9002
- [39] Diamond light source: I21 beamline homepage https://www.diamond.ac.uk/Instruments/ Magnetic-Materials/I21.html
- [40] Max iv: Veritas beamline homepage https: //www.maxiv.lu.se/accelerators-beamlines/ beamlines/veritas/
- [41] Wray L A, Huang S W, Jarrige I, Ikeuchi K, Ishii K, Li J, Qiu Z Q, Hussain Z and Chuang Y D 2015 Frontiers in Physics 3 32 ISSN 2296-424X
- [42] Bambynek W, Crasemann B, Fink R W, Freund H U, Mark H, Swift C D, Price R E and Rao P V 1972 Rev. Mod. Phys. 44 (4) 716-813
- [43] Hal´ asz G B, Perkins N B and van den Brink J 2016 Phys. Rev. Lett. 117 (12) 127203
- [44] Ament L J P, van Veenendaal M, Devereaux T P, Hill J P and van den Brink J 2011 Rev. Mod. Phys. 83 (2) 705-767
- [45] Suzuki H, Gretarsson H, Ishikawa H, Ueda K, Yang Z, Liu H, Kim H, Kukusta D, Yaresko A, Minola M, Sears J A, Francoual S, Wille H C, Nuss J, Takagi H, Kim B J, Khaliullin G, Yavas H and Keimer B 2019 Nature Materials 18 563567 ISSN 1476-4660
- [46] Gretarsson H, Suzuki H, Kim H, Ueda K, Krautloher M, Kim B J, Yava¸ s H, Khaliullin G and Keimer B 2019 Phys. Rev. B 100 (4) 045123
Blair W Lebert Subin Kim Valentina Bisogni Ignace Jarrige Andi M Barbour Young-June Kim
Resonant inelastic x-ray scattering study of $α$-RuCl$_3$: a progress report
Ru M$_3$-edge resonant inelastic x-ray scattering (RIXS) measurements of
RuCl$_3$ with 27 meV resolution reveals a spin-orbit exciton without noticeable
splitting. We extract values for the spin-orbit coupling constant
($\lambda=154\pm2$ meV) and trigonal distortion field energy
($\left|\Delta\right|<65$ meV) which support the $j_{\rm eff}=1/2$ nature of
RuCl$_3$. We demonstrate the feasibility of M-edge RIXS for $4d$ systems, which
allows ultra high-resolution RIXS of $4d$ systems until instrumentation for
L-edge RIXS improves.