Next Article in Journal
Application of Hollow-Core Photonic Crystal Fibers in Gas Raman Lasers Operating at 1.7 μm
Next Article in Special Issue
Ultrahigh Ballistic Resistance of Twisted Bilayer Graphene
Previous Article in Journal
Synthesis, Structure and In Vitro Anticancer Activity of Pd(II) Complex of Pyrazolyl-s-Triazine Ligand; A New Example of Metal-Mediated Hydrolysis of s-Triazine Pincer Ligand
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance of SCAN Meta-GGA Functionals on Nonlinear Mechanics of Graphene-Like g-SiC

1
Physics Department, King Fahd University of Petroleum & Minerals, Dhahran 31261, Saudi Arabia
2
K.A.CARE Energy Research & Innovation Center at Dhahran, Dhahran 31261, Saudi Arabia
Crystals 2021, 11(2), 120; https://doi.org/10.3390/cryst11020120
Submission received: 5 January 2021 / Revised: 20 January 2021 / Accepted: 23 January 2021 / Published: 27 January 2021
(This article belongs to the Special Issue Graphene Mechanics (Volume II))

Abstract

:
Although meta-generalized-gradient approximations (meta-GGAs) are believed potentially the most accurate among the efficient first-principles calculations, the performance has not been accessed on the nonlinear mechanical properties of two-dimensional nanomaterials. Graphene, like two-dimensional silicon carbide g-SiC, has a wide direct band-gap with applications in high-power electronics and solar energy. Taken g-SiC as a paradigm, we have investigated the performance of meta-GGA functionals on the nonlinear mechanical properties under large strains, both compressive and tensile, along three deformation modes using Strongly Constrained and Appropriately Normed Semilocal Density Functional (SCAN) as an example. A close comparison suggests that the nonlinear mechanics predicted from SCAN are very similar to that of Perdew-Burke-Ernzerhof (PBE) formulated functional, a standard Density Functional Theory (DFT) functional. The improvement from SCAN calculation over PBE calculation is minor, despite the considerable increase of computing demand. This study could be helpful in selection of density functionals in simulations and modeling of mechanics of materials.

1. Introduction

Both elements of carbon and silicon can form mono-elemental 2D crystal structures. Carbon atoms of a graphene mono-layer is coplanar, whereas silicene has a buckled structure. The structure of silicon carbide mono-layer is coplanar [1]. Graphene-like SiC mono-layer (i.e., g-SiC ) is a stable 2D honeycomb structure with a wide indirect band gap (∼2.55 eV) [2]. It is the most thermodynamically stable one among over 300 SixCy structures [1]. g-SiC has promising applications as blue and ultraviolet light emitted devices (LEDs), and photo-voltaic devices [2,3,4]. It also has potential applications in sensors, composites, field emitters, super-capacitors, nanoelectromechanical devices, bio-imaging probes, catalysts, field-effect transistors, and wave absorber. g-SiC mono-layers have not been fabricated yet: only a few layers 2D-SiC films have so far been fabricated through solution-based exfoliation [5]. Extensive theoretical studies have been done on the structural stability, thermodynamic stability, electronic, catalytic, and optical properties of g-SiC [3,6,7,8,9,10,11,12,13,14,15,16,17,18].
Strain engineering is a common approach to tune a material’s properties, especially two-dimensional nanomaterials due to their outstanding flexibility [19]. In general, a nanomaterial is vulnerable to strain [20,21,22,23,24,25], which could be induced by thermal vibration, surface adhesion, epitaxial grown, substrate deformation, thermal expansion mismatch, pre-stretched substrate, pressurized blisters, substrate topography modification, tip indentation, and many other process. It was predicted that the indirect-direct band gap transition of the SiC nanosheet can be tuned through changing the lattice constant [3] which is equivalent to biaxial straining. Control of the electronic and optical properties of multi-layer g-SiC sheets by strain engineering up to 9% strain was also reported recently [26]. The knowledge the mechanical properties is crucial to the applications of these nanomaterials [27,28,29].
We have investigated the non-linear mechanical behaviors by means of first-principles calculations within the framework of Density Functional Theory (DFT) [30]. DFT is a feasible first-principles method with less computationally demanding than other computational methods with a similar accuracy. Therefore, DFT is the most widely-used method of electronic structure calculations. The biggest challenge of DFT is the description of the exchange-correlation term, which represents the correction to the kinetic energy and the electron-electron repulsion energy [31]. There are many approximations to the exchange-correlation energy, with different level of accuracy, which are described by various rungs of “Jacob’s Ladder” [32]. The first three rungs are the local density approximation (LDA), generalized gradient approximation (GGA), and meta-GGA approximation (MGGA) [33]. Among them, Perdew-Burke-Ernzerhof (PBE) formulated GGA [34] is widely used as a “standard” of DFT calculations for its combined good performance of accuracy and computation demand. Recently, Strongly Constrained and Appropriately Normed Semilocal Density Functional (SCAN) rises up with promising improvement of accuracy whereas not-much increasing of computing demands [35,36]. There is a concern about the precision of first-principles modeling of the nonlinear elastics at large strain. In general, it is believed that the selection of exchange-correlation functionals affect the prediction precisions in DFT calculations. To which extent, however, is still obscure.
In this work, we aim to assess the performance of the rungs of DFT exchange-correlation approximations on the outcome precision of nonlinear mechanics of 2D materials. We focus on the GGA (second rung) and the MGGA (the third rung) for g-SiC, taking PBE-GGA and SCAN-MGGA as an example. We deform the g-SiC mono-layers to various large strains and compute the high order (up to 5th order) elastic constants. These elastic constants are important for an accurate nonlinear elastic description. Furthermore, the pressure dependence properties are predicted. These properties include sound wave speed (both p-wave and s-wave) and the second order elastic constants. Our method is introduced in the Section 2. Our results and discussion are presented in Section 3, followed by Conclusions in Section 4.

