Next Article in Journal
Friction Determination by Atomic Force Microscopy in Field of Biochemical Science
Previous Article in Journal
Design and Performance Test of an Ocean Turbulent Kinetic Energy Dissipation Rate Measurement Probe
Previous Article in Special Issue
Hard-Baked Photoresist as a Sacrificial Layer for Sub-180 °C Surface Micromachining Processes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Frequency Tuning of Graphene Nanoelectromechanical Resonators via Electrostatic Gating

1
School of Electronic Science and Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China
2
Electrical Engineering, Case School of Engineering, Case Western Reserve University, Cleveland, OH 44106, USA
*
Authors to whom correspondence should be addressed.
Micromachines 2018, 9(6), 312; https://doi.org/10.3390/mi9060312
Submission received: 16 April 2018 / Revised: 11 May 2018 / Accepted: 31 May 2018 / Published: 20 June 2018
(This article belongs to the Special Issue Carbon Based Materials for MEMS/NEMS)

Abstract

:
In this article, we report on a comprehensive modeling study of frequency tuning of graphene resonant nanoelectromechanical systems (NEMS) via electrostatic coupling forces induced by controlling the voltage of a capacitive gate. The model applies to both doubly clamped graphene membranes and circumference-clamped circular drumhead device structures. Frequency tuning of these devices can be predicted by considering both capacitive softening and elastic stiffening. It is shown that the built-in strain in the device strongly dictates the frequency tuning behavior and tuning range. We also find that doubly clamped graphene resonators can have a wider frequency tuning range, while circular drumhead devices have higher initial resonance frequency with same device characteristic parameters. Further, the parametric study in this work clearly shows that a smaller built-in strain, smaller depth of air gap or cavity, and larger device size or characteristic length (e.g., length for doubly clamped devices, and diameter for circular drumheads) help achieve a wider range of electrostatic frequency tunability. This study builds a solid foundation that can offer important device fabrication and design guidelines for achieving radio frequency components (e.g., voltage controlled oscillators and filters) with the desired frequencies and tuning ranges.

1. Introduction

Two-dimensional (2D) materials such as graphene and atomic layer semiconductors have recently attracted tremendous attention and research interest due to their unconventional and extraordinary material properties stemming from their atomically thin layered structures. In addition to excellent and unique properties in optical [1] and electrical domains (e.g., mobility of 1,000,000 cm2/(V⋅s) [2]), graphene also exhibits remarkable mechanical properties such as ultrahigh Young’s modulus EY = 1 TPa, and a large breaking strain limit up to 25% [3,4]. These attributes have made graphene an attractive and promising candidate for highly miniaturized and aggressively scaled resonant-mode nanoelectromechanical systems (NEMS) with unprecedented device performance. To date, mainly doubly clamped [5,6,7,8] membranes and ribbons, and circumference-clamped circular drumhead graphene resonators [9,10,11,12,13,14,15,16] have been prototyped, and fundamental device physics and basic device characteristics have been studied. Further, potential applications of these graphene 2D NEMS resonators in sensing of external stimuli and perturbations [17], and components for radio frequency (RF) signal processing and communications (e.g., oscillators [11]), have been attempted. For these applications, continuous and wide frequency tuning controlled by moderate level electrical signals (i.e., voltage or current) is highly desirable, to render these systems tunable, flexible, or even programmable and reconfigurable. Indeed, strong frequency tunability has been quite heavily pursued in more conventional and state-of-the-art microelectromechanical systems (MEMS) based resonators and oscillators; however, tuning range is often limited up to 5% due to their high stiffness [18]. Thanks to ultra-strong yet highly stretchable crystals and their related material properties, graphene NEMS resonators can exhibit remarkably broad tunability of resonance frequency, with Δfres/fres > 300% [12].
Among various device platforms and resonance excitation and detection methods, including photothermal and electrothermal schemes, electrostatic excitation and control of graphene 2D NEMS devices are always attractive for on-chip integration with mainstream technologies, and are particularly promising toward achieving wide frequency tuning ranges [6,7,8,9,10,11,12]. In the electrostatic scheme, device structures form a capacitor between a suspended 2D material and a bottom (or top) gate, enabling electrostatic control of electrical and mechanical device performance by applying direct current (DC) electric potential (i.e., DC gate voltage). The tunable conductance with respect to the applied DC potential at the gate offers strong coupling between the mechanical motion and the electrical conductance, enabling mechanical motion detection via sensitive vibrating channel field effect transistor (VCFET) signal transduction [6]. Further, the electrostatic scheme only requires ultralow power consumption owing to lack of DC current flowing between the suspended graphene NEMS and the gate. More importantly, in the electromechanical domain, frequency tuning can be easily achieved by applying the gate DC voltage, which modifies the effective spring constants of the devices by electrostatic-mechanical coupling. The reported frequency tuning results in the electrostatic scheme to date, however, exhibit complicated behavior including frequency downshifts and upshifts [6,7,8,9,10,11,12,13,14,16]. Several attempts have been made to develop frequency tuning models that can be used to describe and explain the precise tuning mechanisms of devices [6,9,10,19,20,21,22,23,24]. So far, however, these models are not sufficient to comprehensively explain the complicated frequency tuning behavior and the device coupling mechanisms. Thus, it can be quite challenging to understand contributions from various parameters, including device geometry and built-in tension (strain) on the frequency tuning ranges.
In this work, we present an accurate frequency tuning model for electrostatically actuated doubly clamped membranes and circumference-clamped circular drumhead graphene resonators. This model can quantitatively clarify various mechanisms of frequency tuning. We carefully develop the frequency tuning model via electrostatic gating voltage by considering multiphysics coupling of elastic and electrostatic effects. Our analysis reveals that various device parameters strongly govern the device resonance and its frequency tuning range. Moreover, we find that frequency tuning range and wide-range tuning behavior can be finely engineered by controlling important parameters such as the built-in strain.

