The Effect of Hole Geometry on the Nonlinear Nanomechanics of γ-Graphyne Structures: A Finite Element Analysis

Graphyne is a material that has unique mechanical properties, but little is known about how these properties change when the material has holes. In this work, the effect of hole geometry, considering circular, triangle, and rhombus hole configurations, on the mechanical nonlinear response of γ-graphyne structures is studied. Graphyne, graphdiyne, graphyne-3, and graphyne-4 structures are under investigation. An efficient nonlinear finite element analysis (FEA) method is adequately implemented under large deformations for this purpose. The study varied the size and shape of the holes to understand how these changes affect the nanostructure’s mechanical response. The results indicate that the hole geometry significantly impacts the mechanical nonlinear response of γ-graphyne structures. The holes’ size and shape affect the structures’ elastic behavior, deformation, and strength. The findings can be used to optimize the design of γ-graphyne structures for specific mechanical applications. The study highlights the importance of considering the hole geometries in the design and fabrication of these materials.


Introduction
The introduction of monolayer graphene in 2004 was a scientific breakthrough in the field of materials due to its combination of superior mechanical, electromagnetic, optical, and thermal properties [1,2].In recent years, efforts have been made to discover materials and structures at the nanoscale that are akin to graphene or possess superior physical and structural attributes compared to graphene.Among the most promising recently identified graphene-like nanomaterials is the family of graphyne materials, initially examined theoretically by Balaban and colleagues [3], and subsequently by Narita and Nagai through first-principles analysis [4].Graphyne represents a novel category of twodimensional (2D) substances that have recently garnered interest owing to their distinct physical and chemical characteristics.Comprising sp-and sp 2 -bonded carbon atoms, these substances exhibit a honeycomb configuration akin to graphene, but with extra atoms or functional groups affixed to the carbon atoms.The presence of these additional atoms or groups gives rise to a wide range of properties, including enhanced electrical conductivity, mechanical strength, and chemical stability [5][6][7].
Despite the successful synthesis of graphdiyne by Li et al. in 2010 [8], there are limited experimental data available on the mechanical properties of graphyne structures in the literature.According to Ivanovskii [9], various types of graphyne and graphdiyne were thoroughly investigated for basic structural and mechanical properties, which were determined by experiments or evaluated by ab initio theory.Additionally, as reported by Zhang et al. [10], the Raman spectra of graphyne and graphdiyne were examined in a systematic manner, and the behavior of these materials under mechanical stress and strain was also investigated.A study by Xiao et al. [11] used atomic force microscopy experiments to determine the thermomechanical properties of graphdiyne films.
The effective design of advanced applications using nanostructures requires a thorough understanding of their properties and an ability to predict them using reliable theoretical or computational tools.Numerous molecular dynamics (MD) simulations and ab initio studies have been carried out on structures like graphene, carbon nanotubes, and fullerenes to comprehend a broad spectrum of properties [12][13][14] and behaviors [15][16][17].Additionally, the finite element method has been used to provide a more comprehensive approach to the computational modeling of nanostructures.With regard to nanostructures, several forms of finite elements, such as beams [18][19][20], spring-like elements [21][22][23], shells [24,25], and bar elements [26,27], have been formulated to analyze the mechanical performance of these materials.
Atomistic finite element models simulate interatomic interactions of atoms using appropriate finite element formulations.Couto and Silvestre [28] introduced an FEA modeling approach to estimate the elastic behavior of graphyne.They conducted uniaxial and biaxial tensile simulations using computational models of zigzag and armchair configurations and compared their results with equivalent molecular dynamics and density functional theory models, finding that they were in good agreement.Lee et al. [29] applied an atomistic FE model to analyze the sensitivity of various graphyne-based resonators under different boundary conditions and size configurations.They compared the results with those of graphene-based resonators and concluded that graphyne derives higher sensitivity.Georgantzinos et al. [30,31] investigated the elastic and elastoplastic mechanical properties of γ-graphyne family structures utilizing FE analysis.They conducted various parametric studies, concentrating on the size-dependency of the sheet structure on the mechanical properties, and generalized their results using regression models.Galhofo and Silvestre [32] derived stress-strain curves to examine both the monotonic and hysteretic behavior of γ-graphyne sheets using the same geometric configurations as in [28].Siampanis et al. [33] examined the stress-strain behavior of graphyne structures under shear loading and obtained higher shear strength and stiffness compared to other structures of this family.
The available studies on the mechanical properties of graphyne structures by numerical methods assume ideal sheet structures without any defects or holes.However, in the case of graphene, a particular interest in the modeling of the structure lies in the assessment of the influence of nano-holes on the mechanical behavior.This is in line with the investigation of nanoporous graphene or graphene sheet with holes that have been studied due to its good properties for gas absorption [34], water desalination [35], energy storage [36,37], and nanoelectronics applications [38].Genoese et al. [39] introduced a new interatomic potential, the Dreading potential, in which they added damping functions, to investigate the in-plane and out-plane tensile behavior of monolayer graphene sheets.In their study, they examined the effect of the diameter of a hole located in the middle of the sheet on the stress and strain distributions.Another study on the effect of holes on graphene sheets can be found in [40], where Muraru et al. developed a new software tool to generate graphene-based molecular models.Yoon et al. [41] examined the mechanical resonance properties of porous graphene sheets computationally and experimentally.However, graphyne structures fall behind in terms of published papers regarding the role of holes and their effect on mechanical properties.Recently, Lee et al. [42] used an atomistic FE model to examine the effect of holes on Young's moduli of porous γ-graphyne sheets under tensile loading, concluding that an increase in porosity decreases Young's modulus.
The presence of holes in a graphyne structure can play a practical role in several areas such as electronics, energy storage, and catalysis.In electronics, the holes can be used to create electronic devices with high electron mobility, such as field-effect transistors.In energy storage, graphyne with holes can be used as a high-capacity electrode material in lithium-ion batteries.The holes can also be used as active sites for catalytic reactions in chemical synthesis, such as hydrogenation reactions.The size and shape of the holes can also be used to control the behavior of the graphyne, such as its strength, electronic conductivity, and chemical reactivity.Additionally, by controlling the size and shape of the holes, it may be possible to create new types of graphyne structures with unique properties that can be used in a variety of applications.As a recent addition to the family of carbon allotropes, graphyne structures have not been extensively studied in comparison to other carbon nanostructures.In this study, we investigate the effect of hole geometry on the mechanical nonlinear response of γ-graphyne structures, using an efficient nonlinear finite element method (FEM).Specifically, we examine the mechanical properties of graphyne, graphdiyne, graphyne-3, and graphyne-4 structures with circular, triangle, and rhombus hole configurations.The implementation of the FEA model is conducted under large deformations, and it is an atomistic spring-based finite element model in which a unique spring is used to simulate the interatomic interactions between carbon atoms.

