## Observation of three-state nematicity in the

## triangular lattice antiferromagnet Fe1/3 NbS2

Arielle Little 1,2 , Changmin Lee 2 , Caolan John 1,2 , Spencer Doyle 1,2 , Eran Maniv 1,2 , Nityan L. Nair 1,2 , Wenqin Chen 1,2 , Dylan Rees 1,2 , Jšrn W.F. Venderbos 4,5 , Rafael M. Fernandes 3 , James G. Analytis 1,2 , Joseph Orenstein 1,2*

- 1 Department of Physics, University of California, Berkeley, California 94720, USA 2 Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 3 School of Physics and Astronomy, University of Minnesota, Minneapolis, Minnesota 55455, USA 4 Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104Ð6323, USA 5 Department of Physics and Astronomy, University of Pennsylvania, Philadelphia,
- Pennsylvania 19104Ð6396, USA

These authors contributed equally to this work *Corresponding author (email: jworenstein@lbl.gov)

Nematic  order  is  the  breaking  of  rotational  symmetry  in  the  presence  of translational  invariance.  While  originally  defined  in  the  context  of  liquid crystals, the concept of nematic order has arisen in crystalline matter with discrete rotational symmetry, most prominently in the tetragonal Fe-based superconductors where the parent state is four-fold symmetric. In this case the nematic director takes on only two directions, and the order parameter in  such  ÒIsing-nematicÓ  systems  is  a  simple  scalar.  Here,  using  a  novel spatially-resolved optical polarimetry technique, we show that a qualitatively distinct  nematic  state  arises  in  the  triangular  lattice  antiferromagnet Fe1/3NbS2. The crucial difference is that the nematic order on the triangular lattice is a ℤ", or three-state Potts-nematic order parameter. As a consequence,  the  anisotropy  axes  of  response  functions  such  as  the resistivity tensor can be continuously re-oriented by external perturbations. This discovery provides insight into realizing devices that exploit analogies with nematic liquid crystals.

Currently there is intense focus on incorporating antiferromagnets (AFMs) in spintronic applications, with the promise of faster response, lower threshold current, and scaling to smaller dimensions 1,2 . For example, the observation that stable switching of electric resistance can be induced by a spin-unpolarized current in thin films of antiferromagnetic CuMnAs, Mn2Au, and NiO has attracted considerable attention 3-6 . In closely related work, multi-stable magnetic memory with electric read-out was demonstrated in the hexagonal AFM, MnTe 7,8 . These systems are easy-plane, collinear AFMs in which switching and memory are associated with rotation of the NŽel vector, 𝑳, between stable states in the plane that contains the perturbing magnetic field or current, where 𝑳 ≡ 𝑴𝟏 - 𝑴𝟐, and 𝑴𝟏 , 𝑴𝟐 are the sublattice magnetizations 9,10 .

Recently, both current-induced switching and multistable memory effects were reported in the Fe-intercalated transition metal dichalcogenide (TMD) FexNbS2 with 𝑥 ≈ 1/3 , compounds that undergo a transition to AFM order below approximately 50 K 11 . It was found that current pulses applied parallel to the TMD atomic layers could reversibly switch a 10-micron scale device between stable resistance states with significantly lower threshold current than required for CuMnAs, albeit at lower temperatures. Moreover, it was discovered that cooling through the NŽel temperature ( 𝑇1 ) in an in-plane magnetic field, 𝑩 , induces an in-plane resistivity anisotropy whose symmetry axes continuously follow the direction of 𝑩 and remain stable after the field is turned off. While the anisotropy of the static magnetic susceptibility, 𝜒 , in Fe1/3NbS2 suggests that 𝑳 is oriented primarily perpendicular to the TMD planes 12-14 . It has been suggested that this compound possesses a small in-plane component of 𝑳 , whose rotation by current pulses may give rise to the observed switching phenomena.

Here we report optical measurements which suggest that an in-plane nematic director plays a role in mediating switching and metastable memory in the Fe1/3NbS2 system. Using spatially-resolved optical polarimetry, we show that the onset of AFM order in Fe1/3NbS2 occurs simultaneously with a first order transition that breaks the 6fold (screw) rotational symmetry of the paramagnetic phase. Below 𝑇1 we observe three nematic domains whose directors are rotated by an angle of 2𝜋/3 with respect to each other. We provide a theoretical understanding of these results by showing that the AFM transition in Fe1/3NbS2 is analogous to the magneto-structural transition in Fe-based

superconductors 15-19 , but with a crucial difference Ð the ℤ6 Ising-nematic degree of freedom obtained in a tetragonal structure becomes a ℤ" or three-state Potts-nematic on the triangular lattice 20,21 .

A schematic of the optical set-up is shown in Fig. 1a. A linearly polarized probe beam is focused through a microscope objective at normal incidence onto the sample surface. The reflection amplitude from the surface is characterized by a 2 x 2 matrix, 𝑟 89 . For a material with 3-fold or higher rotational symmetry, 𝑟 89 = 𝛿89 𝑟 , and the polarization of the probe beam remains unchanged upon reflection. Lowering rotational symmetry leads to a reflection matrix with principal optic axes 𝒂 and 𝒃 , and birefringence, 𝛥𝑟 ≡ 𝑟 ? -𝑟@ . In this broken symmetry state, the probe polarization is rotated through an angle given by 𝜙 = 𝑅𝑒{𝑟 ∗ 𝛥𝑟 𝑐𝑜𝑠 2(𝜃 - 𝜃 K )}/|𝑟| 6 , where 𝜃 - 𝜃 K is the difference between the probe beam polarization and principal axis direction. By measuring 𝜙(𝜃) with a balanced optical bridge detector (see Supplementary Information) we obtain the optic axes directions 𝒂 and 𝒃 and birefringence amplitude, 𝛥𝑟 .

To obtain the maps of 𝛥𝑟 and 𝜃K described below we overlap the probe beam with an 800 nm pump laser chopped at 2 kHz, which modulates the sample temperature. Lock-in detection at the pump chopping frequency eliminates long-time-period drifts and enables microradian sensitivity to the modulated polarization rotation, 𝛿𝜙 . Fig. 1b shows 𝛿𝜙 as a function of temperature, 𝑇 , at a single 10-micron spot on the sample and for a fixed probe polarization. The abrupt onset of birefringence at 49 K indicates a firstorder phase transition to a state in which the 6-fold screw symmetry of the paramagnetic state is lowered to at most 2-fold rotational symmetry. The decaying oscillations in 𝛿𝜙 with further lowering of 𝑇 are explained below. Plotted in Fig. 1c is 𝛿𝜙 as a function of probe polarization angle 𝜃 in the vicinity of the phase transition. Below 𝑇1 , 𝛿𝜙 shows the expected cos2(𝜃 - 𝜃K ) dependence, allowing us to determine the orientation of the optic axes in the magnetic phase.

The clue as to the origin of the oscillations in 𝛿𝜙(𝑇) shown in Fig. 1b is that the optical probe detects the onset of order at a higher temperature than bulk probes do. Figure 2a shows a summary of bulk measurements in the temperature range of the transition. A singularity in heat capacity, 𝐶S , coincides with abrupt jumps in 𝜒T and inplane resistivity, 𝜌VV , indicating a single, first-order transition to the AFM state at 43 K.

The oscillatory 𝛿𝜙(𝑇) is plotted in conjunction with 𝐶S in Fig. 2b, showing that there is a large offset between the transition temperatures measured optically (49 K) and by bulk probes (43 K).

