# High-frequency limits of graphene field-effect transistors with velocity saturation

Quentin Wilmart, Mohamed Boukhicha, Holger Graef, David Mele, José Palomo, Michael Rosticher, Takashi Taniguchi, Kenji Watanabe, Vincent Bouchiat, Emmanuel Baudin, et al.

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

## To cite this version:

Quentin Wilmart, Mohamed Boukhicha, Holger Graef, David Mele, José Palomo, et al. High-frequency limits of graphene field-effect transistors with velocity saturation. Applied Sciences, 2020, 10 (2), pp.446. 10.3390/app10020446. hal-02520621

HAL Id: hal-02520621

https://hal.science/hal-02520621

Submitted on 26 Mar 2020

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_2.png" /></figure>

Article

# High-Frequency Limits of Graphene Field-Effect Transistors with Velocity Saturation

Quentin Wilmart<fnref n="1" />, Mohamed Boukhicha <fnref n="1" />,<fnref n="2" /> ID, Holger Graef <fnref n="1" /> ID, David Mele <fnref n="1" /> ID, Jose Palomo <fnref n="1" />, Michael Rosticher <fnref n="1" />, Takashi Taniguchi <fnref n="3" />, Kenji Watanabe <fnref n="3" /> ID, Vincent Bouchiat <fnref n="4" />, Emmanuel Baudin <fnref n="1" />, Jean-Marc Berroir <fnref n="1" /> ID, Erwann Bocquillon <fnref n="1" />, Gwendal Fève <fnref n="1" />, Emiliano Pallecchi <fnref n="5" /> and Bernard Placais <fnref n="1" />,* ID

Received: 29 November 2019; Accepted: 31 December 2019; Published: 8 January 2020

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

check for
updates

**Abstract:** The current understanding of physical principles governing electronic transport in graphene field effect transistors (GFETs) has reached a level where we can model quite accurately device operation and predict intrinsic frequency limits of performance. In this work, we use this knowledge to analyze DC and RF transport properties of bottom-gated graphene on boron nitride field effect transistors exhibiting pronounced velocity saturation by substrate hyperbolic phonon polariton scattering, including Dirac pinch-off effect. We predict and demonstrate a maximum oscillation frequency exceeding 20 GHz. We discuss the intrinsic 0.1 THz limit of GFETs and envision plasma resonance transistors as an alternative for sub-THz narrow-band detection.

**Keywords:** Graphene transistor; velocity saturation; high-frequency; Dirac pinch-off; contact gating

## 1. Introduction

Graphene is considered as a high-potential material for electronics, especially radio-frequency (RF) field effect transistors (GFETs) [1], owing to the high intrinsic mobility of charge carriers [2–5]. The main challenge in this respect is to achieve current saturation at high bias to obtain significant voltage and power gain. Two routes are considered. The first relies on band-gap engineering using, e.g., lateral confinement in graphene nanoribbon devices [6], which is often obtained at the cost of a reduced mobility. As in MOSFETs [7], the saturation results from a potential barrier induced at the drain side. The second route is velocity saturation, which promises high performances [8–10]. It requires high mobility transistors ($\mu \gtrsim 1 \text{ m}^2\text{V}^{-1}\text{s}^{-1}$) and results in large saturation currents ($I_{sat} \sim 20 \text{ Acm}^{-1}$) [11]. Limited by the Fermi energy at low doping, the saturation velocity is determined by optical phonon (OPs) scattering at large doping when $E_F$ exceeds the optical phonon energy $\hbar\Omega_{OP}$ [11,12]. In oxide-supported graphene, relevant OPs are substrate polar phonons (SPPs) [13–15] that have a moderate energy $\hbar\Omega_{SPP} \lesssim 50$ meV. The saturation current being proportional

<footnote n="1">¹ Laboratoire de Physique de l'Ecole Normale Supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université Paris-Diderot, Sorbonne Paris Cité, 75005 Paris, France; quentin.wilmart@gmail.com (Q.W.); med.boukhicha@yahoo.fr (M.B.); holger.graef@lpa.ens.fr (H.G.); david.mele@lpa.ens.fr (D.M.); jose.palomo@ens.fr (J.P.); michael.rosticher@phys.ens.fr (M.R.); emmanuel.baudin@phys.ens.fr (E.B.); jean-marc.berroir@ens.fr (J.-M.B.); Erwann.Bocquillon@phys.ens.fr (E.B.); gwendal.feve@phys.ens.fr (G.F.)</footnote>

<footnote n="2">² Sustainable Energy Technologies Department, Brookhaven National Laboratory, Upton, NY 11973, USA</footnote>

<footnote n="3">³ Advanced Materials Laboratory, National Institute for Materials Science, Tsukuba, Ibaraki 305-0047, Japan; taniguchi.takashi@nims.go.jp (T.T.); WATANABE.Kenji.aml@nims.go.jp (K.W.)</footnote>

<footnote n="4">⁴ Institut Néel, CNRS and Université Grenoble-Alpes, 38042 Grenoble, France; Vincent.Bouchiat@grenoble.cnrs.fr</footnote>

<footnote n="5">5 Institut d'Électronique, de Microélectronique et de Nanotechnologie, CNRS and University of Lille, 59491 Villeneuve d'Asq, France; Emiliano.Pallecchi@iemn.univ-lille1.fr</footnote>

