Performance of high impedance resonators in dirty dielectric environments

High-impedance resonators are a promising contender for realizing long-distance entangling gates between spin qubits. Often, the fabrication of spin qubits relies on the use of gate dielectrics which are detrimental to the quality of the resonator. Here, we investigate loss mechanisms of high-impedance NbTiN resonators in the vicinity of thermally grown SiO2 and Al2O3 fabricated by atomic layer deposition. We benchmark the resonator performance in elevated magnetic fields and at elevated temperatures and find that the internal quality factors are limited by the coupling between the resonator and two-level systems of the employed oxides. Nonetheless, the internal quality factors of high-impedance resonators exceed 103 in all investigated oxide configurations which implies that the dielectric configuration would not limit the performance of resonators integrated in a spin-qubit device. Because these oxides are commonly used for spin qubit device fabrication, our results allow for straightforward integration of high-impedance resonators into spin-based quantum processors. Hence, these experiments pave the way for large-scale, spin-based quantum computers.


Introduction
Understanding the origin of dielectric losses is crucial when exploiting superconducting resonators for quantum information science.But studies of dielectric losses in various dielectrics have so far been based on low-impedance resonators [1][2][3][4][5][6][7][8][9].However, a large resonator impedance is desirable, in-particular in the context of spin-qubits, as the coupling to a weak electric dipole moment scales with the square root of the impedance [10].Previous studies showed that high-impedance resonators in a clean electrostatic environment typically reach quality factors on the order of 10 5 [11][12][13][14] in the absence of dielectrics.
Incorporating high-impedance resonators with dielectrics would allow one to use established fabrication recipes that employ gate dielectrics aiming on electrostatic stability of quantum dot devices.The fact that the quality factor of a high-impedance resonator, incorporated in a real device, is not limited by dielectric losses therefore raises the following question: To what extent can dielectrics be employed when fabricating resonator devices without limiting the quality factor to 10 3 ?
To answer this question, here we investigate high-impedance NbTiN resonators in a variety of dielectric configurations that are commonly used when fabricating double-quantum dots based on semiconductor nanowires [29][30][31][32][33].We demonstrate that despite a reduction of the resonator quality factors due to the additional dielectrics, their quality is sufficient even in a dielectric configuration ideal for nanowire integration.
We describe the crucial parameters during sputtering of the material and investigate the dominant resonator loss mechanism.A crucial criterion for resonators aiming on coupling to spin qubits is their magnetic-field resilience which can be achieved by employing disordered superconductors with a high critical magnetic field [12,13,[34][35][36].Recently the community has started to operate spin qubits at elevated temperatures [37][38][39].Using these arguments as a motivation, we characterize the resonator performance in large magnetic fields up to 5 T and at elevated temperatures.

Experimental setup
We fabricated a total of 16 NbTiN coplanar waveguide resonators with an impedance of Z = √ L/C ∼ 2 k distributed on four different chips with differing dielectric configurations.The sample preparation is described in Sect.3.Each chip hosts a feedline with four notch-type, half-wave resonators as shown in Fig. 1a).From left to right, the different wrapping of the feedline results in coupling quality factors between Q c ≈ 10 3 and Q c ≈ 10 5 (see Table 3 in the Appendix).This large spread of Q c enables us to investigate the film properties accurately for a large range of internal quality factors Q i .Resonance frequencies are in the range between f r ∼ 4.2 GHz and f r ∼ 5.6 GHz.
Figure 1b) shows the four different dielectric configurations in colored, dashed boxes.For illustrative purposes, the corresponding dielectric configuration of a nanowire device is shown as well and we will explain the advantages and disadvantages of each configura-  tion in the context of nanowire integration below.Case I -NbTiN on Si -is the most ideal configuration for the resonator.But fabricating bottom-gate based devices [33,40] directly on top of intrinsic silicon comes with the problem of gate leakage, because of the small gate pitch.Moreover, on the surface of the intrinsic silicon, a native silicon oxide forms under ambient conditions which might result in a poor electrostatic device stability.Therefore, fabrication of nanowire devices on top of intrinsic silicon involves sandwiching the bottom gates with two oxide layers grown by atomic-layer deposition (ALD).To maintain a pure dielectric environment of the resoantor, the oxide has to be wet-etched or the ALD-layers have to be deposited locally by a lift-off process [41].Wet-etching of the oxide might lead to unwanted surface-chemistry on the surface of the NbTiN [42].And, since ALD growth is a conformal processes, the lift-off process might result in irregular, rough edges around the desired structures that may protrude significantly out of plane with respect to the substrate (red crosses in schematic).These edges in turn, may lead to step coverage issues on subsequent metal layers.
The local deposition of oxides for nanowire device integration is alleviated if the whole chip, including the resonator, can be covered with an ALD-grown oxide.We investigate this in case II -Al 2 O 3 on NbTiN on Si.
For device integration, it is desirable to work with electrostatically silent oxides.Therefore, nanowire devices are commonly fabricated on top of thermally grown silicon oxide.We therefore investigate the performance of resonators on top of silicon oxide in case III -NbTiN on SiO 2 on Si.In this case bottom-gate based nanowire devices only require one local oxide deposition step as indicated in the schematic.Additionally, the remaining local oxide deposition is alleviated in case IV -Al 2 O 3 on NbTiN on SiO 2 on Si.
The color codes as introduced in Fig. 1b) are used throughout the rest of the manuscript and denote the dielectric configuration and Table 1 shows the calculated participation ratios of each dielectric for the four investigated dielectric configurations.

