Introduction

Superconducting nanowires represent an important example of a system where low-temperature physics is dominated by both thermal and quantum fluctuations1,2,3,4,5, thus making their properties entirely different from those of bulk superconductors well described by the standard Bardeen–Cooper–Schriffer (BCS) mean-field theory6.

A large part of fluctuation phenomena in such nanowires are attributed to the so-called phase slips1,2 which correspond to temporal local suppression of the superconducting order parameter \({{\Delta }}\exp (i\varphi )\) accompanied by the phase slippage process. At temperatures T close enough to the BCS critical temperature Tc such phase slips are induced by thermal fluctuations7,8,9 whereas at lower temperatures TTc quantum fluctuations of the order parameter take over and generate quantum phase slips (QPS)10,11.

As the phase φ changes in time by 2π during a QPS event, each such event causes a voltage pulse \(V=\dot{\varphi }/2e\) inside the wire. As a result, a current biased superconducting nanowire acquires a non-vanishing electric resistance down to lowest T10,11. This effect received its convincing experimental confirmation12,13,14,15,16,17,18. The same effect is also responsible for voltage fluctuations in superconducting nanowires19,20. Quantum phase slips also cause suppression of persistent currents in uniform superconducting nanorings21,22.

A fundamentally important property of superconducting nanowires is the so-called flux-charge duality. This feature was extensively discussed for ultrasmall Josephson junctions23,24,25,26 implying that under the duality transformation 2e ↔ Φ0 quantum dynamics of Cooper pairs (with charge 2e) should be identical to that of magnetic flux quanta Φ0 = hc/2e. All the same arguments remain applicable for shorter superconducting nanowires27 which properties are dual to those of small Josephson junctions (Fig. 1). The duality considerations can further be extended to longer nanowires5,22.

Fig. 1: Flux-charge duality.
figure 1

Quantum tunneling of a magnetic fluxon Φ0 across a superconducting nanowire and a dual tunneling process of a Cooper pair with charge 2e across a Josephson junction between two superconductors.

Manifestations of flux-charge duality in superconducting nanowires were observed in a variety of experiments thereby opening new horizons for applications of such structures in modern nanoelectronics, information technology, and metrology. These observations include, e.g., coherent tunneling of magnetic flux quanta through superconducting nanowires28,29 and the so-called Bloch steps30. Operations of duality-based single-charge transistor31,32 and charge quantum interference device33 were demonstrated. Superconducting nanowires were also proposed to serve as central elements for QPS flux qubits34 as well as for creating a QPS-based standard of electric current35.

Quantum fluctuations in superconducting nanowires are controlled by two different parameters

$${g}_{\xi }={R}_{q}/{R}_{\xi }\ \ \ \ {\rm{and}}\ \ \ {g}_{Z}={R}_{q}/Z.$$
(1)

Here Rq = h/e2 25.8 kΩ is the quantum resistance unit, Rξ is the normal state resistance of the wire segment of length equal to the superconducting coherence length ξ and \(Z=\sqrt{{\mathcal{L}}/C}\) is the wire impedance determined by the kinetic wire inductance (times length) \({\mathcal{L}}\) and the geometric wire capacitance (per length) C.

The dimensionless conductance gξ accounts for the fluctuation correction to the BCS order parameter9 Δ → Δ − δΔ (with δΔ ~ Δ/gξ) and determines the QPS amplitude (per unit wire length)11\({\gamma }_{{\rm{qps}}}=b({g}_{\xi }{{\Delta }}/\xi )\exp (-a{g}_{\xi })\) (with a ~ 1 and b ~ 1). The dimensionless admittance gZ, in turn, accounts for hydrodynamic (long wavelength) fluctuations of the superconducting phase intimately related to sound-like plasma modes36 propagating along the wire with the velocity \(v=1/\sqrt{{\mathcal{L}}C}\). Different quantum phase slips interact by exchanging such plasmons and, hence, the parameter gZ also controls the strength of inter-QPS interactions. By reducing the wire diameter \(\sqrt{s}\propto {g}_{Z}\) one eventually arrives at the “superconductor-insulator” quantum phase transition10 that occurs at gZ = 16 and T → 0.

