Time-reversal symmetry-breaking superconductivity in epitaxial bismuth/nickel bilayers

See allHide authors and affiliations

Science Advances  31 Mar 2017:
Vol. 3, no. 3, e1602579
DOI: 10.1126/sciadv.1602579


Superconductivity that spontaneously breaks time-reversal symmetry (TRS) has been found, so far, only in a handful of three-dimensional (3D) crystals with bulk inversion symmetry. We report an observation of spontaneous TRS breaking in a 2D superconducting system without inversion symmetry: the epitaxial bilayer films of bismuth and nickel. The evidence comes from the onset of the polar Kerr effect at the superconducting transition in the absence of an external magnetic field, detected by the ultrasensitive loop-less fiber-optic Sagnac interferometer. Because of strong spin-orbit interaction and lack of inversion symmetry in a Bi/Ni bilayer, superconducting pairing cannot be classified as singlet or triplet. We propose a theoretical model where magnetic fluctuations in Ni induce the superconducting pairing of the Embedded Image orbital symmetry between the electrons in Bi. In this model, the order parameter spontaneously breaks the TRS and has a nonzero phase winding number around the Fermi surface, thus making it a rare example of a 2D topological superconductor.

  • Superconductivity
  • MBE
  • 2D
  • time-reversal symmetry
  • Topological


Unconventional superconductors, characterized by spontaneous breaking of additional symmetries other than U(1) gauge symmetry, often have novel properties and reveal new physics. Spontaneous time-reversal symmetry breaking (TRSB) in a superconducting state is one of the most fascinating phenomena. Such superconductors may host unusual particles called Majorana fermions (1), much discussed for topological quantum computing (2). However, up to now, spontaneous TRSB has been observed only in a few three-dimensional (3D) superconductors with bulk inversion symmetry, such as Sr2RuO4, UPt3, and URu2Si2 (35). The possibility of spontaneous TRSB has not been investigated experimentally, so far, for a 2D surface superconductivity at interfaces between nonsuperconducting materials (6), which are currently attracting a lot of interest for possible device applications. These 2D superconducting systems typically lack inversion symmetry because of the difference between the top and bottom surroundings of the layer. Because of the spin-orbit interaction, the spin degeneracy of the electronic states is lifted in these noncentrosymmetric materials, thus making the usual singlet-triplet classification of superconducting pairing inapplicable (7). Recent progress in the fabrication of high-quality heterostructures enables artificial realization of these exotic superconducting states. Here, we report an observation of spontaneous TRSB in the superconducting state of a 2D epitaxially grown bilayer made from the nonsuperconducting elements Bi and Ni.

The bulk rhombohedral crystalline Bi, which is the heaviest nonradioactive element in the periodic table and a semimetal with an extremely long mean free path, is not superconducting at temperatures down to 50 mK (8). Thus, it was a complete surprise when tunneling experiments (9, 10) identified a superconducting transition at Tc = 3.6 K in Bi thin films evaporated on a thin seed layer of ferromagnetic Ni. It is unknown what the symmetry of this superconducting pairing is and, in particular, whether it breaks the time-reversal symmetry (TRS). Because of the recent progress in epitaxial film technology, it became possible to fabricate high-quality superconducting Bi/Ni films (11) suitable for optical measurements of the polar Kerr effect (PKE). The PKE can be used for experimental detection of spontaneous TRSB at a superconducting transition, as it has been demonstrated for Sr2RuO4, UPt3, and URu2Si2 (35). The PKE measures the optical phase difference between two opposite circular polarizations of light reflected on a sample surface, thus giving information about the TRS of the system. For optical frequencies ω, which are much higher than the superconducting gap Δ, the Kerr signal θK is expected (12) to be proportional to (Δ/ω)2 and is generally less than 1 μrad (10−9 rad) (3, 4). To measure such a small θK, we used a custom-made ultrasensitive loop-less fiber-optic Sagnac interferometer (see fig. S1) (12, 13). The geometry and design of the Sagnac interferometer (12) guarantee that all the time reversal–invariant effects (such as optical birefringence, temperature fluctuations, and mechanical vibrations) cancel out, and a nonzero signal θK appears only if the measured sample breaks the TRS. This method (14, 15) probes the sample at an optical penetration depth of a few tens of nanometer and is ideally suited for studying thin films.


