Rotational symmetry breaking in the topological superconductor SrxBi2Se3 probed by upper-critical field experiments.

Recently it was demonstrated that Sr intercalation provides a new route to induce superconductivity in the topological insulator Bi2Se3. Topological superconductors are predicted to be unconventional with an odd-parity pairing symmetry. An adequate probe to test for unconventional superconductivity is the upper critical field, Bc2. For a standard BCS layered superconductor Bc2 shows an anisotropy when the magnetic field is applied parallel and perpendicular to the layers, but is isotropic when the field is rotated in the plane of the layers. Here we report measurements of the upper critical field of superconducting SrxBi2Se3 crystals (Tc = 3.0 K). Surprisingly, field-angle dependent magnetotransport measurements reveal a large anisotropy of Bc2 when the magnet field is rotated in the basal plane. The large two-fold anisotropy, while six-fold is anticipated, cannot be explained with the Ginzburg-Landau anisotropic effective mass model or flux flow induced by the Lorentz force. The rotational symmetry breaking of Bc2 indicates unconventional superconductivity with odd-parity spin-triplet Cooper pairs (Δ4-pairing) recently proposed for rhombohedral topological superconductors, or might have a structural nature, such as self-organized stripe ordering of Sr atoms.

Currently, topological insulators (TIs) are at the focus of condensed matter research, because they offer unprecedented possibilities to study novel quantum states123. 3D TIs are bulk insulators with a non-trivial topology of the electron bands that gives rise to surface states at the edge of the material. The gapless surface states have a Dirac-type energy dispersion with the spin locked to the momentum and are protected by symmetry. This makes TIs promising materials for applications in fields like spintronics and magnetoelectrics12. The concept of a TI can also be applied to superconductors, where the superconducting gap corresponds to the gap of the band insulator45. Topological superconductors are predicted to be unconventional with an odd-parity pairing symmetry67. Much research efforts are devoted to 1D and 2D superconductors, where Majorana zero modes exist as protected states at the edge of the superconductor89. Majorana zero modes with their non-Abelian statistics offer a unique platform for future topological quantum computation devices10. Prominent candidates for 3D topological superconductivity are the Cu intercalated TI Bi2Se31112, the doped topological crystalline insulator Sn1−InTe13 and selected topological half-Heusler compounds141516.

Among the 3D topological superconductors, CuBi2Se3, which has a superconducting transition temperature T = 3 K for x = 0.31112, is the most intensively studied material. ARPES (Angle Resolved PhotoEmission Spectroscopy) experiments conducted to study the bulk and surface states reveal that the topological character is preserved when Bi2Se3 is intercalated with Cu17. By evaluating the topological invariants of the Fermi surface, CuBi2Se3 is expected to be a time-reversal invariant fully-gapped odd-parity topological superconductor67. This was put on a firmer footing by a two-orbital pairing potential model where odd-parity superconductivity is favoured by strong spin-orbit coupling18. Several experiments have been interpreted in line with topological superconductivity. The specific heat shows a full superconducting gap12. The upper critical field exceeds the Pauli limit and has a temperature variation that points to spin-triplet superconductivity19. Much excitement was generated by the observation of a zero-bias conductance peak in point contact spectroscopy, that was attributed to a Majorana surface state20. However, STS (Scanning Tunneling Spectroscopy) showed that the density of states at the Fermi level is fully gapped without any in-gap states21. On the other hand, the superconducting state shows a large inhomogeneity21 and the superconducting volume fraction depends on quenching conditions22. Consequently, the issue of topological superconductivity in CuBi2Se3 has not been settled and further experiments are required, as well as new materials.

Very recently it has been demonstrated that Sr intercalation provides a new route to induce superconductivity in Bi2Se323. Resistivity and magnetization measurements on SrBi2Se3 single crystals with x = 0.06 show T = 2.5 K. The superconducting volume fraction amounts to 90% which confirms bulk superconductivity. By optimizing the Sr content a maximum T of 2.9 K was found for x = 0.1024. The topological character of Bi2Se3 is preserved upon Sr intercalation. ARPES showed a topological surface state well separated from the bulk conduction band2526. Based on the first measurements of the electronic parameters in the normal and superconducting states, and the close analogy to CuBi2Se3, it has been advocated that SrBi2Se3 is a new laboratory tool to investigate topological superconductivity2324.

