High mobility SiMOSFETs fabricated in a full 300mm CMOS process

The quality of the semiconductor-barrier interface plays a pivotal role in the demonstration of high quality reproducible quantum dots for quantum information processing. In this work, we have measured SiMOSFET Hall bars on undoped Si substrates in order to investigate the quality of the devices fabricated in a full CMOS process. We report a record mobility of 17'500 cm2/Vs with a sub-10 nm oxide thickness indicating a high quality interface, suitable for future qubit applications. We also study the influence of gate materials on the mobilities and discuss the underlying mechanisms, giving insight into further material optimization for large scale quantum processors.

The spin of an electron in Silicon has been considered as one of the most promising candidates for largescale quantum computers, due to its long coherence time, compactness, potential to operate at relatively high temperatures, and compatibility with CMOS technology for upscaling [1][2][3]. High fidelity single and two qubit operations have been demonstrated in academic lab-based devices [4,5]. The scale up process to multi-qubit arrays will however entail high quality qubits in large numbers, necessitating a transition to industrial fabrication techniques [6].
While spin qubit-devices can be fabricated using semiconductor manufacturing techniques, a detailed and careful analysis on the impact of different fabrication process steps on the performance of spin qubits is crucial. Hallbar mobility is a widely used metric to characterize the MOS interface quality, which can provide valuable feedback to fabrication process optimization. In this work, we perform transport measurements on SiMOSFETs and investigate the quantum Hall effect. The gate stack for SiMOSFET is the same as that for qubit integration [7]. We extract key characteristics of the MOS gate stack such as critical density and mobility as well as transport and quantum lifetimes. Comparing the obtained results to a transport model, we identify and quantify the leading scattering terms. Additionally, we explore possible mobility limiting factors when changing the top gate material from PolySi to titanium nitride (TiN).
The devices measured in this work are undoped inversion gated Hall bars fabricated in a state-of-the art 300 mm fab line [7]. The starting substrate consists of a (100) silicon wafer with a background doping <5 × 10 15 cm −3 . The Hall bar top gates are defined by electron beam lithography and subsequent dry etching process, and the fanout metal pins are defined by deep UV lithography. After depositing the gate material, the whole wafer is covered with a passivation layer. In order to measure the electronic properties of the investigated samples, a positive voltage has been applied to the top gate to form a 2-dimensional electron gas (2DEG) at the interface of Si and SiO 2 . Throughout this work, we perform four-wire electronic transport measurements using standard lockin techniques. Both the longitudinal and transverse Hall voltages were measured simultaneously with differential voltage amplifiers. Additionally, the current through the sample was measured using an IV transimpeadance amplifier. The measurements were carried out in a dilution refrigerator with a base temperature of 20 mK with a 9 T magnet. Higher temperatures were achieved by heating the mixing chamber, resulting in electron temperatures between 100 mK and 4.2 K. Fig.1 shows the density and mobility study performed on a sample with a PolySi top gate at three different temperatures. First, the density of charge carriers in the 2DEG was extracted from the transverse Hall effect in the linear regime far from quantum oscillations, particularly at the lowest temperature with magnetic field between ±1 T. The density of charge carriers in the 2DEG has an overall linear dependence on the applied top gate voltage V G as expected from a plate-capacitor model. Assuming the relative dielectric constant r (SiO 2 ) = 3.9, we extract the thickness of the oxide and obtain a value of 8.1 ± 0.1 nm. This in good agreement with the oxide thickness of 8 nm defined during the fabrication process. As one can see in the inset of Fig.1a, the extracted density deviates from a linear regime at very low V G near the threshold voltage, particularly at the lowest temperature. A temperature dependent dielectric constant would result in a temperature dependent slope of the n S regardless of the V G [8], which is not observed in this study. Partial localization of charge carriers in the disorder potential at low sheet density, on the other hand, could explain this [9]. Here, the disorder potential is mainly due to bulk charge scattering and scattering from oxide charges. At these rather low densities close to the 2D Metal to Insulator Transition (MIT), partial localization would effectively reduce the sheet density of mobile carriers below a simple capacitor model, as observed, and thus slightly enhance the mobility, which is inversely proportional to the density.
The extracted densities, in combination with the sheet conductance from the longitudinal voltage measurements, allow for the calculation of the mobility µ at various temperatures. The resulting curves are shown in Fig.1b. We measured a maximum mobility of 17.5 × 10 3 cm 2 /Vs, among the highest mobility reported for SiMOSFET devices with sub 10 nm oxide thickness [10][11][12]. We fit the mobility curve with the Kruithof-Klapwijk-Bakker (KKB) model (Ref. [13]) which includes remote impurity scattering and surface roughness terms [13,14]. Solving the model self-consistently for three separate temperatures (0.5, 0.9 and 4.2 K) yields the charged impurity center density N C = 1.53 ± 0.03 × 10 10 cm −2 , the rms height of the amplitude of the interface roughness ∆ = 5.9 ± 0.51Å and the correlation length L = 1.8 ± 0.35 nm. These numbers indicate a smooth and clean interface which is critical for future qubit implementations.
As shown in Fig.1b, there is a slight deviation of the maximum mobility from the theory curve at 0.9 K, and the theoretical curve does not agree with the reported data for 0.1 K in the low density regime. This could be explained by the assumption of a linear relationship between mobility and temperature in the KKB theory. According to our fitting, this assumption seems to be only valid above 1 K and can become a more complex function of temperature below 1 K if the low density regime becomes non-linear. This is due to the fact that the slightly reduced density (compared to the capacitor model) results in a slightly higher mobility, explaining the deviation of our data from the theory. Furthermore, for sheet densities n S close to the depletion charge density n dep , the calculation of the finite extent of the wave function in the KKB theory can become unreliable because of the assumed approximations. The theoretically predicted curve can be recovered by using a linearly extrapolated density during the calculation of the mobility. Further details are discussed in Appendix A.
Another key parameter is the percolation density, which indicates the minimum density required to form a conducting channel. Generally, the percolation density describes the disorder in the channel at low densities, where quantum devices typically operate. In Fig.1d, we employed a MIT model in order to fit the densitydependent conductivity σ xx ∼ (n − n p ) p [15], where n p and p are the percolation density and universal exponent, respectively. However, this model is only valid in the low density regime and deviates above n =8 × 10 11 cm −2 . We fixed the exponent p = 1.31, as expected in a twodimensional system, and allowed n p to vary [15]. Our fits yield n p = 1.86 × 10 11 cm −2 at T = 4.2 K. As seen in the inset of 1d the extracted n p decreases with decreasing temperature, which is in agreement with previous studies [15], and is below 1 × 10 11 cm −2 at the lowest temperature. Such a low percolation density is among the lowest reported for similar devices [12,16]. The small percola-tion density and high peak mobility indicate low interface disorder of the MOS gate stack, which is crucial for large scale quantum dot arrays.
Additionally, we investigate the temperature dependence of mobility at various top gate voltages. For sheet densities less than 2 × 10 12 cm −2 the mobility is a function of temperature, with a linear dependence for temperatures above 1 K (see Fig.1c). The temperature dependence of the mobility is only visible for low electron densities because mobilities at higher densities are limited by surface roughness, which is not temperature dependent. Additionally, the peak mobility shifts to higher densities at higher temperatures due to the activation of background charges limiting the mobility at high temperatures.
In order to extract more information about the transport properties and the dominating scattering mechanism, we investigated the transverse resistance R xy and longitudinal resistivity ρ xx in the quantum Hall regime. The magnetic field B was swept up to 8.5 T at a fixed sheet density of 7.1 × 10 11 cm −2 near the peak mobility. As shown in Fig.2a, we observe Shubnikov-De Haas (SdH) oscillations starting at 0.8 T and begin to resolve spin splitting at 3.2 T. The ρ xx shows well developed zeros for magnetic fields larger than 3.8 T. In combination with the single frequency of the SdH oscillations extracted below 2.5 T and the absence of beatings strongly indicate single subband occupancy. We further observe a four-fold degeneracy in the filling factor which changes to a two-fold degeneracy as we increased the magnetic field (Fig.2b). The remaining two-fold degeneracy is due to the two lowest lying valleys which are split off from the remaining four valleys due to the z-confinement in Silicon (100) [17]. A linear fit of the filling factor ( Fig.2b) with respect to inverse magnetic field yields a sheet density 7.1×10 11 cm −2 in agreement with the extracted Hall density, further confirming clean single subband transport in the investigated samples.
Upon further investigations into the SdH oscillation amplitudes below 2.5 T, we were able to extract an effective mass m * = 0.21 ± 0.02 · m 0 from the temperature dependent data, which corresponds to the transverse mass expected for Silicon samples with similar background doping [18]. Further details are discussed in Appendix B. From fitting the SdH envelope as a function of temperature and density, we can extract the quantum lifetime τ q . By combining τ q , which includes all possible scattering events, and the transport lifetime τ t = e * µ/m * , which includes only the large angle scattering, we extract the so-called Dingle ratio τ t /τ q . The Dingle ratio thus characterizes the dominant scattering mechanism ranging from a ratio of 1 for predominant large angle scattering to values >> 1 for dominant small angle scattering. Fig.2c shows an increase of the Dingle ratio from ≈ 1 to ≈ 2 with increasing temperature, indicating a shift in the type of scattering dominating at different tempera-ture ranges. We can see an inverted trend for the Dingle ratio with respect to density. For a fixed temperature of 0.5 K, we observe a decrease from a ratio of ≈ 1.5 to ≈ 1 by increasing the density, again showing a slight shift towards large angle scattering. This is true for all measured temperatures and densities. In all, the observed Dingle ratio around 1 to 2 shows that short range / large angle scattering is the dominant scattering mechanism in this system, where the scattering centers are near or in the channel itself [18].
Furthermore, we investigate the effect of two different gate metals on the mobility and percolation density. We chose PolySi and titanium nitride (TiN) as top gate materials since both are widely applied gating materials for integrated CMOS device fabrication. While PolySi has generally larger mobilities, which indicates a small disorder potential and therefore a low dot-to-dot variation, TiN exhibits low resistivity and is superconducting at low temperatures, which is desirable for microwave qubit control and readout. The density dependent mobilities at 4.2 K, for both TiN and PolySi are shown in Fig.3a. There is a decrease by a factor of three in the maximum mobility when changing the material from PolySi to TiN while all other parameters in the fabrication remain the same. To investigate the underlying mechanism for the mobility difference, we perform high magnetic field Hall measurements on the TiN devices as well. Unlike PolySi which shows vanishing ρ xx in Fig.2a, ρ xx of TiN does not have well developed zeros. This can be attributed to the lower mobility and could be due to multiple subband transport [19] (see Appendix C for further details). At high electron concentrations and in the presence of compressive stress, multiple subbands can be occupied and thus introduce inter-subband scattering at low temperatures.
To explain the inter-subband scattering, we firstly consider the effect of strain. On silicon (100) surfaces, materials used as top gates with certain coefficient of thermal expansion (CTE) introduce compressive strain [20]. This leads to a decrease in the energy splitting of two subsequent subbands [18,21,22]. Strain simulations, described in Appendix D, indicate a large difference in strain between TiN and PolySi. However, the strain could only decrease the TiN mobility by roughly 10 percent rather than the three fold decrease observed [22]. This difference could be explained by the oxygen scavenging of TiN on SiO 2 [23]. When depositing TiN on a SiO 2 surface, SiO 2 gets reduced to SiO x and the free oxygen combines with TiN to form TiO 2 and titanium oxynitride (TiN x O y ). This process creates charged oxygen vacancies in the SiO 2 layer and because the oxide is very thin in these samples, the charge carriers can then scatter at these defects thus reducing the overall mobility. Indeed, fits of the KKB model to the TiN data indicate a significantly higher N C than PolySi consistent with the scavenging hypothesis. Those two effects com-bined therefore could explain the difference in mobilities between the two gate metals.
In summary, we have characterized SiMOSFET samples using Hall effect measurements. We found a record mobility of 17.5 × 10 3 cm 2 /Vs which is the highest reported mobility to date for SiMOSFET samples with sub-10 nm oxide thickness. Further, we have shown that the dominant mobility limiting factor are charges near the conducting channel, which can be improved by further reducing the impurity levels in the substrate. For example, an epitaxial layer can reach a background doping level of roughly 1×10 12 cm −3 , which is three orders of magnitude lower than the samples measured in this study. Additionally, we show that the choice of top gate material directly affects the mobility and percolation density. Factors such as strain and oxygen scavenging could explain the strong influence on both parameters. This demonstrates that the choice of top gate material is of paramount importance. The MOS interface has been widely identified as the source for charge noise and disorder sites, which limits the qubit performance and upscaling. To further improve the qubit device fabrication, more systematical investigations are needed. The Hallbar mobility and the analysis presented in this study could provide valuable insight into the gate stack over a large device area and direct further fabrication optimizations, which could pave the way for a large scale spin qubit processor integration with a full CMOS process.