Computational Model 2.1. Geometric Characteristics
Graphyne consists of three types of bonds: aromatic, single, and triple.The repeatable unit cell of the structure is the single-triple-single bond interconnection.The number of the unit-cell triple bonds defines the graphyne variation, i.e., one bond yields graphyne, two bonds yield graphdiyne, three bonds yield graphyne-3, and four bonds yield graphyne-4. Figure 1 shows the structural configuration of a graphyne sheet along with the bond variations.The effects of the geometric parameters of the holes are indirectly calculated via the calculation of the area of the hole.The dimensions of all γgraphyne models were approximately 10 × 10 nm.The investigations are carried out for five different aspect ratios, including 5%, 10%, 15%, 20% and 25%.The hole aspect ratio AR is obtained by dividing the area of the hole (Af ) by the total area of each nanosheet (At).
This approach allows for a more generalized analysis of the impact of hole geometry on the mechanical properties of the graphyne structures, as it accounts for variations in hole size and shape by normalizing the hole area to the total area of the nanosheet.
Building on the aforementioned methodology, the determination of the effect of the parameter becomes straightforward when utilizing the provided formulas for specific hole shapes and given values of AR and At.For a circular hole, the diameter (D) can be derived from the equation D = √ 4ARAt/π.Similarly, for a square hole, the side length (A) is given by A = √ ARAt.In the case of a rhombus-shaped hole, where both diagonals are equal (B = C), the diagonal length (B) can be calculated using the formula B = √ 2ARAt.These formulas offer a systematic approach to determine the dimensions of the holes based on their shape, ensuring a consistent and accurate analysis.By employing these equations, researchers can effectively gauge the influence of different hole geometries on the mechanical properties of γ-graphyne structures, further enhancing the comprehensiveness of the study.