2. Computation Method

We adopted the model of a 6-atom conventional unitcell, which has 3 carbon atoms and 3 silicon atoms, as shown in Figure 1. This large unitcell is relatively larger than the primitive unitcell which only contains two atoms (one carbon and one silicon). Our model is capable to capture the soft mode in hexagonal rings [37]. We performed the DFT calculations using the program of Vienna Ab-initio Simulation Package (VASP) [38]. We adopted the PBE-GGA [34] and SCAN-MGGA [35] for exchange-correlation functions in DFT calculations in this study. We denoted them as “PBE” and “SCAN”, respectively, in the rest of text for convenience of description. We have employed projector augmented wave (PAW) and pseudo-potential approach [39] for modeling electrons. For carbon atoms, the 2 s 2 2 p 2 electrons are explicitly considered. For silicon atoms, the 3 s 2 3 p 2 electrons are included. The van der Waals interactions are treated with DFT+q method [40]. The periodic boundary conditions are applied along all three directions to avoid the surface effect. A cutoff of energy of 600 eV for the basis of plane waves is used in all the DFT calculations. The criterion of electronic relaxation is the convergence of total energy within a fluctuation less than 1.0 × 10 5 eV. The ionic relaxation is achieved when the atomic forces are less than 0.001 eV/Å. The irreducible Brillouin Zone of the unitcell is sampled with a Gamma-centered 24 × 24 × 1 k-mesh. A 15 Å thick vacuum region is included to model the single layer structure. The second Piola-Kirchhoff (K-P) stress [27] is used for the 2D forces per length with units of N/m.
The mechanical deformations are applied in thee modes: uniaxial along the armchair direction (i.e., nearest neighbor direction), uniaxial along the zigzag direction (i.e., second nearest neighbor direction), and biaxial. Both tensile (namely positive strain) and compressive (i.e., negative) strains are examined in all the three modes. The fourteen non-zero elastic constants of g-SiC are obtained by a least-squares fit to stress-strain results [41]. This method has been widely used to study the mechanical properties of various graphene-like materials, including g-BN, g-GeC, g-GaN, g-TlN, g-ZnO, g-Boron, g-BNC, g-ZnS, g-MoS2, graphene-oxide, and g-InN [23,24,25,29,42,43,44,45,46,47,48,49,50,51,52,53,54].

3. Results and Discussion

3.1. Atomic Structure

We have fully relaxed the system with all degree of freedom including both the cell geometry and individual atoms. All atoms of a g-SiC mono-layer are coplanar, same as that of graphene. The lattice constant of a = 3.103 Å and 3.096 Å for SCAN and PBE, respectively. These values are in an excellent agreement with previous DFT studies [1,2,3]. The bond length of C-Si bond is 1.792 (1.787) Å from SCAN (PBE) calculations, which agrees well with previous DFT calculations of 1.77 Å [2] and 1.79 Å [1]. The bond angles of C-Si-C and Si-C-Si are 120 . Our result of atomic structure of g-SiC agrees well with literature [1,2,3]. The lattice constant of g-SiC mono-layers is compared to graphene and silicene as summarized in Table 1.

3.2. Strain Energy Profile

The energy profile is the stored potential energy under applied strain. The energy profile illustrates the change of energy with respect to strain, which can be used to describe the mechanical stability. When a strain is mechanically loaded, all the atoms in the supercell are allowed to relax. We have studied both the compression (negative strains) and tension (positive strains). The Lagrangian strains are used due to describe the applied large strains. For PBE calculations, the range is applied strain is −0.1 to 0.4. The tensile strain is up to 0.3 for SCAN calculations because the critical strains are less than 0.3. The step size of applied strain is 0.01 all three deformation modes. The reduced strain energy is defined as E s = ( E t o t E 0 ) / n for convenience of comparison between different systems including graphene and silicene [25], where E t o t is the total energy of the strained system, E 0 is the total energy of the strain-free system, and n = 6 is the number of atoms in the unitcell. The strain energy profile (the strain energy as a function of the strain) of g-SiC are displayed in Figure 2.
We conclude from the comparison in Figure 2 that SCAN predicts similar behaviors of strain energy profiles for uniaxial armchair, uniaxial zigzag and equibiaxial deformations to that of PBE. The critical strains (denoted by dashed lines) in strain energy profile are the ultimate tensile strains, where the second derivative of the strain energy to the strain is zero. The critical strains of the three deformation modes are summarized in Table 1 for comparison. The armchair critical strain from SCAN calculation is identical to that of PBE calculation. However, the zigzag and biaxial critical strain of SCAN is 0.01 larger than that of PBE. The difference is trivial because the error bar in strain is 0.01, same as the minimum increment strain step in this study.

