# NEMS Resonators for Detection of Chemical Warfare Agents Based on Graphene Sheet

Authors: Nikola Anđelić; Zlatan Car; Marko Čanađija
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

Graphene sheets are the basis of nanoelectromechanical (NEMS) mass resonators that are in recent years experimentally used for detection of specific gas molecules in air. The aim of this paper is to theoretically investigate the possibility of application of graphene sheets in detection of Chemical Warfare Agents (CWA’s). By application of the nonlocal theory of elasticity and classical plate theory, equations that describe natural vibrations of the simply supported single-layer graphene sheet (SLGS) and double-layer graphene sheet (DLGS) were derived. In order to detect attached CWA molecule, the frequency shift method was applied. The results indicate that the proposed methodology could be successfully implemented in order to detect an attached CWA molecule. However, only specific mode shapes could be employed for such detection.


---

## Body

## 1. Introduction

Due to its extraordinary mechanical, electrical, and thermal properties graphene sheets have many potential applications such as reinforced materials, molecule sensors, and NEMS mass resonators [1]. For structural applications of graphene sheets, it is important to know their macroscopic properties. Therefore, experimental and theoretical analysis of graphene sheets are necessary in order to design optimal materials and devices. Nowadays a lot of resources and effort is focused on extensive research for possible application of graphene sheets in high-performance hybrid supercapacitors, optoelectronic devices, and various types of high-performance sensors. Although the graphene is material with extraordinary characteristics, studies [2, 3] are conducted to examine their potential impact on health and environment. Nanoelectromechanical systems or shortly NEMS are strong candidates for variety of applications in semiconductor-based technology and fundamental science [4, 5]. NEMS resonators are used as the precision mass sensors which are used to weigh cells, biomolecules, and gas molecules [6–8]. A NEMS mass sensor relies on monitoring how the resonance frequency of the NEMS resonator changes when an additional mass is absorbed onto its surface. Today there is intense interest in implementation of the graphene sheet as the building block of the NEMS resonator based on their material properties such as low specific mass and high stiffness. The research [9] introduced a method for the large-scale production of the high-quality free-standing graphene NEMS in which the graphene top surface is continuously covered by the supporting polymer throughout the fabrication process. Although graphene sheet is highly sensitive material and used is for gas sensing applications in research [10] the effective, single-step method is shown to enhance its sensitivity.Graphene sheets based on their structural features can be divided as Pristine Graphene (PG), Graphene Oxide (GO), and Reduced Graphene Oxide (RGO) [11]. Each of these groups of graphene sheets has distinct gas sensing capabilities. According to the research [12] GO proved to be promising candidate of gas detection under standard humidity conditions. Transducing mechanisms that can be used for detection of nanoparticles are optical [13, 14] (luminescence, absorption, polarization, and fluorescence), mass [15] (acoustic waves, absorption, and polarization), electrochemical [16] (conductometric, amperometric, potentiometric, and ion-sensitive), piezoelectric [17], and thermal. Today most commonly used transducing mechanisms are optical, electrochemical, and mass. The research [18] showed that two-dimensional structure of graphene makes electron transport through graphene highly sensitive to the adsorption of the gas molecules. When certain gas molecule attached itself on the graphene surface, this leads to change in the electrical conductivity. The conductivity change can be attributed to the change in the local carrier concentration induced by the surface adsorbates which act as electron donors or electron acceptors. NEMS with mass transducing mechanisms operate on the frequency shift principle. This means that resonant frequencies are determined for the system without attached nanoparticles and for the system with attached nanoparticles. Based on the difference between these two frequencies the type of chemical compound can be approximately determined. So far, the NEMS resonators based on graphene sheets were being used for investigating detection from theoretical and experimental point of view with simple chemical compounds such as H2, H20, O2, CO2, CO, NO2, Fe, Co, Ni, Ru, Rh, Pd, Os, Ir, and Pt [11]. Although the graphene as a building block of gas sensor has high sensitivity to gas molecules the proof-of-principle comparison studies [19] showed that graphene sensing capabilities could be enhanced using a nanoporous substrate.The tendency among researchers [20–22] is to apply the NEMS resonators in detection of more complex chemical compounds. In the focus of the research at hand is detection of Chemical Warfare Agents (CWA). According to Organization for the Prohibition of Chemical Weapons (OPCW) document [23] the CWAs are classified according to mechanism of toxicity in humans into blister agents, nerve agents, asphyxiants, choking agents, and incapacitating/behavior altering agents. Blister agents are sulfur mustard, nitrogen mustard and lewisite. According to NATO codes nerve agents can be categorized as G series agents, GB (sarin), GD (Soman), GA (Tabun), GF (cyclosarin), and V series agents, VE, VG, VM, and VX. The letter G represents the country of origin which is Germany and letter V represents Venomous. Asphyxiants agents are cyanogen chloride, hydrogen cyanide, and arsine. Chocking agents are chlorine, chloropicrin, and phosgene [24].The goal of this research is to develop a theoretical framework that can be implemented into NEMS in order to detect presence of CWA. To the best of authors’ knowledge, such a framework was not proposed in the literature so far. The concept relies on the small size effects observed in both one-dimensional [25] and two-dimensional structures [26] operating at nano- and microscale. In short, forces that are irrelevant at the macroscale become important at smaller scales. In that way, the stress state at some point does not depend only on the strain state at the same location, but on strains in whole neighborhood of the point in question. This approach in mechanics is known as the nonlocal theory. The theory is based on two main methodologies. Currently predominating one is based on gradients developed for either isothermal [26–28] or nonisothermal environments [29–31]. The newer approach is based on the integral approach [32, 33]. Consequently, in this paper, the vibrations of simply supported graphene sheet (GS) with and without attached CWA molecule are analyzed using the gradient based nonlocal plate theory. In a nutshell, the nonlocal plate theory is an extension of the Kirchhoff-Love plate theory and is used to model mechanical response of single-layer graphene sheets (SLGS) and double-layer graphene sheets (DLGS). Today, the nonlocal theory is extensively applied in order to investigate mechanical and vibrational properties of SLGS and DLGS. In [4, 34] the authors used a continuum model to develop the equation which describes vibrations of multilayered graphene sheets (MLGSs) and in their investigation the formula is derived to predict the van der Waals (vdW) interaction. Large amplitude vibrations of multilayered GS were analyzed in [35] using Hamilton’s principle in order to obtain the coupled nonlinear partial differential equations of motion which are based on the von Karman strain theory and Eringen theory of nonlocal continuum. The circular plate theory with vdW force was applied in [4] to investigate vibrations of the circular DLGS system. In [36] the higher order shear deformation theory is used to investigate buckling properties of single graphene sheet. Using this theory the authors developed closed-form solutions for critical buckling forces of the graphene sheets. In [37] shear deformation theory for buckling analysis of SLGS is developed using nonlocal differential constitutive relations of Eringen. Nonlocal elasticity theory was employed to investigate effects of small scale on buckling of the rectangular nanoplate. The equations of motion of the nonlocal theories are derived and solved via Navier's procedure for all edges simply supported boundary conditions. Obtained results are verified with the known results in the literature. The influences played by effects of nonlocal parameter, length, and thickness of the graphene sheets and shear deformation effect on the critical buckling load are studied. Verification studies show that the proposed theory is not only accurate and simple in solving the buckling nanoplates, but also comparable with the other higher order shear deformation theories which contain more number of unknowns.To summarize, the main novelty of this paper is to investigate the possibility of utilizing SLGS and DLGS as sensing elements of NEMS resonators in detection of CWA molecules based on absolute and relative frequency shift method. In [1] the authors applied the nonlocal theory in order to obtain frequency equations of simply supported SLGS and DLGS systems used as the nanomechanical mass sensor. However, during the derivation of the frequency equation for simply supported DLGS system, the nonlocal parameter was neglected so the final equation for this system is reduced to the classical (local) form. The research in this paper starts from these existing results. After the derivation of the nonlocal equation which describes the vibrations of simply supported SLGS and DLGS, influence of the nonlocal parameter, geometry of SLGS and DLGS, and nanoparticle mass (CWA molecule) attached to the GS on vibrations of SLGS is analyzed and results are graphically represented. In this investigation the center of SLGS and DLGS upper plate is considered as a point of the attached CWA molecules. The influence of attached CWA molecule position on resonant frequencies of SLGS is also investigated. The vibrations of SLGS and DLGS with and without attached nanoparticle are determined and frequency shift as a function of the nonlocal parameter and nanoparticle mass is also analyzed.


## 2. Nonlocal Theory or Elasticity

In the general elasticity model, the stress components at any point depend only on the strain component at the same point. In the nonlocal elasticity theory, the stress field at a certain point is a function of the strain distribution over a certain representative volume of the material centered at that point. This theory is used to study several phenomena related to nanoscale materials. The equation of motion for nonlocal linear elastic solids can be written in the following form:(1)tij+fi=ρui¨,where ρ,fi,ui,andtij are the mass density, applied body forces vector, displacement vector, and the nonlocal stress tensor, respectively. The nonlocal stress tensor can be written in the following form:(2)tij=∫Vαx-x′,e0aσijx′dVx′,where x is the reference point in the body and e0a is the nonlocal parameter or scale coefficient. According to [29] the material parameter e0 is a constant appropriate to each material and it is estimated in a way that relations of the nonlocal elasticity model could provide approximation of atomic dispersion curves of plane waves with those of atomic lattice dynamics. The value of material parameter e0 can be found by comparing results obtained using nonlocal theory with results obtained with molecular dynamics simulation and atomistic based techniques. The parametera can be the lattice parameter, granular size, distance between atoms or similar quantitiy. x-x′ is the Euclidian distance andV is the region occupied by the body. In previous equation the term αx-x′,e0a represents the two-dimensional nonlocal modulus and is given as(3)αx,e0a=12πe0a2K0x∙xe0a,where K0 is modified Bessel function. Integration of (3) over the entire domain will result in unity. If the value e0a in the kernel function α is equal to 0 then the nonlocal elasticity is reduced to classical elasticity. In analysis of graphene sheets and carbon nanotubes according to [38] the most commonly used nonlocal kernel function is the one represented in (3) and will be also used in this paper.A parameterσij in (2) represents the local stress tensor of the classical elasticity theory at any point x′ in the body and follows from the usual constitutive relation: (4)σijx′=λεkkx′δij+2μεijx′,where λ and μ are Lame constants and εij is the strain tensor.


