Molecular Dynamics Modeling of Interfacial Interactions between Flattened Carbon Nanotubes and Amorphous Carbon: Implications for Ultra-Lightweight Composites.

Flattened carbon nanotubes (flCNTs) naturally form in many carbon nanotube-based materials and can exhibit mechanical properties similar to round carbon nanotubes but with tighter packing and alignment. To facilitate the design, fabrication, and testing of flCNT-based composites for aerospace structures, computational modeling can be used to efficiently and accurately predict their performance as a function of processing parameters, such as reinforcement/matrix cross-linking. In this study, molecular dynamics modeling is used to predict the load transfer characteristics of the interface region between the flat region of flCNTs (i.e., bi-layer graphene) and amorphous carbon (AC) with various levels and locations of covalent bond cross-linking and AC mass density. The results of this study show that increasing the mass density of AC at the interface improves the load transfer capability of the interface. However, a much larger improvement is observed when cross-linking is added both to the flCNT-AC interface and between the flCNT sheets. With both types of cross-linking, substantial improvements in interfacial shear strength, transverse tension strength, and transverse tension toughness are predicted. The results of this study are important for optimizing the processing of flCNT/AC composites for demanding engineering applications.

Simulating the effects of carbon nanotube continuity and interfacial bonding on composite strength and stiffness.

Molecular dynamics simulations of carbon nanotube (CNT) composites, in which the CNTs are continuous across the periodic boundary, overestimate the experimentally measured mechanical properties of CNT composites along the fiber direction. Since the CNTs in these composites are much shorter than the composite dimensions, load must be transferred either directly between CNTs or through the matrix, a mechanism that is absent in simulations of effectively continuous CNTs. In this study, the elastic and fracture properties of high volume fraction discontinuous carbon nanotube/amorphous carbon composite systems were compared to those of otherwise equivalent continuous CNT composites using ReaxFF reactive molecular dynamics simulations. These simulations were used to show how the number of nanotube-matrix interfacial covalent bonds affect composite mechanical properties. Furthermore, the mechanical impact of interfacial bonding was decomposed to reveal its effect on the properties of the CNTs, the interfacial layer of matrix, and the bulk matrix. For the composites with continuous reinforcement, it was found that any degree of interfacial bonding has a negative impact on axial tensile strength and stiffness. This is due to disruption of the structure of the CNTs and interfacial matrix layer by the interfacial bonds. For the discontinuous composites, the modulus was maximized between 4%-7% interfacial bonding and the strength continues to increase up to the highest levels of interfacial bonding studied. Areas of low stress and voids were observed in the simulated discontinuous composites at the ends of the tubes, from which fracture was observed to initiate. Experimental carbon nanotube yarn composites were fabricated and tested. The results were used to illustrate knockdown factors relative to the mechanical performance of the tubes themselves.

Comparing the Mechanical Response of Di-, Tri-, and Tetra-functional Resin Epoxies with Reactive Molecular Dynamics.

The influence of monomer functionality on the mechanical properties of epoxies is studied using Molecular Dynamics (MD) with the Reax Force Field (ReaxFF). From deformation simulations, the Young's modulus, yield point, and Poisson's ratio are calculated and analyzed. The results demonstrate an increase in stiffness and yield strength with increasing resin functionality. Comparison between the network structures of distinct epoxies is further advanced by the Monomeric Degree Index (MDI). Experimental validation demonstrates the MD results correctly predict the relationship in Young's moduli. Therefore, ReaxFF is confirmed to be a useful tool for studying the mechanical behavior of epoxies.

Chemical Reactions in Classical Molecular Dynamics.

An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as fix bond/react, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candidate reactants are first identified by interatomic separation, followed by the application of a generalized topology matching algorithm to confirm they match the pre-reaction template. This is followed by a topology conversion to match the post-reaction template and a dynamic relaxation to minimize high energy configurations. Two case studies, the condensation polymerization of nylon 6,6 and the formation of a highly-crosslinked epoxy, are simulated to demonstrate the robustness, stability, and speed of the algorithm. Improvements which could increase its utility are discussed.

Simulation of the Elastic and Ultimate Tensile Properties of Diamond, Graphene, Carbon Nanotubes, and Amorphous Carbon Using a Revised ReaxFF Parametrization.

In light of the enduring interest in using nanostructured carbon materials as reinforcing elements in composite materials, there is a significant need for a reliable computational tool capable to predict the mechanical properties, both elastic properties and properties at the point of fracture, in large-scale atomistic simulations. A revised version of the ReaxFF reactive force field parametrization for carbon, ReaxFFC-2013, was recently published and is notable because of the inclusion of density functional theory (DFT)-derived mechanical data for diamond and graphite in the fitting set. The purpose of the present work is to assess the accuracy of this new force field for predicting the mechanical properties for several allotropes of carbon, both in the elastic regime and during fracture. The initial discussion focuses on the performance of ReaxFFC-2013 for diamond and graphene, the two carbon forms for which mechanical properties were included in the parametrization data set. After it is established that simulations conducted with the new force field yield results that agree well with DFT and experimental data for most properties of interest, its transferability to amorphous carbon and carbon nanotubes is explored. ReaxFFC-2013 is found to produce results that, for the most part, compare favorably with available experimental data for single and multiwalled nanotubes and for amorphous carbon models prepared over a range of densities. Although there is opportunity for improvement in some predicted properties, the ReaxFFC-2013 parametrization is shown to generally perform well for each form of carbon and to compare favorably with DFT and experimental data.

The effect of time step, thermostat, and strain rate on ReaxFF simulations of mechanical failure in diamond, graphene, and carbon nanotube.

As the sophistication of reactive force fields for molecular modeling continues to increase, their use and applicability has also expanded, sometimes beyond the scope of their original development. Reax Force Field (ReaxFF), for example, was originally developed to model chemical reactions, but is a promising candidate for modeling fracture because of its ability to treat covalent bond cleavage. Performing reliable simulations of a complex process like fracture, however, requires an understanding of the effects that various modeling parameters have on the behavior of the system. This work assesses the effects of time step size, thermostat algorithm and coupling coefficient, and strain rate on the fracture behavior of three carbon-based materials: graphene, diamond, and a carbon nanotube. It is determined that the simulated stress-strain behavior is relatively independent of the thermostat algorithm, so long as coupling coefficients are kept above a certain threshold. Likewise, the stress-strain response of the materials was also independent of the strain rate, if it is kept below a maximum strain rate. Finally, the mechanical properties of the materials predicted by the Chenoweth C/H/O parameterization for ReaxFF are compared with literature values. Some deficiencies in the Chenoweth C/H/O parameterization for predicting mechanical properties of carbon materials are observed.

Parametric Study of ReaxFF Simulation Parameters for Molecular Dynamics Modeling of Reactive Carbon Gases.

The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although the Reax Force Field (ReaxFF) can be used to simulate the chemical behavior of carbon-based systems, the simulation settings required for accurate predictions have not been fully explored. Using the ReaxFF, molecular dynamics (MD) simulations are used to simulate the chemical behavior of pure carbon and hydrocarbon reactive gases that are involved in the formation of carbon structures such as graphite, buckyballs, amorphous carbon, and carbon nanotubes. It is determined that the maximum simulation time step that can be used in MD simulations with the ReaxFF is dependent on the simulated temperature and selected parameter set, as are the predicted reaction rates. It is also determined that different carbon-based reactive gases react at different rates, and that the predicted equilibrium structures are generally the same for the different ReaxFF parameter sets, except in the case of the predicted formation of large graphitic structures with the Chenoweth parameter set under specific conditions.