Figure 2c shows the polarization rotation 𝜙(𝑇) measured directly, that is, without photo-thermal modulation. The presence of the same oscillations in 𝜙(𝑇) seen in 𝛿𝜙(𝑇) confirms that they are not an artifact of photo-thermal modulation. Instead, our analysis, presented in detail in the Supplementary Information, shows that oscillations appear in the optical measurements because 𝑇1 decreases with increasing depth below the surface. For temperatures below the surface transition temperature 𝑇1 (0) = 49 K, a buried interface separates a birefringent surface layer from the isotropic bulk, as illustrated in Fig. 2d. As the depth, 𝑧 , of the interface increases with decreasing 𝑇 , the phase difference between the reflections from the surface and the buried interface produces the interference pattern seen in 𝜙(𝑇) . The oscillations are cut off when the buried interface reaches beyond the penetration depth of the probe beam, which is approximately 150 nm 22 . At this depth, 𝑇1 (𝑧) has reached 47.5 K, which is still 4.5 K higher than the bulk 𝑇1 . The mesoscopic scale of the probe wavelength in the medium enables a new method for depth profiling of transition temperatures on a 10-nm scale.

We now turn to the spatial mapping of the amplitude and principal axis directions of the nematic order. Below the transition we detect three orientations of optic axes offset by 2𝜋/3 from each other, as illustrated by polar plots of 𝛿𝜙(𝜃) at three locations on the sample (Fig. 3a). By registration of the x-ray Laue diffraction pattern with the probe polarization angle, we find that the three orientations of optic axes correspond to the three crystallographic symmetry directions of the triangular Fe-lattice (Fig. 4b). We assign a color to each of the orientations, such that [100] is red, [010] is green, and [11 [ 0] [ [ [ is blue. A map of a 900 𝜇 m × 500 𝜇 m region of the sample is shown in Fig. 3b, revealing the presence of all three domains, which can be as large as hundreds of microns or as small as our resolution of 10 microns. The domain distribution is deterministic upon warming and cooling, and does not change in magnetic fields up to 400 G.

The nematic order that accompanies the onset of AFM order in Fe1/3NbS2 can be qualitatively understood to result from the geometric frustration of Ising spins 23-25  on

the triangular ( √3𝑎 × √3𝑎) superlattice 26-29  of Fe atoms. A distortion along one direction of Fe-Fe bonds relieves the frustration, allowing one of three degenerate ordered phases to condense. Both stripe-like order (ordering wavevectors 𝑸8 in the 𝛤 - 𝑀 directions) and zigzag order ( 𝑸8 in the 𝛤 - 𝐾 directions) are possible (Figs. 4b,c), and single-crystal neutron measurements are necessary to fully distinguish between the two 12 . Both types of order give rise to nematicity and the stripe case is discussed here (see Supplementary Information for zigzag).

To gain further insight into the nature of the transition, we construct a phenomenological Landau model assuming stripe order. We introduce three order parameters 𝑳d , 𝑳6 , and 𝑳" , whose wavevectors are parallel to one of the three 𝛤 - 𝑀 directions. Spin-rotational and 6-fold rotational symmetries restrict the corresponding magnetic Landau free-energy 𝐹f to the form,

$$F _ { M } = a \sum _ { \iota = 1, 3 } L _ { i } ^ { 2 } + \nu _ { 0 } \left ( \sum _ { \iota = 1, 3 } L _ { i } ^ { 2 } \right ) ^ { 2 } + \nu _ { 1 } \sum _ { \iota < j } L _ { i } ^ { 2 } L _ { j } ^ { 2 } + \nu _ { 2 } \sum _ { \iota < j } \left ( L _ { i } \cdot L _ { j } \right ) ^ { 2 }.$$

For simplicity, here we assume commensurate wave-vectors (Fig. 4c), but the results can be extended in a straightforward way to incommensurate order. The Landau coefficients 𝑣d and 𝑣6 determine the nature of the magnetic state below the NŽel temperature. In particular, for 𝑣d &gt; 0 and 𝑣6 &gt; -𝑣d a single𝑄 magnetic state is favored in which only one of the 𝑳8 is nonzero and rotational symmetry is broken. In contrast, other parameter values lead to equal amplitude triple𝑄 magnetic ordering, preserving the rotational symmetry of the paramagnetic state. While the particular values of these and the other Landau coefficients a and 𝑣K depend on microscopic considerations related to the mechanism responsible for the magnetic instability, for our purposes a phenomenological approach suffices.

Indeed, our experimental observation of rotational symmetry breaking points uniquely to singleQ AFM order in Fe1/3NbS2 (i.e. 𝑣d &gt; 0 and 𝑣6 &gt; -𝑣d ), in which 〈𝑳8 〉 ≠ 0 for one of the three order parameter components. As seen in Fig. 4d, the six nearestneighbor links are no longer equivalent in the stripe AFM state; four bonds couple antiparallel spins whereas two bonds couple parallel spins. This rotational symmetry

breaking is captured by an order parameter 𝒏 = (𝑛d, 𝑛 6 ) , which is given in terms of the magnetic components as,

$$( n _ { 1 }, n _ { 2 } ) = ( L _ { 1 } ^ { 2 } + L _ { 2 } ^ { 2 } - 2 L _ { 3 } ^ { 2 }, \sqrt { 3 } L _ { 1 } ^ { 2 } - \sqrt { 3 } L _ { 2 } ^ { 2 } ).$$

Here 𝒏 is a nematic director indicative of rotational symmetry breaking and can be parametrized as 𝒏 = 𝑛(cos2𝜃, sin 2𝜃) . The structure of 𝒏 for the triangular lattice contrasts sharply with nematic order in tetragonal systems such as the Fe-pnictides, which is described by a single-component ℤ6 Ising-nematic order parameter 15,18 .

The magnetic free energy (Eq. 1) can be used to derive an effective free energy Fn for nematic order, since the nematic order parameter (Eq. 2) is a composite magnetic order parameter. Integrating out magnetic fluctuations and going beyond mean-field theory, one finds,

$$F _ { n } = \alpha n ^ { 2 } + \beta n ^ { 3 } \, c o s \, 6 \theta + \gamma n ^ { 4 },$$

where the Landau parameters 𝛼 𝛽 , , and 𝛾 depend on the Landau parameters of the magnetic free energy, Eq. 1. The key observation is that Eq. 3 represents the same free energy as that of the three-state Potts model 21,30 , implying that 𝒏 is a ℤ" Potts-nematic order parameter. This follows from the third-order term in Eq. 3, which defines a ℤ" Potts-anisotropy and restricts the director to point along one of the three highsymmetry directions of the lattice. The first-order jump in the nematic order parameter arising from the cubic term triggers a jump in the magnetic correlation length and a simultaneous first-order nematic-AFM transition consistent with the data.

Next we show that the relative population of the three AFM domains can be tuned in response to an external perturbation that couples to the nematic director, in this case uniaxial strain. The birefringence map shown in Fig. 3b was obtained with the sample resting on a Cu plate secured by vacuum grease, a configuration in which the external strain from thermal contraction is negligibly small. The inset presents a histogram illustrating the domain population measured in this configuration, showing that all three domains are represented.

To apply strain the sample was glued onto a piezoelectric stack, which upon cooling applies a uniaxial strain of ≈ 0.1 % 31-33 . Imaging the same region of the sample with compressive strain applied parallel to the bond direction, as illustrated in Fig. 5a, yields the birefringence map shown in Fig. 5b. This map, together with the domain histogram (Fig. 5c), shows that the areal fraction of the red [100] domains, whose principal axes are parallel and perpendicular to the bond direction, is strongly suppressed and redistributed to green and blue domains. The sample was then cleaved, rotated by 90 ∘ and re-mounted on the piezo stack such that now tensile strain is applied in the bond direction (Fig. 5d). In this configuration the areal fraction of the red domains grows at the expense of the other two (Figs. 5e,f). In the Supplementary Information we show that the domain repopulation follows in a straightforward fashion from the coupling of the strain tensor to the nematic director, 𝒏 . Basically, tensile strain acts as a positive nematic conjugate field oriented parallel to the strain direction, favoring the domain corresponding to this direction. On the other hand, compressive strain acts as a negative conjugate field, which suppresses the corresponding domain and favors the other two domains.

The repopulation of ℤ" nematic domains in response to uniaxial strain suggests a mechanism for switching and memory phenomena in response to in-plane perturbations 𝑩 and currents, 𝑱. In Fe1/3NbS2, 𝒏 can couple to such perturbing fields, playing the same role as the in-plane 𝑳 does in the easy-plane systems such as CuMnAs and MnTe. As shown for coupling to 𝑳 in MnTe 6 , and for coupling to 𝒏 in Fe1/3NbS2 (see Supplementary Information), application of in-plane 𝑩 or 𝑱 unbalances the domain population, resulting in a global resistivity tensor whose principal axes continuously follow 34  the direction of 𝑩 or 𝑱 . We note that this effect cannot occur in a ℤ6 Isingnematic system, where the anisotropy is locked to the crystal axes, regardless of the direction of the perturbation 15,32 .

In conclusion, we have shown that the onset of AFM order in Fe1/3NbS2 occurs via a first-order transition that lowers the rotational symmetry of the triangular lattice from 6 to at most 2-fold. Below TN , maps of local birefringence reveal three domains, characterized by in-plane nematic directors with relative angles 2𝜋/3 . The first-order nature of this transition may be understood in terms of a three-state Potts model where

the symmetry-breaking order parameter is a ℤ" nematic, which in turn is a composite magnetic order parameter. Furthermore, we demonstrated that the relative population of ℤ" nematic domains can be controlled by uniaxial strain. In principle, the coupling of external perturbations to a nematic director could be sufficiently strong to enable low switching thresholds in nematic-magnetic systems.

## Methods

Imaging the direction of the optic axes is performed by mounting the sample on a cryogenic piezo-actuated xyz stage in a Montana Instruments 4 K cryostat with a lowworking-distance cryogenic window. As the sample is rastered beneath the probe beam focus, a polarization scan of 𝛿𝜙 is taken at each position on a grid. Polarization scanning is enabled by co-rotating half wave plates, a scheme described in more detail in the Supplementary Information. Whereas dichroism-based imaging 35,36  is limited to the contrast between two fixed orthogonal polarizations, here the continuous control of the probe polarization measures the precise local orientation of the optic axes.

## Acknowledgements:

We thank D.H. Lee and J. Moore for useful discussions and N. Tamura for support at the Advanced Light Source. Optical measurements were performed at the Lawrence Berkeley National Laboratory in the Quantum Materials program supported by the supported by the Director, Office of Science, Office of Basic Energy Sciences, Materials Sciences and Engineering Division, of the U.S. Department of Energy under Contract No. DE-AC02- 05CH11231. A.L. and J.O. received support for optical measurements from the Gordon and Betty Moore FoundationÕs EPiQS Initiative through Grant No. GBMF4537 to J. O. at UC Berkeley.  Synthesis of Fe1/3NbS2 was supported by Laboratory Directed Research and Development Program of Lawrence Berkeley National Laboratory under Contract No. DE-AC02-05CH11231 . J.G.A. and N.L.N. received support from the Gordon and Betty Moore FoundationÕs EPiQS Initiative Grant No. GBMF4374 to J. G.A. at UC Berkeley. R.M.F. is supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, under Award DE-SC0012336. X-ray diffraction to register crystal orientation was carried out at beam line 12.3.2 at the

Advanced Light Source, which is a Department of Energy User Facility under Contract No. DE-AC02-05CH11231.

## Author contributions:

A.L. and C.L. performed and contributed equally to the birefringence microscopy measurements and data analysis. C.J., S.D. and E.M. grew and characterized the crystals. E.M., N.L.N., and J.G.A. discovered the switching effect which motivated this project. J.W.F.V. and R.M.F. developed the theoretical model. W.C. performed the simulation of depth profiling. D.R. assisted with optical measurements. J.O., A.L., C.L., J.W.F.V., and R.M.F. wrote the manuscript. All authors commented on the manuscript.

## Competing Interests:

The authors declare no competing financial interests.

## References

- 1 Jungwirth, T., Marti, X., Wadley, P. &amp; Wunderlich, J. Antiferromagnetic spintronics. Nat. Nanotechnol. 11 , 231-241, (2016).
- 2 Baltz, V. et al. Antiferromagnetic spintronics. Rev. Mod. Phys. 90 , 015005, (2018).
- 3 Wadley, P. et al. Electrical switching of an antiferromagnet. Science 351 , 587590, (2016).
- 4 Bodnar, S. Y. et al. Writing and reading antiferromagnetic Mn2Au by NŽel spinorbit torques and large anisotropic magnetoresistance. Nat. Commun. 9 , 348, (2018).
- 5 Moriyama, T., Oda, K., Ohkochi, T., Kimata, M. &amp; Ono, T. Spin torque control of antiferromagnetic moments in NiO. Sci. Rep. 8 , 14167, (2018).
- 6 Saidl, V. et al. Optical determination of the NŽel vector in a CuMnAs thin-film antiferromagnet. Nat. Photon. 11 , 91, (2017).
- 7 Kriegner, D. et al. Multiple-stable anisotropic magnetoresistance memory in antiferromagnetic MnTe. Nat. Commun. 7 , 11623, (2016).
- 8 Kriegner, D. et al. Magnetic anisotropy in antiferromagnetic hexagonal MnTe. Phys. Rev. B 96 , 214418, (2017).

|   9 | Ž elezn ý , J. et al. Relativistic NŽel-Order Fields Induced by Electrical Current in Antiferromagnets. Phys. Rev. Lett. 113 , 157201, (2014).                                                                             |
|-----|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
|  10 | Ž elezn ý , J. et al. Spin-orbit torques in locally and globally noncentrosymmetric crystals: Antiferromagnets and ferromagnets. Phys. Rev. B 95 , 014403, (2017).                                                         |
|  11 | Nair, N. L. et al. Electrical switching in a magnetically intercalated transition metal dichalcogenide, Preprint at: https://arxiv.org/abs/1907.11698> (2019).                                                             |
|  12 | Van Laar, B., Rietveld, H. M. &Ijdo, D. J. W. Magnetic and crystallographic structures of Me x NbS 2 and Me x TaS 2 . J. Solid State Chem. 3 , 154-160, (1971).                                                            |
|  13 | Gorochov, O., Blanc-soreau, A. L., Rouxel, J., Imbert, P. & Jehanno, G. Transport properties, magnetic susceptibility and Mšssbauer spectroscopy of Fe 0.25 NbS 2 and Fe 0.33 NbS 2 . Philos. Mag. B 43 , 621-634, (1981). |
|  14 | Yamamura, Y. et al. Heat capacity and phase transition of Fe x NbS 2 at low temperature. J. Alloys Compd. 383 , 338-341, (2004).                                                                                           |
|  15 | Chu, J.-H. et al. In-Plane Resistivity Anisotropy in an Underdoped Iron Arsenide Superconductor. Science 329 , 824-826, (2010).                                                                                            |
|  16 | Johnston, D. C. The puzzle of high temperature superconductivity in layered iron pnictides and chalcogenides. Adv. Phys. 59 , 803-1061, (2010).                                                                            |
|  17 | Paglione, J. & Greene, R. L. High-temperature superconductivity in iron-based materials. Nat. Phys. 6 , 645, (2010).                                                                                                       |
|  18 | Fernandes, R. M., Chubukov, A. V. & Schmalian, J. What drives nematic order in iron-based superconductors? Nat. Phys. 10 , 97, (2014).                                                                                     |
|  19 | Si, Q., Yu, R. &Abrahams, E. High-temperature superconductivity in iron pnictides and chalcogenides. Nat. Rev. Mater. 1 , 16017, (2016).                                                                                   |
|  20 | Fernandes, R. M., Orth, P. P. & Schmalian, J. Intertwined Vestigial Order in Quantum Materials: Nematicity and Beyond. Annu. Rev. Condens. Matter Phys. 10 , 133-154, (2019).                                              |
|  21 | Hecker, M. & Schmalian, J. Vestigial nematic order and superconductivity in the doped topological insulator Cu x Bi 2 Se 3 . npj Quantum Mater. 3 , 26, (2018).                                                            |
|  22 | Fan, S. et al. Electronic chirality in the metallic ferromagnet Fe 1/3 TaS 2 . Phys. Rev. B 96 , 205119, (2017).                                                                                                           |
|  23 | Wannier, G. H. Antiferromagnetism. The Triangular Ising Net. Phys. Rev. 79 , 357-364, (1950).                                                                                                                              |

|   24 | Korshunov, S. E. Nature of phase transitions in the striped phase of a triangular- lattice Ising antiferromagnet. Phys. Rev. B 72 , 144417, (2005).                                                        |
|------|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
|   25 | Smerald, A., Korshunov, S. & Mila, F. Topological Aspects of Symmetry Breaking in Triangular-Lattice Ising Antiferromagnets. Phys. Rev. Lett. 116 , 197201, (2016).                                        |
|   26 | Friend, R. H., Beal, A. R. & Yoffe, A. D. Electrical and magnetic properties of some first row transition metal intercalates of niobium disulphide. Philos. Mag. A 35 , 1269-1287, (1977).                 |
|   27 | Parkin, S. S. P. & Friend, R. H. 3d transition-metal intercalates of the niobium and tantalum dichalcogenides. I. Magnetic properties. Philos. Mag. B 41 , 65-93, (1980).                                  |
|   28 | Horibe, Y. et al. Color theorems, chiral domain topology, and magnetic properties of Fe(x)TaS2. J AmChem Soc 136 , 8368-8373, (2014).                                                                      |
|   29 | Doyle, S. et al. Tunable Giant Exchange Bias in an Intercalated Transition Metal Dichalcogenide, Preprint at: https://arxiv.org/abs/1904.05872 (2019).                                                     |
|   30 | Straley, J. P. & Fisher, M. E. Three-state Potts model and anomalous tricritical points. J. Phys. A 6 , 1310-1326, (1973).                                                                                 |
|   31 | Simpson, A. M. & Wolfs, W. Thermal expansion and piezoelectric response of PZT Channel 5800 for use in low - temperature scanning tunneling microscope designs. Rev. Sci. Instrum. 58 , 2193-2195, (1987). |
|   32 | Chu, J.-H., Kuo, H.-H., Analytis, J. G. & Fisher, I. R. Divergent Nematic Susceptibility in an Iron Arsenide Superconductor. Science 337 , 710-712, (2012).                                                |
|   33 | Hicks, C. W., Barber, M. E., Edkins, S. D., Brodsky, D. O. & Mackenzie, A. P. Piezoelectric-based apparatus for strain tuning. Rev. Sci. Instrum. 85 , 065003, (2014).                                     |
|   34 | Fradkin, E., Kivelson, S. A., Lawler, M. J., Eisenstein, J. P. & Mackenzie, A. P. Nematic Fermi Fluids in Condensed Matter Physics. Annu. Rev. Condens. Matter Phys. 1 , 153-178, (2010).                  |
|   35 | Sapozhnik, A. A. et al. Direct imaging of antiferromagnetic domains in Mn 2 Au manipulated by high magnetic fields. Phys. Rev. B 97 , 134429, (2018).                                                      |
|   36 | Bodnar, S. Y. et al. Imaging of current induced NŽel vector switching in antiferromagnetic Mn 2 Au. Phys. Rev. B 99 , 140409, (2019).                                                                      |

Fig. 1 | Photo-thermal modulated birefringence measurements. a, Pump and probe

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00657v1-figure-6077cf7c-a538-4a1e-9cfa-8ed6f743cfa6.jpg)

beams are spatially overlapped onto the sample surface. The 633 nm beam probes the thermally modulated polarization rotation induced by 24 𝜇 W of 800 nm pump that is modulated at 2 kHz with an optical chopper. The dependence on probe beam polarization is measured using a pair of co-rotating half-wave plates. b , Photo-thermally induced polarization rotation 𝛿𝜙 exhibits a sharp peak at 𝑇1 followed by decaying oscillations. c , Polarization rotation 𝛿𝜙 as a function of input probe polarization angle at temperatures close to 𝑇1 . The input polarization angle with maximal 𝛿𝜙 indicates the direction of the optic fast axis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00657v1-figure-a639150b-50f9-402c-9174-50bc4fb4233b.jpg)

Temperature (K)

Fig. 2 | First-order phase transitions at surface and bulk. a, Heat capacity (red), resistivity (green), and magnetic susceptibility (blue) measurements all reveal a sharp first-order phase transition at 𝑇1 = 43 K. b , Optical birefringence data (black) shows a higher 𝑇1 (= 49 K) compared to heat capacity (red) measurement. c, Unmodulated birefringence data also exhibits oscillations at temperatures slightly below the surface 𝑇1 . d, Illustration of surface and bulk in two temperature regimes, showing the presence of a buried interface when 𝑇 is between the surface and bulk transition temperatures.

Data

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00657v1-figure-faafd8e7-8536-41f7-b093-9273b1100c94.jpg)

Fig. 3 | Crystal structure and birefringence map. a, Polar plots of 𝛿𝜙 at different sample locations reveal breaking of rotational symmetry. b , Birefringence map across a 900 𝜇 m × 500 𝜇 m area of the sample shows three distinct domains whose optic fast axes are offset by 120¡ from one another. This region was sampled in 50 um steps, and the positions of the foci are indicated by the white arrows. The three nematic domains in c are color-coded as red, green, and blue. Each is associated with a high-symmetry direction of the Fe-Fe triangular superlattice, as shown in Fig. 4b.  Scale bar: 100 µm. Inset: histogram showing the distribution of domain population in the birefringence map.

Fig 4| Crystal structure and nematic order. a, Cut along the c-direction of Fe1/3NbS2 illustrating staggered intercalation of Fe-ions between the TMD planes. b, View of triangular Fesuperlattice in the ab plane. Colored arrows indicate the high-symmetry directions and three domain orientations shown in Fig. 3. c, Hexagonal Brillouin zone for a single Fe-layer with highsymmetry directions shown. The three Γ - 𝑀 ordering wave-vectors ( Q 1, Q 2, Q 3), related by 2 𝜋 /3 rotations, corresponding to the three domains. d, Real space cartoon of Q 1 stripe order. Solid black circles indicate spin up, white circles with black XÕs indicate spin down. Light blue bonds connect parallel spins, pink bonds connect anti-parallel spins.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00657v1-figure-a47976a8-4b75-4862-b52b-badcf0ed70de.jpg)