### 2.1. Single-Layer Graphene Sheets

As already mentioned, the dynamic behavior of a single-layer graphene sheet (SLGS) with the attached concentrated massmc located at an arbitrary position x0,y0 is considered using nonlocal continuum mechanics. The SLGS is assumed to be simply supported at the boundary. The schematic illustration of the problem at hand is shown in Figure 1.Figure 1
Schematic illustration of SLGS with attached mass at arbitrary position.The origin of coordinate system is positioned in the lower left corner of the mid-plane of SLGS whilex andy axes are oriented along the length La and Lb, respectively. The out of plane axisz is taken along the thicknessh of SLGS.By applying the nonlocal elasticity theory [38] on SLGS the two-dimensional nonlocal constitutive equations are obtained and can be written in the following form:(5a)σxx-e0a2∂2σxx∂x2+∂2σxx∂y2=E1-v2εxx+vεyy,(5b)σyy-e0a2∂2σyy∂x2+∂2σyy∂y2=E1-v2εyy+vεxx,(5c)τxy-e0a2∂2σxy∂x2+∂2σxy∂y2=Gγxy,whereE,G, andv are the elastic modulus, the shear modulus, and Poisson’s ratio of the SLGS, respectively. The internal characteristic lengtha corresponds to the distance between two atoms in a C-C bond and is equal to 0.142 nm. Flexural moments of SLGS are obtained using the following expressions:(6)Mxx=∫-h/2h/2zσxxdz,Myy=∫-h/2h/2zσyydz,Mxy=∫-h/2h/2zτxydz.If in-plane displacements of the middle surface are neglected in thex andy direction then the relationship between strain and displacement fields is expressed as (7)εxx=-z∂2w∂x2,εyy=-z∂2w∂y2,γxy=-2z∂2w∂x∂y,where w is displacement in direction ofz axis of GS. With substitution of (5a), (5b), and (5c) into (6) and using (7) the following expressions are obtained:(8a)Mxx-e0a2∂2Mxx∂x2+∂2Mxx∂y2=-D∂2w∂x2+v∂2w∂y2(8b)Myy-e0a2∂2Myy∂x2+∂2Myy∂y2=-D∂2w∂y2+v∂2w∂x2,(8c)Mxy-e0a2∂2Mxy∂x2+∂2Mxy∂y2=-D1-vG∂2w∂x∂y,whereD is the flexural rigidity of SLGS, expressed as (9)D=Eh3121-v2.The governing equation that describes the flexural vibration of SLGS carrying a nanoparticle at arbitrary position(x0,y0) can be written in the following form: (10)∂2Mxx∂x2+2∂2Mxy∂x∂y+∂2Myy∂y2=ρh+mcδx-x0δy-y0∂2w∂t2,where ρ,mc,t, and δ are the mass density of SLGS, mass of nanoparticle, time, and the Dirac delta function, respectively.With substitution of (8a), (8b), and (8c) into (10), the following governing equation is obtained: (11)D∂4w∂x4+2∂4w∂x2∂y2+∂4w∂y4+1-e0a2∂2w∂x2+∂2w∂y2ρh+mcδx-x0δy-y0∂2w∂t2=0.The harmonic solution of previous equation can be expressed as (12)wx,y,t=Yx,yeiωt,where Yx,y represents a shape function of deflection and ω is the resonant frequency of the SLGS.Substituting (12) into (11), the governing equation can be expressed as (13)∂4Y∂x4+2∂4Y∂x2∂y2+∂4Y∂y4-ω2D1-e0a2∂2Y∂x2+∂2Y∂y2ρh+mcδx-x0δx-y0∂2Y∂t2=0.Since the GS is simply supported, the following boundary conditions are applied:(14)w=0,∂2w∂x2=0,∂2w∂y2=0onx=0,La,andy=0,Lb.The shape function from (12) can be expressed as (15)Yx,y=Amnsin⁡mπxLasin⁡nπyLb,where parameter Amn represents the vibration amplitude and indexesm andn represent mode numbers in the periodic directions, respectively.In order to obtain the equation for determining frequencies of SLGS, (15) is substituted in (13). After the substitution, the equation is multiplied by x=sin⁡(mπx/La)sin⁡(nπy/Lb) and integrated over the whole region with respect tox andy in limitsx = 0 toLa andy = 0 toLb. With some simplifications the following frequency equation is obtained:(16)∫0Lb∫0LaAmnDπ4m2La2+n2Lb22sin2⁡mπxLasin2⁡nπyLbdxdy-ω2∫0Lb∫0LaAmn1+e0a2π2m2La2+n2Lb2ρh+mcδx-x0δy-y0sin2⁡mπxLasin2⁡nπyLbdxdy=0.The roots of the previous equation are the desired resonant frequencies and each resonant frequency corresponds to a given shape function. After further manipulation of the previous equation, the expression for determining resonant frequency is obtained in the following form:(17)ωmn2=Dπ4m2/La2+n2/Lb221+e0a2π2m2/La2+n2/Lb2ρh+4mc/LaLbsin2⁡mπξsin2⁡nπη,where ξ=x0/La and η=y0/Lb are the nondimensional parameters of nanoparticle position. In (17) in denominator sin2⁡mπξsin2⁡nπη is multiplied by the nonlocal parameter (e0a) so the nonlocal parameter is influencing the mass localization on SLGS. The frequency is obtained from (18)fmn=ωmn2π.It is worth mentioning that when the nonlocal parametere0a is equal to zero, the resonant frequency of SLGS with attached nanoparticles collapses to the classical plate theory.


### 2.2. Double-Layer Graphene Sheets

DLGS are composed of two SLGS and they interact with each other by vdW forces. In these simulations, the nanoparticle is attached on the upper layer while no nanoparticles are attached to the lower level of GS. The DLGS is assumed to be simply supported at the boundary. The DLGS system with attached nanoparticle at arbitrary position is shown in Figure2.Figure 2
Schematic illustration of DLGS with attached mass at arbitrary position.The governing equation that describes vibration of this system can be written in the following form:(19a)D∇4w1+1-e0a2∇2ρh+mcδx-x0δy-y0∂2w1∂t2=p12,(19b)D∇4w2+ρh1-e0a2∇2∂2w2∂t2=p21,where wjx,y,t,j=1,2 are the flexural deflections of upper j=1 and lower j=2 sheet and p12p21=-p12 is the transverse pressure between layers of DLGS caused by vdW forces. As in the case of SLGS the variableD is flexural rigidity given by (9) while ∇4 is the biharmonic operator and in the Cartesian coordinate system it can be written in the following form: (20)∇4.=∂4∂x4+∂4∂x2∂y2+∂4∂y4.The distributed transverse pressure acting on the upper and lower layer of DLGS system can be written as(21)p12=cw2-w1,where c is the vdW interaction coefficient between upper and lower layers [4, 34, 35, 40–42]. It can be obtained from the Lennard-Jones potential as (22)c=439a224εσ2σa83003π256∑k=05-1k2k+15kσa61z¯1-z¯212-35π8∑k=02-1k2k+12k1z¯1-z¯26,where ε=2.968meV and σ=0.34nm are parameters that are chosen to fit the physical properties of graphene sheets. The parameter z-j can be determined from the following equation: (23)z-j=zja,j=1,2,where zj is the coordinate ofj-th layer in the direction of the thickness with the origin in the mid-plane of the graphene sheets, and parameter a=0.142nm is the C-C bond length.The second term in (19a) is given by (24)μx,y=ρh+mcδx-x0δy-y0,where μx,y is the variable mass distribution function of the upper sheet.As in the case of SLGS, to obtain the resonant frequencies of DLGS system described by governing equations (19a) and (19b) the following function is introduced:(25)wjx,y,t=Yjx,yeiωt,j=1,2,where Yjx,y,j=1,2 represents a shape function of deflection in the upper and the lower sheet of DLGS system and ω is the resonant frequency.(26)wjx,y,t=Yjx,yeiωt,j=1,2,Substituting (26) into differential equations (19a) and (19b) of DLGS system the following equations are obtained in the matrix form: (27)D∇4-μω2e0a2∇2+c-μω2-c-cD∇4-ρhω2e0a2∇2+c-ρhω2Y1Y2=00.With further algebraic manipulation of the previous equation it is obtained that(28)∇8Yj+μ+ρhω2e0a2D2∇6Yj+2c-ρhω2-μω2D+μρhω4e0a4D2∇4Yj+cρh+μ-2ρhω2ω2e0a2D2∇2Yj+ρhμω4-cρh+μω2D2Yj=0,j=1,2.Again, the DLGS system is simply supported which means that the following boundary conditions are applied: (29)w=0,∂2w∂x2=0,∂2w∂y2=0onx=0,La,andy=0,Lb.To fulfill previously defined boundary conditions, the following shape function of deflection is proposed: (30)Yjx,y=Amnsin⁡mπxLasin⁡nπyLb,j=1,2,where Amn is the vibrational amplitude of oscillations, while indexesm andn indicate the mode numbers. With substitution of (30) into (28), the following equation is obtained: (31)π8Lb2m2+La2n24La8Lb8Yj-π6Lb2m2+La2n23La6Lb6μ+hρω2e0a2DYj+π4Lb2m2+La2n22La4Lb4D2c-ω2μ+hρ+ρhμω4e0a2D2Yj+π2Lb2m2+La2n2La2Lb2e0a2ω22hμρω2-cμ+ρhD2Yj+ω2hμρω2-cμ+ρhD2Yj=0.Next step is to substitute (24) into the previous equation:(32)
π
8
L
b
2
m
2
+
L
a
2
n
2
4
L
a
8
L
b
8
Y
j
-
π
6
L
b
2
m
2
+
L
a
2
n
2
3
L
a
6
L
b
6
h
ρ
ω
2
e
0
a
2
D
Y
j
+
π
4
L
b
2
m
2
+
L
a
2
n
2
2
L
a
4
L
b
4
2
c
-
ω
2
ρ
h
D
Y
j
+
π
2
L
b
2
m
2
+
L
a
2
n
2
L
a
2
L
b
2
e
0
a
2
ω
2
ρ
h
D
2
Y
j
+
ω
2
c
ρ
h
D
2
Y
j
=
π
6
L
b
2
m
2
+
L
a
2
n
2
3
L
a
6
L
b
6
ω
2
e
0
a
2
D
+
π
4
L
b
2
m
2
+
L
a
2
n
2
2
L
a
4
L
b
4
ω
2
-
ρ
h
ω
4
e
0
a
4
D
2
-
π
2
L
b
2
m
2
+
L
a
2
n
2
L
a
2
L
b
2
e
0
a
2
ω
2
2
h
ρ
ω
2
-
c
D
2
Y
j
-
ω
2
h
ρ
ω
2
-
c
D
2
m
c
δ
x
-
x
0
δ
y
-
y
0
Y
j
.After multiplication of the previous equation with the term sin⁡mπx/Lasin⁡nπy/Lb from both sides and with integration over the whole region with respect tox andy the following solution is obtained: (33)r0ω4+r1ω2+r2=0,where the coefficients r0,r1, and r2 are (34a)r0=18D2La4Lb4hLa2Lb2+n2La2+m2Lb2e0a2π22ρ2mc+hLaLbρ+mc-2cos⁡2nπη-2cos⁡2mπξ+cos⁡2πnη-mξ+cos⁡2πnη+mξ,(34b)r1=-12D2La6Lb6cijLa4Lb4+Dn2La2+m2Lb22π4m2Lb2e0a2π2+La2Lb2+n2e0a2π2hLaLbρ+2mcsin⁡nπη2sin⁡mπξ2,(34c)r2=n2La2+m2Lb22π42cijLa4Lb4+Dn2La2+m2Lb22π44DLa7Lb7.ξ and η, as in the SLGS case, represent the nondimensional parameter of the attached nanoparticle position. r0 and r1 depend on the nonlocal parameter so the nonlocal parameter influences natural frequencies as well. The solution of (33) gives a term for determining natural frequencies of DLGS system with attached nanoparticle and can be expressed as (35)ωmn2=-r1+r12-4r0r22r0,And, finally, the frequency is obtained from(36)fmn=ωmn2π.


