Weak spin-flip scattering in Pd89Ni11 interlayer of NbN-based ferromagnetic Josephson junctions.

We studied niobium nitride (NbN)-based π-junctions with a diluted ferromagnetic Pd89Ni11 interlayer (NbN/PdNi/NbN junctions). In the NbN/PdNi/NbN junctions with various PdNi thicknesses, we observed a non-monotonic dependence of the critical currents on PdNi thickness, indicating the effects of the exchange interaction on the superconducting order parameter. From theoretical fitting of the experimental data, we found that the NbN/PdNi/NbN junctions showed a significantly smaller degree of spin-flip scattering in the PdNi interlayer than in the CuNi interlayer of NbN/CuNi/NbN junctions reported previously. The weak spin-flip scattering leads to a longer decay length of the Josephson critical current, so the critical currents were observed over a wide range of PdNi thicknesses (10-40 nm). We also fabricated superconducting quantum interference devices (SQUIDs) including the NbN/PdNi/NbN junction, using a PdNi thickness in which the π-state was expected. A half-flux-quantum shift, as evidence of the π-state, was observed in the magnetic field-dependent critical currents of the SQUIDs. This result represents an important step towards the practical application of NbN-based π-Josephson junctions.

Introduction

In a Josephson junction, two superconductors are coupled by insulating or normal metal interlayers (SIS or SNS junctions), and the superconducting order parameter shows an exponential decay through the interlayers[1]. On the other hand, when the interlayer consists of ferromagnetic materials (SFS junctions), the superconducting order parameter shows spatial oscillation in the ferromagnetic layer due to the interplay between superconductivity and ferromagnetism[2]. Thus, in contrast to conventional junctions, SFS junctions can be in a so-called π-state over a certain range of ferromagnetic layer thicknesses, in which the phase difference between the two superconducting layers becomes π in the ground state. The π-state in an SFS junction (π-junctions) is an interesting topic, not only due to its intriguing physics but also its potential applications.

In the field of quantum computing, superconducting quantum bits (qubits) are proving to be one of the best candidates for building a practical quantum computer[3]. Flux qubits with long coherence times and high anharmonicity have been reported, to facilitate the scaling-up of quantum processors[4]. Conventional flux qubits always require a half-flux-quantum bias to operate at the optimum point for achieving the longest coherence time. However, the necessity of having a half-flux-quantum bias reduces scalability[5]. One way to solve this issue is to insert a π-junction into a superconducting loop of the flux-qubit. A coherent two-level state can be created with a zero external magnetic field, so that the flux qubit can be operated optimally without half-flux-quantum bias[6-8]. Recent research into superconducting single-flux-quantum logic circuits has demonstrated the possibility of improving operating margins and energy consumption by introducing π-junctions into these circuits[9-11]. However, despite extensive study of π-junctions for superconducting quantum and classical circuits[2,12-14], there have been few reports on device implementation of the π-junction to date[15-17]. Therefore. further research to advance the development of the π-junction toward practical application are necessary.

Niobium (Nb) is often used as a superconducting material for π-junctions. From the viewpoint of superconducting logic circuits application, niobium nitride (NbN) is an attractive material, as it provides a high operating frequency of up to 1.2 THz and an operating temperature of 10 K, which cannot be achieved with Nb electrodes[18]. In addition, its large superconducting gap of 5.2 meV can suppress noise caused by quasiparticle excitation in the junctions, which is necessary for enhancing the coherence time of qubits[19]. For NbN-based ferromagnetic Josephson junctions, the ferromagnetic CuNi interlayer has been studied, and a π-phase shift in NbN/CuNi/NbN junctions has been demonstrated[20,21]. However, the CuNi-alloy contains Ni magnetic clusters and shows strong spin-flip scattering; thus, the Josephson critical current of the junctions with the CuNi interlayer decays quickly with increasing CuNi thickness[14,20]. On the contrary, in PdNi interlayers, Josephson critical currents have been observed over a wider range of PdNi thicknesses in Nb-based junctions, even in those using a 100-nm-thick PdNi interlayer[22-24]. From the viewpoint of device applications, the PdNi interlayer also provides advantages for NbN-based ferromagnetic Josephson junctions, including easier control over the critical currents and ensuring the π-state of junctions in large-scale circuits, even if there are spatial and/or run-to-run variations in PdNi thickness.