The Bi(110)/Ni(001) samples were fabricated using molecular beam epitaxy on MgO(001) substrate (see fig. S2) and then mounted in an optical cryostat, with the optical probe beam incident normally on the sample surface for the PKE measurements, as shown in Fig. 1A. Four gold wires were attached to the sample surface with silver epoxy for transport measurements. A study of thickness dependence (see fig. S4) shows that an increase in the thickness of the face-centered cubic (fcc) Ni layer suppresses Tc, whereas an increase in Bi thickness restores Tc back to 4.1 K. These observations (11) suggest that although the presence of Ni is necessary for providing the pairing interaction in Bi, superconductivity does not develop at the direct Bi/Ni interface because of the strong ferromagnetism of Ni but, instead, develops at the opposite surface of Bi facing the vacuum. A further increase of Bi thickness starts to suppress Tc (11), as expected for going toward the limit of bulk Bi. Also, there is an optimal Bi layer thickness for maximal Tc. In this experiment, we choose a Bi thickness of about 20 nm to minimize the Kerr signal from the underlying Ni layer while maintaining the optimal Tc = 4.1 K.

Fig. 1 Structure and Kerr signal of a Bi (25 nm)/Ni (2 nm) sample.

(A) Left: Side view of the sample structure with a transmission electron microscopy image. Right: Schematic diagram of Kerr rotation measurement. (B) Kerr angle (purple) and resistance (green) measured on the Bi side for ZF cooldown, showing the onset of θK at Tc = 4.1 K. The pink dashed line is a guide to the eye with the form of 1 − (T/Tc)2. (C) Kerr signal measured on the Ni side at 1600-Oe applied perpendicular magnetic field (red) or at ZF (blue). In either case, no change of the Kerr signal was observed across Tc = 4.1 K.

Our experimental setup allows us to perform optical measurements on both the Bi and Ni sides of the bilayer, as shown in Fig. 1 (B and C). The main results of the paper are the measurements on the Bi side shown in Fig. 1B. The Bi (25 nm)/Ni (2 nm) sample was cooled at zero magnetic field (<0.1 Oe) to the base temperature of the cryostat and then measured upon warming in zero field (ZF). The onset of the ZF Kerr signal is observed at 4.1 K, which coincides with the superconducting transition, as signaled by the initial drop of resistance, thus demonstrating that spontaneous TRSB is due to superconductivity. The Kerr signal θK becomes larger at lower temperatures, reaching 120 nrad at T = 3 K. The dashed line is a guide to the eye of the form of 1 − (T/Tc)2 (16). Each Kerr angle θK data point represents an average of 1000 s, with the error bar denoting the statistical uncertainty of less than 10 nrad. The systematic source of error/uncertainty is mainly the long-term drift of the interferometer, which is estimated to be 5 nrad. The Sagnac interferometer operates at an optical power of 2 μW to prevent any optical heating (3, 4) of the sample.