2. Analytical Model and Computational Methods

2.1. Development of Frequency Tuning Model and Analysis Procedure

We start to develop our model by calculating the deflection of the suspended graphene as a thin membrane subject to external forces. We assume that the suspended 2D graphene membrane has negligible flexural rigidity, and it can be stretchable by transversely pulling the membrane using the electrostatic force induced by a gate voltage. Figure 1a shows the side-view schematic of the deflection of the graphene resonator, and Figure 1b,c illustrate the schematics of a doubly clamped membrane and the circumference-clamped circular drumhead device structure, respectively. The deflection of the 2D material is precisely computed by following the procedure outlined in the flowchart shown in Figure 1d. The equilibrium displacement (ze) on the center of the devices is solved by minimizing the total energy (sum of elastic energy Uel and electrostatic energy Ues) upon the application of electrostatic force. To make the model accurate, we perform several iterative calculations of the capacitance between the gate and suspended membrane Cg, equilibrium displacement ze, and total strain ε after deflection. Substituting the modified strain ε and Cg into the effective spring constant keff, the resonance frequency fres can then be obtained using a relation of fres = (keff/meff)1/2/2π, with the calculated effective mass of resonance mode, meff.

2.2. Frequency Tuning of Doubly Clamped Graphene Resonator

Frequency tuning in a doubly clamped graphene resonator under electrostatic coupling can be calculated by examining the elastic energy of the suspended graphene and electrostatic energy in the device capacitor. Consider the doubly clamped device has a static deflection profile of ud(x) = 4ze(Lxx2)/L2, (0 < x < L), where L is the length of the suspended graphene membrane (or ribbon), and it has the maximum static deflection of ze at the midpoint (x = L/2) due to the center of the symmetric structure. Under uniformly distributed force, we neglect the Poisson’s ratio in doubly clamped case since the resonance frequency is independent on width of the membrane. For the doubly clamped device structure shown in Figure 1b, the elastic energy stored in the stretched membrane Uel,d can be described as [24]
U el , d = 0 w 0 L γ 2 u 2 ( x ) d x d y = 0 w 0 L u 2 ( x ) [ E Y t ε 0 2 + E Y t 4 L 0 L u 2 ( χ ) d χ ] d x d y
where w, EY, t, and ε0 are the width, Young’s modulus, thickness, and built-in strain of membrane, respectively. The total tension after deflection is γ = EYd, where εd is the total strain that is composed of both the built-in strain and the added strain according to the deflection. By introducing deflection profile ud(x) into Equation (1), the elastic energy is then
U el , d = 64 E Y t w z e 4 9 L 3 + 8 E Y t w ε 0 z e 2 3 L .
We consider the suspended graphene as a capacitor where the graphene is a deformable electrode, the bottom back gate is a fixed plate electrode, and the air gap is the dielectric layer. Considering deflection profile, the device capacitance Cg,d is
C g , d = 0 w 0 L є 0 z 0 z e 4 L 2 ( L x x 2 ) d x d y ,
where є 0 is the permittivity, and z0 is the depth of the air gap. The electrostatic energy stored in the device structure is then
U es = 1 2 C g , d V g 2 ,
where Vg is the gate voltage. The equilibrium displacement ze can be obtained by minimizing the total energy thus finding the deflection position where the elastic and electrostatic forces are equal:
( U el , d + U es ) z e = 256 E Y t w z e 3 9 L 3 + 16 E Y w t ε 0 z e 3 L 1 2 C g V g 2 = 0 .
Although Equation (3) can be further simplified by expanding it to Cg,d = є 0 wL/z0 + 2 є 0 wLze/3 z 0 2 + 8є0wLze2/15 z 0 3 + …, and the initial three terms of the series can be used for calculating ze using Equation (5), this calculation produces an approximated value of ze and it may give quite large error, especially when we calculate the exponent of ze. Instead, here we perform several iterations to find accurate device Cg by repeating a process of plugging ze obtained using Equation (5) back into Equation (3) till it is convergent. The total strain εd in the graphene includes built-in strain and the added strain upon deformation. The added strain after deflection is obtained by comparing the extended length from the deformed graphene sheet and the initial graphene length. The total strain is
ε d = 1 + ε 0 L 0 L 1 + ( u x ) 2 d x 1 = 1 + ε 0 2 1 + 16 z e 2 L 2 + ( 1 + ε 0 ) L 8 z e ln ( 1 + 16 z e 2 L 2 + 4 z e L ) 1 .
Next, we consider the dynamic motion of the doubly clamped graphene membrane vibrating at its resonance frequency. The mode shape of the fundamental mode resonance of doubly clamped devices can be expressed by a sinusoidal function. Considering the clamped boundary conditions, we obtain the mode shape of fundamental resonance uf,d(x) = δzsin(πx/L), where δz is the dynamic displacement at the midpoint. We assume that at the small static defection ze of the graphene membrane, the vibration mode shape remains to be in the sinusoidal form. The elastic energy δUel for the fundamental mode resonance of the strained graphene resonator is
δ U el , d = 0 w 0 L γ 2 u f , d 2 ( x ) d x d y = π 2 E Y t w ε d 4 L ( δ z ) 2 .
The effective spring constant keff,d can be given by the second order differentiation of the total energy:
k eff , d = 2 δ U el , d δ z 2 + 2 U es z e 2 = π 2 E Y w t ε d 2 L 8 є 0 w L 15 z 0 3 V g 2 .
Resonance frequency is determined by fres = (1/2π)(keff/meff)1/2, where meff is effective mass which can be calculated from kinetic energy of membrane at resonance. The peak kinetic energy Ekin,d is
E kin , d = 1 2 m eff , d ( δ z ) 2 = 1 2 ρ t 0 w 0 L ( δ z sin ( π L x ) ) 2 d x d y = 1 4 ρ t w L ( δ z ) 2 .
From Equation (9), the effective mass of the fundamental mode is meff,d = 0.5ρtwL in doubly clamped case, where ρ is the mass density. The frequency tuning for the doubly clamped resonator can be given by
f res , d = 1 2 π k eff , d m eff , d = 1 2 π π 2 E Y ε d ρ L 2 16 є 0 15 ρ t z 0 3 V g 2 .