3.3. Stress-Strain Curves

The second P-K stress is used due to the two-dimension nature of the system. The stress-strain relationships of g-SiC are shown in Figure 3 for three deformation modes, which are uniaxial strains along the armchair and zigzag directions and biaxial strains. The ultimate tensile strength is the maximum stress of the stress-strain curve, which means the upper limit that a material can withstand in a tensile test. The corresponding strain is named the ultimate tensile strain. Therefore, the system under strains beyond the ultimate strains is unstable. The ultimate strain serves as an upper limit of elastic stability. As a result, the ultimate tensile strength is a key parameter for strain engineering and applications. As listed in Table 1, the ultimate tensile stress at the three deformation modes are compared to that of graphene [45] and silicene [49].
The ultimate strains of SCAN are similar to that of PBE. However, the ultimate tensile stress of SCAN are larger than that of PBE, which is 6.7%, 9.0%, and 9.9% for armchair, zigzag, and biaxial deformation, respectively.

3.4. High Order Elastic Constants

The stress-strain relationship is the central part of the elastic theory of continuum mechanics, where the elastic constants are the key parameters. These elastic constants are critical to predict the mechanical behaviors at macro-scale. Among them, the high order elastic constants are essential in nonlinear elastic mechanics. We have predicted the 14 non-zero elastic constants from SCAN calculations, which includes normal and high order elastic constants. There are 2 SOECs (second order elastic constants), 3 TOECs (third order elastic constants), 4 FOECs (fourth order elastic constants), and 5 FFOECs (fifth order elastic constants) as listed in Table 1. The SOECs is for linear elastic mechanics and the high order elastic constants (TOECs, FOECs, and FFOECs) are for nonlinear elastic mechanics. Our results of elastic constants and the modulus from SCAN calculations are compared to that of PBE calculations as summarized in Table 1, along with that of graphene and silicene. It is worth noting that the in-plane Young’s modulus Y s are defined as Y s = ( C 11 2 C 12 2 ) / C 11 and Poisson ratio is computed as ν = C 12 / C 11 .

3.5. Pressure Effect on Elastic Constants

One of the applications of the high order (>order 2) elastic constants is to describe the pressure effects on mechanical properties. The second-order elastic moduli of g-SiC are not constants when the in-plane pressure is introduced. Instead, second-order elastic moduli are a function of in-plane pressure. The pressure-dependent second-order elastic moduli ( C 11 ˜ , C 12 ˜ , C 22 ˜ ) can be obtained from second-order elastic moduli at zero pressure ( C 11 , C 12 , C 22 , Y s , and ν ), and third order elastic constants ( C 111 , C 112 , and C 222 ) as [30,41]:
C 11 ˜ = C 11 ( C 111 + C 112 ) 1 ν Y s P ,
C 22 ˜ = C 11 C 222 1 ν Y s P ,
C 12 ˜ = C 12 C 112 1 ν Y s P .
The results of the pressure effects are illustrated in Figure 4 for both SCAN and PBE calculations. The predicted second-order elastic moduli, in-plane Young’s modulus, and Poisson ratio as a function of the in-plane pressure from SCAN are similar to that from PBE calculations.

3.6. Pressure Effect on the Soundwave Speed

One of the applications of the high order elastic constants is to predict the pressure effects on the speed of soundwaves propagating in the g-SiC. We have investigated the speed of anharmonic sound waves as a function of the in-plane pressure. The pressure-dependent soundwaves are interesting in many applications. Depending on the deformation mode during propagation, there are two main classes of soundwaves: compressive waves or p-waves, and shear waves or s-waves. With the pressure-dependent elastic constants, the speed of these two types of soundwave can be computed as
u p = Y s ˜ ( 1 ν ˜ ) ρ m ( 1 + ν ˜ ) ( 1 2 ν ˜ ) ,
u s = C 12 ˜ ρ m .
The dependences of both u p and u s on in-plane pressures, are plotted in Figure 5. It is well demonstrated that the speed of both p-wave and s-wave predicted from SCAN calculation is very similar to that of PBE.
Besides the speed along, the velocity ratio of the two, u p / u s , is a key parameter in the mechanics of materials [41]. The velocity ratio has been used for many purposes, such as a lithology indicator, identifying pore fluid, determining degree of consolidation, and predicting velocities [55]. The velocity ratio relates to the Poisson’s ratio as
u p u s = 1 ν ˜ ( 1 + ν ˜ 2 1 2 ν ˜ ) .
It is well demonstrated in Figure 5 that the velocity ratio predicted from SCAN is similar to that of PBE calculations.