In this work, we studied NbN-based Josephson junctions with a diluted ferromagnetic Pd89Ni11 interlayer (NbN/PdNi/NbN junctions) and demonstrated their π-phase shift for the first time. We observed a non-monotonic dependence of critical currents (Ic) on the PdNi layer thickness, indicating phase transition between the 0- and π-states owing to oscillation of the superconducting order parameter in the PdNi ferromagnetic interlayers. From the fitting of the experimental results based on an expression derived from the Usadel equation, we found that the spin-flip scattering degree in the PdNi interlayer was much smaller than that in the CuNi interlayer, leading to a longer decay length of the Josephson critical current. Based on these results, superconducting quantum interference devices (SQUIDs) were fabricated with SIS and SFS junctions. A clear half-flux-quantum shift was observed in the dependence of the critical current on the external magnetic field, thus providing evidence of the π-state.

Fabrication and experimental setup

The MgO (100) was used as a substrate for NbN (100) film growth[25,26]. First, an NbN/PdNi/NbN trilayer was deposited in a multi-chamber sputtering system. The system was equipped with two separate sputtering chambers for NbN and PdNi films, in which the chambers were connected through a load lock chamber. A 200-nm-thick NbN was deposited as a base layer via reactive direct current (dc)-sputtering of an Nb target under an argon (Ar) and nitrogen gas atmosphere[20]. Then, the sample was transferred to the other chamber without breaking the vacuum, for deposition of the PdNi barrier by dc-sputtering with a Pd89Ni11 target. The thickness of the PdNi layer ranged from 10 to 40 nm by controlling the sputtering time. Next, the sample was sent back to the sputtering chamber to deposit the 200-nm-thick counter electrode NbN layer. All the sputtering processes are conducted at room temperature without heating, so the crystal or magnetic structure of PdNi film will not be affected by the high-temperature deposition.

The junctions were patterned in a square shape (size: 10 × 10 µm2) using conventional photolithography, and the counter electrode NbN layer was etched via reactive ion etching (RIE) with CF4 gas. The PdNi layer was fabricated by physical etching with Ar gas. Then, the base NbN electrode was patterned and etched by RIE. A 300-nm-thick SiO2 layer was deposited for isolating the base and NbN counter electrodes. The contact hole was created in the SiO2 layer and etched using RIE with CHF3 gas. Finally, a 450-nm-thick wiring NbN layer was deposited and patterned. Figure 1a,b show a schematic cross-sectional view and microscope image of the junction, respectively.

Figure 1

(a) Schematic cross-sectional view of the NbN/PdNi/NbN junction on an MgO substrate. (b) Microscope image of 10 × 10 µm2 junction.

In addition to the single NbN/PdNi/NbN SFS junctions, we fabricated SQUID structures including two SIS junctions with an SFS junction, to determine the fabricated junction was in the π-state. First, the NbN-based SIS (NbN/AlN/NbN) junctions were fabricated. Then, the SFS junction was constructed to form the SQUID structure[21]. Although the fabrication process of the SFS junction was nearly the same as that for the single junction, the junction part was defined by the lift-off technique to avoid damage to the underlying structure, including the SIS junctions.

The current–voltage (I–V) characteristics of the junctions and magnetic field dependence of the critical currents of the SQUIDs were measured at 4.2 K using a magnetically shielded cryostat filled with liquid helium.

Results and discussion

To investigate the magnetic properties of the PdNi interlayer, we deposited a 35-nm-thick PdNi film on a Si substrate and measured the temperature and magnetic field dependences of the magnetization using a commercial magnetic properties measurement system (MPMS; Quantum Design). Figure 2a shows the temperature dependence of the magnetization (M–T) for in-plane and out-of-plane magnetic fields at 10,000 Oe. To eliminate the signal of the substrate, we measured an Si substrate having the same size as that used for PdNi growth and subtracted the measurements from the data of the PdNi/Si sample. We confirmed that the jump of the magnetization around 10,000 Oe for the in-plane field was originated from the Si substrate. The Curie temperature (TCurie) was approximately 163 K both for in-plane and out-of-plane magnetic fields (Fig. 2a). Figure 2b shows the magnetic field dependences of the magnetization (M–H). From the M–H curves, the saturation magnetization (Msat) of the PdNi film for the out-of-plane field was determined as 120 emu/cm3. The obtained values of TCurie and Msat indicated that the Ni concentration of the deposited PdNi film was around 11%, which was the same as that of the PdNi target[24,27]. In a previous study of NbN/CuNi/NbN junctions, the Ni concentration obtained from TCurie, CuNi = 200 K was ~ 59%, which was higher than that of the CuNi target (53%)[20]. The M–H curve showed clear out-of-plane magnetic anisotropy. For ferromagnetic thin films, the shape anisotropy usually dominates the magnetocrystalline anisotropy, resulting in in-plane anisotropy. However, in very thin films, the magnetic surface energy may be dominant, leading to out-of-plane anisotropy[28]. To date, the out-of-plane anisotropy of PdNi thin films with a thickness of up to 60 nm has been reported, as well as that for CuNi films with a thickness of 34 nm[23,29]. It is expected that the stray field from the PdNi interlayer could be reduced in SFS junctions compared to those related to in-plane anisotropy films[30].