<footnote>* Correspondence: placais@phys.ens.fr or bernard.placais@phys.ens.fr</footnote>
to the OP energy, power applications favor non-oxide dielectrics such as SiC ($\hbar\Omega_{SiC} \approx 116$ meV), or hexagonal boron nitride (h-BN), which has two Reststrahlen bands, $\hbar\Omega_I \approx 95–100$ meV and $\hbar\Omega_{II} \approx 170–200$ meV, associated respectively with out-of-plane and in-plane optical phonons. Uniaxial dielectrics, such as h-BN, actually sustain hyperbolic phonon polariton modes (HPhPs) that differ from SPPs (and intrinsic OPs), insofar as they propagate deep in the dielectric bulk, providing a unique radiative cooling pathway [11,16]. This paper discusses velocity saturation GFETs, their optimal design and operating conditions, their intrinsic frequency limit, and envisions the possibility to overcome this limit using plasma resonance devices [17,18]. It relies on experimental data taken on two series of h-BN-supported GFETs based on CVD graphene, which is appropriate for wafer scale applications.

## 2. Methods

Samples used in this work are made of CVD graphene on a thin h-BN flake deposited on top of a metallic bottom gate. The substrate consists in a 280 nm thermal SiO₂ insulating barrier on top of high-resistivity Si suitable for RF measurements. Two types of gate electrodes are used: a thin tungsten electrode deposited on SiO₂ (type A samples) for low-frequency devices and DC characterization, and a thick gold dual gate buried in SiO₂ (type B samples) for RF devices. Several samples of both types have been measured presenting similar characteristics, but for simplicity, we focus on one illustrative sample of each type.

Sample-A (Figure 1a) is made out of a 20-nm-thin tungsten film that is etched to define the gate whose dimensions are $L = 0.5$ µm for the length and $W = 1.5$ µm for the width. In addition, we have included local contact gates (also tungsten) controlling the carrier density of the contacted graphene to reduce the contact resistance (see [19]). An h-BN flake, of thickness $t_{BN} = 16$ nm, is deposited on top of the gate by a dry transfer technique [2]. The h-BN surface is cleaned by furnace annealing at 300 °C under an Ar/H₂ flow. The CVD graphene is grown on copper following the *pulsed CVD* method at Institut Néel [3] and transferred with a wet technique on top of h-BN. Again, the polymer residues are removed by furnace annealing. Finally, 50 nm of Pd is evaporated for the electrical contacts, and devices are embedded in a 50 Ω-matched coplanar waveguide (CPW) used for DC and RF characterization. A schematic of a type-A device can be seen in Figure 1b. In very-high mobility devices such as the Zener-Klein transistor of [11] radiative cooling by hyperbolic phonon polaritons in the can be achieved, which is optimized using thicker ($t_{BN} \gtrsim 100$ nm) h-BN substrates. Note also that alternative WSe₂ substrates can be used that provide even higher mobility at room temperature [20].

The fabrication steps of sample-B (see Figure 2a below) are identical, except for the gate electrode: First a trench of SiO₂ is realized with e-beam lithography and dry etching with a carefully controlled depth of 200 nm. Another e-beam lithography step defines a gate electrode within the SiO₂ trench, and gold is evaporated up to a nearly perfect alignment with the SiO₂ surface (an excess of a few nm of gold is preferable to avoid any air gap between the gate and h-BN). The gate dimensions are $L = 0.8$ µm and $W = 22$ µm, and h-BN thickness is 18 nm.

