Nonreciprocal charge transport at topological insulator/superconductor interface.

Topological superconductor is attracting growing interest for its potential application to topological quantum computation. The superconducting proximity effect on the topological insulator surface state is one promising way to yield topological superconductivity. The superconductivity realized at the interface between Bi2Te3 and non-superconductor FeTe is one such candidate. Here, to detect the mutual interaction between superconductivity and topological surface state, we investigate nonreciprocal transport; i.e., current-direction dependent resistance, which is sensitive to the broken inversion symmetry of the electronic state. The largely enhanced nonreciprocal phenomenon is detected in the Bi2Te3/FeTe heterostructure associated with the superconducting transition. The emergent nonreciprocal signal at low magnetic fields is attributed to the current-induced modulation of supercurrent density under the in-plane magnetic fields due to the spin-momentum locking. The angular dependence of the signal reveals the symmetry of superconductivity and indicates the existence of another mechanism of nonreciprocal transport at high fields.

Introduction

Topological superconductivity (TSC) is extensively sought for in contemporary condensed matter physics[1-3]. In a conventional superconductor, spin-up and spin-down electrons form a spin-singlet Cooper pair, and even-parity s-wave superconductor is realized. However, in a spinless superconductor, i.e., paired systems with only one active fermionic species, an odd-parity p + ip-wave superconductor is induced because of Pauli’s exclusion principle. Such a superconductor is topologically non-trivial and supports Majorana zero modes at the boundary or at the vortex core, which can act as building blocks of topological quantum computation owing to their non-abelian statistics[1-3]. To realize a non-degenerate spinless mode, one has to fight against the Kramers degeneracy: One way to deal with the Kramars degeneracy is time-reversal symmetry breaking as discussed in superfluid 3He-A phase[4] and superconducting Sr2RuO4[5]. Another strategy is to utilize spin-splitting by inversion symmetry breaking as proposed in the synthetic topological superconductors, such as superconducting–proximity coupled Rashba wires[6-8] and topological insulator (TI) surface states[9-14]. Spin splitting caused by relativistic spin–orbit coupling and the inversion symmetry breaking yields a spin non-degenerate band, as exemplified by the spin–momentum-locked topological surface state. In these systems, the understanding of the superconducting property on the spin-splitted band structure is a prerequisite for the realization of TSC.

Recently, nonreciprocal charge transport measurement has been developed as a new probe to detect the spin-splitting in inversion symmetry-broken systems[15-24]. Reflecting the asymmetric nature of the electronic band structure with spin splitting, the electrical resistivity is expected to vary depending on the current and magnetic field direction. For example, nonreciprocal resistance is experimentally observed in materials with Rashba-type spin splitting, such as a bulk Rashba semiconductor[15], interfaces[16-19], and TI surfaces[20-22]. Under the in-plane magnetic field perpendicular to the current, the resistance value depends on the spin accumulation direction, either parallel or antiparallel to the external magnetic field[15,18-22]. Hence, this method works as an alternative to angle-resolved photoemission spectroscopy (ARPES)[13] and scanning tunneling microscope (STM)[14] to discuss the effect of spin splitting via transport measurement. Furthermore, nonreciprocal transport is also observed in non-centrosymmetric superconductors, such as a chiral nanotube[23] and a transition metal dichalcogenide[24], demonstrating that this probe is also applicable to inversion symmetry broken superconductors.

In this study, we target the emergent superconductivity at the interface between Bi2Te3 and FeTe[25]. The interface superconductivity is induced when a three-dimensional TI Bi2Te3[26] is grown on a parent compound of iron-based superconductor FeTe as schematically shown in Fig. 1a. Since the topological surface state locates at the interface, the interaction between them makes this system a strong candidate of TSC. However, it has remained unclear how the spin–momentum locking of the surface state affects the superconducting property[25]. To study the interplay between them, we here employ nonreciprocal charge transport. The nonreciprocal transport is especially effective to this system because surface-sensitive methods, such as ARPES and STM is not applicable to such a buried interface. The large nonreciprocal transport is observed to be associated with the superconducting transition of Bi2Te3/FeTe. This is revealed to originate from the current-induced-modulation of supercurrent density due to spin–momentum locking, which represents the close connection between superconductivity and topological surface state.