Force Field Description
The aggregate potential energy of a graphyne molecular arrangement can be depicted as the cumulative energies from interatomic interactions.Considering that nonbonded interactions are negligible for the molecular systems under investigation, and only inplane attributes are being assessed, the potential energy of a γ-graphyne configuration is articulated as the cumulative energies of bond stretching and bond angle bending: where U i r represents the bond stretching term and U i θ denotes the bond angle bending term, respectively.The parameter i denotes the kind of bond, assuming the values s, a, and t correspond to the single, aromatic, and triple bond, respectively.According to the Morse potential field [31], the potential energy terms can be written as In Equation ( 2), r i and r i 0 represent the changed and original bond lengths, respectively.D i e and B i are constants that vary based on the type of bond (s, a, t).In Equation (3), θ i,j and θ i,j 0 are the angles between the deformed and original state of two adjacent bonds i and j, respectively.k i,j θ and k i,j sextic are constants that control the linear and nonlinear force-deformation relationship of bond angle bending interactions.It is important to note that nonbonded interatomic interactions and resulting nonlocality typically have minimal impact on the overall mechanical behavior of graphene-like nanostructures and can be omitted in the mathematical formulation.
These formulas can be differentiated to mimic the force field inside the nanostructure using a mix of translational and rotational springs.The longitudinal rigidity of the translational springs is determined by differentiating Equation (2) concerning the variation in bond length (∆r i = r i − r i 0 ), and the necessary rigidity of the rotational spring can be obtained by differentiating Equation (3) concerning the variation in the bending angle (∆θ i,j = θ i,j − θ i,j 0 ), where the developed axial force between two bonded atoms F i r ∆r i is a result of a change in their interatomic distance, while a bending moment M i,j θ ∆θ i,j arises from a change in the angle between two linked bonds.

Finite Element Model
The simulation of potential energies resulting from bond length and bending angle variations is consistent with the finite element formulation outlined in [30,32].This is achieved by utilizing two nodes, spring-based line elements that connect two bonded atoms.The longitudinal rigidity of the translational springs is determined by differentiating Equation (4) concerning the variation in bond length, while the rotational spring is derived by differentiating Equation ( 5) with respect to bending angle variation.
The stiffness matrix of the finite elements and their force-displacement behavior are deduced by using a local Cartesian coordinate system (x, y).This results in a 2 × 2 elemental matrix that incorporates two fundamental coefficients in the x and y directions: where the symbolization el of the element receives distinct string values: sas, ast, sts, and tst, indicating that the current formulation necessitates the implementation of four kinds of spring-like elements with varying longitudinal and transverse stiffness coefficients, as illustrated in [30,32].For example, the notation 'sts' corresponds to a spring-based element used to describe a triple bond located between two single bonds.The nonlinear axial stiffness coefficients of these elements can be determined using the subsequent equation.
The notation el of the element is assigned one of four different string values: sas, ast, sts, and tst.This indicates that the current formulation necessitates the implementation of four types of spring-like elements, each with distinct longitudinal and transverse stiffness coefficients, as outlined in [21,31].For example, the notation sts represents a spring-like element used to describe a triple bond between two single bonds.The nonlinear axial stiffness coefficients of these elements can be determined using the following equation.
In accordance with the simplified modeling technique, detailed in [42], the bond angle bending variation can be effectively addressed by utilizing an suitable rigidity coefficient in the y-direction for each element el = sas,ast,sts,tst.This is provided by: To investigate the elastoplastic behavior of a nanostructure, the equilibrium equation for each nonlinear spring-based element can be written as: The equilibrium equation is for each hypothetical nonlinear spring-based element, where u el and f el represent the respective force and displacement vectors for each element.These vectors have the following notation if i and j are the two nodes of the finite element: The system of nonlinear equations can be assembled into its final form, through the transformation of the global coordinate system to the elemental stiffness equation for each finite element, in accordance with the constraints of nodal connectivity: where F el , U el , and K el (U) el are the assembled force vector, assembled displacement vector, and assembled deformation-dependent stiffness matrix, respectively.
To simulate the nonlinear mechanical tensile tests, appropriate supports and loads are applied as boundary conditions.The next step is to use an established incremental-iterative method based on the Newton-Raphson algorithm to address the nonlinearity of the global stiffness matrix numerically and determine the response of the structure as a result.