## 2.1. Single-Layer Graphene Sheets

As already mentioned, the dynamic behavior of a single-layer graphene sheet (SLGS) with the attached concentrated massmc located at an arbitrary position x0,y0 is considered using nonlocal continuum mechanics. The SLGS is assumed to be simply supported at the boundary. The schematic illustration of the problem at hand is shown in Figure 1.Figure 1
Schematic illustration of SLGS with attached mass at arbitrary position.The origin of coordinate system is positioned in the lower left corner of the mid-plane of SLGS whilex andy axes are oriented along the length La and Lb, respectively. The out of plane axisz is taken along the thicknessh of SLGS.By applying the nonlocal elasticity theory [38] on SLGS the two-dimensional nonlocal constitutive equations are obtained and can be written in the following form:(5a)σxx-e0a2∂2σxx∂x2+∂2σxx∂y2=E1-v2εxx+vεyy,(5b)σyy-e0a2∂2σyy∂x2+∂2σyy∂y2=E1-v2εyy+vεxx,(5c)τxy-e0a2∂2σxy∂x2+∂2σxy∂y2=Gγxy,whereE,G, andv are the elastic modulus, the shear modulus, and Poisson’s ratio of the SLGS, respectively. The internal characteristic lengtha corresponds to the distance between two atoms in a C-C bond and is equal to 0.142 nm. Flexural moments of SLGS are obtained using the following expressions:(6)Mxx=∫-h/2h/2zσxxdz,Myy=∫-h/2h/2zσyydz,Mxy=∫-h/2h/2zτxydz.If in-plane displacements of the middle surface are neglected in thex andy direction then the relationship between strain and displacement fields is expressed as (7)εxx=-z∂2w∂x2,εyy=-z∂2w∂y2,γxy=-2z∂2w∂x∂y,where w is displacement in direction ofz axis of GS. With substitution of (5a), (5b), and (5c) into (6) and using (7) the following expressions are obtained:(8a)Mxx-e0a2∂2Mxx∂x2+∂2Mxx∂y2=-D∂2w∂x2+v∂2w∂y2(8b)Myy-e0a2∂2Myy∂x2+∂2Myy∂y2=-D∂2w∂y2+v∂2w∂x2,(8c)Mxy-e0a2∂2Mxy∂x2+∂2Mxy∂y2=-D1-vG∂2w∂x∂y,whereD is the flexural rigidity of SLGS, expressed as (9)D=Eh3121-v2.The governing equation that describes the flexural vibration of SLGS carrying a nanoparticle at arbitrary position(x0,y0) can be written in the following form: (10)∂2Mxx∂x2+2∂2Mxy∂x∂y+∂2Myy∂y2=ρh+mcδx-x0δy-y0∂2w∂t2,where ρ,mc,t, and δ are the mass density of SLGS, mass of nanoparticle, time, and the Dirac delta function, respectively.With substitution of (8a), (8b), and (8c) into (10), the following governing equation is obtained: (11)D∂4w∂x4+2∂4w∂x2∂y2+∂4w∂y4+1-e0a2∂2w∂x2+∂2w∂y2ρh+mcδx-x0δy-y0∂2w∂t2=0.The harmonic solution of previous equation can be expressed as (12)wx,y,t=Yx,yeiωt,where Yx,y represents a shape function of deflection and ω is the resonant frequency of the SLGS.Substituting (12) into (11), the governing equation can be expressed as (13)∂4Y∂x4+2∂4Y∂x2∂y2+∂4Y∂y4-ω2D1-e0a2∂2Y∂x2+∂2Y∂y2ρh+mcδx-x0δx-y0∂2Y∂t2=0.Since the GS is simply supported, the following boundary conditions are applied:(14)w=0,∂2w∂x2=0,∂2w∂y2=0onx=0,La,andy=0,Lb.The shape function from (12) can be expressed as (15)Yx,y=Amnsin⁡mπxLasin⁡nπyLb,where parameter Amn represents the vibration amplitude and indexesm andn represent mode numbers in the periodic directions, respectively.In order to obtain the equation for determining frequencies of SLGS, (15) is substituted in (13). After the substitution, the equation is multiplied by x=sin⁡(mπx/La)sin⁡(nπy/Lb) and integrated over the whole region with respect tox andy in limitsx = 0 toLa andy = 0 toLb. With some simplifications the following frequency equation is obtained:(16)∫0Lb∫0LaAmnDπ4m2La2+n2Lb22sin2⁡mπxLasin2⁡nπyLbdxdy-ω2∫0Lb∫0LaAmn1+e0a2π2m2La2+n2Lb2ρh+mcδx-x0δy-y0sin2⁡mπxLasin2⁡nπyLbdxdy=0.The roots of the previous equation are the desired resonant frequencies and each resonant frequency corresponds to a given shape function. After further manipulation of the previous equation, the expression for determining resonant frequency is obtained in the following form:(17)ωmn2=Dπ4m2/La2+n2/Lb221+e0a2π2m2/La2+n2/Lb2ρh+4mc/LaLbsin2⁡mπξsin2⁡nπη,where ξ=x0/La and η=y0/Lb are the nondimensional parameters of nanoparticle position. In (17) in denominator sin2⁡mπξsin2⁡nπη is multiplied by the nonlocal parameter (e0a) so the nonlocal parameter is influencing the mass localization on SLGS. The frequency is obtained from (18)fmn=ωmn2π.It is worth mentioning that when the nonlocal parametere0a is equal to zero, the resonant frequency of SLGS with attached nanoparticles collapses to the classical plate theory.


## 2.2. Double-Layer Graphene Sheets