Fig. 5 | Strain tuning of birefringent domains. a, Compressive strain is applied along the [100] direction of the lattice, parallel to the Fe-Fe bond direction. b, The birefringence map and c, histogram showing that the green and blue domains are favored over the red. d, Conversely, tensile strain along the [100] direction results in dominance of red domains, as illustrated by the (e) birefringence map and (f) histogram.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00657v1-figure-ae2096b5-96cc-407e-980f-6278a0ec5ff5.jpg)

## Supplementary Information

ÒObservation of three-state nematicity in the triangular lattice antiferromagnet Fe1/3NbS2Ó

## I. Sample Growth

Single crystals of Fe1/3NbS2 were grown using vapor transport techniques. 1  The samples studied are hexagonal crystals 1-3 mm in size and 20 Ð 100 microns thick. Samples were grown with a nominal Fe-concentration of x = 0.37. Energy-dispersive xray spectroscopy (EDX) and inductively coupled plasma (ICP) measurements were performed to verify the chemical composition, yielding an Fe-intercalation value of x = 0.34. X-ray diffraction and transmission electron microscopy measurements 1  confirm the crystal structure is the nominal x =1/3 non-centrosymmetric space group P 6322.

## II. Birefringence Microscopy and Balanced Optical Bridge Detector

Fig. S1 | Schematic of photo-thermal modulated birefringence microscopy and optical bridge detector.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00657v1-figure-cb262e04-86f8-4917-8ca1-73dfaab67b77.jpg)