Results and Discussion
The size of the structure and the number of aromatic rings are indeed inversely related when considering the length of the acetylenic linkage in graphyne structures.Specifically, graphyne, with the shortest acetylenic linkage, accommodates more aromatic rings within a given size compared to its counterparts.As we progress to graphdiyne, graphyne-3, and graphyne-4, the length of the acetylenic linkage increases, which results in a reduced number of aromatic rings for structures of the same size (Figure 1).This relationship has significant implications for the mechanical properties of these materials.Aromatic rings inherently contribute to the rigidity and stability of the structure due to their conjugated and electron-delocalized nature.Therefore, materials with a higher number of aromatic rings, like graphyne, tend to exhibit superior mechanical properties [30,31].In essence, the more aromatic rings present in the structure, the better the mechanical strength and stability.
The proposed finite element model has been previously validated through comparisons with other available results in the open literature [31] regarding pristine graphyne structures.However, a number of comparisons were performed between the outcomes of this computational approach and related data taken from various literature sources, as provided in Table 1, in order to validate the numerical calculation of the key mechanical properties of γ-graphynes.This comparison serves as further verification of the validity and accuracy of the proposed model.These sources employed different approaches to extract the data.The results indicate a reasonable agreement with the findings of previous studies and demonstrate that the model can be used to predict the mechanical properties of graphyne structures with holes with a high degree of accuracy, providing confidence in the accuracy of the computational scheme, and thus can be useful for the design and optimization of graphyne-based materials and devices.The behavior of the elastic modulus of graphyne structures (graphyne, graphdiyne, graphyne-3, graphyne-4) with increasing aspect ratio of the hole is an important aspect of understanding the mechanical properties of these materials.Figure 3 provides insights on this behavior, illustrating the change in elastic modulus for different types of holes.Figure 3a,b show the behavior of the elastic modulus in the x and y axis for the case of a central circular hole for the four structures.Figure 3b,c depict the case of a square hole, while Figure 3c,d describe the case of a rhombus-type hole.The elastic modulus is seen to decrease as the size of the hole increases in all circumstances.Additionally, the rate of decrease in the elastic modulus with an increase in hole size is similar for all cases.This suggests that the graphyne structures have a relatively consistent response to the introduction of holes in their structure.However, small fluctuations in the elastic modulus may be attributed to the non-symmetric geometric effects caused by the hole creation, depending on the number of atoms and bonds cut off.As the length of the acetylene series increases, the elastic modulus is generally lower in all cases.This highlights the importance of the number of atoms and bonds in determining the mechanical properties of graphyne structures.Close values are observed concerning the elastic modulus in the y-direction between the graphyne and graphdiyne specifically, indicating that these two structures may have similar mechanical properties in certain directions.
The behavior of the elastic modulus of graphyne structures (graphyne, graphdiyne, graphyne-3, graphyne-4) with increasing aspect ratio of the hole is an important aspect of understanding the mechanical properties of these materials.Figure 3 provides insights on this behavior, illustrating the change in elastic modulus for different types of holes.Figure 3a,b show the behavior of the elastic modulus in the x and y axis for the case of a central circular hole for the four structures.Figure 3b,c depict the case of a square hole, while Figure 3c,d describe the case of a rhombus-type hole.The elastic modulus is seen to decrease as the size of the hole increases in all circumstances.Additionally, the rate of decrease in the elastic modulus with an increase in hole size is similar for all cases.This suggests that the graphyne structures have a relatively consistent response to the introduction of holes in their structure.However, small fluctuations in the elastic modulus may be attributed to the non-symmetric geometric effects caused by the hole creation, depending on the number of atoms and bonds cut off.As the length of the acetylene series increases, the elastic modulus is generally lower in all cases.This highlights the importance of the number of atoms and bonds in determining the mechanical properties of graphyne structures.Close values are observed concerning the elastic modulus in the y-direction between the graphyne and graphdiyne specifically, indicating that these two structures may have similar mechanical properties in certain directions.The mechanical tensile strength of graphyne structures (graphyne, graphdiyne, graphyne-3, graphyne-4) seems to be considerably affected by the aspect ratio of the hole.Figure 4 illustrates the change in ultimate strength for different types of holes.Figure 4a,b show the ultimate strength in the x and y axis for the case of a central circular hole for the four structures.Figure 4b,c depict the case of a square hole, while Figure 4c,d describe the case of a rhombus-type hole.As the size of the hole increases, the structure can experience a significant reduction in strength, with losses of up to almost 50% observed for aspect Figure 3.The elastic modulus of graphyne structures vs. aspect ratio in (a) x-direction for the circular hole, (b) y-direction for the circular hole, (c) x-direction for the square hole, (d) y-direction for the square hole, (e) x-direction for the rhombus-type hole, and (f) y-direction for the rhombus-type hole.
The mechanical tensile strength of graphyne structures (graphyne, graphdiyne, graphyne-3, graphyne-4) seems to be considerably affected by the aspect ratio of the hole.Figure 4 illustrates the change in ultimate strength for different types of holes.Figure 4a,b show the ultimate strength in the x and y axis for the case of a central circular hole for the four structures.Figure 4b,c depict the case of a square hole, while Figure 4c,d describe the case of a rhombus-type hole.As the size of the hole increases, the structure can experience a significant reduction in strength, with losses of up to almost 50% observed for aspect ratios near 0.25.The rate of decrease in the ultimate strength with an increase in hole size is similar for all cases.Small fluctuations in the ultimate strength may be observed for the same reasons as previously.The mechanical tensile strength of graphyne structures (graphyne, graphdiyne, graphyne-3, graphyne-4) seems to be considerably affected by the aspect ratio of the hole.Figure 4 illustrates the change in ultimate strength for different types of holes.Figure 4a,b show the ultimate strength in the x and y axis for the case of a central circular hole for the four structures.Figure 4b,c depict the case of a square hole, while Figure 4c,d describe the case of a rhombus-type hole.As the size of the hole increases, the structure can experience a significant reduction in strength, with losses of up to almost 50% observed for aspect ratios near 0.25.The rate of decrease in the ultimate strength with an increase in hole size is similar for all cases.Small fluctuations in the ultimate strength may be observed for the same reasons as previously.As the length of the acetylene series increases, the ultimate strength is generally lower in all cases.This emphasizes the significance of the number of atoms and bonds in establishing the mechanical properties of graphyne structures.
In Figure 5, we observe the variation in fracture strain for different types of holes.Figure 5a,b present the fracture strain in the x and y axis for the case of a central circular hole for the four structures.Figure 5b,c depict the case of a square-type hole, while Figure 5c,d describe the case of a rhombus-type hole.It is noticeable that in all cases, the rate of As the length of the acetylene series increases, the ultimate strength is generally lower in all cases.This emphasizes the significance of the number of atoms and bonds in establishing the mechanical properties of graphyne structures.
In Figure 5, we observe the variation in fracture strain for different types of holes.Figure 5a,b present the fracture strain in the x and y axis for the case of a central circular hole for the four structures.Figure 5b,c depict the case of a square-type hole, while Figure 5c,d describe the case of a rhombus-type hole.It is noticeable that in all cases, the rate of decrease in the fracture strain with an increase in hole size is similar for all cases.Slight variations in the fracture strain can be attributed to non-symmetric geometric effects caused by the hole creation, depending on the number of atoms and bonds cut off.Similarly, here, as the length of the acetylene series increases, the fracture strain is generally lower in all cases.We conducted a nonlinear regression analysis to extract analytical equations that capture the hole-dependent variations in mechanical properties of the four types of graphyne Figure 5.The fracture strain of graphyne structures vs. aspect ratio in (a) x-direction for the circular hole, (b) y-direction for the circular hole, (c) x-direction for the square hole, (d) y-direction for the square hole, (e) x-direction for the rhombus-type hole, and (f) y-direction for the rhombus-type hole.
We conducted a nonlinear regression analysis to extract analytical equations that capture the hole-dependent variations in mechanical properties of the four types of graphyne with central holes, aiming to enable a more systematic and simplified prediction of their mechanical properties.The following fitting 3-parameter function of hole aspect ratio was used for all of the attempted property approximations: where Y is the suggested mathematical model to fit the data computed for all examined properties E x , E y , σu x , σu y , εf x , and εf y whereas, a, b, and c are constant parameters that must be found using a different regression procedure for each mechanical property variation.Y 0 is the corresponding property of pristine graphyne structure.Equation (20) from [31] is a valuable tool for predicting the mechanical properties of graphyne structures, by considering the effect of size variations.This equation can be used to estimate Y o , an important parameter for understanding the mechanical behavior of these structures.By utilizing Equation (17), which relates Y o to the corresponding mechanical properties Y, predictions of any size can be made, providing a comprehensive understanding of the mechanical response of graphyne structures with holes.By incorporating the FEM results and performing a regression analysis using Equation (17), the optimized values of the parameters can be calculated.Appendix A presents the regression curves for each graphyne structure type, hole type, and mechanical property in Figures A1-A3.These figures also include the average curve obtained by considering a single equation for all graphyne structure types.Table 2 provides the parameter values for the unique equation applied to all graphyne structure types.For a more accurate prediction, Table A1 contains the values of the parameters for Equation ( 17) that were obtained from the fitting process for mechanical properties data of graphyne and each hole type.The fitting process for the mechanical properties data based on Equation ( 17) yielded the parameter values shown in Table A2, regarding the graphdiyne.Tables A3 and A4 provide the parameter values obtained from fitting Equation (17) to the mechanical properties data for each hole type in graphyne-3 and graphyne-4, respectively.The nonlinear FEA method was implemented under large deformations to examine the mechanical properties of γ-graphyne structures with circular, triangle, and rhombus hole configurations.The study used an atomistic spring-based finite element model, simulating the interatomic interactions between carbon atoms.This approach provided a comprehensive understanding of the mechanical response of graphyne structures with holes, enabling a more systematic and simplified prediction of their mechanical properties.
In addition to the FEA used in this study, the nonlinear response of γ-graphyne structures can also be investigated using several other robust numerical modeling methods.For instance, the finite difference method [46], Bezier multi-step method [47], and differential quadrature method [48] can provide alternative approaches to understanding the nonlinear response of γ-graphyne structures and can be considered for future studies or comparisons.
The study methodically investigates the nonlinear mechanical response of γ-graphyne structures, including variations such as graphyne, graphdiyne, graphyne-3, and graphyne-4, under the influence of different hole geometries (circular, square, and rhombus) and aspect ratios (5%, 10%, 15%, 20%, 25%).The investigation reveals that hole geometry significantly impacts the elastic behavior, deformation, and strength of the structures, with the elastic modulus, ultimate strength, and fracture strain decreasing as hole size increases.This is attributed to the variations in stress distribution and deformation patterns caused by the holes.The study also presents analytical equations derived from regression analysis to capture hole-dependent variations in mechanical properties, providing a valuable tool for systematic and simplified prediction of mechanical properties.This comprehensive analysis underscores the importance of considering hole geometries in the design and optimization of graphyne-based materials and devices, ultimately aiding in the development of materials with desired mechanical properties for specific applications.
The study of hole geometry, encompassing circular, triangular, and rhombus configurations, is pivotal for comprehending the nonlinear mechanical behavior of γ-graphyne structures.Different hole geometries lead to variations in the stress distribution, deformation patterns, and ultimately, the mechanical properties of the material.Specifically, the hole geometry influences the initiation and propagation of cracks, affecting the material's strength and failure mechanisms.This understanding is crucial for optimizing the design and fabrication of γ-graphyne structures, as it can lead to enhanced mechanical properties, making the material more suitable for various applications, such as in nanoelectronics and nanomechanical systems.

