Effect of Mass on the Dynamic Characteristics of Single- and Double-Layered Graphene-Based Nano Resonators

Graphene has been widely and extensively used in mass sensing applications. The present study focused on exploring the use of single-layer graphene (SLG) and double-layer graphene (DLG) as sensing devices. The dynamic analysis of SLG and DLG with different boundary conditions (BDs) and length was executed using the atomistic finite element method (AFEM). SLG and DLG sheets were modelled and considered as a space–frame structure similar to a 3D beam. Spring elements (Combin14) were used to identify the interlayer interactions between two graphene layers in the DLG sheet due to the van der Waals forces. Simulations were carried out to visualize the behavior of the SLG and DLG subjected to different BDs and when used as mass sensing devices. The variation in frequency was noted by changing the length and applied mass of the SLGs and DLGs. The quantity of the frequency was found to be highest in the armchair SLG (6, 6) for a 50 nm sheet length and lowest in the chiral SLG (16, 4) for a 20 nm sheet length in the bridged condition. When the mass was 0.1 Zg, the frequency for the zigzag SLG (20, 0) was higher in both cases. The results show that the length of the sheet and the various mass values have a significant impact on the dynamic properties. The present research will contribute to the ultra-high frequency nano-resonance applications.


Introduction
Graphene was discovered in 2004, and its field of research now deals with various aspects of graphene. Graphene possesses amazing features, such as high electrical and thermal conductivity, and this drives the research [1][2][3][4]. It is a suitable material for flexible electronics and nanomechanical systems due to its higher stiffness, reduced mass per unit area, strength, and improved electrical conductivity nature [5][6][7][8][9]. A one-atom-thick planar sheet of sp2-bonded carbon atoms densely packed in a two-dimensional honeycomb crystal lattice is referred to as SLG. DLG is a stack of graphene sheets separated by 0.34 nanometres. Nanorings are made up of distorted graphite sheets, and they are the basic structural constituent of these structures. As one of the most widely used methods for producing graphene sheets, electrochemical exfoliation was the first synthetic process. Graphene sheets are thin sheets of graphene [10][11][12][13]. Microbial detections and diagnosis devices are all examples of nanodevices with outstanding properties in engineering and medicine [14][15][16][17]. The theoretical modelling of nanostructured materials can be conducted in three ways: (a) atomistic modelling, (b) hybrid atomistic-continuum mechanics, and (c) continuum mechanics, in addition to experimental studies. Classic molecular dynamics (MD) is an example of the atomic modelling approach [18]. MD simulations for the free vibration of multiple SLGs with varied lengths were utilized to find the suitable values of the nonlocal parameter. The focus of this research is to investigate the variation in the frequency of SLGs and DLGs for different values of the mass and length of the sheets using an atomistic model and MD simulations [19]. The resonant frequencies of the graphene sheets, the effects of the sheet sizes, the boundary conditions, and the number of layers were investigated. The vibrations of DLGs with various BDs between two layers were also considered. The modelling of the SLGs and DLGs is conducted at the atomic level [20]. The carbon atoms of the graphene sheets are covalently linked together to form a hexagonal lattice. Individual atoms cannot move due to the external forces and extreme bonds. Graphene can be described as a space-frame structure by seeing the bonds as connecting load-carrying elements and the atoms as the joints of the connecting elements. The mechanical behavior of graphene can be studied using standard structural mechanics methods if it is considered as a space-frame structure [21][22][23][24][25][26][27]. MD is an efficient technique for modeling graphene's complete mechanical performance. However, it has a significant computing cost, which could be very costly for large-scale problems, particularly in the case of vibration analysis of graphene single-and double-layer sheets. The present work focused on the vibration analysis of SLG and DLG using zeptogram mass. For both designs, the results show a decrease in resonance frequencies as the associated mass increases.