A schematic of the measurement apparatus is shown in Fig. S1. The polarization of the incident probe laser is set by a rotatable half-wave plate, HWP1. For the photothermal modulation measurements, the pump is modulated at 2 KHz by an optical chopper and both beams are focused onto the sample surface with a microscope objective. The measured 𝛿𝜙(𝑇) is independent of the pump polarization. The pump reflection is rejected via a color filter while the reflected probe is sent through a Wollaston prism which spatially separates the horizontal ( x ) and vertical ( y ) polarization components of the beam. The two orthogonally polarized beams are sent to an optical bridge detector consisting of two unbiased photodiodes connected in parallel but with opposite polarity. The net photocurrent will be zero, or ÔbalancedÕ, when the x and y components have equal intensity. Balancing is achieved by warming the sample above the NŽel temperature and adjusting a second half-wave plate (HWP2).

In the antiferromagnetic (AFM) phase the sample becomes birefringent, rotating the polarization of the reflected probe and unbalancing the detector. The intensity admitted to each of the two photodiodes can be calculated using Jones calculus,

$$\left [ \begin{matrix} V _ { x } \\ V _ { y } \end{matrix} \right ] = \left [ \begin{matrix} \cos \theta & - \sin \theta \\ \sin \theta & \cos \theta \end{matrix} \right ] \left [ \begin{matrix} r _ { a } & 0 \\ 0 & r _ { b } \end{matrix} \right ] \left [ \begin{matrix} \cos \theta & \sin \theta \\ - \sin \theta & \cos \theta \end{matrix} \right ] \left [ \begin{matrix} 1 \\ 1 \end{matrix} \right ],$$