DATA AVAILABILITY STATEMENT
The data supporting the plots of this paper are available at the Zenodo repository at https://doi.org/10. 5281/zenodo.xxxxx.

A. Mobility as a function of linear density and MIT fitting range
During the calculation of the mobility at the lowest temperature, we observed a peak mobility higher than the KKB model (see Fig.1b). We associated this deviation to the non-linear density dependence in the small top gate voltage regime at the lowest temperature as shown in the inset of Fig.1a. This assumption is further supported because the theoretical curve can be recovered using a linear extrapolated density when calculating the mobility (see Fig.4). This behaviour can be explained within the KKB theory, where the depletion charge density n dep is implicitly assumed to be significantly smaller than the sheet density n S especially in the critical regime of maximum mobility [13]. If the density becomes non-linear, this violation of the linear density assumption has a significant impact when calculating the finite extent of the wave function b within the Coulomb scattering term |U C (q)| 2 and the dielectric response function ε. Further research in this regime is needed in order to capture the behaviour found in this work, but it exceeds the scope of this work. Nevertheless, we emphasize the mobility of 17.5 × 10 3 cm 2 /Vs is the properly measured mobility, and the slightly lower peak mobility of 16.5 × 10 3 cm 2 /Vs in Fig.4 is a modeled mobility assuming a linear density model. The previous mobility record was ≈10 × 10 3 cm 2 /Vs for a comparable sample although the oxide was slightly thicker in that work [12].
The metal to insulator transition model (MIT) we used in this work is only valid at very low densities (below 1 × 10 12 cm −2 and therefore is rather sensitive to the fitting range. In order to find the appropriate fitting range, we varied both the minimum and maximum value independently. As shown in Fig.4, we determined that fitting the data between 2.7 × 10 11 cm −2 and 6.0 × 10 11 cm −2 yields the most reliable outcome because variations in the fitting range do not significantly change the result.

B. Extraction of effective mass and quantum lifetime
In order validate our evaluation approach and our experiments, we extract the effective mass from the temperature dependence of the SdH oscillations. For that, we subtract a polynominal background from the ρ xx data and use the following formula as described in [24,25] ∆ρ xx (T ) with χ(T ) = 2π 2 k B T /hω c and ω c = eB/m * . Here, T is the temperature, T 0 the lowest temperature, k B is the Boltzmann constant,h the reduced Plank constant, e the electric charge, B the magnetic field and m * the effective mass. Using the effective mass as the only free fit parameter, we extracted m * = 0.21 ± 0.02 · m 0 , which is close to the expected value of m * = 0.19 · m e . Further, the envelope of the SdH oscillations allows for the extraction of the quantum lifetime τ q . For that, we use the following formula Therefore, by analysing the slope of the magnetic field dependent logarithmic term we can extract τ q for all measured densities and temperatures. In combination with the transport lifetime τ t we can then further extract the dingle ratio as described in the main text.

C. TiN data at high fields
In comparison to the reported data in Fig.2a, TiN does not show well developed zeroes in the longitudinal resistivity ρ xx at high magnetic fields at a temperature of 500 mK (see Fig.7), as expected for a much lower mobility. Further, the transversal Hall resistance R xy shows oscillating features instead of the expected Hall plateaus. This could be due to multiple subbands contributing to the overall transport or with the increased density of background charges due to the above mentioned oxygen scavenging.

D. Strain simulations
We quantify the stain induced in the 2DEG due to the different thermal expansion coefficients of both gate materials. The strain was simulated using a commercially available software [26]. The silicon was simulated isotropically. The temperature dependence of the CTE for PolySi up to room temperature is assumed to be linear because there is no reported data above 200 K [27]. Further, we neglect the intrinsic strain induced during the fabrication process, which is beyond the scope of this work. Thermal contraction is considered from room temperature down to the cryogenic measurement temperatures. All simulations were performed in three dimensions and the components of strain x , y , z are calculated 5 nm below the oxide due to the finite extension of the wavefunction. As shown in Fig.7b and   Linear relationship between filling factor ν (blue scatter) and inverse magnetic field. The solid line is a fit from which the sheet density is calculated. (c) and (d) Extracted dingle ratio τt/τq as a function of temperature T and sheet density nS, respectively. The dingle ratio is slightly temperature and density dependent and is ≈ 1, which indicates that large angle scattering is most dominant scattering mechanism in this sample. FIG. 4. Mobility as a function of linear density and MIT fitting range (a) When using a linearly extrapolated (instead of the carefully extracted) density in order to calculate the mobility of our sample (blue), the peak mobility becomes slightly smaller than previously demonstrated in Fig.1b. Additionally, we can fit the KKB theory (black) selfconsistently with the other temperature data highlighting the importance of the density values. (b) By varying ∆ns for extracting the percolation density np within the metal to insulator transition (MIT) model, we determine the optimal fitting range. We chose ∆ns =3.2 × 10 11 cm −2 (indicated by the red box) because the extracted np does not change significantly while changing ∆ns.