Here we report a study of unusual basal-plane anisotropy effects in the upper critical field, B, of SrBi2Se3. Bi2Se3 crystallizes in a rhombohedral structure with space group . It is a layered material and Sr is intercalated in the Van der Waals gaps between the quintuple Bi2Se3 layers23. For a standard BCS (Bardeen, Cooper, Schrieffer) layered superconductor the anisotropy of B is expressed by the parameter , where and are measured with the B-field parallel and perpendicular to the layers, respectively27. Whereas is normally isotropic, SrBi2Se3 presents a unique exception. Field-angle-dependent magnetotransport experiments demonstrate a large two-fold basal-plane anisotropy of B, with T and T for x = 0.15 at T/T = 0.1 (T = 3.0 K), where a and a* are orthogonal directions in the basal plane. This large effect cannot be explained with the anisotropic effective mass model2728 or the variation of B caused by flux flow29. The rotational symmetry breaking of B indicates unconventional superconductivity3031, or might have a structural nature, such as preferential ordering of Sr atoms.

Results

The resistivity, ρ(T), of our SrBi2Se3 crystals with x = 0.10 and x = 0.15 shows a metallic temperature variation with superconducting transition temperatures T of 2.8 K and 3.0 K, respectively, see Fig. S4 in the Supplementary Information 32. The superconducting volume fractions of the crystals measured by ac-susceptibility amount to 40% and 80%, respectively32. In Fig. 1 we show the angular variation of the resistance, R(θ), measured in a fixed field B = 0.4 T directed in the basal plane (aa*-plane), in the temperature range 2–3 K around T (T = 2.8 K at B = 0 T), for x = 0.10. Rather than attaining a constant value, the curves show a pronounced angular variation which demonstrates that B(T) (or T(B)) is field-angle dependent. For instance, at 2.5 K and 0.4 T (violet symbols) the sample is in the normal state at θ = 3° and superconducts (R = 0) at 93°. By raising the temperature from 2 K to 3 K superconductivity is smoothly depressed for all field directions. The data show a striking two-fold symmetry, which is most clearly demonstrated in a polar plot (Fig. 2). We remark, the same two-fold anisotropy is observed in crystals with x = 0.15. In the top panel of Fig. 1 we show R(θ) in the normal state measured in 8 T for x = 0.10. The data have been symmetrized after measuring R(θ) for opposite field polarities to eliminate a small Hall component. R(θ) in the normal state shows the same two-fold symmetry as in Fig. 1a. The variation in R(θ) is small and amounts to 3% in 8 T. The data follow a sin θ dependence, which tells us the variation is due to the classical magnetoresistance related to the Lorentz force F = BI sin θ, where I is the transport current that flows in the basal plane. R(θ) is minimum in the longitudinal case (B || I) and maximum in the transverse case (B ⊥ I).

Figure 1

Angular variation of the resistance of Sr0.10Bi2Se3.

Lower panel: Resistance of Sr0.10Bi2Se3 as a function of angle θ at B = 0.4 T and temperatures between 2.0 K (bottom) and 3.0 K (top). The angle θ = 3° corresponds to B ⊥ I and θ = 93° to B || I as indicated by arrows. The current direction is along the a-axis, with a precision of several degrees. The data are measured with increasing angle, and reproduce when the rotation direction is reversed, apart from a small backlash in the rotator of 2°. Upper panel: R(θ) in the normal state at T = 3.0 K and B = 8 T. The solid line shows R(θ) can be described by a sin θ function.

Figure 2

Polar plot of the resistance of Sr0.10Bi2Se3.

Resistance of Sr0.10Bi2Se3 as a function of angle θ in a magnetic field of 0.4 T and temperatures ranging from 2.0 K to 3.0 K presented in a polar plot. The angle θ = 3° corresponds to B || a* ⊥ I, while θ = 93° corresponds to .

