Mechanical properties of defective single-layered graphene sheets via molecular dynamics simulation

In this paper, the effects of two main types of structural defects, i.e. Stone–Wales and single vacancy, on the mechanical properties of single-layered graphene sheets (SLGSs) are investigated. To this end, molecular dynamics simulations based on the Tersoff–Brenner potential function and Nose–Hoover thermostat technique are implemented. The results obtained have revealed that the presence of defects significantly reduces the failure strain and the intrinsic strength of SLGSs, while it has a slight effect on Young’s modulus. Furthermore, the examination of loading in both armchair and zigzag directions demonstrated that SLGSs are slightly stronger in the armchair direction and defects have lower effect in this direction. Considering the fracture mechanism, the failure process of defective and perfect graphene sheets is also presented.     

Vibration analysis of defective graphene sheets using nonlocal elasticity theory

Many papers have studied the free vibration of graphene sheets. However, all this papers assumed their atomic structure free of any defects. Nonetheless, they actually contain some defects including single vacancy, double vacancy and Stone-Wales defects. This paper, therefore, investigates the free vibration of defective graphene sheets, rather than pristine graphene sheets, via nonlocal elasticity theory. Governing equations are derived using nonlocal elasticity and the first-order shear deformation theory (FSDT). The influence of structural defects on the vibration of graphene sheets is considered by applying the mechanical properties of defective graphene sheets. Afterwards, these equations solved using generalized differential quadrature method (GDQ). The small-scale effect is applied in the governing equations of motion by nonlocal parameter. The effects of different defect types are inspected for graphene sheets with clamped or simply-supported boundary conditions on all sides. It is shown that the natural frequencies of graphene sheets decrease by introducing defects to the atomic structure. Furthermore, it is found that the number of missing atoms, shapes and distributions of structural defects play a significant role in the vibrational behavior of graphene. The effect of vacancy defect reconstruction is also discussed in this paper.     

Effects of structural defects on the compressive buckling of boron nitride nanotubes

In this paper, we examined the buckling of perfect and defective armchair boron nitride nanotubes with three types of vacancy defects, i.e. B- and N- single vacancy defects and B–N- double vacancy defect, using molecular dynamics simulations. To this end, all systems were modeled with a Tersoff-type potential, which is able to accurately describe covalent bonding of BN systems. We applied external uniaxial compressive forces to the nanotubes in vacuum and derived the critical buckling loads and strains, at room temperature in an NVT-ensemble. Our results showed significant differences between the critical buckling strengths of pristine and defective nanotubes. The resistance to axial buckling decreased with the introduction of one vacancy defect, and the B–N- double vacancy was the most seriously damaged structure, followed by B-vacancy and N-vacancy defects. Furthermore, the B-vacancy was shown to have the most significant effect on the decrease of the critical buckling strain. This can be attributed to the excessive asymmetries and perturbations induced in the structure of the nanotube and the local deformations around the defective site around the B-vacancy, even before loading. Moreover, results show that reduction in the buckling strength of the nanotube due to the presence of more than one B-vacancy defect depends on their distribution. If the two or three defects are close to each other, they act as a single point of weakness and the critical buckling load is only slightly reduced (similar to the existence of only one vacancy defect). However, if the defects are at more distant points, the critical buckling load may experience a higher decrease. Results show that vacancy defects play a critical role in the compressive buckling performance of boron nitride nanotubes and special attention must be paid to the presence of structural defects when designing members against buckling, especially for micro- and nano-electro-mechanical systems. On the other hand, defect engineering is a great means for tailoring the buckling strength of boron nitride nanotubes, in cases where the nanotube is expected to absorb energy through compressive buckling deformation and is not designed against, but for buckling.     

Atomistic finite element model for axial buckling and vibration analysis of single-layered graphene sheets

In this article, an atomistic model is developed to study the buckling and vibration characteristics of single-layered graphene sheets (SLGSs). By treating SLGSs as space-frame structures, in which the discrete nature of graphene sheets is preserved, they are modeled using three-dimensional elastic beam elements for the bonds. The elastic moduli of the beam elements are determined via a linkage between molecular mechanics and structural mechanics. Based on this model, the critical compressive forces and fundamental natural frequencies of single-layered graphene sheets with different boundary conditions and geometries are obtained and then compared. It is indicated that the compressive buckling force decreases when the graphene sheet aspect ratio increases. At low aspect ratios, the increase of aspect ratios will result in a significant decrease in the critical buckling load. It is also indicated that increasing aspect ratio at a given side length results in the convergence of buckling envelops associated with armchair and zigzag graphene sheets. The influence of boundary conditions will be studied for different geometries. It will be shown that the influence of boundary conditions is not significant for sufficiently large SLGSs.     

Diffusion, coalescence, and reconstruction of vacancy defects in graphene layers

