Dynamic properties of high-Tc superconducting nano-junctions made with a focused helium ion beam.

The Josephson junction (JJ) is the corner stone of superconducting electronics and quantum information processing. While the technology for fabricating low Tc JJ is mature and delivers quantum circuits able to reach the "quantum supremacy", the fabrication of reproducible and low-noise high-Tc JJ is still a challenge to be taken up. Here we report on noise properties at RF frequencies of recently introduced high-Tc Josephson nano-junctions fabricated by mean of a Helium ion beam focused at sub-nanometer scale on a YBa2Cu3O7 thin film. We show that their current-voltage characteristics follow the standard Resistively-Shunted-Junction (RSJ) circuit model, and that their characteristic frequency fc = (2e/h)IcRn reaches ~300 GHz at low temperature. Using the "detector response" method, we evidence that the Josephson oscillation linewidth is only limited by the thermal noise in the RSJ model for temperature ranging from T ~ 20 K to 75 K. At lower temperature and for the highest He irradiation dose, the shot noise contribution must also be taken into account when approaching the tunneling regime. We conclude that these Josephson nano-junctions present the lowest noise level possible, which makes them very promising for future applications in the microwave and terahertz regimes.

Introduction

The astonishing recent evolution of Information and Communication Technologies (ICTs) is based on an accurate control of quantum properties of semiconductors at sub-micron scales. As some limitations appear, new paradigms emerge to further improve the performances of ICT devices, based on coherent quantum states and nano-scale engineering. Superconductivity is a very interesting platform which provides robust quantum states that can be entangled and controlled to realize quantum computation and simulation[1], or classical computation at very high speed using the so called SFQ (Single Flux Quantum) logic[2]. This platform can also be used to make detectors of electromagnetic fields and photons operating at the quantum limit, i.e. with unsurpassed sensitivity and resolution. These quantum sensors can be used for classical or quantum communications[3], THz waves detection and imaging[4], sensitive high frequency magnetic fields measurements[5,6]. Impressive results have been achieved in the recent years with devices based on Low critical Temperature (T) Superconductors (LTS) working at liquid helium temperature, and well below for Quantum Computing.

The main building block of this superconducting electronics is the Josephson Junction (JJ), a weak link between two superconducting reservoirs. While the technology for LTS JJ of typically 1 μm in size required for complex systems is mature[2], other ways are explored to downsize the JJ using Carbon Nano-Tubes[7], Copper nanowires[8] or LaAlO3/SrTiO3 interfaces[9] for examples. The complexity and the cost of the needed cryogenic systems are clearly obstacles for large scale applications of such devices. High-T superconductors (HTS) operating at moderate cryogenic temperature (≤40 K) appear as an interesting alternative solution, provided reliable JJ are available.

Different routes to make HTS JJ with suitable and reproducible characteristics are explored[10,11]. One of them relies on the extreme sensitivity of HTS materials such as YBa2Cu3O7 (YBCO) to disorder, which first reduces T and then makes it insulating. High energy ion irradiation (HEII) have been used to introduce disorder in YBCO thin films through e-beam resist masks with apertures at the nanometric scale (20–40 nm), to make JJ[12-14] and arrays[15-17] with interesting high frequencies properties, from microwaves to THz ones. Recently, Cybart et al. successfully used a Focused Helium Ion Beam (He FIB) to locally disorder YBCO thin films and make JJ[18]. In this technique, a 30 keV He+ ion beam of nominal size 0.7 nm is scanned onto a thin film surface to induce disorder. It has been used to engineer nanostructures in two-dimensional (2D) materials[19-21], magnetic ones[22] or to make plasmonic nano-antennas for instance[23]. Superconducting nano-structures and JJ have been fabricated with cuprate superconductors[24-28], MgB2[29] and pnictides[30]. While mainly DC properties of HTS JJ made by this technique have been reported to date, the present work aims at exploring their dynamic behavior, by studying the Josephson oscillation linewidth in the tens of GHz frequency range. This characterization, which gives access to the intrinsic noise of the JJ, is essential for high frequency applications of JJ such as mixers and detectors[14,16,17] and to study unconventional superconductivity[31].

