Molecular mechanics-based finite element analysis of graphene sheet and carbon nanotubes using the rebo potential

Molecular mechanics-based finite element (FE) models of graphene sheet and single-walled zigzag and armchair carbon nanotubes (CNTs) are developed on the basis of the assumption that the carbon nanostructures, when loaded, behave like frame structures. The behavior of carbon–carbon bonds, which are represented by beam elements, is simulated using the many-body second generation reactive empirical bond order (REBO) potential. By means of the FE models, the tensile behavior of carbon nanostructures is simulated. The FE models are verified against molecular dynamics simulations. The computed results in terms of tensile stress–strain curves and fracture patterns are compared with results obtained using the pairwise modified-Morse potential. Different tensile properties and fracture patterns are predicted using the two potentials. This is mainly attributed to the deviations in the force–bond length curves and to the contribution of bond angle variation which is present in REBO. The present work is the first attempt to implement the REBO potential into a continuum model of carbon nanostructures and paves the way for a more systematic incorporation of atomistic simulation methods into continuum models.

Elastic moduli of boron nitride, aluminium nitride and gallium nitride nanotubes using second generation reactive empirical bond order potential

Purpose – The purpose of this paper is to study elastic properties of III-V nitride nanotubes (NNTs) using second generation (REBO) potential. Design/methodology/approach – In the present research paper elastic properties of BN, AlN and GaN nanotubes have been investigated, using the second generation REBO potential by Brenner and co-workers, which is a bond order potential earlier used for carbon nanostructures successfully. In the present calculation, the same form of potential is used with adjusted parameters for h-BN, h-AlN and h-GaN. In all these cases the authors have considered graphite like network and strongly polar nature of these atoms so electrostatic forces are expected to play an important role in determining elastic properties of these nanotubes. The authors generate the coordinates of nanotubes of different chirality’s and size. Each and every structure thus generated is allowed to relax till the authors obtain minima of energy. The authors then apply the requisite compressions, elongations and twists to the structures and compute the elastic moduli. Young’s Modulus, Shear Modulus and Poisson’s ratio for single-walled armchair and zigzag tubes of different chirality’s and size have been calculated. The computational results show the variation of Young’s Modulus, Poisson’s ratio and Shear Modulus for these NNTs with nanotube diameter. The results have been compared with available data, experimental as well as theoretical. Findings – The authors have calculated bond length, cohesive energy/bond, Strain energy, Young’s Modulus, Shear Modulus and Poisson’s ratio. Originality/value – To the best of the knowledge this work is the first attempt to study elastic properties of III-V NNTs using second generation REBO potential

Study of the Nanomechanics of CNTs under Tension by Molecular Dynamics Simulation Using Different Potentials

At four different strain rates, the tensile stress strain relationship of single-walled 12-12 CNT with aspect ratio 9.1 obtained by Rebo potential (Brenner, 1990), Airebo potential (Stuart et al., 2000), and Tersoff potential (Tersoff, 1988) is compared with that of Belytschko et al. (2002) to validate the present model. Five different empirical potentials such as Rebo potential (Brenner, 1990), Rebo potential (Brenner et al., 2002), Inclusion LJ with Rebo potential (Brenner, 1990), Airebo potential (Stuart et al., 2000), and Tersoff potential (Tersoff, 1988) are used to simulate CNT subjected to axial tension differing its geometry at high strain rate. In Rebo potential (Mashreghi and Moshksar, 2010) only bond-order term is used and in Rebo potential (Brenner et al., 2002) torsional term is included with the bond-order term. At high strain rate the obtained stress strain relationships of CNTs subjected to axial tension differing its geometries using five different potentials are compared with the published results and from the comparison of the results, the drawback of the published results and limitations of different potentials are evaluated and the appropriate potential is selected which is the best among all other potentials to study the elastic, elastic-plastic properties of different types of CNTs. The present study will help a new direction to get reliable elastic, elastic-plastic properties of CNTs at different strain rates. Effects of long range Van der Waals interaction and torsion affect the elastic, elastic-plastic properties of CNTs and why these two effects are really needed to consider in bond-order Rebo potential (Brenner, 1990) to get reliable elastic, elastic-plastic properties of CNTs is also discussed. Effects of length-to-diameter ratio, layering of CNTs, and different empirical potentials on the elastic, elastic-plastic properties of CNTs are discussed in graphical and tabular forms with published results as a comparative manner to understand the nanomechanics of CNTs under tension using molecular dynamics simulation.

Molecular Dynamics Simulation to Study the Effect of Empirical Potential Functions on Modeling of Diamond-Like Carbon

Molecular dynamics (MD) simulations is an effective method to investigate the mechanical and tribological properties of amorphous carbon since the coordinates of all atoms can be calculated as a function of time. Several empirical potentials can be used to model the interatomic interactions of carbon atoms, including the Tersoff potential, the Reactive Bond Order (REBO) potential and its revised versions, and the Reactive Force Field (ReaxFF) potential. The choice of empirical potential is one of the fundamental and important assumptions in the MD approach since it can affect the properties of amorphous carbon during the MD simulations. In this study, liquid quenching method is used to model amorphous carbon for computational efficiency. We will study the influence of the three types of potentials, specifically the Tersoff potential, the 2nd REBO potential and the ReaxFF potential on DLC parameters. These parameters include the sp3 content as a function of density, the arrangement of the amorphous carbon atoms, hybridization and the radial distribution functions G(r).