DLGS are composed of two SLGS and they interact with each other by vdW forces. In these simulations, the nanoparticle is attached on the upper layer while no nanoparticles are attached to the lower level of GS. The DLGS is assumed to be simply supported at the boundary. The DLGS system with attached nanoparticle at arbitrary position is shown in Figure2.Figure 2
Schematic illustration of DLGS with attached mass at arbitrary position.The governing equation that describes vibration of this system can be written in the following form:(19a)D∇4w1+1-e0a2∇2ρh+mcδx-x0δy-y0∂2w1∂t2=p12,(19b)D∇4w2+ρh1-e0a2∇2∂2w2∂t2=p21,where wjx,y,t,j=1,2 are the flexural deflections of upper j=1 and lower j=2 sheet and p12p21=-p12 is the transverse pressure between layers of DLGS caused by vdW forces. As in the case of SLGS the variableD is flexural rigidity given by (9) while ∇4 is the biharmonic operator and in the Cartesian coordinate system it can be written in the following form: (20)∇4.=∂4∂x4+∂4∂x2∂y2+∂4∂y4.The distributed transverse pressure acting on the upper and lower layer of DLGS system can be written as(21)p12=cw2-w1,where c is the vdW interaction coefficient between upper and lower layers [4, 34, 35, 40–42]. It can be obtained from the Lennard-Jones potential as (22)c=439a224εσ2σa83003π256∑k=05-1k2k+15kσa61z¯1-z¯212-35π8∑k=02-1k2k+12k1z¯1-z¯26,where ε=2.968meV and σ=0.34nm are parameters that are chosen to fit the physical properties of graphene sheets. The parameter z-j can be determined from the following equation: (23)z-j=zja,j=1,2,where zj is the coordinate ofj-th layer in the direction of the thickness with the origin in the mid-plane of the graphene sheets, and parameter a=0.142nm is the C-C bond length.The second term in (19a) is given by (24)μx,y=ρh+mcδx-x0δy-y0,where μx,y is the variable mass distribution function of the upper sheet.As in the case of SLGS, to obtain the resonant frequencies of DLGS system described by governing equations (19a) and (19b) the following function is introduced:(25)wjx,y,t=Yjx,yeiωt,j=1,2,where Yjx,y,j=1,2 represents a shape function of deflection in the upper and the lower sheet of DLGS system and ω is the resonant frequency.(26)wjx,y,t=Yjx,yeiωt,j=1,2,Substituting (26) into differential equations (19a) and (19b) of DLGS system the following equations are obtained in the matrix form: (27)D∇4-μω2e0a2∇2+c-μω2-c-cD∇4-ρhω2e0a2∇2+c-ρhω2Y1Y2=00.With further algebraic manipulation of the previous equation it is obtained that(28)∇8Yj+μ+ρhω2e0a2D2∇6Yj+2c-ρhω2-μω2D+μρhω4e0a4D2∇4Yj+cρh+μ-2ρhω2ω2e0a2D2∇2Yj+ρhμω4-cρh+μω2D2Yj=0,j=1,2.Again, the DLGS system is simply supported which means that the following boundary conditions are applied: (29)w=0,∂2w∂x2=0,∂2w∂y2=0onx=0,La,andy=0,Lb.To fulfill previously defined boundary conditions, the following shape function of deflection is proposed: (30)Yjx,y=Amnsin⁡mπxLasin⁡nπyLb,j=1,2,where Amn is the vibrational amplitude of oscillations, while indexesm andn indicate the mode numbers. With substitution of (30) into (28), the following equation is obtained: (31)π8Lb2m2+La2n24La8Lb8Yj-π6Lb2m2+La2n23La6Lb6μ+hρω2e0a2DYj+π4Lb2m2+La2n22La4Lb4D2c-ω2μ+hρ+ρhμω4e0a2D2Yj+π2Lb2m2+La2n2La2Lb2e0a2ω22hμρω2-cμ+ρhD2Yj+ω2hμρω2-cμ+ρhD2Yj=0.Next step is to substitute (24) into the previous equation:(32)
π
8
L
b
2
m
2
+
L
a
2
n
2
4
L
a
8
L
b
8
Y
j
-
π
6
L
b
2
m
2
+
L
a
2
n
2
3
L
a
6
L
b
6
h
ρ
ω
2
e
0
a
2
D
Y
j
+
π
4
L
b
2
m
2
+
L
a
2
n
2
2
L
a
4
L
b
4
2
c
-
ω
2
ρ
h
D
Y
j
+
π
2
L
b
2
m
2
+
L
a
2
n
2
L
a
2
L
b
2
e
0
a
2
ω
2
ρ
h
D
2
Y
j
+
ω
2
c
ρ
h
D
2
Y
j
=
π
6
L
b
2
m
2
+
L
a
2
n
2
3
L
a
6
L
b
6
ω
2
e
0
a
2
D
+
π
4
L
b
2
m
2
+
L
a
2
n
2
2
L
a
4
L
b
4
ω
2
-
ρ
h
ω
4
e
0
a
4
D
2
-
π
2
L
b
2
m
2
+
L
a
2
n
2
L
a
2
L
b
2
e
0
a
2
ω
2
2
h
ρ
ω
2
-
c
D
2
Y
j
-
ω
2
h
ρ
ω
2
-
c
D
2
m
c
δ
x
-
x
0
δ
y
-
y
0
Y
j
.After multiplication of the previous equation with the term sin⁡mπx/Lasin⁡nπy/Lb from both sides and with integration over the whole region with respect tox andy the following solution is obtained: (33)r0ω4+r1ω2+r2=0,where the coefficients r0,r1, and r2 are (34a)r0=18D2La4Lb4hLa2Lb2+n2La2+m2Lb2e0a2π22ρ2mc+hLaLbρ+mc-2cos⁡2nπη-2cos⁡2mπξ+cos⁡2πnη-mξ+cos⁡2πnη+mξ,(34b)r1=-12D2La6Lb6cijLa4Lb4+Dn2La2+m2Lb22π4m2Lb2e0a2π2+La2Lb2+n2e0a2π2hLaLbρ+2mcsin⁡nπη2sin⁡mπξ2,(34c)r2=n2La2+m2Lb22π42cijLa4Lb4+Dn2La2+m2Lb22π44DLa7Lb7.ξ and η, as in the SLGS case, represent the nondimensional parameter of the attached nanoparticle position. r0 and r1 depend on the nonlocal parameter so the nonlocal parameter influences natural frequencies as well. The solution of (33) gives a term for determining natural frequencies of DLGS system with attached nanoparticle and can be expressed as (35)ωmn2=-r1+r12-4r0r22r0,And, finally, the frequency is obtained from(36)fmn=ωmn2π.


## 3. Results and Discussion

As already described, the considered nanomechanical resonator consists of simply supported SLGS or DLGS with the attached nanoparticle. Unless otherwise stated, the dimensionless parametersξ and ηare equal to 0.5. This means that the nanoparticle is at the center of the graphene sheet. Regarding mechanical properties of the graphene sheet in general Young’s modulus is not yet precisely determined. To this day researchers used different approaches such as density function theory (DFT) [43], molecular dynamics [44, 45] or MD, continuum mechanics approach [46, 47], finite element modeling [48], and experimental investigation [49] in order to determine its value. All these analyses gave Young’s modulus of graphene sheet around 1 TPa which will be used as Young’s modulus in our simulations. Geometrical and mechanical data of the simply supported graphene sheet is shown in Table 1.Table 1
Geometrical and mechanical characteristics of SLGS.
Description Value Length [nm] 30  Width [nm] 20  Thickness [nm] 0.345  Internal characteristic length [nm] 0.142  Modulus of elasticity [TPa] 1.00  Density [kg/m3] 2300  Poisson’s number 0.149The attached nanoparticles in these analyses are CWA molecules and the data for each CWA molecule is given in Table2.Table 2
Technical data of CWA molecules [39].
Name Molecular Formula Average Mass [Da = g/mol ] Mass of one molecule [zg] Type of Chemical Warfare Agent: Blister agent  Sulfur Mustard C4H8Cl2S 159.077 264.16  Nitrogen Mustard HN-1 C6H13Cl2N 170.08 282.431  Nitrogen Mustard HN-2 C5H11Cl2N 156.053 259.138  Nitrogen Mustard HN-3 C6H12Cl3N 204.53 337.977  Lewisite C2H2AsCl3 207.318 344.268  Type of Chemical Warfare Agent: Nerve gases  G series  Sarin C4H10FO2P 140.093 232.635  Tabun C5H11N2O2P 162.127 262.224  Soman C7H16FO2P 182.173 302.512  Cyclosarin C7H14FO2P 180.16 299.17  V series  VE C10H24NO2Ps 253.342 420.694  VG C10H24NO3Ps 269.341 447.261  VM C10H22NO2Ps 239.315 397.401  VX C11H26NO2Ps 267.368 443.985  Type of Chemical Warfare Agent: Asphyxiants  Cyanogen chloride CClN 61.470 102.076  Hydrogen Cyanide CHN 27.025 44.8771  Arsine H3As 77.9454 129.434  Type of Chemical Warfare Agent: Choking Agent  Chlorine Cl2 70.906 117.745  Chloropicrin CCl3NO2 164.375 272.957  Phosgene CCl2O 98.916 164.258  Type of Chemical Warfare Agent: Behavioral Agents/Incapacitating Agents  3-Quinuclidinyl benzilate (BZ) C21H23NO3 337.41 560.296


### 3.1. Influence of Different Parameters on Frequencies of SLGS

In order to investigate how an attached CWA molecule affects resonant frequencies of SLGS, influence of parameters such as the nonlocal parameter, geometry, and mass of attached nanoparticle needs to be investigated. The influence of the nonlocal parameter on natural frequencies will be investigated first. For the geometrical and mechanical parameters listed in Table1 the influence of the nonlocal parameter on natural frequencies of simply supported SLGS for different mode shapes is shown in Figure 3.Figure 3
Variation of the natural frequency of SLGS versus nonlocal parameter for different mode shapes.As seen in Figure3, when the value of the nonlocal parameter corresponds to 0 the frequency of all the mode shapes is at its highest. In other words, the highest natural frequencies are obtained in the case of local solutions. As the value of the nonlocal parameter increases the natural frequency starts to decrease. Looking at the range of the nonlocal parameter from 0 to 20 nm, it is obvious that the value of the nonlocal parameter is strongly influencing the resonant frequency. In range from 20 to 50 nm, the resonant frequency drop is much slower. In order to achieve natural frequencies of SLGS as in experimental investigations that are lower than 1 GHz without considering any dissipative mechanisms the value of the nonlocal parameter must be greater than 20 nm. Since the focus of this paper is not to determine the true value of the nonlocal parameter but to examine the possibility of using graphene in CWA molecule detection, the nonlocal parameter in further analyses will range from 1 to 4 nm.It should be emphasized that experimental investigations [19, 50, 51] showed that resonant frequencies of graphene sheets are lower (MHz range) when compared to results obtained using the nonlocal theory of elasticity (GHz range) [1, 4, 31, 34]. This difference in frequency is due to the fact that the effect of all dissipative mechanisms on the experimental system is not known and because the exact value of the nonlocal parameter is still unknown, which represents the main disadvantage of the nonlocal theory for smaller values of the nonlocal parameter.In order to investigate influence of geometry on the resonant frequency the mechanical characteristics from Table1 are taken into account and for this analysis the value of nonlocal parameter is 1 nm and values of nondimensional parameters of shape function ofm andn are either 1 or 2.As seen from Figure4 when the lengthsLa andLb are in range from 5 to 20 nm then the natural frequencies are in range of 200 to 20 GHz. As the length of graphene sheet gets larger the frequency starts to drop down and is in range from 0 to 2 GHz.Figure 4
Variation of the natural frequency of the simply supported SLGS versus geometry (e0a = 1 nm).To investigate the influence of the attached nanoparticle mass on resonant frequencies of the simply supported SLGS, the particle is placed at the center of SLGS,ξ=0.5,η=0.5. Since the mass of CWA molecules is in range from 0 to 600 zg (Table 2), the mass of attached nanoparticle is varied in that range, respectively. For this simulation the mechanical and geometrical data are again taken from Table 1 and results of simulation for different mode shapes are shown in Figure 5.Figure 5
Influence of attached mass on resonant frequencies of SLGS (ξ=0.5,η=0.5).From Figure5 it is obvious that the nonlocal parameter has greater influence on the resonant frequency than the attached nanoparticle. Analyzed mode shapes points out that if the mass of attached nanoparticle and the nonlocal parameter is low the resonant frequency is around 8 - 30 GHz depending on the mode shape. As the mass of the attached nanoparticle and the nonlocal parameter gets larger the resonant frequency drops down to a MHz range which is closer to experimentally obtained results [19, 50, 51].


