We fabricated a resonator structure on a high resistivity (> 5k ohms. cm) silicon wafer (see experimental section for more details), as shown in Fig. 1A. We characterized the structural properties and chemical composition of the film before and after the device fabrication: a crystalline orientation using an X-ray diffraction (XRD) pattern, surface morphology by atomic force microscopic (AFM) image, and the chemical composition investigated with X-ray photoelectron spectroscopy (XPS). Figure S1A shows XRD patterns of the as-deposited Nb, UV-ozone exposed (Nb-O), and oxygen plasma exposed (Nb-OP) samples. The results reveal two peaks at 38.47◦ and 82.45◦, corresponding to the [110] and [220] crystal orientations, indicating that the Nb film contains a plane of [110] structure as the dominant orientation compared to the [220] plane structure. The inset Figure S1A displays an expanded view of the [110] pattern, indicating that the fabrication and plasma cleaning processes do not affect the crystal orientation of the native Nb film. However, there is a small change in the crystallite size, which increases from 21.54 nm to 22.23 nm after the fabrication, indicating that the lithography followed by the dry etching process does not change the lattice strain significantly. The oxygen plasma-exposed samples exhibited a crystallite size of ∼21.82 nm. These results further confirm that the structural characteristics of Nb film are retained after the oxygen plasma exposure.
The surface morphology of the Nb thin film is characterized using tapping mode AFM. Figure 1B&C shows the AFM image of ozone and oxygen plasma exposed. The elongated crystal structure obtained (Fig. 1B) resembles previously reported highly crystalline niobium resonator structures fabricated on sapphire and silicon surfaces.9 The Nb thin film with ozone exposure has an RMS value of 0.98 nm over a 1×1 µm² area. After 6 minutes of oxygen plasma exposure, the elongated, densely packed crystal structure disappeared (Fig. 1C), and the Nb surface exhibited an atomically smooth surface with an RMS value of 0.31 nm over the same area. This observation is similar to what was previously reported for the Nb film exposed to an argon milling process,35 which led to smoothening the Nb film surface by increasing valley and peak areas by a factor of 5.6, simultaneously, argon milling causes damage to the surface of niobium.35 Here, the Oxygen plasma exposure partially etches the Nb film surface, making it smoother and significantly reducing the average pinhole depth from 4 nm to 1.5 nm (see Fig. 1D). C.T. Earnest et al.36 demonstrated that surface roughness, specifically horizontal peak-to-valley height, has been linked to the increased loss due to the enhanced participation ratio of the electric fields at the surface while other parameters remaining same. This indicates that it is essential to reduce the overall horizontal peak-to-valley height of the surface, including pinholes, to decrease the electric fields at the surface. We studied the time-dependent oxygen exposure on a Nb film to understand the structural changes (e.g. RMS roughness) of the surface without affecting the internal structure of the film. Figure S1B illustrates changes in the RMS value relative to exposure time, with the RMS value reaching a minimum at approximately 6 minutes of exposure. Upon further exposure time, the roughness value increased, and pinholes were formed on the surface of the Nb film (see Figure S1B). Thus, we confirm that controlled oxygen plasma exposure only alters the surface morphology, resulting in minimized RMS roughness (see Figure S1C&D) and a reduced pinhole on the surface without impacting the Nb film’s crystal structure.
To better understand how oxygen plasma exposure affects the Nb resonator fabricated on a silicon substrate, we investigated the surface chemical composition of the Nb resonator with ozone and oxygen plasma treated chip using the XPS surface characterization. Figure 2A displays a C 1s spectrum with different post-fabrication treatments, clearly showing similar peak positions and shapes. This indicates that any carbon contamination is only attributed to environmental factors and not leftover polymer contamination from the fabrication process. The Si 2p spectrum (Fig. 2B) demonstrates an increase in SiO2 composition compared to ozone treatment alone, indicating that oxygen plasma increases the SiO2 concentration at the substrate-air interface. In Fig. 2C, the O 1s spectrum shows an additional peak at 529.85 eV for the Nb-O sample, corresponding to Nb metal hydroxide. This indicates a significant change in the surface metallic oxide composition due to oxygen plasma treatment. Furthermore, Fig. 2D depicts a decrease in the intensity of the low binding energy peak corresponding to the metallic Nb at 201 to 204 eV, indicating that oxygen plasma exposure increases the oxide thickness of the fabricated Nb film. These observations collectively suggest that while oxygen plasma treatment changes the chemical composition of metal surface oxide, it does not have a notable impact on the oxide composition of the bulk Nb film.
A cross-sectional TEM image provides a clear view of the effect of oxygen plasma exposure on the Nb film surface native oxide and bulk interfaces. Figures 3A&B show cross-section images of Nb surfaces treated with ozone and oxygen plasma. These images demonstrate that the native oxide thickness increases by 2.0 nm after exposure to oxygen plasma for 6 minutes, consistent with the XPS results presented earlier. The ozone-treated Nb film showed an interplanar spacing of 2.22 Å, which is smaller than the standard out-of-plane orientation of Nb crystal [110] interplanar spacing of 2.34 Å. This suggests that there is an internal compressed stress present on the patterned Nb-O resonator structure. However, the Nb film sample exposed to 6 minutes of oxygen plasma had an interplanar spacing reaching the standard value of 2.35 Å, indicating that the oxygen plasma treatment slightly released the compressed stress in the Nb film. To further understand the effect of oxygen plasma treatment, we investigated the low-temperature superconducting Tc, as shown in Fig. 4A. The Tc is extracted from the typical phase transition signature in the resistance vs temperature curve. The measured Tc is 9.02 K for Nb-O and 8.98 K for Nb-OP, close to the bulk Tc of 9.35 K,37 and similar to the previously reported high-quality Nb film deposited using the PVD method.38 The observed close and similar Tc for both patterned and exposed ozone and oxygen plasma further confirms that quasiparticle formations are expected to be low, as mentioned in the introduction. We observe a slight decrease in the Tc after exposure to oxygen plasma treatment. This decrease is attributed to a 2.0 nm increase in surface oxide thickness. Additionally, the relative resistance ratio of the Nb-OP films (4.83) slightly increases compared to the Nb-O film (4.67). This confirms that as the film surface roughness decreases, the low-temperature resistance decreases due to less surface scattering. These results confirm that the Tc value depends on the material properties, precisely the difference in oxide thickness. In contrast, the relative resistance ratio may depend on the surface quality of the film.
Figure 4B illustrates the microwave transmission response of the Nb-O sample at 1.2 K, revealing six distinct dips corresponding to the designed resonators at a photon number 109 on the measured chip. However, some chips exhibit inconsistencies in resonator presence, which we attribute to environmental factors such as dust particles adsorbed near the resonator conductor lines during assembly or wire-bonding processes. The yield of working resonators varied for each chip but always exceeded > 90%, the previously reported yield of ≥ 50% by Drimmer et al.32 Before the assembly process, it is crucial to ensure the absence of dust particles and fabrication defects near or around the resonator conductor lines. Through careful investigation and optical inspection of the fabricated chip, we aim to either improve or avoid missing resonators on a single chip. Figure 4C presents a high-resolution view of a single sampled resonator scan and a fitted line using a previously reported fitting procedure39 (see the detailed description in Annexure 1&2 and corresponding Equation S1) to extract each resonator's Qi. To comprehend the variation in quality factor across a single chip, we measured three randomly chosen chips sourced from separate wafers and diverse locations on those wafers. The plot of all the Qi’s, combined with their respective standard deviation obtained from the three independent chips at each specific frequency, reveals a change in the Qi as a function of the frequency. To understand the high-frequency drive couples to the defects on the thin film surface, which are area-dependent, we plotted the Qi as a function of the resonator interface surface area (Fig. 4D). In particular, we are focused on the Nb vacuum/air interface, henceforth called the top interface, while keeping all other interfaces unchanged. The resonators in this study are designed for frequencies between 4 to 6 GHz, corresponding to length variation from 16 to 9 mm, respectively. In the low-frequency (correspondingly larger top-interface area), the smoother samples (Nb-OP) always show higher Qi than their rougher (Nb-O) counterparts. As the top-metal/air interface area becomes smaller (higher frequency), the Qis are independent of the smoothness of the interface. This could be explained by considering the total loss, which is the inverse of the internal quality factor. This is contributed by three components of resistance arising from TLS (RTLS), QP (RQP), and a residual (Rres) in an additive manner (see detailed discussion in Annexure 3). The TLS is temperature dependent and only becomes significant for temperatures in the mK range, and at low power, the RQP is substantial at our operating temperature (1.2 K); however, it is related to the bulk property and hence volume. On the other hand, the Rres originates from the surface resistance originating from different sources, including fundamental residual surface resistivity from boundary states as given by BCS theory, surface defects/grain boundaries, and, more recently, postulated edge boundary states.40 The first one is a fundamental limit, but the other two are technical in origin, and it is not easy to distinguish their contribution as they impact the Rres in the same way. The samples, Nb-OP and Nb-O, show different roughness, leading to significantly different residual surface resistances and, consequently, different internal quality factors. Furthermore, for the rough Nb-O samples, we show an area dependence as expected from theory and also observed in our experiment.