## Polarization-resolved Raman spectroscopy of α -RuCl 3 and evidence of room temperature two-dimensional magnetic scattering

Thuc T. Mai, 1, 2, ∗ A. McCreary, 2 P. Lampen-Kelley, 3, 4 N. Butch, 5, 6 J. R. Simpson, 7, 2

J.-Q. Yan, 4, 3 S. E. Nagler, 8 D. Mandrus, 3, 4 A. R. Hight Walker, 2 and R. Vald´ es Aguilar 1, †

1 Department of Physics, The Ohio State University. Columbus, OH 43210, USA

2 Nanoscale Device Characterization Division, Physical Measurement Laboratory, NIST, Gaithersburg, MD

3 Materials Science and Engineering, University of Tennessee, Knoxville, TN

4 Materials Science and Technology Division Oak Ridge National Laboratory, Oak Ridge, TN

5 NIST Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA

6 Center for Nanophysics and Advanced Materials, Department of Physics,

University of Maryland, College Park, MD 20742 USA

7 Physics, Astronomy, and Geosciences, Towson University, Towson, MD

8 Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN

(Dated: August 6, 2019)

Polarization-resolved Raman spectroscopy was performed and analyzed from large, high quality, mono-domain single crystal of α -RuCl , a proximate Kitaev quantum spin liquid. 3 Spectra were collected with laser polarizations parallel and perpendicular to the honeycomb plane. Pairs of nearly degenerate phonons were discovered and show either a 4-fold or 2-fold polarization angle dependence in their Raman intensity, thereby providing evidence to definitively assign the bulk crystal point group as C 2 h . The low frequency continuum that is often attributed to scattering from pairs of Majorana fermions was also examined and found to disappear when the laser excitation and scattered photon polarizations were perpendicular to the honeycomb plane. This disappearance, along with the behavior of the phonon spectrum in the same polarization configuration, strongly suggests that the scattering continuum is 2-dimensional. We argue that this scattering continuum originates from the Kitaev magnetic interactions that survives up to room temperature, a scale larger than the bare Kitaev exchange energy of approximately 50 K.

## I. INTRODUCTION

in both α -RuCl 3 and β and γ -Li IrO 2 3 10-12 .

In the last decade, quantum spin liquids (QSL) research has developed into a burgeoning field in condensed matter physics 1-6 . Due to the lack of spin ordering and the massively entangled ground state, a QSL is predicted to host a variety of exotic phenomena . 3 Because of the exact analytical solution of its ground state and excitations, the Kitaev model 2 is one of the most studied theoretical models of a QSL. Recent experimental and theoretical studies have converged onto a few real systems that could host this model 7-17 , including α -RuCl 3 and the iridate family of A IrO . In these systems, strong 2 3 spin-orbit coupling, which leads to J eff = 1 2, and the / honeycomb lattice of the edge-sharing IrO 6 and RuCl 6 octahedra give rise to effective bond-dependent exchange interaction of the Kitaev model 18,19 .

Although long range magnetic order does occur in α -RuCl 3 below 7 K 20 , inelastic neutron scattering (INS) revealed contributions from both Heisenberg and Kitaev exchange . 7 A magnetic scattering continuum measured by various INS studies 8,9,21 was suggested to arise from excitations of Majorana fermions, a signature of the Kitaev QSL phase 2,4,5 . Others have argued that this continuum comes from magnons that decay rapidly as a consequence of frustrated interactions 22 . Interestingly, theoretical studies have shown that inelastic scattering with light (such as Raman scattering) can also couple to the fractionalized spin excitations 13,14,23 . The predicted broad Raman scattering signal has been measured

The space group symmetry of bulk crystal α -RuCl 3 is C /m 2 (point group C 2 h ) at room temperature 24,25 . The a-b plane consists of Ru 3+ ions that form a 2-dimensional (2D) honeycomb lattice. The Ru ions reside inside Cl -octahedra, where due to crystal field splitting and spin orbit coupling, have pseudo-spin values of 1/2 18 . These honeycomb planes are monoclinically stacked along the c axis (Fig. 1(a) and (b)). The c axis makes an approximately 108 ° angle with the a axis while being perpendicular to the b axis, and the van der Waals coupling weakly holds the 2D layers together. The monoclinic stacking breaks the 3-fold rotational symmetry of the honeycomb layers, leaving the 2-fold rotational symmetry around the b -axis as the highest symmetry of the bulk crystal. Refinement of diffraction data from previous studies revealed a small distortion of the honeycomb plane, in the order of ∼ 0.01 ˚ , effectively also breaking A the individual layer's 3-fold symmetry 24,25 .