### 3.2. Influence of Different Parameters on Frequencies of DLGS

For analysis of influence of different parameters on resonant frequencies of the simply supported DLGS the mechanical and geometrical parameters are used from Table1. Firstly, the influence of the nonlocal parameter on natural frequencies is investigated. In this analysis the value of nonlocal parameter ranges from 0 to 50 nm, the mass of attached nanoparticle is equal to 0 and the nondimensional shape function parametersm andn are set in range from 1 to 4. The result of calculation is shown in Figure 6.Figure 6
Natural frequencies of DLGS versus nonlocal parameter for different mode shapes.As expected, the natural frequency of DLGS system is highest when the value of nonlocal parameter is equal to 0 nm, i.e., for the local formulation. From Figure6 it is obvious that resonant frequency starts to drop as the nonlocal parameter increases.The next step is to investigate the influence of geometry on natural frequencies of DLGS. For this analysis the nonlocal parameter is set to 1 nm and mechanical parameters are taken from Table1. Results of the performed analysis are shown in Figure 7.Figure 7
Natural frequency of DLGS versus its dimensions.In Figure7 the influence of geometry of simply supported DLGS system on natural frequency is shown. For each mode shape when the valuesLa andLb are small (in range from 5 to 40 nm) the natural frequency of DLGS system is in range from 20 to 600 GHz depending on the mode shape. WhenLa andLb are in range from 40 to 160 nm the resonant frequency of DLGS is in range from 0 to 20 GHz. So as the dimensions of DLGS system increase, the frequencies of DLGS system decrease.The next step is to show how the mass of attached nanoparticle and the nonlocal parameter together is influencing the resonant frequency of the simply supported DLGS system. As in the case of simply supported SLGS system the mass of attached nanoparticle is taken in range from 0 to 600 zg since the mass CWA molecules from Table2 are in that range. For this analysis the nonlocal parameter is taken in range from 0 to 100 nm. The position of attached CWA nanoparticle is at the center of upper plate of DLGS (ξ=0.5,η=0.5). Nondimensional values of shape function ofm andn are set to 1 and 2, respectively. The result of the analysis is shown in Figure 8.Figure 8
Resonant frequency of DLGS versus mass of the attached nanoparticle and nonlocal parameter (e0a=1nm,ξ=0.5,η=0.5).From Figure8 it is easy to conclude that mass of attached nanoparticle and nonlocal parameter lowers the resonant frequency of simply supported DLGS system. As in the case of simply supported SLGS the nonlocal parameter has greater influence on resonant frequency than the mass of attached nanoparticle.


### 3.3. Absolute and Relative Frequency Shift of Single-Layer and Double-Layer Graphene Sheets

In order to detect the attached nanoparticle by the absolute frequency shift method first the resonant frequencies of the SLGS or DLGS system without attached nanoparticle must be determined and then resonant frequencies with attached nanoparticle must be determined, respectively. Resonant frequencies of the free system and system with attached nanoparticle have to be determined since the frequency shift is defined as the difference between the resonant frequency of a SLGS (or DLGS) with and without attached nanoparticles. Thus, it can be written in the following form:(37)Δf=f-f0,where f0 is the resonant frequency of SLGS or DLGS without attached nanoparticle and f the resonant frequency of SLGS or DLGS with attached nanoparticle. According to [32] relative frequency shift can be written in the following form:(38)frelative=Δff0=ff0-1.From (38) it can be concluded that the nonlocal parameter has no influence on the relative frequency shift. In particular, since the term involving the nonlocal parameter in denominator 1+e0a2π2m2/La2+n2/Lb2 (17) can be canceled in division of frequencies with and without attached mass f/f0, the relative frequency shift has to be independent of the nonlocal parameter. The same conclusion can be made for DLGS case, eq. (36). The importance of this observation is raised once more later in this section by means graphical representation.


#### 3.3.1. Frequency Shift for Simply Supported SLGS System

As described above, frequency shifts for SLGS for different mode shapes are calculated and the results are provided in Figure9.Figure 9
Variation of the frequency shift of simply supported SLGS versus attached mass for different mode shapes (ξ=0.5,η=0.5,e0a=1.0nm).From Figure9 it can be seen that the mass of the attached CWA is not influencing all mode shapes. Attached CWA mass affects only first mode shape (m = 1,n = 1), fifth mode shape (m = 3,n = 1), and seventh mode shape (m = 3,n = 3). In these mode shapes as the mass of attached nanoparticle increases the frequency shift increases too. The lack of influence of the frequency shift on other analyzed mode shapes can be explained by the position of attached mass. Since the attached nanoparticle is at the center of simply supported SLGS the molecule is located at the nodal line. The nodal line is the line which remains stationary during the vibration cycle. Since the attached nanoparticle is located on this nodal line there is no change in the resonant frequency or in the frequency shift. Figure 10 shows the difference between the mode shape without nodal line and mode shape with one nodal line.Figure 10
Two different mode shapes of simply supported SLGS: (a) mode shape without nodal lines (m = 1,n =1) and (b) mode shape with one nodal line (m = 2,n = 1).
(a) (b)In Figure10 two mode shapes are shown. The first mode shape (m = 1,n = 1) does not have nodal lines so the attached CWA molecule is influencing the resonant frequency. In second mode shape (m = 2,n = 1) the attached CWA molecule is located on the nodal line which means that this molecule is not influencing the resonant frequency since the nodal line is not moving during the vibration cycle.Since the mode shapes that reacted specific position of attached CWA molecule are found, the next step is to see how the nonlocal parameter is influencing the frequency shift of simply supported SLGS. The results of this analysis are shown in Figures11-12.Figure 11
Variation of frequency shift of simply supported SLGS as a function of attached mass for different nonlocal parameters (m=1,n=1,ξ=0.5,η=0.5).Figure 12
Variation of frequency shift of simply supported SLGS as a function of attached mass for different nonlocal parameters (m=3,n=1,ξ=0.5,η=0.5).As seen from Figures11-12 the nonlocal parameter lowers the absolute frequency shift. When mass of nanoparticle is low the difference between resonant frequencies of SLGS without nanoparticle and SLGS with the attached nanoparticle is negligible. If a mass of attached nanoparticle becomes larger the difference in resonant frequencies between SLGS with and without attached nanoparticle becomes larger.


#### 3.3.2. Influence of Attached CWA Molecule Location on the Frequency Shift for SLGS System

So far, the location of attached nanoparticle or CWA molecule was at a center of the simply supported SLGS system. The results of performed simulations showed that SLGS could detect a CWA molecule based on frequency shift method only for specific mode shapes. In this analysis the influence of attached CWA molecule location on frequency shift will be investigated for simply supported SLGS. This analysis was performed for mode shape with nondimensional mode shape parametersm = 2 andn =1. This mode shape is a mode shape with one nodal line and the goal is to show that other mode shapes are also good for detection of CWA molecules but only when the CWA molecule is not attached on the nodal line. For this analysis two different CWA molecules are chosen, that is, hydrogen cyanide (CHN) and 3-quinuclidinyl benzilate (BZ). The mass of these two molecules represent the boundaries in the mass range of CWA molecules. Looking at Table 2 the CHN has lowest mass value of one molecule while BZ has the highest mass value. The results of the analysis are shown in Figures 13 and 14.Figure 13
Influence of the nondimensional location parametersξ,η on the frequency shift for the case of attached hydrogen cyanide on the simply supported SLGS.Figure 14
Influence of the nondimensional location parametersξ,η on the frequency shift for the case of attached 3-quinuclidinyl benzilate (BZ) on the simply supported SLGS.From Figures13 and 14 it can be concluded that nondimensional location parameters also have the influence on the frequency shift and resonant frequency of simply supported SLGS system. As the values of nondimensional location parameters increase the value of frequency shift also increases. The specific mode shape used in this analysis showed that it can be used for detection of CWA molecules but not when these molecules are attached to nodal lines. As seen from Figures 13 and 14 the highest frequency shift is at the center of SLGS (ξ=0.5,η=0.5).


#### 3.3.3. Frequency Shift of Simply Supported DLGS System