In Fig. 3 we report B(T) for two single crystals measured with the B-field along the orthogonal directions in the hexagonal unit cell. The data points are obtained by measuring the superconducting transition in R(T) in fixed fields, where T is identified by the 50% drop of R with respect to its value in the normal state32. In determining the values of B we did not correct for demagnetization effects, since the demagnetization factors calculated for our crystals are small32. As expected from the data in Fig. 1, we observe a large difference between and , with an in-plane anisotropy parameter of 6.8 (at 1.9 K) and 2.6 (at 0.3 K) for x = 0.10 and x = 0.15, respectively. For both crystals . Obviously, the B ratio γ for the field || and ⊥ to the layers now depends on the field angle and ranges from 1.2 to 3.2 for x = 0.15. In ref. 24 a value for γ of 1.5 is reported, whereas from the data in ref. 23 we infer a value of 1. In the top panels of Fig. 3 we show ρ(B) measured along the a, a* and c axis at T = 2.0 K and T = 0.3 K for x = 0.10 and x = 0.15, respectively. The B(T) values are determined by the midpoints of the transitions to the normal state, and are indicated by open symbols in the lower panels. The agreement between both methods (field sweeps and temperature sweeps) is excellent. For the x = 0.15 sample we see a remarkable broadening for B || a. The initial small increase of ρ(B) between 4 and 6 T is most likely related to a sample inhomogeneity, because a similar tail is also observed in the R(T) data32.

Figure 3

Upper critical field of SrBi2Se3.

Panel (a,b): Resistance of SrBi2Se3 as a function of a, a* and c, for x = 0.10 and 0.15, respectively. The open symbols indicate the midpoints of the transitions to the normal state. Panel (c,d): B obtained for a, a* and c, for x = 0.10 and 0.15, respectively. Solid symbols from midpoints of R(T)-curves in fixed B32. Open symbols from ρ(B) at fixed T. In the experiments for x = 0.15 the crystal was not mounted on the rotator but oriented by eye, which adds some inaccuracy as regards field alignment. The current direction was always along the a-axis, with a precision of several degrees.

In Fig. 4 we show the angular variation of the upper critical field, B(θ). For this experiment the crystals are placed on the rotator and the field is oriented in the basal plane. The data points are obtained as the midpoints of the transitions to the normal state of the R(B) curves measured at temperatures of 2 K for x = 0.10 and of 0.3 K and 2 K for x = 0.15 (see Fig. S732). All data sets show the pronounced two-fold basal-plane anisotropy of B, already inferred from Figs 1 and 2.

Figure 4

Angular variation of B of SrBi2Se3 in the basal plane.

Panel (a,b): Angular variation of B in the basal plane (aa*-plane) for x = 0.10 and 0.15, respectively, at temperatures as indicated. The data are obtained from R(B) measurements at fixed T. The angle θ = 0° corresponds to B || a* ⊥ I and θ = 90° to . The solid black line in panel (b) represents B(θ) for an anisotropic effective mass model with two-fold symmetry and Γ = 3.2 (see text). The a and a* directions in the hexagonal basal plane are defined as in the figure in the upper right corner of panel (b).

Discussion

Having conclusively established the two-fold anisotropy of B in the basal plane, we now turn to possible explanations. A first explanation could be a lowering of the symmetry caused by a crystallographic phase transition below room temperature. However, the powder X-ray diffraction patterns measured at room temperature and T = 10 K are identical (see Fig. S2 in ref. 32). Moreover, the resistivity traces (T = 2–300 K, Fig. S4) and the specific heat (T = 2–200 K, Fig. S8) all show a smooth variation with temperature and do not show any sign of a structural phase transition32. We therefore argue our crystals keep the space group at low temperatures.