3.7. Elastic Stability

We have also examined the elastic stability according to the Born elastic stability criteria [56]. The necessary and sufficient elastic stability conditions are illustrated in the following equation reflecting that the strain energy density is positive for a strain [57]:
Φ ( η ) > 0 ,   η 0 .
Our results of SCAN calculation of strain energy profiles shown in Figure 2 implies the elastic stabilities of g-SiC mono-layers. The same conclusion can be drawn from PBE calculations.
In terms of elastic stiffness tensor, furthermore, we examine elastic stability from the aspect of elastic stiffness tensor. The criteria are expressed in the following equations for a 2D hexagonal lattice:
C 11 > 0 ; C 11 > | C 12 | .
Both the SCAN and PBE results of the second order elastic constants summarized in Table 1 suggest the satisfactory of the stiffness tensor criteria (Equation (8)) for g-SiC mono-layers.

4. Conclusions

We have assessed the performance of the SCAN calculations on the nonlinear mechanical properties and mechanical stabilities of g-SiC under large strains with point-by-point comparison to that of PBE calculations, including ultimate strength, ultimate strain, 14 non-trivial elastic constants up to 5th order, strain energy profile, stress-strain behaviors, and pressure effect on normal modulus, shear modulus, Poisson ratio, p-wave, s-wave, and velocity ratio, under large deformations in three modes. We demonstrate that results from SCAN are similar to that of PBE calculations. The armchair ultimate strain from SCAN is identical to that of PBE. The difference in zigzag and biaxial ultimate strain is only 0.01, within the precision of strain we applied in this study. The differences in ultimate stress are all less than 10%, which fall in the precision of a DFT calculations.
The agreement between SCAN and PBE calculations is strikingly good in prediction of the pressure on the normal modulus, shear modulus, Poisson ratio, p-wave, s-wave, and velocity ratio. The results of the elastic stability from SCAN and PBE are also identical. From these comparisons, we could conclude that the SCAN calculations have very similar performance in predicting the nonlinear mechanics. In other words, the improvement in precision in predicting nonlinear mechanics from SCAN calculation over PBE calculation is minor. It is worth noting that the computing demand of SCAN calculation is about double that of PBE calculation. Considering the considerable increase of computing demand, the improvement of the precision seems not worthwhile. The PBE results are still reliable in predicting mechanical properties as the standard of DFT calculations.

Funding