One may ask whether the observed onset of the Kerr signal at 4.1 K on the Bi side of the sample could somehow originate from the ferromagnetic Ni on the other side. We believe this is not the case because of the following reasons. First, the 400-K ferromagnetic Curie temperature of Ni is two orders of magnitude greater than the 4.1-K onset temperature of PKE on the Bi side; thus, it is unlikely that the two effects are related. Second, the magnetic moment of Ni is almost parallel to the interface surface, whereas the PKE is sensitive only to the perpendicular component. Third, the opacity and thickness of the Bi layer greatly suppress an optical pickup of a signal originating from the Ni side. Nevertheless, to rule out the possibility that a yet-unknown magnetic transition in the underlying Ni layer at 4.1 K generates the Kerr signal on the Bi side, we have directly measured the PKE on the Ni side by flipping the sample, as shown in Fig. 1C. Thin Ni films are known to have an almost in-plane magnetic anisotropy; hence, the magnetic moment of Ni is almost in the plane, resulting in near-zero polar Kerr signal. This is shown by the blue curve in Fig. 1C, where a temperature-independent θK of 0.75 μrad is observed without any sign of transition at 4.1 K. When a strong magnetic field of 1600 Oe is applied perpendicularly to the interface, the majority of the magnetic moment of Ni is out of the plane. In addition, the corresponding perpendicular component produces a Kerr signal of 130 μrad, as shown by the red curve in Fig. 1C. These indicate a magnetic anisotropy angle of less than 0.33° from the in-plane direction, which is common for thin Ni films (10). During the Kerr measurement on the Bi side, because only a tiny fraction of the light could penetrate into the underlying Ni layer, the Kerr signal is less than 5 μrad, or 4% of the Kerr signal from the Ni side. Therefore, if the 120-nrad Kerr signal in Fig. 1B was due to an unknown magnetic transition in Ni at 4.1 K, it would necessarily require an increase of the Kerr signal on the Ni side by 3 μrad at a temperature below 4.1 K, which is nonexistent in the blue curve in Fig. 1C up to a precision of 0.1 μrad.

Because of spontaneous symmetry breaking, the sign of the Kerr signal observed below Tc may be positive or negative, but it can be controlled by applying a small external “training” magnetic field when cooling through the superconducting transition. Such a training effect has been observed in other TRSB superconductors (35), and we also find it in Bi/Ni. As shown in Fig. 2 (A and B), a Bi (20 nm)/Ni (2 nm) sample was cooled down to the lowest temperature in a magnetic field of +190 or −70 Oe perpendicular to the film, and then the Kerr signal was measured on warm-up in zero magnetic field. We observed that the sign of the Kerr signal on warm-up corresponds to the sign of the training field applied during cooldown, thus indicating that the TRSB order parameter couples to the external magnetic field. We have repeated these measurements with training fields up to 230 Oe and found the absolute values of θK at 2.8 K to be identical within the 20-nrad measurement uncertainty. The fact that |θK| is independent of the cooling field magnitude implies that the trapped superconducting flux is not responsible for the observed θK, because the number of trapped vortices would scale with the cooling field. The applied training magnetic fields are also too small to alter the small out-of-plane component of an Ni magnetic moment.

Fig. 2 Chirality training and domain size estimation of a Bi (20 nm)/Ni (2 nm) sample.

(A) ZF warm-up data after cooling down in +190-Oe magnetic field. (B) ZF warm-up data after cooling in −70-Oe magnetic field. Pink dashed lines in (A) and (B) are guides to the eye with the form of 1 − (T/Tc)2. (C) Kerr effect θK measured during ZF warm-up, after cooling down in ZF. Inset: Random (2D) chiral domains under the optical spot. (D) SD σ (θK) between experiments that contain two random contributions: σ0 due to chiral domains and σapp from the apparatus.

For a cooldown in zero magnetic field, we expect that chiral domains with positive (+) and negative (−) signs of the Kerr signal should coexist in the sample. If the optical spot is positioned within a single chiral domain, the full Kerr signal would be observed, which we believe is the case for the measurement performed without a training field (Fig. 1B). However, more often, the measured θK is an average of chiral domains under the focused optical spot of the diameter w ≈ 2 μm, with a value between +θ0 and −θ0, where θ0 = 60 nrad is the value for a single domain (that is, after being trained). This is observed in 11 measurements of the Bi (20 nm)/Ni (2 nm) sample (Fig. 2C). The standard deviation σ (θK) between these measurements has two contributions: σ0 due to random chiral domains and σapp from the apparatus uncertainty. Above Tc, there is no TRSB, and σ (θK) remains constant at σapp = 18 nrad. Below Tc, the TRSB occurs, and σ (θK) increases toward σ3K = 30 nrad at T = 3 K (Fig. 2D). Therefore, σ0 can be estimated to be 24 nrad at T = 3 K. By assuming 2D domains with the average size d, we estimate that d ≈ w × σ00 ≈ 0.8 μm.