Temperature-dependent Raman scattering experiments and theory provided additional evidence that this scattering continuum follows a fermionic behavior at low temperature, as opposed to bosonic 10,12,26,27 . Scattering of a lattice vibrational mode in the proximity of this continuum exhibits a highly asymmetric lineshape 10,11 . However, a determination of the symmetry and frequency of all of the bulk crystal's Raman phonons has not been uniquely made. Armed with a high quality single crystal sample of α -RuCl , our Raman scattering experiment 3 is able to definitively determine the symmetry of the

TABLE I. Amplitude of the Raman response of the a-b and a-c plane of α -RuCl , derived as 3 e i · R S · e s , where e i and e s are the incident and scattered polarization directions, respectively, and the raman tensor R S are given in Eq. 1. Here the crystalline b axis is the tensor y axis, and a is the x axis, while c is not the z axis. N.A. means that no Raman response should occur for that configuration. θ and φ are defined in the text.

|     |     | Parallel                                       | Crossed                             |
|-----|-----|------------------------------------------------|-------------------------------------|
| a-b | A g | ( a + b ) +( a - b ) cos(2 θ )                 | ( b - a ) sin(2 θ )                 |
|     | B g | e sin(2 θ )                                    | e cos(2 θ )                         |
| a-c | A g | ( a + c ) +( a - c ) cos(2 φ ) - 2 d sin(2 φ ) | ( a - c ) sin(2 φ ) - 2 d cos(2 φ ) |
| a-c | B g | N.A.                                           | N.A.                                |

Raman-active phonons. In addition, multiple phonon peaks, not observed in previous studies, with ≈ 2 cm -1 separation in frequency were resolved, including the one affected by the continuum. Then, we systematically study the change in the asymmetry of this phonon, with polarization-resolved Raman spectroscopy. We find that the continuum and the phonon asymmetry are not observed when the photon polarizations are perpendicular to the honeycomb plane. This strongly suggests that these effects are of 2D magnetic origin, even at room temperature, a scale larger than the bare Kitaev exchange of ≈ 5-8 meV ∼ 58-90 K 7,19,26-28 .

## II. EXPERIMENT

The single crystal of α -RuCl 3 used here was grown as described elsewhere 24 . Similar crystals have been studied by neutron scattering 7-9 , terahertz 29 and Raman 16,27 spectroscopies. Laue diffraction patterns were taken on various spots across the surface of the sample ( a-b plane) at room temperature to confirm the large domain, single crystal nature of our sample (Fig. 1(c), (d), (e)). Raman spectra with 514.5 nm laser excitation were measured in ambient conditions using a single-grating spectrometer (grating with lines density of 1800 mm -1 , and 1.27 cm CCD detector) in the 180 ° backscattering configuration. A combination of half wave plates and linear polarizers was used to measure spectra in either parallel or crossed polarization configurations, where the polarization orientations were confirmed using MoS 2 as a reference 30 . The α -RuCl 3 sample, where the a and b axes were known from Laue diffraction, was placed on a rotating stage to access different crystal orientations. Enabled by the large thickness of our single crystal, we also measured Raman spectra from the a-c plane by mounting the crystal on its side. Integration times ranged between 5 minutes ( a-b plane) and 12 minutes ( a-c plane), and the laser power was kept below 200 µ W through the objective (50x, numerical aperture 0.75 to probe a-b plane, 100x long working distance objective, numerical aperture 0.75 to probe a-c plane) to avoid any local heating of the sample.

Given the point group symmetry of the bulk crystal is C 2 h , we expect 12 Raman active phonons: 6 A g + 6 B g .

Group theory predicts the Raman response tensor of A g and B g to be the following:

$$\mathbf R _ { A _ { g } } = \begin{bmatrix} a & 0 & d \\ 0 & b & 0 \\ d & 0 & c \end{bmatrix} ; \mathbf R _ { B _ { g } } = \begin{bmatrix} 0 & e & 0 \\ e & 0 & f \\ 0 & f & 0 \end{bmatrix} \quad \quad ( 1 )$$

where a, b, c, d, e, f are independent coefficients. The 3 × 3 matrices in Eq. 1 are written in a coordinate system

FIG. 1. (color online) Crystal structure and Laue X-ray diffraction pattern . Crystal structure of α -RuCl 3 projected onto the a-c plane (a) and a-b plane (b). The solid black lines represent the unit cell. The orange balls represent Cl atoms, while the larger blue balls represent Ru atoms. (c) Picture of studied sample; the red circles are the approximate location of where the Laue patterns in (d) and (e) were taken. The well defined circular dots indicate the single-crystallinity of our sample. This pattern remained unchanged over several millimeters on the a-b plane surface. The relationship between the crystallographic axes (a,b,c) and the laboratory (Raman tensor) coordinate (x,y,z) is shown in (f). Notice that c and z are not parallel.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01021v1-figure-d550885c-8b00-4dfb-8645-9c69b40626f9.jpg)