Authors would like to acknowledge the support provided by the Deanship of Scientific Research (DSR) at King Fahd University of Petroleum and Minerals (KFUPM) through project No. SR191013 and K.A.CARE Research Fellowship.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated during and/or analyzed during the current study are available from the corresponding author on request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shi, Z.; Zhang, Z.; Kutana, A.; Yakobson, B.I. Predicting Two-Dimensional Silicon Carbide Mono layers. ACS Nano 2015, 9, 9802–9809. [Google Scholar] [CrossRef]
  2. Şahin, H.; Cahangirov, S.; Topsakal, M.; Bekaroglu, E.; Akturk, E.; Senger, R.T.; Ciraci, S. Monolayer honeycomb structures of group-IV elements and III-V binary compounds: First-principles calculations. Phys. Rev. B 2009, 80, 155453. [Google Scholar] [CrossRef] [Green Version]
  3. Lu, T.Y.; Liao, X.X.; Wang, H.Q.; Zheng, J.C. Tuning the indirect-direct band gap transition of SiC, GeC and SnC monolayer in a graphene-like honeycomb structure by strain engineering: A quasiparticle GW study. J. Mater. Chem. 2012, 22, 10062–10068. [Google Scholar] [CrossRef]
  4. Ganguly, G.; De, S.C.; Ray, S.; Barua, A.K. Polycrystalline silicon-carbide films deposited by low-power radiofrequency plasma decomposition of SiF4-CF4-H2 Gas-mixtures. J. Appl. Phys. 1991, 69, 3915–3923. [Google Scholar] [CrossRef]
  5. Lin, S.S. Light-Emitting Two-Dimensional Ultrathin Silicon Carbide. J. Phys. Chem. C 2012, 116, 3951–3955. [Google Scholar] [CrossRef]
  6. Wu, R.; Zhou, K.; Yue, C.Y.; Wei, J.; Pan, Y. Recent progress in synthesis, properties and potential applications of SiC nanomaterials. Prog. Mater Sci. 2015, 72, 1–60. [Google Scholar] [CrossRef]
  7. Lin, X.; Lin, S.; Xu, Y.; Hakro, A.A.; Hasan, T.; Zhang, B.; Yu, B.; Luo, J.; Li, E.; Chen, H. Ab initio study of electronic and optical behavior of two-dimensional silicon carbide. J. Mater. Chem. C 2013, 1, 2131–2135. [Google Scholar] [CrossRef]
  8. Hsueh, H.C.; Guo, G.Y.; Louie, S.G. Excitonic effects in the optical properties of a SiC sheet and nanotubes. Phys. Rev. B 2011, 84, 085404. [Google Scholar] [CrossRef] [Green Version]
  9. Wu, I.J.; Guo, G.Y. Optical properties of SiC nanotubes: An ab initio study. Phys. Rev. B 2007, 76, 035343. [Google Scholar] [CrossRef] [Green Version]
  10. Sun, L.; Li, Y.; Li, Z.; Li, Q.; Zhou, Z.; Chen, Z.; Yang, J.; Hou, J.G. Electronic structures of SiC nanoribbons. J. Chem. Phys. 2008, 129, 174114. [Google Scholar] [CrossRef]
  11. Bekaroglu, E.; Topsakal, M.; Cahangirov, S.; Ciraci, S. First-principles study of defects and adatoms in silicon carbide honeycomb structures. Phys. Rev. B 2010, 81, 075433. [Google Scholar] [CrossRef] [Green Version]
  12. Zhou, L.J.; Zhang, Y.F.; Wu, L.M. SiC2 Siligraphene and Nanotubes: Novel Donor Materials in Excitonic Solar Cells. Nano Lett. 2013, 13, 5431–5436. [Google Scholar] [CrossRef] [PubMed]
  13. Yang, C.; Xie, Y.; Liu, L.M.; Chena, Y. Versatile electronic properties and exotic edge states of single-layer tetragonal silicon carbides. Phys. Chem. Chem. Phys. 2015, 17, 11211–11216. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Zhang, P.; Xiao, B.B.; Hou, X.L.; Zhu, Y.F.; Jiang, Q. Layered SiC Sheets: A Potential Catalyst for Oxygen Reduction Reaction. Sci. Rep. 2014, 4, 3821. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Lin, X.; Lin, S.; Xu, Y.; Chen, H. Electronic structures of multilayer two-dimensional silicon carbide with oriented misalignment. J. Mater. Chem. C 2015, 3, 9057–9062. [Google Scholar] [CrossRef]
  16. Yu, M.; Jayanthi, C.S.; Wu, S.Y. Geometric and electronic structures of graphitic-like and tubular silicon carbides: Ab-initio studies. Phys. Rev. B 2010, 82, 075407. [Google Scholar] [CrossRef] [Green Version]
  17. Guilhon, I.; Teles, L.K.; Marques, M.; Pela, R.R.; Bechstedt, F. Influence of structure and thermodynamic stability on electronic properties of two-dimensional SiC, SiGe, and GeC alloys. Phys. Rev. B 2015, 92, 075435. [Google Scholar] [CrossRef]
  18. Wang, N.; Tian, Y.; Zhao, J.; Jin, P. CO oxidation catalyzed by silicon carbide (SiC) monolayer: A theoretical study. J. Mol. Graph. Modell. 2016, 66, 196–200. [Google Scholar] [CrossRef]
  19. Deng, S.; Sumant, A.V.; Berry, V. Strain engineering in two-dimensional nanomaterials beyond graphene. Nano Today 2018, 22, 14–35. [Google Scholar] [CrossRef]
  20. Ariga, K.; Mori, T.; Hill, J.P. Mechanical Control of Nanomaterials and Nanosystems. Adv. Mater. 2012, 24, 158–176. [Google Scholar] [CrossRef]
  21. Dai, Y.; Li, S.; Gao, H.; Wang, W.; Sun, Q.; Peng, Q.; Gui, C.; Qian, Z.; Liu, S. Stress evolution in AlN and GaN grown on Si(111): Experiments and theoretical modeling. J. Mater. Sci. 2015, 27, 2004–2013. [Google Scholar] [CrossRef]
  22. Dai, Y.; Wang, W.; Gui, C.; Wen, X.; Peng, Q.; Liu, S. A first-principles study of the mechanical properties of AlN with Raman verification. Comput. Mater. Sci. 2016, 112, 342–346. [Google Scholar] [CrossRef]
  23. Peng, Q.; Chen, X.J.; Ji, W.; De, S. Chemically Tuning Mechanics of Graphene by BN. Adv. Eng. Mater. 2013, 15, 718–727. [Google Scholar] [CrossRef]
  24. Peng, Q.; Liang, C.; Ji, W.; De, S. A First-principles Study of the Mechanical Properties of g-GeC. Mech. Mater. 2013, 64, 135–141. [Google Scholar] [CrossRef]
  25. Peng, Q.; Sun, X.; Wang, H.; Yang, Y.; Wen, X.D.; Huang, C.; Liu, S.; De, S. Theoretical prediction of a graphene-like structure of Indium Nitride: A promising excellent material for optoelectronics. Appl. Mater. Today 2017, 7, 169. [Google Scholar] [CrossRef]
  26. Xu, Z.; Li, Y.; Liu, Z. Controlling electronic and optical properties of layered SiC and GeC sheets by strain engineering. Mater. Des. 2016, 108, 333–342. [Google Scholar] [CrossRef]
  27. Peng, Q.; Zamiri, A.R.; Ji, W.; De, S. Elastic Properties of Hybrid Graphene/Boron Nitride Monolayer. Acta Mech. 2012, 223, 2591–2596. [Google Scholar] [CrossRef] [Green Version]
  28. Peng, Q.; Ji, W.; De, S. Mechanical Properties of Graphyne Monolayer: A First-Principles Study. Phys. Chem. Chem. Phys. 2012, 14, 13385–13391. [Google Scholar] [CrossRef] [Green Version]
  29. Peng, Q.; Dearden, A.K.; Chen, X.J.; Huang, C.; Wen, X.; De, S. Peculiar pressure effect on Poisson ratio of graphone as a strain damper. Nanoscale 2015, 7, 9975–9979. [Google Scholar] [CrossRef]
  30. Peng, Q. High-order nonlinear mechanical properties of g-SiC. Mech. Mater. 2020, 148, 103473. [Google Scholar] [CrossRef]
  31. Sousa, S.F.; Fernandes, P.A.; Ramos, M.J. General Performance of Density Functionals. J. Phys. Chem. A 2007, 111, 10439–10452. [Google Scholar] [CrossRef] [PubMed]
  32. Perdew, J.P.; Schmidt, K. Jacob’s ladder of density functional approximations for the exchange-correlation energy. AIP Conf. Proc. 2001, 577, 1–20. [Google Scholar]
  33. Perdew, J.P.; Ruzsinszky, A.; Tao, J.; Staroverov, V.N.; Scuseria, G.E.; Csonka, G.I. Prescription for the design and selection of density functional approximations: More constraint satisfaction with fewer fits. J. Chem. Phys. 2005, 123, 062201. [Google Scholar] [CrossRef] [PubMed]
  34. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Sun, J.; Ruzsinszky, A.; Perdew, J.P. Strongly Constrained and Appropriately Normed Semilocal Density Functional. Phys. Rev. Lett. 2015, 115, 036402. [Google Scholar] [CrossRef] [Green Version]
  36. Sun, J.; Remsing, R.C.; Zhang, Y.; Sun, Z.; Ruzsinszky, A.; Peng, H.; Yang, Z.; Paul, A.; Waghmare, U.; Wu, X.; et al. Accurate first-principles structures and energies of diversely bonded systems from an efficient density functional. Nat. Chem. 2016, 8, 831–836. [Google Scholar] [CrossRef]
  37. Marianetti, C.A.; Yevick, H.G. Failure Mechanisms of Graphene under Tension. Phys. Rev. Lett. 2010, 105, 245502. [Google Scholar] [CrossRef] [Green Version]
  38. Kresse, G.; Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 1993, 47, 558. [Google Scholar] [CrossRef]
  39. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef] [Green Version]
  40. Peng, Q.; Wang, G.; Liu, G.R.; De, S. Van der Waals Density Functional Theory vdW-DFq for Semihard Materials. Crystals 2019, 9, 243. [Google Scholar] [CrossRef] [Green Version]
  41. Peng, Q.; Ji, W.; De, S. Mechanical properties of the hexagonal boron nitride monolayer: Ab initio study. Comput. Mater. Sci. 2012, 56, 11–17. [Google Scholar] [CrossRef] [Green Version]
  42. Peng, Q.; Liang, C.; Ji, W.; De, S. A First Principles Investigation of the Mechanical Properties of g-TlN. Model. Numer. Simul. Mater. Sci. 2012, 2, 76–84. [Google Scholar] [CrossRef] [Green Version]
  43. Peng, Q.; Liang, C.; Ji, W.; De, S. A First Principles Investigation of the Mechanical Properties of g-ZnO: The Graphene-like Hexagonal Zinc Oxide Monolayer. Comput. Mater. Sci. 2013, 68, 320–324. [Google Scholar] [CrossRef]
  44. Peng, Q.; Ji, W.; De, S. First-Principles study of the Effects of Mechanical Strains on the Radiation Hardness of Hexagonal Boron Nitride Monolayers. Nanoscale 2013, 5, 695–703. [Google Scholar] [CrossRef] [PubMed]
  45. Peng, Q.; Liang, C.; Ji, W.; De, S. A Theoretical Analysis of the Effect of the Hydrogenation of Graphene to Graphane on Its Mechanical Properties. Phys. Chem. Chem. Phys. 2013, 15, 2003–2011. [Google Scholar] [CrossRef] [PubMed]
  46. Peng, Q.; Chen, X.J.; Liu, S.; De, S. Mechanical Stabilities and Properties of Graphene-like Aluminum Nitride Predicted from First-principles Calculations. RSC Adv. 2013, 3, 7083–7092. [Google Scholar] [CrossRef]
  47. Peng, Q.; Liang, C.; Ji, W.; De, S. Mechanical Properties of g-GaN: A First Principles Study. Appl. Phys. A 2013, 13, 483–490. [Google Scholar] [CrossRef]
  48. Peng, Q.; Chen, Z.; De, S. A density functional theory study of the mechanical properties of graphane with van der Waals corrections. Mech. Adv. Mater. Struct. 2015, 22, 717–721. [Google Scholar] [CrossRef]
  49. Peng, Q.; Wen, X.; De, S. Mechanical stabilities of Silicene. RSC Adv. 2013, 3, 13772–13781. [Google Scholar] [CrossRef]
  50. Peng, Q.; De, S. Outstanding mechanical properties of monolayer MoS2 and its application in elastic energy storage. Phys. Chem. Chem. Phys. 2013, 15, 19427–19437. [Google Scholar] [CrossRef]
  51. Peng, Q.; De, S. Mechanical properties and instabilities of ordered graphene oxide C6O monolayer. RSC Adv. 2013, 3, 24337–24344. [Google Scholar] [CrossRef]
  52. Peng, Q.; Han, L.; Wen, X.; Liu, S.; Chen, Z.; Lian, J.; De, S. Mechanical properties and stabilities of α-Boron monolayers. Phys. Chem. Chem. Phys. 2015, 17, 2160–2168. [Google Scholar] [CrossRef] [PubMed]
  53. Peng, Q.; Han, L.; Wen, X.; Liu, S.; Chen, Z.; Lian, J.; De, S. Mechanical properties and stabilities of g-ZnS monolayers. RSC Adv. 2015, 5, 11240–11247. [Google Scholar] [CrossRef]
  54. Peng, Q.; Han, L.; Lian, J.; Wen, X.; Liu, S.; Chen, Z.; Koratkar, N.; De, S. Mechanical degradation of graphene by epoxidation: Insights from first-principles calculations. Phys. Chem. Chem. Phys. 2015, 17, 19484–19490. [Google Scholar] [CrossRef] [PubMed]
  55. Lee, M.W. Velocity Ratio and Its Application to Predicting Velocities; US Department of the Interior, US Geological Survey: Reston, VA, USA, 2003.
  56. Mouhat, F.; Coudert, F.X. Necessary and sufficient elastic stability conditions in various crystal systems. Phys. Rev. B 2014, 90, 224104. [Google Scholar] [CrossRef] [Green Version]
  57. Mazdziarz, M. Comment on ‘The Computational 2D Materials Database: High-throughput modeling and discovery of atomically thin crystals’. 2D Mater. 2019, 6, 048001. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Geometry and orientation. The atomistic structure of Graphene-like SiC (g-SiC ) is displayed in (a) top-view the plane with the unit cell marked. (b) The atomic configuration of the conventional six-atom unit-cell and the coordinate orientations. The large (light green) and small (yellow) ball represents the Si and C atoms, respectively.
