Granular Aluminum kinetic inductance nonlinearity

Granular Aluminum is a superconductor known for more than eighty years, which recently found its application in qubits, microwave detectors and compact resonators, due to its high kinetic inductance, critical magnetic field and critical current. Here we report on the nonlinear dependence of granular Aluminum inductance on current, which hints towards parametric amplification of the microwave signal in granular Aluminum films. The phase shift of the microwave signal reached 4 radians at a frequency of 7 GHz, which makes it possible to estimate the nonlinearity of the system as 1.4 % and the potential gain of the order of 17 dB


I. INTRODUCTION
Superconductors are one of the most promising platforms for quantum processor realisation.Transmonbased processors have been used in first experiments on quantum supremacy demonstration [1,2], quantum error-correcting surface codes realizations [3,4] and solving quantum chemistry problems [5,6].Such processors require an amplification of a readout signal.Superconducting parametric amplifiers are the best in terms of signal-to-noise ratio, having an additional noise at the quantum limit level: 168 mK at 7 GHz [7].
Traveling-wave parametric amplifiers (TWPA), used for a wide-band parametric amplification, employ either an array of Josephson junctions [8][9][10], or a high kinetic inductance materials like NbTiN [11][12][13], TiN [14] or NbN [15].Josephson junctions fabrication consists of several complicated technological operations implementing a multi-layered structure of the junctions, which makes fabrication of a 1000 junction array [8,[16][17][18] a rather complex task.On the other hand TWPA based on high kinetic inductance films [19,20] reduce fabrication to one technological cycle and minimize production risks.Nevertheless a nonlinearity of such systems is concomitant with local heating and hot spots emerging, which leads to the structure's critical current decrease [21,22].
Here we consider grAl as an alternative to traditional in TWPA realizations NbTiN [11][12][13]19].The microstructure of grAl consists of pure aluminum grains separated by thin aluminum oxide barriers, forming a 3D-network of Josephson junctions, which results in high kinetic inductance [42].Such a microstructure is obtained by evaporating of Aluminum in Oxygen atmosphere.An * zhdanova.mar@gmail.comadditional advantage of this material is a critical temperature on the order of 2-3 K (achievable by pumping He 4 and significantly higher then Aluminum T Al c ≈ 1.2 K), high critical fields [43], and handy stable fabrication technology.

II. FABRICATION
By adjusting an oxygen flow during Aluminum evaporation, one can control a resistance per square of a grAl film [44].Depending on the flow a resistivity varies between 10 and 10 4 µΩ • cm and critical temperature goes up to 3 K [44].
The samples for this study were fabricated in an ISO-7 cleanroom using thin film planar technology methods.Films were deposited by electron beam evaporator Plassys MEB 550 S. The samples topology was formed by an optical lithography using a photoresist MICROPOSIT S1805.After forming the protective mask, the substrate was placed back in Plassys MEB 550 S, where uncoated grAl areas were etched by an ion beam.This stages of sample fabrication are schematically shown in Fig. 1a.Finally, substrates were diced in 5x5 mm peaces and all organic residuals were washed off.Using test structures of etched bridges 8 µm wide and 200 µm long, the resistance per square for each oxygen flow value was measured by four-probe method resulting in a calibration curve shown in Fig. 1b.For this study we choose resistance of R □ n ≈ 200 Ω, corresponding to 1.6 sccm O 2 flow, for optimizition of the phase velocity and characteristic impedance of transmission line.

III. EXPERIMENTAL STUDY OF GRANULAR ALUMINUM NONLINEARITY
Considering thin grAl films we have a kinetic inductance L k dominating over geometric L g .Moreover, for high resistances R kinetic inductance can be several orders greater then geometric [45].A superconductor kinetic inductance can be found in the frame of the Bardeen-Cooper-Schrieffer (BCS) theory.In the limit of low frequencies (hf ≪ k B T c ) according to Mattis-Bardeen formula for dirty superconductors the complex conductivity can be written in terms of the ratio of the imaginary conductivity σ 2 to the normal state conductivity σ n as: where ∆ is the superconducting energy gap depending in general on temperature and bias current.Introducing a kinetic inductance as L k = 1/2πf σ 2 , one can write it down for a strip of length l and width w where N = l/w is a number of squares in a strip and R □ n is its sheet resistance in the normal state.In the limit of low temperatures (T ≪ T c ), BCS formula for the energy gap ∆ = 1.76πkB T c can be used and the standard expression [39] can be derived Following a method given by [46] a kinetic inductance dependence on bias current at T ≪ T c is considered as where L k (0) is given by Eq. ( 3) and ζ = Dk 2 2∆(I) is a combination of the diffusion coefficient D, the gradient of the phase along the length of the superconducting strip k and the enrgy gap ∆(I), corresponding to bias current I.For each bias current value ζ can be found by solving following equation where I dep is the deparing current from BCS. Eqs. ( 4) and ( 5) are obtained by solving Gorkov equation in the dirty superconductor limit [47].On the other hand a kinetic inductance in general nonlinearly depends on a bias current due to the Cooper pairs breaking and can be estimated from Pippard's equation as where L k (0) is given by Eq. ( 3), and I * is of the order of the deparing current [11,46].This estimation is valid at T ≪ T c , taking into account that the linear term must vanish due to symmetry considerations.A nonlinearity of grAl opens up a possibility to use a superconducting transmission line made of grAl for parametric interaction of microwave signals.For this study we fabricated a coplanar transmission line 78 mm long with a central strip width of 3.5 µm and a gap of 5 µm, on an oxidized silicon substrate (see Fig. 2c).We used a film of thickness t = 30 nm with normal state sheet resistance R □ n = 200 Ω and critical temperature T c = 2.15 K. Permittivity of our substrates was ϵ r = 11.9.The input and output of the transmission line were connected to the contact pads via impedance converters.A characteristic impedance of a coplanar is where L kl and L gl are geometric and kinetic inductances per unit length and C l is a capacitance per unit length.A phase velocity in a transmission line can be written as Using an experimental setup shown in Fig. 2a transport properties of our samples were measured.The grAl transmission lines were placed in a dilution refrigerator Oxford Instruments Triton DR200 at temperature of 35 mK.A sample (Fig. 2c) was placed in a Copper sample holder (Fig. 2b) and bonded by Aluminum wires with a diameter of 25 µm.A sample holder was thermally anchored to a mixing chamber plate of a refrigerator.Microwave signals were generated by Vector Network Analyzer Agilent Technologies PNA-X N5242A.An input signal was attenuated for reduction of thermal noises.In order to measure a nonlinearity of grAl we used a Mini-Circuits Bias-Tee ZX85-12G-S+ for injection of a dc bias current from a current source Keithley 6221.

