# Analysis of Local Properties and Performance of Bilayer Epitaxial Graphene Field Effect Transistors on SiC

Dalal Fadil, Wlodek Strupinski, Emiliano Pallecchi, Henri Happy

► To cite this version:

Dalal Fadil, Wlodek Strupinski, Emiliano Pallecchi, Henri Happy. Analysis of Local Properties and Performance of Bilayer Epitaxial Graphene Field Effect Transistors on SiC. Materials, 2024, 17 (14), 3553, 12 p. 10.3390/ma17143553. [hal-04673340](https://hal.science/hal-04673340v1)

HAL Id: hal-04673340

https://hal.science/hal-04673340v1

Submitted on 20 Aug 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Analysis of Local Properties and Performance of Bilayer Epitaxial Graphene Field Effect Transistors on SiC

Dalal Fadil <fnref n="1" />,<fnref n="2" />,* ID, Wlodek Strupinski <fnref n="3" />, Emiliano Pallecchi <fnref n="1" /> and Henri Happy <fnref n="1" />,* ID

**Abstract:** Epitaxial bilayer graphene, grown by chemical vapor deposition on SiC substrates without silicon sublimation, is crucial material for graphene field effect transistors (GFETs). Rigorous characterization methods, such as atomic force microscopy and Raman spectroscopy, confirm the exceptional quality of this graphene. Post-nanofabrication, extensive evaluation of DC and high-frequency properties enable the extraction of critical parameters such as the current gain ($f_{max}$) and cut-off frequency ($f_t$) of hundred transistors. The Raman spectra analysis provides insights into material property, which correlate with Hall mobilities, carrier densities, contact resistance and sheet resistance and highlights graphene's intrinsic properties. The GFETs' performance displays dispersion, as confirmed through the characterization of multiple transistors. Since the Raman analysis shows relatively homogeneous surface, the variation in Hall mobility, carrier densities and contact resistance cross the wafer suggest that the dispersion of GFET transistor's performance could be related to the process of fabrication. Such insights are especially critical in integrated circuits, where consistent transistor performance is vital due to the presence of circuit elements like inductance, capacitance and coplanar waveguides often distributed across the same wafer.

**Keywords:** bilayer graphene; silicon carbide; field effect transistors; nanofabrication; Raman spectra analysis; DC and RF characterizations

## 1. Introduction

Graphene is considered one of the most famous two-dimensional (2D) materials of this century [1–4]. This honeycomb carbon atom lattice exhibits extraordinary electrical [5,6], optical [7], thermal [8] and mechanical [9] properties and attracts huge attention for a wide panel of device applications, especially for radio-frequency transistors on rigid or on flexible substrate [4,10–12]. Since 2004, researchers continue improving the technique of growth, either by exfoliation [13–15] or by chemical vapor deposition (CVD). The most common technique used to grow graphene is using CVD on copper and transferring the carbon monolayer to a host substrate [16,17]. However, efforts remain to be made regarding the reliability and the reproducibility of the quality after transfer due to the cracks, wrinkles and residues, as reported by Smith et al. [18]. Another alternative is to grow graphene by CVD directly on SiC substrate by graphitization or by epitaxial growth without SiC sublimation or by hydrogen intercalation [19–25]. Besides the large-scale fabrication, the main advantage of CVD on SiC is to avoid the transfer from copper and improve the electronic properties of the material by using hydrogen intercalation, as reported by Ciuk et al. [26,27]. One of the advantages of the bilayer graphene on SiC substrate is overcoming the lack of a band gap; other possibilities reported in the literature include introducing defects, doping, strain and chemical bounding to the substrate [28–30]. Furthermore, the bilayer graphene field effect transistors (GFETs) on SiC have demonstrated better performances than monolayer

<figure><img src="image_2.png" /></figure>

### check for updates

**Citation:** Fadil, D.; Strupinski, W.; Pallecchi, E.; Happy, H. Analysis of Local Properties and Performance of Bilayer Epitaxial Graphene Field Effect Transistors on SiC. *Materials* **2024**, *17*, 3553. https://doi.org/10.3390/ma17143553

Academic Editors: Marilena Vivona and Mike Jennings

Received: 5 June 2024

Revised: 8 July 2024

Accepted: 11 July 2024

Published: 18 July 2024

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ University of Lille—IEMN CNRS UMR 8520, Avenue Poincaré, CS 60069, 59652 Villeneuve d'Ascq, France; emiliano.pallecchi@univ-lille.fr</footnote>

<footnote n="2">² Departament d'Enginyeria Electrònica, Universitat Rovira I Virgili, 43007 Tarragona, Spain</footnote>

<footnote n="3">³ Faculty of Physics, Warsaw University of Technology, Koszykowa 75 Str., 00-662 Warsaw, Poland; wlodek.strupinski@pw.edu.pl</footnote>

<footnote>* Correspondence: dalal.fadil@urv.cat (D.F.); henri.happy@univ-lille.fr (H.H.)</footnote>
GFETs [31]. The 100 nm gate length ($L_g$) and 2 × 4 µm channel length ($W_g$) transistor exhibits 60 GHz of extrinsic current gain cut-off frequency $f_{T_extr}$ and 25 GHz of maximum frequency of oscillation $f_{max_extr}$, and by annealing the bilayer graphene in hydrogen, the performance can be enhanced and reach $f_{T_extr}$ = 70 GHz and $f_{max_extr}$ = 120 GHz for $L_g$ = 60 nm and $W_g$ = 2 × 8 µm [32]. To analyze graphene properties, the standard non-destructive techniques used are (i) the atomic force microscopy AFM that give information about surface morphology including defects like holes, wrinkles or grain boundaries (ii) Raman spectroscopy which provides information about the number of layers, doping, strain and disorders in the graphene surface [33,34]. A new specific technique is non-contact Terahertz time domain spectroscopy. This complementary technique gives large scale information about mobility and carrier density [35].

The motivation behind this study goes beyond evaluating the RF performance of a single transistor. It also emphasizes the reproducibility of transistor performance across a large area, investigating whether the growth was uniform and the process reproducible. Understanding the material's properties and analyzing the local properties, mobility, contact resistance and their impact on device performance is crucial for enhancing the reliability and uniformity of transistors in integrated circuits.

In this work, we report the fabrication process of bilayer graphene field effect transistors. Initially, we verify the quality of the graphene samples before fabrication using AFM and Raman spectroscopy. Post-fabrication, we conduct further analysis using Raman spectroscopy, Hall mobility and TLM measurements to assess the local properties and the homogeneity of the samples. To evaluate device performance across the entire wafer, we performed DC and RF characterizations on hundreds of transistors. The RF performance data were then compared with the local properties of the samples, supporting a conclusion about device variability across the wafer. Our findings demonstrate significant performance variation across the wafer. The average $<f_t>$/<f<sub>max</sub>>$ ranged from 3/0.5 GHz in cell 4 to 16/11 GHz in cell 7 for devices with $L_g$ = 200 nm and $W_g$ = 2 × 30 µm. Notably, the devices in cell 7 and 8, which exhibited the best performance, also had the lowest contact resistance. This work highlights the important role of the fabrication process in the RF performance variability of graphene FETs fabricated from high-quality epitaxial bilayer graphene on a SiC substrate.

## 2. Materials and Methods

Epitaxial bilayer graphene was synthesized by chemical vapor deposition on a 500 µm thick, high-resistivity 6H-SiC (0001) substrate using a commercial horizontal CVD hot wall Aixtron VP508 reactor (Aixtron, Herzogenrath, Germany) equipped with an RF generator for heating. Prior to growth, in situ etching of the SiC surface was performed under a hydrogen atmosphere at 1600 °C and a chamber pressure of 100 mbar. The carbon films were deposited using propane as the carbon precursor. Our method employs high-temperature and low-argon-pressure CVD, creating laminar argon flow dynamics to protect the SiC substrate from silicon sublimation and facilitate propane mass transport, thereby enabling graphene epitaxy, as detailed in references [22,25,26]. The growth process was followed by in situ hydrogen intercalation at 1000 °C in a 900 mbar Ar atmosphere. The optimization of growth parameters aimed at achieving uniform bilayer graphene. Before fabrication, initial carrier density and the mobility was provided by Ciuk et al. to be around +8.3 × 10¹² cm⁻² and 850 cm²·V⁻¹s⁻¹, respectively. In contrast to the methodology described by P. Wehrfritz et al. [36], which utilized a similar SiC substrate, our process involves unique conditions including high-temperature in situ hydrogen etching and controlled low-argon pressure to prevent silicon sublimation. These optimizations contribute to the superior uniformity and controlled thickness of our bilayer graphene compared to prior studies. Notably, our approach introduces an innovative in situ hydrogen intercalation step at 1000 C, enhancing graphene quality, which was not emphasized in previous literature [36].
The AFM Bruken Icon model was used on mode tapping to determine the properties, the material and image of the surface. Image analysis was performed with WSxM5.0 Develop 8.3, a free software. High-resolution electron beam lithography VISTEC EBPG5000Plus was used for device fabrication. To analyze the local properties of the material, we used a HORIBA Jobin-Yvon lab system for Raman spectroscopy at a laser wavelength of 473 nm, using a 1 µm laser spot size and filters to deliver power less than 0.1 mW and ×100 objective lens to measure the different positions of the sample. HL5500PC was used to carry out Hall measurement. To measure the DC and RF performance of hundreds of devices, we used a standard probe station with Microtech's probes, Semiconductor Analyzer HP4155A, the Vector Network Analyzer HP4155A and the vector Network Analyzer Rohde & Schwarz ZVA67.

A. Before fabrication

Primary characterization was performed before fabrication. Figure 1a represents atomic force microscopy (AFM) images of a 60 × 60 µm² surface area of graphene observed at room temperature. The AFM images reveal well-oriented, parallel atomic steps with SiC terraces approximately 10 µm wide, separated by steps estimated to be a few nanometers in height. This phenomenon, known as step bunching, occurs during the growth-preceding in situ hydrogen etching of the SiC surface, as described in reference [27]. Within these terraces, Figure 1b showcases two distinct surface morphologies. The graphene roughness within individual terraces has a root mean square (RMS) value of 0.273 nm (Figure 1b, left) and 0.2406 nm (Figure 1b, right), significantly smoother than the 1 nm roughness reported in previous studies [36]. Figure 1c shows the SEM image of epitaxial graphene on the SiC substrate. We observed the graphene surface nucleation on the SiC steps. An example of a Raman spectroscopy of graphene on SiC and after extracting SiC peaks is presented in Figure 1d. We observed small intensity of the D peak compared to the G and 2D peaks, indicating a small disorder and defects in the material. Previous work about Raman studies in graphene on SiO₂/Si show that the shape of the 2D peak is an indicator of the number of layers [37]. For Raman analysis, the nature of the growth (exfoliation, CVD graphitization, CVD without graphitization) and the type of the substrates (SiO₂ or SiC) can affect the peaks' position and the value of the full width at half maximum (FWHM), which enhance the need to have more results about bilayer graphene on SiC [26,38,39]. Here, the full width at the half maximum of the 2D peak is around 59 cm⁻¹, comparable to the 41–62 cm⁻¹ reported in bilayer graphene on SiC and other substrates [25,37]. In bilayer graphene, the 2D peak is typically broader and upshifted compared to monolayer graphene. This broadening is due to the presence of an additional phonon mode in bilayer graphene, resulting in a more complex 2D band shape [37]. Before fabrication, Raman spectra was established in three different locations randomly selected on the SiC wafer and presented in Figure 1e. The Raman spectrum does not change with changing locations.

B. Device fabrication.

The field effect transistors based on graphene were fabricated on a 15 × 15 mm² SiC wafer. Figure 2a represents the layout of the device. Each level of the layout is represented by a different color or contrast and represents a fabrication step of the process. First, the process is fixed by defining the alignment marks. It follows the etching of the graphene channel and contacts as a hole for improving contact resistance, as previously reported in [41,42]. The source and drain contacts were obtained by the standard lift-off process after evaporating 1.5 nm of nickel and 30 nm of gold metals. Here, a thin layer of nickel (1.5 nm) was deposited before in order to improve the metal adhesion on the surface. The dual T-gate with gate length ($L_g$) were defined by using three layers of poly-meta-methacrylate (PMMA), as shown in Figure 2b, where three different thicknesses were defined: 160 nm thick at the bottom part of the T-gate, 720 nm thick at the top part of the T-gate and 130 nm at the resist followed by electron beam (e-beam) lithography. After the development of these multilayers' resists, the gate oxide is deposited using 2 nm of evaporated aluminum four times, followed by oxidation in ambient air for 24 h. Finally, the coplanar access
Ni (50 nm)/Au (300 nm) are deposited, followed by lift-off. The cross-section schema of the final transistor is presented in Figure 2c. The top view of the active part of the final dual-T-gate-transistor is illustrated in the scanning electron microscopy (SEM) image in Figure 2d (left) while the cross-section of the gate part was illustrated by the focused ion beam (FIB) technique in Figure 2d (right) and shows the shape of the final T-gate of the GFET. A picture of the 15 × 15 mm² final wafer where there are 8 cells and 458 transistors is shown in Figure 2e.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** AFM images of bilayer graphene surface on 6H-SiC substrate [40]. (a) 60 × 60 µm² image and the blue and red dash square 10 × 10 µm² is represented in (b). (c) SEM image with 10 µm scale bar. (d) Raman spectroscopy of graphene on SiC substrate. The black trace is the spectrum of graphene and SiC, the red trace is the graphene spectrum once SiC Raman peaks are subtracted, the blue trace is the smoothed spectrum. (e) Raman spectra at different locations (blue, red, and black) on the SiC wafer.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** (a) A picture of the layout of the device. “+” is source-drain contact layer. (b) A schematic image representing the step of the realization of the T gate process with the three layers’ resists. (c) A schematic of the side view of the T gate transistor. (d) (d, left) An SEM image of the transistor in the end of the process. (d, right) An FIB image of the transistor showing the T-gate structure, $L_g$ = 200 nm. (e) A photograph of the final devices fabricated on the 15 mm × 15 mm SiC substrate, including 458 transistors and height cells numbered from 1 to 8.</figcaption></figure>

**3. Results and Discussion**

After fabrication, we conducted an in-depth characterization of the electrical properties of graphene. The wafer was tested within patterned structures, including Hall measurements and transmission line measurement (TLM), as illustrated in Figure 3a, across various cells. Contact and sheet resistance, which represent the resistance as a function of the distance between contact points (0.5, 1, 2, 4, 16, 24 µm), were extracted from a linear fit using the TLM method. The y-axis intercept provides the contact resistance (2Rc), while the slope of the linear fit represents the Rsh ratio to the width of the graphene channel [41,43]. All extracted values of contact resistance and sheet resistance for each cell are presented in Figure 3d.

To complete the analysis, we acquired seven Raman spectra from each cell of the wafer and extracted an average Raman spectrum. Figure 3c summarizes the average Raman spectra representing each cell. Mobility values, carrier densities and Raman peaks are provided in the table in Figure 3d. It shows a dispersion of values for contact resistance, sheet resistance, Hall mobility and carrier density across the wafer. The highest Hall
mobility (626–832 cm²/V·s), carrier densities 18–16 (×10¹² cm⁻²), lowest values of contact resistance (678–650 Ω·m) and sheet resistance (323–286 Ω/sq) were obtained in the cells 7–8. However, the low value of mobility observed in cells 1–6 is possibly due to the high contact resistance. The ohmic contact between graphene and metal contacts exhibits large contact resistance, significantly reducing the apparent mobility of contacted graphene and hindering its potential in high-frequency applications [42].

<figure><img src="image_6.png" /></figure>

<table><caption>**Figure 3.** (a) An optical image of the Hall and TLM pattern. (b) The transmission line measurement (TLM) and the linear fit to extract the correspondent contact resistance and sheet resistance. (c) A recapitulation of the averaged Raman spectra of graphene on SiC measured in the eight cells of the wafer. (d) A table summarizing the values of the Hall mobility, carrier densities and Raman peaks' positions and the full width at the half maximum of the G and 2D peaks. Missing values are due to faulty components.</caption><thead><tr><th>Cell No</th><th>Carrier densities x10<sup>12</sup> (cm<sup>-2</sup>)</th><th>Hall Mobility (cm²/V·s)</th><th>Contact resistance (Ω·μm)</th><th>Sheet Resistance (Ω/sq)</th><th>G peak position (cm<sup>-1</sup>)</th><th>2D peak position (cm<sup>-1</sup>)</th><th>FWHM (G)</th><th>FWHM (2D)</th></tr></thead><tbody><tr><td>1</td><td>4</td><td>360</td><td>1689</td><td>5175</td><td>1591</td><td>2750</td><td>33</td><td>65</td></tr><tr><td>2</td><td>5</td><td>439</td><td>1335</td><td>1582</td><td>1590</td><td>2751</td><td>19</td><td>65</td></tr><tr><td>3</td><td>4</td><td>367</td><td>-</td><td>-</td><td>1590</td><td>2750</td><td>24</td><td>57</td></tr><tr><td>4</td><td>-</td><td>-</td><td>1971</td><td>2845</td><td>1593</td><td>2748</td><td>16</td><td>57</td></tr><tr><td>5</td><td>10</td><td>281</td><td>-</td><td>-</td><td>1590</td><td>2744</td><td>16</td><td>54</td></tr><tr><td>6</td><td>6</td><td>104</td><td>-</td><td>-</td><td>1592</td><td>2749</td><td>17</td><td>55</td></tr><tr><td>7</td><td>18</td><td>625</td><td>678</td><td>323</td><td>1590</td><td>2738</td><td>14</td><td>57</td></tr><tr><td>8</td><td>16</td><td>832</td><td>650</td><td>286</td><td>1591</td><td>2740</td><td>14</td><td>57</td></tr></tbody></table>

From the Raman analysis, the position of the G peak remains consistent across the wafer, while the full width at the half maximum of the G peak FWHM (G) varies. The thinnest FWHM (G) is observed in cells 7 and 8. The 2D peak exhibits a blueshift of approximately 20 cm⁻¹, from 2751 cm⁻¹ in cell 2 to 2731 cm⁻¹ in cell 7. This shift could be attributed to interlayer coupling effects in bilayer graphene. Ciuk et al. [25] investigated bilayer CVD graphene on SiC and observed the blueshift of the 2D peak and redshift of the
G peak due to the strain at the edges of the steps. Moreover, Das et al. [44] demonstrated that in graphene on SiO$_{2}$/Si, the position of the G peak increased, while the 2D peak and the FWHM (G) decreased with increased dopant concentration. Ferrari et al. [37] showed that the electron and hole doping upshifts and sharpens the G peak. Additionally, graphene can be doped chemically during fabrication, affecting carrier densities, 2D-peak width and position. Previous studies have shown that doping increases the 2D position due to the modification of lattice parameters, which modifies the total number of charges and leads to a stiffening/softening of the phonons [45]. The analysis of the table in Figure 3d clearly shows that cells 7–8 present the highest mobility, carrier densities, lowest contact and sheet resistance and smallest FWHM (G), suggesting favorable prospects for effective device and circuit applications.

To evaluate the performance of our devices, DC and RF measurements were performed using a standard probe station with Microtech's probes, Semiconductor Analyzer HP4155A, the Vector Network Analyzer HP4155A and the vector Network Analyzer Rohde & Schwarz ZVA67. A common calibration procedure of Line-Reflect-Reflect-Match (LRRM) for RF measurements was established before measurements. In total, 195 transistors of the GFETs were measured.

Figure 4 represents an example of the electrical characterization of a transistor in cell 7 (Figure 2e), with a 200 nm gate length and a 30 × 2 µm channel width. The transfer characteristic $I_{DS}$ as function of $V_{GS}$ at $V_{DS} = 1.5$ V and the transconductance ($g_m = dI_{DS}/dV_{GS}$) are shown in Figure 4a. The $I_{DS}$–$V_{GS}$ curve is non-monotonic. As $V_{GS}$ changes from –3 V to +3 V, $I_{DS}$ shows a change in slopes initially decreasing, then reversing direction. This non-monotonic behavior is accurately reflected in the $g_m$ curve, calculated as the derivative of $I_{DS}$ as a function of $V_{GS}$. The characteristic of decreasing and then increasing $g_m$ is consistent with our observations and not indicative of bipolar behavior as observed in [46]. The on/off extracted is approximately 1.4, within the $V_{GS}$ range of ±3 V. The peak $g_m$ reaches 4.6 mS (76 mS/mm) at $V_{GS} = 1.3$ V. The Dirac point is located at $V_{GS}$ higher than 3 V, suggesting p-type doping in the graphene channel. Figure 4b shows the output characteristics, sweeping the gate voltage from –2 V to +2 V in 0.5 V steps. The maximum current is around 53 mA (0.88 mA/µm) at $V_{DS} = +1.5$ V and $V_{GS} = –2$ V, with no observed saturation current. This work investigates graphene's potential for high-frequency analog applications, focusing on its high-frequency characteristics rather than switching capability and looking for it to reach the saturation current. Graphene's semi-metallic nature results in a low on/off ratio, which is less critical for RF amplifiers compared to synaptic transistors or for digital applications [12,46]. Due to its high carrier mobility and conductivity, graphene excels in GHz frequency operations, maintaining steady current mobility flow at specific bias points [47].

<figure><img src="image_7.png" /></figure>

**Figure 4.** Example of DC and RF characteristics of dual T-gate graphene transistor in cell 7 having gate channel length ($L_g$) equal to 200 nm and dual-gate width ($w$) 2 × 30 µm. (a) DC measurement $I_{DS}$–$V_{GS}$ and $g_m$–$V_{GS}$. (b) The voltage transfer characteristics as a function of $V_{GS}$ varying from –2 V to +2 V by 0.5 V steps. (c) RF characteristic includes the as-measured values of the current gain $H_{21}$ and the unilateral power gain $U$ as a function of the frequency at $V_{DS} = 1.5$ V and $V_{GS} = 1.3$ V. The cut off frequency $f_t$ and the maximum oscillation frequency $f_{max}$ have been extracted [40].
The extrinsic RF characterization is reported in Figure 4c. It shows the current gain H21 and the Mason's gain U of the device in a frequency range of 0.6 to 67 GHz. The maximum frequency of oscillation (fmax) and the current gain cut-off frequency (ft) are, respectively, the frequency at which the power gain (U) and the current gain (H21) are equal to 1. The on-probe measurement of the cut-off frequency and the maximum oscillation frequency reveal, respectively, ft = 21 GHz and fmax = 18 GHz at VGS = +1 V and the gain curves follow the slope of 20 dB/decade, as expected. The extrinsic value of the performance is obtained after the capacitances related to the length of the transmission line of transistor access are removed [48]. They are, respectively, fT-extr = 62 GHz and fmax-extr remains the same.

These values are comparable to the recent value achieved by [32] in the bilayer graphene on SiC.

To understand the evolution of RF performance across the eight cells of the full wafer shown in Figure 5a, we analyzed the average values of $f_t$, $f_{max}$ and the optimum gate voltage for each individual cell. The optimum gate voltages for RF measurement were determined from the $g_m-V_{GS}$ characteristic, where the transconductance is at its maximum. The average gate voltage $<V_{GS}>$ is shown in Figure 5b and represents the mean gate voltage computed in each cell. The average cut-off frequency $<f_t>$ and maximum oscillation frequency $<f_{max}>$ are reported in Figure 5c,d, respectively. They clearly show that the best performance is achieved in cells 7 and 8, with $<f_t>$ of 16–14 GHz and $<f_{max}>$ of 11–9 GHz, respectively. Figure 5e summarizes the $f_t$ and $f_{max}$ values for each cell as a function of the optimum gate voltage. Interestingly, the best transistor performances are obtained in cells 7 and 8 for gate voltages of at least 1 V. This variation in gate voltage may be attributed to the variation of the dopant across the wafer as mentioned in the Raman analysis. Additionally, differences in contact resistance across various cells contribute to the observed variability. To mitigate this issue, implementing a robust cleaning process is essential to eliminate polymer contamination during fabrication. This approach will reduce contact resistance, enhance uniformity, and minimize performance variability. It is important to acknowledge that both contact resistance and graphene quality are important for overall device performance. However, our study suggests that the RF performance is immediately limited by contact resistance related to the process of fabrication.

The findings highlight the critical role of the fabrication process in the RF performance variability of graphene FETs fabricated from high-quality epitaxial bilayer graphene on a SiC substrate. Despite the material's exceptional quality, confirmed through AFM and Raman spectroscopy, the observed dispersion in values of Hall mobility, carrier densities and contact resistance across the wafer suggest that these performance discrepancies result from the fabrication stages rather than inherent material inconsistencies. This understanding underscores the need for refined manufacturing techniques to ensure consistent transistor behavior, which is essential for the reliable operation of integrated circuits that rely on uniform performance of GFET in wafer scale. The insights gained from this study are essential for advancing GFET technology, particularly in high-frequency applications where reproducible device characteristics are crucial.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** (a) Photography of the 15 mm × 15 mm graphene wafer presented previously in Figure 2e. (b) $V_{GS}$ is the bias gate voltage related to the maximum transconductance and where the best performance of the transistor is expected. $<V_{GS}>$ is the average value of all the gate voltage measured in each cell. (c,d) Representation of the average values of the cut-off frequency $<f_t>$ and the maximum oscillation frequency $<f_{max}>$ computed for each cell. (e) Graph summarizing the evolution of the on-probe values of $f_t$ and $f_{max}$ as function of the biased gate voltage $V_{GS}$ of the transistors in each cell.</figcaption></figure>

4. Conclusions

In this paper, we fabricated hundreds of bilayer graphene field effect transistors on a SiC substrate. We provided Raman spectrum analysis combined with the Hall mobility, carrier densities, contact resistance and sheet resistance toward a 15 mm × 15 mm graphene wafer. The analyses of local Raman have shown that the 2D peak and the FWHM (G) changed, as well as the Hall mobility, carrier densities and contact resistance in the wafer. The analysis of the RF performances was performed and compared to the local properties of the wafer. It revealed dispersion of the performances and correlation between the RF performance, contact resistance and Hall mobilities. This work completes the few works about local Raman investigation on bilayer graphene on SiC and highlights the importance of local analyses of the properties of the material to evaluate the performance of electronic devices where large-scale, homogenous and high-quality fabrication of the material is needed for high-performance circuit.
**Author Contributions:** D.F. is the main author of this work. She assumed the design and the fabrication of the devices, the DC and RF characterizations, as well as the preparation of the graphs of the manuscript. W.S. provided the graphene material used in this work. The CVD bilayer graphene used was grown in his research group. E.P. and H.H. initiated this work and contributed to high-frequency characterizations and interpretations of the results. All authors reviewed and commented on the manuscript. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Horizon 2020 Framework Program: GrapheneCore2 785219 and GrapheneCore3 881603.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data is contained within the article.

**Acknowledgments:** The authors gratefully acknowledge financial support from the EU FP7-ICT-2013-FET-F GRAPHENE Flagship project (GrapheneCore2 no. 785219 and GrapheneCore3 n° 881603). This work was also partly supported by the French RENATECH network. D.F. is currently funded from the European Union's Horizon 2020 research and innovation program under the Marie Skłodowska Curie grant agreement (no. 101025770).

**Conflicts of Interest:** The authors declare no conflicts of interest.

## References

1. Geim, A.K. Graphene: Status and Prospects. *Science* **2009**, *324*, 1530–1534. [CrossRef] [PubMed]

2. Kaul, A.B. Two-dimensional layered materials: Structure, properties, and prospects for device applications. *J. Mater. Res.* **2014**, *29*, 348–361. [CrossRef]

3. Ferrari, A.C.; Bonaccorso, F.; Fal'Ko, V.; Novoselov, K.S.; Roche, S.; Bøggild, P.; Borini, S.; Koppens, F.H.L.; Palermo, V.; Pugno, N.; et al. Science and technology roadmap for graphene, related two-dimensional crystals, and hybrid systems. *Nanoscale* **2015**, *7*, 4598–4810. [CrossRef] [PubMed]

4. Zhang, H.; Chhowalla, M.; Liu, Z. 2D nanomaterials: Graphene and transition metal dichalcogenides. *Chem. Soc. Rev.* **2018**, *47*, 3015–3017. [CrossRef] [PubMed]

5. Bolotin, K.; Sikes, K.; Jiang, Z.; Klima, M.; Fudenberg, G.; Hone, J.; Kim, P.; Stormer, H. Ultrahigh electron mobility in suspended graphene. *Solid State Commun.* **2008**, *146*, 351–355. [CrossRef]

6. Pizzocchero, F.; Gammelgaard, L.; Jessen, B.S.; Caridad, J.M.; Wang, L.; Hone, J.; Bøggild, P.; Booth, T.J. The hot pick-up technique for batch assembly of van der Waals heterostructures. *Nat. Commun.* **2016**, *7*, 11894. [CrossRef] [PubMed]

7. Koppens, F.H.L.; Mueller, T.; Avouris, P.; Ferrari, A.C.; Vitiello, M.S.; Polini, M. Photodetectors based on graphene, other two-dimensional materials and hybrid systems. *Nat. Nanotechnol.* **2014**, *9*, 780–793. [CrossRef] [PubMed]

8. Balandin, A.A. Thermal properties of graphene and nanostructured carbon materials. *Nat. Mater.* **2011**, *10*, 569–581. [CrossRef]

9. Lee, C.; Wei, X.; Kysar, J.W.; Hone, J. Measurement of the elastic properties and intrinsic strength of monolayer graphene. *Science* **2008**, *321*, 385–388. [CrossRef]

10. Schwierz, F. Graphene transistors. *Nat. Nanotechnol.* **2010**, *5*, 487–496. [CrossRef]

11. Wu, Y.; Jenkins, K.A.; Valdes-Garcia, A.; Farmer, D.B.; Zhu, Y.; Bol, A.A.; Dimitrakopoulos, C.; Zhu, W.; Xia, F.; Avouris, P.; et al. State-of-the-Art Graphene High-Frequency Electronics. *Nano Lett.* **2012**, *12*, 3062–3067. [CrossRef] [PubMed]

12. Schwierz, F. Graphene Transistors: Status, Prospects, and Problems. *Proc. IEEE* **2013**, *101*, 1567–1584. [CrossRef]

13. Novoselov, K.S.; Geim, A.K.; Morozov, S.V.; Jiang, D.; Zhang, Y.; Dubonos, S.V.; Grigorieva, I.V.; Firsov, A.A. Electric field effect in atomically thin carbon films. *Science* **2004**, *306*, 666–669. [CrossRef] [PubMed]

14. Novoselov, K.S.; Fal'Ko, V.I.; Colombo, L.; Gellert, P.R.; Schwab, M.G.; Kim, K. A roadmap for graphene. *Nature* **2012**, *490*, 192–200. [CrossRef] [PubMed]

15. Tielrooij, K.-J.; Hesp, N.C.H.; Principi, A.; Lundeberg, M.B.; Pogna, E.A.A.; Banszerus, L.; Mics, Z.; Massicotte, M.; Schmidt, P.; Davydovskaya, D.; et al. Out-of-plane heat transfer in van der Waals stacks through electron–hyperbolic phonon coupling. *Nat. Nanotechnol.* **2018**, *13*, 41–46. [CrossRef] [PubMed]

16. Kim, K.S.; Zhao, Y.; Jang, H.; Lee, S.Y.; Kim, J.M.; Kim, K.S.; Ahn, J.-H.; Kim, P.; Choi, J.-Y.; Hong, B.H. Large-scale pattern growth of graphene films for stretchable transparent electrodes. *Nature* **2009**, *457*, 706–710. [CrossRef] [PubMed]

17. Deokar, G.; Avila, J.; Razado-Colambo, I.; Codron, J.-L.; Boyaval, C.; Galopin, E.; Asensio, M.-C.; Vignaud, D. Towards high quality CVD graphene growth and transfer. *Carbon* **2015**, *89*, 82–92. [CrossRef]

18. Smith, A.D.; Wagner, S.; Kataria, S.; Malm, B.G.; Lemme, M.C.; Ostling, M. Wafer-Scale Statistical Analysis of Graphene FETs—Part I: Wafer-Scale Fabrication and Yield Analysis. *IEEE Trans. Electron. Devices* **2017**, *64*, 3919–3926. [CrossRef]

19. Berger, C.; Song, Z.M.; Li, T.B.; Li, X.B.; Ogbazghi, A.Y.; Feng, R.; Dai, Z.N.; Marchenkov, A.N.; Conrad, E.H.; First, P.N.; et al. Ultrathin Epitaxial Graphite: 2D Electron Gas Properties and a Route toward Graphene-based Nanoelectronics. *J. Phys. Chem. B* **2004**, *108*, 19912–19916. [CrossRef]
20. Wu, Y.Q.; Ye, P.D.; Capano, M.A.; Xuan, Y.; Sui, Y.; Qi, M.; Cooper, J.A.; Shen, T.; Pandey, D.; Prakash, G.; et al. Top-gated graphene field-effect-transistors formed by decomposition of SiC. *Appl. Phys. Lett.* **2008**, *92*, 092102. [CrossRef]

21. Ferrer, F.J.; Moreau, E.; Vignaud, D.; Godey, S.; Wallart, X. Atomic scale flattening, step formation and graphitization blocking on 6H- and 4H-SiC[0 0 0 1] surfaces under Si flux. *Semicond. Sci. Technol.* **2009**, *24*, 125014. [CrossRef]

22. Strupinski, W.; Grodecki, K.; Wysmolek, A.; Stepniewski, R.; Szkopek, T.; Gaskell, P.E.; Grüneis, A.; Haberer, D.; Bozek, R.; Krupka, J.; et al. Graphene Epitaxy by Chemical Vapor Deposition on SiC. *Nano Lett.* **2011**, *11*, 1786–1791. [CrossRef] [PubMed]

23. Guo, Z.; Dong, R.; Chakraborty, P.S.; Lourenco, N.; Palmer, J.; Hu, Y.; Ruan, M.; Hankinson, J.; Kunc, J.; Cressler, J.D.; et al. Record Maximum Oscillation Frequency in C-Face Epitaxial Graphene Transistors. *Nano Lett.* **2013**, *13*, 942–947. [CrossRef] [PubMed]

24. Pallecchi, E.; Lafont, F.; Cavaliere, V.; Schopfer, F.; Mailly, D.; Poirier, W.; Ouerghi, A. High Electron Mobility in Epitaxial Graphene on 4H-SiC(0001) via post-growth annealing under hydrogen. *Sci. Rep.* **2014**, *4*, 4558. [CrossRef] [PubMed]

25. Ciuk, T.; Cakmakyapan, S.; Ozbay, E.; Caban, P.; Grodecki, K.; Krajewska, A.; Pasternak, I.; Szmidt, J.; Strupinski, W. Step-edge-induced resistance anisotropy in quasi-free-standing bilayer chemical vapor deposition graphene on SiC. *J. Appl. Phys.* **2014**, *116*, 123708. [CrossRef]

26. Ciuk, T.; Caban, P.; Strupinski, W. Charge carrier concentration and offset voltage in quasi-free-standing monolayer chemical vapor deposition graphene on SiC. *Carbon* **2016**, *101*, 431–438. [CrossRef]

27. Ciuk, T.; Strupinski, W. Statistics of epitaxial graphene for Hall effect sensors. *Carbon* **2015**, *93*, 1042–1049. [CrossRef]

28. Zhou, S.; Gweon, G.-H.; Fedorov, A.; First, P.N.; De Heer, W.A.; Lee, D.-H.; Guinea, F.; Neto, A.H.C.; Lanzara, A. Substrate-induced bandgap opening in epitaxial graphene. *Nat. Mater.* **2007**, *6*, 770–775, Corrigendum in **2007**, *6*, 916. [CrossRef]

29. Zhang, Y.; Tang, T.-T.; Girit, C.; Hao, Z.; Martin, M.C.; Zettl, A.; Crommie, M.F.; Shen, Y.R.; Wang, F. Direct observation of a widely tunable bandgap in bilayer graphene. *Nature* **2009**, *459*, 820–823. [CrossRef]

30. Xu, X.; Liu, C.; Sun, Z.; Cao, T.; Zhang, Z.; Wang, E.; Liu, Z.; Liu, K. Interfacial engineering in graphene bandgap. *Chem. Soc. Rev.* **2018**, *47*, 3059–3099. [CrossRef]

31. He, Z.-Z.; Yang, K.-W.; Yu, C.; Liu, Q.-B.; Wang, J.-J.; Song, X.-B.; Han, T.-T.; Feng, Z.-H.; Cai, S.-J. Comparative Study of Monolayer and Bilayer Epitaxial Graphene Field-Effect Transistors on SiC Substrates. *Chin. Phys. Lett.* **2016**, *33*, 086801. [CrossRef]

32. Yu, C.; He, Z.Z.; Song, X.B.; Liu, Q.B.; Han, T.T.; Dun, S.B.; Wang, J.J.; Zhou, C.J.; Guo, J.C.; Lv, Y.J.; et al. Improvement of the Frequency Characteristics of Graphene Field-Effect Transistors on SiC Substrate. *IEEE Electron. Device Lett.* **2017**, *38*, 1339–1342. [CrossRef]

33. Neumann, C.; Reichardt, S.; Venezuela, P.; Drögeler, M.; Banszerus, L.; Schmitz, M.; Watanabe, K.; Taniguchi, T.; Mauri, F.; Beschoten, B.; et al. Raman spectroscopy as probe of nanometre-scale strain variations in graphene. *Nat. Commun.* **2015**, *6*, 8429. [CrossRef] [PubMed]

34. Long, F.; Yasaei, P.; Sanoj, R.; Yao, W.; Král, P.; Salehi-Khojin, A.; Shahbazian-Yassar, R. Characteristic Work Function Variations of Graphene Line Defects. *ACS Appl. Mater. Interfaces* **2016**, *8*, 18360–18366. [CrossRef] [PubMed]

35. Buron, J.D.; Pizzocchero, F.; Jepsen, P.U.; Petersen, D.H.; Caridad, J.M.; Jessen, B.S.; Booth, T.J.; Bøggild, P. Graphene mobility mapping. *Sci. Rep.* **2015**, *5*, 12305. [CrossRef] [PubMed]

36. Wehrfritz, P.; Fromm, F.; Malzer, S.; Seyller, T. Quasi-freestanding epitaxial graphene transistor with silicon nitride top gate. *J. Phys. D Appl. Phys.* **2014**, *47*, 305103. [CrossRef]

37. Ferrari, A.C. Raman spectroscopy of graphene and graphite: Disorder, electron-phonon coupling, doping and nonadiabatic effects. *Solid State Commun.* **2007**, *143*, 47–57. [CrossRef]

38. Wang, Y.Y.; Ni, Z.H.; Yu, T.; Shen, Z.X.; Wang, H.M.; Wu, Y.H.; Chen, W.; Wee, A.T.S. Raman Studies of Monolayer Graphene: The Substrate Effect. *J. Phys. Chem. C* **2008**, *112*, 10637–10640. [CrossRef]

39. Röhrl, J.; Hundhausen, M.; Emtsev, K.V.; Seyller, T.; Graupner, R.; Ley, L. Raman spectra of epitaxial graphene on SiC(0001). *Appl. Phys. Lett.* **2008**, *92*, 201918. [CrossRef]

40. Fadil, D.; Wei, W.; Deng, M.; Fregonese, S.; Strupinski, W.; Pallecchi, E.; Happy, H. 2D-Graphene Epitaxy on SiC for RF Application: Fabrication, Electrical Characterization and Noise Performance. In Proceedings of the 2018 IEEE/MTT-S International Microwave Symposium—IMS, Philadelphia, PA, USA, 10–15 June 2018; pp. 228–231.

41. Passi, V.; Gahoi, A.; Ruhkopf, J.; Kataria, S.; Vaurette, F.; Pallecchi, E.; Happy, H.; Lemme, M.C. Contact resistance Study of “edge-contacted” metal-graphene interfaces. In Proceedings of the ESSDERC 2016—46th European Solid-State Device Research Conference, Lausanne, Switzerland, 12–15 September 2016; pp. 236–239.

42. Anzi, L.; Mansouri, A.; Pedrinazzi, P.; Guerriero, E.; Fiocco, M.; Pesquera, A.; Centeno, A.; Zurutuza, A.; Behnam, A.; A Carrion, E.; et al. Ultra-low contact resistance in graphene devices at the Dirac point. *2D Mater.* **2018**, *5*, 025014. [CrossRef]

43. Watanabe, E.; Conwill, A.; Tsuya, D.; Koide, Y. Low contact resistance metals for graphene based devices. *Diam. Relat. Mater.* **2012**, *24*, 171–174. [CrossRef]

44. Das, A.; Pisana, S.; Chakraborty, B.; Piscanec, S.; Saha, S.K.; Waghmare, U.V.; Novoselov, K.S.; Krishnamurthy, H.R.; Geim, A.K.; Ferrari, A.C.; et al. Monitoring dopants by Raman scattering in an electrochemically top-gated graphene transistor. *Nat. Nanotechnol.* **2008**, *3*, 210–215. [CrossRef]

45. Bruna, M.; Ott, A.K.; Ijäs, M.; Yoon, D.; Sassi, U.; Ferrari, A.C. Doping Dependence of the Raman Spectrum of Defected Graphene. *ACS Nano* **2014**, *8*, 7432–7441. [CrossRef]

46. Ni, Y.; Ma, M.; Wei, H.; Gong, J.; Han, H.; Liu, L.; Xu, Z.; Xu, W. Multiplexed neurotransmission emulated for emotion control. *Nano Energy* **2021**, *86*, 106038. [CrossRef]
47. Asad, M.; Bonmann, M.; Yang, X.; Vorobiev, A.; Jeppson, K.; Banszerus, L.; Otto, M.; Stampfer, C.; Neumaier, D.; Stake, J. The Dependence of the High-Frequency Performance of Graphene Field-Effect Transistors on Channel Transport Properties. *IEEE J. Electron. Devices Soc.* **2020**, *8*, 457–464. [CrossRef]

48. Nougaret, L.; Happy, H.; Dambrine, G.; Derycke, V.; Bourgoin, J.-P.; Green, A.A.; Hersam, M.C. 80 GHz field-effect transistors produced using high purity semiconducting single-walled carbon nanotubes. *Appl. Phys. Lett.* **2009**, *94*, 243505. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.