Figure 1. Geometry and orientation. The atomistic structure of Graphene-like SiC (g-SiC ) is displayed in (a) top-view the plane with the unit cell marked. (b) The atomic configuration of the conventional six-atom unit-cell and the coordinate orientations. The large (light green) and small (yellow) ball represents the Si and C atoms, respectively.
Crystals 11 00120 g001
Figure 2. Energy profile of g-SiC Strongly Constrained and Appropriately Normed (SCAN) (square) and Perdew-Burke-Ernzerhof (PBE) (solid line) results of the reduced strain-energy as a function of applied strain. Three deformation modes are applied: uniaxial strain in armchair (a) and zigzag (b) directions, and equibiaxial (c). The vertical dashed lines indicate the critical strain.
Figure 2. Energy profile of g-SiC Strongly Constrained and Appropriately Normed (SCAN) (square) and Perdew-Burke-Ernzerhof (PBE) (solid line) results of the reduced strain-energy as a function of applied strain. Three deformation modes are applied: uniaxial strain in armchair (a) and zigzag (b) directions, and equibiaxial (c). The vertical dashed lines indicate the critical strain.
Crystals 11 00120 g002
Figure 3. Comparison of stress-strain relationship. SCAN (left panels) calculation of g-SiC for the stress-strain responses under the three deformation modes of armchair (a), zigzag (b), and biaxial (c) strain. The counterparts of PBE calculations (right panels) are in (df), respectively. Σ 1 and Σ 2 stands for the x and y component of stress, respectively. The Density Functional Theory (DFT) calculations and the fitting for continuum mechanics are marked by “DFT” and “Cont”, respectively. The regions of strain are reflected by different colors. The tensile region (positive strain) consists of two parts of elastic region (light green) and unstable region (yellow). The compressive region (negative strain) is displayed in light blue. The elastic tensile region is divided into linear elastic (pink) and nonlinear (dark violet) regions.
Figure 3. Comparison of stress-strain relationship. SCAN (left panels) calculation of g-SiC for the stress-strain responses under the three deformation modes of armchair (a), zigzag (b), and biaxial (c) strain. The counterparts of PBE calculations (right panels) are in (df), respectively. Σ 1 and Σ 2 stands for the x and y component of stress, respectively. The Density Functional Theory (DFT) calculations and the fitting for continuum mechanics are marked by “DFT” and “Cont”, respectively. The regions of strain are reflected by different colors. The tensile region (positive strain) consists of two parts of elastic region (light green) and unstable region (yellow). The compressive region (negative strain) is displayed in light blue. The elastic tensile region is divided into linear elastic (pink) and nonlinear (dark violet) regions.
Crystals 11 00120 g003
Figure 4. Pressure effect on elastic constants. The SCAN calculation predicts in-plane Young’s modulus, the second-order elastic moduli, and Poisson ratio with respect to the applied in-plane pressure. The comparison of SCAN (marked) and PBE (solid lines) is illustrated in (a) the normal modulus of C 11 ˜ , C 22 ˜ and in-plane Young’s modulus, (b) shear modulus of C 12 ˜ , (c) Poisson ratio C 12 ˜ .
Figure 4. Pressure effect on elastic constants. The SCAN calculation predicts in-plane Young’s modulus, the second-order elastic moduli, and Poisson ratio with respect to the applied in-plane pressure. The comparison of SCAN (marked) and PBE (solid lines) is illustrated in (a) the normal modulus of C 11 ˜ , C 22 ˜ and in-plane Young’s modulus, (b) shear modulus of C 12 ˜ , (c) Poisson ratio C 12 ˜ .
Crystals 11 00120 g004
Figure 5. Velocities of sound waves. The SCAN calculation predicts the speed of (a) p-waves and (b) s-waves with respect to the applied in-plane pressure, compared with that of PBE calculation. (c) The velocity ratio u p / u s with respect to the applied in-plane pressure. The red squares are for results of SCAN calculations and the blue solid lines are for that of PBE calculations.
Figure 5. Velocities of sound waves. The SCAN calculation predicts the speed of (a) p-waves and (b) s-waves with respect to the applied in-plane pressure, compared with that of PBE calculation. (c) The velocity ratio u p / u s with respect to the applied in-plane pressure. The red squares are for results of SCAN calculations and the blue solid lines are for that of PBE calculations.
Crystals 11 00120 g005
Table 1. Mechanical properties The results of SCAN calculations of lattice constant a, Poisson’s ration ν , in-plane stiffness Y s , ultimate tensile strength in three modes ( Σ u a , Σ u z and Σ u b ), ultimate tensile strains ( η u a , η u z , η u b ), the fourteen independent components of g-SiC, are compared to that of PBE calculations and that of silicene and graphene.
Table 1. Mechanical properties The results of SCAN calculations of lattice constant a, Poisson’s ration ν , in-plane stiffness Y s , ultimate tensile strength in three modes ( Σ u a , Σ u z and Σ u b ), ultimate tensile strains ( η u a , η u z , η u b ), the fourteen independent components of g-SiC, are compared to that of PBE calculations and that of silicene and graphene.
g-SiC (SCAN)g-SiC (PBE) [30]Silicene [49]Graphene [45]
a (Å)3.1033.0963.8652.468
Y s 174160.163.8340.8
ν 0.3040.3140.3250.178
Σ u a 17.616.56.028.6
η u a 0.20.20.170.19
Σ u z 18.216.75.930.4
η u z 0.260.250.210.23
Σ u b 17.816.26.232.1
η u b 0.200.190.170.23
SOEC C 11 191.3177.671.3352.0
C 12 58.255.723.262.6
TOEC C 111 −1487.6−1390.5−397.6−3089.7
C 112 −398.8−391.9−14.1−453.8
C 222 −1261.5−1180.9−318.9−2928.1
FOEC C 1111 99909165−83021,927
C 1112 36943604−3092731
C 1122 9421162−50913888
C 2222 34133000−62918,779
FFOEC C 11111 −41,367−37,19620,641−118,791
C 11112 −5129−60756923−19,173
C 11122 −23,635−19,30011,681−15,863
C 12222 −26,812−21,523−7593−27,463
C 22222 −26,130−14,127−29,735−134,752
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peng, Q. Performance of SCAN Meta-GGA Functionals on Nonlinear Mechanics of Graphene-Like g-SiC. Crystals 2021, 11, 120. https://doi.org/10.3390/cryst11020120

AMA Style

Peng Q. Performance of SCAN Meta-GGA Functionals on Nonlinear Mechanics of Graphene-Like g-SiC. Crystals. 2021; 11(2):120. https://doi.org/10.3390/cryst11020120

Chicago/Turabian Style

Peng, Qing. 2021. "Performance of SCAN Meta-GGA Functionals on Nonlinear Mechanics of Graphene-Like g-SiC" Crystals 11, no. 2: 120. https://doi.org/10.3390/cryst11020120

APA Style

Peng, Q. (2021). Performance of SCAN Meta-GGA Functionals on Nonlinear Mechanics of Graphene-Like g-SiC. Crystals, 11(2), 120. https://doi.org/10.3390/cryst11020120

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

Article Metrics

Back to TopTop