Van Der Waal Interaction in Graphene
At distances greater than 1 nm, the VdW force is significant in the interactions of atoms and molecules with carbon nanostructures [28]. Lifshitz and Pitaevskii developed a comprehensive VdW force theoretical framework based on the resonance dielectric permittivity principle, which is based on quantum statistical physics [29]. The intermolecular and atomic interaction with a microbody's flat surface is described by this theory. Figure 1 shows two layers of the double-layer graphene sheet; the spring element Combin 14 is used to represent the weak VdW force and the atomic mass at the node. The FEM model of the DLG is indicated in Figure 2.
the effects of the sheet sizes, the boundary conditions, and the number of layers we vestigated. The vibrations of DLGs with various BDs between two layers were als sidered. The modelling of the SLGs and DLGs is conducted at the atomic level [20 carbon atoms of the graphene sheets are covalently linked together to form a hexa lattice. Individual atoms cannot move due to the external forces and extreme bonds phene can be described as a space-frame structure by seeing the bonds as connecting carrying elements and the atoms as the joints of the connecting elements. The mech behavior of graphene can be studied using standard structural mechanics methods considered as a space-frame structure [21][22][23][24][25][26][27]. MD is an efficient technique for mod graphene's complete mechanical performance. However, it has a significant comp cost, which could be very costly for large-scale problems, particularly in the case of tion analysis of graphene single-and double-layer sheets. The present work focus the vibration analysis of SLG and DLG using zeptogram mass. For both designs, t sults show a decrease in resonance frequencies as the associated mass increases.

Van Der Waal Interaction in Graphene
At distances greater than 1 nm, the VdW force is significant in the interactio atoms and molecules with carbon nanostructures. [28]. Lifshitz and Pitaevskii deve a comprehensive VdW force theoretical framework based on the resonance dielectri mittivity principle, which is based on quantum statistical physics [29]. The intermol and atomic interaction with a microbody's flat surface is described by this theory. Figure 1 shows two layers of the double-layer graphene sheet; the spring ele Combin 14 is used to represent the weak VdW force and the atomic mass at the nod FEM model of the DLG is indicated in Figure 2.

Interatomic Modelling of SLG and DLG
The beam elements in the atomistic model accurately represent the covalent interactions between the carbon atoms. The elastic characteristics of the beam can be calculated by establishing links between the molecular and structural potential energy [30].
At low stresses, the total steric potential energy of a graphene can be calculated by adding the energies owing to the valence of the bonded contacts or the connected and non-bonded interactions [31].
In the present study, potential energy has been used to analyse linear nanospring stiffness by applying the finite element method. The sum of the forces exerted by the electrons and the electrostatic forces exerted between the positively charged nuclei equals the overall force exerted on each atomic nucleus.
where Ur is the bond stretch interaction energy, Uθ is the bending (bond angle variation energy), Uθ is the dihedral angle torsion energy, Uω is the out-of-plane torsion energy, and Uvdw is the non-bonded van der Waals interaction energy.
where kr, kθ, and k τ are the bond stretching, bond bending, and torsional resistance force constants, respectively, while Δr, Δθ and ΔФ represent the bond stretching increment, the bond angle variation, and the angle variation of the bond twisting, respectively. The second derivatives of the potential energy terms in Equations (2)-(4) with respect to bond

Interatomic Modelling of SLG and DLG
The beam elements in the atomistic model accurately represent the covalent interactions between the carbon atoms. The elastic characteristics of the beam can be calculated by establishing links between the molecular and structural potential energy [30].
At low stresses, the total steric potential energy of a graphene can be calculated by adding the energies owing to the valence of the bonded contacts or the connected and non-bonded interactions [31].
In the present study, potential energy has been used to analyse linear nanospring stiffness by applying the finite element method. The sum of the forces exerted by the electrons and the electrostatic forces exerted between the positively charged nuclei equals the overall force exerted on each atomic nucleus.
where U r is the bond stretch interaction energy, U θ is the bending (bond angle variation energy), U φ is the dihedral angle torsion energy, U ω is the out-of-plane torsion energy, and U vdw is the non-bonded van der Waals interaction energy.
where k r , k θ , and k τ are the bond stretching, bond bending, and torsional resistance force constants, respectively, while ∆r, ∆θ and ∆φ represent the bond stretching increment, the bond angle variation, and the angle variation of the bond twisting, respectively. The second derivatives of the potential energy terms in Equations (2)-(4) with respect to bond length, bond angle, and twisting bond angle variations produce the spring stiffness coefficients k r , k θ , and k τ according to Castigliano's theorem. The elements representing the bond are assumed to be elastic beams with Young's modulus E, length L, cross-sectional area A, and moment of inertia I.
The strain energy under pure tension N is given by The strain energy of the beam element under the pure bending moment M can be expressed as: where α is the rotational angle of the beam ends. Similarly, the strain energy of the beam element under pure twisting moment T is given by where J, G, and ∆β are the polar moment of inertia, the shear modulus, and the relative rotations of the beam ends, respectively [20].
The elastic properties of the beam element are given as [21] where d, L, E, and G represent the diameter, length, Young's modulus, and shear modulus of the beam element. The stiffness K s of the special spring is defined by the following equation: For the Finite Element Model of DLG, as indicated in Figure 2, the DLGs are modeled considering the weak van der Waals force of attraction between the upper and lower layers as a spring element. The van der Waals force field between the interfacial layers is represented by the spring element Combin 14.
The xand y-axes are aligned with the SLG's length Lp and width Lq, respectively, and the z-axis is aligned with the SLG's thickness h [50].
The two-dimensional nonlocal constitutive equations of the SLGs are written using nonlocal elasticity theory, as follows: σ yy − (e 0a ) 2 σ 2 yy σx 2 + σ 2 yy where E, G, and ν are the elastic modulus, the shear modulus, and the Poisson's ratio of the GSs, respectively. The internal characteristic length a is the distance between two atoms in a C-C bond, which is 0.142 nm.  The x-and y-axes are aligned with the SLG's length Lp and width Lq, respectively, and the z-axis is aligned with the SLG's thickness h [50].
The two-dimensional nonlocal constitutive equations of the SLGs are written using nonlocal elasticity theory, as follows: where E, G, and ν are the elastic modulus, the shear modulus, and the Poisson's ratio of the GSs, respectively. The internal characteristic length a is the distance between two atoms in a C-C bond, which is 0.142 nm.
The link between the strain and the displacement fields is represented as being when the middle surface displacements in the x and y directions are ignored. The link between the strain and the displacement fields is represented as being when the middle surface displacements in the x and y directions are ignored. (13) where w is the displacement along the GS' thickness. We obtain Equation (15) by substituting Equation (13) into Equation (14) where D is the flexural rigidity of SLGs, expressed as The governing equation for the flexural vibration of SLGs carrying a nanoparticle can be given as where m c is the mass of nanoparticles connected at point (x 0 , y 0 ) and is the Dirac delta function indicated by Substituting Equation (17) into Equation (18), the governing equation can be written in terms of w as The harmonic solution of Equation (19) can be expressed as where Y (x, y) is the shape function of deflection, and ω is the resonant frequency of the SLGs. Substituting Equation (20) into Equation (19), we obtain Note that the boundary conditions of the SLGs with simply supported edges are On x = 0, L p and y = 0 L q . Therefore, the shape function (Y) in Equation (20) can be expressed as where A mn denotes the oscillation's vibration amplitude, and m and n denote the mode numbers in the periodic directions. After some optimizations, we obtain the following frequency equation by substituting the shape function of Equation (23) into Equation (21), multiplying both sides of Equation (17) by sinmx/L p sinny/L q , and integrating over the entire region with respect to x and y, with the limits x = 0 to x = L p and y = 0 to y The required resonant frequencies corresponding to a given form function are all roots of the above equation. For the non-trivial solution, the A mn coefficient should be zero. As a result, the resonant frequency of a mass sensor may be calculated.
The resonance frequency of an SLG with connected nanoparticles may also be calculated using classical elasticity theory when the nonlocal parameter (e 0a ) is considered to be zero. Graphene is an enormously robust nanomaterial, with a Young's module near to 1 TPa [51-53].

Validation of the Model
We obtained an analytical technique from Natsuki et al. [50] to compare with the FEA model with the objective of validating the existing model results. As shown in Table 1, the analytical and FEA results are in close proximity. The following material attributes were considered by the authors when validating the SLG model: E = 3792.47739 nN/nm 2 , ν = 0.17, and ρ = 8362.714 kg/m 3 . The results are presented in Figure 4. The frequency variation between the analytical model and the simulated FEM proximity has been identified. However, the current model for the application as a mass sensing device is reflected in the results of the FEA technique and analytical results.
The results are presented in Figure 4. The frequency variation between the analytical model and the simulated FEM proximity has been identified. However, the current model for the application as a mass sensing device is reflected in the results of the FEA technique and analytical results.

Results and Discussion
The SLG and DLG frame structures were modelled. Nine different varieties of SLGs and DLGs were evaluated for analysis, including three different types of armchairs (8,8), (6,6), (10,10); zigzag (20, 0); (16, 0), (12, 0); and chiral (16,4), (12,4), (8,4). The influence of the diverse boundary conditions, such as cantilever and bridged, on the mass sensing ability of the DLG-and SLG-based mass sensors was investigated. Table 2 shows the fundamental frequency of several types of cantilever DLGs with various mass additions at the tip of sheet. Table 3 demonstrates the frequency of several types of bridged DLGs with various mass additions in the sheet's centre point. When the

Results and Discussion
The SLG and DLG frame structures were modelled. Nine different varieties of SLGs and DLGs were evaluated for analysis, including three different types of armchairs (8,8), (6,6), (10,10); zigzag (20, 0); (16, 0), (12, 0); and chiral (16,4), (12,4), (8,4). The influence of the diverse boundary conditions, such as cantilever and bridged, on the mass sensing ability of the DLG-and SLG-based mass sensors was investigated. Table 2 shows the fundamental frequency of several types of cantilever DLGs with various mass additions at the tip of sheet. Table 3 demonstrates the frequency of several types of bridged DLGs with various mass additions in the sheet's centre point. When the mass is 0.1 Zg, the frequency for the zigzag DLG (20, 0) is higher in both cases. Furthermore, when the associated mass on the tip and centre of the DLG increases, the fundamental frequency decreases. Table 4 shows the fundamental frequency for several cantilever SLGs with various mass additions at the sheet's tip. Table 5 demonstrates the frequency of several types of bridged SLGs with various mass additions at the centre of the sheet. When the mass is 0.1 Zg, the frequency for the zigzag SLG (20, 0) is higher in both cases. Moreover, as the associated mass on the tip and centre of SLGs increases, the fundamental frequency decreases.    6 3.51 × 10 6 1.14 × 10 6 3.60 × 10 5 1.14 × 10 5 The frequency of the cantilever SLG fluctuates with length and chirality, as shown in Figure 5. The quantity of the frequency was found to be larger in the zigzag SLG (6, 6) for a sheet length of 50 nm and lower in the armchair SLG (20 0) for a sheet length of 50 nm in the cantilever condition. The quantity of frequency fluctuates for the different lengths and chirality of the bridged SLGs; this is shown in Figure 6. The quantity of frequency was found to be highest in the armchair SLG (6, 6) for a 50 nm sheet length and lowest in the chiral SLG (16, 4) for a 20 nm sheet length in the bridged condition. The graph clearly shows that as chirality increases, the frequency drops. In a smaller graphene, a higher frequency can be attained. In the bridged state of single-layer graphene, there is a huge variation in frequency. It shows that the bridged condition has a higher frequency. The graph in Figure 6 shows that a considerable fluctuation in frequency is attained at different lengths of the double-layered graphene sheet. As the chirality of a frequency grows, its value eventually falls. With a smaller graphene sheet, the armchair cantilever SLG (6,6), and a sheet length of 40, we obtain the highest frequency value. Its frequency begins to decrease as its chirality increases. The greatest frequency is attained for the smallest cantilever armchair DLG (6,6). DLG (20, 0) achieves the lowest frequency value at a length of 20 nm. The graph in Figure 6 shows that a considerable fluctuation in frequency is attained at different lengths of the double-layered graphene sheet. As the chirality of a frequency grows, its value eventually falls. With a smaller graphene sheet, the armchair cantilever SLG (6,6), and a sheet length of 40, we obtain the highest frequency value. Its frequency begins to decrease as its chirality increases. The greatest frequency is attained for the smallest cantilever armchair DLG (6,6). DLG (20, 0) achieves the lowest frequency value at a length of 20 nm. The bridged armchair DLG (6, 6) with a length of 50 nm achieves the highest frequency. We can see that as the chirality increases, the value of the frequency drops in the bridged-conditioned double-layer graphene as well. The minimum frequency for DLG (20,0) at length 20 nm is obtained under the bridged state.
We discovered the highest frequency value for bridged chiral DLG (12,4) at length 50 nm and mass attached 1 × 10 −20 using the frequency shift graph given in Figure 7. For the shortest length of 10 nm of bridged DLG, we obtain the lowest frequency shift (12,4). As the length of the graph grows, so does the frequency shift.  The graph in Figure 6 shows that a considerable fluctuation in frequency is attained at different lengths of the double-layered graphene sheet. As the chirality of a frequency grows, its value eventually falls. With a smaller graphene sheet, the armchair cantilever SLG (6,6), and a sheet length of 40, we obtain the highest frequency value. Its frequency begins to decrease as its chirality increases. The greatest frequency is attained for the smallest cantilever armchair DLG (6,6). DLG (20, 0) achieves the lowest frequency value at a length of 20 nm. The bridged armchair DLG (6, 6) with a length of 50 nm achieves the highest frequency. We can see that as the chirality increases, the value of the frequency drops in the bridged-conditioned double-layer graphene as well. The minimum frequency for DLG (20,0) at length 20 nm is obtained under the bridged state.
We discovered the highest frequency value for bridged chiral DLG (12,4) at length 50 nm and mass attached 1 × 10 −20 using the frequency shift graph given in Figure 7. For the shortest length of 10 nm of bridged DLG, we obtain the lowest frequency shift (12,4). As the length of the graph grows, so does the frequency shift. The bridged armchair DLG (6, 6) with a length of 50 nm achieves the highest frequency. We can see that as the chirality increases, the value of the frequency drops in the bridgedconditioned double-layer graphene as well. The minimum frequency for DLG (20, 0) at length 20 nm is obtained under the bridged state.
We discovered the highest frequency value for bridged chiral DLG (12,4) at length 50 nm and mass attached 1 × 10 −20 using the frequency shift graph given in Figure 7. For the shortest length of 10 nm of bridged DLG, we obtain the lowest frequency shift (12,4). As the length of the graph grows, so does the frequency shift. Figure 8 shows the frequency shift graph of the armchair single-layer cantilever graphene SLG (8,8). The maximum value frequency was attained for the longest length of 50 nm and the shortest length of 10 nm of a double-layer graphene sheet. The mass value of 1 × 10 −20 was found to give the highest frequency.  Figure 8 shows the frequency shift graph of the armchair single-layer cantilever graphene SLG (8,8). The maximum value frequency was attained for the longest length of 50 nm and the shortest length of 10 nm of a double-layer graphene sheet. The mass value of 1 × 10 −20 was found to give the highest frequency.     Figure 8 shows the frequency shift graph of the armchair single-layer cantilever graphene SLG (8,8). The maximum value frequency was attained for the longest length of 50 nm and the shortest length of 10 nm of a double-layer graphene sheet. The mass value of 1 × 10 −20 was found to give the highest frequency.     (8,8), (16,0), (18,0), and (12,4). It can be shown that at lower mass values, there are far fewer frequency shifts in the SLG, which significantly rises as the attached mass grows. For an attached mass of more than 0.1 zg, the amount of frequency shift is observed to be larger in SLG (8,8) at a length of 50 nm. This clearly indicates that the frequency fluctuation is much less at lower mass values. It suggests that SLGs with a mass sensitivity of 0.1 zg can be impacted by changes in the sheet length.  Figure 10 depicts the frequency shift caused by mass fluctuation at different cantilever SLGs (8,8), (16,0), and (12,4). It can be shown that at lower mass values, there will be far fewer frequency shifts in the SLG, which significantly rises as the attached mass increases. For an attached mass of more than 0.1 zg, the amount of frequency shift is found to be larger in the armchair SLG (8,8) at a length of 50 nm, and the minimum frequency is found in the chiral SLG (12,4). This clearly indicates that the frequency fluctuation is much less at lower mass and chirality values. It suggests that SLGs with a mass sensitivity of 0.1 zg can be impacted by changes in sheet length.   Figure 10 depicts the frequency shift caused by mass fluctuation at different cantilever SLGs (8,8), (16,0), and (12,4). It can be shown that at lower mass values, there will be far fewer frequency shifts in the SLG, which significantly rises as the attached mass increases. For an attached mass of more than 0.1 zg, the amount of frequency shift is found to be larger in the armchair SLG (8,8) at a length of 50 nm, and the minimum frequency is found in the chiral SLG (12,4). This clearly indicates that the frequency fluctuation is much less at lower mass and chirality values. It suggests that SLGs with a mass sensitivity of 0.1 zg can be impacted by changes in sheet length.  Figure 10 depicts the frequency shift caused by mass fluctuation at different cantilever SLGs (8,8), (16,0), and (12,4). It can be shown that at lower mass values, there will be far fewer frequency shifts in the SLG, which significantly rises as the attached mass increases. For an attached mass of more than 0.1 zg, the amount of frequency shift is found to be larger in the armchair SLG (8,8) at a length of 50 nm, and the minimum frequency is found in the chiral SLG (12,4). This clearly indicates that the frequency fluctuation is much less at lower mass and chirality values. It suggests that SLGs with a mass sensitivity of 0.1 zg can be impacted by changes in sheet length. According to the graph in Figure 11, at a length of 20 nm the frequency was highest in the cantilever chiral SLG (16,4). As the length of the sheet grows longer, it steadily diminishes. As demonstrated in the graph, the zigzag SLG (20, 0) has the lowest frequency According to the graph in Figure 11, at a length of 20 nm the frequency was highest in the cantilever chiral SLG (16,4). As the length of the sheet grows longer, it steadily diminishes. As demonstrated in the graph, the zigzag SLG (20, 0) has the lowest frequency at a length of 20 nm.
OR PEER REVIEW 15 of 22 Figure 11. Graph of frequency (Hz) vs. different lengths of chiral, armchair, and zigzag cantilever DLG.
The graph shows that when the chirality grows, the frequency drops progressively. We can determine from this graph that the maximum frequency can be reached with reduced chirality and length. The frequency drops as the mass on the sheet increases. Figure 12 illustrates the frequency shift due to mass fluctuation at various lengths of SLG (8,8), (16,0), (18,0), and (12,4). It can be shown that at lower mass values, there were much smaller frequency shifts in the SLG, which significantly rises as the attached mass grows. For an attached mass of greater than 0.1 zg, the highest amount of frequency shift was discovered in SLG (8,8), at a length of 50 nm. This clearly indicates that the frequency fluctuation is much less at lower mass values. It suggests that SLGs with a mass sensitivity of 0.1 zg can be impacted by changes in the sheet length.  The graph shows that when the chirality grows, the frequency drops progressively. We can determine from this graph that the maximum frequency can be reached with reduced chirality and length. The frequency drops as the mass on the sheet increases. Figure 12 illustrates the frequency shift due to mass fluctuation at various lengths of SLG (8,8), (16,0), (18,0), and (12,4). It can be shown that at lower mass values, there were much smaller frequency shifts in the SLG, which significantly rises as the attached mass grows. For an attached mass of greater than 0.1 zg, the highest amount of frequency shift was discovered in SLG (8,8), at a length of 50 nm. This clearly indicates that the frequency fluctuation is much less at lower mass values. It suggests that SLGs with a mass sensitivity of 0.1 zg can be impacted by changes in the sheet length. The graph shows that when the chirality grows, the frequency drops progressively. We can determine from this graph that the maximum frequency can be reached with reduced chirality and length. The frequency drops as the mass on the sheet increases. Figure 12 illustrates the frequency shift due to mass fluctuation at various lengths of SLG (8,8), (16,0), (18,0), and (12,4). It can be shown that at lower mass values, there were much smaller frequency shifts in the SLG, which significantly rises as the attached mass grows. For an attached mass of greater than 0.1 zg, the highest amount of frequency shift was discovered in SLG (8,8), at a length of 50 nm. This clearly indicates that the frequency fluctuation is much less at lower mass values. It suggests that SLGs with a mass sensitivity of 0.1 zg can be impacted by changes in the sheet length. The graph shows in Figure 13. that in single-layer graphene at maximum length, maximum frequency may be reached. The graph clearly demonstrates that as the length of the sheet rises, the frequency increases as well, with the highest frequency at 50 nm. As The graph shows in Figure 13. that in single-layer graphene at maximum length, maximum frequency may be reached. The graph clearly demonstrates that as the length of the sheet rises, the frequency increases as well, with the highest frequency at 50 nm. As the length lowers and the mass rises, the frequency decreases gradually. As a result, the frequency variation will be much less with higher mass.

ER REVIEW
16 of 22 Figure 13. Frequency vs. mass graph for bridged armchair SLG (8,8). Figure 14 shows the first ten mode shapes of vibration for the bridged SLG (16, 4) with the length of 50 nm. The vibration of the bridged chiral SLG (16,4) rises from the second mode, as can be seen from the mode shape above. Mode 5 shows a half-sine-wave shape of response.  Figure 14 shows the first ten mode shapes of vibration for the bridged SLG (16,4) with the length of 50 nm. The vibration of the bridged chiral SLG (16,4) rises from the second mode, as can be seen from the mode shape above. Mode 5 shows a half-sine-wave shape of response.
FEA results show areas of high stress with red color, lower than that with green and yellow color and lowest stress with blue. MN is stress at the fixed end and MX stress point at the span of layer.
As shown in Figure 15. The vibrations increase with the second mode of the cantilever single-layer graphene chiral graphene (16,4) and continue to increase until the last mode of vibration. A nonlinear pattern is observed for the different mode shapes. Figure 16 shows the first ten mode shapes of bridged DLG (20, 0) with a length of 50 nm. The variation in the vibration of the sheet from the second mode can be seen from the mode shape of the bridged DLG (20, 0). The last mode shows a non-coaxial vibration of response.
The variation in the vibration of the sheet as shown in Figure 17. starts from the second mode, as shown in the mode shape of the cantilever DLG (20, 0). FEA results show areas of high stress with red color, lower than that with green and yellow color and lowest stress with blue. MN is stress at the fixed end and MX stress point at the span of layer.
As shown in Figure 15. The vibrations increase with the second mode of the cantilever single-layer graphene chiral graphene (16,4) and continue to increase until the last mode of vibration. A nonlinear pattern is observed for the different mode shapes.  Figure 16 shows the first ten mode shapes of bridged DLG (20, 0) with a length of 50 nm. The variation in the vibration of the sheet from the second mode can be seen from the mode shape of the bridged DLG (20, 0). The last mode shows a non-coaxial vibration of response.  The variation in the vibration of the sheet as shown in Figure 17. starts from the second mode, as shown in the mode shape of the cantilever DLG (20, 0).

Conclusions
To investigate the vibrational characteristics of the cantilever and bridged doublelayer zigzag, armchair, and chiral graphene sheets, an atomistic FEA model for the van der Waal interaction between the upper and lower graphene sheets and the concentrated

Conclusions
To investigate the vibrational characteristics of the cantilever and bridged double-layer zigzag, armchair, and chiral graphene sheets, an atomistic FEA model for the van der Waal interaction between the upper and lower graphene sheets and the concentrated masses was prepared.
The frequency shift graph of SLG and DLG for the bridging and cantilever conditions concludes that as the length of the sheet increases, so does its frequency, and as the mass of the sheet increases, the frequency also increases. The effect of different lengths of singleand double-layer graphene sheets on frequency was studied using an atomistic modelling technique. It was discovered that the longer the graphene sheet, the higher the frequency. The impacts of mass attached to the tip of the single-layer graphene and the double-layer graphene for the cantilever conditions, as well as to the centre of a bridging graphene sheet, were investigated. In both cases, maximum frequency was reached for mass values greater than 1 zg. The vibrations began with the second mode shape as the SLG and DLG vibrate in different planes, as observed from the mode shapes of the SLGs and DLGs for the cantilever and bridged states. In the bridged DLG (6, 6) with a 1.00 × 10 −16 mass attached at the centre of the sheet, the maximum frequency achieved was 7.76 × 10 6 , and in bridged SLG (6, 6) with a similar mass attached, the highest frequency obtained was 7.50 × 10 5 . The results show that the graphene sheet (6, 6) with a length of 50 nm achieves the maximum frequency for both the bridged SLG and the DLG. This research will help in future applications of graphene in advanced nano-resonator applications.