where the tensor's x and y axes are parallel to the crystallographic a and b axes respectively. The selection rules allow for both A g and B g phonons to be measured from the a-b plane. Furthermore, A g and B g phonons can be measured simultaneously in both crossed and parallel polarization configurations when the polarizations are not coincident with either the a or b crystallographic axis.

The predicted angular dependence of the Raman response in the backscattering geometry is shown in Table I, where θ is the angle between the incoming polarization and the crystal's b axis for the measurement in the a-b plane (inset in Fig. 2), and φ is the angle away from the a axis for the a-c plane measurement (inset in Fig. 5). The Raman intensity is given by the square of the expressions given in Table I: I Raman = | e i · R e · s | 2 , where e i and e s are the incident and scattered polarization directions, respectively. We expect the intensity of the A g and B g phonons to be 4-fold modulated, and out of phase with each other as θ changes within the a-b plane. In contrast in the a-c plane, B g phonons are not allowed and the A g phonons should only show 2-fold modulated intensities.

## III. DATA

We obtained multiple Laue X-ray diffraction patterns across our mm sized sample, where two examples are shown in Fig. 1 (d) and (e). The consistency of the diffraction patterns at various locations confirmed that our sample is not only of single crystal nature, but also that the domain size is in the order of a few mm. The Laue diffractometer determined the highest symmetry axis of our crystal, the b axis with 2-fold rotational symmetry, in Fig. 1(d), (e). With the a-b plane defined as the basal plane, we can uniquely determine the other two crystal axes, a and c . The unique b axis places constraints on any rank 2 property tensor (i.e. Raman response tensor) by Neumann's principle 31 . From this constraint, the tensor y-coordinate (in Eq. 1) is identically the crystal's b axis. We then define the x axis to be the a axis. This leaves the laboratory frame's z axis to be different from the crystal's c axis (Fig. 1(f)).

Using this knowledge of the crystallographic axes, our initial measurement of the a-b plane collected with incident polarization parallel to b is presented in Fig. 2. The labels ¯( z bb z ) and ¯( z ab z ) represent the experimental configurations of the direction of propagation of the incoming beam, the polarization of incident and scattered photons, and the propagation direction of the scattered beam. At first glance, we observed similar spectra to previous Raman studies on α -RuCl 3 10,11 . When we carefully compare the two traces in Fig. 2, between the parallel and crossed polarization configuration around 270 cm -1 and 295 cm -1 , there are clear differences in the phonon peak frequencies of approximately 2 cm -1 . Similarly, two distinct and weaker phonons can be seen around 119 cm -1 . Our measurements show what appears to be two asymmetric phonon peaks around 164 cm -1 (Fig. 2) and a broad background. An asymmetric Fano lineshape typically arises when a discrete resonance, a phonon in this case, is coupled to a broad continuum 32 . This scattering continuum in α -RuCl 3 has been subjected to multiple low temperature studies 10,11 , and its origin at room temperature will be discussed in the next section.

To determine the phonon frequencies, the peaks in the Raman spectra were fit with Lorentzian lineshapes, with the exception of the asymmetric peaks near 164 cm -1 , where the Fano lineshape:

$$I ( \omega ) = I _ { 0 } \frac { 1 } { ( 1 + q ^ { 2 } ) } \frac { ( q + ( \frac { \omega - \omega _ { 0 } } { \Gamma } ) ) ^ { 2 } } { 1 + ( \frac { \omega - \omega _ { 0 } } { \Gamma } ) ^ { 2 } } \quad \quad ( 2 )$$

was used instead. Here, q represents the asymmetry parameter. A Lorentzian lineshape is recovered as q goes to infinity (or 1 / q | | = 0). I 0 is an overall intensity scale, Γ is the intrinsic width (inverse lifetime) of the phonon, and ω 0 is the natural frequency of the mode. We label the phonons with the irreducible representations A g and B g of the point group C 2 h in Tab. II (and Fig. 2) in contrast to previous reports. The Fano parameters (1 / q | | ) found for parallel ( A 2 g ) and crossed ( B 2 g ) configuration are 0 064 . and 0 091, respectively. .

With fixed polarizers and rotating the sample, we collected Raman spectra as a function of the angle between the laser polarizations and the crystalline axes. In the a-b plane, oscillations in the scattered intensity in both