Diffusion, coalescence, and reconstruction of vacancy defects in graphene layers are investigated by tight-binding molecular dynamics (TBMD) simulations and by first principles total energy calculations. It is observed in the TBMD simulations that two single vacancies coalesce into a 5-8-5 double vacancy at the temperature of 3000 K, and it is further reconstructed into a new defect structure, the 555-777 defect, by the Stone-Wales type transformation at higher temperatures. First principles calculations confirm that the 555-777 defect is energetically much more stable than two separated single vacancies, and the energy of the 555-777 defect is also slightly lower than that of the 5-8-5 double vacancy. In TBMD simulation, it is also found that the four single vacancies reconstruct into two collective 555-777 defects which is the unit for the hexagonal haeckelite structure proposed by Terrones et al. [Phys. Rev. Lett. 84, 1716 (2000)].     

单、双原子空位缺陷对扶手椅型单层碳纳米管屈曲性能的不同影响

采用分子动力学方法,对完善和含缺陷扶手椅型单层碳纳米管进行轴向压缩的数值模拟,对比研究三种不同的温度环境下单、双原子空位缺陷对碳纳米管轴压变形性能的特殊影响.研究结果表明管壁缺陷显著降低了纳米管低温时的承载能力,由于单原子空位缺陷造成的特殊应力集中效应会引发纳米管过早发生局部屈曲,单原子缺陷管的屈曲强度反而小于双原子管的屈曲强度. 关键词： 分子动力学 碳纳米管 屈曲 缺陷     

Deformation and fracture in graphene with divacancies of the 555–777 type

The deformation and fracture of graphene sheets containing 555–777 defects have been investigated by molecular dynamics simulations. Each such defect is a divacancy forming a localized configuration of three pentagonal and three septangular cells of carbon atoms in a hexagonal graphene lattice. An emphasis is made on the influence of 555–777 defects in graphene on its mechanical characteristics (stress–strain curve, uniaxial tensile strength, and maximum elastic strain).     

Atomistic evaluation of the stress concentration factor of graphene sheets having circular holes

Stress concentration factor concept has been developed for single-layered graphene sheets (SLGSs) with circular holes through an atomistic point of view by the application of molecular structural mechanics (MSM) approach. In this approach the response of SLGSs against unidirectional tensile loading is matched to the response of a frame-like macro structure containing beam elements by making an equivalence between strain energies of beam elements in MSM and potential energies of chemical bonds of SLGSs. Both chirality and size effects are considered and the atomistic evaluation of stress concentration factor is performed for different sizes of circular holes. Also, molecular dynamics simulations are implemented to verify the existence and location of the predicted stress concentration. The results reveal that size effects and the diameters of circular holes have a significant influence on the stress concentration factor of SLGSs and armchair SLGSs show a larger value of stress concentration than zigzag ones.     

Influence of line defects on relaxation properties of graphene: A molecular dynamics study

The relaxation properties of single layer graphene sheets containing line defects were investigated using molecular dynamics simulation with AIROBE bond-order interatomic potential. The dynamic evolution of graphene sheets during relaxation condition was analyzed. The simulation results show that the single layer graphene sheets are not perfectly flat in an ideal state, and the graphene sheet shows a significant corrugations at the verge of sheet. The graphene sheet is bent with the line defects at the end of the sheet, and the extent of this bend also increases with the increase of the defect number. Furthemore, the graphene sheet transforms into a paraboloid with the line defects at the middle of the sheet.     

First-principles study of lithium intercalated bilayer graphene

Lithium intercalated bilayer graphene has been investigated using first-principles density functional theory calculations. Results show that there exist AB and AA stacking sequences for bilayer graphene in which the latter is more favorable for the Li storage and the former will evolve into the latter with the intercalation of Li ions. The relationship between the interlayer distance of two graphene sheets and the intercalated capacity of Li ions is discussed. It is found that structural defect is identified to store Li ions more favorably than pristine bilayer graphene and an isolated C atom vacancy in bilayer graphene can capture three Li ions between two graphene sheets.     

Nonlocal plate model for free vibrations of single-layered graphene sheets

Vibration analysis of single-layered graphene sheets (SLGSs) is investigated using nonlocal continuum plate model. To this end, Eringens's nonlocal elasticity equations are incorporated into the classical Mindlin plate theory for vibrations of rectangular nanoplates. In contrast to the classical model, the nonlocal model developed in this study has the capability to evaluate the natural frequencies of the graphene sheets with considering the size-effects on the vibrational characteristics of them. Solutions for frequencies of the free vibration of simply-supported and clamped SLGSs are computed using generalized differential quadrature (GDQ) method. Then, molecular dynamics (MD) simulations for the free vibration of various SLGSs with different values of side length and chirality are employed, the results of which are matched with the nonlocal model ones to derive the appropriate values of the nonlocal parameter relevant to each boundary condition. It is found that the value of the nonlocal parameter is independent of the magnitude of the geometrical variables of the system.     