To explain the experimental results, we propose a minimal model involving electronic surface states of Bi that interact with fluctuating magnetic moments of Ni (technical details are given in the Supplementary Materials). There are indications (17) that electronic states in the interior of a Bi film are gapped, and the only metallic states are the surface states. Thus, we only consider the Bi surface exposed to vacuum, because the surface exposed to Ni does not contribute to superconductivity, according to Gong et al. (11). The (110) surface of Bi may contain multiple Fermi pockets (1822); however, we only consider the largest pocket enclosing the surface Brillouin zone center Embedded Image. Because of the strong Rashba-type spin-orbit coupling, the electron spin and momentum are locked on this pocket, as shown in Fig. 3A. To simplify the presentation, we assume rotational symmetry in the (x, y) plane of the surface (anisotropy is discussed at the end) and use the standard Rashba Hamiltonian Embedded Image, where k and σ are the momentum and spin operators for 2D electrons, respectively. The Zeeman term −Mσ due to the in-plane magnetization M of the Ni can be eliminated by the gauge transformation Ψ = e−iwr Ψ′ of the electron fields, which is equivalent to a redefinition of the electron momentum p = kw and a shift of the Fermi surface center from Embedded Image to Embedded Image by the vector Embedded Image, as shown in Fig. 3A. Then, we consider the superconducting pairing of the electrons with the opposite momenta p and −p relative to Embedded Image (Fig. 3B), which corresponds to a pairing with the total momentum 2w relative to Embedded Image.

Fig. 3 Shifted Fermi surface and pairing symmetry of Bi electrons.

(A) The Fermi circle of the electron surface states in Bi. The original dashed blue Fermi circle centered at Embedded Image is shifted by the vector w to the solid blue circle (shaded area) centered at Embedded Image because of the in-plane magnetization M produced by Ni. The electron momentum p is measured from Embedded Image and characterized by the azimuthal angle φp. The green arrows show spin polarization locked to the momentum direction. (B) Superconducting pairing of the electrons with opposite spins and opposite momenta p and p. The TRSB condensate has the total angular momentum Jz = ±2, corresponding to the d ± id pairing, as indicated by the red double-curved arrows. A weak training magnetic field can select one of the two degenerate states.

Because of the strong spin-orbit interaction and the lack of inversion symmetry (signified by the direction perpendicular to the surface), the electron states are nondegenerate for each momentum p. Symmetry classification of 2D superconducting pairing in this case (23) differs substantially from the conventional, spin-degenerate case. In particular, classification into singlet and triplet pairing is not appropriate (24). With the spin quantization axis taken along the direction, the electron state of the momentum p is the spinor Embedded Image, where ϕp is the azimuthal angle of the vector p in Fig. 3A. The time-reversed state Embedded Image produced by the time-reversal operator Embedded Image, where Embedded Image is the complex conjugation and iσy is the spinor metric tensor, has the momentum –p; thus, Embedded Image, where Embedded Image is a phase factor. Because Θ2 = −1 for spin 1/2, the phase factor has the property ηp = − ηp (which can also be checked using ϕp = ϕp + π) and cannot be eliminated by any gauge transformation. The corresponding second-quantized operators also involve this phase factor: Embedded Image = Embedded Image. Following the study by Samokhin (23), the superconducting condensate Embedded Image is introduced for the fermion fields that are time-reversal partners of each other. An important observation is that f(p) is an even function of p: Embedded Image, where fermion anticommutation was used. Thus, only even values of m are permitted in the expansion Embedded Image over circular harmonics for a rotationally invariant system (because ϕp = ϕp + π). The simplest case m = 0, studied by Gor’kov and Rashba (24), represents a pairing that is rotationally and time reversal–invariant. The px ± ipy pairing with the odd m = ±1, often discussed for Sr2RuO4, is forbidden in the nondegenerate case.