Conclusions
In this study, we investigated the effect of hole geometry on the mechanical nonlinear response of γ-graphyne structures, including graphyne, graphdiyne, graphyne-3, and graphyne-4 structures, under large deformations using a finite element analysis (FEA) method.We found that the size and shape of the holes significantly impacted the structures' elastic behavior, deformation, and strength, highlighting the importance of considering hole geometries in the design and fabrication of graphyne-based materials and devices.The results indicated that the size and shape of the holes significantly impacted the structures' elastic behavior, deformation, and strength.Specifically, the elastic modulus of the graphyne structures decreased with an increase in the size of the hole, and the ultimate strength and fracture strain experienced a similar rate of decrease with an increase in hole size.Small fluctuations in these properties were observed due to non-symmetric geometric effects caused by the hole creation.Our study provides a valuable tool for predicting the mechanical properties of graphyne structures with holes by incorporating the effect of size variations, enabling a more systematic and simplified prediction of their mechanical properties.The analytical equations extracted from the regression analysis, which capture the hole-dependent variations in mechanical properties, were presented.This study highlights the importance of considering the hole geometries in the design and fabrication of these materials, providing a comprehensive understanding of the mechanical response of graphyne structures with holes, which can be useful for the design and optimization of graphyne-based materials and devices.Overall, our findings can aid in the optimization of the design of graphyne structures for specific mechanical applications.