Geometric stability and adsorption property of hydroxyl group on graphene sheets

The stable geometrics and adsorption behaviors of hydroxyl (OH) groups on graphene sheets are investigated using the first-principles calculations. The single hydroxyl adatom has small adsorption energy and diffusion barrier on pristine graphene. The binding strength of the hydroxyl group increases with the coverage, and the aggregations of the hydroxyl groups reduce the structural bucking of graphene sheet. On the graphene with single vacancy (SV-graphene), the large trapping zones mean the adsorbed OH would be easily trapped at the vacancy site. The hydroxyl groups prefer to aggregate on graphene surfaces and form the water molecule, leaving the epoxy group on pristine graphene or oxygen dopant in SV-graphene, which is used to constitute the structural model of oxidized graphene. These results would provide us a useful reference to understand the atomic structure and adsorption property of functional groups on graphene sheets.     

Hydrogen adsorption and storage of Ca-decorated graphene with topological defects: A first-principles study

As a candidate for hydrogen storage medium, geometric stability and hydrogen capacity of Ca-decorated graphene with topological defects are investigated using the first-principle based on density functional theory (DFT), specifically for the experimentally realizable single carbon vacancy (SV), 585 double carbon vacancy (585 DCV) and 555–777 double carbon vacancy (555–777 DCV) defects. It is found that Ca atom can be stabilized on above defective graphenes since Ca׳s binding energy on vacancy defect is much larger than its cohesive energy. Up to six H₂ molecules can stably bind to a Ca atom on defective graphene with the average adsorption energies of 0.17–0.39 eV/H₂. The hybridization of the Ca-3d orbitals with H₂-σorbitals and the electrostatic interaction between the Ca cation and the induced H₂ dipole both contribute to the H₂ molecules binding. Double-side Ca-decorated graphene with 585 DCV and 555–777 DCV defects can theoretically reach a gravimetric capacity of 5.2 wt% hydrogen, indicating that Ca-decorated defective graphene can be used as a promising material for high density hydrogen storage.     

Density functional theory study on the adsorption of alkali metal ions with pristine and defected graphene sheet

The graphene-based materials along with the adsorption of alkali metal ions are suitable for energy conversion and storage applications. Hence in the present work, we have investigated the structural and electronic properties of pristine and defected graphene sheet upon the adsorption of alkali metal ions (Li⁺, Na⁺, and K⁺) using density functional theory (DFT) calculations. The presence of vacancies or vacancy defects enhances the adsorption of alkali ions than the pristine sheet. From the obtained results, it is found that the adsorption energy of Li⁺ on the vacancies defected graphene sheet is higher (3.05?eV) than the pristine (2.41?eV) and Stone–Wales (2.50?eV) defected sheets. Moreover, the pore radius of the pristine and defected graphene sheets are less affected by metal ions adsorption. The increase in energy gap upon the adsorption of metal ions is found to be high in the vacancy defected graphene than that of other sheets. The metal ions adsorption in the defective vacancy sheets has high charge transfer from metal ions to the graphene sheet. The bonding characteristic between the metal ions and graphene sheet are analysed using QTAIM analysis. The influence of alkali ions on the electronic properties of the graphene sheet is examined from the Total Density of States (TDOS) and Partial Density of States (PDOS).     

The characterization of electronic defect states of single and double carbon vacancies in graphene sheets using molecular density functional theory

ABSTRACT

A detailed picture of the electronic states manifolds of single- and double-vacancy defects in molecular models of graphene based on polycyclic aromatic hydrocarbons (PAHs) is presented. DFT calculations using various density functionals including long-range corrected ones have been performed for pyrene, circumpyrene and 7a,7z-periacene. It has been found for pyrene defect models that DFT results reproduced well the set of closely-spaced singlet and triplet states predicted by the CCSD(T) and previous MRCI?+?Q calculations, indicating the applicability of DFT for accessing the excited states manifolds also for larger graphene models. For the single-carbon vacancy defect, all structures have a triplet ground state. As expected, in the largest system, 7a,7z-periacene-1C, the lowest lying states are much closer in energy. For all double-vacancy defect structures, a significant rearrangement of the electronic states with increasing size of the sheet is observed. The closed-shell ¹A_g state in the smallest systems is destabilised in the extended 7a,7z-periacene system, which has the ³B_2u state as the ground state. As observed for the single-vacancy defect, the lowest lying states are closer in energy for the larger systems, since there are more π orbitals close in energy available. For all states, the formation of the bridging bonds for the double vacancy leads to distances shorter than for the single vacancy defect indicating a larger rigidity of the former structure which does not allow stronger distortions.     

石墨烯/BiOI纳米复合物电子结构和光学性质的第一性原理模拟计算