To explain spontaneous TRSB in Bi/Ni bilayers, we propose that the superconducting pairing in this system takes place with m = ±2, where the sign is either selected spontaneously or controlled by a small training magnetic field, as shown in Fig. 3B. The pairing states with m = ±2, which can be also labeled as Embedded Image, are the time-reversal partners of each other. As shown in the study by Samokhin (23), these states have two chiral Majorana edge modes propagating around the sample edge, either clockwise or counterclockwise. In the conventional spinor basis, the superconducting condensate Fαβ(p) = 〈ΨpαΨpβ〉 is a 2 × 2 matrix labeled by spin indices α, β = ↑, ↓. Although the spin Sz and orbital angular momentum Lz have different values for different components of this matrix, representing a mix of singlet and triplet pairing, the total angular momentum has the same value Jz = Sz + Lz = m = ±2 for all matrix components. Similar conclusions were obtained by Vafek and Wang (25) using renormalization group for repulsive interaction between electrons. In a less-convenient representation used in the studies by Farrell and Pereg-Barnea (26) and Sergienko and Curnoe (27), the pairing function 〈ψpψp〉 = ηpf(p) has the winding number m − 1, with the values 1 and −3 (corresponding to p- and f-wave pairing) for the time-reversed pairings with the total angular momentum Jz = m = ±2.

It was shown in the study by Scheurer (28) that phonons cannot possibly produce the TRSB superconductivity. Thus, we propose that the superconducting pairing of the surface electron states in the Bi layer is mediated by magnetic fluctuations in the Ni layer. This scenario is supported by the observation (11) that superconductivity exists only for very thin Ni layers and is suppressed with the increase of Ni thickness, which would reduce magnetic fluctuations. We particularly focus on the out-of-plane fluctuations of the Ni magnetic moments [the effect of in-plane fluctuations is discussed in the study by Kargarian et al. (29)]. The effective interaction in the Cooper channel between electrons with the momenta p and p′ is given by the properly symmetrized product of the magnon propagator D(pp′) and the spin form factor Embedded Image (30). The dominant contribution to the m = ±2 pairing channel comes from the product of the first circular harmonics of D and Λs but is somewhat reduced by the product of the second harmonic of D and the zeroth harmonic of Λs.

Finally, we discuss the effects of crystal anisotropy. The (110) surface of Bi has a low-symmetry group D1 with a mirror plane perpendicular to the surface (18). Consequently, there are only two symmetry representations, A1 and A2, for the superconducting pairing (23), which are symmetric and antisymmetric with respect to reflection. The antisymmetric representation A2 has nodes and is even in p; thus, its wave function is pxpy of the dxy type. The symmetric representation A1 can be of the s or Embedded Image type, which belongs to the same representation. Because two different representations, A1 and A2, are available, it is possible to make two degenerate time-reversal breaking combinations Embedded Image out of them. However, because the A1 and A2 representations are not degenerate, their transition temperatures are generally different. Thus, we come to the conclusion that there must be two superconducting transitions in the Bi/Ni bilayers, as in UPt3 (4). We expect that the nodal superconducting state of the dxy type would develop at a higher transition temperature and that the Embedded Image state would develop at a lower temperature, spontaneously breaking the TRS by the Embedded Image combination and eliminating the nodes. If the crystal anisotropy is weak, the difference between transition temperatures may be small, which may be the reason why the two superconducting transitions have not been resolved so far, but it is an important direction for future experimental studies.


In conclusion, we report observation of the PKE in the Bi/Ni epitaxial bilayers, which indicates spontaneous TRSB in the superconducting state. A combination of strong spin-orbit coupling and noncentrosymmetry makes this 2D interfacial system distinctly different from the bulk TRSB superconductors (35). To describe the TRSB state, we propose a model where magnetic fluctuations in Ni induce superconducting pairing in Bi, with the superconducting phase winding number around the Fermi surface equal to ±2. Thus, Bi/Ni is likely a rare example of a 2D topological superconductor with two chiral edge states moving in the same direction (23).


Setup of the Sagnac interferometer

