Benjamin D. Jensen

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.

Simulation of mechanical performance limits and failure of carbon nanotube composites

The mechanical properties of carbon nanotube (CNT) fiber composites are steadily approaching those of traditional carbon fiber composites. This work is focused on establishing a plausible upper bound on these properties by modeling the elastic deformations, yield, and fracture of idealized CNT composites using reactive molecular dynamics. Amorphous carbon (AC) was used for the matrix material because of its structural simplicity and physical compatibility with the CNT fillers. Three different arrangements of CNTs in the simulation cell were investigated: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. The SWNT and MWNT array systems are clearly idealizations, but the SWNT bundle system is a step closer to real systems in which individual tubes aggregate into large assemblies. Chemical crosslinking was modeled by adding bonds between the CNTs and AC to explore the balance between weakening the CNTs and improving fiber-matrix load transfer. The simulation results reported here clarify the impact of CNT dispersion, the extent of crosslinking, and CNT-templated matrix structuring on the mechanical properties of CNT composites.

Atomistic Modeling of Cross-linked Epoxy Polymer

Molecular Dynamics simulations are used to study cross-linking of an epoxy polymer. OPLS force field parameters are used for modeling a 2:1 stoichiometric mixture of epoxy resin and the cross-linking agent. The model has 17,928 united atoms and a static cross-linking method is used along with molecular minimization and molecular dynamics techniques to achieve two different cross-link densities. The crosslinked models can be used for understanding various phenomenon occurring in cross-linked epoxy resins at the atomic scale. Glass-transition temperature ranges of two differently cross-linked samples have been predicted using the models. These models will be used for studying aging behavior at the atomic level in epoxy materials and understanding the influence of aging on mechanical properties.