The next step is to determine if a DLGS system could be possible candidate for nanoparticle detection or detection of CWA molecules based on the absolute and relative frequency shift calculation. Like in case of simply supported SLGS, the mode shapes useful for detection of attached CWA molecules have to be identified. Results of the analysis are shown in Figure15.Figure 15
Variation of the frequency shift as a function of the attached nanoparticle mass on DLGS system for different mode shapes (e0a=1nm,ξ=0.5,η=0.5).As seen in Figure15 the frequency shift can be obtained for specific mode shapes of simply supported DLGS system and these are first mode shape (m = 1,n = 1), fifth mode shape (m = 3,n = 1), and seventh mode shape (m = 3,n = 3). As in the case of simply supported SLGS the reason why some mode shapes could not detect the attached CWA molecule is because the attached CWA molecule is located on the nodal line. If the location of CWA molecule was changed then other mode shapes could be used for its detection. Now, mode shapes of the simply supported DLGS system that could detect the CWA molecules are found. The influence of the nonlocal parameter on the frequency shift is shown in Figures 16-17.Figure 16
Variation of the frequency shift of the simply supported DLGS as a function of attached mass for different nonlocal parameters (m=1,n=1,ξ=0.5,η=0.5).Figure 17
Variation of the frequency shift of the simply supported DLGS as a function of attached mass for different nonlocal parameters (m=3,n=1,ξ=0.5,η=0.5).As in the case of simply supported SLGS system, the nonlocal parameter lowers the frequency shift for analyzed mode shapes of simply supported DLGS system. From previous figures it can also be concluded that for each of the following mode shapes the nonlocal parameter has greater influence on values of the frequency shift. For example, in first mode shape the influence of nonlocal parameter on the frequency shift is very small since the distance between the frequency shift curves is very small, Figure16. Looking at Figure 17 it can be concluded that the distance between each of the frequency shift curves is larger when compared to the first mode shape so the influence of the nonlocal parameter is not the same for different mode shapes.


### 3.4. Detection of Chemical Warfare Agents with Single-Layer and Double-Layer Graphene Sheets Based on the Frequency Shift Method

In this analysis the simply supported SLGS is modelled based on the mechanical and geometrical properties in Table1. Mass of the attached nanoparticle is now replaced with the CWA molecules which are presented in Table 2. The results of the analysis are shown in Figures 18-19.Figure 18
Frequency shift of simply supported SLGS versus mass, attached at SLGS center (ξ=0.5,η=0.5,m=1,n=1).Figure 19
Frequency shift of simply supported SLGS versus mass, attached at SLGS center (ξ=0.5,η=0.5,m=3,n=1).From Figures18-19 it can be concluded that simply supported SLGS system is capable of detecting all of the CWA molecules given in Table 2 by the application of the frequency shift method.The simply supported DLGS system is also modeled based on the mechanical and geometrical data in Table1. The CWA molecule is attached at the upper plate of DLGS system at its center. The analysis is conducted for first (m = 1,n = 1) and fifth (m = 3,n =1) mode shape since they are sensitive to the frequency shift. The results of the performed analysis are shown in Figures 20-21.Figure 20
Frequency shift of the simply supported DLGS versus mass of attached mass at the center of upper plate (ξ=0.5,η=0.5,m=1,n=1).Figure 21
Frequency shift of the simply supported DLGS versus mass of attached mass at the center of upper plate (ξ=0.5,η=0.5,m=3,n=1).Based on results shown in Figures20-21 the specific mode shapes of simply supported DLGS system can detect the CWA molecules given in Table 2 based on the frequency shift method.


### 3.5. Detection of Chemical Warfare Agents with Single-Layer and Double-Layer Graphene Sheets Based on Relative Frequency Shift Method

In this analysis the possibility of the SLGS and DLGS system in detection of CWA will be investigated using relative frequency shift method. As in the case of frequency shift (absolute) analysis the geometrical and material data of the SLGS and DLGS system are taken from Table2. The position of CWA molecule is at center of SLGS and center of upper plate of DLGS. The investigation will be conducted only for the first mode shape (m = 1,n = 1) since the values of relative frequency shift are the same for other two modes and there is no influence of the nonlocal parameter on relative frequency shift as described before. The results of analysis are shown in Figures 22-23.Figure 22
Relative frequency shift of simply supported SLGS versus mass, attached at SLGS center (ξ=0.5,η=0.5,m=1,n=1).Figure 23
Relative frequency shift of simply supported SLGS versus mass, attached at DLGS center (ξ=0.5,η=0.5,m=1,n=1).In Figures22 and 23 the variation of relative frequency is shown as a function of attached nanoparticle at the center of the SLGS and upper DLGS plate. These figures also show that all CWA molecules line on the relative frequency shift curve which means they could be detected. The results obtained for the absolute and relative frequency shift showed that SLGS and DLGS plate can detect CWA molecules.


## 3.1. Influence of Different Parameters on Frequencies of SLGS

In order to investigate how an attached CWA molecule affects resonant frequencies of SLGS, influence of parameters such as the nonlocal parameter, geometry, and mass of attached nanoparticle needs to be investigated. The influence of the nonlocal parameter on natural frequencies will be investigated first. For the geometrical and mechanical parameters listed in Table1 the influence of the nonlocal parameter on natural frequencies of simply supported SLGS for different mode shapes is shown in Figure 3.Figure 3
Variation of the natural frequency of SLGS versus nonlocal parameter for different mode shapes.As seen in Figure3, when the value of the nonlocal parameter corresponds to 0 the frequency of all the mode shapes is at its highest. In other words, the highest natural frequencies are obtained in the case of local solutions. As the value of the nonlocal parameter increases the natural frequency starts to decrease. Looking at the range of the nonlocal parameter from 0 to 20 nm, it is obvious that the value of the nonlocal parameter is strongly influencing the resonant frequency. In range from 20 to 50 nm, the resonant frequency drop is much slower. In order to achieve natural frequencies of SLGS as in experimental investigations that are lower than 1 GHz without considering any dissipative mechanisms the value of the nonlocal parameter must be greater than 20 nm. Since the focus of this paper is not to determine the true value of the nonlocal parameter but to examine the possibility of using graphene in CWA molecule detection, the nonlocal parameter in further analyses will range from 1 to 4 nm.It should be emphasized that experimental investigations [19, 50, 51] showed that resonant frequencies of graphene sheets are lower (MHz range) when compared to results obtained using the nonlocal theory of elasticity (GHz range) [1, 4, 31, 34]. This difference in frequency is due to the fact that the effect of all dissipative mechanisms on the experimental system is not known and because the exact value of the nonlocal parameter is still unknown, which represents the main disadvantage of the nonlocal theory for smaller values of the nonlocal parameter.In order to investigate influence of geometry on the resonant frequency the mechanical characteristics from Table1 are taken into account and for this analysis the value of nonlocal parameter is 1 nm and values of nondimensional parameters of shape function ofm andn are either 1 or 2.As seen from Figure4 when the lengthsLa andLb are in range from 5 to 20 nm then the natural frequencies are in range of 200 to 20 GHz. As the length of graphene sheet gets larger the frequency starts to drop down and is in range from 0 to 2 GHz.Figure 4
Variation of the natural frequency of the simply supported SLGS versus geometry (e0a = 1 nm).To investigate the influence of the attached nanoparticle mass on resonant frequencies of the simply supported SLGS, the particle is placed at the center of SLGS,ξ=0.5,η=0.5. Since the mass of CWA molecules is in range from 0 to 600 zg (Table 2), the mass of attached nanoparticle is varied in that range, respectively. For this simulation the mechanical and geometrical data are again taken from Table 1 and results of simulation for different mode shapes are shown in Figure 5.Figure 5
Influence of attached mass on resonant frequencies of SLGS (ξ=0.5,η=0.5).From Figure5 it is obvious that the nonlocal parameter has greater influence on the resonant frequency than the attached nanoparticle. Analyzed mode shapes points out that if the mass of attached nanoparticle and the nonlocal parameter is low the resonant frequency is around 8 - 30 GHz depending on the mode shape. As the mass of the attached nanoparticle and the nonlocal parameter gets larger the resonant frequency drops down to a MHz range which is closer to experimentally obtained results [19, 50, 51].


## 3.2. Influence of Different Parameters on Frequencies of DLGS

For analysis of influence of different parameters on resonant frequencies of the simply supported DLGS the mechanical and geometrical parameters are used from Table1. Firstly, the influence of the nonlocal parameter on natural frequencies is investigated. In this analysis the value of nonlocal parameter ranges from 0 to 50 nm, the mass of attached nanoparticle is equal to 0 and the nondimensional shape function parametersm andn are set in range from 1 to 4. The result of calculation is shown in Figure 6.Figure 6
Natural frequencies of DLGS versus nonlocal parameter for different mode shapes.As expected, the natural frequency of DLGS system is highest when the value of nonlocal parameter is equal to 0 nm, i.e., for the local formulation. From Figure6 it is obvious that resonant frequency starts to drop as the nonlocal parameter increases.The next step is to investigate the influence of geometry on natural frequencies of DLGS. For this analysis the nonlocal parameter is set to 1 nm and mechanical parameters are taken from Table1. Results of the performed analysis are shown in Figure 7.Figure 7
Natural frequency of DLGS versus its dimensions.In Figure7 the influence of geometry of simply supported DLGS system on natural frequency is shown. For each mode shape when the valuesLa andLb are small (in range from 5 to 40 nm) the natural frequency of DLGS system is in range from 20 to 600 GHz depending on the mode shape. WhenLa andLb are in range from 40 to 160 nm the resonant frequency of DLGS is in range from 0 to 20 GHz. So as the dimensions of DLGS system increase, the frequencies of DLGS system decrease.The next step is to show how the mass of attached nanoparticle and the nonlocal parameter together is influencing the resonant frequency of the simply supported DLGS system. As in the case of simply supported SLGS system the mass of attached nanoparticle is taken in range from 0 to 600 zg since the mass CWA molecules from Table2 are in that range. For this analysis the nonlocal parameter is taken in range from 0 to 100 nm. The position of attached CWA nanoparticle is at the center of upper plate of DLGS (ξ=0.5,η=0.5). Nondimensional values of shape function ofm andn are set to 1 and 2, respectively. The result of the analysis is shown in Figure 8.Figure 8
Resonant frequency of DLGS versus mass of the attached nanoparticle and nonlocal parameter (e0a=1nm,ξ=0.5,η=0.5).From Figure8 it is easy to conclude that mass of attached nanoparticle and nonlocal parameter lowers the resonant frequency of simply supported DLGS system. As in the case of simply supported SLGS the nonlocal parameter has greater influence on resonant frequency than the mass of attached nanoparticle.