Fig. 1

Nonreciprocal transport in a Bi2Te3/FeTe heterostructure. a The illustration of a Bi2Te3/FeTe heterostructure (Bi2Te3: blue, FeTe: green). The interface of Bi2Te3 and FeTe shows superconductivity (yellow region). The nonreciprocal transport is measured with current I along x-direction under the in-plane magnetic field B (along y-direction) perpendicular to the current. The mutual interaction between spin–momentum-locked topological surface state and superconductivity (as represented by ΔSC) is expected at the interface as expressed in momentum space in the right circle. σ represents the spin direction of the electron. b Temperature dependence of resistance in FeTe(18 nm)/CdTe(100) (green), Bi2Te3(15 nm)/InP(111) (blue) and Bi2Te3(15 nm)/FeTe(18 nm)/CdTe(100) (red) thin films. The superconductivity and topological surface state appear at the interface between Bi2Te3 and FeTe as denoted by the yellow allow. c The temperature dependence of resistance in Bi2Te3/FeTe heterostructure under in-plane magnetic field B = 0, 2, 4, 6, 8, 10, 12, and 14 T. The resistance is measured under I = 1 μA. d The magnetic field dependence of R2 in Bi2Te3/FeTe heterostructure at T = 12 K (blue, normal state) and at T = 9.5 K (light blue, below the superconducting onset temperature) measured with a current magnitude of I = 200 μA

Result

Sample and basic characterization

We grew the heterostructure of Bi2Te3/FeTe using molecular-beam epitaxy (MBE) on CdTe(100) substrate. Van der Waals nature of Bi2Te3 and FeTe allows the sharp interface despite the difference in the lattice symmetry. The scanning transmission electron microscopy and energy-dispersive x-ray spectroscopy images demonstrate the high quality of the interface without discernable interdiffusion (see Supplementary Fig. 1 and Supplementary Note 1). Figure 1b shows the temperature dependence of resistance for the respective samples, FeTe(18 nm), Bi2Te3(15 nm), and Bi2Te3(15 nm)/FeTe(18 nm) bilayer thin films. FeTe thin film shows a semi-metallic behavior with a kink at around 50 K associated with the transition to the bicollinear antiferromagnetic state[27], while Bi2Te3 shows a metallic behavior due to its electron-doped nature. The resistance of Bi2Te3/FeTe behaves as parallel conduction of these two materials except for the resistance drop to zero at around Tc0 = 10.7 K (Fig. 1c), indicating the appearance of superconductivity at the interface[25]. The fitting of the temperature dependence of resistivity with the Berezinskii–Kosterlitz–Thouless (BKT) transition[28-30] and the jump in the power law of current–voltage characteristics (Supplementary Fig. 2 and Supplementary Note 2) verifies the two-dimensional nature of superconductivity[25], where the binding of the vortex–antivortex pairs realizes the zero-resistance state. Because of the two-dimensional nature of superconductivity, the superconducting order parameter appears at Tc0, but finite resistance exists in the intermediate state because of the current-induced motion of vortex and antivortex. When vortex and antivortex form a pair at TBKT, zero-resistance state, i.e., the superconducting state is realized.

Nonreciprocal transport in Bi2Te3/FeTe

If inversion symmetry of the electronic system is broken, the electrical voltage V under the in-plane magnetic field is phenomenologically described up to a linear order in the magnetic field B and second order of current I as[16]Here, we define the direction normal to the plane as z-axis and γ is a coefficient representing the nonlinear and nonreciprocal transport. When the current is perpendicular to the in-plane magnetic field, Eq. (1) becomes V = R0I + γR0BI2. If the second term is finite, the resistance value depends on the current direction, namely the nonreciprocal transport appears. Since the second term is proportional to I2, measurements of the second harmonic resistance under ac current can detect the nonreciprocal transport (see Supplementary Note 3 for the derivation). We fabricated a Hall bar of 100 μm width and measured R2 with a lock-in amplifier[18-24] under the ac current of I = 200 μA (i = 2 A/m in current density) unless otherwise noted. Figure 1d displays the magnetic field dependence of R2 in the normal and the intermediate states. In the normal state (T = 12 K), R2 is almost zero. On the other hand, B-linear R2 appears at T = 9.5 K in the intermediate region as the superconducting order parameter develops.