Sample preparation
Since this work aims to investigate resonator losses due to the choice of the dielectric configuration, the intrinsic Q of the resonators must not be limited by the NbTiN film quality.Here we summarize the steps taken to optimize the fabrication of the used films.
1.As a substrate, we select two undoped Si wafers with a resistivity larger than 10 k cm; one with only a layer of native SiO 2 and the other with ∼100 nm of thermally grown SiO 2 .
2. In order to minimize the impurity density at the metal-substrate interface, the wafer with only native oxide undergoes the following etching steps: (i) a Piranha etch to oxidize the top ∼10 nm that may contain contamination, (ii) an HF bath to remove this oxide layer, and (iii) a second Piranha etch, followed by (iv) a second HF bath seconds before loading the wafer into the sputtering chamber.The second wafer hosts ∼100 nm of thermally grown SiO 2 .In order to remove organic residues but keeping the oxide layer intact, we consecutively use ultrasonic cleaning of the wafer in an aqueous solution of tripotassium orthophosphat, 1 distilled water, acetone and isopropanol before loading the wafer into the sputtering chamber.
3. The vacuum quality in the sputtering chamber plays a vital role.We perform Ti presputtering, resulting in a significant reduction of the chamber base pressure.
4. We pre-sputter the NbTi target to remove the top, potentially contaminated or oxidized layer [43].
5. The sputtering rate has to be maximized by choosing an ideal set of sputtering parameters.See Appendix A for details.Because the impingement rates of oxygen and water decreases with increasing growth rates, higher sputtering rates result in a purer film and accordingly lower loss tangents of the resonators. 2  6.We perform sputtering as close as possible to stoichiometry of NbTiN [44].See Appendix A for details.
7. The resonators are dry-etched using argon/chlorine, offering a higher selectivity against silicon etching compared to the more widely used fluorine based etching recipes [43,45].This makes it easier to prevent over-etching.We note that etch-induced losses [45] might limit internal Q factors to ∼ 10 5 , beyond the regime of our interest.
8. After fabrication, each film is characterized in dc measurements by measuring the critical temperature T dc c and the sheet resistance R sq,dc close to T dc c using etched reference structures.This allows us to estimate the sheet kinetic inductance [46][47][48] which we use to design the resonator geometry.The resonance frequency is designed using analytical equations of coplanar waveguide resonators [49] and the coupling quality factor is estimated by simulating the structure using the electromagnetic simulation software Sonnet.