FIG.
2. An experimental setup scheme (a).In a dilution refrigerator a Copper sample holder (b) with a long coplanar line represented by a meander (c) is cooled down.We measure the transmission coefficient with a vector network analyzer (VNA) as a function of the dc current injected to the transmission line by bias-tees.An input signal is attenuated by 30 dB and an output signal is amplified by a commercial wideband amplifier.

IV. RESULTS AND DISCUSSION
In order to demonstrate the nonlinear properties of superconducting films, the phase of the microwave signal is measured as a function of the amplitude of the bias current.In approximation of low currents the nonlinear phase shift is proportional to the total phase length of the transmission line [48]: In the limit of a high kinetic inductance L gl ≪ L kl , knowing a current dependence of a phase shift ∆ϕ(I) we can find a current dependence of a kinetic inductance using Eqs.( 8) and ( 9): where L k (0) ≈ 3.63 H/m is a kinetic inductance per unit length at low current defined by Eq. ( 3), l is a total physical length of a transmission line, C l = 1.25 × 10 −10 F/m is a capacitance per unit length.Figure 3 summarizes our results of the grAl nonlinear properties study.Red dots in Fig. 3a show the dependence of the microwave signal phase shift on the current at a frequency of 7 GHz.At the bias current of 80 µA the phase shift reaches its maximum ∆ϕ = 4 radian with the total line phase length of 284 radians.This results in the nonlinearity ∆ϕ/ϕ of about 1.4 %, which is comparable to values previously obtained on NbTiN films [11,48].In Fig. 3b red dots show the kinetic inductance calculated from the measured phase shift using Eq. ( 10).
An analytical dependence using an approximation (6) is shown in Fig. 3b by solid blue line, here I * ≈ 318 µA.At the same figure by dashed orange line an analytical dependence from Eqs. (4) and ( 5) is shown, where I dep ≈ 144 µA.Measurement results are in a good agreement with both analytical curves.A measured critical current I c ≈ 80 µA, which is 56 % of the theoretical limit.We believe that the measured value of the critical current could be reduced by defects and not optimal structure geometry.At this frequency the largest phase shift is 4 radians, with total line phase length of 284 radians.Experimentally obtained kinetic inductance (red dots) is in a good agreement with both analytical curves: BCS approximation given by Eqs. ( 4) and ( 5) -a dashed orange line, and an approximation of Eq. ( 6)a solid blue line.
For a trivial long dispersionless line, a gain G can be estimated as [11] G ≈ 1 + (∆ϕ) 2 . (11) In our experiments, the phase shift is caused by a dc bias current passing through the transmission line, which also occurs when passing a microwave signal of higher power.Using Eq. ( 11) we estimate a microwave signal amplification in a dispersionless transmission line to be expected of the order of 17 dB at a frequency of 7 GHz.
In conclusion, we experimentally demonstrated the nonlinear dependence of the grAl thin films kinetic inductance on the dc current, which, to the best of the authors' knowledge, has not been previously studied.On a line with a phase length of 284 radians, a phase shift of 4 radians is obtained at a frequency of 7 GHz.This value allows us to expect a significant amplification of the microwave signal, opening up the possibility of creating parametric amplifiers based on grAl.

FIG. 1 .
FIG. 1.(a) Stages of a sample fabrication: a grAl film is deposited by an electron beam evaporator (2) on an oxidized silicon substrate (1).After spin-coating a photoresist (3), a mask is formed by an optical lithography (4).Then by an ion gun etching (5) a final topology is created (6).(b) A calibration curve: an oxygen flow dependence of a grAl resistance per square.A vertical line indicates a sample with resistance per square of 200 Ω, considered in this work.

FIG. 3 .
FIG.3.The microwave signal's phase shift ∆ϕ (a) and the kinetic inductance (b) versus dc current at a frequency of 7 GHz.At this frequency the largest phase shift is 4 radians, with total line phase length of 284 radians.Experimentally obtained kinetic inductance (red dots) is in a good agreement with both analytical curves: BCS approximation given by Eqs.(4) and (5) -a dashed orange line, and an approximation of Eq. (6)a solid blue line.