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.

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.

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.

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.

Evolution of Glassy Carbon Derived from Pyrolysis of Furan Resin.

Glassy carbon (GC) material derived from pyrolyzed furan resin was modeled by using reactive molecular dynamics (MD) simulations. The MD polymerization simulation protocols to cure the furan resin precursor material are validated via comparison of the predicted density and Young's modulus with experimental values. The MD pyrolysis simulations protocols to pyrolyze the furan resin precursor is validated by comparison of calculated density, Young's modulus, carbon content, sp2 carbon content, the in-plane crystallite size, out-of-plane crystallite stacking height, and interplanar crystallite spacing with experimental results from the literature for furan resin derived GC. The modeling methodology established in this work can provide a powerful tool for the modeling-driven design of next-generation carbon-carbon composite precursor chemistries for thermal protection systems and other high-temperature applications.

In situ mechanical property measurements of amorphous carbon-boron nitride nanotube nanostructures.

To understand the mechanical properties of amorphous carbon (a-C)/boron nitride nanotube (BNNT) nanostructures, in situ mechanical tests are conducted inside a transmission electron microscope equipped with an integrated atomic force microscope system. The nanotube structure is modified with amorphous carbon deposited by controlled electron beam irradiation. We demonstrate multiple in situ tensile, compressive, and lap shear tests with a-C/BNNT hybrid nanostructures. The tensile strength of the a-C/BNNT hybrid nanostructure is 5.29 GPa with about 90 vol% of a-C. The tensile strength and strain of the end-to-end joint structure with a-C welding is 0.8 GPa and 5.2% whereas the lap shear strength of the side-by-side joint structure with a-C is 0.25 GPa.

Voxelized Atomic Structure Potentials: Predicting Atomic Forces with the Accuracy of Quantum Mechanics Using Convolutional Neural Networks.

This paper introduces voxelized atomic structure (VASt) potentials as a machine learning (ML) framework for developing interatomic potentials. The VASt framework utilizes a voxelized representation of the atomic structure directly as the input to a convolutional neural network (CNN). This allows for high-fidelity representations of highly complex and diverse spatial arrangements of the atomic environments of interest. The CNN implicitly establishes the low-dimensional features needed to correlate each atomic neighborhood to its net atomic force. The selection of the salient features of the atomic structure (i.e., feature engineering) in the VASt framework is implicit, comprehensive, automated, scalable, and highly efficient. The calibrated convolutional layers learn the complex spatial relationships and multibody interactions that govern the physics of atomic systems with remarkable fidelity. We show that VASt potentials predict highly accurate forces on two phases of silicon carbide and the thermal conductivity of silicon over a range of isotropic strain.

Multifunctional Electroactive Nanocomposites Based on Piezoelectric Boron Nitride Nanotubes.

Space exploration missions require sensors and devices capable of stable operation in harsh environments such as those that include high thermal fluctuation, atomic oxygen, and high-energy ionizing radiation. However, conventional or state-of-the-art electroactive materials like lead zirconate titanate, poly(vinylidene fluoride), and carbon nanotube (CNT)-doped polyimides have limitations on use in those extreme applications. Theoretical studies have shown that boron nitride nanotubes (BNNTs) have strength-to-weight ratios comparable to those of CNTs, excellent high-temperature stability (to 800 °C in air), large electroactive characteristics, and excellent neutron radiation shielding capability. In this study, we demonstrated the experimental electroactive characteristics of BNNTs in novel multifunctional electroactive nanocomposites. Upon application of an external electric field, the 2 wt % BNNT/polyimide composite was found to exhibit electroactive strain composed of a superposition of linear piezoelectric and nonlinear electrostrictive components. When the BNNTs were aligned by stretching the 2 wt % BNNT/polyimide composite, electroactive characteristics increased by about 460% compared to the nonstretched sample. An all-nanotube actuator consisting of a BNNT buckypaper layer between two single-walled carbon nanotube buckypaper electrode layers was found to have much larger electroactive properties. The additional neutron radiation shielding properties and ultraviolet/visible/near-infrared optical properties of the BNNT composites make them excellent candidates for use in the extreme environments of space missions.

Toward high performance thermoset/carbon nanotube sheet nanocomposites via resistive heating assisted infiltration and cure.

Thermoset/carbon nanotube (CNT) sheet nanocomposites were successfully fabricated by resistive heating assisted infiltration and cure (RHAIC) of the polymer matrix resin. Resistive heating takes advantage of the electrical and thermal conductivity of CNTs to rapidly and uniformly introduce heat into the CNT sheet. Heating the CNT sheet reduces the viscosity of the polymer resin due to localized temperature rise in close proximity to the resin, which enhances resin flow, penetration, and wetting of the CNT reinforcement. Once the resin infusion process is complete, the applied power is increased to raise the temperature of the CNT sheet, which rapidly cures the polymer matrix. Tensile tests were used to evaluate the mechanical properties of the processed thermoset/CNT sheet nanocomposites. The improved wetting and adhesion of the polymer resin to the CNT reinforcement yield significant improvement of thermoset/CNT nanocomposite mechanical properties. The highest specific tensile strength of bismaleimide(BMI)/CNT sheet nanocomposites was obtained to date was 684 MPa/(g/cm(3)), using 4 V (2 A) for resin infiltration, followed by precure at 10 V (6 A) for 10 min and post curing at 240 °C for 6 h in an oven. The highest specific Young's modulus of BMI/CNT sheet nanocomposite was 71 GPa/(g/cm(3)) using resistive heating infiltration at 8.3 V (4.7 A) for 3 min followed by resistive heating cure at 12.5 V (7 A) for 30 min. In both cases, the CNT sheets were stretched and held in tension to prevent relaxation of the aligned CNTs during the course of RHAIC.

Polyaniline/carbon nanotube sheet nanocomposites: fabrication and characterization.

Practical approaches are needed to take advantage of the nanometer-scale mechanical properties of carbon nanotubes (CNTs) at the macroscopic scale. This study was conducted to elucidate the salient factors that can maximize the mechanical properties of nanocomposites fabricated from commercially available CNT sheets. The CNT sheets were modified by stretching to improve CNT alignment and in situ polymerization using polyaniline (PANI), a π-conjugated conductive polymer, as a binder. The resulting CNT nanocomposites were subsequently postprocessed by hot pressing and/or high temperature treatment to carbonize the PANI as a means to improve mechanical properties. The PANI/CNT nanocomposites demonstrated significant improvement in mechanical properties compared to pristine CNT sheets. The highest specific tensile strength of PANI/stretched CNT nanocomposite was 484 MPa/(g/cm3), which was achieved in a sample with ∼42 wt % of PANI. This specimen was fabricated by in situ polymerization followed by hot pressing. The highest specific Young's modulus of 17.1 GPa/(g/cm3) was measured on a sample that was hot-pressed and carbonized. In addition, the highest DC-electrical conductivity of 621 S/cm was obtained on a sample prepared by in situ polymerization of PANI on a stretched CNT sheet.

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.

Ultrafast dynamics of single-walled carbon nanotubes dispersed in polymer films.

Wavelength-resolved femtosecond transient absorption spectroscopy is used to follow the electronic dynamics of single-walled carbon nanotubes in polymers following visible photoexcitation. Electron-hole (e-h) pairs give rise to sharp features in the transient spectra that decay in amplitude and exhibit rapid spectral shifts. The decay reflects (e-h) recombination on both short (1.3 ps) and long (35 ps) time scales. Transient spectra also exhibit a broad photobleach at early times that arises from the cooling of a hot electron gas created via excitation at the red edge of a pi-plasmon band.