Figure 2

(a) Temperature-dependent magnetization (M–T) of a 35-nm-thick PdNi film at 10,000 Oe. (b) Magnetic field-dependent magnetization of the 35-nm-thick PdNi film at 4.2 K.

Next, we measured the I–V characteristics of the NbN/PdNi/NbN junctions at 4.2 K. Figure 3a shows the I–V characteristics of the junction with a 25-nm-thick PdNi interlayer. Overdamped I–V characteristics, which are typical for SFS junctions, were observed[20]. The normal resistance of the junction (R) is 256 µΩ with 25-nm-thick PdNi. This value is higher than that extracted from the resistivity of a single Pd89Ni11 film (~ 90 µΩ)[31], which may be due to the interface resistance arising in the non-epitaxial NbN counter-layer grown on PdNi. However, this interface resistance is much smaller than that reported for NbN/CuNi/NbN junctions[32], in which the interface resistance was seven-fold larger than the resistance extracted from the resistivity of a single CuNi film.

Figure 3

(a) Current–voltage characteristics of a 10 × 10 µm2 junction with a 25-nm-thick PdNi interlayer measured at 4.2 K. (b) Dependence of critical currents on the thickness of the PdNi interlayer. The solid line denotes the theoretical fitting of the data, where 0 and π indicate the expected 0 and π regions based on the fitting curve, respectively.

Figure 3b shows the dependence of Josephson critical current (Ic) on the PdNi thickness. The obtained Ic (symbols) showed non-monotonic dependence, which may indicate a phase transition between 0- and π-states. To analyze the experimental data, we considered the following theoretical expression of the critical currents in diffusive SFS junctions, derived from the Usadel equation[14,20]. where I is a constant prefactor, T = 15 K is the critical temperature of NbN, T = 4.2 K is the measurement temperature, and ξ is the coherence length in the ferromagnetic interlayer. n is an integer, = , = (i + )/(i +  +  ), and = /[ + ]2. Δ(T) is the superconducting gap, = /E = π(2n + 1)kT/Eex, and α = ℏ/(τsE), where E is the exchange energy and τs is the spin-flip scattering time in the ferromagnetic layer. k is Boltzmann’s constant and ℏ is the reduced Plank’s constant. The parameter α indicates the degree of spin-flip scattering. In Eq. (1), the free fitting parameters are E, ξ, α, and I. The solid line in Fig. 3b indicates the curve fitted to the experimental data. From the fitting, we obtained the values as E/k = 175 K (E 15 meV), ξ = 3.5 nm, α = 0.1, and I = 1.8 × 108 µA. The values of E 15 meV and Δ = 5 meV indicate our junctions fulfill the diffusive limit of Eq. (1), since E, Δ ˂˂ ℏ/τe = 78 meV, where τe ≈ 8.4 × 10–15 s is the electron collisions time estimated from the resistivity and electron mean-free path (le) of PdNi[23].

In order to evaluate the relevance of Eq. (1) when apply to our data, we estimated the interface transparency Ƭ of PdNi junctions using the relation[33]: ϒB = . Here, the parameter ϒB = (RBS/ρFξ*), with RB is interface resistance per unit area, S is the SFS junction area, ρF is the resistivity of PdNi layer and ξ* = . In our junction, RB = 83 µΩ for 10 × 10 µm2 junction (estimated from junction’s resistance and resistance extracted from resistivity of PdNi film), ρF = 36 µΩ cm, ξ* = 4.8 nm (calculated from ξ*/ξ ratio), so ϒB = RBS/ρFξ* ≈ 4.8. The electron mean-free path le ≈ 4.7 nm was estimated using the relation ρFle = 1.7 fΩ m2, which has been reported for PdNi[23]. Finally, we obtained the interface transparency Ƭ = 0.12. The condition Ƭ = 0 corresponds to a completely reflecting interface, whereas Ƭ = 1 corresponds to a perfect transparent boundary. Experimentally, it is not practical to get the perfect transparency with Ƭ = 1. Our result of Ƭ = 0.12, however, equal to the transparency value calculated for CuNi junctions (with ϒB = 0.52, le = 1 nm, ξ* = 9.4 nm) in the reference[34], where Eq. (1) has been used with good approximation. Therefore, our PdNi junctions fulfill the condition of interface transparency for this expression.