## 3.3. Absolute and Relative Frequency Shift of Single-Layer and Double-Layer Graphene Sheets

In order to detect the attached nanoparticle by the absolute frequency shift method first the resonant frequencies of the SLGS or DLGS system without attached nanoparticle must be determined and then resonant frequencies with attached nanoparticle must be determined, respectively. Resonant frequencies of the free system and system with attached nanoparticle have to be determined since the frequency shift is defined as the difference between the resonant frequency of a SLGS (or DLGS) with and without attached nanoparticles. Thus, it can be written in the following form:(37)Δf=f-f0,where f0 is the resonant frequency of SLGS or DLGS without attached nanoparticle and f the resonant frequency of SLGS or DLGS with attached nanoparticle. According to [32] relative frequency shift can be written in the following form:(38)frelative=Δff0=ff0-1.From (38) it can be concluded that the nonlocal parameter has no influence on the relative frequency shift. In particular, since the term involving the nonlocal parameter in denominator 1+e0a2π2m2/La2+n2/Lb2 (17) can be canceled in division of frequencies with and without attached mass f/f0, the relative frequency shift has to be independent of the nonlocal parameter. The same conclusion can be made for DLGS case, eq. (36). The importance of this observation is raised once more later in this section by means graphical representation.


### 3.3.1. Frequency Shift for Simply Supported SLGS System

As described above, frequency shifts for SLGS for different mode shapes are calculated and the results are provided in Figure9.Figure 9
Variation of the frequency shift of simply supported SLGS versus attached mass for different mode shapes (ξ=0.5,η=0.5,e0a=1.0nm).From Figure9 it can be seen that the mass of the attached CWA is not influencing all mode shapes. Attached CWA mass affects only first mode shape (m = 1,n = 1), fifth mode shape (m = 3,n = 1), and seventh mode shape (m = 3,n = 3). In these mode shapes as the mass of attached nanoparticle increases the frequency shift increases too. The lack of influence of the frequency shift on other analyzed mode shapes can be explained by the position of attached mass. Since the attached nanoparticle is at the center of simply supported SLGS the molecule is located at the nodal line. The nodal line is the line which remains stationary during the vibration cycle. Since the attached nanoparticle is located on this nodal line there is no change in the resonant frequency or in the frequency shift. Figure 10 shows the difference between the mode shape without nodal line and mode shape with one nodal line.Figure 10
Two different mode shapes of simply supported SLGS: (a) mode shape without nodal lines (m = 1,n =1) and (b) mode shape with one nodal line (m = 2,n = 1).
(a) (b)In Figure10 two mode shapes are shown. The first mode shape (m = 1,n = 1) does not have nodal lines so the attached CWA molecule is influencing the resonant frequency. In second mode shape (m = 2,n = 1) the attached CWA molecule is located on the nodal line which means that this molecule is not influencing the resonant frequency since the nodal line is not moving during the vibration cycle.Since the mode shapes that reacted specific position of attached CWA molecule are found, the next step is to see how the nonlocal parameter is influencing the frequency shift of simply supported SLGS. The results of this analysis are shown in Figures11-12.Figure 11
Variation of frequency shift of simply supported SLGS as a function of attached mass for different nonlocal parameters (m=1,n=1,ξ=0.5,η=0.5).Figure 12
Variation of frequency shift of simply supported SLGS as a function of attached mass for different nonlocal parameters (m=3,n=1,ξ=0.5,η=0.5).As seen from Figures11-12 the nonlocal parameter lowers the absolute frequency shift. When mass of nanoparticle is low the difference between resonant frequencies of SLGS without nanoparticle and SLGS with the attached nanoparticle is negligible. If a mass of attached nanoparticle becomes larger the difference in resonant frequencies between SLGS with and without attached nanoparticle becomes larger.


### 3.3.2. Influence of Attached CWA Molecule Location on the Frequency Shift for SLGS System

So far, the location of attached nanoparticle or CWA molecule was at a center of the simply supported SLGS system. The results of performed simulations showed that SLGS could detect a CWA molecule based on frequency shift method only for specific mode shapes. In this analysis the influence of attached CWA molecule location on frequency shift will be investigated for simply supported SLGS. This analysis was performed for mode shape with nondimensional mode shape parametersm = 2 andn =1. This mode shape is a mode shape with one nodal line and the goal is to show that other mode shapes are also good for detection of CWA molecules but only when the CWA molecule is not attached on the nodal line. For this analysis two different CWA molecules are chosen, that is, hydrogen cyanide (CHN) and 3-quinuclidinyl benzilate (BZ). The mass of these two molecules represent the boundaries in the mass range of CWA molecules. Looking at Table 2 the CHN has lowest mass value of one molecule while BZ has the highest mass value. The results of the analysis are shown in Figures 13 and 14.Figure 13
Influence of the nondimensional location parametersξ,η on the frequency shift for the case of attached hydrogen cyanide on the simply supported SLGS.Figure 14
Influence of the nondimensional location parametersξ,η on the frequency shift for the case of attached 3-quinuclidinyl benzilate (BZ) on the simply supported SLGS.From Figures13 and 14 it can be concluded that nondimensional location parameters also have the influence on the frequency shift and resonant frequency of simply supported SLGS system. As the values of nondimensional location parameters increase the value of frequency shift also increases. The specific mode shape used in this analysis showed that it can be used for detection of CWA molecules but not when these molecules are attached to nodal lines. As seen from Figures 13 and 14 the highest frequency shift is at the center of SLGS (ξ=0.5,η=0.5).


### 3.3.3. Frequency Shift of Simply Supported DLGS System

The next step is to determine if a DLGS system could be possible candidate for nanoparticle detection or detection of CWA molecules based on the absolute and relative frequency shift calculation. Like in case of simply supported SLGS, the mode shapes useful for detection of attached CWA molecules have to be identified. Results of the analysis are shown in Figure15.Figure 15
Variation of the frequency shift as a function of the attached nanoparticle mass on DLGS system for different mode shapes (e0a=1nm,ξ=0.5,η=0.5).As seen in Figure15 the frequency shift can be obtained for specific mode shapes of simply supported DLGS system and these are first mode shape (m = 1,n = 1), fifth mode shape (m = 3,n = 1), and seventh mode shape (m = 3,n = 3). As in the case of simply supported SLGS the reason why some mode shapes could not detect the attached CWA molecule is because the attached CWA molecule is located on the nodal line. If the location of CWA molecule was changed then other mode shapes could be used for its detection. Now, mode shapes of the simply supported DLGS system that could detect the CWA molecules are found. The influence of the nonlocal parameter on the frequency shift is shown in Figures 16-17.Figure 16
Variation of the frequency shift of the simply supported DLGS as a function of attached mass for different nonlocal parameters (m=1,n=1,ξ=0.5,η=0.5).Figure 17
Variation of the frequency shift of the simply supported DLGS as a function of attached mass for different nonlocal parameters (m=3,n=1,ξ=0.5,η=0.5).As in the case of simply supported SLGS system, the nonlocal parameter lowers the frequency shift for analyzed mode shapes of simply supported DLGS system. From previous figures it can also be concluded that for each of the following mode shapes the nonlocal parameter has greater influence on values of the frequency shift. For example, in first mode shape the influence of nonlocal parameter on the frequency shift is very small since the distance between the frequency shift curves is very small, Figure16. Looking at Figure 17 it can be concluded that the distance between each of the frequency shift curves is larger when compared to the first mode shape so the influence of the nonlocal parameter is not the same for different mode shapes.


## 3.3.1. Frequency Shift for Simply Supported SLGS System

As described above, frequency shifts for SLGS for different mode shapes are calculated and the results are provided in Figure9.Figure 9
Variation of the frequency shift of simply supported SLGS versus attached mass for different mode shapes (ξ=0.5,η=0.5,e0a=1.0nm).From Figure9 it can be seen that the mass of the attached CWA is not influencing all mode shapes. Attached CWA mass affects only first mode shape (m = 1,n = 1), fifth mode shape (m = 3,n = 1), and seventh mode shape (m = 3,n = 3). In these mode shapes as the mass of attached nanoparticle increases the frequency shift increases too. The lack of influence of the frequency shift on other analyzed mode shapes can be explained by the position of attached mass. Since the attached nanoparticle is at the center of simply supported SLGS the molecule is located at the nodal line. The nodal line is the line which remains stationary during the vibration cycle. Since the attached nanoparticle is located on this nodal line there is no change in the resonant frequency or in the frequency shift. Figure 10 shows the difference between the mode shape without nodal line and mode shape with one nodal line.Figure 10
Two different mode shapes of simply supported SLGS: (a) mode shape without nodal lines (m = 1,n =1) and (b) mode shape with one nodal line (m = 2,n = 1).
(a) (b)In Figure10 two mode shapes are shown. The first mode shape (m = 1,n = 1) does not have nodal lines so the attached CWA molecule is influencing the resonant frequency. In second mode shape (m = 2,n = 1) the attached CWA molecule is located on the nodal line which means that this molecule is not influencing the resonant frequency since the nodal line is not moving during the vibration cycle.Since the mode shapes that reacted specific position of attached CWA molecule are found, the next step is to see how the nonlocal parameter is influencing the frequency shift of simply supported SLGS. The results of this analysis are shown in Figures11-12.Figure 11
Variation of frequency shift of simply supported SLGS as a function of attached mass for different nonlocal parameters (m=1,n=1,ξ=0.5,η=0.5).Figure 12
Variation of frequency shift of simply supported SLGS as a function of attached mass for different nonlocal parameters (m=3,n=1,ξ=0.5,η=0.5).As seen from Figures11-12 the nonlocal parameter lowers the absolute frequency shift. When mass of nanoparticle is low the difference between resonant frequencies of SLGS without nanoparticle and SLGS with the attached nanoparticle is negligible. If a mass of attached nanoparticle becomes larger the difference in resonant frequencies between SLGS with and without attached nanoparticle becomes larger.