2.3. Circumference-Clamped Circular Membrane

We now turn to analyze the frequency tuning behavior of the circular drumhead graphene resonator. For convenience, we use a polar coordinate for the circular drumhead that closely matches the device geometry. The elastic energy stored in the stretched membrane Uel,c can be obtained by [19,25]
U el , c = 0 2 π 0 R { E Y t 2 [ ε 0 , r ( u r ) 2 + ε 0 , θ r 2 ( u θ ) 2 ] + E Y t 8 ( 1 v 2 ) ( 2 u r 2 ) 2 } r d r d θ ,
where R, ν, ε0,r, ε0,θ are radius, Poisson’s ratio, initial radial strain, initial tangential strain, respectively. Similar to the doubly clamped case, we assume that the curvature of the static deflection forms the parabolic shape, and it has the maximum static deflection at its center due to symmetry. The solution of deflection can also be uc(r) = ze(R2r2)/R2, where ze is the static defection at the center of the drumhead. The elastic energy induced by electrostatic deformation is
U el , c = 2 π E Y t z e 4 3 ( 1 v 2 ) R 2 + π E Y t ε 0 , r z e 2 + π E Y t ε 0 , r R 2 2 ( 1 v 2 ) .
The device capacitance between the suspended graphene and the bottom gate is
C g , c = 0 2 π 0 R є 0 z 0 z e 1 R 2 ( R 2 r 2 ) r d r d θ ,
and the equilibrium displacement ze can be obtained by minimizing the total energy
( U el , c + U es ) z e = 8 π E Y t z e 3 3 ( 1 v 2 ) R 2 + 2 π E Y t ε 0 , r z e 1 2 C g V g 2 = 0 .
To find accurate ze, we conduct iterations using Equations (12)–(14) till it is convergent. Based on the deflection curvature and estimated ze, the total radial strain from stretching of the drumhead is estimated by calculating the radial elongation of the membrane:
ε r = 1 + ε 0 , r 2 R R R 1 + ( u r ) 2 d r 1 = 1 + ε 0 , r 2 1 + 4 z e 2 R 2 + ( 1 + ε 0 , r ) R 4 z e 2 ln ( 1 + 4 z e 2 R 2 + 2 z e R ) 1 .
The fundamental mode resonance and its mode shape of circular membrane can be expressed using the 0th-order Bessel function of the first kind, J0. Considering the circumference-clamped boundary conditions, the mode shape of the fundamental resonance is uf,c(r) = δzJ0(2.405r/R). The elastic energy of the strained circular drumhead at resonance is then modified by substituting the mode shape uf,c(r) into Equation (11):
δ U el , c = 0.271 2.405 2 π E Y t ε r 2 ( δ z ) 2 + π E Y t ε r R 2 2 ( 1 v 2 ) .
With capacitive softening, the effective spring constant for the circular drumhead is given by
k eff , c = 2 δ U el , c δ z 2 + 2 U es z e 2 = 4.924 E Y t ε r є 0 π R 2 3 z 0 3 V g 2 .
The effective mass of the device can be calculated from the peak kinetic energy,
E kin , c = 1 2 m eff , c ( δ z ) 2 = 1 2 ρ t 0 2 π 0 R ( δ z J 0 ( 2.4 R r ) ) 2 r d r d θ = 0.136 π ρ t R 2 ( δ z ) 2 .
We obtain the effective mass meff = 0.271ρtπR2 for the fundamental mode resonance in the circular drumhead case. Thus, the resonance frequency for the circular drumhead under electrostatic gating is given by
f res , c = 1 2 π k eff , c m eff , c = 1 2 π 2.4 2 E Y ε r ρ R 2 є 0 0.813 ρ t z 0 3 V g 2 .