Appendix A
Figure A1 shows the fitting curves based on Equation ( 17) that fit the FEM data for each structure and hole type, as well as the corresponding average fitting curve, which represents the effect of hole magnitude on the Young's modulus of the graphyne structures.Figure A2 displays the fitting curves for the ultimate tensile strength, while Figure A3 presents the corresponding fitting curves for the tensile fracture strain.
. Fitting curves for the normalized fracture strain of graphyne structures vs aspect ratio in (a) x-direction for the circular hole, (b) y-direction for the circular hole, (c) x-direction for the square hole, (d) y-direction for the square hole, (e) x-direction for the rhombus-type hole, and (f) y-direction for the rhombus-type hole.
Tables A1-A3 provide the values of the parameters obtained from the fitting process for Equation (17) applied to the normalized Young's modulus, ultimate tensile strength, and tensile fracture strain data, respectively.
The presented tables, namely Tables A1-A4, provide crucial parameter values for Equation (17), which is an analytical model used to predict the mechanical properties of different types of graphyne materials with circular, square, and rhombic holes.Specifically, Tables A1-A4 contain the parameter values for graphyne, graphdiyne, graphyne-3, and graphyne-4, respectively.These tables provide valuable information for designing and engineering materials with specific mechanical properties by controlling the hole shape and the material type.The R 2 values in the tables indicate the goodness of fit of the analytical model and can be used to evaluate the reliability of the predicted mechanical properties.Overall, these tables offer useful insights into the mechanical properties of graphyne materials with different hole shapes, and their practical applications in material science and engineering can be significant.Tables A1-A3 provide the values of the parameters obtained from the fitting process for Equation (17) applied to the normalized Young's modulus, ultimate tensile strength, and tensile fracture strain data, respectively.
The presented tables, namely Tables A1-A4, provide crucial parameter values for Equation (17), which is an analytical model used to predict the mechanical properties of different types of graphyne materials with circular, square, and rhombic holes.Specifically, Tables A1-A4 contain the parameter values for graphyne, graphdiyne, graphyne-3, and graphyne-4, respectively.These tables provide valuable information for designing and engineering materials with specific mechanical properties by controlling the hole shape and the material type.The R 2 values in the tables indicate the goodness of fit of the analytical model and can be used to evaluate the reliability of the predicted mechanical properties.Overall, these tables offer useful insights into the mechanical properties of graphyne materials with different hole shapes, and their practical applications in material science and engineering can be significant.