In this work, we experimentally and theoretically investigate both global and local ground-state properties of superconducting nanowires in the “insulating” regime gZ < 16. We demonstrate that quantum fluctuations of magnetic flux in long nanowires yield effective localization of Cooper pairs at a fundamental length scale Lc that essentially depends on both parameters (1). We also show that nominally uniform nanowires exhibit a non-trivial mixture of superconducting-like features at shorter length scales and resistive long-scale behavior which should actually tend to insulating at T → 0. This state of matter can thus be named as a superconducting insulator.

Results

In order to accomplish our goal, we fabricated long and thin titanium nanowires having the form of narrow strips overlapping a relatively wide aluminum electrode through a tunnel barrier (aluminum oxide), as it is shown in Fig. 2. The normal state resistance of these wires RN measured above the BCS critical temperature Tc ≈ 400 mK was found in the range RN ~ 25–70 kΩ. The length L 20 μm and thickness d 35 nm remain the same for all Ti samples, whereas their width w varies in the range 30–60 nm within which quantum phase slips usually proliferate in Ti nanowires15,37. The zero-temperature superconducting coherence length in our Ti samples is estimated to be ξ ~ 140–150 nm and, hence, the quasi-one-dimensional limit condition d, wξL holds for all samples. With these parameters, one obtains the dimensionless admittance gZ ≈ 1–3, i.e., the desired condition gZ < 16 is well satisfied in all our nanowires. The dimensions of the aluminum strip are large enough, enabling one to ignore fluctuation effects.

Fig. 2: Schematics of experiment and sample layout.
figure 2

a Long and thin titanium nanowires having the form of narrow strips overlap a relatively wide aluminum electrode through a tunnel barrier (aluminum oxide). The structure enables one to carry out both pseudo-four-terminal measurements of the total resistance for all nanowires and local measurements of the current–voltage characteristics for all Al–AlOx–Ti tunnel junctions. b Scanning electron microscope image of our typical structure. Inset: Zoom of the junction region taken with atomic force microscope. Fake color corresponds to variation of the sample height from 0 (substrate, dark blue) to 80 nm (overlapping titanium, orange).

Nanowire resistance

The results of our measurements of a total resistance R(T) for five different nanowires are displayed in Fig. 3a. With the values gZ 16, in the low-temperature limit all these samples should remain deep in the insulating regime. We observe, however, that two thicker samples with nominal widths w ≈ 62 nm (sample Ti1) and w ≈ 46 nm (sample Ti2), demonstrate a pronounced resistive behavior with R(T) ≈ RN only at temperatures not far below the bulk titanium critical temperature Tc ≈ 400 mK followed by a rather sharp resistance drop by ~2 orders of magnitude at temperatures T ~ 300 mK (sample Ti1) and T 200 mK (sample Ti2). The remaining samples Ti3, Ti4, and Ti5 with nominal widths just slightly below that for Ti2 (respectively, w ≈ 41, 40, and 30 nm) show no sign of superconductivity down to the lowest T and only very weak dependence R(T), in particular for the thinnest samples Ti4 and Ti5.

Fig. 3: Resistance data and localization of Cooper pairs.
figure 3

a Resistance versus temperature R(T) measured for five Ti nanowires of length L = 20 μm, thickness d 35 nm, and nominal width values w indicated in brackets for each of the samples. Solid lines represent fits of the data for samples Ti1 and Ti2 to the theory9 within its validity range. Resistance saturation observed in these two samples at low T is due to finite voltage sensitivity of about few nV corresponding to residual resistance ~100 Ω measured using ac bias current ~10 pA rms. Error bars are smaller than data points. b Localization of Cooper pairs (with charge 2e) generated by quantum tunneling of magnetic fluxons Φ0 across the nanowire. This phenomenon explains the low-temperature behavior of R(T) observed in samples Ti3, Ti4, and Ti5.

At temperatures not far below Tc the system behavior should be dominated by thermally activated phase slips which contribution to the wire resistance Rtaps(T)9 indeed provides very accurate fits for the resistance of two of the above samples (see Fig. 3a) and allows to extract effective values gξ 37.4 and gξ 9.0, respectively, for samples Ti1 and Ti2 (see Supplementary Note 1 for more details). These values are smaller than the nominal ones, most likely indicating certain non-uniformity of our nanowires.