In contrast to the first harmonic resistance R, which is independent of the current magnitude, R2 is proportional to the current as shown in Fig. 2a, b, which is consistent with Eq. (1). Thus, the γ value can be derived from the slope of the magnetic field dependence of R2/R as shown in Fig. 2c. To see the correspondence between the temperature dependence of R and γ value (Fig. 2d, e), we categorize them into three regions: normal, intermediate, and superconducting region. In the normal region, R2 is not observed within the measurement noise level. In the superconducting region, R2 vanishes with the disappearance of R (see Supplementary Fig. 3 and Supplementary Note 4). On the other hand, γ becomes finite in the intermediate region. Especially, γ value shows a divergent behavior towards the BKT transition temperature, TBKT, since R2 drops to zero more slowly than R. The maximum of γ is 6.5 × 10−3 T−1 A−1 m at T = 6.9 K, which is more than five orders of magnitude larger than the reported value of γ = 6 × 10−8 T−1 A−1 m in TI without superconductivity[22]. The large enhancement of γ associated with the superconducting transition can be qualitatively explained as follows: The energy scale governing the system changes from the Fermi energy EF (~a few hundred meV) to the much smaller superconducting gap ΔSC (~1 meV) along with the superconducting transition. Consequently, the relative energy scale of spin–orbit interaction and magnetic field gets larger, resulting in the enhancement of nonreciprocal transport[24].

Fig. 2

Current-magnitude and temperature dependence of nonreciprocal transport. a The current-magnitude dependence of second harmonic resistance R2ω at T = 9.5 K, and B = 0.5 K deduced from the data shown in b. The black dotted line is the fitting line. b The magnetic field dependence of R2ω at T = 9.5 K measured under I = 40, 80, 120, 160, and 200 μA. c The magnetic field dependence of R2ω/Rω measured under I = 200 μA at T = 6.9, 7.2, 7.5, 7.8, 8.1, 8.5, 9, 9.5, and 10 K. d The temperature dependence of resistance measured under I = 200 μA. The black curve is the fitting of the Berezinskii–Kosterlitz–Thouless (BKT) transition using Halperin–Nelson formula, , where R0 and b are material parameters. Tc0 is the temperature at which the finite amplitude of the order parameter develops and TBKT is the BKT transition temperature. The fitting gives the values, Tc0 = 10.7 K and TBKT = 6.0 K. The blue, green, and red regions correspond to normal, intermediate, and superconducting regions, respectively. e The temperature dependence of γ-value measured under I = 200 μA derived from c. The red points are the measurement on Bi2Te3(15 nm)/FeTe(18 nm) sample (denoted as BT(15 nm)/FT). The green point are the measurement on Bi2Te3(1.5 nm)/FeTe(18 nm) sample (denoted as BT(1.5 nm)/FT). Note that all the measurements, except for Fig. 1b, and the green curve of Fig. 2e, are done on the BT(15 nm)/FT sample. The purple curve is the fitting of the red points with the formula γ = β, where β = 5.3 × 10−3 T−1 A−1 m. Note that the BKT model and the fitting is valid only at around TBKT, which is represented by the solid purple curve. The purple dotted curve is out of the applicable range of theory

Surface state origin of nonreciprocal transport