Depending on the ion dose, HTS JJs made by the He FIB technique behave as Superconductor/Normal metal/Superconductor (SNS) JJs or Superconductor/Insulator/Superconductor (SIS) ones[18]. Müller et al.[28] evidenced scaling relations obeyed by the characteristic parameters I (critical current) and R (normal resistance), which are typical of highly disordered materials and known for a long time in HTS Grain-Boundary (GB) JJ[32] for example. The large density of localized electronic defect states at the origin of this behavior is a source of 1/f low-frequency noise[33,34], which broadens the Josephson oscillation linewidth at much higher frequency[35]. To assess the potential performances of HTS Josephson devices made with He FIB, we directly measured the Josephson linewidth up to 40 GHz.

Results

To fabricate HTS JJ, we begin with a commercial 50 nm-thick c-oriented YBCO thin film on a sapphire substrate[36] capped in-situ with 250 nm of gold for electrical contacts. After removing the Au layer by Ar+ ion etching everywhere except at contact pads, we structure 4 μm wide and 20 μm long channels using the HEII technique[12]. An e-beam resist mask protects the film from a 70 keV oxygen ion irradiation at a dose of 1 × 1015 ions/cm2 to keep it superconducting. The unprotected part becomes insulating. In a second step, samples are loaded into a Zeiss Orion NanoFab Helium/Neon ion microscope and the 30 keV He+ beam (current ~1.15 pA) was scanned across the 4 μm-wide superconducting bridges to form JJs. A single line is used in these experiments, whose trace can be imaged directly in the microscope[28] (Fig. 1(a)). Imaging with the He+ beam creates disorder, which adds to the one used to make the JJ. This is why we did not image the channels that we measured in the present study. On the same YBCO chip, we irradiated different channels with different doses ranging from 200 to 1000 ions/nm. The samples were then measured in a cryogen-free cryostat with a base temperature of 2 K, equipped with filtered DC lines. The RF illumination is performed via a broadband spiral antenna placed 1 cm above the chip, and connected to a generator in a Continuous Wave (CW) mode at frequency f. To measure the “detector response” described below, the RF signal is electrically modulated at low frequency (f = 199 Hz). The output signal V is measured via a lock-in amplifier synchronized on this frequency.

Figure 1

(a) Image of a JJ using the imaging mode of the He-FIB microscope. The light grey area in between dashed lines is the superconducting channel defined by HEII The horizontal line is the 600 ions/nm dose irradiated zone, which corresponds to the barrier of the JJ. (b) Sketch of the circuit used to measure the detector response signal V. A RF signal whose amplitude is modulated at the frequency f is sent onto the JJ via an antenna. V is measured with a lock-in amplifier at f. The JJ is described according to the RSJ model, as a junction in parallel with a resistance R. (c) R vs T curves (solid lines) for JJ made using different irradiation doses: 200 ions/nm (blue), 400 ions/nm (green) and 600 ions/nm (red). Same color code for panels (c,d). Below T, R (symbols) is extracted from RSJ fits. Dashed lines show the linear decrease of R(T) curves below T. (d) I vs T for different irradiation doses. Dashed lines are quadratic fits. inset: I − V characteristics of a 200 ions/nm JJ. Blue line are data and black line is the RSJ fit. (e) IR product vs T for different irradiation doses. Colored areas are calculated from the dashed lines in (c,d), and correspond to the Josephson regimes.