Figure 2
Figure 2 depicts the geometric failures, which are the circular hole, the square-type hole and rhombus-type hole.The diameter of the circular hole is denoted as D, the length of the side of the square-type hole is A, and the lengths of the rhombus-type hole are B and C. In the case studied, the interior angles of the rhombus are 90 degrees.The holes are located at the geometrical center of the structure.The effects of the geometric parameters of the holes are indirectly calculated via the calculation of the area of the hole.The dimensions of all γgraphyne models were approximately 10 × 10 nm.The investigations are carried out for five different aspect ratios, including 5%, 10%, 15%, 20% and 25%.The hole aspect ratio AR is obtained by dividing the area of the hole (Af ) by the total area of each nanosheet (At).

Figure 3 .
Figure3.The elastic modulus of graphyne structures vs. aspect ratio in (a) x-direction for the circular hole, (b) y-direction for the circular hole, (c) x-direction for the square hole, (d) ydirection for the square hole, (e) x-direction for the rhombus-type hole, and (f) y-direction for the rhombus-type hole.

Figure 3 .
Figure3.The elastic modulus of graphyne structures vs. aspect ratio in (a) x-direction for the circular hole, (b) y-direction for the circular hole, (c) x-direction for the square hole, (d) ydirection for the square hole, (e) x-direction for the rhombus-type hole, and (f) y-direction for the rhombus-type hole.