A second explanation for breaking the symmetry in the basal plane could be the measuring current itself. Since the current flows in the basal plane it naturally breaks the symmetry when we rotate the field in the basal-plane. Indeed B is largest for B || I and smallest for B ⊥ I. In the latter geometry, and for large current densities, the Lorentz force may cause flux lines to detach from the pinning centers, which will lead to a finite resistance, a broadened R(B)-curve and a lower value of B29. This effect has been observed for instance in the hexagonal superconductor MgB2 by rotating B with respect to I in the basal plane33. For a current density 30 A/cm2, the two-fold anisotropy obtained just below T = 36 K is small, ~8%33. In our transport experiments the current densities are ≤0.4 A/cm2 and we did not detect a significant effect on the resistance when the current density was varied close to T (see Fig. S9 32). Also, when flux flow has a significant contribution, one expects the R(B)-curves for B ⊥ I to be broader than the curves for B || I. However, we observe the reverse (see Fig. 3a,b). Moreover, the anisotropy is still present at T/T = 0.1 and is much larger (of the order of 300%, see Fig. 4) than can be expected on the basis of flux flow. In order to further rule out the influence of the current direction we have investigated B(θ) in the basal plane with the transport current perpendicular to the layers (I || c) and thus keeping B ⊥ I (see Fig. S11, ref. 32). The angular variation of the resistance, measured in this geometry using a two-probe method, is similar to that reported in Fig. 1. Thus the two-fold anisotropy in B is also present for the B-field in the aa*-plane and the current along the c-axis.

Next we address whether the variation of B in the basal plane can be attributed to the anisotropy of the effective mass. Within the Ginzburg-Landau model2734 the anisotropy of B is attributed to the anisotropy of the superconducting coherence length, ξ, which in turn relates to the anisotropy of the effective mass. For a layered superconductor the anisotropy ratio 28. Here m and M are the effective masses || and ⊥ to the layers. In the rhombohedral structure and M = m, where the subscripts a, a* and c refer to the effective masses for the energy dispersion along the main orthogonal crystal axes (i.e. in the hexagonal unit cell). For a field rotation in the aa*-plane is in general isotropic, since . For a 3D anisotropic superconductor the angular variation B(θ) in a principal crystal plane can be expressed as , where Γ = B(90°)/B(0°). To provide an estimate of Γ for Sr0.15Bi2Se3, we compare in Fig. 4b the measured B(θ) with the angular variation in the anisotropic effective mass model (solid line). We obtain B(0°) = 2.3 T, B(90°) = 7.4 T and Γ = 3.2. The effective mass ratio 34 would then attain the large value of 10.2. As we show below, this is not compatible with the experimental Fermi-surface determination.

The Fermi surface of n-doped Bi2Se3, with a typical carrier concentration n ~ 2 × 1019 cm−3 representative for the superconducting SrBi2Se3 crystals2324, has been investigated by the Shubnikov - de Haas effect233536. It can be approximated by an ellipsoid of revolution with the longer axis along the k-axis. A trigonal warping of the Fermi surface due to the rhombohedral symmetry has been detected, but the effect is small: the variation of the effective mass in the basal plane amounts to a few % only35. This also explains why R(θ) in the normal state (Fig. 1a), does not show a 2π/3 periodicity superimposed on the two-fold symmetry induced by the current. Clearly, the two-fold symmetry (Fig. 4), while three fold is expected, and the calculated large ratio using the Ginzburg-Landau model are at variance with the experimental Fermi-surface determination35 and we discard this scenario.

Having excluded these conventional explanations for the rotational symmetry breaking we now proceed to a more exciting scenario. Nagai (ref. 30) and Fu (ref. 31) recently proposed a model for odd parity spin-triplet superconductivity developed in the context of CuBi2Se3, and investigated the experimental consequences of Δ4 pairing in the two-orbital model18. Here, superconductivity is described by an odd-parity two-dimensional representation, E, where the attractive potential pairs two electrons in the unit cell to form a spin triplet, i.e. a vectorial combination of c1↑c2↑ and c1↓c2↓. The indices 1, 2 refer to the two orbitals and the arrows to the spin. The Δ4 state has zero-total spin along an in-plane direction n = (n, n) that is regarded as a nematic director and breaks rotational symmetry. By taking into account the full crystalline anisotropy in the Ginzburg-Landau model, it can be shown that n is pinned to a direction in the basal plane. For , point nodes in the superconducting gap are found along , whereas for two gap minima occur at 31. Our B-data can be interpreted as reflecting a strongly anisotropic superconducting gap function. The superconducting coherence length, ξ, along the main axes can be evaluated from the Ginzburg-Landau relations , and . Here Φ0 is the flux quantum. With the experimental B-values, taken at T/T = 0.1 in Fig. 3d for x = 0.15, we calculate ξ = 19.6 nm, nm and ξ = 5.4 nm. Interpreting ξ as the Cooper-pair size, this implies that the pairing interaction is strongest along the a* and c-axis, and weakest along the a-axis. The observation that can naively be translated to the gap structure consistent with the one predicted for . More recent calculations show that B for the two-dimensional E representation retains the hexagonal symmetry of the crystal lattice, but its symmetry can be lowered to two-fold in the presence of a symmetry breaking field3738. As regards SrBi2Se3 the origin of the symmetry breaking is not clear yet. Possible candidates are sample shape, residual strain and local ordering of Sr atoms. We remark that rotational symmetry breaking in the spin system has been observed by Nuclear Magnetic Resonance (NMR) in the related superconductor CuBi2Se3, which is considered to provide solid evidence for a spin-triplet state39.