Figure 1(c) shows the resistance R as a function of temperature T for samples irradiated with a dose of 200, 400 and 600 ions/nm. Below the T of the reservoirs (T = 84 K), a resistance plateau develops till a transition to a zero-resistance state takes place, corresponding to the Josephson coupling through the irradiated part of the channel. Let T be this coupling temperature, which decreases as the dose is increased as already reported[18,28]. The resistance above T increases with disorder as expected, from 1 Ω (200 ions/nm) to ~3 Ω (600 ions/nm). For irradiation dose higher than 1000 ions/nm, an insulating behavior is observed down to the lowest temperature. For the samples studied here (doses between 200 and 600 ions/nm), we measured the current-voltage (I − V) characteristics below T. The inset of Fig. 1(d) shows the I − V curve of the 200 ions/nm sample recorded at T = 73 K, which can be accurately fitted with the Resistively-Shunted-Junction (RSJ) model including thermal noise (black line), as already reported[18,28].

This model accounts for Josephson weak links and Superconductor-Normal Metal-Superconductor (SNS) junctions, where the quasiparticle current is in parallel with the superconducting one, in the limit of small junction capacitance[37,38]. Finite temperature effect is introduced by mean of a noise current whose power spectral density corresponds to the Johnson noise of the normal state resistance R[35,39] (see Methods section for more detail and numerical calculation).

The RSJ fits are still valid when the dose and the temperature are varied, as proved by extended fits shown in Fig. 2. From these fits, we extracted the temperature-dependent normal-state resistance R(T) and the critical current I(T) presented in Fig. 1(c,d), respectively, with an uncertainty of typically a few percents, indicated by error bars in the figures. The former roughly follows the R(T) curve measured above T, decreases linearly with temperature (dashed lines) and goes to zero at the superconducting temperature of the irradiated part where the Josephson regime ends. The latter has a quadratic temperature dependence (dashed lines) as expected for SNS JJ[12,40]. Its absolute value can exceed 1 mA (200 and 400 ions/nm doses), which corresponds to critical current densities larger than 500 kA/cm2 [28]. We show in Fig. 1(e) the IR product as a function of temperature for the different irradiation doses.

Figure 2

I − V characteristics of (a) a 200 ions/nm JJ, (b) a 400 ions/nm JJ and (c) a 600 ions/nm JJ, with (red) and without (blue) 10 GHz irradiation. Solid lines are data and dashed lines RSJ fits. Color-plot of R as a function of I and V at f = 10 GHz for (d) a 200 ions/nm JJ and (e) a 400 ions/nm JJ (Color-scale at the bottom). (f) Same plot for a 600 ions/nm JJ (Color scale on the right).

At low doses, it shows a maximum, characteristic of SS’S junctions (where S’ is a superconductor with a T lower than the one of S) as observed with HEII HTS JJ[14]. However, for the highest dose (600 ions/nm), it raises monotonically as the temperature is lowered. Its maximum value (IR ~ 600 μV at 4 K) lies in between the values reported by Cybart et al.[18] and Müller et al.[28]. The corresponding characteristic frequency GHz ( the flux quantum) is higher than the one obtained by the HEII technique, which is promising for operations up to the THz frequency range[14].

We now focus on properties of these JJ at frequencies f much lower than f, and more specifically on the Shapiro steps which appear on I − V characteristics at voltages (n is an integer). Figure 2(a) shows the I − V curves of the 200 ions/nm JJ measured at T = 72 K without (blue) and with (red) RF irradiation at f = 10 GHz, where we observe clear Shapiro steps. Both curves are well fitted with the RSJ model (dashed lines) with the following parameters: R = 0.5 Ω, I = 133 μA and for the RF curve: I = 133 μA (the RF current). For this temperature close to T = 75 K, R does not depend on the bias current I, contrary to the HEII HTS JJ[15,41]. Sweeping both the RF voltage V and the bias current I, we recorded the I − V curves from which we computed numerically the differential resistance . The result is presented in color-scale in Fig. 2(d). The observation of well-defined and high-index (up to n = 12) Shapiro steps attests the quality of this SNS JJ. Similar measurements were performed on the other JJs. The results are shown in Fig. 2(b,c) for the 400 ions/nm JJ and the 600 ions/nm, respectively. In each case, the measurement temperature (T = 60 K and T = 30 K) are close to their respective T. In this regime, all the curves are very well fitted with the RSJ model with the following parameters: R = 0.68 Ω, I = 205 μA and I = 143.5 μA for the 400 ions/nm JJ, and R = 2.9 Ω, I = 32.5 μA and I = 26.2 μA for the 600 ions/nm one. It is important to note that the RSJ fits were performed while taking a noise temperature equals to the bath temperature. Figure 2(e,f) show color-plot for the corresponding samples. In both cases, pronounced oscillations with RF voltage corresponding to high order Shapiro steps are observed.