Because of the following two reasons, it is natural to guess that topological surface state plays a vital role in nonreciprocal transport: (i) spin–orbit coupling is an essential source of nonreciprocal transport[31] and (ii) topological surface state with strong spin–orbit coupling coexists with superconductivity at the interface where nonreciprocal transport appears. To substantiate this argument, we study the nonreciprocal transport of Bi2Te3(1.5 nm)/FeTe(18 nm) with much thinner Bi2Te3 thickness. Although the normal state resistance and the BKT transition temperature is almost the same in the two samples (Supplementary Figs 4 and 5, and Supplementary Notes 5 and 6), the γ-value for Bi2Te3(1.5 nm)/FeTe(18 nm) is about an order of magnitude smaller than that of Bi2Te3(15 nm)/FeTe(18 nm) as shown in the green points in Fig. 2e. This dramatic change is attributed to the suppressed spin polarization of the surface state in the ultrathin limit because of the hybridization between the top and surface states[32,33]. Thus, the large suppression of the signal points to the surface state origin of nonreciprocal transport.

On the basis of the above observation, we theoretically treat the superconducting-proximity-coupled surface state to discuss the temperature dependence of γ. A recent theoretical study[34] suggests that spin–momentum locking of the surface state causes the renormalization of supercurrent density under the current in superconducting-proximity-coupled TI. As a result, the BKT transition temperature is modulated as T′BKT = TBKT(1 + αBI). Interestingly, this equation means that resistance in this system can decrease for a specific current direction, different from the conventional superconductors, where resistance increases by the excitation with large current. By substituting T′BKT to the standard BKT model which is valid at around BKT transition temperature, we getHere, we employed perturbative expansion about αBI, which is justified because R2 is several orders of magnitude smaller than R as shown in Fig. 2c. Comparing the result with Eq. (1), the γ-value is expected to diverge as (T − TBKT)−1.5 [34]. In fact, the fitting of the γ-value by the formula (T − TBKT)−1.5 in Fig. 2e well reproduces the divergent behavior toward TBKT (see Supplementary Fig. 7 and Supplementary Note 8 for the reproducibility in another sample). The good consistency between the theory and the experiment supports that the enhancement of nonreciprocal transport can be understood in terms of the current-induced modulation of supercurrent density in the superconducting-proximity-coupled surface state.

High field behavior and angular dependence of nonreciprocal transport

In the above, we have considered the second harmonic resistance at a low magnetic field region, we next discuss R2 at higher fields. Figure 3a, b show the magnetic field dependence of R and R2 for T = 9.5 and 10 K. The breaking of superconductivity with the application of the magnetic field is clearly discerned in Fig. 3a. At low magnetic fields, the linear-positive slope of R2 with respect to the magnetic field is observed. The positive slope can be explained in terms of the modulation of supercurrent density as discussed earlier. At high magnetic fields, however, R2 shows a sign reversal. With the further application of the magnetic field, R2 goes to zero along with the breaking of superconductivity, which is clearly discerned at T = 10 K at around B = 14 T. The B–T diagrams for R and R2 are displayed as contour plots in Fig. 3c, d, respectively. We can again see that R2 appears positively at low fields and negatively at high fields at the intermediate temperature region. To discuss the origin of the negative component, we investigate the magnetic field direction dependence of R2. Figure 3e, f show the xy plane (in-plane) and zy plane (out-of-plane) magnetic-field directional dependence of R2, respectively. R2 becomes almost zero at B||x (φ = 0°, 180°, 360°) and B||z (θ = 0°, 180°, 360°), consistent with the symmetry of a superconducting-proximity-coupled TI[20-22] (see Supplementary Fig. 9 and Supplementary Discussion 1 for the discussion on the disappearance of R2 under B||z). In xy plane (Fig. 3e), R2 shows sin φ dependence both at 2 and 9 T, meaning that the signal is scaled with the y component of the magnetic field, consistent with Eq. (1). Similarly, in the zy plane (Fig. 3f), R2 shows sinθ behavior at 2 T, and the signal at 9 T basically follows sinθ behavior with the same sign as the 2 T signal. Remarkably, however, a sudden sign reversal of the signal is observed at θ = 90° and 270°, namely the negative component appears only when B is almost exactly aligned to in-plane (see Supplementary Figs. 6 and 8, and Supplementary Notes 7 and 8 for further analysis and reproducibility). In fact, the negative component almost vanishes at θ = 70°, leaving only the positive component, as shown in Fig. 3g. These results imply that the positive and the negative components are governed by the different origins. The sudden appearance of the negative component under the in-plane magnetic field indicates that it originates from the in-plane magnetic vortex along the y-direction[35]. The current-induced renormalization of supercurrent density may result in the modulation of the nucleation/annihilation energy barrier of the in-plane vortex, which can appear as nonreciprocal resistance[35]. Further theoretical and experimental investigation is required to fully reveal the origins of the negative component.

Fig. 3

Nonreciprocal transport at high magnetic fields. a The magnetic field dependence of first harmonic resistance Rω at T = 9.5 K (light blue) and T = 10 K (gray) measured under I = 200 μA. b The magnetic field dependence of second harmonic resistance R2. at T = 9.5K and T = 10 K measured under I = 200 μA. The red and blue triangles represent the positive and negative peaks at around B = 2 T and B = 9 T, respectively. c The contour plot of R in the plane of magnetic field and temperature at θ = 90° (B || y, defined in f) measured under I = 200 μA. The superconducting onset and R = 1 Ω are shown in circle and triangle, respectively. d The contour plot of R2 in the plane of magnetic field and temperature at θ = 90° measured under I = 200 μA. e The in-plane magnetic-field direction dependence of R2 at B = 2 T (red) and B = 9 T (blue) within xy plane measured under I = 200 μA. φ is defined as an angle in xy plane measured from the x-axis. f The out-of-plane magnetic-field direction dependence of R2 at B = 2 T (red) and B = 9 T (blue) within zy plane measured under I = 200 μA. θ is defined as an angle in zy plane measured from the z-axis. g The contour plot of R2 in the plane of magnetic field and temperature at θ = 70° measured under I = 200 μA

Discussion

The large nonreciprocal charge transport in the superconducting interface of Bi2Te3/FeTe can be utilized as a magnetically controllable superconducting diode[24] and as a rectenna at low temperatures[35]. The nonreciprocal transport can be well-explained by the modulation of supercurrent density due to spin–momentum locking, exemplifying the close connection between superconductivity and topological surface state. Thus, tuning the Fermi level to the surface state[36] will make the present system a desirable platform to study TSC and the associated formation of Majorana fermion. The nonreciprocal transport measurement is also applicable to other TSC candidates, such as LaAlO3/SrTiO3 interface[37,38], superconducting-proximity coupled Rashba wire[7,8], and topological surfaces[10-14], which will help to discuss the effect of spin splitting on the superconducting properties.

During the final revision of the manuscript, we noticed the independent ARPES study of Bi2Te3/FeTe, which suggests that charge transfer from FeTe to Bi2Te3 induces the interfacial superconductivity in hole-doped FeTe[39]. The interfacial superconductivity induces the proximity effect on the surface state of Bi2Te3, which is in line with the discussion of our present study.

Methods

MBE thin film growth

Thin films were grown by MBE on insulating CdTe(100) and InP(111) substrates. The growth temperatures for FeTe and Bi2Te3 were 400 and 180 °C, respectively. CdTe(100) substrates were etched with bromine–methanol (bromine 0.01%) for 5 min before the deposition. For FeTe, the flux ratio was fixed at Fe:Te = 1:30 and the growth rate was about 0.7 nm/min. For Bi2Te3, the flux ratio was fixed at Bi:Te = 1:2 and the growth rate was about 0.2 nm/min. To suppress the degradation of the film, AlO capping layer with a thickness of about 3 nm was deposited at room temperature with an atomic layer deposition system immediately after taking out the samples from the MBE vacuum chamber.

Device fabrication

The Hall-bar device pattern was defined by a photolithography technique and wet etching processes with HCl:H3PO4:H2O = 1:1:8. For Ohmic-contact electrodes, 5 nm Ti/45 nm Au were deposited with an electron-beam evaporator.

Transport measurements

The first and second harmonic resistances, R and R2, were measured using a current source (Keithley: Model 6221) and lock-in amplifiers (SRS: SR830). The measurement current and frequency were fixed at 200 μA (root mean square) unless otherwise noted and at 13 Hz, respectively. The measurements were done in Physical Property Measurement System (Quantum Design: PPMS). R2 was anti-symmetrized with respect to B.