where 𝜃 , set by the two half-wave plates, represents the angle between the incident probe beam polarization and the sample principal optic axes, 𝒂 and 𝒃 . The reflection coefficients along each of the sample principle axes are written as 𝑟 @ and 𝑟 ? . The balanced diode detector measures the intensity difference between the two polarization components, |𝑉 V | 6 -GLYPH&lt;140&gt;𝑉 GLYPH&lt;136&gt; GLYPH&lt;140&gt; 6 . To first order in Δ𝑟/𝑟̅ , where,

$$\bar { r } = \frac { r _ { a } + r _ { b } } { 2 }, \ \Delta r = r _ { a } - r _ { b },$$

we obtain,

$$| V _ { 1 } | ^ { 2 } - | V _ { 2 } | ^ { 2 } = 2 R e ( \bar { r } ^ { * } \Delta r ) \cos 2 \theta.$$

To measure the unmodulated polarization rotation, 𝜙(𝑇) , for a fixed HWP1 angle as in the main text Fig. 2c, the pump is removed entirely. The optical chopper is moved to modulated the probe beam and detection method remains the same.

## III. Depth Profiling of the AFM Transition

As discussed in the main text, the AFM phase sets in at the sample surface at a higher temperature than in the bulk. The temperature oscillations in 𝜙(𝑇) shown in the main text Fig 2(b), (c) are well-described by a model wherein a birefringent surface layer propagates into the bulk with decreasing temperature.

Fig. S2 | a, Illustration of the sharp, shifting interface between the AFM phase, where Δ𝑛 = 0.1 , and the paramagnetic phase, where Δ𝑛 = 0 . b, The result of simulation. The ratio between the first maximum point and the first minimum point is approximately 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00657v1-figure-62c9697d-5eb5-4c93-9f4e-c7c95779b499.jpg)

We first consider a sharp interface between the AFM region and the paramagnetic region, as illustrated in Fig. S2a. In the paramagnetic region we assume an isotropic index of refraction given by = 2.6 + 0.6𝑖 , the value measured for the similar compound Fe1/4TaS2 by Fan et al 2 . The AFM surface region has a birefringence characterized by the difference in index of refraction along the optic axes, Δ𝑛 = 𝑛@ - 𝑛? . The simulation result shown in Fig. S2b assumes a relative index birefringence of 0.04 and confirms the sharp interface model provides a good qualitative explanation of the data. However, for this simple case the oscillation amplitude is larger than in the data shown in main text Fig. 2c. In the measurement, the ratio between the first maximum point and the first minimum point is approximately 1.5, but this sharp interface calculation gives a ratio of 2.

This discrepancy may be remedied by considering a smooth interface with a finite width. We model the smooth interface by discretizing the interface into multiple layers, or ÔslabsÕ, with varying Δ𝑛 , as illustrated in Fig. S3a. We use a modified Fermi-Dirac distribution function to describe the broadened interface,

$$\Delta n ( z ) = \frac { \Delta n _ { 1 } } { e ^ { \frac { z - z _ { 0 } } { \alpha } } + 1 }.$$

with three parameters: Δ𝑛d , the difference between refractive indices along two principal optical axes in the ordered phase; 𝑧K , the depth of interface as a function of temperature; and 𝛼 , the interface width.

Our simulation is carried out using the scattering matrix method 3 , which allows us to calculate the transmitted and reflected electric field through a series of birefringent slabs. First, we calculate the 𝐼 8 matrices that describe transmitted electric field amplitude at the surface of a given slab and the 𝑇8 matrices that describe the phase shift of the electric field as it is transmitted through the slab,

$$I _ { i } & = \begin{bmatrix} \frac { 1 } { t _ { i } } & \frac { r _ { i } } { t _ { i } } \\ r _ { i } & \frac { 1 } { t _ { i } } \end{bmatrix}, \quad & T _ { i } & = \begin{bmatrix} e ^ { i \delta _ { i } } & 0 \\ 0 & e ^ { - i \delta _ { i } } \end{bmatrix}, \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \ \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \	 & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad &\quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad \\ \alpha & \quad \cdot & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad &
 & \quad \cdot & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & & \quad & \quad \\ \alpha & \quad \cdot & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & & \quad & \quad & \quad & \quad & \quad \\ \alpha & \quad \cdot & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad \\ \alpha & \quad \cdot & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad \\ \alpha & \quad \cdot & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad &
 & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \ \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad \\ \alpha & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \ \quad & \ \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \ \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \ \quad & \quad & \ \quad & \quad & \quad & \ \quad & \quad & \ \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \ \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \ \quad & \quad & \quad & \ \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \quad & \ \quad & \quad & \ \quad & \quad & \quad & \quad & \ \quad & \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \ & & \ \quad & \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \ \\ \alpha & \quad & \quad & \quad & \ \quad & \ \quad & \quad & \ \quad & \ \quad & \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \ & & \ \quad & \ \quad & \ \ & & \ \ & & \ \ \\ \alpha & \ \quad & \ \quad & \ \ & \ \ } & \ \quad & \ \ & & \ \ \\ \alpha & \quad & \quad & \quad & \quad & \quad & \ \quad & \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \quad & \ \ & \quad & \ \quad & \ \ \ & & \ \ \\ \alpha & \quad & \ \quad & \ \ & \ \ \\ \alpha & \quad & \quad & \quad & \ \quad & \ \quad & \quad & \ \quad & \ \quad & \ \quad & \ \ & \ \ & \ \ & & \ \ & & \ \ & & \ \ \\ \alpha & \quad & \quad & \quad & \quad & \ \ & \quad & \ \ \ & & \ \ & \ \ & & \ \ & & \ \quad & \ \ & & \ \ & \ \ \ & & \ \ & & \ \ \\ \alpha & \ \ & \ \ & \ \ & \ \ & & \ \ & \ \ & \ \ & \ \ & \ \ & & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & & \ \ & \ \ & & \ \ & \ \ & \ \ & & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & & \ \ & \ \ & \ \ & & \ \ & & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & & \ \ & \ \ & \ \ & & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & & \ \ & \ \ & \ \ & \ \ & & \ \ & \ \ & \ \ & & \ \ & & \ \ & & \ \ & & \ \ & & \ \ & \ \ & \ \ & \ \ & & \ \ & \ \ & \ \ & & \ \ & \ \ & & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & & \ \ & \ \ & \ \ & \ \ & \ \ & \ \ & & \ \ & \ \ & \ & \ \ & \ & \ \ & \ & \ \ & \ & \ \ & \ \ & \ \ & & \ \ & & \ \ & & \ \ & \ \ & & \ \ & \ \ & \ \ & & \ \ & & \ \ & & \ \ & \ \ & & \ \ & & \ \ & & \ \ & \ \ & & \$$