FIG. 2. Raman intensity of the a-b plane . The parallel configuration data are offset for clarity. Here z is the lab frame direction in which the incident and scattered beams propagate (¯ = z -z ) and a and b are the crystallographic axes. A clear difference is observed between some of the phonon frequencies in the crossed and parallel configuration. We labeled these phonons with A g and B g , of the bulk point group symmetry C 2 h . The * mode seen in ¯( z ab z ) is the A 6 g phonon signal leakage due to imperfect polarizers. The inset shows a schematic of the a-b plane measurement, where θ is defined as the angle between the incident polarization and the b axis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01021v1-figure-5a09756a-42dc-454e-89e7-29242eea8e1e.jpg)

TABLE II. List of A g and B g phonon frequencies in units of cm -1 . The symmetries of the phonons were determined by their intensities when e i or e s were along a crystalline axis. In such polarization configuration, A g modes appear in parallel and B g in crossed. The † denotes pairs of nearly degenerate phonon modes, not observed in previous publications. B 6 g was deliberately labeled out of order to make the nearly degenerate phonon pairs clearer. N.O. means not observed.

| Symmetry Phonon   |   1 † |   2 † | 3     |   4 † |   5 † |     6 |
|-------------------|-------|-------|-------|-------|-------|-------|
| A g               | 119.4 | 164.2 | 222.0 | 271.9 | 295.8 | 313.1 |
| B g               | 115.8 | 165.8 | N.O.  | 269.3 | 297.6 | 291.5 |

crossed and parallel configuration (Fig. 3), around 164 cm -1 , 270 cm -1 , and 295 cm -1 were observed. In contrast, the A 6 g phonon remained constant in frequency. Each peak frequency oscillation can be explained by a pair of A g and B g phonons. The phonons in each pair change in intensity out of phase with each other as the sample is rotated which results in the apparent oscillation of a phonon peak frequency. We employed a fitting model that contains Lorentzians with frequencies fixed by the A g and B g phonon frequencies. Three example spectra (from crossed polarizations configuration) and their fitted curves are plotted in Figure 4. The spectral weight of the A 4 g , B 4 g , A 5 g , B 5 g phonons are also plotted as a function of the angle (Fig. 4). The solid lines are group theory predicted responses, Table I, which are in excellent agreement with our data. In order to measure only the A g or B g phonons in the a-b plane, the incident and scattered polarizations have to be parallel to either the a or b axis of the crystal.

with ≈ 2 cm -1 difference in their frequencies between crossed and parallel configuration. In these studies, the authors assumed the point group D 3 d of a single perfect honeycomb layer to label the α -RuCl 3 Raman phonon spectra with 4 E g +2 A 1 g . We postulate that the samples used in those studies either contain multiple domains in the volume probed by their laser, or the laser polarizations are not parallel to a crystalline axis. In either case, broader, averaged phonon peaks would appear in place of the separate A g and B g modes as seen in Fig. 4 (b). It is interesting to see how the conclusions from these

Our angular dependent measurements of the a-c plane confirmed the group theoretical prediction, showing no oscillating phonon peaks. In the backscattering geometry, only A g phonons are active in this plane (see eq. 1). Surprisingly, the A 2 g phonon showed a dramatic shift in the frequency of its maximum intensity as a function of sample orientation φ (Fig. 5). Additionally, its asymmetry changed dramatically, becoming nearly symmetric at certain angles, e.g. 1 / q | | → 0 when φ = 90 ° (see inset of Fig. 5). This angle corresponds to the incident and scattered laser polarizations being perpendicular to the a-b plane. Also at this angle, the background/scattering continuum below the phonon frequency vanishes.

## IV. DISCUSSION

## A. The effects of symmetry reduction on Raman phonons

The observations above (the frequency separation of the nearly degenerate phonons, and the polarizationdependent details that allow us to determine the point group of the crystal), are made possible only by the large, high quality, single crystal sample with crystallographic domain size larger than a few mm. Previous studies 10,11,16,27,33,34 did not report multiple phonons

FIG. 3. Raman spectra in the a-b plane as a function of rotation angle θ . False color maps of the Raman intensity as a function of frequency and angle θ , shown for the crossed (a), and parallel (b), configurations. θ = 0 corresponds to e i , e s ‖ b in the parallel configuration, and e i ‖ a and e s ‖ b in the crossed configuration.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01021v1-figure-662a0e19-1bee-4b7c-9e13-de6e4e2b7708.jpg)