Here we discuss the obtained fitting parameters. First, the exchange energy E of 15 meV is in good agreement with that for the Pd89Ni11 alloy in a previous report[27]. The value of ξ = 3.5 nm is also consistent with values from the literature: ξ = 3.3–4.0 nm for Nb-based SFS junctions with a PdNi alloy with similar Ni concentrations[23,27,35]. Next, we compared the obtained parameters with those for the NbN/Cu40Ni60/NbN junction[20]. The value of ξ (3.5 nm) was larger than that of 1.9 nm reported for NbN/CuNi/NbN junctions, due to the smaller E and larger diffusion constant D in for the NbN/PdNi/NbN junctions[13]. Regarding the degree of spin-flip scattering, we found that the value of α was 0.1, which is much smaller than that of 0.8 for the NbN/CuNi/NbN junction[20]. We calculated the decay length (ξ) and oscillation length (ξ) of the Josephson critical current Ic(d) using the equation ξ (T) = ξ/, and obtained the following[20]: ξ(0) = 3.3 nm and ξ(0) = 3.9 nm. Owing to the smaller α and larger ξ, the decay length ξ in the NbN/PdNi/NbN junction became much larger than that of 1.3 nm in the CuNi junctions. Furthermore, ξ(0) and ξ(0) were comparable due to the small value of α, and the oscillation length ξ(0) was also larger than that of the 2.8 nm reported for NbN/CuNi/NbN junctions[20]. This small value of α, i.e., weak spin-flip scattering, in the PdNi interlayer may originate from better magnetic uniformity with fewer magnetic impurities (e.g., clusters). It is known that, in CuNi films, ferromagnetism originates from clusters of at least eight Ni atoms, due to statistical concentration fluctuations[36]. On the other hand, in the PdNi film, the ferromagnetism is caused not only by the diluted Ni atoms but also by paramagnetic Pd atoms. This leads to homogeneous ferromagnetism in PdNi rather than a matrix with ferromagnetic clusters in CuNi. From a device application perspective, the longer ξ and ξ obtained in the PdNi junction are attractive, as it becomes easier to fabricate a π-junction with the desired critical current, even under the conditions of process variability and/or spatial variation of the PdNi thickness. Lastly, the obtained value of 1.8 × 108 µA for I is also comparable to that calculated from the actual physical parameters of the fabricated junctions (T and R) using the expression obtained from the Usadel equation[14,20].

It is necessary to discuss the role of the spin–orbit scattering in PdNi interlayer since Pd is a heavy metal. In SFS junctions, the spin–orbit and spin-flip scattering have a similar effect on ξ and ξ. If the spin–orbit scattering rate is larger than Eex/ℏ, the oscillations of the junction’s critical currents disappear completely[23]. Our results in Fig. 3b, however, shows a clear oscillation of critical currents, so the spin–orbit scattering rate was expected to smaller than the limit above.

From the fitting of the experimental data, we can expect a π-state region with a PdNi thickness in the range of 8–20 nm, as shown in Fig. 3b. To determine whether the π-state actually emerged, we selected SFS junctions with 10- and 15-nm-thick PdNi interlayers, and fabricated and characterized SQUID structures with these junctions as follows.

On a chip, we fabricated two SQUIDs: a conventional SQUID with two SIS junctions (dc-SQUID) and a SQUID with two SIS junctions and one SFS junction (an SFS-SQUID), as shown in Fig. 4a,b, respectively. The SIS junctions were NbN/AlN/NbN junctions with critical currents of 15 µA and a size of 5 × 5 µm2, giving a critical current density of 60 A/cm2, which is typically used for superconducting qubits[19]. The SFS junctions were circular NbN/PdNi/NbN junctions with a diameter of 5 µm. The critical currents of the SFS junctions were 3.0 and 2.9 mA for 10- and 15-nm-thick PdNi, respectively. The critical current density of the SFS junctions with both 10- and 15-nm-thick PdNi layers was around 15 kA/cm2, which was about three-fold higher than that of the SIS junction. Therefore, the SFS junction worked as a passive π-phase shifter, even with its much smaller size[10,21]. For the dc-SQUID, when the inductance of the loop is negligible, the magnetic field-dependent critical currents can be described by the equation: = 2 │cos(π)│. Here, Ф is the externally applied magnetic field and Ф0 is the single flux quantum. In the case of the SFS-SQUID, when the π-junction is integrated into the loop having a critical current higher than that of the SIS junctions ( > ), it plays the role of a π-phase shifter, resulting in the π-shifted relation[37,38]: = 2 │sin(π )│.