Yet another interesting possibility is a self-organized structural stripiness in the optimum for superconductivity due to ordering of Sr atoms in the Van der Waals gaps. This could naturally lead to an anisotropy of B when measured for a current in the basal plane, because of an effective reduced dimensionality. The higher B-values will then be found for B || I along the stripes. On the other hand, for I perpendicular to the layers the basal-plane anisotropy of B is found as well32. This calls for a detailed compositional and structural characterization of SrBi2Se3 by techniques such as Electron Probe Microprobe Analysis (EPMA) or Transmission Electron Microscopy (TEM). Notice that in CuBi2Se3 crystals EPMA has revealed that the Cu concentration shows variations on the sub-mm scale, which gives rise to superconducting islands40. Moreover, a STM study reports an oscillatory behaviour of the Cu pair distribution function due to screened Coulomb repulsion of the intercalant atoms41.

In conclusion, we have investigated the angular variation of the upper critical field of superconducting crystals of SrBi2Se3. The measurements reveal a striking two-fold anisotropy of the basal-plane B. The large anisotropy cannot be explained with the anisotropic effective mass model or the variation of B caused by flux flow. We have addressed two alternative explanations: (i) unconventional superconductivity, with an odd-parity triplet Cooper-pair state (Δ4 pairing), and (ii) self-organized striped superconductivity due to preferential ordering of Sr atoms. The present experiments and results provide an important benchmark for further unraveling the superconducting properties of the new candidate topological superconductor SrBi2Se3.

After completion of this work we learned that rotational symmetry breaking has been observed in two related superconductors, namely in CuBi2Se3 by means of specific heat experiments42 and in NbBi2Se3 by means of torque magnetometry43.

Methods

Sample preparation

Single crystals SrBi2Se3 with x = 0.10 and x = 0.15 were prepared by melting high-purity elements at 850 °C in sealed evacuated quartz tubes, followed by slowly cooling till 650 °C at the rate of 3 °C/hour. Powder X-ray diffraction confirms the space group (see Supplementary Information 32). Laue back-scattering diffraction confirmed the single-crystallinity and served to identify the crystal axes a and a*. Thin bar-like samples with typical dimensions 0.3 × 1.5 × 3 mm3 were cut from the bulk crystal for the transport measurements.

Magnetotransport experiment

Magnetotransport experiments were carried out in a PPMS-Dynacool (Quantum Design) in the temperature range from 2 K to 300 K and magnetic fields up to 9 T and in a 3-Helium cryostat (Heliox, Oxford Instruments) down to 0.3 K and fields up to 12 T. The resistance was measured with a low-frequency ac-technique in a 4-point configuration with small excitation currents, I, to prevent Joule heating (I = 0.5–1 mA in the PPMS and 100 μA in the Heliox experiments). The current was applied in the basal plane along the long direction of the sample. For in-situ measurements of the angular magnetoresistance the crystals were mounted on a mechanical rotator in the PPMS and a piezocrystal-based rotator (Attocube) in the Heliox. The samples were mounted such that the rotation angle corresponds to B ⊥ I. Care was taken to align the a-axis with the current direction, but a misorientation of several degrees can not be excluded.

Additional Information

How to cite this article: Pan, Y. et al. Rotational symmetry breaking in the topological superconductor SrBi2Se3 probed by upper-critical field experiments. Sci. Rep. 6, 28632; doi: 10.1038/srep28632 (2016).