# Introduction raphene is an outstanding material which has a number of multifunctional properties that repeatedly gross it into the title of "wonder material" which is a road map on the way to guide the community toward the development of products [1].The remarkable mechanical behavior and properties of graphene-based material's has concerned with important study concern in recent years, in line for to their encouraging forecasts, now adaptable divisions for example micromechanics [2], microelectronics [3], and thermal [4] application with desired mechanical properties, and electrical conductivities [2][3][4]. The trial and hypothetical revision of graphene, two-dimensional (2D), is a tremendously growing field of today's condensed matter research. The causes for this massive methodical attention were diverse; on the other hand, one might highlight some key inspirations. Keeping given of the science-based interest generated via graphene and its promising upcoming contribution toward electronic engineering and sensing applications, so a group of research effort is steadfastly hooked on considering the configuration and properties of graphene in this paper. Outstanding toward its excellent mechanical behavior, thermal and electrical conductivities of graphene could also use for more conventional purposes as compared with carbon nanotubes, which was quit restricted to aerospace industries and graphene is also known to have veryhigh stiffness in addition strength until now an extensive scatter have been witnessed in the mechanical properties [1][2][3][4].In this effort, we present molecular dynamics model simulation for the initiation of defects and the influence of different defects(vacancy defects) and pristine one on mechanical strength of graphene sheets were observed and, the fracture strength was predicted from the numerical simulation and the properties of graphene in table 1 and investigated young's modulus displayed in table 2 below. Table 1: Properties for a Single Sheet of Graphene [1]. # Property Value Young's modulus [1] 1.0TPa Rupture strength [1] 130GPa Tensile strength [2] 100GPa Thermal conductivity [3][4] 5000w/mK Shear modulus [5][6] 280GPa Longitudinal sound velocity [5,[7][8][9] 20km/s Melting temperature [5,10] 4900K Specific surface area [11] 2630m 2 /g Optical transmittance [12] 97.70% High electron mobility [13] 250,000cm # Modeling and Methodology a) Molecular Dynamics based Simulation Molecular dynamics-based simulations were performed to study the effect of non-bonded interactions on the mechanical behavior and failure morphology of defective graphene sheet. The success of any molecular dynamics-based simulations entirely depends on the interatomic potentials chosen for simulating the atomic interactions. A Significant amount of advancement in conjunction with computational techniques has already been made by the researchers in developing potentials for capturing the realistic properties for the range of materials. In this study, AIBO (adaptive intermolecular reactive bond order) potential was used to compute the interatomic forces between carbon atoms in graphene. Simulations were performed with a single cutoff distance of 1.95Å as proposed in the work of [25]. AIREBO potential consists of a summation of pair potential REBO (E ij REBO ), non-bonded Lennard Jones potential (E ij LJ ) and torsional component between carbon atoms (E ijk tors ), also described with the help of mathematical expressions in equation (1). , , ,1 2AIREBO REBO LJ ltors ij ij kij i j i k i j l i j k E E E E ? ? ? ? ? = + + ? ? ? ? ? ? ? ?(1) Here, i, j, k, and l refers to individual atoms, E is the total potential energy of the system estimated with the help of AIBO potential. To perform this study, a graphene sheet consisting of 800 atoms was generated in the simulation box. The dimensions of a single sheet of graphene was kept fixed at A=46.599Å and B=49.19Å (as shown in Fig. 1) along the zig-zag and arm chair direction respectively. In-plane periodic boundary conditions were imposed on the simulation box. The interlayer spacing between the sheets of graphene in bilayer graphene was kept constant at 3.45Å. During the simulations, the NPT (isothermalisobaric) ensemble in conjunction with an integration time step of 1fs was enforced. After achieving a minimum energy configuration of graphene, atoms at a temperature of 1K, tensile loading was applied at a strain rate of 0.005ps -1 . To avoid thermal effects on the failure mechanism of graphene, simulations were Global Journal of Researches in Engineering (A ) Volume XIx X Issue III Version I performed at such a low temperature of 1K. Stressstrain response was estimated in this study with the help of the virial stress component [26,27], which can be calculated with the help of mathematical expression given in equation (2). 1, 1 1 ( ); 2 i j ij i j n v f m v v r ? ? ? ? ?? ?? ? ? = = + ? (2) Here, i and j denote indices in Cartesian coordinates system; ?? and ?? are the atomic indices; ?? ?? and ?? ?? are mass and velocity of atom ?;?? ???? is the distance between ?? ?????? ?? atoms and V is the surrounding volume of atom ??. Figure 2: Snapshots on the way to confirm the mathematical method, the fracture strength of a pristine graphene sheet was initially calculated. Stress-strain bends of pristine graphene sheet under same tension along the zigzag way (black color) and armchair way (red color) at 300K. Now the direction of validating the mathematical method, the rupture stress of pure graphene sheet was initially designed. The Consequence of minimal stressstrain bend next to the temperature of 300 K, subjected to tension load alongside both armchair and zigzag directions shown above Fig. 1, was revealed, that fracture stress beside the armchair and the zigzag way are calculated as 91 and 106 GPa, separately. In Cauchy stress; the rupture stiffness was 100GPa and 126 GPa, and the rupture strain is 0.13 and 0.22 correspondingly. These results were promising new examination, i.e., ?? f ?130 GPa and ?? f ?0:25 [28] as well as previous numerical simulation [29], verifying dynamism and exactness of our mathematical approach. Also, graphene can be subjected to a higher temperature at the production stage as well as when graphene-based devices operate at the higher temperature. As we discussed above Chemical vapor deposition (CVD) is one of the most commonly used methods of graphene manufacture; that products graphene at a temperature of around 800 K. Therefore, understanding the temperature behavior of graphene helps to fabricate best excellence graphene founded devices. Studying the effect of high temperature on mechanical properties of a substantial armchair and zigzag is presented. In the temperature range of 200K, 300K, and 450K, the breakage stress with a vacancy III. # Results and Discussion Molecular dynamics-based simulations were performed to capture the failure morphology of pristine graphene either as a single or in bi-layer sheet configuration. These simulations were performed with the help of three models to study the effects of nonbonded interactions on the mechanical behavior of pristine graphene. Stress and strain response estimated along the zig-zag and arm chair directions of pristine single sheet graphene were plotted in Fig. 4. It can be observed from Fig. 4 that the mechanical properties of pristine graphene along with the zig-zag and arm chair directions are quite different because of edge defects. In direction to get a better insight on the failure mechanism of the pristine form of graphene under the influence of tensile loading, snapshots of the simulation box were taken at the time of initiation of the failure as provided in Fig. 4. It is observed that the failure morphology of graphene sheet inferred from the snapshots provided in Fig. 4 is almost independent of the non-bonded interactions. A brittle nature of failure can be observed in zig-zag as well as arm chair directions of graphene sheets under the influence of tensile loading. Stress and strain response estimated along with the zig-zag and arm chair directions of pristine single sheet graphene & bi-layer with (LJ-On) & (LJ-Off) were plotted in Fig. 5 below. It can be observed from Fig. 5 that the mechanical properties of pristine graphene single & bi-layer along the zig zag and arm chair direction. where SG (single graphene sheet), BG (LJ-On) (bilayer graphene sheet with non-bonded interactions) and BG (LJ-Off) (bilayer graphene without non-bonded interactions) [24]. It can be inferred from Fig. 6 that non-bonded interactions as well as stiffness of pristine graphene have an impact on the failure morphology of defective graphene sheet containing single vacancy defects. Snapshots of the simulation box provided in Fig. 6 (c3) for defective graphene sheet accompanied by a pristine sheet of graphene connected with non-bonded interactions showthat the failure initiates at two different regions subsequently and helps in achieving higher failure strength. This initiation of failure at two different defects helps in distributing the energy among these points, which can be attributed to the higher failure strength for defective graphene sheets in bilayer configuration connected with non-bonded interactions. In the way to investigate the reasons behind the improvement in the fracture strength and strain of defective graphene in bilayer configuration of graphene, snapshots at the time of initiation of failure are provided in Fig. 7. It can be observed in Fig. 6 (b3 and c3) that at the higher concentration of single vacancy defects failure triggers from the vacancies at two separate locations. Distribution of loading with the help of nonbonded interactions as well as pristine graphene sheet accompanied the defective graphene can be attributed to the higher strength of defective graphene in bi-layer sheets of graphene. This subsection of the molecular dynamics based simulation helps in concluding that at higher percentage of single vacancy defects, bilayer sheets of graphene shows higher strength and strain values for the failure of defective graphene sheet. Improvement in the strength of defective sheet was observed in the presence of another pristine graphene connected with non-bonded interactions, but no transition from brittle behavior was observed in any of the failure morphology. # a) Result of single, double and multiple vacancy defects Failure morphology of the single graphene with uniformly distributed vacancies during strain failure vs vacancy defect ratio was displayed in Fig. 7. A very instance concentrated stress occurred near unperfected; at that moment breakages happen to open from were vacancy defect started then growth in the (a3) (b3) (c3) defects. direction of nearby defects where fracture starts randomly from the defect of vacancies exist. We now turn to analyze the mechanical properties at the failure point for defective graphene. It should be noted that the ultimate strength is the maximum stress in the stressstrain curves, while the fracture strain is determined from the spontaneous large drop of the total energy increment curves. Without defect, the ultimate tensile strength is 91GPa and 106 GPa intended for armchair and zigzag graphene separately. On behalf of through evenly concentrated defects, the correlation among stress, strains besides defects are revealed below & (b). Obviously, the stress decreases with the increase in vacancy defect, and the strain failure decreases with increase vacancy defect. On or after this we decided that in contrast, stiffness to some extent drops by the rising in vacancy Fig. 8 (b) defect; because lack of an atom implies vacancy defect that graphene is more sensitive to vacancy where carbon bond breakage is happens at the time. This study revealed that fracture stress in zig zag direction with different single, double, and multiple vacancy defects are much better in Pristine single graphene than bilayer di-vacancy, single bilayer vacancy (dangling bond because of odd vacancy defect) and multi-vacancy defect in bilayer single graphene defects are also shown in this bar graph below Fig. 9. # Figure 9: Fracture stress in zig zag direction with different single, double and multiple vacancy defects. Here, the pristine, BG, SG, refers to pristine single graphene sheet, bilayer graphene, and single graphene respectively; Whereas, SV, DV and MV refer to single, di-and multi-vacancy defects. IV. # Conclusions Molecular dynamics-based simulations were performed to predict the effect of non-bonded interactions on the mechanical behavior and failure morphology of defective graphene sheet. Simulations were performed with an isolated defective sheet of graphene or defective sheet of graphene accompanied by a pristine sheet of graphene. Atomistic modeling with single as well as bilayer configuration of graphene was performed with different defect concentrations as well as geometries of vacancy defects such as single, double, and multiple vacancy defects. Di-vacancy defects have predicted higher strength in zig-zag configuration, whereas lower strength in arm chair configuration while compared with the single vacancy defects. A Shift in the failure morphology of graphene along the arm chair direction was observed in bi-layer configuration of defective graphene containing di-vacancy defects. It can be concluded that non-bonded interaction helps in achieving a uniform distribution of load around the defects which triggers the failure simultaneously from different regions & initiating of failure simultaneously from two different points help in achieving a higher strength. # Pristine BG (DV) BG (SV) BG (MV) SG (SV) SG (MV) Fracture Stress (GPa) Graphene Configuration1![Figure 1: Single graphene sheet: same tension alongside; (a) Zigzag direction and, (b) armchair direction; Simulation models of single graphene sheet where the dimension is given by Lx and Ly: uniaxial tension along.](image-2.png "Figure 1 :") ![Journal of Researches in Engineering (A ) Volume XIx X Issue III Version I was evaluated subjected to load; along with armchair and zigzag way. Modeling was held at a temperature of 200K, 300K, and 450K as shown in Fig.3. (a), and (b) Shows the fracture strength ?? f for the graphene sheets without defect at temperatures of 200 K, 300 K, and 450 K.](image-3.png "Global") 3![Figure 3: Stress versus strain curves; (a) Zigzag direction graphene sheets at different temperatures, (b) Armchair direction graphene sheets at different temperatures.](image-4.png "Figure 3 :") 4![Figure 4: Failure morphology of pristine single and bilayer sheets (a &d) along zig-zag and arm chair directions respectively (b & e) along zig-zag and arm chair directions with non-bonded interactions between bi-layer configuration of graphene containing one defective and another pristine sheet (c & f) along zig-zag and arm chair directions with (bi-layer configuration of graphene containing one defective and another pristine that without nonbonded interactions) [24].](image-5.png "Figure 4 :") 5![Figure 5: Stress-strain curves of pristine graphene under (a) zigzag and (b) armchair direction;where SG (single graphene sheet), BG (LJ-On) (bilayer graphene sheet with non-bonded interactions) and BG (LJ-Off) (bilayer graphene without non-bonded interactions)[24].](image-6.png "Figure 5 :") 6![Figure6: Failure morphology of defective graphene with varying concentration of single vacancy defects in arm-chair direction (a1) initiation of failure in isolated single defective sheet of graphene at 0.5% single vacancy defects (a2) initiation of failure in isolated single defective sheet of graphene at 1.0% single vacancy defects (a3) initiation of failure in isolated single defective sheet of graphene at 1.5% single vacancy defects (b1) initiation of failure in defective sheet of graphene in bilayer sheets connected with non-bonded interactions at 0.5% single vacancy defects (b2) initiation of failure in defective sheet of graphene in bilayer sheets connected with non-bonded interactions at 1.0% defects (b3) initiation of in defective sheet of graphene in bilayer sheets connected with non-bonded interactions at 1.5% single vacancy defects (c1) initiation of failure in defective sheet of graphene in bilayer sheets without non-bonded interactions at 0.5% single vacancy defects (c2) initiation of failure in defective sheet of graphene in bilayer sheets without non-bonded interactions at 1.0% single vacancy defects (c3) initiation of failure in defective sheet of graphene in bilayer sheets without non-bonded interactions at 1.5% single vacancy](image-7.png "Figure 6 :") 7![Figure 7: Steps of breakage development in graphene having evenly concentrated vacancies; the vacancy defect ratio was (a) 0.5%, (b) 1%, (c) 1.5%, where, UDV refers to uniformly distributed vacancy defect ratios.](image-8.png "Figure 7 :") 8![Figure 8: (a) Breakage strong point of unperfected graphene sheet on stress and strain against the number of vacancy defect ratio in Zigzag direction where both lines remain the outcomes of quantized fracture mechanics (QFM). (b) Breakage strong point of an unperfected graphene sheet on stress and strain against the number of vacancy defect ratio in Armchair direction where both lines remain the outcomes of quantized fracture mechanics (QFM).](image-9.png "Figure 8 :") 2Studied byConditions/ Types of DefectsMethods AdoptedYoung's Modulus (TPa)Poisson's RatioJiang et al. 14T = 100-500 KMolecular Dynamics0.95 -1.10.17Shen et al . 15T = 300-700 KMolecular Dynamics0.905Lee et al. 16Pristine grapheneExperiment1 ± 0.1Tsai et al .17NPT ensembleMolecular Dynamics0.9120.261Sakhaee-Pour 18Pristine grapheneFinite Element Method1.025Georgantzinos et al. 19Pristine grapheneFinite Element Method1.367Kvashnin et al. 20Vacancy defectsMolecular Mechanics1.08Neek-Amal et al. 21randomly distributed vacancy defectsSTW defects0.501 ± 0.032Shokrieh et al. 22Pristine grapheneContinuum Mechanics1.04R.Ansary et al. 23STW defectsMolecular Dynamics60% reductionMuse Degefe & Avinash Parashar et al. 24Vacancy bi-layer T=300K grapheneMolecular Dynamics0.91II. © 2019 Global Journals * Effect of Point and Line Defects on Mechanical and Thermal Properties of Graphene: A Review GRajasekeran AvinashPrarthana Narayanan Parashar Critical Reviews in Solid State and Materials Sciences 2015 * Rotational actuators based on carbon nanotubes AMFennimore TDYuzvinsky Wei-QiangHan MSFuhrer JCumings AZettl Nature 424 6947 2003 * Morphology and in-plane thermal conductivity of hybrid graphene sheets BoLiu CDReddy JinwuJiang JuliaABaimova VSergey AyratADmitriev KunNazarov Zhou Applied Physics Letters 101 21 211909 2012 * Vibrational properties of graphene nanoribbons by first-principles calculations RolandGillen MarcelMohr ChristianThomsen Janinamaultzsch Physical Review B 80 15 155418 2009 * Phonon dispersion of graphite by inelastic x-ray scattering MMohr JMaultzsch EDobard?i ? SReich IMilo?evi? MDamnjanovi? ABosak MKrisch CThomsen Physical Review B 76 3 35439 2007 * SVDmitrriev JABaimoya AVSavin YSKivshar Comput. Mater. Sci 53 2012 * Scaling properties of flexible membranes from atomistic simulations: application to graphene JHLos IMikhail OVKatsnelson KVYazyev AnnalisaZakharchenko Fasolino Physical Review B 80 12 121405 2009 * Two-dimensional atomic crystals KSNovoselov DJiang FSchedin TJBooth VVKhotkevich SVMorozov AKGeim Proceedings of the National Academy of Sciences of the United States of America 102 30 2005 * Multiscale model to investigate the effect of graphene on the fracture characteristics of graphene/polymer nanocomposites AvinashParashar PierreMertiny Nanoscale research letters 7 1 2012 * Lee structures of thermally and chemically reduced graphenes HMTu SHHuh SHChoi HL Mater-lett 64 2010 * Graphene-DNA hybrids: self-assembly and electrochemical detection performance WeiLv MinGuo Ming-HuiLiang Feng-MinJin LanCui LinjieZhi Quan-HongYang J. Mater. Chem 20 32 2010 * Charge Transport in Polycrystalline Graphene: Challenges and Opportunities AWCummings DLDuong VLNguyen DVTuan JKotakoski JE BVarga Adv. Mater 26 * Characterizing mechanical properties of graphite using molecular dynamics simulation JLTsai JFTu Mater. Des 31 2010 * Anisotropic mechanical properties of graphene sheets from molecular dynamics ZNi HBu MZou HYi KBi YChen Physica B 405 2010 * Representative volume element to estimate buckling behavior of graphene/polymer nanocomposites PParashar Mertiny Nanoscale Res. Lett 7 515 2012 * Two-dimensional atomic crystals KSNovoselov DJiang FSchedin TJBooth VVKhotkevich SVMorozov AKGeim Proceedings of the National Academy of Sciences of the United States of America 102 30 2005 * Multiscale model to investigate the effect of graphene on the fracture characteristics of graphene/polymer nanocomposites AvinashParashar PierreMertiny Nanoscale research letters 7 1 2012 * Lee structures of thermally and chemically reduced graphenes HMTu SHHuh SHChoi HL Mater-lett 64 2010 * Observation of electron-antineutrino disappearance at Daya Bay FPAn JZBai ABBalantekin HRBand DBeavis WBeriguete MBishai Physical Review Letters 108 17 171803 2012 * Measurement of the elastic properties and intrinsic strength of monolayer graphene CLee XWei JWKysar JHone Science 321 2008 * Mechanical properties of suspended graphene sheets WFrank DMTanenbaum AMVan Der Zande PLMc Euen J. Vac. Sci. Technol. B 25 2007 * Young's modulus of Graphene: a molecular dynamics study WJiang JSWang BLi Phys. Rev. B Condens. Matter 80 113405 2009 * Temperaturedependent elastic properties of single layer graphene sheets LShen HSShen CLZhang Mater. Des 31 2010 * Effect of nonbonded interactions on failure morphology of a defective graphene sheet MDegefe AParashar Materials Research Express 3 4 45009 2016 * Measurement of the elastic properties and intrinsic strength of monolayer graphene CLee XWei JWKysar JHone Science 321 2008 * Elastic properties of single-layered graphene sheet Sakhaee-Pour Solid State Commun 149 2009 * Linear reduction of stiff-ness and vibration frequencies in defected circular monolayer graphene M FMPeters Phys. Rev. B Condens. Matter 81 235437 2010 * Prediction of Young's modulus of graphene sheets and carbon nanotubes using nanoscale continuum mechanics approach MShokrieh RRafiee Mater. Des 31 2010 * Mechanical properties of defective single-layered graphene sheets via molecular dynamics simulation, Super lattices Microstructure RAnsari SAjori BMotevalli 2012 51