Z3-vestigial nematic order due to superconducting fluctuations in the doped topological insulators NbxBi2Se3 and CuxBi2Se3

A state of matter with a multi-component order parameter can give rise to vestigial order. In the vestigial phase, the primary order is only partially melted, leaving a remaining symmetry breaking behind, an effect driven by strong classical or quantum fluctuations. Vestigial states due to primary spin and charge-density-wave order have been discussed in iron-based and cuprate materials. Here we present the observation of a partially melted superconductivity in which pairing fluctuations condense at a separate phase transition and form a nematic state with broken Z3, i.e., three-state Potts-model symmetry. Thermal expansion, specific heat and magnetization measurements of the doped topological insulators NbxBi2Se3 and CuxBi2Se3 reveal that this symmetry breaking occurs at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{T}}_{\mathrm{nem}} \simeq 3.8\,K$$\end{document}Tnem≃3.8K above \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{{\mathrm{c}}} \simeq 3.25\,K$$\end{document}Tc≃3.25K, along with an onset of superconducting fluctuations. Thus, before Cooper pairs establish long-range coherence at Tc, they fluctuate in a way that breaks the rotational invariance at Tnem and induces a crystalline distortion.

N ematic electronic phases with vestigial order 1,2 are known from iron-based superconductors and cuprates, where it has been suggested that the nematic phase and the nearby spin-and charge-density wave states are not independent competing but intertwined electronic phases [1][2][3][4][5][6][7][8][9][10] . The density-wave states are the primary electronic phases and characterized by a multi-component order parameter. The nematic phase is a fluctuation-driven phase and characterized by a composite order parameter. Then the spin-or charge density-wave order melts partially, but leaves an Ising, i.e., Z 2 -nematic state as a vestige. Vestigial order whose primary order is superconductivity has not been observed. Such partially molten superconductivity requires a material with unconventional, multi-component order parameter, and strong pairing fluctuations.
When the topological insulator Bi 2 Se 3 is doped with electrons, e.g., by intercalation of Cu, Sr, Nb, or other metal ions in its layered structure, a superconducting state is formed 11,12 . The presence of a strong spin-orbit coupling, which also manifests itself in a topological surface state of the parent insulator 13 , led to the proposal of unconventional pairing with an odd-parity symmetry and topological superconductivity 14 . The low carrier concentration, the layered structure, and the low ratio ξ/λ F of the superconducting coherence length and the Fermi wavelength 11,12 , strongly enhance fluctuation effects. In addition, numerous experiments have shown that the superconducting state is accompanied by a spontaneous breaking of rotational symmetry with a pronounced twofold anisotropy within the Bi 2 Se 3 basal plane [15][16][17][18][19][20][21][22] ; see Ref. 23 for a recent review. The twofold symmetry can be observed in field-angle resolved experiments where a magnetic field is rotated in the plane with respect to the crystalline axes and the corresponding physical quantity (e.g., spin susceptibility, specific heat, magneto-resistance, upper critical field, magnetization, magnetic torque) is represented as a function of angle [15][16][17][18][19][20][21][22][23] . This behavior directly reflects the anisotropy of the superconducting state. Thus, doped Bi 2 Se 3 is an unconventional nematic superconductor with a pairing wave function in either the two-component E u or E g point group representation, the only pairing states that spontaneously break the trifold crystal symmetry within the basal plane. The temperature dependence of the penetration depth of Ref. 24 supports point nodes, consistent with E u odd-parity pairing.
In this article we report on high-resolution thermal expansion experiments on a superconducting mono-crystalline Nb-doped Bi 2 Se 3 sample in combination with electrical transport, DC magnetization, and specific heat data demonstrating a Z 3 -vestigial nematic phase with enhanced superconducting fluctuations. We have measured the linear thermal expansion in three different crystalline directions in the Bi 2 Se 3 basal plane and observed a strong anisotropic expansion occurring at a temperature of~0.5 K above the superconducting transition. Our high-resolution magnetization, electrical resistivity, and specific heat data-after zooming near T c -show that an anomaly with increasing superconducting fluctuations occurs at the nematic transition. As we will explain below, these observations are perfectly consistent with a vestigial nematic phase of symmetry-breaking pairing fluctuations, recently predicted in Ref. 25 . This observation of a genuine symmetry breaking of pairing fluctuations above T c is qualitatively distinct from the gradual onset of order-parameter fluctuations in the disordered phase 26 or the crossover to Bose-Einstein condensation of pairs 27 . It corresponds to a sharply defined state of matter that might, e.g., undergo a separate quantum phase transition when a magnetic field is applied in the plane at low temperatures. Qualitatively similar results are obtained both on another Nb-doped Bi 2 Se 3 single crystal and on a Cu-doped Bi 2 Se 3 single crystal, thus, demonstrating the reproducibility and universality of the observed features.