光催化材料在解决能源短缺和环境污染等问题方面具有广泛的应用前景,本文通过构建BiOI纳米薄膜并将其与石墨烯复合起来,得到具有较高的比表面积和良好的光催化活性的纳米复合物光催化材料.采用基于密度泛函理论的第一性原理方法分别计算了单层和双层BiOI纳米片及其与石墨烯复合结构的电子结构和光学性质,并考虑了BiOI中的Bi,O,I三种空位缺陷对电子结构和光学特性的影响.计算结果表明,由于BiOI和石墨烯之间的相互作用,在石墨烯和BiOI界面处自发发生电荷转移,形成电子-空穴对,且石墨烯衬底可有效提高BiOI对可见光的光吸收,提高其光催化活性.对空位缺陷的计算表明,Bi空位缺陷可促进石墨烯和BiOI之间的电荷转移,形成更多的层间电子-空穴对;相反,O和I空位缺陷则抑制层间电荷转移,减少电子-空穴对的生成.     

温度、应变率和空位缺陷对氮化硼纳米管压缩性能的影响

采用分子动力学方法,分别模拟了完好的和含有缺陷的氮化硼纳米管的轴向压缩过程。原子间的相互作用采用Tersoff多体势函数来描述。结果表明,同尺寸的锯齿型氮化硼纳米管的临界轴向压缩强度高于扶手型氮化硼纳米管,这与碳纳米管的研究结果一致。发现纳米管的压缩强度,如临界轴向内力在低温下受温度影响明显,并且和应变率的大小有关。然而,应变率对纳米管的弹性变形没有影响。另外,还发现空位缺陷降低了纳米管的力学性能。与完好的纳米管相比,含有缺陷的纳米管轴向压缩强度对于温度的影响并不敏感。     

单空位缺陷对石墨纳米带电子结构和输运性质的影响

基于第一性原理电子结构和输运性质计算,研究了单空位缺陷对单层石墨纳米带(包括zigzag型和armchair型带)电子性质的影响.研究发现,单空位缺陷使石墨纳米带在费米面上出现一平直的缺陷态能带;单空位缺陷的引入使zigzag型半导体性的石墨纳米带变为金属性,这在能带工程中有重要的应用价值;奇数宽度的armchair型石墨纳米带表现出金属特性,有着很好的导电性能,同时,偶数宽度的armchair型石墨带虽有金属性的能带结构,但却有类似半导体的伏安特性;单空位缺陷使得奇数宽度的armchair石墨纳米带导电 关键词： 石墨纳米带 单空位缺陷 电子结构 输运性质     

空位结构缺陷对C纳米管弹性性质的影响

用分子动力学方法对不同空位缺陷的扶手椅型与锯齿型单壁C纳米管杨氏弹性模量进行了计算和分析. 结果表明：扶手椅型（5, 5）, （10,10）和锯齿型（9, 0）, （18, 0） 纳米管在无缺陷时其杨氏模量分别为948,901和804,860 GPa. 随管径的增大,扶手椅型和锯齿型单壁C纳米管弹性模量分别减小和增大,表现出完全不同的变化规律. 随着C纳米管中单点空位缺陷的均匀增加,杨氏模量下降,当缺陷比率增加到一定程度时,杨氏模量下降骤然趋缓,形成一下降平台;双空位缺陷对C纳米管杨氏模量的影响与其分布方向有关;随单点空位缺陷间原子数的增加,在轴向上,杨氏模量下降到某一值小幅波动,而在周向上杨氏模量先下降,然后上升到某一稳定值. 随两单点空位缺陷的空间距离进一步增大,杨氏模量又呈微降趋势. 通过分子间σ键与π键特征及缺陷间近程电子云耦合作用规律与空位缺陷内部5-1DB缺陷的形成特点等理论对上述规律进行了分析. 关键词： 空位缺陷 C纳米管 分子动力学 杨氏模量     

Effect of triangular vacancy defect on thermal conductivity and thermal rectification in graphene nanoribbons

We investigate the thermal transport properties of armchair graphene nanoribbons (AGNRs) possessing various sizes of triangular vacancy defect within a temperature range of 200–600 K by using classical molecular dynamics simulation. The results show that the thermal conductivities of the graphene nanoribbons decrease with increasing sizes of triangular vacancy defects in both directions across the whole temperature range tested, and the presence of the defect can decrease the thermal conductivity by more than 40% as the number of removed cluster atoms is increased to 25 (1.56% for vacancy concentration) owing to the effect of phonon–defect scattering. In the meantime, we find the thermal conductivity of defective graphene nanoribbons is insensitive to the temperature change at higher vacancy concentrations. Furthermore, the dependence of temperatures and various sizes of triangular vacancy defect for the thermal rectification ration are also detected. This work implies a possible route to achieve thermal rectifier for 2D materials by defect engineering.