Determining loss due to two-level fluctuators
To benchmark the performance of the resonators in the different dielectric configurations, we probe each notch-type half-wave resonator by measuring the transmission S 21 through the feedline at the base temperature of a dilution refrigerator T base ≈ 30 mK. Figure 1c) shows an exemplary resonance which is described by [50,51] Here ) is the loaded quality factor and describes a small resonance asymmetry due to interference with a standing-wave background [50].
We identify 4 resonances on every investigated chip and assign every measured resonance frequency f r to a physical resonator.Using the physical dimensions of the resonator for calculating its geometric inductance and capacitance [49], we deduce its kinetic inductance from the measured resonance frequency.Thereby, we use the center conductor width which we measure for every resonator by means of scanning electron microscopy after having performed the experiments presented in this work.The so obtained averaged square kinetic inductance L sq,rf k is given in Table 2 where the error bar represents the root variance.Values for individual resonators are shown in Table 3 in the Appendix.We find that the values of L sq,dc k and L sq,rf k are consistent for either film. 1 The used solution has the brand name deconex® 12 BASIC 2% solution.
2 Discussion with Mihai Gabureac.is independently inferred from the 8 measured resonance frequencies of either film, where the error represents the root variance.Q TLS , Q other , n c and β are fit parameters of Eq. ( 3).The paramagnetic impurities leading to the feature in Fig. 4a) are characterized by the minimum internal Q-factor Q s i and feature width s .The values of D, Q TLS , Q other , n c , β, Q s i , and s represent weighted averages over the data sets of 4 resonators of either dielectric configuration with weights proportional to the inverse of the error bar of the fit, resulting in a maximum weight, if Q c ∼ Q i for which the resonance is most pronounced.g s is the Landé g-factor extracted from Fig. 4b).The color code corresponds to Fig. 1b) The differences in R sq,dc and T dc c between the two films such as the large variance of L sq,rf k is attributed to the small film thickness giving rise to a large effect of film thickness inhomogeneities.
In order to quantify the loss due to two-level systems (TLS) residing in the differing dielectric structures, we measure a resonance trace for every resonator in every dielectric configuration and extract Q i as a function of power applied on the feedline by fitting Eq. ( 1) in a circular fit [52].Figure 2 shows the fitted internal quality factor Q i for two resonators of each configuration.We convert the applied power P in on the feedline to an average photon number in the resonator using [13] where Q i , Q c and ω r are extracted from fitting the resonance curve.TLS residing in the oxides close to the resonator give rise to a power dependent dielectric loss which is usually modeled by [5,[53][54][55][56][57][58][59] 1 In the low power limit, and at low temperatures, Q i is approximately given by Q TLS due to TLS.When increasing n ph above a critical value n c , Q i increases with a characteristic scaling β until eventually saturating at Q other .We fit Eq. ( 3) to the data (solid lines in Fig. 2) and extract Q TLS , Q other , n c and β as fit parameters.The weighted average of these fit parameters for each film are specified in Table 3.
We find that Eq. ( 3) fits well to our data in all four dielectric configurations implying that in the limit of low photon numbers, all resonators are limited by their coupling to TLS.However, the quantitative behavior for the different dielectric configurations differs by a lot.Let us first consider the low-photon limit in Fig. 2. At low photon numbers, the internal quality factor is determined by the coupling to TLS, In case I -NbTiN on Si, Q i saturates at the largest value as the number of photons in the resonator approaches zero.This implies a low abundance of TLS at the interface between the intrinsic silicon and the NbTiN.In case II -Al 2 O 3 on NbTiN on Si, Q i saturates at values approximately an order of magnitude lower which we attribute to the larger abundance of TLS stemming from the ALD-grown oxide on top of the metal and on top of the dielectric.For case III -NbTiN on SiO 2 -and for case IV -Al 2 O 3 on NbTiN on SiO 2 , the saturation of Q i in the low-photon limit happens another order of magnitude lower than for case II.We attribute this decrease to the larger participation ratio (compare Table 1 of the layers below the center conductor compared to the ones above it due to the larger dielectric constant of silicon as compared to the vacuum dielectric constant.The larger importance of the oxides below the center conductor is confirmed by the negligible difference of Q TLS in case III and case IV (with additional oxide on top of the resonator).
In all cases, once the average number of photons n ph exceeds a critical value n c , Q i increases, because the TLS are increasingly saturated and no longer open a photon leakage path [5,53,55,60].
In the high power limit in case I and case II, all TLS saturate, and Q i asymptotically approaches Q other which originates from a power independent source of loss.The origin of Q other potentially lies in the interaction with phonons or quasiparticles.In case III and case IV, Q i does not saturate even at photon numbers on the order of 10 7 underlining the importance of losses due to TLS in these cases.
Despite the TLS being the dominant source of loss for these resonators, we highlight that Q TLS well exceeds 10 3 even for the configuration where the resonator is sandwiched between SiO 2 and Al 2 O 3 .This result is a central point of this manuscript as it allows for easier integration of semiconductor nanowires into a resonator architecture maintaining a good resonator quality.Moreover, we stress that Q TLS is larger by almost an order of magnitude when oxides are only grown on top of the metal and not below.The inset is a zoom in onto the peak that is observed at ∼ 0.5 K.In all sub-figures, the color encodes the dielectric configuration according to Fig. 1b) and legend in Fig. 2. The data was acquired at an average photon number of ∼ 10 4

Resonator stability at elevated temperatures and fields
After having determined the quality of the resonators in each dielectric configuration, quantified by Q TLS , we aim on benchmarking the resonator stability at elevated temperatures and magnetic fields in regimes relevant for spin-qubit operation.Figure 3a) shows Q i as a function of temperature for all dielectric configurations where the color code corresponds to the one introduced in Fig. 1b).For all curves, we measure an increase in Q i for increasing temperature peaking at ∼ 0.8 K.We attribute this increase in the quality factor to an increasing saturation of the TLS with temperature.When the temperature exceeds ∼ 1 K, a decline in Q i is observed which is attributed to an increasing quasiparticle population because of the closing of the superconducting gap [12,61].
Simultaneously to measuring the quality factor, we also measure the shift in resonance frequency δf r and plot it in Fig. 3b).We find that δf r peaks at a temperature T 0 ∼ 0.5 K (see inset of figure), with the effect being most pronounced for the resonators fabricated on top of SiO 2 (blue points in Fig. 3).Moreover, the positive frequency shift at increased temperature exceeds the positive frequency shift in the case of TLS saturation due to a large photon population (see Fig. 6 in the Appendix).This effect can be understood by noting that the latter only saturates TLS in a narrow band around the resonance frequency ω r while the elevated temperature saturates TLS in a much wider frequency range [5].The temperature of maximum positive frequency shift corresponds to a frequency of ω thermal /2π = k B T 0 /h ∼ 10 GHz.As ω thermal > ω r , saturation of TLS in this frequency range explains the postitive resonance frequency shift due to the dispersive interaction between resonator and TLS.For larger temperature the resonance frequency starts to decrease due to the closing of the superconducting gap resulting in a larger kinetic inductance, explaining the decrease of δf r .
In order to benchmark the stability of the NbTiN resonators in an external magnetic field, we plot Q i as a function of in-plane field B in Fig. 4a) and as a function of out-ofplane field B ⊥ in Fig. 4b).During the latter measurement, the magnetic field B was applied with an angle of 49 • with respect to the sample plane as indicated by the second longitudinal axis in Fig. 4b).In both cases, we observe a dip in Q i at an absolute field strength B s .The dip is attributed to a resonant interaction with paramagnetic impurities in the substrate [12].Remarkably, the dip is much less pronounced for the resonator fabricated  4).In all sub-figures, the color encodes the dielectric configuration according to Fig. 1b) without any additional oxides, indicating that the paramagnetic impurities mainly reside within the oxides.We quantitatively analyze the influence of the paramagnetic impurities by fitting a lorentzian to the dip of each resonator.Therefrom, we extract the minimum quality factor Q s i , the feature width s and the field B s around which the dip is observed (see Table 3 in the Appendix).Figure 4(b) shows B s versus the resonance frequency of the corresponding resonator ω r .This allows to extract the Landé g-factor g s by fitting a straight line given by B s = ω r /g s μ B , where μ B is the Bohr magnetron.The resulting values of g s for the different dielectric configurations are shown in Table 2.
Besides this dip-like feature, for the in-plane field, no noteworthy decline in Q i is observed up to the very largest applied field strengths of 5 T, confirming a magneticfield resilience for in-plane fields.This is expected because the penetration depth λ ∼ 260 nm [48] is much larger than the thickness of the NbTiN film (∼ 10 nm).Out-of-plane, Q i declines monotonously for increasing field-strengths which we attribute to an increasing quasiparticle density in the film.However, Q i remains larger than 10 3 up to B ⊥ ∼ 100 mT for all dielectric configurations, once again confirming their suitability for spin qubit integration.Finally, Fig. 4c) shows the resonance frequency versus out-of-plane field.The data is well fitted by which is deduced from BCS theory [12,46].Here, the width w of every resonator center conductor is measured by electron beam microscopy.The average width is w = 390 ± 120 nm where the error bar is the root variance.D denotes the diffusion constant which is a fit parameter.
The weighted average of the fitted diffusion constants are given in Table 2 and the values for individual resonators are given in Table 3.We note that the diffusion constants of the two investigated films are similar and slightly lower than the one in Reference [12].