Localization of Cooper pairs

In order to understand drastic difference in the low-temperature behavior of our samples with various cross-sections it is necessary to account for the effect of quantum phase slips. The dual Hamiltonian for superconducting nanowires in the presence of QPS reads5,22

$$\hat{H}=\int_{0}^{L}dx\left[\frac{{\hat{{{\Phi }}}}^{2}}{2{\mathcal{L}}}+\frac{{({\partial }_{x}\hat{Q})}^{2}}{2C}-{\gamma }_{{\rm{qps}}}\cos \left(\frac{\pi \hat{Q}}{e}\right)\right],$$
(2)

where \(\hat{{{\Phi }}}\) and \(\hat{Q}\) are canonically conjugate flux and charge operators obeying the commutation relation \([\hat{{{\Phi }}}(x),\hat{Q}(x^{\prime} )]=-i\hslash \delta (x-x^{\prime} )\). Employing this Hamiltonian one can demonstrate5,22 that in the “insulating” phase, i.e., for gZ < 16, the wire ground-state properties are controlled by a non-perturbative correlation length \({L}_{c}\propto {\gamma }_{{\rm{qps}}}^{-\alpha }\) with 1/α = 2 − gZ/8 or, equivalently,

$${L}_{c} \sim \xi \exp \left(\frac{a{g}_{\xi }-{\mathrm{ln}}\,b}{2-{g}_{Z}/8}\right).$$
(3)

Physically the appearance of this QPS-induced fundamental length scale can be viewed as a result of spontaneous tunneling of magnetic flux quanta Φ0 back and forth across the wire, as it is illustrated in Fig. 3b. These quantum fluctuations of magnetic flux wipe out phase coherence at distances ~Lc and yield effective localization of Cooper pairs at such length scales. Accordingly, samples with LLc may still exhibit superconducting properties also in the presence of QPS, whereas in the limit LLc the supercurrent gets disrupted by quantum fluctuations and such nanowires remain non-superconducting even at T → 0.

This is exactly what the data in Fig. 3a demonstrate. Indeed, the value Lc (3) for the sample Ti1 with gξ ≈ 37 obviously exceeds L by several orders of magnitude, and hence, this sample should remain superconducting at low enough T. In order to estimate the length scale (3) for sample Ti2 with ξ ~ 140 nm, gξ 9.0, and gZ 2.5 it is desirable to explicitly determine the prefactors a and b. The data analysis for this sample yields a lower bound for the combination \(a{g}_{\xi }-{\mathrm{ln}}\,b\;\gtrsim\; 7.5\), see Supplementary Note 2 for details. With this in mind, Eq. (3) allows to estimate Lc 12 μm, i.e., in this case, Lc ~ L and the sample Ti2 should also remain superconducting at low T in accordance with our observations. By contrast, three thinner nanowires Ti3, Ti4, and Ti5 with lower effective values gξ and Lc significantly smaller than L exhibit a non-superconducting behavior down to lowest T.

In order to interpret this behavior let us recall that for gZ < 16 quantum phase slips are no longer bound in pairs. According to the exact solution for the sine-Gordon model38, in this case, an effective minigap in the spectrum \(\tilde{{{\Delta }}}\propto {\gamma }_{{\rm{qps}}}^{\alpha }\) develops implying that at T → 0 samples Ti3, Ti4, and Ti5 should behave as insulators. In line with these arguments, our resistance data in Fig. 3a demonstrate that the supercurrent in these samples is fully blocked by QPS down to the lowest available temperatures, and hence, their insulating behavior should indeed be expected at \(T\;<\;\tilde{{{\Delta }}}\). The absence of any visible resistance upturn at low T most likely implies that the latter condition is not yet reached and/or the inequality LLc is not satisfied well enough for these samples. In any event, here superconductivity is totally wiped out by quantum fluctuations in accordance with our theoretical arguments.

Note that the resistance data similar to those of Fig. 3a were also reported previously12,13,39 for a large number of MoGe nanowires with shorter values of ξ and L. In some of these samples, the resistance upturn at lower T indicating the insulating behavior was observed. Reanalyzing the data12,13,39 we conclude that they are also consistent with the above physical picture involving the correlation length Lc (3), i.e., the superconducting MoGe samples obey the condition LLc, whereas the non-superconducting ones typically have the length L exceeding Lc. Hence, retrospectively the observations12,13,39 also receive a natural explanation which was not yet available at that time.