FIG. 4. Angular dependence of A 4 g , B 4 g , A 5 g , B 5 g phonons . Panels (a), (b), (c) show the raw data (dots) and fitted Lorentzian (lines) at certain angles in the crossed polarization configuration. The spectral weight outputs from the data fit are shown in (d) and (e) ((f) and (g)) for A 4 g and B 4 g ( A 5 g and B 5 g ) for crossed ((d) and (f)) and parallel ((e) and (g)) polarization configurations. The solid lines in ((d)-(g)) are obtained by fitting the spectral weights with the modulus squared of the corresponding expressions found in Table I. All polar plots are on the same radial scale. The error bars reprensent one standard deviation, obtained by the fitting procedure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01021v1-figure-f8de34de-a289-4117-ae20-461f07239546.jpg)

previous experiments would change with the full details of the Raman phonon spectrum.

## B. Scattering continuum and the Fano lineshape

By comparing our result to other Raman studies 10,16,33,34 , we noticed that the pairs of phonons A 1 g and B 1 g , A 2 g and B 2 g , A 4 g and B 4 g , A 5 g and B 5 g were labeled as E g modes. Factor group analysis reveals that the doubly degenerate E g phonons of D 3 d , the often assumed symmetry of α -RuCl 3 crystals, would split into the A g and B g modes of C 2 h , the point group symmetry of bulk α -RuCl 3 24,25 . Our observation of the nearly degenerate A g and B g phonons confirms this picture. It also suggests that the energy scale of the symmetry breaking effect that splits these modes is small glyph[lessorsimilar] 1 meV.

Similar splitting of the E g modes in the phonon spectra were calculated by density functional theory 35 and confirmed experimentally 36 in CrI , a material that shares 3 the same point group symmetry with α -RuCl 3 at room temperature. Larson and Kaxiras 35 conclude that the energy differences between the A g and B g pairs of phonons come from the weak interlayer coupling. This is consistnet with the small frequency differences that we have measured. We also expect the distortion of each honeycomb layer, where two of the bonds of the hexagon differ from the other four, to play a similar role in the phonon splitting. Although the magnitude of this lattice distortion is minute ( ≈ 0 2%) . 24,25 , it is not inconsistent with our results. Further studies are necessary to clarify the contribution from each mechanism responsible for the observed splitting.

Previous reports have characterized the Raman intensity around 164 cm -1 with a Fano peak and a broad scattering continuum. Our measurements of the a-b plane show that there are two Raman phonons near 164 cm -1 , with two different asymmetry parameters 1 / q | | as shown in Figure 3. It was argued that at low temperature, such scattering continuum comes from the magnetism in α -RuCl 3 10,11 . In this picture, it is expected that the Raman scattering process measures the excitation of the itinerant Majorana fermions in the Kitaev QSL, and can appear as a broad feature in the spectrum 4,13,14,17,23,26 . Our observation in the a-c plane (Fig. 5) where only the A 2 g phonon is allowed, showed three dramatic changes in the phonon Raman response as the sample was rotated around the b axis: the shift of the phonon peak frequency, the vanishing of the continuum and of 1 / q | | when the laser polarization is aligned perpendicular to the honeycomb plane. The Fano line shape is expected to show a large shift in the measured peak frequency when the energy of the discrete state is comparable to the energy of the background 32 . Indeed, we observed a difference of about 2 5 cm . -1 in the bare phonon frequency of A 2 g from fitting with eq. 2 at φ = 0 ° and 90 . ° Therefore, the Fano line shape and the continuum share the same origin, as they disappear together when φ = 90 . ° These phenomena are consistent with the idea that the continuum cannot be excited with e i,s ⊥ a-b plane, and thus suggests that the scattering continuum seen in our data is a 2D effect.

FIG. 5. Changes in the asymmetry of the A 2 g phonon . Here, φ is defined as the angle the incoming polarization makes with the a axis, as shown in the left inset. When the laser polarizations are perpendicular to the a-b plane ( φ =90 ), the ° A 2 g phonon appears to be fully symmetric. The dashed arrow highlights both the change in asymmetric lineshape and the disappearance of the scattering continuum. The phonon frequencies obtained by fitting φ = 0 ° and 90 ° differs by about 2 5 cm . -1 . The inset shows the Fano parameter 1 /q . We expect 1 /q to vanish when the phonon becomes fully symmetric.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01021v1-figure-2756e6fe-4208-4512-98dc-949d47fd50cc.jpg)

While the continuum disappears, the A g phonons are still present. Combined with the presence of the continuum when B g phonons are not active ( ¯ aa b b ( ) ), we conclude that the continuum cannot come from any phononic contribution of the A g or B g modes. We can also rule out the single-particle excitation of free electrons contributing to the continuum. We would expect the scattering of electrons between the weakly coupled van der Waals layers to be weak, leading to a 2D-like quasi-elastic response. However, this material is a Mott insulator, with an optical gap of about 1 eV 37-40 . Temperature dependent electrical resistivity measurement further confirm the insulating behavior of α -RuCl 3 33 . The lack of free carriers is also evident in various far-infrared and terahertz transmission spectroscopy experiments, especially through ∼ mm thick single crystals of α -RuCl 3 29,40-42 . In addition, another candidate to host the Kitaev QSL phase, β -Li IrO , showed similar broad scattering back2 3 ground at low frequency along with multiple asymmetric phonons, even at room temperature 12 .