3. Results and Discussions

We first calculate frequency scaling of the single layer (1L) graphene resonators without the external gate voltage to understand effects of device dimension and built-in strain on the resonance frequency. Figure 2 shows the simulated frequency scaling of doubly clamped and circular drumhead graphene resonators using Equations (10) and (19) at zero gate voltage Vg = 0 V as a function of device characteristic dimension (length for the doubly clamped structures, and diameter for the circular drumhead devices) with different built-in strain levels of 0.002%, 0.01%, 0.05%, and 0.25%, respectively. Due to possible photoresist residue on the device during fabrication and surface adsorbates, we assume the mass is higher than the intrinsic device mass estimated by device dimensions and mass density, and we take an effective mass ratio of 2 (i.e., mdevice = 2 × mgraphene) as a typical value. With strain levels from 0.002% to 0.25%, the resonance frequency depends on the length for the doubly clamped case and the diameter for the circular drumhead case, with fres~L−1, and fres~D−1 power laws, at zero bias condition. It is clearly shown that in both doubly clamped and circular drumhead cases, resonance frequency increases as the characteristic device dimension decreases and as built-in strain increases. For the same characteristic length and built-in strain, the circular device gives higher resonance frequency than the doubly clamped device does, since its circumference-clamped structure provides higher spring constant and smaller effective mass. These correlations for both cases agree with the results from previous models [6,9,10,19,20,21,22,23,24] at gate voltage Vg = 0 V. To enable resonance frequency above gigahertz, fres > 1 GHz, built-in strain ≈0.25% with length smaller than 0.4 μm or built-in strain ≈0.05% with length smaller than 0.18 μm are required for doubly clamped graphene resonators. For circular drumhead resonators to achieve fres > 1 GHz, it needs a built-in strain ≈0.05% with a diameter smaller than 0.22 μm or a built-in strain ≈0.25% with a diameter smaller than 0.6 μm. Accordingly, fabricating smaller size membranes, larger built-in strain, and circular drumhead structure are preferred to attain higher resonance frequency.
Now we focus on electrostatic frequency tuning with respect to gate voltage Vg. At different built-in strain levels, we find three frequency tuning behaviors: resonance frequency increases monotonically with |Vg|, it decreases monotonically with |Vg|, and it first decreases, then increases with increasing |Vg|. All these cases are observed in existing experiments [5,6,7,8,9,10,11,12,13,14,15,16] for both of doubly clamped and circular drumhead graphene resonators. Figure 3 shows simulated frequency tuning of the doubly clamped and circular drumhead single-layer (1L) graphene resonators, with typical effective mass ratio = 2, air gap z0 = 300 nm, length L = 1 μm and diameter D = 1 μm, respectively. The results show that a big difference between our model and the previous modeling. For the device with relatively small built-in strain (ε0 = 0.01% in Figure 3a,d), applying |Vg| leading to elastic stiffening which is much bigger than that of capacitive softening, causing frequency elevation (“U” shape). Whereas for the device with large built-in strain (ε0 = 0.25% in Figure 3c,f), the capacitive softening dominates as |Vg| increases, which keeps reducing the resonance frequency. For intermediate built-in strain (ε0 = 0.05% in Figure 3b,e), initially the capacitive softening dominates, then the spring constant stiffening dominates as |Vg| increases, showing “W” shape frequency tuning. These tuning behaviors also affect frequency tuning range of the graphene resonators. For smaller built-in strain (ε0 = 0.01%), the resonance frequency of the doubly clamped device (L = 1 μm, z0 = 300 nm) shifts from 75 MHz to 98 MHz when |Vg|= 0 to 10 V (Figure 3a), showing the frequency tuning range of Δfres/fres ≈ 31%. For the device with the higher built-in strain (ε0 = 0.25%), although it leads to higher initial resonance frequency, the frequency shift is very small (378.9 MHz to 378.3 MHz), offering very limited frequency tunability, Δfres/fres ≈ 0.15%. These results suggest that there is an important trade-off between achieving high initial resonance frequency and wide frequency tuning range, which should be considered for designing the desired device performance for the specific applications. For example, devices with large initial strain are preferred for high frequency signal processing applications, and small initial strain facilitates realizing voltage controlled tunable devices that could be useful for making systems with tunable or programmable functions.
We also study frequency tuning with varying device parameters such as the depth of the air gap and the number of layers in the doubly clamped and circular drumhead resonators. Figure 4a,b shows that the device (L = 1 μm, ε0 = 0.01%) with the smaller air gap can provide much larger frequency tuning range for the doubly clamped and circular drumhead resonators (e.g., Δfres/fres ≈ 24.4% for z0 = 350 nm, and Δfres/fres ≈ 57.8% for z0 = 150 nm, for doubly clamped devices). Interestingly, with the same device parameters, such as the characteristic length, initial strain, depth of air gap, doubly clamped devices show larger tuning ranges when the gate voltage is changed (|Vg| = 0–10 V) (Δfres/fres ≈ 57.8% for L = 1 μm, z0 = 150 nm, ε0 = 0.01%), compared to those of the circular drumhead resonators (Δfres/fres ≈ 48.2% for D = 2R = 1 μm, z0 = 150 nm, ε0 = 0.01%). This is attributed to the fact that the fully clamped circular drumhead devices have larger stiffness at the beginning, making it less responsive to incremental external electrostatic force.
Further, we have investigated the dependence of resonance characteristics on a number of layers. Figure 4c,d shows frequency tuning in single-layer (1L) and few-layer graphene doubly clamped and circular drumhead resonators. In the same device parameters (L = 1 μm, ε0 = 0.01%, z0 = 300 nm) without an applied gate voltage, the resonance frequencies are independent of the number of layers. With the gate voltage, the resonance frequency for the 1L graphene resonators exhibits a much higher frequency tuning capability, compared to those in the few-layer graphene resonators (Δfres/fres ≈ 30.2% for 1L, Δfres/fres ≈ 17.5% for 2L, and Δfres/fres ≈ 12.0% for 3L doubly clamped graphene resonators), and we also find a similar trend in the circular drumhead resonators. According to Equations (9) and (19), the elastic stiffening is independent of the thickness t, while the capacitance softening increases with the increasing thickness of the membrane. For the built-in strain ≈0.01% in these cases, the softening effects of the few-layer are stronger compared to the softening in the single layer. Therefore, 1L graphene resonators give a much higher frequency at the same voltage and thus a wider tuning capability, compared to those in the few-layer graphene resonators.
We plot the frequency tuning range to clearly display tuning capability with respect to important parameters, including characteristic dimension, depth of air gap, and strain (Figure 5). Figure 5 shows the frequency tuning at gate voltage from Vg = 0 V to 10V as a function of the characteristic dimension and air gap with high built-in strain (0.2%) and low built-in strain (0.002%), for both doubly clamped and circular drumhead graphene resonators, respectively. With an air gap from 200 nm to 350 nm, length or diameter from 0.2 μm to 2 μm, and lower built-in strain (0.002%), this gives a larger frequency tuning range, while a higher built-in strain (0.2%) leads to a smaller frequency tuning range for both types of geometries for the graphene resonators. We find that a smaller depth of air gap, smaller initial strain, and longer characteristic dimension help achieve a wide frequency tuning range.
Figure 6 shows the experimental frequency tuning data from [6] and fitted results using our model (blue solid lines) and previous modeling (red dashed lines) [24]. Due to the trapped charges in the membrane, the symmetry axis of the frequency tuning curve is shifted from zero gate voltage. The results show that our model agrees well with experimental data and it is much more accurate compared to that of the previous modeling, particularly at a high voltage regime. In a previous model [24], the equilibrium displacement ze, which is a positive correlate to the strain, depends on the DC gate voltage Vg, with ze~ V g 2 power laws [6,10]. This is the only correct form for sufficiently small Vg; while we expect ze~ V g 2 / 3 for large Vg. This is the reason that the estimated resonance frequency from the previous modeling is much larger than the experimental results. In our model, we obtain ze by performing several iterative calculations, to make the model accurate at a wide voltage range. Based on the iterative calculations, we solve the effective spring constant by considering both static deflection and the fundamental mode shape, which are largely unexplored in the previous model.