Figure 4

Microscope images of (a) a direct-current superconducting quantum interference device (dc-SQUID) and (b) a superconductor/ferromagnetic/superconductor (SFS)-SQUID with superconductor/insulator/superconductor (SIS) and SFS junctions, as denoted by the dotted squares and circles, respectively.

Figure 5a,b show the dependences of the positive and negative critical currents on the magnetic coil currents for the dc- and SFS-SQUIDs with a 10-nm-thick PdNi interlayer, respectively. The dashed lines represent the border for the vortex state and intersection between dotted lines (two connected peaks of ), and the horizontal axis indicates the quantum number of the SQUID. Our results show that = 2 at zero magnetic coil current for the dc-SQUID, as expected. When the coil current increased (or decreased), local minima appeared in , corresponding to the half-flux quantum value (± Ф/2) or a π-phase difference between the two arms of the SQUID loop. The peaks and dips in the modulation pattern can be understood as constructive and destructive interference of the currents in the two arms of the SQUID, respectively. At the half-flux quantum value of the magnetic fields, the drops to around zero due to the small inductance of the SQUID loop. We also calculated the inductance of the SQUID loop from the geometric inductance (L = 1.25 µd ≈ 15.7 pH) and kinetic inductance (L = µ0λ2l/wh ≈ 0.3 pH). Here, µ0 is the vacuum permeability, d is arm’s length of SQUID, λ is the penetration depth of NbN, l is total loop’s length, w is arm’s width, and h is thickness of superconducting film[37]. The estimated value of the loop inductance was L = 16.0 pH, giving a screening parameter =  = 0.2. There is a relationship between β and the depth of the modulation pattern, at β = 0.2, only 10% of the maximum occurs, corresponding to around 3 µA for this dc-SQUID[37]. This very small critical current is usually suppressed by external noise, thus, giving the results shown in Fig. 5a.

Figure 5

Magnetic field-dependent critical currents of (a) the dc-SQUID and (b) SFS-SQUID. The arrows indicate the field value at a half-flux-quantum.

On the other hand, for the SFS-SQUID shown in Fig. 5b, reached its minimum at zero magnetic coil current. As the coil current increased (or decreased), the local maxima of = 2 appeared at Ф =  ± Ф0/2. The intersections between the horizontal axis and dotted lines were equally shifted from the origin by a half period of the modulation pattern. Compared to the dc-SQUID, a clear half-flux-quantum shift was observed in the modulation pattern of the SFS-SQUID[21,38]. A half-flux-quantum shift was also observed in the SFS-SQUID with a 15-nm-thick PdNi interlayer. These results provide evidence of the π-state in SFS junctions. In Fig. 5b, a reduction of the modulation depth in the SFS-SQUID compared to the dc-SQUID can be seen, due to the increased inductance of the SQUID loop caused by the kinetic inductance of the SFS junction.

We also fabricated the SFS-SQUIDs with Nb/AlO/Al/Nb junctions. In these samples, the Nb/AlO/Al/Nb junctions were fabricated on Nb base electrode by using the standard fabrication process of clean room for analog and digital superconductivity (CRAVITY) in the National Institute of Advanced Industrial Science and Technology (AIST) Japan. Then, we fabricated the NbN/PdNi/NbN junctions with 15- and 35-nm-thick PdNi on the polycrystalline Nb base layer of the SQUIDs. A π-phase shift was observed in the modulation pattern of these SFS-SQUIDs. The π-phase shift demonstrated by the NbN-based SFS-SQUIDs and Nb four-layer structure SFS-SQUIDs expands its application fields.

Conclusions

In conclusion, we studied NbN/Pd89Ni11/NbN junctions with a PdNi thickness of 10–40 nm and observed non-monotonic behavior in the thickness-dependent critical currents of the junctions, indicating the oscillation of the superconducting order parameter due to the exchange interaction in the PdNi ferromagnetic layer. Spin-flip scattering in the PdNi junctions was significantly weaker than that in CuNi junctions, leading to a longer decay length. This provides advantages for the fabrication and performance of the device, such as ease in controlling the critical currents and ensuring working π-junctions for large-scale circuits, even with spatial variation in the PdNi thickness. Lastly, we fabricated and measured the magnetic field-dependent critical currents of dc-SQUIDs and SFS-SQUIDs for demonstrating the π-state of junctions. A half-flux quantum shift was observed in the modulation pattern of SFS-SQUIDs compared to the dc-SQUID, indicating the presence of π junctions. These results not only provide an important understanding of the physics of ferromagnetic Josephson junctions but also an initial step toward practical application.