There have been several studies of Raman scattering in the Kitaev model. They predict a scattering continuum coming from the excitation of the Majorana fermions, in both the gapped, and gapless Kitaev QSL ground states 4,13,14,17,23,26 . In particular, Perreault et al. 23 showed that Raman scattering from Kitaev exchange is only absent for incoming and outgoing polarizations perpendicular to the honeycomb plane. At the same time, the response in the a -c plane is not zero if both the 3-fold symmetry of the perfect honeycomb lattice is bro- ken, and both direct and mediated superexchange are considered in the Kitaev interaction. Our observations support this picture, as the phonon asymmetry and the continuum are present in all polarization combinations, except when they both point perpendicular to the a -b plane, and by the fact that our A g -B g phonon splitting means that the 3-fold symmetry is broken. Furthermore, it is expected that the Kitaev scattering continuum persist up to temperatures of at least 7 × J K , where J K is the isotropic Kitaev exchange constant 26,27 . The value of J K has been estimated to be 5 meV to 8 meV (58 K to 93 K) 7,10,19,28 , which places the continuum persistence temperature well above 300 K. Combined with the polarization-resolved results described above, it strongly suggests that the scattering continuum measured in bulk α -RuCl 3 comes from the Kitaev QSL physics at room temperature.

## V. SUMMARY

Using polarization-sensitive Raman scattering of the phonon spectrum in α -RuCl , we showed that the point 3 group symmetry of the crystal is C 2 h , and not the D 3 d point group of its perfect honeycomb layer. Because of the large, high-quality single crystal used in this experiment, we correctly identified and labeled the full Raman active phonon spectrum. We observed a broad scattering continuum and two Fano line shape resonances near 164 cm -1 . The scattering continuum and the Fano asymmetry were found to vanish with laser polarizations perpendicular to the honeycomb plane. We argued that the Kitaev exchange is the most likely cause for the continuum and the asymmetry even at room temperature. Thus, we provide an example that shows the power of polarization-resolved Raman scattering in the study of Kitaev materials.

## ACKNOWLEDGEMENTS

Funding for this work was provided by OSU's Center for Emergent Materials, an NSF MRSEC, under award number DMR-1420451. T.T.M., A.M., and A.R.H.W. would like to acknowledge the National Institute of Standards and Technology (NIST)/National Research Council Postdoctoral Research Associateship Program and NIST-STRS (Scientific and Technical Research Services) for funding. P.L.K. and D.M. acknowledge support from the Gordon and Betty Moore Foundations EPiQS Initiative through grant GBMF441. J.Y. was supported by the DOE Office Basic Energy Sciences, Materials Sciences and Engineering Division. S.E.N. was supported by the US DOE, Division of Scientific User Facilities, Office of Basic Energy Sciences. We acknowledge useful discussions with N.P. Armitage, Y-M. Lu and N. Trivedi.