A schematic of the Sagnac interferometer used in this experiment is shown in fig. S1. In a fiber-optic Sagnac interferometer (13), light from a 1550-nm broadband source passed through a polarization-maintaining (PM) circulator and was polarized. A half-wave plate was then used to rotate the polarization axis to 45° with respect to the axis of a phase shift modulator. This modulator added a time-varying phase shift to the light. However, the amplitude of the phase shift and the optical path length were different for the in-plane and out-of-plane polarizations of light. After passing through the modulator, the in-plane and out-of-plane polarizations were no longer coherent with each other, and each polarization beam had a time-varying phase shift of different amplitudes. The light was then routed via a PM fiber to optics mounted on a piezo stage capable of submicrometer step size in both the x and y directions. After exiting the PM fiber, the light was passed through a quarter-wave plate that converted the orthogonal, linear polarizations into left- and right-circular polarizations. Finally, it was focused onto the sample surface. The sample rested inside an optical cryostat with a base temperature of 2.5 K, which contained an optical viewport on the top, allowing the light to pass from the focusing element to the sample surface. Upon reflection from the sample, there was a phase shift between the two polarizations that was twice the size of the Kerr signal. After passing through the quarter-wave plate, the polarization axis of the two beams had been swapped. On the return trip, each beam traveled through the system with a polarization orthogonal to that of its outgoing trip. The net result is that each beam traveled along the exact same optical path but in opposite direction over the round trip. Once each beam passed through the phase shift modulator again, the two beams were once again coherent but had phase differences due to the Kerr effect (2θK) as well as the phase shift modulation. The two beams interfered, resulting in an elliptically polarized beam, the in-plane component of which was routed to a photodetector. If ωm was chosen so that τωm = π, lock-in detection could be used to determine θK by comparing the detected signal at ωm and 2ωm. The dc power received by the photodetector is a measure of the reflectivity of the sample, and the ratio between the signal at 2ωm and the dc signal is a measure of optical anisotropy, as described in the study by Xia et al. (13).

Procedure of epitaxially growing samples

All the samples were epitaxially grown in ultrahigh vacuum (UHV), with a base pressure 6 × 10−8 Pa. We first chemically cleaned the MgO(001) substrates before putting them into the UHV chamber. After transfer, we annealed them to 750 K and incubated them at this temperature for 70 min. Figure S2A shows the reflection high-energy electron diffraction (RHEED) patterns for the MgO substrate. The Ni layer was deposited at 300 K, whereas the Bi layer was deposited at 110 K. Figure S2B shows the RHEED patterns of the topmost Bi layer taken at room temperature. The sharp lines indicate that the sample is single crystalline and that the quality is very high, which can also be verified by the cross-sectional scanning transmission electron microscopy image in Fig. 1A.


Supplementary material for this article is available at

Experimental methods and supplementary data

Theory of superconducting pairing in epitaxial bismuth/nickel bilayers

fig. S1. Schematic of Sagnac interferometer.

fig. S2. RHEED patterns of the substrate and the sample.

fig. S3. Supplementary data of the Bi (25 nm)/Ni (2 nm) sample.

fig. S4. Dependence of the critical temperature Tc on the thickness of the Bi and Ni layers.

fig. S5. Kerr signal of the Bi (40 nm)/Ni (2 nm) sample.

fig. S6. Optical signal from the MgO substrate.

References (3141)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: V.M.Y. thanks K. Samokhin for the discussion. Funding: This study was supported by the NSF (grant DMR-1350122). X.G., D.Y., H.Z., and X.J. acknowledge support from the National Natural Science Foundation of China (grant nos. 11434003 and 11421404) and the National Basic Research Program of China (grant no. 2015CB921402). M.K. and V.M.G. acknowledge support from the U.S. Department of Energy–Basic Energy Sciences (DESC0001911) and the Simons Foundation. Author contributions: X.G. took the responsibility of conducting the measurements and analyzing the data. A.S. assisted in the experiment. D.Y. and H.Z. prepared the samples. X.J. developed the Sagnac interferometer. X.J. and J.X. designed the project. M.K. and V.M.Y. conducted the theoretical calculation. J.X., V.M.Y., and V.M.G. wrote the paper. All of the authors contributed to the analysis of the results and the preparation of the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.
View Abstract

Navigate This Article