Results
Magnetoresistance and nematicity. The advantage of the Nbdoped Bi 2 Se 3 system is that single crystals with a high superconducting volume fraction and a complete zero resistance can be found, as the results presented here show. All our bulk thermodynamic data (thermal expansion, magnetization, and specific heat) show relatively large anomalies at the superconducting transition. Figure 1a shows magnetoresistance data recorded at 0.35 K with the magnetic field applied strictly parallel to the Bi 2 Se 3 basal plane for different directions in the plane with respect to the yields Γ % 3:32Ω and H c2 (0) = 1.42 T as a measure of the anisotropy in the basal plane. It should be noted that the normal state resistance well above the upper critical field has no variation for the different orientations of the magnetic field in the plane, indicating an isotropic normal state within the trigonal basal plane. We have previously found that the orientation of this nematic superconducting order parameter for this sample always appears to be pinned along the same of the three equivalent crystal directions in the Bi 2 Se 3 basal plane, even if the sample is warmed to room temperature between different experiments 19 . The origin of this preference for a particular direction is unknown, but likely associated with microscopic details of the sample morphology, such as internal strain or microcracks (see Supplementary Discussion and Supplementary Fig. 6 for more details).
Thermal expansion and vestigial order. Through thermal expansion experiments, we have a highly sensitive bulk thermodynamic probe that is not only sensitive to the anharmonicity of phononic contributions but also to electronic degrees of freedom including nematic and superconducting order. We focus here on ΔLðT; HÞj μ =L 0 measured along different directions μ in the Bi 2 Se 3 basal plane, which directly represents the change in length ΔL of the sample as a function of temperature or magnetic field, normalized to the length L 0 at ambient temperature. This quantity is directly related to the linear thermal expansion coefficient α μ (T,H) = 1/L 0 dL μ (T)/dT. Figure 2a shows the linear thermal expansion ΔL| μ /L 0 measured along the three directions within the Bi 2 Se 3 basal plane of 90°, 155°, and 215°. All data fall perfectly on each other in the normal state, but begin to deviate gradually from each other below 3.8 K. In the following, we will refer to this characteristic temperature as T nem , because here the onset of a twofold crystalline distortion and thus nematicity occurs. The twofold crystalline distortion with a relative length change ΔL/L 0 = 2 × 10 −7 amounts to a distortion of less than 0.1 femtometers within the unit cell. Still, these minute changes smaller than the size of the proton, are clearly resolvable in our measurements. The distortion is correlated with the upper critical field, with a small negative length change along 90°where the H c2 minimum occurs (Fig. 1b), and large positive anomalies at 155°and 215°, both near the mean H c2 value. At 3.25 K, much smaller anomalies are visible that can be identified as the superconducting transition, as the comparison with the specific heat (shown in the same graph) reveals. In the specific heat, the anomaly at T nem is obscured by the phonon background, but becomes visible in the temperature derivative of C/T (Fig. 2b), where a small step-shaped anomaly occurs. The anomalies in ΔL/L 0 at T nem show up as a somewhat broadened step. As a first-order derivative of the free energy, a step-like transition in ΔL/L 0 is the characteristic signature of a first-order transition, while a second-order transition would appear as a kink.
In distinction, the superconducting transition at T c remains the standard second-order transition, as evidenced by the jump in the specific heat. Figure 2c shows the Meissner signal in the zero-field cooled and field-cooled DC magnetization, which agrees with T c ≈ 3.25 K obtained from the specific heat. We also show the same data, but with a magnification of 10 5 , to illustrate that an enhanced diamagnetic response, signaling superconducting fluctuations, already sets in at T nem . The electrical resistivity in Fig. 2d shows a similar trend with a drop in resistivity consistent with paraconductivity, i.e., superconducting fluctuations, well above the main transition. Our X-ray diffraction results (see section "Methods") show that the doped Bi 2 Se 3 phase of the R-3m space group is the majority phase responsible for T c at 3.25 K. As minority phases Bi 2 Se 3 of space group P-3m1 and NbBiSe 3 were found. The latter occupy far too little volume to explain such large anomalies in thermal expansion. From this, we conclude that the observed crystalline distortion below T nem = 3.8 K is caused by a transition separate from the main superconducting transition, but linked to the -lin. backgr. occurrence of superconducting fluctuations, which cause a weak Meissner effect and decrease in resistivity below this temperature.
Two stage transition. Our data thus show that nematic superconductivity occurs in the form of a two-stage transition, see also Fig. 3 for an illustration. The distortion forms near the higher onset temperature T nem , where superconducting fluctuations in the magnetization and the first-order derivative of the specific heat are visible. The superconducting transition occurs at a lower temperature T c , which corresponds to the formation of a global phase-coherent superconducting state. The signs and magnitudes of the length changes for the three measuring directions are consistent with the indicated distortions, see "Methods" section.
It could be argued that thermal expansion reveals a separate structural transition at T nem that has nothing to do with superconductivity 28 . Such a sequence of independent or competing transitions would be allowed within the Landau theory of phase transitions. Obviously, the anomaly in the diamagnetic response at T nem is already strong evidence that this is not the case. Furthermore, in Fig. 4 we show magnetostriction data for the 155°direction, along which we observed the greatest change in length in Fig. 2a. Here, ΔL/L 0 was measured at a fixed temperature of 350 mK as a function of the magnetic field. A broad step-like transition occurs with a total length change of ΔL/L 0 ≈ 0.22 × 10 −7 with onset at~1.0 T. Figure 1 shows that the resistively determined H c2 for this direction occurs at 0.7 T. Magnetostriction shows that the crystalline distortion is constant up to this field where a kink occurs, then it is gradually removed up to 1.0 T. The overall anomaly represents a broad step in ΔL/L 0 , which is the expected characteristics of a first-order transition. A small hysteresis can be seen in the data measured upon sweeping the field up and down. Experimental artifacts as source for the hysteresis have been carefully excluded, especially since the field scanning speed was kept very slow at 0.02 T/min, which typically does not cause any hysteresis-like effects in reversible samples. Given the layered structure of the sample with the field aligned in parallel, this hysteresis is most likely a consequence of flux pinning effects, which are typically strong at such a low temperature. Therefore, the hysteresis should not necessarily be regarded as evidence of a first-order H c2 nature. The field-induced length change at low temperature shown in Fig. 4 corresponds to the temperature-induced length change at zero field shown in Fig. 2a. This observation provides further evidence that the nematic distortion is closely linked to the superconducting state, with a separate nematic transition H nem occurring above the main superconducting H c2 transition, as shown in Fig. 4.
Further data on thermal expansion measured on a second sample of the same batch are shown in Supplementary Figs. 1 and 2 and are discussed in Supplementary Note 1. They show a similar behavior, although the crystalline distortion occurring below T nem is weaker due to a multi-domain structure. For this sample we have also measured thermal expansion in fixed magnetic fields, and it can be seen that T nem is suppressed by the magnetic field together with T c , further confirming that the two transitions are closely related. We also show similar data of a Cu 0.2 Bi 2 Se 3 single crystalline sample (Sample 3) in Supplementary  Fig. 3-5 and the data are discussed in Supplementary Note 2.