Shapiro steps unveil the internal Josephson oscillation that is produced when a JJ is biased beyond its critical current. The width of the steps is the linewidth of the Josephson oscillation[35]. Within the RSJ model, Likharev and Semenov[35,39] calculated the voltage power spectral density S(f) and the resulting Josephson oscillation linewidth Δf as follows:

This thermal Δf is the minimum Josephson linewidth which can be measured, as any other noise source will increase this intrinsic linewidth. Divin et al. showed that Δf can be measured experimentally from the Shapiro steps by mean of the “detector response” method[42-44]. The JJ is DC biased while the RF illumination is modulated at low frequency (f)[45]. The “detector response” signal V is measured with a lock-in amplifier synchronized at f, and plotted as a function of the DC voltage V converted into a frequency f through the Josephson relation . Centered on the Josephson frequency, i.e. on the Shapiro step, an odd-symmetric structure appears, whose width (distance between the extrema) corresponds to Δf. To be more precise, Divin et al.[42] showed that the inverse Hilbert transform of the quantity is directly S(f), a Lorentzian of width Δf centered at the Josephson frequency f. This procedure, successfully used in LTS[43] and HTS[44] materials, allows to extract the Josephson linewidth accurately.

We measured the Josephson oscillation linewidth of the different JJ irradiated at f = 40 GHz using the “detector response” method. Figure 3(a) (bottom panel) shows V (left axis) as a function of the frequency () for the 200 ions/nm JJ. Around 40 GHz, the characteristic double-peak structure predicted by Divin[42] is observed, from which we extracted the Josephson linewidth (Δf = 4.87 GHz). It is worth noticing that this method is highly sensitive, since the Shapiro steps cannot be seen in the I − V curve simultaneously recorded (top panel). We then computed S through the above explained procedure, and plotted it on the same graph (Fig. 3(a) (bottom panel, right axis). Two peaks are observed, corresponding to the first an second Shapiro steps, that can be fitted with Lorentzian to extract the corresponding Josephson oscillation linewidths. For the first step (index n = 1), the value is exactly the same as the one calculated above. Depending on the experimental conditions, we could accurately measure the Josephson linewidth of the first two Shapiro steps, or only of one of them.

Figure 3

(a) I − V characteristics of the 200 ions/nm JJ under f = 40 GHz irradiation. V vs . The distance Δf between the extrema corresponds to the Josephson oscillation linewidth. S extracted from the inverse Hilbert transform of the the normalized response g(V) (pink), whose width is Δf (grey). (b) V vs V for a 400 ions/nm JJ measured at different temperatures under f = 40 GHz irradiation. (c) Δf vs T for different irradiation doses. Solid symbols correspond to the first Shapiro step, open symbols to the second one. Solid lines (blue, green, red) are calculated from the RSJ model for n = 1 (200 and 400 ions/nm) and n = 2 (600 ions/nm), the black one from the tunneling one. (d) Δf vs Δf (left panel) and vs Δf (right panel, 600 ions/nm JJ (n = 2) for ). The slope of the dashed lines is one. inset: Δf vs I for the 600 ions/nm JJ (n = 2) for . Dashed lines are best fits with a power-law exponent 1 (purple) and 1.35 (orange).

We measured V as a function of V at different temperatures, as for example reported in Fig. 3(b) for the 400 ions/nm JJ. The odd-symmetric structure at V = 82.7 μV (corresponding to 40 GHz, dashed line) widens with increasing temperature as expected for thermal noise. We extracted Δf as a function of temperature for the different samples. The result is shown in Fig. 3(c). Open (respectively solid) symbols correspond to measurements on the first (respectively second) Shapiro step. On the same graph, we added the linewidth Δf calculated for the thermal noise in the RSJ model using Eq. 1[35,39], with no adjustable parameter. The agreement is excellent for the 200 and 400 ions/nm JJ at all temperatures. For the 600 ions/nm, data are well reproduced at high temperature, but strongly depart from the calculation below T = 20 K. In Fig. 3(d), we made a parametric plot of the same data: the experimental Δf as a function of the calculated Δf in the RSJ model with thermal noise (left panel). All data align along the dashed line of slope 1, which means that noise in He FIB JJ is purely thermal, except for 600 ions/nm JJ at low temperature.

This indicates that an extra source of noise takes place below T = 20 K in this JJ. We notice that this temperature corresponds to an up-turn in the R(T) curve (Fig.1(c)).

This thermally activated electronic transport, characteristic of a disorder-induced Anderson insulator where charge carriers hop between localized states[46], is well known for ion-irradiated cuprates[47]. It has been reported by Cybart et al.[18] in YBCO JJ made by the He FIB technique for a sample slightly more irradiated than our 600 ions/nm one. They showed that a SIS junction is formed, and they observed a structure in the conductance related to the superconducting gap, as expected in tunnel junctions where the differential conductance is proportional to the Density of States of the reservoirs in first approximation. It is worth noting that in this regime, and contrary to the SS’S one, both I and R increase as the temperature is lowered, and so does the IR product (see Fig. 1e)) to reach an interesting high value (~600 μV). In that case, the tunneling approach proposed by Dahm et al.[48] is more appropriate than the RSJ one to calculate the Josephson oscillation linewidth, which includes the non-linear superposition of thermal and shot noises in these JJ at intermediate damping[35,48]. They show that:

We calculated Δf for the 600 ions/nm JJ with this expression, and obtained a very good agreement with the data as shown in Fig. 3(c)(black line), once again with no adjustable parameter. The excess noise comes therefore from the shot noise contribution when approaching the tunneling limit. The parametric plot of the experimental Δf as a function of the calculated Δf including the shot noise (Fig. 3(d) right panel) clearly shows that there is no additional noise source in our JJ.

Discussion

This seems quite surprising, given the recent result from Müller et al. pointing towards highly disordered JJ[28], which are usually associated with strong 1/f Flicker noise[33]. They measured the low-frequency noise of SQUIDs made with low-irradiation dose JJ (230 ions/nm). A clear 1/f noise component is observed up to ~100 kHz. These measurements have been performed on JJ fabricated on the so-called “LSAT” substrate, which have much lower I and IR products than the others for the same dose. The role of the substrate on the JJ characteristics is not understood yet, but it is clear that the microstructure of the film matters for final JJ performances, and more specifically for noise properties related to defects. This has been evidenced long ago on YBCO Grain-Boundary JJ by annealing experiments[49]. Our samples grown on sapphire may have therefore less fluctuating centers at the origin of Flicker noise than others. Past studies showed that this noise in HTS JJ is often induced by enhanced critical current fluctuations in inhomogeneous barriers[50-52], and that the maximum 1/f noise power in the vicinity of I scales with it ()[52,53]. This noise translates into a broader Josephson linewidth at high frequency, especially when f < f, i.e. for DC bias current close to I[44]. It has been evaluated by Hao et al.[52] as:where f0 is a low frequency cut-off of the 1/f noise (typically f0 ~ 1 Hz). Through the above mentioned scaling relation, Δf should thus increase as or so. The data on our 600 ions/nm sample below T = 20 K are compatible with this relation (inset Fig. 3(d)), but we cannot make a quantitative fit since we do not know the value of S(f0). Moreover, the same data can be fitted with the shot noise model as well (inset Fig. 3(d)), since the latter states that Δf ~ I close to I (derived from the above equation). This model fits quantitatively the evolution of Δf with temperature with no adjustable parameter, and qualitatively the one with the critical current. We therefore conclude that the shot noise contribution fully explains the low temperature data of the most irradiated JJ, and that there is no evidence of strong Flicker noise in the present study.