Local properties

Measurements of the total resistance R(T) alone are not yet sufficient to obtain complete information about the quantum mechanical ground state of superconducting nanowires. In order to probe their local properties, we performed measurements of the IV curves for tunnel junctions between Ti nanowires and bulk Al electrodes (with the BCS gap ΔAl 190 μV), see Fig. 2. The corresponding results for all five samples are displayed in Fig. 4. In these samples, the differential conductance for Ti–Al tunnel junctions has a peak which position varies slightly from sample to sample. As the peak is expected to occur at eV = Δ + ΔAl, we immediately reconstruct the local gap value ranging between Δ ≈ 50 μeV and Δ ≈ 37 μeV depending on the sample. Hence, quantum fluctuations tend to reduce Δ in superconducting nanowires below its bulk value ΔTi 60 μeV and this effect appears more pronounced for thinner samples. On the other hand, a non-zero local superconducting gap Δ remains clearly observable in all our samples.

Fig. 4: Local differential conductance and electron density of states.
figure 4

a Differential conductance dI/dV as a function of voltage V measured in Ti−Al tunnel junctions at T 21 mK for five samples Ti1 to Ti5. A sharp peak is observed at eV = Δ + ΔAl. Inset: The same data for sample Ti3 at different temperatures. b Fit of the data for sample Ti3 at T 21 mK to the theory40. Inset: The density of states ν (in units of the normal density of states at the Fermi energy) as a function of energy E reconstructed for the same sample at the same T. Error bars are smaller than data points.

As compared to the standard BCS-like IV curve, systematic broadening of this peak in dI/dV with decreasing wire cross-section is observed. This broadening increases with T (cf. inset in Fig. 4a) and it can be explained40,41 if we bear in mind that electrons exchange energies with an effective dissipative environment formed by Mooij–Schön plasmons propagating along the wire. As a result, in our Ti nanowires the singularity in the electron density of states (DOS) ν(E) at E = Δ and T → 0 gets weaker with decreasing wire cross-section and becomes washed out by quantum fluctuations at gZ ≤ 2.

This is exactly what we observe in our experiment. By fitting the corresponding IV data for Ti−Al tunnel junctions to theoretical predictions40 (see Supplementary Note 3) we reconstruct the energy-dependent DOS ν(E) for our Ti nanowires, as displayed in Fig. 4b. The best fit for sample Ti3 yields the value gZ 1.50 just slightly below our theoretical estimate gZ 2.26. In contrast to the standard BCS dependence \({\nu }_{{\rm{BCS}}}(E)={\rm{Re}}| E| /\sqrt{{E}^{2}-{{{\Delta }}}^{2}}\), here the gap singularities are totally smeared due to electron-plasmon interactions. Nevertheless the superconducting gap in DOS ν(E) remains clearly visible. At non-zero T and subgap energies DOS decays exponentially with decreasing E as \(\nu (E)\propto \exp (-({{\Delta }}-| E| )/T)\) (cf. inset in Fig. 4b) which is also due to the interaction between electrons and Mooij–Schön plasmons40.

A pronounced superconducting gap in DOS is not the only feature indicating that locally superconducting properties remain preserved despite the effect of quantum fluctuations. In Fig. 5, we display the IV curves measured at different T for yet one more Ti nanowire with local DOS also showing a pronounced superconducting gap (Fig. 5a) and total resistance R(T) behaving qualitatively similar to that of the samples Ti3, Ti4, and Ti5 (cf. inset in Fig. 5b). Zooming at the origin of these IV curves we clearly observe the Josephson current ~5 pA at T 75 mK, see Fig. 5b. Quite naturally, due to strong fluctuation effects inside our Ti nanowire42,43 this current value is orders of magnitude smaller than the nominal maximum Josephson current of few nA estimated from the standard Ambegaokar–Baratoff formula6. Fluctuation effects become even stronger with increasing temperature and totally wash out the Josephson current already at T 150 mK.

Fig. 5: Current–voltage characteristics and Josephson current.
figure 5