- ∗ thuc.mai@nist.gov
- † valdesaguilar.1@osu.edu
- 1 L. Balents, Nature 464 , 199 EP (2010).
- 2 A. Kitaev, Annals of Physics 321 , 2 (2006), january Special Issue.
- 3 L. Savary and L. Balents, Reports on Progress in Physics 80 , 016502 (2016).
- 4 M. Hermanns, I. Kimchi, and J. Knolle, Annual Review of Condensed Matter Physics 9 , 17 (2018), https://doi.org/10.1146/annurev-conmatphys-033117053934.
- 5 J. Knolle and R. Moessner, Annual Review of Condensed Matter Physics 10 , 451 (2019), https://doi.org/10.1146/annurev-conmatphys-031218013401.
- 6 H. Takagi, T. Takayama, G. Jackeli, G. Khaliullin, and S. E. Nagler, Nature Reviews Physics 1 , 264 (2019).
- 7 A. Banerjee, C. A. Bridges, J.-Q. Yan, A. A. Aczel, L. Li, M. B. Stone, G. E. Granroth, M. D. Lumsden, Y. Yiu, J. Knolle, S. Bhattacharjee, D. L. Kovrizhin, R. Moessner, D. A. Tennant, D. G. Mandrus, and S. E. Nagler, Nature Materials 15 , 733 EP (2016), article.
- 22 S. M. Winter, K. Riedl, P. A. Maksimov, A. L. Chernyshev, A. Honecker, and R. Valent´ ı, Nature communications 8 , 1152 (2017).
- 23 B. Perreault, J. Knolle, N. B. Perkins, and F. J. Burnell, Phys. Rev. B 94 , 104427 (2016).
- 24 H. B. Cao, A. Banerjee, J.-Q. Yan, C. A. Bridges, M. D. Lumsden, D. G. Mandrus, D. A. Tennant, B. C. Chakoumakos, and S. E. Nagler, Phys. Rev. B 93 , 134423 (2016).
- 25 R. D. Johnson, S. C. Williams, A. A. Haghighirad, J. Singleton, V. Zapf, P. Manuel, I. I. Mazin, Y. Li, H. O. Jeschke, R. Valent´ ı, and R. Coldea, Phys. Rev. B 92 , 235119 (2015).
- 26 J. Nasu, J. Knolle, D. L. Kovrizhin, Y. Motome, and R. Moessner, Nature Physics 12 , 912 (2016).
- 27 Y. Wang, G. B. Osterhoudt, Y. Tian, P. Lampen-Kelley, A. Banerjee, T. Goldstein, J. Yan, J. Knolle, J. Nasu, Y. Motome, S. Nagler, D. Mandrus, and K. S. Burch, 'High temperature fermi statistics from majorana fermions in an insulating magnet,' (2018), arXiv:1809.07782.
- 28 R. Yadav, N. A. Bogdanov, V. M. Katukuri, S. Nishimoto, J. Van Den Brink, and L. Hozoi, Scientific reports 6 , 37925 (2016).
- 8 A. Banerjee, J. Yan, J. Knolle, C. A. Bridges, M. B. Stone, M. D. Lumsden, D. G. Mandrus, D. A. Tennant, R. Moessner, and S. E. Nagler, Science 356 , 1055 (2017), http://science.sciencemag.org/content/356/6342/1055.full.pdf.
- 9 A. Banerjee, P. Lampen-Kelley, J. Knolle, C. Balz, A. A. Aczel, B. Winn, Y. Liu, D. Pajerowski, J. Yan, C. A. Bridges, A. T. Savici, B. C. Chakoumakos, M. D. Lumsden, D. A. Tennant, R. Moessner, D. G. Mandrus, and S. E. Nagler, npj Quantum Materials 3 , 8 (2018).
- 10 L. J. Sandilands, Y. Tian, K. W. Plumb, Y.-J. Kim, and K. S. Burch, Phys. Rev. Lett. 114 , 147201 (2015).
- 11 A. Glamazda, P. Lemmens, S.-H. Do, Y. S. Kwon, and K.-Y. Choi, Phys. Rev. B 95 , 174429 (2017).
- 12 A. Glamazda, P. Lemmens, S.-H. Do, Y. S. Choi, and K.-Y. Choi, Nature Communications 7 , 12286 EP (2016).
- 13 J. Knolle, G.-W. Chern, D. L. Kovrizhin, R. Moessner, and N. B. Perkins, Phys. Rev. Lett. 113 , 187201 (2014).
- 14 B. Perreault, J. Knolle, N. B. Perkins, and F. J. Burnell, Phys. Rev. B 92 , 094439 (2015).
- 15 K. Ran, J. Wang, W. Wang, Z.-Y. Dong, X. Ren, S. Bao, S. Li, Z. Ma, Y. Gan, Y. Zhang, J. T. Park, G. Deng, S. Danilkin, S.-L. Yu, J.-X. Li, and J. Wen, Phys. Rev. Lett. 118 , 107203 (2017).
- 16 B. Zhou, Y. Wang, G. B. Osterhoudt, P. Lampen-Kelley, D. Mandrus, R. He, K. S. Burch, and E. A. Henriksen, Journal of Physics and Chemistry of Solids (2018), https://doi.org/10.1016/j.jpcs.2018.01.026.
- 17 I. Rousochatzakis, S. Kourtis, J. Knolle, R. Moessner, and N. B. Perkins, Phys. Rev. B 100 , 045117 (2019).
- 18 G. Jackeli and G. Khaliullin, Phys. Rev. Lett. 102 , 017205 (2009).
- 19 S. M. Winter, Y. Li, H. O. Jeschke, and R. Valent´ ı, Phys. Rev. B 93 , 214431 (2016).
- 20 J. A. Sears, Y. Zhao, Z. Xu, J. W. Lynn, and Y.-J. Kim, Phys. Rev. B 95 , 180411 (2017).
- 21 S.-H. Do, S.-Y. Park, J. Yoshitake, J. Nasu, Y. Motome, Y. Kwon, D. T. Adroja, D. J. Voneshen, K. Kim, T.-H. Jang, J.-H. Park, K.-Y. Choi, and S. Ji, Nature Physics 13 , 1079 (2017).
- 29 A. Little, L. Wu, P. Lampen-Kelley, A. Banerjee, S. Patankar, D. Rees, C. A. Bridges, J.-Q. Yan, D. Mandrus, S. E. Nagler, and J. Orenstein, Phys. Rev. Lett. 119 , 227201 (2017).
- 30 X. Lu, X. Luo, J. Zhang, S. Y. Quek, and Q. Xiong, Nano Research 9 , 3559 (2016).
- 31 R. Birss, Symmetry and magnetism , Selected topics in solid state physics (North-Holland Pub. Co., 1964).
- 32 U. Fano, Phys. Rev. 124 , 1866 (1961).
- 33 S. Mashhadi, D. Weber, L. M. Schoop, A. Schulz, B. V. Lotsch, M. Burghard, and K. Kern, Nano Letters 18 , 3203 (2018), pMID: 29635914, https://doi.org/10.1021/acs.nanolett.8b00926.
- 34 G. Li, X. Chen, Y. Gan, F. Li, M. Yan, F. Ye, S. Pei, Y. Zhang, L. Wang, H. Su, J. Dai, Y. Chen, Y. Shi, X. Wang, L. Zhang, S. Wang, D. Yu, F. Ye, J.-W. Mei, and M. Huang, Phys. Rev. Materials 3 , 023601 (2019).
- 35 D. T. Larson and E. Kaxiras, Phys. Rev. B 98 , 085406 (2018).
- 36 S. Djurdji ´ c Mijin, A. Solaji´, ˇ c J. Peˇ si´, c M. ˇ ´epanovi´, Sc c Y. Liu, A. Baum, C. Petrovic, N. Lazarevi´, c and Z. V. Popovi´ c, Phys. Rev. B 98 , 104307 (2018).
- 37 K. W. Plumb, J. P. Clancy, L. J. Sandilands, V. V. Shankar, Y. F. Hu, K. S. Burch, H.-Y. Kee, and Y.-J. Kim, Phys. Rev. B 90 , 041112 (2014).
- 38 S. Sinn, C. H. Kim, B. H. Kim, K. D. Lee, C. J. Won, J. S. Oh, M. Han, Y. J. Chang, N. Hur, H. Sato, B.-G. Park, C. Kim, H.-D. Kim, and T. W. Noh, Scientific Reports 6 , 39544 (2016).
- 39 L. J. Sandilands, C. H. Sohn, H. J. Park, S. Y. Kim, K. W. Kim, J. A. Sears, Y.-J. Kim, and T. W. Noh, Phys. Rev. B 94 , 195156 (2016).
- 40 S. Reschke, F. Mayr, S. Widmann, H.-A. Krug von Nidda, V. Tsurkan, S.-H. Eremin, M. V. Do, K.-Y. Choi, Z. Wang, and A. Loidl, J. Phys.: Condens. Matter 30 , 475604 (2018).
- 41 Z. Wang, S. Reschke, D. H¨vonen, S.-H. Do, K.-Y. Choi, u M. Gensch, U. Nagel, T. R˜ om, o and A. Loidl, Phys. Rev.