Conclusion

We fabricated HTS JJ by the He FIB technique, and studied their DC and RF properties in the 10 to 40 GHz range. Their IR product reaches 300 GHz at low temperature, which is higher than for HEII JJ. We showed that Shapiro steps in the I − V characteristics that appear under RF irradiation are well described by the RSJ model for SNS junctions with thermal noise. Using the “detector response” method, we determined the Josephson oscillation linewidth, and showed that it corresponds to the sole Johnson-Nyquist thermal noise in the RSJ model for the low-dose irradiated JJ. Below T = 20 K, the high-dose irradiated sample has a SIS character. We demonstrated that the associated enhanced noise is due to shot noise when approaching the tunneling regime. We did not evidenced any Flicker noise component, which means that the barrier is rather homogeneous in these JJ. This study paves the way for using He FIB JJ in high frequency applications[53].

Methods

Resistively shunted junction (RSJ) model

Calculation of the I–V curve at finite temperature

The Resistively Shunted Junction (RSJ) model describes the equivalent circuit of a Josephson Junction (JJ) as two elements in parallel (the junction described by the two Josephson equations written below and its normal state resistance R as sketched in Fig. 1(b)), biased with a current I[37,38]. In this “overdamped” limit, the capacitance of the junction is neglected. The Josephson equations state:where I is the bias current and V the voltage of the JJ, I its critical current, φ the quantum phase difference across it, and the superconducting flux quantum.

The time evolution of the current is therefore:

The voltage is given by Eq. 5.

These equations are valid in the limit of zero-temperature. At finite temperature, the Johnson noise of the resistance must be added. The power spectral density of the current fluctuations at temperature T is[35,39]:

We introduce a noise current δI(t) whose power spectral density is given by Eq. 7. It has therefore a Gaussian variation in time with a variance:where Δt is the time interval considered.

The time evolution of the current is now:

To get the DC I-V curve, one needs to time average this equation.

I–V curve in a presence of RF irradiation

If the JJ is submitted to a time varying current (RF irradiation) at finite temperature, the time evolution of the JJ is given by:and Eq. 5. After time averaging, the I-V curves present current (Shapiro) steps at voltages , where n is an integer. The width of the transition from one step to the next one is given by the thermal noise.

I–V curve simulation

In practice, the simulation of an IV curve consists in solving the Eqs. 5 and 10 by numerical integration using the Euler method. Hence the system to be numerically solved is:where the bracket notation means discrete time steps of pace δτ, and n is the step index. For each current bias I, a voltage vector is thus found by iteration, for each step δτ from τ = 0 to τ = τ, starting with a random initial phase and V[0] = 0. The Gaussian noise δI is a random variable changed at every step, whose variance is given by Eq. 8 (Δt = δτ).

The system is numerically heavy to solve: first because one needs a sufficiently small δτ to account for the rapid variation of the voltage oscillations, especially at low bias, and at the same time one needs a sufficiently high τ in order to have enough oscillations to average. In practice, δτ must be much smaller than 1/v, and τ should be sufficiently high to average enough oscillations. We typically have vectors of 200000 points, and δτ ~ 1 · 10−12 s. Second, because the presence of the (actually pseudo random) noise also requires to average the calculation of each V over several iterations of the same IV curve, typically 10 times.