Figure 4 .
Figure 4.The ultimate tensile strength of graphyne structures vs. aspect ratio in (a) x-direction for the circular hole, (b) y-direction for the circular hole, (c) x-direction for the square hole, (d) ydirection for the square hole, (e) x-direction for the rhombus-type hole, and (f) y-direction for the rhombus-type hole.

Figure 4 .
Figure 4.The ultimate tensile strength of graphyne structures vs. aspect ratio in (a) x-direction for the circular hole, (b) y-direction for the circular hole, (c) x-direction for the square hole, (d) y-direction for the square hole, (e) x-direction for the rhombus-type hole, and (f) y-direction for the rhombus-type hole.

Figure 5 .
Figure5.The fracture strain of graphyne structures vs. aspect ratio in (a) x-direction for the circular hole, (b) y-direction for the circular hole, (c) x-direction for the square hole, (d) y-direction for the square hole, (e) x-direction for the rhombus-type hole, and (f) y-direction for the rhombus-type hole.

Figure A1 .
Figure A1.Fitting curves for the normalized Young's modulus of graphyne structures vs. aspect ratio in (a) x-direction for the circular hole, (b) y-direction for the circular hole, (c) x-direction for the square hole, (d) y-direction for the square hole, (e) x-direction for the rhombus-type hole, and (f) y-direction for the rhombus-type hole.

Figure A1 .
Figure A1.Fitting curves for the normalized Young's modulus of graphyne structures vs. aspect ratio in (a) x-direction for the circular hole, (b) y-direction for the circular hole, (c) x-direction for the square hole, (d) y-direction for the square hole, (e) x-direction for the rhombus-type hole, and (f) y-direction for the rhombus-type hole.

Figure A2 .
Figure A2.Fitting curves for the normalized ultimate tensile strength of graphyne structures vs. aspect ratio in (a) x-direction for the circular hole, (b) y-direction for the circular hole, (c) xdirection for the square hole, (d) y-direction for the square hole, (e) x-direction for the rhombus-type hole, and (f) y-direction for the rhombus-type hole.

Figure A3 .
Figure A3.Fitting curves for the normalized fracture strain of graphyne structures vs aspect ratio in (a) x-direction for the circular hole, (b) y-direction for the circular hole, (c) x-direction for the square hole, (d) y-direction for the square hole, (e) x-direction for the rhombus-type hole, and (f) y-direction for the rhombus-type hole.

Table 1 .
Comparison of the mechanical characteristics of pristine graphene predicted by the finite element method with equivalent findings from previous investigations.

Table 2 .
(17)meter values of Equation(17)for the analytical prediction of mechanical properties.

y σu x σu y εf x εf y
Fitting curves for the normalized ultimate tensile strength of graphyne structures vs. aspect ratio in (a) x-direction for the circular hole, (b) y-direction for the circular hole, (c) x-direction for the square hole, (d) y-direction for the square hole, (e) x-direction for the rhombus-type hole, and (f) y-direction for the rhombus-type hole.Fitting curves for the normalized ultimate tensile strength of graphyne structures vs. aspect ratio in (a) x-direction for the circular hole, (b) y-direction for the circular hole, (c) xdirection for the square hole, (d) y-direction for the square hole, (e) x-direction for the rhombus-type hole, and (f) y-direction for the rhombus-type hole.

Table A1 .
Parameter values of Equation (17) for the analytical prediction of mechanical properties of graphyne.

Table A1 .
Parameter values of Equation (17) for the analytical prediction of mechanical properties of graphyne.

Table A2 .
Parameter values of Equation (17) for the analytical prediction of mechanical properties of graphdiyne.

Table A3 .
Parameter values of Equation (17) for the analytical prediction of mechanical properties of graphyne-3.

Table A4 .
Parameter values of Equation (17) for the analytical prediction of mechanical properties of graphyne-4.