a Current–voltage characteristics for Ti–Al tunnel nanojunction corresponding to the Ti nanowire with L = 20 μm, d = 35 nm, and w = 38 nm recorded at T = 26 mK. b Zoom of the current versus voltage dependencies taken at various temperatures. At T 75 mK one clearly observes the Josephson current which gradually disappears at higher T. Inset: The total resistance R measured for this nanowire as a function of temperature. Error bars are smaller than data points.

Note that the same Josephson current feature is detected in other Ti nanowires at low T and V → 0, cf., e.g., Fig. 4a (inset) and Fig. 4b. These observations of dc Josephson effect in Ti–Al tunnel junctions further support our conclusion suggesting the presence of local superconductivity in all investigated Ti samples, including the most resistive ones.

Discussion

We arrive at the following physical picture describing ultrathin superconducting wires in the “insulating” regime gZ < 16 at low enough temperatures. In this regime, QPS proliferate while TAPS effects can already be neglected. In thicker nanowires with LcL (samples Ti1 and Ti2) quantum phase slips alone cannot disrupt phase coherence across the wire. Such samples then behave to a large extent similarly to effectively zero-dimensional objects, such as, e.g., small-size Josephson junctions with the fluctuating phase26 embedded in a low resistive external circuit. Depending on the experimental realization22,42,43, these nanowires may either stay superconducting or become resistive, albeit typically with rather small \(R\propto {\gamma }_{{\rm{qps}}}^{2}\). In contrast, thinner samples with LcL remain highly resistive with R ~ RN even at TTc and should turn insulating in the limit T → 0. This behavior is due to QPS that suppresses long-range phase coherence in such nanowires.

Remarkably, the superconducting gap Δ in the energy spectrum of all our Ti nanowires, including highly resistive ones, is reduced but not destroyed by quantum fluctuations. In addition, this spectrum is also affected by the interaction between electrons and soft phase fluctuation modes (Mooij–Schön plasmons) which washes out the BCS gap singularity in DOS of ultrathin (gZ < 2) nanowires and produces a weak subgap tail in ν(E) at non-zero T. We have demonstrated that the wire segments of length Lc retain their superconducting properties. On the other hand, longer nanowires composed of many such superconducting segments exhibit effective localization of Cooper pairs at lengths ~Lc and loose their ability to sustain any measurable supercurrent. These nanowires demonstrate a resistive behavior with R(T) ~ RN even at TTc and should turn insulating in the limit of large L and T → 0.

It is well-known that under certain conditions granular superconducting arrays and Josephson junction chains may also become resistive and even insulating44,45,46,47. In that regime, superconductivity is well preserved only inside grains while dissipativeless charge transfer across the system is prohibited due to Coulomb blockade of Cooper pair tunneling between such grains. Here, in contrast, we are dealing with nominally uniform nanowires which do not contain any grains and dielectric barriers at all. Nevertheless, such nanowires may exhibit both resistive and insulating behavior as long as their length L strongly exceeds the typical size of a “superconducting domain” \({L}_{c}\propto {\gamma }_{{\rm{qps}}}^{-\alpha }\). Similarly to normal metallic structures1,48,49, this non-trivial feature can be interpreted as weak Coulomb blockade of Cooper pairs that—as it is illustrated by our results—may occur even in the absence of tunnel barriers.

In summary, we have demonstrated—both experimentally and theoretically—that long and uniform superconducting nanowires in the so-called “insulating” regime actually exhibit a more complicated behavior characterized by superposition of local superconductivity and effective global localization of Cooper pairs. This fundamental property of superconducting nanowires needs to be accounted for while designing various nanodevices with novel functionalities.

Methods

E-beam lift-off process, vacuum deposition of metals and in situ oxidation were used to fabricate tunnel junctions between aluminum electrodes and titanium nanostripes. Each structure enables one to carry out both pseudo-four-terminal measurements of the total resistance R(T) for all Ti nanowires and local measurements of the IV curve for all Al−Ti tunnel junctions (Fig. 2). Differential conductances dI/dV were obtained by modulation technique using lock-in amplification. All experiments were made inside 3He4He dilution refrigerator with carefully filtered50 input/output lines connecting sample to laboratory digital electronics through battery-powered analog pre-amplifiers (see Supplementary Note 4 for details).