Conclusion
We have investigated superconducting, high-impedance resonators based on NbTiN in four different dielectric configurations.The largest internal quality factor in the lowphoton limit is found for the resonator fabricated on intrinsic silicon.Nevertheless, all other dielectric configurations result in internal quality factors > 10 3 , a value which has proven sufficient for performing key experiments in the context of spin qubits.Moreover, we benchmark the resonator performance at elevated temperatures and magnetic-field strengths.Since the resonators are compatible with existing fabrication protocols, our results allow for straightforward integration of these types of resonators with various kinds of spin qubits defined in semiconductor nanowires.
We acknowledge very fruitful discussions with Mihai Gabureac, Sergey Amitonov and Alessia Pally.Furthermore, we thank Dario Marty for the support in wafer etching in the facilities of the Paul Scherrer Institute.
Figure 5 Sputtering parameters.a) growth rate γ as function of nitrogen flow J N 2 at different back-ground pressures ρ bg and plasma power P. b) dc voltage V dc between plasma source and target as function of J N 2 for two valued of ρ bg and P. A kink in the curve is obersed in both a) and b) at the same values of J N 2 .This kink corresponds to sputtering at a stoichiometric ratio.c) growth rate γ as a function of background pressure ρ bg showing a monotonous decay as lower ρ bg correspond to larger mean-free paths and therefore to a smaller scattering of the sputtered material.d) dc voltage V dc between the plasma source and the target as a function of background pressure ρ bg .When lowering ρ bg < 2 mTorr, V dc increases drastically and the plasma becomes unstable.All traces were measured at a fixed argon flow rate Q Ar = 50 sccm Figure 6 Power dependence.Relative frequency shift δf r = f r -f r (n → 0) as function of average number of photons in the resonator n ph .We do not observe a noteworthy positive frequency shift in comparison with the Fig. 3 in the main text.For n ph 10 5 , an increasingly negative relative frequency shift is observed which we attribute to the finite non-linearity of the resonators For very large drive powers with n ph 10 5 , a negative relative frequency shift is observed in Fig. 6.We attribute this negative shift to the onset of the bifurcation of the resonator due to a finite non-linearity.