where 𝑟 8 is the reflection coefficient and 𝑡 8 is the transmission coefficient for the i th slab. The phase shift over a propagation length Δ𝑧8 is given by 𝛿8 = 2𝜋𝑛8𝛥𝑧8/𝜆 . Second, we calculate the scattering matrix by multiplying all the 𝐼 8 and 𝑇8 matrices in sequence,

$$M = I _ { o } T _ { 1 } I _ { 1 } T _ { 2 } I _ { 2 } \dots I _ { m - 1 } T _ { m } I _ { m }.$$

We can now calculate 𝐸 £/ and 𝐸 ¥ƒ§ at the sample surface separately for the polarization components along the principle optic axes,

$$\left [ \begin{matrix} E _ { a } ^ { \text{in} } \\ E _ { a } ^ { \text{out} } \end{matrix} \right ] = M _ { a } \left [ \begin{matrix} 1 \\ 0 \end{matrix} \right ], \quad \left [ \begin{matrix} E _ { b } ^ { \text{in} } \\ E _ { b } ^ { \text{out} } \end{matrix} \right ] = M _ { b } \left [ \begin{matrix} 1 \\ 0 \end{matrix} \right ].$$

Finally, we can calculate the detector output, 𝑅𝑒(𝑟̅ ∗ Δ𝑟) ,

$$r _ { a } = \frac { E _ { a } ^ { o u t } } { E _ { a } ^ { i n } } \, \quad r _ { b } = \frac { E _ { b } ^ { o u t } } { E _ { b } ^ { i n } } \,,$$

$$R e ( \bar { r } ^ { * } \Delta r ) = R e \left [ \left ( \frac { r _ { a } + r _ { b } } { 2 } \right ) ^ { * } ( r _ { a } - r _ { b } ) \right ].$$

Fig. S3 | a, Illustration of the Fermi-Dirac-like distribution of Δ𝑛 and multilayer slabs. The gradient of blue shows the smooth transition from the birefringent region to isotropic region. b, Fitting parameter 𝛼 determines how broad the interface is. c, The result of calculation for different values of 𝛼 . There is a significant decrease in the oscillation amplitude in birefringence if the interface is broader.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00657v1-figure-f773459c-c633-4f52-99ea-267f21ed5b07.jpg)

The result of the calculation shown in Fig. S3(c) for different values of 𝛼 . If the interface is broader, i.e. 𝛼 is larger, there is a significant overall decrease in amplitude, indicating that a broadened distribution of Δ𝑛 at the interface provides a better fit to the data. This broadened-interface model in combination with the experimental data demonstrates a novel technique for profiling depth dependence of the transition temperature.

## IV. Derivation of the nematic free energy