## 3.3.2. Influence of Attached CWA Molecule Location on the Frequency Shift for SLGS System

So far, the location of attached nanoparticle or CWA molecule was at a center of the simply supported SLGS system. The results of performed simulations showed that SLGS could detect a CWA molecule based on frequency shift method only for specific mode shapes. In this analysis the influence of attached CWA molecule location on frequency shift will be investigated for simply supported SLGS. This analysis was performed for mode shape with nondimensional mode shape parametersm = 2 andn =1. This mode shape is a mode shape with one nodal line and the goal is to show that other mode shapes are also good for detection of CWA molecules but only when the CWA molecule is not attached on the nodal line. For this analysis two different CWA molecules are chosen, that is, hydrogen cyanide (CHN) and 3-quinuclidinyl benzilate (BZ). The mass of these two molecules represent the boundaries in the mass range of CWA molecules. Looking at Table 2 the CHN has lowest mass value of one molecule while BZ has the highest mass value. The results of the analysis are shown in Figures 13 and 14.Figure 13
Influence of the nondimensional location parametersξ,η on the frequency shift for the case of attached hydrogen cyanide on the simply supported SLGS.Figure 14
Influence of the nondimensional location parametersξ,η on the frequency shift for the case of attached 3-quinuclidinyl benzilate (BZ) on the simply supported SLGS.From Figures13 and 14 it can be concluded that nondimensional location parameters also have the influence on the frequency shift and resonant frequency of simply supported SLGS system. As the values of nondimensional location parameters increase the value of frequency shift also increases. The specific mode shape used in this analysis showed that it can be used for detection of CWA molecules but not when these molecules are attached to nodal lines. As seen from Figures 13 and 14 the highest frequency shift is at the center of SLGS (ξ=0.5,η=0.5).


## 3.3.3. Frequency Shift of Simply Supported DLGS System

The next step is to determine if a DLGS system could be possible candidate for nanoparticle detection or detection of CWA molecules based on the absolute and relative frequency shift calculation. Like in case of simply supported SLGS, the mode shapes useful for detection of attached CWA molecules have to be identified. Results of the analysis are shown in Figure15.Figure 15
Variation of the frequency shift as a function of the attached nanoparticle mass on DLGS system for different mode shapes (e0a=1nm,ξ=0.5,η=0.5).As seen in Figure15 the frequency shift can be obtained for specific mode shapes of simply supported DLGS system and these are first mode shape (m = 1,n = 1), fifth mode shape (m = 3,n = 1), and seventh mode shape (m = 3,n = 3). As in the case of simply supported SLGS the reason why some mode shapes could not detect the attached CWA molecule is because the attached CWA molecule is located on the nodal line. If the location of CWA molecule was changed then other mode shapes could be used for its detection. Now, mode shapes of the simply supported DLGS system that could detect the CWA molecules are found. The influence of the nonlocal parameter on the frequency shift is shown in Figures 16-17.Figure 16
Variation of the frequency shift of the simply supported DLGS as a function of attached mass for different nonlocal parameters (m=1,n=1,ξ=0.5,η=0.5).Figure 17
Variation of the frequency shift of the simply supported DLGS as a function of attached mass for different nonlocal parameters (m=3,n=1,ξ=0.5,η=0.5).As in the case of simply supported SLGS system, the nonlocal parameter lowers the frequency shift for analyzed mode shapes of simply supported DLGS system. From previous figures it can also be concluded that for each of the following mode shapes the nonlocal parameter has greater influence on values of the frequency shift. For example, in first mode shape the influence of nonlocal parameter on the frequency shift is very small since the distance between the frequency shift curves is very small, Figure16. Looking at Figure 17 it can be concluded that the distance between each of the frequency shift curves is larger when compared to the first mode shape so the influence of the nonlocal parameter is not the same for different mode shapes.


## 3.4. Detection of Chemical Warfare Agents with Single-Layer and Double-Layer Graphene Sheets Based on the Frequency Shift Method

In this analysis the simply supported SLGS is modelled based on the mechanical and geometrical properties in Table1. Mass of the attached nanoparticle is now replaced with the CWA molecules which are presented in Table 2. The results of the analysis are shown in Figures 18-19.Figure 18
Frequency shift of simply supported SLGS versus mass, attached at SLGS center (ξ=0.5,η=0.5,m=1,n=1).Figure 19
Frequency shift of simply supported SLGS versus mass, attached at SLGS center (ξ=0.5,η=0.5,m=3,n=1).From Figures18-19 it can be concluded that simply supported SLGS system is capable of detecting all of the CWA molecules given in Table 2 by the application of the frequency shift method.The simply supported DLGS system is also modeled based on the mechanical and geometrical data in Table1. The CWA molecule is attached at the upper plate of DLGS system at its center. The analysis is conducted for first (m = 1,n = 1) and fifth (m = 3,n =1) mode shape since they are sensitive to the frequency shift. The results of the performed analysis are shown in Figures 20-21.Figure 20
Frequency shift of the simply supported DLGS versus mass of attached mass at the center of upper plate (ξ=0.5,η=0.5,m=1,n=1).Figure 21
Frequency shift of the simply supported DLGS versus mass of attached mass at the center of upper plate (ξ=0.5,η=0.5,m=3,n=1).Based on results shown in Figures20-21 the specific mode shapes of simply supported DLGS system can detect the CWA molecules given in Table 2 based on the frequency shift method.


## 3.5. Detection of Chemical Warfare Agents with Single-Layer and Double-Layer Graphene Sheets Based on Relative Frequency Shift Method

In this analysis the possibility of the SLGS and DLGS system in detection of CWA will be investigated using relative frequency shift method. As in the case of frequency shift (absolute) analysis the geometrical and material data of the SLGS and DLGS system are taken from Table2. The position of CWA molecule is at center of SLGS and center of upper plate of DLGS. The investigation will be conducted only for the first mode shape (m = 1,n = 1) since the values of relative frequency shift are the same for other two modes and there is no influence of the nonlocal parameter on relative frequency shift as described before. The results of analysis are shown in Figures 22-23.Figure 22
Relative frequency shift of simply supported SLGS versus mass, attached at SLGS center (ξ=0.5,η=0.5,m=1,n=1).Figure 23
Relative frequency shift of simply supported SLGS versus mass, attached at DLGS center (ξ=0.5,η=0.5,m=1,n=1).In Figures22 and 23 the variation of relative frequency is shown as a function of attached nanoparticle at the center of the SLGS and upper DLGS plate. These figures also show that all CWA molecules line on the relative frequency shift curve which means they could be detected. The results obtained for the absolute and relative frequency shift showed that SLGS and DLGS plate can detect CWA molecules.


## 4. Conclusions

In this paper, the resonant frequencies of simply supported SLGS and DLGS were analyzed by the nonlocal plate theory. The nonlocal plate theory, being an extension of the classical plate theory (Kirchhoff) can be used to tackle problems at nanoscale. The results showed the following.(i)
The nonlocal parameter has a huge influence on resonant frequencies of simply supported SLGS and DLGS. As demonstrated, maximal natural frequencies are obtained for the local solutions. When the nonlocal parameter gets larger than zero, the natural frequency is lowered. The nonlocal parameter generally tends to lower natural or resonant frequencies of SLGS and DLGS for every mode shape.(ii)
For small dimensions of SLGS or DLGS the natural or resonant frequencies in MHz range (obtained in experimental investigations [19, 50, 51]) could be obtained by the nonlocal theory in the absence of dissipation mechanisms, only if the nonlocal parameter is greater than 20 nm.(iii)
For the large values of the nonlocal parameter and the attached nanoparticle mass the resonant frequencies of the SLGS and DLGS systems are in the MHz range as in [19, 50, 51],(iv)
The influence of geometry of simply supported SLGS or DLGS system showed that as dimensions of plates increase the frequency drops to a certain value and after that it remains constant.(v)
The investigation also showed that the mass of attached nanoparticle has huge influence on resonant frequencies of simply supported SLGS and DLGS. As the mass of the attached nanoparticle increases the resonant frequencies become lower.The central investigation in this paper was calculation of the frequency shift. In order to build fully functional NEMS where GS is a building block for detection of CWA molecules, the resonant frequency of SLGS or DLGS with and without CWA molecules needs to be determined. The results showed the following.(i)
Based on the difference of resonant frequencies free SLGS or DLGS and resonant frequency of SLGS or DLGS with attached CWA, it can be determined which CWA molecule is attached to the system.(ii)
The position of attached CWA on SLGS and DLGS system is also important. In this paper the position of CWA molecule was at the center of simply supported SLGS and DLGS system and the analysis showed that only certain mode shapes could detect CWA molecules. The main reason for this anomaly stems from the stationarity of nodal lines during the vibration cycle, so the mass of CWA molecule has no influence on the resonant frequency or the frequency shift of system. So, in order to use SLGS or DLGS system as a building block for NEMS mass resonators analysis of mode shapes which are capable of detection of CWA molecules must be conducted and the frequency shift method must be applied. This will provide information about the specific CWA molecule which is attached to GS. In this paper the position of CWA molecules on frequency shift method was also investigated and results showed that other nonactive mode shapes could be utilized. It is shown that both SLGS and DLGS systems can be highly effective in such detection.(iii)
The nonlocal parameter does not affect relative frequency shift regardless of mode shape. The result showed that the mode shapes that showed change in absolute frequency shift versus added mass have equal relative frequency changes.



---

*Source: 6451861-2019-06-17.xml*