4. Conclusions

Using iterative computational modeling of device capacitance, we have developed comprehensive models for frequency tuning behaviors in doubly clamped and circular drumhead graphene resonators. We have also examined the effect of various parameters such as built-in tension, characteristic dimension, depth of air gap, number of layers, and device structure on frequency tuning, which provides useful guidelines for future design and fabrication to achieve new tunable graphene devices with the desired resonance frequencies and tuning ranges.

Author Contributions

T.M. performed the modeling and computation and analyzed the data with J.L. and P.X.-L.F. T.M., J.L., Y.X. and P.X.-L.F. wrote the manuscript. P.X.-L.F. and Y.X. conceived and supervised the project.

Acknowledgments

We thank the support from National Science Foundation CCSS Award (Grant ECCS-1509721), Fundamental Research Funds for the Central Universities (Grant ZYGX2016J036), the National Natural Science Foundation of China (Grant Number 61474020), Special Innovation Zone of NDST (Grant Number: 1716313-ZD00100103), and the China Scholarship Council (File Number 201706070021).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mak, K.F.; Sfeir, M.Y.; Wu, Y.; Lui, C.H.; Misewich, J.A.; Heinz, T.F. Measurement of the optical conductivity of graphene. Phys. Rev. Lett. 2008, 101, 196405. [Google Scholar] [CrossRef] [PubMed]
  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, 614–617. [Google Scholar] [CrossRef] [PubMed]
  3. 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. [Google Scholar] [CrossRef] [PubMed]
  4. Lee, G.-H.; Cooper, R.C.; An, S.J.; Lee, S.; van der Zande, A.; Petrone, N.; Hammerberg, A.G.; Lee, C.; Crawford, B.; Oliver, W.; et al. High-strength chemical-vapor–deposited graphene and grain boundaries. Science 2013, 340, 1073–1076. [Google Scholar] [CrossRef] [PubMed]
  5. Bunch, J.S.; van der Zande, A.M.; Verbridge, S.S.; Frank, I.W.; Tanenbaum, D.M.; Parpia, J.M.; Craighead, H.G.; McEuen, P.L. Electromechanical resonators from graphene sheets. Science 2007, 315, 490–493. [Google Scholar] [CrossRef] [PubMed]
  6. Chen, C.Y.; Rosenblatt, S.; Bolotin, K.I.; Kalb, W.; Kim, P.; Kymissis, I.; Stormer, H.L.; Heinz, T.F.; Hone, J. Performance of monolayer graphene nanomechanical resonators with electrical readout. Nat. Nanotechnol. 2009, 4, 861–867. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Xu, Y.H.; Chen, C.Y.; Deshpande, V.V.; DiRenno, F.A.; Gondarenko, A.; Heinz, D.B.; Liu, S.; Kim, P.; Hone, J. Radio frequency electrical transduction of graphene mechanical resonators. Appl. Phys. Lett. 2010, 97, 243111. [Google Scholar] [CrossRef] [Green Version]
  8. Singh, V.; Sengupta, S.; Solanki, H.S.; Dhall, R.; Allain, A.; Dhara, S.; Pant, P.; Deshmukh, M.M. Probing thermal expansion of graphene and modal dispersion at low-temperature using graphene nanoelectromechanical systems resonators. Nanotechnology 2010, 21, 165204. [Google Scholar] [CrossRef] [PubMed]
  9. Weber, P.; Güttinger, J.; Tsioutsios, I.; Chang, D.E.; Bachtold, A. Coupling graphene mechanical resonators to superconducting microwave cavities. Nano Lett. 2014, 14, 2854–2860. [Google Scholar] [CrossRef] [PubMed]
  10. Miao, T.; Yeom, S.; Wang, P.; Standley, B.; Bockrath, M. Graphene nanoelectromechanical systems as stochastic-frequency oscillators. Nano Lett. 2014, 14, 2982–2987. [Google Scholar] [CrossRef] [PubMed]
  11. Chen, C.; Lee, S.; Deshpande, V.V.; Lee, G.-H.; Lekas, M.; Shepard, K.; Hone, J. Graphene mechanical oscillators with tunable frequency. Nat. Nanotechnol. 2013, 8, 923–927. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Ye, F.; Lee, J.; Feng, P.X.-L. Electrothermally tunable graphene resonators operating at very high temperature up to 1200 K. Nano Lett. 2018, 18, 1678–1685. [Google Scholar] [CrossRef] [PubMed]
  13. Mathew, J.P.; Patel, R.N.; Borah, A.; Vijay, R.; Deshmukh, M.M. Dynamical strong coupling and parametric amplification of mechanical modes of graphene drums. Nat. Nanotechnol. 2016, 11, 747–751. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Alba, R.D.; Massel, F.; Storch, I.R.; Abhilash, T.S.; Hui, A.; McEuen, P.L.; Craighead, H.G.; Parpia, J.M. Tunable phonon-cavity coupling in graphene membranes. Nat. Nanotechnol. 2016, 11, 741–746. [Google Scholar] [CrossRef] [Green Version]
  15. Davidovikj, D.; Alijani, F.; Cartamil-Bueno, S.J.; van der Zant, H.S.J.; Amabili, M.; Steeneken, P.G. Nonlinear dynamic characterization of two-dimensional materials. Nat. Commun. 2017, 8, 1253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Güttinger, J.; Noury, A.; Weber, P.; Eriksson, A.M.; Lagoin, C.; Moser, J.; Eichler, C.; Wallraff, A.; Isacsson, A.; Bachtold, A. Energy-dependent path of dissipation in nanomechanical resonators. Nat. Nanotechnol. 2017, 12, 631–636. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Lee, J.; Krupcale, M.J.; Feng, P.X.-L. Effects of γ-ray radiation on two-dimensional molybdenum disulfide (MoS2) nanomechanical resonators. Appl. Phys. Lett. 2016, 108, 023106. [Google Scholar] [CrossRef]
  18. He, R.; Feng, X.L.; Roukes, M.L.; Yang, P. Self-transducing silicon nanowire electromechanical systems at room temperature. Nano Lett. 2008, 8, 1756–1761. [Google Scholar] [CrossRef] [PubMed]
  19. Lee, J.; Wang, Z.; He, K.; Yang, R.; Shan, J.; Feng, P.X.-L. Electrically tunable single- and few-layer MoS2 nanoelectromechanical systems with broad dynamic range. Sci. Adv. 2018, 4, eaao6653. [Google Scholar] [CrossRef] [PubMed]
  20. Xu, Y.H.; Li, O.P.; Xu, R.M. Graphene resonant channel transistor. In Proceedings of the IEEE International Wireless Symposium (IWS), Beijing, China, 14–18 April 2013; pp. 1–6. [Google Scholar]
  21. Lekas, M.; Lee, S.; Cha, W.; Hone, J.; Shepard, K. Noise modeling of graphene resonant channel transistors. IEEE Trans. Electron Devices 2015, 62, 1276–1283. [Google Scholar] [CrossRef]
  22. Mei, T.D.; Xu, Y.H.; Li, O.P.; Lan, Y.; Wu, Y.Q.; Xu, R.M.; Chen, Y.F.; Li, Y.R. Accurate multi-bias equivalent circuit model for graphene resonant channel transistors. In Proceedings of the IEEE International Microwave Symposium (IMS), San Francisco, CA, USA, 22–27 May 2016; pp. 1–4. [Google Scholar]
  23. Mei, T.D.; Xu, Y.H.; Lan, Y.; Li, O.P.; Sander, M.R.; Xu, R.M.; Li, Y.R. A high-frequency compact model for graphene resonant channel transistors including mechanical nonlinear effects. IEEE Trans. Microw. Theory Tech. 2017, 65, 4063–4072. [Google Scholar] [CrossRef]
  24. Chen, C.Y. Graphene nanoElectroMechanical Resonators and Oscillators. Ph.D. Thesis, Columbia University, New York, NY, USA, 2013. [Google Scholar]
  25. Schomburg, W.K. Introduction to Microsystem Design; Springer: New York, NY, USA, 2011. [Google Scholar]