Discussion
The two thermodynamic properties, the specific heat and the linear thermal expansion coefficient α μ (T) = 1/L 0 dL μ (T)/dT are closely related in the vicinity of a phase transition through the Clausius Clapeyron (2) and Ehrenfest (3) relation for first and second order phase transitions, respectively. The proportionality is the uniaxial pressure dependence dT c /dp μ of the transition (T c is the critical temperature at which the phase transition occurs, p μ is the uniaxial pressure applied along a certain crystalline direction μ, V mol is the molar volume, ΔS is the jump in entropy at a SC + Nematic Nematic Metallic ig. 3 Lattice distortion and real space image of the superconducting order parameter. Purple dashed lines below T nem indicate the high-temperature atomic positions. The lattice distortion is strongly exaggerated. Arrows in the nematic phase indicate the directions of the thermal expansion measurements. In the latter, strong superconducting fluctuations break the discrete lattice symmetry without broken U(1) symmetry and superconducting coherence. A globally coherent superconducting state only sets in at the superconducting transition temperature T c . The data (Sample 1) was measured along the 155°direction at a fixed temperature of 350 mK. The crystalline distortion is removed at H nem = 1 T, well above the resistively determined H c2 , which coincides with the kink at 0.7 T for this field orientation (marked as H c2 ). The magnetic field was applied parallel to the measured sample direction. A weak linear normal state background was subtracted for clarity. The inset shows the magnetic phase diagram.
first order transition and ΔC p is the jump in the specific heat at constant pressure at a second order phase transition). dT c dp μ ¼ A large anomaly in thermal expansion and a small anomaly in specific heat means in both cases that T nem is strongly dependent on uniaxial pressure and the electronic nematic order is strongly coupled to the crystalline lattice. The strong crystalline distortion observed here using linear thermal expansion therefore means that the nematic transition is strongly dependent on pressure or strain. Such a behavior can also be observed, for example, in iron based superconductors, where a nematic transition occurs in the vicinity to a spin density wave transition and causes large anomalies in thermal expansion 29 .
Our findings can be explained in terms of vestigial order due to superconducting fluctuations. In fact, recently it has been suggested that such vestigial order should emerge from the superconducting phase in doped Bi 2 Se 3 25 . On the one hand, the superconducting order parameter of either the E g or the E u representation has two components 14,30-32 that are characterized by the overall amplitude Δ 0 , the global U(1) phase φ and three distinct values of the angle θ ¼ π 6 ; π 2 ; 5π 6 È É that select a specific crystalline axis. Superconducting fluctuations will then induce a phase transition to a vestigial nematic state at a temperature T nem above T c . While superconductivity is signaled by a finite expectation value of Δ x and/or Δ y , the nematic phase is characterized by a finite expectation value of the composite order parameter Upon increasing the temperature, superconducting fluctuations continue to break the rotational symmetry, even after restoration of global U(1) symmetry at T c . The composite order parameter Q μν is made up of combinations of the superconducting order parameter, similar to charge-4e superconductivity proposed within the context of pair-density wave order in cuprate superconductors 33,34 or proton-electron superconducting condensate in liquid hydrogen 35 . As a traceless second-rank tensor, hQ μν i ¼ Q 0 n μ n ν À 1 2 δ μν behaves, however, like a nematic order parameter with director n ¼ cos θ; sin θ ð Þ 36 and strongly couples to the strain tensor ε μν via κtrðQεÞ with nemato-elastic coupling constant κ. A nonzero Q 0 then induces a lattice distortion ε μν / κQ μν , see Fig. 3. Thus, the lattice can be utilized to detect this unconventional electronic order. The point group analysis further yields a first-order transition at T nem into a state with Q 0 ≠ 0, since it is in the three-state Potts model, i.e., the Z 3 universality class. The superconducting transition continues to be of second order, all in agreement with our experimental findings. In Fig. 5a we show the nematic (Q 0 ) and superconducting (Δ 0 ) order parameters and in Fig. 5b the diamagnetic susceptibility obtained within the theory of Ref. 25 . The susceptibility is compared with the data of Fig. 2c, where the logarithmic axis is used to illustrate the rapid growth of diamagnetic fluctuations below T nem . While the in-plane anisotropy of the susceptibility is only finite below the nematic transition, the crystal symmetry allows χ zz to be distinct already above T nem . The magnitude of the out-of-plane anisotropy is determined by the ratio of the electron velocities in the corresponding directions. The anisotropy of H c2 ϕ ð Þ shown in Fig. 1b (orange line) was also obtained within the same theory and is compared with the behavior without nematic order (Q 0 = 0) where H c2 should have sixfold symmetry 24,25 . Without nematic phase above T c , the superconducting order parameter directly at H c2 ϕ ð Þ is infinitesimal and no twofold rotational symmetry breaking should be visible, in clear contrast to experimental observations. In Ref. 31 , the twofold symmetric behavior of H c2 ϕ ð Þ only occurred after an additional symmetry-breaking strain was added. Vestigial nematic order offers a natural explanation for this strain field.
After this work was completed, we learned about Ref. 37 , where a twofold symmetry breaking above T c is reported. Our results agree with those of Ref. 37 and make evident that the high temperature phase is separated by an actual first order transition where superconducting fluctuations are enhanced. Furthermore, Refs. 38,39 reported on the control of nematic superconductivity by uniaxial strain, which is consistent with our observation of a coupling of the nematic order parameter to the crystal lattice.
To summarize, our data demonstrate that a separate nematic transition occurs in the doped topological insulator Nb 0.25 Bi 2 Se 3 at T nem = 3.8 K, i.e., about 0.5 K above T c , with a distinct crystalline distortion occurring in the Bi 2 Se 3 basal plane. T nem coincides with the onset temperature of superconducting fluctuations. The direction of the crystalline distortion is correlated with the direction of the twofold symmetry of the superconducting order parameter and is removed together with the superconductivity at or near the upper critical field H c2 . The two transitions are thus interconnected. Our observations are perfectly consistent with vestigial nematic order and a sequential restoration of U(1) and rotational symmetry. The new nematic phase is a state of matter in which Cooper pairs have lost their off-diagonal long-range order, yet fluctuate in a way that breaks the rotational symmetry of the crystalline lattice.