We start by discussing the symmetry properties of Fe1/3NbS2. The point group embedded in its space group ( 𝑃6"22 ) is 𝐷GLYPH&lt;176&gt; , which is generated by a six-fold rotation 𝐶GLYPH&lt;176&gt;GLYPH&lt;144&gt; and a perpendicular twofold rotation 𝐶6V . In space group 𝑃6"22 the sixfold rotation becomes a six-fold screw 𝑆GLYPH&lt;176&gt;GLYPH&lt;144&gt; = {𝐶GLYPH&lt;176&gt;GLYPH&lt;144&gt; |00 d 6 } . As a result, the space group is generated by 𝑆GLYPH&lt;176&gt;GLYPH&lt;144&gt; and {𝐶 6V |000} . Nematic order lowers the six-fold rotation 𝐶GLYPH&lt;176&gt;GLYPH&lt;144&gt; to twofold 𝐶6GLYPH&lt;144&gt; .

Since the magnetic ordering originates from the Fe sites, we focus on the Fe layers of the crystal structure. Two inequivalent Fe layers should be distinguished, which are exchanged under the screw 𝑆GLYPH&lt;176&gt;GLYPH&lt;144&gt; . If we label the two layers 𝐴 and 𝐵 then we obtain an effective honeycomb geometry, where the 𝐴, 𝐵 layers play the role of the honeycomb sublattices. Introducing a layer index α = 𝐴, 𝐵 a general spin configuration can be expressed in terms of Fourier components as

$$S _ { \alpha } ( { \mathbf R } ) = \sum _ { q } S _ { \alpha } ( { \mathbf q } ) \, e ^ { i { \mathbf q } \cdot { \mathbf R } }.$$

Here we consider magnetic ordering with three ordering vectors 𝐐8hd,6," related by 2𝜋/3 rotations. In the case of stripe order, they correspond to the three Γ - 𝑀 directions of the hexagonal Brillouin zone:

$$\mathbf Q _ { 1 } = \frac { \pi } { \sqrt { 3 } } ( \sqrt { 3 }, 1 ), \ \mathbf Q _ { 2 } = \frac { \pi } { \sqrt { 3 } } ( - \sqrt { 3 }, 1 ), \ \mathbf Q _ { 3 } = \frac { \pi } { \sqrt { 3 } } ( 0, - 2 ).$$

Magnetic ordering at these ordering vectors implies that only 𝐒GLYPH&lt;147&gt; (𝐐8 ) and 𝐒GLYPH&lt;147&gt; (-𝐐8 ) are nonzero. Since 𝐐8 = -𝐐8 modulo a reciprocal lattice vector, we only need to consider 𝐒GLYPH&lt;147&gt; (𝐐8 ) .

Clearly, 𝐒' (𝐐 8 ) and 𝐒' (𝐐 8 ) transform into each other under the sixfold screw 𝑆GLYPH&lt;176&gt;GLYPH&lt;144&gt; . We can form irreducible magnetic order parameters by taking the sum and difference, which are respectively even and odd under the screw rotation. For concreteness, in what follows we focus on the even magnetic order parameters 𝐋8 :

$$\mathfrak { L } _ { i } = S _ { A } ( \mathbf Q _ { i } ) + S _ { B } ( \mathbf Q _ { i } ),$$

but the results are the same also for the odd case. Symmetry constrains the magnetic free energy density to:

$$F = a \sum _ { i } \mathbf L _ { i } ^ { 2 } & + \frac { u } { 2 } \left ( \sum _ { i } \mathbf L _ { i } ^ { 2 } \right ) ^ { 2 } - \frac { g } { 2 } [ ( \mathbf L _ { 1 } ^ { 2 } + \mathbf L _ { 2 } ^ { 2 } - 2 \mathbf L _ { 3 } ^ { 2 } ) ^ { 2 } ] \\ & + \frac { w } { 2 } [ ( \mathbf L _ { 1 } \cdot \mathbf L _ { 2 } ) ^ { 2 } + ( \mathbf L _ { 2 } \cdot \mathbf L _ { 3 } ) ^ { 2 } + ( \mathbf L _ { 1 } \cdot \mathbf L _ { 3 } ) ^ { 2 } ].$$

Comparing to the expression (1) in the main text, here we have 𝑢 = 2𝑣K + 6`´ " , 𝑔 = `´ GLYPH&lt;176&gt; , and 𝑤 = 2𝑣6 . A stripe ground state, in which only one of the 𝐋8 becomes non-zero, takes place when 𝑔 &gt; 0 and 𝑤 &gt; -12𝑔 . This state not only breaks spin-rotational symmetry, but also the 𝐶GLYPH&lt;176&gt;GLYPH&lt;144&gt; rotational symmetry, since one among the three 𝐐8 wave-vectors is chosen. Within a mean-field approach, as long as 𝑢 is large enough, both the magnetic (spin-rotational symmetry-breaking) and nematic ( 𝐶GLYPH&lt;176&gt;GLYPH&lt;144&gt; to 𝐶6GLYPH&lt;144&gt; symmetry-breaking) transitions happen at the same temperature 𝑇K as a second-order transition when 𝑎 ∝ 𝑇 - 𝑇 K &lt; 0 .

Fluctuations, however, completely change this picture. Formally, fluctuations can be included self-consistently by decoupling the first two quartic terms in Eq. (13) via three Hubbard-Stratonovich fields:

$$\psi & = u ( \mathbf L _ { 1 } ^ { 2 } + \mathbf L _ { 2 } ^ { 2 } + \mathbf L _ { 3 } ^ { 2 } ) \\ n _ { 1 } & = g ( \mathbf L _ { 1 } ^ { 2 } + \mathbf L _ { 2 } ^ { 2 } - 2 \mathbf L _ { 3 } ^ { 2 } ) \\ n _ { 2 } & = g \sqrt { 3 } ( \mathbf L _ { 1 } ^ { 2 } - \mathbf L _ { 2 } ^ { 2 } ).$$

Physically, 𝜓 is simply the amplitude of the magnetic fluctuations. As for 𝑛d and 𝑛6 , they transform as 𝑑 V ¨ GLYPH&lt;145&gt;GLYPH&lt;136&gt; ¨ quadrupolar order and as 𝑑VGLYPH&lt;136&gt; quadrupolar order. These are precisely the transformation properties of the two components of a generic twodimensional electronic nematic order parameter. In terms of the 𝐷GLYPH&lt;176&gt; point group, the vector 𝐧 = (𝑛d, 𝑛 6 ) transforms as the two-dimensional 𝐸6 irreducible representation. In

terms of the Hubbard-Stratonovich fields in Eq. (14), the free energy density (13) is given by:

$$F [ \mathbf L _ { i }, \psi, \phi _ { i } ] = ( a + \psi ) \sum _ { i } \mathbf L _ { i } ^ { 2 } - \frac { \psi ^ { 2 } } { 2 u } + \frac { n ^ { 2 } } { 2 g } - n _ { 1 } ( \mathbf L _ { 1 } ^ { 2 } + \mathbf L _ { 2 } ^ { 2 } - 2 \mathbf L _ { 3 } ^ { 2 } ) - n _ { 2 } \sqrt { 3 } ( \mathbf L _ { 1 } ^ { 2 } - \mathbf L _ { 2 } ^ { 2 } )$$

with 𝜙 6 = 𝜙d 6 +𝜙6. 6 The partition function is given by:

$$Z = \int \, \mathcal { D } \mathbf L _ { i } \, \mathcal { D } \psi \, \mathcal { D } n _ { i } \exp \left ( - \int _ { q } F [ \mathbf L _ { i }, \psi, \phi _ { i } ] \, / T _ { 0 } \right ).$$

Here, ∫ 	≡	 ∫ GLYPH&lt;209&gt; GLYPH&lt;211&gt; GLYPH&lt;212&gt; GLYPH&lt;209&gt; (6GLYPH&lt;213&gt;) GLYPH&lt;212&gt; and fluctuations of the 𝐋8 fields are taken into account by including the gradient term in Eq. (15), yielding:

$$F [ \mathbf L _ { i }, \psi, \phi _ { i } ] = \sum _ { i } \chi _ { i i } ^ { - 1 } ( \mathbf q ) \, \mathbf L _ { i, q } \, \mathbf L _ { i, - q } \, - \frac { \psi ^ { 2 } } { 2 u } + \frac { n ^ { 2 } } { 2 g }.$$

The diagonal inverse magnetic susceptibility matrix 𝜒 88 GLYPH&lt;145&gt;d is given by:

$$\chi _ { i i } ^ { - 1 } ( \mathbf q ) = \begin{pmatrix} a + \psi + q ^ { 2 } - n _ { 1 } - \sqrt { 3 } n _ { 2 } & 0 & 0 \\ 0 & a + \psi + q ^ { 2 } - n _ { 1 } + \sqrt { 3 } n _ { 2 } & 0 \\ 0 & a + \psi + q ^ { 2 } + 2 n _ { 1 } \end{pmatrix} \quad \ \ ( 1 8 )$$

In the paramagnetic state, magnetic fluctuations can be integrated out. It is convenient to change variables to 𝑟 ≡ 𝜓 + 𝑎 ;  the resulting partition function becomes:

$$Z = \, \Big | \, \mathcal { D } r \, \mathcal { D } n _ { i } \exp ( - f _ { \text{eff} } [ r, n _ { i } ] ),$$

where we absorbed the constant 𝑇K in the coupling constants 𝑢 and 𝑔 and derived the free energy:

$$f _ { \text{eff} } [ r, n _ { i } ] = \frac { N } { 2 } \sum _ { i } \int _ { q } \log ( \chi _ { i i } ( \mathfrak { q } ) ) - \frac { ( r - a ) ^ { 2 } } { 2 u } + \frac { n ^ { 2 } } { 2 g }.$$

Here, 𝑁 is the number of components of the 𝐋8 order parameters. Instead of working with 𝑛d and 𝑛6 , it is convenient to parametrize 𝑛d = 𝑛 cos2𝜃 and 𝑛6 = 𝑛 sin 2𝜃 . The magnetic susceptibility matrix becomes, in this parametrization:

$$& \chi _ { i i } ^ { - 1 } ( \mathfrak { q } ) \\ & = \begin{pmatrix} r + q ^ { 2 } - 2 n \cos \left ( 2 \theta - \frac { \pi } { 3 } \right ) & 0 & 0 \\ 0 & r + q ^ { 2 } - 2 n \cos \left ( 2 \theta - \frac { 5 \pi } { 3 } \right ) & 0 \\ 0 & 0 & r + q ^ { 2 } - 2 n \cos ( 2 \theta - \pi ) \end{pmatrix}.$$

Up to this point, all results are exact. To proceed, we must perform some approximation. It is instructive to expand the free energy in powers of 𝑛 only. We find:

$$\frac { f _ { e f f } } { N } = \frac { n ^ { 2 } } { 2 } \left [ \frac { 1 } { g } - \int _ { q } \frac { 3 } { ( r + q ^ { 2 } ) ^ { 2 } } \right ] + \left [ \int _ { q } \frac { 1 } { ( r + q ^ { 2 } ) ^ { 3 } } \right ] n ^ { 3 } \cos 6 \theta - \left [ \int _ { q } \frac { 9 } { ( r + q ^ { 2 } ) ^ { 4 } } \right ] n ^ { 4 } + ( \cdots ),$$

where, as customary, we renormalized (𝑢, 𝑔) → (𝑢, 𝑔)/𝑁 . This is precisely the free energy for the 3-state Potts model discussed in the main text. The key point is that, from Eq. (21), 𝑟 can be interpreted as the inverse squared magnetic correlation length 𝜉 GLYPH&lt;145&gt;6 in the absence of Potts-nematic order. Thus, 𝑟 must vanish at the bare magnetic transition. This implies that the quadratic coefficient of 𝑛 6 changes sign from positive (when 𝑟 → ∞ ) to negative (when 𝑟 → 0 ) before the magnetic transition is reached. As a result, Pottsnematic order appears already in the paramagnetic phase. Note that, because of the positive cubic coefficient, the free energy is minimized by one of the three values 𝜃 = 𝜋/6,𝜋/2,5𝜋/6 . Note also that in principle 𝑟 can also be a function of 𝑛 , which will alter the quartic coefficient in front of 𝑛 GLYPH&lt;128&gt; .

The existence of the cubic coefficient implies that, in three dimensions, the Pottsnematic transition is mean-field and first-order. This opens up the possibility of a simultaneous first-order Potts-nematic and magnetic transition, depending on how large the jump in 𝑛 is. This can be seen directly from the inverse susceptibility, Eq. (21): it is clear that the largest component of the 𝑞 = 0 magnetic susceptibility is given by 𝜒88 = 1/(𝑟 - 2𝑛) , where we used the fact that 𝜃 assumes one of the values 𝜃 =

𝜋/6,𝜋/2,5𝜋/6 . Now, 𝑛 jumps by ∆𝑛 for a finite value of 𝑟 = 𝑟 GLYPH&lt;238&gt; . If however ∆𝑛 &gt; 𝑟 GLYPH&lt;238&gt; /2 , this implies a diverging renormalized magnetic susceptibility, and thus a simultaneous firstorder magnetic transition.

Similar results hold also if the three ordering vectors 𝐐8 are incommensurate and aligned along the Γ - 𝐾 and Γ - 𝐾′ directions. For instance, a zigzag phase is obtained when the ordering vectors are given by

$$\mathbf Q _ { 1 } = \frac { \pi } { 2 } ( 1, \sqrt { 3 } ), \ \mathbf Q _ { 2 } = \frac { \pi } { 2 } ( - 2, 0 ), \ \mathbf Q _ { 3 } = \frac { \pi } { 2 } ( 1, - \sqrt { 3 } ).$$

These ordering vectors are located on the Γ - 𝐾 and Γ - 𝐾′ high symmetry lines, are related by threefold rotation, and have the property that 𝐐8 ≠ -𝐐8 . The latter implies that, in contrast to Eq. (12), the corresponding order parameters 𝐋8 are complex. As a result, the nematic order parameter is now given by:

$$& n _ { 1 } \in ( | L _ { 1 } | ^ { 2 } + | L _ { 2 } | ^ { 2 } - 2 | L _ { 3 } | ^ { 2 } ) \\ & n _ { 2 } \in \sqrt { 3 } \left ( | L _ { 1 } | ^ { 2 } - | L _ { 2 } | ^ { 2 } \right )$$

## V. Domain Repopulation

We model the repopulation of the ℤ" -nematic domains under uniaxial strain, demonstrated in the main text Fig. 5, with a free energy, 𝐹 , that contains two terms,

$$F = [ ( x _ { 1 } - x _ { 2 } ) ^ { 2 } + ( x _ { 2 } - x _ { 3 } ) ^ { 2 } + ( x _ { 3 } - x _ { 1 } ) ^ { 2 } ] - b \sum _ { n = 1 } ^ { 3 } x _ { n } \epsilon \cos 2 ( \theta - \theta _ { n } ).$$

The first term favors equal domain populations whereas the second term breaks the degeneracy of the three domains. Here 𝑥{ is the areal fraction of the 𝑛 §GLYPH&lt;242&gt; domain, 𝑏 is the strain-nematic coupling parameter, 𝜖 is the amplitude of the uniaxial strain, 𝜃 - 𝜃 { is angle of the strain with respect to the director of domain 𝑛 , and 𝜃{ = 2𝜋𝑛/3 . Using the constraint Σ𝑥{ = 1 and linearizing with respect to 𝛿{ ≡ 𝑥8 -1/3 , we find that 𝐹 is minimized for 𝛿{ = (3𝑏𝜖/4)cos 2(𝜃 - 𝜃 { ) . The model predicts that for tensile strain along a bond the domain aligned with the strain grows according to 𝛿K = 3𝑏𝜖/4 , while

the other two shrink, 𝛿d = 𝛿6 = -3𝑏𝜖/8 , and the signs of the 𝛿{ are reversed for tensile strain perpendicular to the bond direction.

The domain repopulation in the ℤ" nematic enables the global symmetry axes to continuously follow the direction of in-plane 𝑩 𝑱 , , or uniaxial strain fields, despite the locking of the Neel vector perpendicular to the plane. As an indicator of the global symmetry, we consider the resistivity anisotropy tensor, 𝝆 ⃗ ⃖ { ≡ 𝑅GLYPH&lt;247&gt; ł Δ𝜌𝜎GLYPH&lt;144&gt; where 𝑅GLYPH&lt;247&gt; ł is the rotation operator and Δ𝜌 and 𝜎GLYPH&lt;144&gt; are the resistivity anisotropy and the Pauli 𝑧 matrix, respectively. We assume that the global symmetry axes align with the spatial average tensor, 〈𝝆 ⃗〉 ⃖ = Δ𝜌∑𝛿{,𝝆 ⃗ { ⃖ . At equilibrium, all three domains are equally populated ( 𝑥8 = 1/3 ) and 〈𝝆 ⃗〉 ⃖ is a null tensor, corresponding to a globally isotropic system. However, an external field applied at an angle 𝜃 unbalances the domain population, in which case,

$$F ( \stackrel { \_ } { \rho } ) \, \varpropto \, \Delta \rho \begin{bmatrix} \cos 2 \theta & \sin 2 \theta \\ - \sin 2 \theta & \cos 2 \theta \end{bmatrix},$$

which shows that rotational symmetry is broken and the global symmetry axis continuously follows the direction of the perturbation. This is in stark contrast with the ℤ6 -Ising-nematic system where the global anisotropy axes are the locked to the crystal axes, regardless of the direction of an external perturbation.

## References

- 1. Doyle, S. et al. Tunable Giant Exchange Bias in an Intercalated Transition Metal Dichalcogenide. Preprint at https://arxiv.org/abs/1904.05872 (2019).
- 2. Fan, S. et al. Electronic chirality in the metallic ferromagnet Fe1/3TaS2. Phys. Rev. B 96 , 1Ð6 (2017).
- 3.        Guenther, R.D. Modern Optics. Wiley, (1990).
Arielle Little Changmin Lee Caolan John Spencer Doyle Eran Maniv Nityan L. Nair Wenqin Chen Dylan Rees Jörn W. F. Venderbos Rafael Fernandes James G. Analytis Joseph Orenstein
Observation of three-state nematicity in the triangular lattice antiferromagnet Fe$_{1/3}$ NbS$_2$
Nematic order is the breaking of rotational symmetry in the presence of
translational invariance. While originally defined in the context of liquid
crystals, the concept of nematic order has arisen in crystalline matter with
discrete rotational symmetry, most prominently in the tetragonal Fe-based
superconductors where the parent state is four-fold symmetric. In this case the
nematic director takes on only two directions, and the order parameter in such
"Ising-nematic" systems is a simple scalar. Here, using a novel
spatially-resolved optical polarimetry technique, we show that a qualitatively
distinct nematic state arises in the triangular lattice antiferromagnet
Fe$_{1/3}$NbS$_2$. The crucial difference is that the nematic order on the
triangular lattice is a Z$_3$, or three-state Potts-nematic order parameter. As
a consequence, the anisotropy axes of response functions such as the
resistivity tensor can be continuously re-oriented by external perturbations.
This discovery provides insight into realizing devices that exploit analogies
with nematic liquid crystals.