Figure 1. Device structures and modeling procedure of electrostatic tuning of resonance frequency in graphene nanoelectromechanical systems (NEMS). (a) Deflection of graphene resonator under electrostatic force. Without gate voltage, the membrane is flat (thin magenta dashed line) and suspended over the trench at equilibrium. With applied DC gate voltage, the membrane deflects to new equilibrium displacement ze (thick black solid line). With AC actuation, the membrane vibrates with amplitude δz (thick black dashed lines). (b,c) Schematics of doubly clamped and circular drumhead graphene resonators, respectively. (d) Flowchart for computational implementation of the frequency tuning model.
Figure 1. Device structures and modeling procedure of electrostatic tuning of resonance frequency in graphene nanoelectromechanical systems (NEMS). (a) Deflection of graphene resonator under electrostatic force. Without gate voltage, the membrane is flat (thin magenta dashed line) and suspended over the trench at equilibrium. With applied DC gate voltage, the membrane deflects to new equilibrium displacement ze (thick black solid line). With AC actuation, the membrane vibrates with amplitude δz (thick black dashed lines). (b,c) Schematics of doubly clamped and circular drumhead graphene resonators, respectively. (d) Flowchart for computational implementation of the frequency tuning model.
Micromachines 09 00312 g001
Figure 2. Calculated resonance frequencies of single-layer (1L) graphene resonators without applying gate voltage. Due to photoresist residue in fabrication and surface absorption, we assume the mass is higher than the intrinsic device mass (of carbon atoms only), and we take an effective mass ratio of 2 as the typical value. Frequency scaling of (a) doubly clamped membranes and (b) circular drumhead graphene resonators. The labels represent built-in strain levels in the devices.
Figure 2. Calculated resonance frequencies of single-layer (1L) graphene resonators without applying gate voltage. Due to photoresist residue in fabrication and surface absorption, we assume the mass is higher than the intrinsic device mass (of carbon atoms only), and we take an effective mass ratio of 2 as the typical value. Frequency scaling of (a) doubly clamped membranes and (b) circular drumhead graphene resonators. The labels represent built-in strain levels in the devices.
Micromachines 09 00312 g002
Figure 3. Three typical frequency tuning behaviors for both a doubly clamped graphene resonator (upper row) and a circular drumhead graphene resonator (lower row) with varying built-in strain of (a,d) 0.01%, (b,e) 0.05%, and (c,f) 0.25%. The red dashed lines show previous modeling and blue solid lines show this model by assuming that effective mass ratio = 2, air gap z0 = 300 nm, length L = 1 μm and diameter D = 1 μm, respectively.
Figure 3. Three typical frequency tuning behaviors for both a doubly clamped graphene resonator (upper row) and a circular drumhead graphene resonator (lower row) with varying built-in strain of (a,d) 0.01%, (b,e) 0.05%, and (c,f) 0.25%. The red dashed lines show previous modeling and blue solid lines show this model by assuming that effective mass ratio = 2, air gap z0 = 300 nm, length L = 1 μm and diameter D = 1 μm, respectively.
Micromachines 09 00312 g003
Figure 4. Dependence of resonance characteristics on depth of air gap and number of layers. Frequency tuning of doubly clamped graphene resonators with varying (a) air gap and (c) number of layers. Simulated resonance frequencies and tuning characteristics for circular drumhead resonators with (b) different depth of air gap and (d) number of layers.
Figure 4. Dependence of resonance characteristics on depth of air gap and number of layers. Frequency tuning of doubly clamped graphene resonators with varying (a) air gap and (c) number of layers. Simulated resonance frequencies and tuning characteristics for circular drumhead resonators with (b) different depth of air gap and (d) number of layers.
Micromachines 09 00312 g004
Figure 5. 3D plots of computed frequency tuning with varying characteristic dimension and air gap with high built-in strain (0.2%) and low built-in strain (0.002%). (a) Frequency tunability of doubly clamped single-layer (1L) graphene resonators, and (b) circular drumhead 1L graphene resonators.
Figure 5. 3D plots of computed frequency tuning with varying characteristic dimension and air gap with high built-in strain (0.2%) and low built-in strain (0.002%). (a) Frequency tunability of doubly clamped single-layer (1L) graphene resonators, and (b) circular drumhead 1L graphene resonators.
Micromachines 09 00312 g005
Figure 6. Fitting frequency tuning model to experimental data. The measured data from [6] and fitted data with our modeling (blue solid lines) and previous modeling (red dashed lines) [24] by assuming that built-in strain ε0, effective mass ratio, length L, air gap z0 are (a) 0.003%, 3.5, 1.1 μm 250 nm and (b) 0.23%, 5.7, 1.8 μm 250 nm, respectively.
Figure 6. Fitting frequency tuning model to experimental data. The measured data from [6] and fitted data with our modeling (blue solid lines) and previous modeling (red dashed lines) [24] by assuming that built-in strain ε0, effective mass ratio, length L, air gap z0 are (a) 0.003%, 3.5, 1.1 μm 250 nm and (b) 0.23%, 5.7, 1.8 μm 250 nm, respectively.
Micromachines 09 00312 g006

Share and Cite

MDPI and ACS Style

Mei, T.; Lee, J.; Xu, Y.; Feng, P.X.-L. Frequency Tuning of Graphene Nanoelectromechanical Resonators via Electrostatic Gating. Micromachines 2018, 9, 312. https://doi.org/10.3390/mi9060312

AMA Style

Mei T, Lee J, Xu Y, Feng PX-L. Frequency Tuning of Graphene Nanoelectromechanical Resonators via Electrostatic Gating. Micromachines. 2018; 9(6):312. https://doi.org/10.3390/mi9060312

Chicago/Turabian Style

Mei, Tengda, Jaesung Lee, Yuehang Xu, and Philip X.-L. Feng. 2018. "Frequency Tuning of Graphene Nanoelectromechanical Resonators via Electrostatic Gating" Micromachines 9, no. 6: 312. https://doi.org/10.3390/mi9060312

APA Style

Mei, T., Lee, J., Xu, Y., & Feng, P. X. -L. (2018). Frequency Tuning of Graphene Nanoelectromechanical Resonators via Electrostatic Gating. Micromachines, 9(6), 312. https://doi.org/10.3390/mi9060312

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop