High-quality III–V narrow bandgap semiconductor materials with strong spin–orbit coupling and large Landé g-factor provide a promising platform for next-generation applications in the field of high-speed electronics, spintronics, and quantum computing. Indium antimonide (InSb) offers a narrow bandgap, high carrier mobility, and small effective mass and, thus, is very appealing in this context. In fact, this material has attracted tremendous attention in recent years for the implementation of topological superconducting states supporting Majorana zero modes. However, high-quality heteroepitaxial two-dimensional (2D) InSb layers are very difficult to realize owing to the large lattice mismatch with all commonly available semiconductor substrates. An alternative pathway is the growth of free-standing single-crystalline 2D InSb nanostructures, the so-called nanoflags. Here, we demonstrate fabrication of ballistic Josephson-junction devices based on InSb nanoflags with Ti/Nb contacts that show a gate-tunable proximity-induced supercurrent up to 50 nA at 250 mK and a sizable excess current. The devices show clear signatures of subharmonic gap structures, indicating phase-coherent transport in the junction and a high transparency of the interfaces. This places InSb nanoflags in the spotlight as a versatile and convenient 2D platform for advanced quantum technologies.
Today, great interest revolves around the possibility to create and manipulate new states of matter with topological properties. This stems mostly from the intrinsic robustness of topological states against local perturbation and the ensuing relevance for quantum computing architectures.1,2 Hybrid superconductor/semiconductor heterostructures represent a promising platform in which topological properties can emerge.3–6
In this context, indium antimonide (InSb) has attracted much attention recently. InSb has a narrow bandgap ( eV).6–8 It also has a very high bulk electron mobility [ cm2/(V s)]9,10 and a small effective mass (),8,9,11–14 which are both important requirements for high-speed and low-power electronic devices.10,15 Finally, it also exhibits a strong spin–orbit interaction and a large Landé g-factor (9,14), and thus, it is useful for spintronics applications8,15 and for the creation of hybrid structures hosting topological states like Majorana zero-modes. Indeed, the first signatures compatible with Majorana bound states were reported in InSb nanowires coupled to a superconductor,3,16 which has triggered strong efforts to improve the quality of hybrid systems based on InSb nanowires.17–24
Aside from one-dimensional nanowires, two-dimensional (2D) InSb structures also attract great attention, owing to their inherent design flexibility.6,8,14 Indeed, InSb 2D electron gases12 and the related ternary compound InSbAs6 have recently been proposed as a platform for topological superconductivity,25 and ballistic superconductivity was demonstrated in InSb quantum wells.12 However, the growth of high-quality heteroepitaxial 2D InSb layers is still a challenge owing to their large lattice mismatch with common semiconductor substrates. In addition, such quantum wells are reported to suffer from instabilities due to the Si dopants.14 A possible strategy to circumvent these problems consists in growing free-standing 2D InSb nanostructures on nanowire stems, because nanowires yield efficient relaxation of elastic strain along the nanowire sidewalls, when lattice-mismatched semiconductor systems are integrated. To emphasize their free-standing 2D shape, such nanostructures are often referred to as nanosails, nanosheets, nanoflakes, or nanoflags. However, until today, only a few studies were reported on the growth and the electrical transport properties of such InSb nanoflags (NFs).7,9,13,19,26–30
InSb NFs were first reported in 2016 by de la Mata et al.9 their growth being based on molecular beam epitaxy (MBE). There, the authors attributed the 2D geometry of the NFs to a single twinning event in the otherwise pure zinc blende structure of the InSb sample, and four-terminal electrical measurements revealed an electron mobility greater than 12 000 cm2/(V s). Pan et al. used Ag-assisted MBE to grow free-standing 2D single-crystalline InSb NFs.26 Hall-bar devices were then fabricated that showed ambipolar behavior and an electron mobility of 18 000 cm2/(V s). The same group also demonstrated functional InSb NF Josephson junctions (JJs) with Al and Nb contacts.27 Furthermore, in the devices with Nb contacts, the authors reported evidence of the coexistence of the quantum Hall effect and proximity-induced superconductivity in the InSb NFs.28 Quantum transport and quantum-dot geometries in such nanostructures were also demonstrated very recently.7,31–33 Gazibegovich et al. combined selective-area growth with the vapor-liquid-solid mechanism in metal organic vapor phase epitaxy, leading to the formation of InSb NFs thanks to the development of a twin-plane boundary.19,34 The same group reported evidence for crossed Andreev reflections in JJs made from such flakes.13
Also our group demonstrated the growth of single-crystal, free-standing InSb NFs, initially on InAs nanowires, using a directional growth technique based on chemical beam epitaxy (CBE).29 The NF size was limited by the flexibility of the InAs nanowires, which led to a bending of the stem and the resulting loss of the orientation for the asymmetric 2D growth. In order to overcome this issue, we recently optimized the growth of InSb NFs.30 In particular, the InSb NFs were grown on sturdy tapered InP nanowires, which did not bend and allowed to grow larger NFs with the same directional-growth approach. This strategy allowed us to obtain InSb NFs of μm length, nm width, and nm thickness.30
The resulting NFs are large enough to fabricate Hall bars with length-to-width ratios enabling precise electrical characterization. An electron mobility of 29 500 cm2/(V s) was measured at a carrier concentration cm−2 at 4.2 K.30 The electron elastic mean free path reached values of 500 nm, which favorably compares with the literature.9,13,34
Here, we report on the fabrication and characterization of JJ devices based on these InSb NFs and provide evidence of ballistic superconductivity. We employ Ti/Nb contacts in InSb JJ devices and show both gate-tunable proximity-induced supercurrent at 250 mK and sizable excess current. The devices also show clear signatures of subharmonic gap structures, indicating phase-coherent transport in the junction and highly transparent interfaces. Our results indicate InSb NF as a promising platform for the study of topological superconductivity.
The upper left inset of Fig. 1 shows a SEM image of the device investigated in this work (for consistent results from other devices, see the supplementary material). In brief, a 100 nm-thick InSb NF was transferred mechanically on a SiO2/Si substrate and contacted with 10/150 nm Ti/Nb. (More details on device fabrication are provided in the supplementary material.) The interelectrode spacing between the two superconductors, i.e., the length of the normal (N) region, is L = 200 nm, while its width is W = 700 nm. Standard transport characterization yields a mean-free path of nm30 for the N region, greater than the junction length L. These numbers place the device in the ballistic regime.
DC voltage drop Vsd as a function of the bias current Isd. A supercurrent of nA is observed. The lower right inset shows the differential resistance dV/dI measured simultaneously by the lock-in technique. Vbg = 30 V, T = 250 mK. B = 6 mT applied to compensate for the residual magnetization of the cryostat. The upper left inset shows an SEM image of the device structure. Scale bar 2 μm.
DC voltage drop Vsd as a function of the bias current Isd. A supercurrent of nA is observed. The lower right inset shows the differential resistance dV/dI measured simultaneously by the lock-in technique. Vbg = 30 V, T = 250 mK. B = 6 mT applied to compensate for the residual magnetization of the cryostat. The upper left inset shows an SEM image of the device structure. Scale bar 2 μm.
The critical temperature of the superconducting leads was determined to be K, from which the bulk gap can be computed using the theory by Bardeen, Cooper, and Schrieffer (BCS):35 meV, consistent with values of Nb superconducting contacts previously reported in the literature.28,36–41 The induced superconducting coherence length13,27,42,43 is , with vF being the Fermi velocity in the N region ( m/s) and Δ being the gap in the superconductor. Here, nm , so the device is in the short junction regime. Equivalently, the Thouless energy42,44 meV .
Figure 1 shows typical voltage–current (V – I) characteristics obtained at T = 250 mK and Vbg = 30 V. The device displays well-developed dissipationless transport, thus demonstrating proximity-induced superconductivity in the InSb NF. As the bias current exceeds the critical value of nA, a sudden jump in the measured voltage to a dissipative quasiparticle branch is observed, indicating that the JJ switches from the superconducting to the normal state with a resistance of Ω. Current sweeps in opposite directions show negligible hysteresis, i.e., switching and retrapping current are the same, so that in the following, we shall use the switching current and critical current as synonyms. Consistently, the switching current is larger than the intrinsic thermal current noise of the junction27,45 ; here, nA. The lower right inset to Fig. 1 shows the differential resistance dV/dI measured using a lock-in amplifier together with the V – I curve. Data clearly show that the differential resistance is zero in the supercurrent branch of the device. Zhi et al. report a supercurrent of 20 nA at 10 mK in Nb/InSb nanoflag SNS junctions.27 We attribute the improved numbers reported here mainly to a higher mobility of the nanoflags and progress in device fabrication.
Superconducting quantum interference was observed in the dependence of the supercurrent on a magnetic field applied perpendicularly to the sample plane (Fig. 2). Supercurrent maximum is obtained for mT instead of the expected maximum at zero B field. This small offset can be attributed to a residual magnetization in the cryostat. Applying higher or lower magnetic fields, the suppercurrent symmetrically decreases, until for mT it is suppressed. The modulation of the critical current by quantum interference is one of the hallmarks of the Josephson effect. The shape of the curve resembles a Fraunhofer pattern with only the central lobe, i.e., without side-lobes. We have verified that in a magnetic field range of mT, no side-lobes appear. de Vries et al. have studied similar InSb NF JJs and report an even-odd Fraunhofer pattern.13 Thus, the intensity of the first side-lobes might be anomalously small, which precludes their observation in our experiment. On the other hand, such anomalous magnetic interference patterns, with a monotonous decay, were reported previously in similar geometries36,37,46,47 and were attributed to geometric factors.36,46,48–53 In addition, the magnetic flux through the junction is , with A the junction area, .54 Here, λL is the London penetration depth of Nb. Thus, the smaller the channel width W, the higher the value of B required to reach . According to Rohlfing et al.,36 larger values of B more strongly suppress via a dramatic reduction of the amplitude of Andreev reflections. Future measurements on devices with larger W/L ratio might help us to clarify whether the missing lobes in Fig. 2 are due to geometric factors.
2D colormap of differential resistance dV/dI, obtained by numerical derivation of measured DC V – I curves, plotted vs current bias Isd and magnetic field B, measured at Vbg = 40 V and T = 250 mK.
2D colormap of differential resistance dV/dI, obtained by numerical derivation of measured DC V – I curves, plotted vs current bias Isd and magnetic field B, measured at Vbg = 40 V and T = 250 mK.
InSb NFs are n-type semiconductors, and the carrier concentration in the NF can be tuned by an applied gate voltage.30 Figure 3(a) shows that the global gate modulation can also be employed to control the magnitude of the supercurrent maximum. The figure shows the differential resistance dV/dI of the device as a function of the current bias Isd and the back gate voltage Vbg. The central white region represents the zero-resistance supercurrent branch. In a range of gate voltages from 20 to 50 V, the supercurrent is approximately constant at nA. Decreasing the gate voltage below 20 V, the supercurrent decreases, until it disappears for V. This demonstrates the gate voltage control of the supercurrent magnitude and that the device implements a Josephson field effect transistor (JoFET).55–59 The corresponding V – I curves as a function of the back gate voltage are shown in Fig. 3(b) as line plots. While the V – I curve for Vbg = 10 V still shows a zero slope at the origin, already for Vbg = 7 V the V – I curve is essentially linear with Ohmic behavior. The gate voltage dependence of the switching current is shown in Fig. 3(b) as red line and confirms that below Vbg = 20 V the switching current decreases. The normal resistance (the slope of the V – I curves in the normal branch) displays an opposite behavior: below Vbg = 20 V, it increases significantly from Ω to above 10 kΩ. The product of the switching current and normal resistance, , is approximately constant at 15 μV in a wide range from Vbg = 10 V to 50 V (see Fig. S1 in the supplementary material), while it drops to zero when the switching current becomes zero. Similar results were reported by Zhi et al.27
(a) Color-scale plot of differential resistance dV/dI, measured by the lock-in technique, as a function of the current bias Isd and back gate voltage Vbg. (b) 3D plot shows the trend of simultaneously measured DC V – I curves at different back gate voltages. The first curve (at lowest back gate voltage) is highlighted in yellow. Red lines indicate the transition between the superconductive and the dissipative regime. The supercurrent increases with increasing the back gate voltage above pinch off. T = 250 mK. B = 6 mT applied to compensate for the residual magnetization of the cryostat.
(a) Color-scale plot of differential resistance dV/dI, measured by the lock-in technique, as a function of the current bias Isd and back gate voltage Vbg. (b) 3D plot shows the trend of simultaneously measured DC V – I curves at different back gate voltages. The first curve (at lowest back gate voltage) is highlighted in yellow. Red lines indicate the transition between the superconductive and the dissipative regime. The supercurrent increases with increasing the back gate voltage above pinch off. T = 250 mK. B = 6 mT applied to compensate for the residual magnetization of the cryostat.
Next, we characterize the dissipative regime. Figure 4(a) shows subharmonic gap structures in the differential conductance that can be attributed to multiple Andreev reflections (MARs). The peak present at Vsd = 0 V corresponds to the superconductive state. On the other hand, above mV, the differential conductance becomes constant and is equal to the inverse of the normal resistance, . Between these two extrema, the differential conductance dI/dV displays characteristic singularities (minima and maxima), which represent the subharmonic gap structures.60–63 Their presence is a signature of the high transparency of the interfaces between S and N regions. The positions of these MAR singularities follow the equation with and being the induced gap in the N region. Most commonly, the position of the maxima in the differential conductance has been analyzed,6,12,13,18,27,28,36,41,64,65 but recently it was pointed out that for highly transparent junctions, the MAR resonances appear as minima in the differential conductance.66
(a) Differential conductance dI/dV, measured by the lock-in technique, as a function of the source-drain voltage Vsd for several back gate voltages Vbg, as indicated in the legend. The dashed lines indicate minima in differential conductance caused by multiple Andreev reflections ( to ). T = 250 mK. B = 6 mT applied to compensate for the residual magnetization of the cryostat. (b) Conductance line trace obtained from a coherent scattering model vs source-drain voltage Vsd. The theoretical fit of the data for Vbg = 15 V yields a transmission Tr = 0.94 and a value of the induced gap of μeV.
(a) Differential conductance dI/dV, measured by the lock-in technique, as a function of the source-drain voltage Vsd for several back gate voltages Vbg, as indicated in the legend. The dashed lines indicate minima in differential conductance caused by multiple Andreev reflections ( to ). T = 250 mK. B = 6 mT applied to compensate for the residual magnetization of the cryostat. (b) Conductance line trace obtained from a coherent scattering model vs source-drain voltage Vsd. The theoretical fit of the data for Vbg = 15 V yields a transmission Tr = 0.94 and a value of the induced gap of μeV.
In order to estimate the junction transparency and the induced gap, we used a simple scattering model that assumes fully coherent transport across a multimode JJ (see Ref. 67 and the supplementary material) and that has been applied to reproduce MAR traces of nanowire junctions.68 Thus, the experimental curves were compared to optimized theoretical MAR conductance traces. One example is shown in Fig. 4(b). The best agreement between the experiment and theory is obtained for a junction model with 40 modes of transparency Tr = 0.94 and an induced gap of μeV.69 The dashed vertical lines in Fig. 4 highlight the series from to , showing an excellent agreement with the predicted behavior for the minima in differential conductance in a highly transparent junction.66 Results of the analysis for the other back gate voltages are very similar. Consistently, while the transparency of the modes and the induced gap do not change with back gate voltage, the number of modes does: from less than 30 at Vbg = 10 V to more than 50 at Vbg = 30 V. The obtained value of the induced gap, μeV, is smaller than the value extracted from the measurement of the critical temperature. We note that similar values for the induced gap have been reported for InSb nanowires proximitized by Nb64 and by NbTiN.16
Aside from that differential conductance data, we recorded also DC traces (I – V curves). Figure 5 shows a representative example obtained with Vbg = 40 V. From a linear fit of the part of the curve at high bias, with , where Andreev reflections are completely suppressed, we obtain an excess current of nA and a normal resistance of Ω. We add that the results do not change significantly if we consider only the voltage range meV. These numbers result in a product μV, an important figure of merit for weak links.70 From Vbg = 40 to 20 V, the excess current is almost constant, while it decreases for smaller back gate values, to nA at Vbg = 5 V. The normal resistance displays the opposite behavior: below Vbg = 20 V, it increases significantly to about 3.1 kΩ at Vbg = 5 V. As shown in the supplementary material, the product of the excess current and normal resistance, , remains approximately constant at 137 ± 19 μV over the whole range of back gate voltages explored, i.e., from Vbg = 5 V to 40 V.
I – V characteristics measured in DC. The linear fit (dashed line) of the I – V curve for gives the excess current Ie as the intercept at zero bias voltage. The inset displays the excess current in a smaller bias voltage range. T = 250 mK and Vbg = 40 V. B = 6 mT applied to compensate for the residual magnetization of the cryostat.
I – V characteristics measured in DC. The linear fit (dashed line) of the I – V curve for gives the excess current Ie as the intercept at zero bias voltage. The inset displays the excess current in a smaller bias voltage range. T = 250 mK and Vbg = 40 V. B = 6 mT applied to compensate for the residual magnetization of the cryostat.
We should like to analyze the superconducting-gap values as obtained from the critical temperature of the superconductor and the observed multiple Andreev-reflection features. The measured critical temperature K is close to the reported value for bulk Nb, and the resulting value of the gap meV is in good agreement with values reported for JJs with Nb contacts.28,36–41,65 Thus, we attribute the observed critical temperature to a switching of the Nb film from the superconducting state to the normal state.
On the other hand, several groups reported gap values extracted from an analysis of MAR features that were smaller than the BCS-like gap of the superconducting leads.12,24,64 Kjaergaard et al. investigated multiple Andreev reflections in an InAs quantum-well heterostructure with epitaxial Al.66 From an analysis of the MAR features, authors obtained a gap value smaller than the BCS-like gap Δ of Al and showed that this is due to an induced gap in the quantum well covered by Al. Andreev reflections of particles in the uncovered region occur at the (vertical) effective interface to the covered region with gap in the quantum well. Since these reflections occur within the InSb crystal, the transparency of the process is high (here, Tr = 0.94).
Another relevant effect is the proximity-effect transfer from Nb into InSb via the thin Ti film. Drachmann et al. studied the proximity-effect transfer from a NbTi film into an InAs quantum well via a thin epitaxial Al layer.71 They found that the induced gap was increased by the NbTi film compared to samples with just the Al film, but it was still smaller than the BCS-like gap of NbTi. This implies that the proximity effect transfer can be weakened by an intermediate superconducting layer with smaller Tc. We recall that the reported Tc value for bulk Ti is K,72,73 i.e., much smaller than the Tc of Nb. The combination of both effects is likely able to explain the observed value of the induced gap in our experiments.
The transparency of the (horizontal) interfaces between the superconducting electrodes and InSb can be estimated using the theory of Aminov et al.,74 which measures the transparency via the dimensionless parameter γB ( for perfectly transparent interfaces, larger γB for larger barriers). Using the BCS gap of Nb, we obtain . Considering that the presence of the Ti film will slightly reduce the gap by an (unknown) amount, . This indicates a small transparency of the interface, consistent with the fact that the induced gap is much smaller than the BCS gap of Nb. On the other hand, Kjaergaard et al.66 and Baumgartner et al.75 report , consistent with their use of epitaxial Al/InAs heterostructures, which are known to have highly transparent interfaces.
Theory predicts for JJ at T = 0 that the product is constant proportional to the gap, , with the prefactor α being constant of order unity.45,76–80 Here, μV is only about 10% of μV. Such a reduction is frequently observed in experiment18,27,81 and has been attributed to a premature switching of the junction due to thermal activation.76,82 On the other hand, excess current is due to Andreev reflections and, thus, depends primarily on the transparency of the (vertical) interface between the covered and uncovered parts of the semiconductor,80 which is high. Consequently, a large product is observed, close to the theoretical value of 8/3 for ballistic junctions.63,80
In summary, we have fabricated JJ devices with InSb NFs as the normal region and Ti/Nb as superconducting contacts. The high electron mobility and large mean free path of the InSb NFs yielded ballistic transport across the normal region of the junction. We showed Josephson coupling between the superconductor and semiconductor, as demonstrated by the zero-resistance supercurrent of nA and the observation of MARs. Analysis of the MAR traces indicates a very high transparency of the interfaces. We also observe a sizable excess current. Our results show that free-standing 2D InSb NFs on InP stems, thanks to their defect-free zinc blende crystal structure,30 are a suitable material platform for fabrication of quantum devices. Considering also their strong spin–orbit interaction and their large Landé g-factor, we envision the use of these structures in future studies toward topological superconductivity.
See the supplementary material for extended methods and additional data.
We gratefully acknowledge helpful discussions with Hervé Courtois. We thank Daniele Ercolani for his help with the growth of the NF. This research activity was partially supported by the SUPERTOP project, QUANTERA ERA-NET Cofound in Quantum Technologies (H2020 Grant No. 731473) and by the FET-OPEN project AndQC (H2020 Grant No. 828948). M.P.N. acknowledges support by the National Science Centre (NCN), Agreement No. UMO-2020/38/E/ST3/00418.
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
I.V., V.Z., and L.S. grew the InSb NF. S.S. fabricated the devices. S.S. and S.H. carried out the experiments. M.P.N. carried out the numerical simulations. S.S., M.C., and S.H. analyzed the data and wrote the manuscript with strong input from all coauthors.
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request.