The low contact resistance of sample-A ($R_c \approx 125$ Ωµm) is minimized by contact gating implying that the total resistance is dominated by the channel. Fully ambipolar contact gating was demonstrated in Ref. [19]; in this work the channel length was small so that the transistor was controlled by the contact junction resistance. Here we work with long channels, longer than the optical phonon scattering length, to realize velocity saturation transistors with large voltage gain; contact gating is used instead to minimize the contact resistance. This geometry is well suited for the identification of such mechanisms and their modeling. The sample-B geometry is best suited for maximizing the oscillation frequency $f_{max}$, as dual gates provide a better impedance match of the device and the thick gold allows a reduction of the gate electrode resistance, which is a limiting factor for $f_{max}$. Our devices fulfill the filed-effect mobility requirement $\mu \gtrsim 1$ m²V⁻¹s⁻¹ at low temperature. Experiments are carried out in a Janis cryogenic probe station at 8 K in the band 0–40 GHz. A vector network analyzer (VNA) is used to measure the S-parameters of our transistors.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Low frequency field effect transistors (GFETs). (a) SEM image of sample-A. Note the absence of ungated area. (b) Electrical sketch of samples used in this study. For sample-A we have *L* = 500 nm, *W* = 1.5 µm and *t*<sub>BN</sub> = 16 nm. (c) I–V characteristic of sample A in the p-doping regime. (d) I–V characteristic of sample-A in the n-doping regime. (e) DC resistance of sample-A. (f) RF small signal resistance of sample-A. The global increase of resistance with *V*<sub>ds</sub> is due to HPhP scattering only for negative gate voltages and both HPhP scattering and Dirac pinch-off for positive gate voltages. The resistance peak of *dV*<sub>ds</sub>/*dI* is the best working point for an RF transistor.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Radio-frequency GFETs. (a) Optical image of sample-B before graphene and contact deposition. Dimensions: L = 800 nm, W = 22 µm and t<sub>BN</sub> = 18 nm. (b) Small signal lumped elements circuit of an RF transistor. (c) RF small signal resistance *r* and transconductance *g*<sub>m</sub> of sample-B measured at 2 GHz for biases in the 0–2 V range. (d) Logarithmic plot of the current H(f) and Mason U(f) gains measured between 1 GHz and 40 GHz after de-embedding. (e) *V*<sub>g</sub> dependence of the corresponding measured cut-off frequencies (black circles) together with the model prediction (red lines).</figcaption></figure>

**3. Characterization of Low-Frequency GFETs**

Figure 1c,d shows the I–V characteristic of sample-A in the p-doping and n-doping regime respectively. Both regimes show a trend toward saturation, with larger currents on the p-side due to a lower contact resistance (the contacted graphene is p-doped at ~0.7 × 10¹² cm⁻² by the contact gate kept at a constant voltage here). Figure 1e shows sample-A’s DC resistance for increasing bias voltages *V*<sub>ds</sub>. In the low doping regime, there is a broadening of the resistance peak as *V*<sub>ds</sub> is increased, which we attribute to a thermal effect. More interestingly, in the highly doped regime, we observe an increase of the resistance which is generally attributed to optical phonon scattering. However the relevant quantity is rather the small signal resistance *dV*<sub>ds</sub>/*dI*, shown in Figure 1f, which is directly obtained from the S-parameters and averaged over a 10 MHz to 4 GHz bandwidth. Indeed, on one hand, the small signal resistance is more sensitive to current saturation and on the other hand, it is used to describe the operation of an RF transistor [21]. In addition to the increase of resistance in the doped regime, we observe a pronounced shift of the resistance peak towards positive gate voltages (the resistance peak shifts towards negative *V*<sub>g</sub> when *V*<sub>ds</sub> is negative). This peculiar shape is characteristic of locally gated devices [11] and not observed when using a remote Si gate. This effect, already introduced in Refs. [8,22], is a consequence of finite quantum capacitance effects and geometry which lead to a non-uniform carrier density in the channel at large *V*<sub>ds</sub> and eventually to a Dirac pinch-off. The next section is devoted to the development of a comprehensive model taking into account both mechanisms.
# 4. Simple Model of GFETs

## 4.1. Substrate Phonon Scattering Model

The scattering by optical phonons has been identified as an important phenomenon of high-bias transport in high mobility materials like carbon nanotubes [23]. However, the influence of the substrate cannot be neglected, and the impact of surface polar phonons (SPP) on transport properties has been invoked in several studies to explain the current saturation [8,12,14,15,24]. Those studies have been conducted with graphene on oxide samples, and to our knowledge, the effect of HPhPs has only recently been demonstrated using high-mobility exfoliated graphene on h-BN transistors [11]. It was found that velocity saturation is provided by the lower Ω_I-HPhP band and hot electron cooling by the upper Ω_II-HPhP band. The out-of-plane optical phonons of the Ω_I band lift the sublattice degeneracy of graphene leading to an (unscreened) electron-phonon coupling reminiscent of the mass disorder mechanism reported in Ref. [25]. In Refs. [11,16] exfoliated graphene was used and the mobility was high enough (μ > 3 m²V⁻¹s⁻¹) to access interband transport at very high bias, the so-called Zener-Klein regime, which is not reached in the present study due to limited mobility.

The heuristic model used to describe HPhP scattering is inspired by the optical phonon-electron scattering model used by Yao et al. for carbon nanotubes [23]. One considers that electrons are accelerated by an electric field $E = V_{ds}^0/L$ up to the HPhP energy (here $ħΩ_I \simeq 0.1$ eV) to emit a phonon. In this notation, $V_{ds}^0$ is the voltage across the channel without potential drop at the contacts. This interaction can be expressed as a scattering length for the electron-HPhP interaction $l_{HPhP} = \frac{ħΩ_I}{eE}$, which is density-independent as represented by the dashed lines in Figure 3a for different electric fields. Our samples are not totally ballistic, and we need to account for a finite impurity scattering length $l_{imp} = μħ\sqrt{πn}/e$, shown as a black solid line (the effective mobility is $μ = 1$ m²/V/s). The resulting scattering length is given by the Matthiessen rule $l_{scatt}^{-1} = l_{imp}^{-1} + l_{HPhP}^{-1}$, which is shown at $E = 1$ V µm⁻¹ in Figure 3a as a green solid line. As seen in the figure, at high bias, the transport is indeed dominated by HPhP scattering. The total scattering length $l_{scatt}$ is translated into conductance $G^0 = \frac{4e^2 k_F W}{ħπ} \frac{l_{scatt}}{L+l_{scatt}}$. It can be more conveniently expressed as a resistance:

$$R = \frac{h\pi}{4e^2 k_F W} + \frac{L}{Wne\mu} + \frac{h\pi}{4e^2 k_F W} \frac{eV_{ds}^0}{\hbar\Omega_I} \quad (1)$$

The first term of Equation (1) is the ballistic resistance, the second is the well-known diffusive term, and the third corresponds to the HPhP scattering contribution; it is responsible for saturation and proportional to $V_{ds}^0$. Note that the linear $V_{ds}$ dependence of the DC resistance (negative quadratic deviation in the current) depends only on the carrier density and the HPhP energy. That is how we can identify the $ħΩ_I \simeq 0.1 \pm 0.01$ eV energy of h-BN HPhPs as the main source of scattering and rule out the scattering by intrinsic graphene optical phonons or the Ω_II-HPhP band, as $ħΩ_{OP} \sim ħΩ_{II} \sim 0.2$ eV. Note that the role of intrinsic OPs or Ω_II-HPhPs in energy relaxation can be identified using noise thermometry respectively in suspended [26] and h-BN supported graphene [11,16]. This illustrates the difference between momentum relaxation (the resistance) and energy relaxation (the temperature or electronic distribution) mechanisms (Another example of such a difference, is acoustic phonon scattering which gives rise to a linear temperature dependence for resistivity but a cubic one for relaxation, reflecting single phonon and inelastic supercollision scattering respectively [27]).

In the full saturation regime all terms are negligible compared to the HPhP term and one obtains the saturation current:

$$I_{sat} = \frac{4e^2 k_F W \hbar \Omega_I}{\pi e} \quad (2)$$
Physically, the current saturation is due to the carrier velocity saturation $v_{sat}$:

$$v_{sat} = \frac{I_{sat}}{neW} \approx \frac{2 \hbar \Omega_{I}}{\pi E_{F}} \times v_{F} \quad (3)$$

Equations (1)–(3) strictly hold at large doping when $E_F \gtrsim \hbar\Omega_I$; the crossover to Fermi energy limited current at low doping can be taken into account by substituting $(E_F^{-2} + \hbar\Omega_I^{-2})^{-1/2}$ for $\hbar\Omega_I$ [11] (The $2/\pi$ pre-factor, which stems from the angular projection of electron momentum along the current axis, is sometimes omitted in the literature [8]). In Ref. [11] the saturation current was deduced by fitting the differential conductance by the law $R^{-1}(V_{ds}) = R^{-1}(0)V_{ds}/[1 + (V_{ds}^0/R(0)I_{sat})]^2$, which can be deduced from Equation (1). The saturation is not complete in our samples due to limited mobility and we cannot extract $I_{sat}$ and $v_{sat}$ directly. In Figure 3b, we show the calculated drift velocity $v_d$ as a function of $\mathcal{E}$ (solid lines) for different densities from the current saturation model with $\hbar\Omega_I = 0.1$ eV. Further, we show the measured $v_d$ at large density ($n = 5 \times 10^{12} \text{cm}^{-2}$), where channel density is mostly homogeneous. The very good agreement between theory and experiment, in a regime where HPhP scattering dominates, confirms the relevance of our model.

As a conclusion of this section, we note that the effect of high mobility is mainly to decrease the threshold voltage of saturation and thus to clarify the transition toward the current saturation regime. This is why full current saturation is only observed in high-mobility samples.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** (a) Scattering length associated with the impurity scattering (black solid line) responsible for the low field mobility; here $\mu = 10^4 \text{cm}^2/\text{V/s}$, and with the HPhP scattering (dashed lines) for different in plane electric fields $\mathcal{E} = V_{ds}^0/L$. The green thick line is the total scattering length at the largest bias obtained by the Matthiessen’s law. (b) Carrier drift velocity simulated using $v_d = V_{ds}^0/(RneW)$ with R from Equation (1), and the measured drift velocity at large doping ($n = 5 \times 10^{12} \text{cm}^{-2}$) in black circles.</figcaption></figure>

## 4.2. Dirac Pinch-Off Model

In a remote gate transistor, for example using a global silicon gate, the gate uniformly controls the charge of the semiconducting channel. However, considering a local gate transistor, where local means that the quantum capacitance $C_Q = 2e^2E_F/(\pi\hbar^2v_F^2)$ and the geometrical capacitance per unit area $C_g = \epsilon_0\epsilon_r/t_{BN}$ are on the same order of magnitude, the charge in the channel is not set by the gate voltage only and it can be influenced by the electrostatic potential in the channel. Applying a voltage on the drain electrode results in a linear electrochemical potential along the channel if we assume diffusive conduction, which is the case at large bias due to HPhP scattering. We write the electrochemical balance at the position x in the channel:

$$\Delta E_F(x) - eV_{gr}(x) = -eV_{ds}^0 \times \frac{x}{L} \quad (4)$$
where $V_{gr}(x)$ is the electrostatic potential of graphene. This equation is combined with the electrostatic equilibrium $n(x)e = C_g(V_{gr}(x) - V_g)$, which allows calculating $\Delta E_F(x)$:

$$ \Delta E_F = \mathrm{sign}(\epsilon_W) \times \epsilon_c \left( \sqrt{1 + 2\frac{|\epsilon_W|}{\epsilon_c}} - 1 \right), \quad (5) $$

with $\epsilon_c = C_g \frac{\pi\hbar^2 v_F^2}{2e^2}$ and $\epsilon_W(x) = eV_g - eV_{ds}^0 \times \frac{x}{L}$. It results in an inhomogeneous $\Delta E_F$ within the channel. This effect is sketched in Figure 4a, and the calculated $\Delta E_F(x)$ is shown in Figure 4b for increasing $V_{ds}^0$. The corresponding carrier density in the channel as a function of $V_g$ and $x$ is shown in the color plot of Figure 4c. Starting from an n-doped channel and increasing $V_{ds}^0$ leads to a decrease of the carrier density on the drain side. Figure 4b shows that the neutrality can be reached and that for larger $V_{ds}^0$ there is even a localized neutral point whose location is shifted toward the source electrode. At the onset of neutrality on the drain side, increasing $V_{ds}^0$ induces an increase of the resistance, which amplifies the sublinear behavior of the I–V characteristic and results in enhanced current saturation. This is the Dirac pinch-off. On the contrary, starting from a p-doped channel, an increase of $V_{ds}^0$ leads to an increase of the carrier density on the drain side which contributes to decreasing the channel resistance. This explains the observed asymmetric shape of the high bias resistance curve in Figure 1f.

<figure><img src="image_7.png" /><figcaption>Figure 4. (a) Sketch of the doping profile across the graphene channel when $V_{ds}$ is increased. Starting from n-doping, $\Delta E_F(x)$ can be decreased down to neutrality at the drain side. This is the Dirac pinch-off. (b) Simulated $E_F(x)$ along the channel for increasing in-plane electric field. (c) Color plot of the local carrier density in the channel as a function of the x-position along the channel and the gate voltage $V_g$ (d) Color plot of the local resistivity $r_{ds}(V_g, x)$ at $V_{ds}^0 = 0.5$ V. The total resistance of the device is the integration over x of this resistivity.</figcaption></figure>

### 4.3. Hot Electron Effects

To accurately describe the behavior of our transistor, we have to take into account the thermal broadening that arises at large $V_{ds}^0$. Indeed, for high biased samples, the electronic temperature can
rise up to 800 K [28], which significantly modifies the energy distribution. We take into account this thermal broadening by writing the carrier density with a residual term $n_0$ that depends on $V_{ds}$:

$$n(x) = \sqrt{\left(\frac{\Delta E_F(x)^2}{\pi \hbar^2 v_F^2}\right)^2 + n_0(V_{ds}^0)^2} \quad (6)$$

where $n_0(V_{ds}^0)$ is the sum of an impurity term ($n_{imp} = 2 \times 10^{11} \text{cm}^{-2}$ at zero bias) and a hot electron term: $n_{th}(V_{ds}^0) \sim \frac{\Delta(k_B T_e)^2}{\pi \hbar^2 v_F^2}$. Here, $T_e$ is governed by the interaction between electrons and acoustic phonons of graphene. It has been shown in previous studies [27] that, at low doping and moderate impurity scattering, electron cooling is controlled by electron-phonon-impurity supercollisions and that for a given Joule power input one has $V_{ds}^2/(RLW) = AT_e^3$ (In very high mobility h-BN supported/encapsulated samples, supercollisions are suppressed and the leading relaxation mechanism becomes black-body and luminescent $\Omega_{II}$-HPhP emission). In this equation we have $A \propto n$ and $R \propto 1/n$ so that AR is constant. Finally, the hot electron residual density reads as follows,

$$n_{th} = \alpha(V_{ds}^0)^{4/3}, \quad (7)$$

where $\alpha \approx 9 \times 10^3/(WL)^{2/3} \text{cm}^{-2}\text{V}^{-4/3}$ is a constant that depends on sample dimensions. This is used to calculate the local resistance $r_{ds}(V_g, x)$, which is plotted in Figure 4d. for $V_{ds}^0 = 0.5$ V. As expected, the maximum resistance line follows the minimum carrier density line, with a width related to $V_{ds}^0$.

## 4.4. Complete Model and Discussion

Let us now calculate the total resistance of the graphene channel taking into account HPhP scattering, the non-uniform carrier density in the channel and the hot electron effect. We compute the total scattering length $l_{scatt}(x)^{-1} = l_{imp}(V_g, x)^{-1} + l_{HPhP}(x)^{-1}$ and neglect the ballistic conductance to obtain the local resistance

$$r_{ds}(x)dx = \frac{\pi h}{4e^2 W \sqrt{n(x)\pi}} \frac{dx}{l_{scatt}(x)} \quad (8)$$

The total channel resistance is then obtained by integrating over the channel length:

$$R_{ds} = \int_{0}^{L} r_{ds}(x)dx \quad (9)$$

Finally, in short devices with high-mobility graphene, the contact resistance $R_c$ is not negligible. Most of the contact resistance comes from the junction between contacted graphene and the channel. This issue has been overcome in a previous study with the transistor geometry of sample-A [19]. It allows us to calculate the $V_g$ dependent contact resistance. The total transistor resistance is then $R(V_g, V_{ds}) = R_{ds}(V_g, V_{ds}) + 2R_c(V_g)$. For sample-A, we have $R_c \approx 125 \text{ } \mu\text{m}$ on the p-doping side (where there are p-p junctions at the two contact edges) and $R_c \approx 275 \text{ } \mu\text{m}$ on the n-doping side (where there are p-n junctions at the two contact edges).

The relevant quantity for an RF transistor is the small signal resistance $dV_{ds}/dI_{ds}$ that we calculate from the current through the transistor $I_{ds} = V_{ds}/R$ and the applied bias voltage $V_{ds}$ across the device. The simulated $dV_{ds}/dI_{ds}$ is plotted in Figure 5a where the mobility $\mu = 10^4 \text{ cm}^2\text{V}^{-1}\text{s}^{-1}$ and $R_c$ are determined from the zero bias transfer curve and the thermal broadening coefficient $\alpha$ (from [28]) is in principle universal. The very good agreement between measurement and simulation shows that we are able to accurately predict the behavior of a local gate transistor driven at high bias knowing the graphene mobility, the contact resistance and the substrate phonon energy. It validates our transistor model and confirms the role of the h-BN HPhPs at $\hbar\Omega = 0.1$ eV in velocity saturation. To emphasize the effect of the local gate on the saturation we plot in Figure 5b, a simulation for a substrate thickness of $t_{BN} = 300$ nm where all parameters are conserved (apart from $V_g$ which is adjusted to conserve the
zero bias density). There is a significant difference between local and remote gate on the shape of the transfer curve. In Figure 5b, the slight asymmetry between p and n doping comes from the contact resistance difference.

<figure><img src="image_8.png" /><figcaption>Figure 5. (a) Small signal resistance of the transistor obtained from the model (solid lines) and from the measured scattering parameters of sample-A at 2 GHz (circles). (b) Small signal resistance simulated for a remote gate transistor (h-BN dielectric thickness of 300 nm). The saturation effect is enhanced for the local gate device in the n-doping regime.</figcaption></figure>

## 5. Characterization of RF-GFETs

The performance of Sample-A as an RF transistor is assessed at high frequency. Measuring the S-parameters gives a direct access to the channel small signal resistance $r$ and the transconductance $g_m = dI_{ds}/dV_g$ shown in the schematic of Figure 2b. We also measure the current gain $H$ and the Mason gain $U$ (equivalent to the power gain). Those gains are sensitive to extrinsic parameters such as stray capacitances and, for $U$ only, the resistance of the gate electrode. After removing the effect of stray capacitances by a de-embedding procedure, we obtain the current gain cut-off frequency $f_T = 9$ GHz and the Mason’s gain cut-off frequency $f_{max} = 7$ GHz. The fact that $f_{max} \sim f_T$ is a distinctive feature of velocity saturation GFETs, as $f_{max} \ll f_T$ in mobility-limited GFETs [29]. The value of $f_{max}$ is limited here by the large resistance $R_G \sim 10$ Ohms of the tungsten gate electrode. This is why we built a thick gold gate device in sample-B (shown in Figure 2a).

In sample-B, the channel width is also further increased to $W = 22$ µm, to increase the total current in the device. In terms of $f_{max}$, however, a compromise has to be found between a large current which means a large transconductance and a larger gate resistance which strongly limits $f_{max}$. The small signal characteristics of sample-B, i.e., the resistance $r$ and transconductance $g_m$, are plotted in Figure 2c. One can notice the strong saturation effect on the n-side where the high bias $r$ is larger than the zero bias resistance peak. In this sample, the n-side is also the low contact resistance side which results in a larger $g_m$ on the n-side. Both effects contribute to a larger voltage gain for positive gate voltages. This is confirmed by measuring $H$ and $U$ at the maximum bias and $V_g = 2$ V plotted in Figure 2d after de-embedding. In Figure 2e, we plot the $V_g$-dependence of $f_T$ and $f_{max}$ at $V_{ds} = 2.4$ V. Our model allows estimating $r$ and $g_m$ which is useful to determine the transistor’s figures of merit from $f_T = g_m/(2\pi C_g)$ and

$$f_{max} = \frac{f_T}{2\sqrt{(R_G + R_c)/r + 2\pi f_T C_g R_G}} \quad (10)$$

The predicted $f_T$ and $f_{max}$ in Figure 2e compare fairly well with our measurement. Both figures of merit present two peaks for p- and n-doping. In the case of $f_T$, the peak asymmetry comes from the
contact resistance asymmetry. For $f_{max}$ the peak asymmetry is stronger because it is also sensitive to a stronger saturation on the n-doping side.

The availability of a reliable model allows fine tuning parameters to optimize $f_{max}$. It appears that in our sample the large gate resistance is the main limitation. Engineering the transistor extrinsic elements can improve $f_{max}$ by a large factor, for example, by using a T-shaped gate [9,30]. Another important limiting factor is the contact resistance.

## 6. Perspectives on High-Frequency Graphene Devices

In an attempt to assess the intrinsic frequency limits of GFETs, we include here an analysis of Zener-Klein transistors where contact resistance is low enough, mobility is large enough and channel length long enough ($L = 4$ µm), to access a fully developed inter-band Zener-Klein tunneling regime [11]. A quasi perfect current saturation ($I_{sat} \gtrsim 20$ Acm$^{-1}$) is achieved thanks to the combined intraband velocity saturation, interband Zener-Klein tunneling and Dirac pinch-off (Figure 1 of Ref. [11]). The small contact resistance $R_c \simeq 150$ Ωµm$^{-1}$, large DC transconductance $g_m \gtrsim 2.5$ Scm$^{-1}$, and saturation resistance $R_{sat}$ lead to a large voltage gain $G \simeq 8$. From these numbers, we estimate $f_{max} \approx f_T \sqrt{G/4g_m R_c} \simeq 7f_T$ (taking $R_G \ll R_c$ in Equation (10)) translating to $f_{max} \gtrsim 100$ GHz (taking $f_T = g_m/2\pi C_g \simeq 16$ GHz with $L = 2$ µm and $t_{BN} = 20$ nm). Considering the optimization of saturation velocity, contact resistance and channel length, these figures constitute the intrinsic frequency limit of GFETs. It still remains to be demonstrated with edge-contacted, h-BN encapsulated, GFETs at room temperature.

This intrinsic limitation of GFETs can be circumvented relying on collective modes such as plasmons which benefit from an enhanced velocity $v_p \gg v_F$ appropriate for sub-THz operation. This corresponds to plasma resonance transistors (PRTs) [18], which are multi-mode narrow-band detectors. Intrinsic Dirac-fermion plasma modes can be studied in plasma resonance capacitors (PRCs), which behave as quarter-mode resonators (fundamental at $f_0 = v_p/4L$ and harmonics $f_k = (2k+1)f_0$). In a recent study [17], an h-BN encapsulated PRC working at $f_0 = 40$ GHz ($v_p \simeq 4v_F$) was demonstrated, showing the relevance of long-wavelength plasma resonance devices in fully characterizing high-mobility graphene, offering a full analysis of the temperature and density dependence of the quantum capacitance $C_Q(T,n)$, kinetic inductance $L_K(T,n)$ and damping or mean-free path $l_{scat}(T,n)$.

## 7. Conclusions

We have fabricated, modeled, and measured high-performance velocity-saturation GFETs showing promising perspectives for microwave detection and amplification. Our modeling involves physical insight into microscopic mechanisms of velocity saturation and a proper account of Dirac pinch-off effects. Our work allows assessing the intrinsic frequency limits of GFETs using velocity saturation associated with Zener-Klein tunneling and opens perspectives “beyond-GFETs”, using plasma resonance techniques for sub-THz detection. The figure of merit is the $f_{max}/f_T$ ratio which increases from $f_{max}/f_T \lesssim 0.2$ in the mobility-limited regime, over $f_{max}/f_T \sim 1$ in the velocity-saturation regime, up to $f_{max}/f_T \gtrsim 5$ in the Zener-Klein regime. The latter offers promising perspectives in RADAR and GSM applications.

**Author Contributions:** Q.W., M.B., E.P. and B.P. conceived the experiment. T.T. and K.W. provided the hBN crystals. V.B. provided the CVD graphene. Q.W., M.B., H.G., J.P. and M.R. did the nanofabrication. Q.W., M.B., H.G. and B.P. did the electrical measurements. Q.W., M.B., H.G., D.M., E.B. (Emmanuel Baudin), J.-M.B., G.F., E.B. (Erwann Bocquillon) and B.P. participated to the data analysis and theoretical modelling. Q.W. and B.P. wrote the manuscript with contributions from all coauthors. All authors have read and agreed to the published version of the manuscript.

**Funding:** The research leading to these results have received partial funding from the European Union Horizon 2020 research and innovation programme under grant No. 785219 "Graphene Core 2", and from the ANR-14-CE08-018-05 "GoBN". The work of Q.W. was supported by a DGA-MRIS scholarship.

**Acknowledgments:** We thank L. Becerra for providing the tungsten films.
**Conflicts of Interest:** The authors declare no conflicts of interest.

**References**

1. 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]

2. Wang, L.; Meric, I.; Huang, P.Y.; Gao, Q.; Gao, Y.; Tran, H.; Taniguchi, T.; Watanabe, K.; Campos, L.M.; Muller, D.A.; et al. One-Dimensional Electrical Contact to a Two-Dimensional Material. *Science* **2013**, *342*.
[CrossRef] [PubMed]

3. Han, Z.; Kimouche, A.; Kalita, D.; Allain, A.; Arjmandi-Tash, H.; Reserbat-Plantey, A.; Marty, L.; Pairis, S.; Reita, V.; Bendiab, N.; et al. Homogeneous Optical and Electronic Properties of Graphene Due to the Suppression of Multilayer Patches During CVD on Copper Foils. *Adv. Funct. Mater.* **2014**, *24*, 964–970.
[CrossRef]

4. Banszerus, L.; Schmitz, M.; Engels, S.; Dauber, J.; Oellers, M.; Haupt, F.; Watanabe, K.; Taniguchi, T.; Beschoten, B.; Stampfer, C. Ultrahigh-mobility graphene devices from chemical vapor deposition on reusable copper. *Sci. Adv.* **2015**, *1*, 2375–2548.
[CrossRef] [PubMed]

5. Banszerus, L.; Schmitz, M.; Engels, S.; Goldsche, M.; Watanabe, K.; Taniguchi, T.; Beschoten, B.; Stampfer, C. Ballistic Transport Exceeding 28 m in CVD Grown Graphene. *Nano Lett.* **2016**, *16*, 1387. [CrossRef]

6. Meng, N.; Fernandez, J.F.; Vignaud, D.; Dambrine, G.; Happy, H. Fabrication and Characterization of an Epitaxial Graphene Nanoribbon-Based Field-Effect Transistor. *IEEE Trans. Electron Dev.* **2016**, *58*, 1594.
[CrossRef]

7. Sze, S.M.; Ng, K. *Physics of Semiconductor Devices*, 3rd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2007; p. 293.

8. Meric, I.; Han, M.Y.; Young, A.F.; Ozyilmaz, B.; Kim, P.; Shepard, K.L. Current saturation in zero-bandgap, top-gated graphene field-effect transistors. *Nat. Nanotechnol.* **2008**, *3*, 654–659. [CrossRef]

9. 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]

10. Feng, Z.H.; Yu, C.; Li, J.; Liu, Q.B.; He, Z.Z.; Song, X.B.; Wang, J.J.; Cai, S.J. An ultra clean self-aligned process for high maximum oscillation frequency graphene transistors. *Carbon* **2014**, *75*, 249–254. [CrossRef]

11. Yang, W.; Berthou, S.; Lu, X.; Wilmart, Q.; Denis, A.; Rosticher, M.; Taniguchi, T.; Watanabe, K.; Fève, G.; Berroir, J.M.; et al. A graphene Zener–Klein transistor cooled by a hyperbolic substrate. *Nat. Nanotechnol.* **2018**, *13*, 47. [CrossRef]

12. Barreiro, A.; Lazzeri, M.; Moser, J.; Mauri, F.; Bachtold, A. Transport Properties of Graphene in the High-Current Limit. *Phys. Rev. B* **2009**, *70*, 076601. [CrossRef] [PubMed]

13. Fratini, S.; Guinea, F. Substrate-limited electron dynamics in graphene. *Phys. Rev. B* **2008**, *77*, 195415.
[CrossRef]

14. Freitag, M.; Steiner, M.; Martin, Y.; Perebeinos, V.; Chen, Z.; Tsang, J.C.; Avouris, P. Energy Dissipation in Graphene Field-Effect Transistors. *Nano Lett.* **2009**, *9*, 1883. [CrossRef] [PubMed]

15. Perebeinos, V.; Avouris, P. Inelastic scattering and current saturation in graphene. *Phys. Rev. B* **2010**, *81*, 5442.
[CrossRef]

16. Baudin, E.; Voisin, C.; Plaçais, B. Hyperbolic Phonon Polariton Electroluminescence as an Electronic Cooling Pathway. *Adv. Funct. Mater.* **2019**. [CrossRef]

17. Graef, H.; Mele, D.; Rosticher, M.; Stampfer, C.; Taniguchi, T.; Watanabe, K.; Bocquillon, E.; Fève, G.; Berroir, J.M.; Teo, E.T.H.; et al. Ultra-long wavelength Dirac plasmons in graphene capacitors. *emph. J. Phys. Mater.* **2018**, *1*, 01LT02. [CrossRef]

18. Bandurin, D.A.; Svintsov, D.; Gayduchenko, I.; Xu, S.G.; Principi, A.; Moskotin, M.; Tretyakov, I.; Yagodkin, D.; Zhukov, S.; Taniguchi, T.; et al. Resonant terahertz detection using graphene plasmons. *Nat. Commun.* **2018**, *9*, 5392. [CrossRef]

19. Wilmart, Q.; Inhofer, A.; Boukhicha, M.; Yang, W.; Rosticher, M.; Morfin, P.; Garroum, N.; Fève, G.; Berroir, J.-M.; Placais, B. Contact gating at GHz frequency in graphene. *Sci. Rep.* **2016**, *6*, 21085. [CrossRef]
20. Banszerus, L.; Sohier, T.; Epping, A.; Winkler, F.; Libisch, F.; Haupt, F.; Watanabe, K.; Taniguchi, T.; Müller-Caspary, K.; Marzari, N.; et al. Extraordinary high room-temperature carrier mobility in graphene-WSe2 heterostructures. *arXiv* **2019**, arXiv:1909.09523v1.

21. Schwierz, F. Graphene transistors. *Nat. Nanotechnol.* **2010**, *5*, 487. [CrossRef]

22. Bai, J.; Liao, L.; Zhou, H.; Cheng, R.; Liu, L.; Huang, Y.; Duan, X. Top-Gated Chemical Vapor Deposition Grown Graphene Transistors with Current Saturation. *Nano Lett.* **2011**, *11*, 2555–2559. [CrossRef] [PubMed]

23. Yao, Z.; Kane, C.L.; Dekker, C. High-Field Electrical Transport in Single-Wall Carbon Nanotubes. *Phys. Rev. Lett.* **2000**, *84*, 2941–2944. [CrossRef] [PubMed]

24. Dorgan, V.E.; Bae, M.-H.; Pop, E. Mobility and saturation velocity in graphene on SiO₂. *Appl. Phys. Lett.* **2010**, *97*, 082112. [CrossRef]

25. Pallecchi, E.; Betz, A.C.; Chaste, J.; Fève, G.; Huard, B.; Kontos, T.; Berroir, J.-M.; Plaçais, B. Transport scattering time probed through rf admittance of a graphene capacitor. *Phys. Rev. B* **2011**, *83*, 125408. [CrossRef]

26. Laitinen, A.; Kumar, M.; Oksanen, M.; Plaçais, B.; Virtanen, P.; Hakonen, P. Coupling between electrons and optical phonons in suspended bilayer graphene. *Phys. Rev. B* **2015**, *91*, 121414. [CrossRef]

27. Betz, A.C.; Jhang, S.H.; Pallecchi, E.; Ferreira, R.; Fève, G.; Berroir, J.-M.; Plaçais, B. Supercollision cooling in undoped graphene. *Nat. Phys.* **2013**, *9*, 109–112. [CrossRef]

28. Betz, A.C.; Vialla, F.; Brunel, D.; Voisin, C.; Picher, M.; Cavanna, A.; Madouri, A.; Fève, G.; Berroir, J.-M.; Plaçais, B.; et al. Hot electron cooling by acoustic phonons in graphene. *Phys. Rev. Lett.* **2012**, *109*, 056805. [CrossRef]

29. Pallecchi, E.; Benz, C.; Betz, A.C.; Löhneysen, H.v.; Plaçais, B.; Danneau, R. Graphene microwave transistors on sapphire substrates. *Appl. Phys. Lett.* **2011**, *99*, 113502. [CrossRef]

30. Mele, D.; Mehdhibi, S.; Fadil, D.; Wei, W.; Ouerghi, A.; Lepilliet, S.; Happy, H.; Pallecchi, E. Graphene FETs Based on High Resolution Nanoribbons for HF Low Power Applications. *Electron. Mater. Lett.* **2018**, *14*, 133. [CrossRef]

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

© 2020 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 (http://creativecommons.org/licenses/by/4.0/).