Additional data
All extracted data for each of the 16 individual resonators is shown in Table 3.
Table 3 Parameters for all 16 resonators.The shown parameters are the resonance frequency f r , the sheet kinetic inductance L sq,rf k , the coupling Q factor Q c , the diffusion constant D, the low-power internal Q factor Q TLS .n c is the critical photon number and β the scaling parameter from Eq. (3).Q S i is the internal Q factor on resonance with the paramagnetic impurities (compare Fig. 4(a)) and S is the width of this resonance.Values are missing where the data could not be extracted

Figure 1
Figure 1 Device overview.a) Top-view design of the four high-impedance (∼2 k ) half-wave resonators, capacitively coupled to a much wider 50 matched feedline.b) Side-view schematics of four different dielectric configurations that were investigated.The color of the dashed boxes corresponds to their respective configurations throughout the paper.The right part of each schematic illustrates how a nanowire device can be realized, given the dielectric configuration on its left.c) Amplitude and IQ-excursion of exemplary resonance.The orange line is a fit to the data

Figure 2
Figure 2 Power dependence.Internal quality factor Q i as a function of average photon number n ph in the resonator.The solid lines are fits to Eq. (3), assuming two level fluctuators as the dominating leakage mechanism at low photon numbers.For each dielectric configuration, as encoded by the color according to Fig. 1b), data for the two resonators with minimal |Q TLS -Q c | are plotted.Different symbols correspond to different resonators

Figure 3
Figure 3 Temperature dependence of resonator properties.a) Internal quality factor Q i as a function of temperature.b) Relative frequency shift δf r = (f r (T) -f r (0)) as a function of temperature.The inset is a zoom in

Figure 4
Figure 4 Magnetic field dependence.a) Internal quality factor Q i as a function of in-plane field amplitude B .A dip is observed at a field B s which is attributed to resonant paramagnetic impurities.b) B s extracted from data similar to the one shown in a) for different resonators as a function of resonance frequency ω r .The solid lines are fits to B s = ω r /g s μ B from which we extract the Landé g-factor g s given in Table 2. c) Internal quality factor Q i as a function of out-of-plane field.The field B is applied with an angle of 49 • with respect to the substrate and the perpendicular component B ⊥ is indicated.Once again a dip is observed around B = hf r /2μ B , being less pronounced for the resonator fabricated on intrinsic silicon.c) Relative frequency shift δf r = f r (B ⊥ ) -f r (0) as a function of out-of-plane field component B ⊥ .The solid lines are fits to Eq. (4).In all

Table 1
Participation ratios.Fraction of electric field energy stored in each dielectric layer for the four investigated films.Values are obtained from performing a dc finite-element simulation using COMSOL

Table 2
Resonator properties.Extracted parameters for the two films with varying dielectric configurations.The sheet resistance R sq,dc and critical temperature T dc c are obtained from a dc measurement.The rf sheet kinetic inductance L