- 42 S. Reschke, F. Mayr, Z. Wang, S.-H. Do, K.-Y. Choi, and A. Loidl, Phys. Rev. B 96 , 165120 (2017).
Thuc T. Mai A. McCreary P. Lampen-Kelley N. Butch J. R. Simpson J. -Q. Yan S. E. Nagler D. Mandrus A. R. Hight Walker R. Valdes Aguilar
Polarization-resolved Raman spectroscopy of α-RuCl3 and evidence of room temperature two-dimensional magnetic scattering
Polarization-resolved Raman spectroscopy was performed and analyzed from
large, high quality, mono-domain single crystal of {\alpha}-RuCl3, a proximate
Kitaev quantum spin liquid. Spectra were collected with laser polarizations
parallel and perpendicular to the honeycomb plane. Pairs of nearly degenerate
phonons were discovered and show either a 4-fold or 2-fold polarization angle
dependence in their Raman intensity, thereby providing evidence to definitively
assign the bulk crystal point group as C2h. The low frequency continuum that is
often attributed to scattering from pairs of Majorana fermions was also
examined and found to disappear when the laser excitation and scattered photon
polarizations were perpendicular to the honeycomb plane. This disappearance,
along with the behavior of the phonon spectrum in the same polarization
configuration, strongly suggests that the scattering continuum is
2-dimensional. We argue that this scattering continuum originates from the
Kitaev magnetic interactions that survives up to room temperature, a scale
larger than the bare Kitaev exchange energy of approximately 50 K.