Methods
Sample characterization. The monocrystalline Nb 0.25 Bi 2 Se 3 sample used in this study was selected because of its particularly large T c anomalies in the specific heat, which indicates a high superconducting volume fraction, and because of its particularly large nematic in-plane H c2 anisotropy 19 . Our previous work also demonstrated that it forms one large nematic domain comprising~90% of the superconducting volume fraction in which the nematic order parameter is pinned  Supplementary Fig. 3-5. These samples show a qualitatively similar behavior as Sample1.
A laboratory X-ray Laue equipped with CCD camera (Photonic Science) was used to characterize the crystal quality and to determine the crystalline directions in the basal Bi 2 Se 3 plane. The Laue images (Fig. 6) were taken on the shiny surface of the sample after cleaving off a thin layer with shooting X-ray beam along the caxis, revealing the crystal orientation and proving the hexagonal structure. We took data on different spots on the sample surface and found that the change in the crystalline direction was less than 0.2°over a distance of 0.8 mm, proving the sufficiently good single crystalline quality of the sample.
Powder X-Ray Diffraction pattern was collected at room temperature in the Bragg-Brentano geometry using a Bruker AXS D8 Advance diffractometer equipped with a Ni-filtered Cu Kα radiation and a 1D LynxEye PSD detector (Fig. 7). The reflections were indexed to Bi 2 Se 3 (SG R-3m; No 166) as a main phase. In addition, some minority phases were found as following: Bi 2 Se 2 and NbBiSe 3 with corresponding space groups: P-3m1 (No 164) and P2 1 2 1 2 1 (No 19), respectively. This demonstrates that, while most of the Nb must be intercalated between the Bi 2 Se 3 layers, some Nb is incorporated into the layers on the Bi sites. This agrees with literature data 40 . The Rietveld refinement 41 of the diffraction patterns was done by the package FULLPROF SUITE 42 (version July-2019) using a previously determined instrument resolution function (based on the small line width polycrystalline sample Na 2 Ca 3 Al 2 F 14 measurements 43 ). Refined parameters were: scale factor, zero displacement, lattice parameters, atomic positions, isotropic Debye-Waller factors, and peak shape parameters as a Thompson-Cox-Hastings pseudo-Voigt function. Determined lattice parameters of the rhombohedral Bi 2 Se 3 are equal a = b = 4.1854(4) Å, and c = 28.4633(7) Å. Because of the habit of the powdered crystal, a preferred orientation as a March-Dollase multi-axial phenomenological model was implemented in the analysis.
Experimental techniques. The high-resolution linear thermal expansion was measured with a capacitive technique using a dilatometer, in which the sample is pressed by a fine screw mechanism against a cantilever forming one of the two plates of a capacitor. A change in the sample length leads to a change in the separation of the capacitor plates, which can be determined with a General Radio 1615A capacitance bridge in combination with a Stanford Research SR830 digital lock-in amplifier. Before the experiments we have carefully checked that the empty dilatometer does not show any anomalies in the temperature range of interest ( Supplementary Fig. 7). We have measured the thermal expansion as a function of temperature in three different directions within the basal Bi 2 Se 3 plane (90°, 155°, and 215°). 0°is the a direction normal to the mirror plane and corresponds to the magnetic-field direction in the plane providing the maximum upper critical field, while 90°corresponds to the a* direction parallel to the mirror plane. The crystalline directions have been obtained by Laue X-ray diffraction as explained above. The other directions were chosen as representative of other characteristic directions in the plane, but largely dictated by the crystal shape, which allowed a stable mounting of the sample only for certain directions. All directions were within 5°f rom the three different a* directions. For technical reasons, all data was taken at a slow rate of 0.02 K/min upon increasing temperature. Figure 8 shows a photograph of the mounted crystal for the three different orientations in the dilatometer. The dilatometer 44 was well characterized using separate calibration measurements: e.g., a measurement with a Cu sample of 1 mm length of the same material as the dilatometer body shows only a very small temperature dependence without significant anomalies in the temperature range of interest. The absolute value of the change in length was calibrated using a 2 mm long undoped silicon sample, which gave a linear thermal expansion coefficient that agrees well with literature 45 (Supplementary Fig. 8).
The specific heat was measured with a home-made calorimeter, which can be used either in AC modulated temperature mode or in long relaxation mode. The long relaxation mode provides high accuracy in the absolute value of 1% precision, while the AC mode provides high relative resolution with a high density of data points of 1000 points per K. The data presented in this letter has been acquired using the AC technique, but the absolute value has been calibrated using the relaxation technique.
DC magnetization was measured using a commercial Quantum Design Vibrating Sample SQUID magnetometer and the electrical resistance was measured using a standard 4-probe technique with a Keithley 6221 AC current source combined with a SR830 digital lock-in amplifier. For the latter, a low temperature piezo rotator was used to precisely align the magnetic field along the different crystalline directions and to study the H c2 anisotropy that reflects the nematic superconductivity. The rotator allowed milli-degree precision for relative changes in orientation. However, a systematic error of less than 5°can occur with respect to the measured crystalline axes.
Theory: q-state Potts model. A q-state Potts model describes a spin-like variable s = 1, 2, …, q that can attain q different values. In our case q = 3 and the three values label the three axes of the crystal along which one assumes a local displacement. The energy of two neighboring variables s and s' that are the same is then lower than for distinct variables, i.e., the bond energy goes like À 1 2 J 0 ðqδ s;s 0 À 1Þ. In general, the model has to be distinguished from a q-states clock model where an angle φ 2 0; 2π ½ can take q distinct values φ s ¼ 2πs q and two sites interact via an energy proportional to ÀJ 0 cosðφ s À φ s0 Þ. For q = 3 both models are equivalent though. Using the Potts language one can define the fraction n s of the lattice in the s-th state, which implies n 1 +n 2 +n 3 = 1. At high temperatures one expects m s = n s À 1 3 to have zero expectation value for all s. Below the transition temperature, one of the m s becomes positive and the two others negative. Because of the condition on the sum of the n s the three m s are not independent: m 1 +m 2 +m 3 = 0. Furthermore, the free energy of the system must be symmetric under m s ↔ m s′ . These two conditions lead to the unique free energy expansion up to quartic order Since there are really only two independent variables one can introduce the parametrization m 1 ¼ 2 ffiffi 3 p Q 1 and m 2;3 ¼ 1 2 À 1 ffiffi 3 p Q 1 ± Q 3 . The expansion in terms of the Q 1,2 leads up to constants to This is precisely the free energy expansion obtained in Ref. 25 in terms of the quadrupolar order parameter of Eq. (2) with Q 1 ¼ Δ x j j 2 − Δ y 2 and This demonstrates that the problem at hand is indeed in the universality class of the q = 3 Potts or clock models.
The sketch of the distorted unit cell in Fig. 3 is deduced from the linear coupling term κtr(Qε) between the strain tensor and the composite order parameter. Thus, the nematic order parameter acts in the same way as an applied external stress field in the E g symmetry channel and hence distorts the unit cell. Moreover, the linear coupling term does not entail a change of the unit cell volume, which is assumed to be unaltered in the following. Figure 9 shows the relative length changes in the three directions 90°, 155°, and 215°caused by an E g unit cell deformation as a function of the lattice parameter a/a 0 . For the strongly exaggerated value a/a 0 = 1.1 Fig. 6 Laue X-ray diffraction image of the Nb x Bi 2 Se 3 single crystal (Sample 1). This data allowed us to define the crystalline orientation in the Bi 2 Se 3 basal plane. we also show the corresponding distorted hexagon in the inset, which is quantitatively similar to the distorted unit cell in Fig. 3. While the computed relative length changes qualitatively capture the measured thermal expansion behavior, the calculated magnitude of the 90 o direction is slightly larger than that observed in experiment when compared with the other two directions. This could be due to the higher-order coupling that gives rise to a change of the unit cell volume at T nem .