Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The origin of passivity in aluminum-manganese solid solutions

In this work, the effects of manganese on the aqueous corrosion of aluminum-manganese alloys were investigated by experiments and atomistic simulations. Electrochemical measurements, x-ray photoelectron spectroscopy, and atom probe tomography analysis indicate that manganese addition enhanced the corrosion resistance of aluminum without participating in the surface oxidation. Furthermore, the selective dissolution of manganese was believed to increase the free volume at the metal/oxide interface to facilitate the formation of a denser, thinner oxide layer with closer to stoichiometry composition. Atomistic simulations showed that the oxide layer compactness increased, and defect density decreased with increasing free volume content in Al, resulting in enhanced barrier characteristics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

ORNL_AISD_NiPt_108atoms

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for crystal structures with 108 atoms with 1,900 configurations. The data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system 'MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001' [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in a compressed zipped folders atoms108.zip. The zipped folder contains the data structured in the following way: - Ni_ground_state.cfg --> atomic configuration for the pure nickel - Pt_ground_state.cfg --> atomic configuration for the pure platinum - Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder - Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms108.zip is 66,132. This dataset is an extension to the dataset ORNL_AISD_NiPt [5] that has been previously released with crystal structures of 256 atoms, 864 atoms, and 2,048 atoms, with the same methodology for data collection. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D [5] M. Karabin, M. Lupo Pasini, and M. Eisenbach. ORNL_AISD_NiPt. United States: N. p., 2023. Web. doi:10.13139/OLCF/1958172.

36 MATERIALS SCIENCE↗
DOE OSTI2023

ORNL_AISD_NiPt

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for three different sizes of the crystal structure: 256 atoms, 864 atoms, and 2,048 atoms, each of which contains 1900 configurations. For each size of the crystal structure, the data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in three compressed zipped folders: atoms256.zip, atoms864.zip, atoms2048.zip Each zipped folder contains the data that describes crystals of size 256 atoms, 864 atoms, and 2,048 atoms respectively. Each one of the three zipped folders contains the data structured in the following way: -Ni_ground_state.cfg --> atomic configuration for the pure nickel -Pt_ground_state.cfg --> atomic configuration for the pure platinum -Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder. Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms256.zip is 65,046. The total number of atomic configurations contained in atoms864.zip is 63,936. The total number of atomic configurations contained in atoms2048.zip is 61,997. The total number of atomic configurations spanned by the entire dataset is 190,979. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büccchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D

36 MATERIALS SCIENCE↗
DOE OSTI2024

Coordinate-based simulation of pair distance distribution functions for small and large molecular assemblies: implementation and applications

X-ray scattering has become a major tool in the structural characterization of nanoscale materials. Thanks to the widely available experimental and computational atomic models, coordinate-based X-ray scattering simulation has played a crucial role in data interpretation in the past two decades. However, simulation of real-space pair distance distribution functions (PDDFs) from small- and wide-angle X-ray scattering, SAXS/WAXS, has been relatively less exploited. This study presents a comparison of PDDF simulation methods, which are applied to molecular structures that range in size from β-cyclodextrin [1 kDa molecular weight (MW), 66 non-hydrogen atoms] to the satellite tobacco mosaic virus capsid (1.1 MDa MW, 81 960 non-hydrogen atoms). The results demonstrate the power of interpretation of experimental SAXS/WAXS from the real-space view, particularly by providing a more intuitive method for understanding of partial structure contributions. Furthermore, the computational efficiency of PDDF simulation algorithms makes them attractive as approaches for the analysis of large nanoscale materials and biological assemblies. The simulation methods demonstrated in this article have been implemented in stand-alone software, SolX 3.0 , which is available to download from https://12idb.xray.aps.anl.gov/solx.html.

Zuo, Xiaobing (ORCID:0000000201344804)↗
DOE OSTI2020

Polarizabilities of neutral atoms and atomic ions with a noble gas electron configuration

Atomic polarizabilities play an important role in the development of force fields for molecular simulations, as well as for the development of qualitative concepts of atomic and molecular behavior. Coupled cluster theory at the coupled cluster singles doubles triples level with very large correlation-consistent basis sets with extended diffuse functions has been used to predict the polarizabilities of the atomic neutrals, mono-cations and mono-anions with a noble gas configuration. Additional corrections for scalar relativistic and spin–orbit effects were also included for the electron configurations of Kr, Xe, and Rn. The results are in excellent agreement with experiment or with other high level calculations where available. The current results for most of these species represent the best available values for the polarizabilities. The results show that the polarizability of H - is very difficult to calculate without extremely diffuse functions. The polarizability of H - is the largest value, 34.05 Å 3 , calculated for all species in the current study. The polarizabilities of the remaining halogen anions are also the best available values. The polarizabilities of the halogen anions (excluding F - ) and H - have a linear correlation with the electron affinity of the neutral atom. Spin– orbit effects, even for closed shell species, cannot be ignored for quantitative accuracy, and the inclusion of spin–orbit effects for Fr + , Rn, and At - increases the polarizability by 4%, 6%, and 15%, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

$\mathrm{SIM}$-$\mathrm{STEM}$ $\mathrm{Lab}$: Incorporating Compressed Sensing Theory for Fast STEM Simulation

We report recently it has been shown that precise dose control and an increase in the overall acquisition speed of atomic resolution scanning transmission electron microscope (STEM) images can be achieved by acquiring only a small fraction of the pixels in the image experimentally and then reconstructing the full image using an inpainting algorithm. In this paper, we apply the same inpainting approach (a form of compressed sensing) to simulated, sub-sampled atomic resolution STEM images. We find that it is possible to significantly sub-sample the area that is simulated, the number of g-vectors contributing the image, and the number of frozen phonon configurations contributing to the final image while still producing an acceptable fit to a fully sampled simulation. Here we discuss the parameters that we use and how the resulting simulations can be quantifiably compared to the full simulations. As with any Compressed Sensing methodology, care must be taken to ensure that isolated events are not excluded from the process, but the observed increase in simulation speed provides significant opportunities for real time simulations, image classification and analytics to be performed as a supplement to experiments on a microscope to be developed in the future

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Interpretable Machine Learning Models for Practical Antimonate Electrocatalyst Performance

Computationally predicting the performance of catalysts under reaction conditions is a challenging task due to the complexity of catalytic surfaces and their evolution in situ, different reaction paths, and the presence of solid-liquid interfaces in the case of electrochemistry. We demonstrate here how relatively simple machine learning models can be found that enable prediction of experimentally observed onset potentials. Inputs to our model are comprised of data from the oxygen reduction reaction on non-precious transition-metal antimony oxide nanoparticulate catalysts with a combination of experimental conditions and computationally affordable bulk atomic and electronic structural descriptors from density functional theory simulations. From human-interpretable genetic programming models, we identify key experimental descriptors and key supplemental bulk electronic and atomic structural descriptors that govern trends in onset potentials for these oxides and deduce how these descriptors should be tuned to increase onset potentials. Here, we finally validate these machine learning predictions by experimentally confirming that scandium as a dopant in nickel antimony oxide leads to a desired onset potential increase. Macroscopic experimental factors are found to be crucially important descriptors to be considered for models of catalytic performance, highlighting the important role machine learning can play here even in the presence of small datasets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Perspective: Are glasses really frozen?

A glass is obtained by cooling a liquid fast enough to avoid crystallization. It is a solid with elastic constants and density similar to those of a crystal with the same composition. It is customarily assumed that the structure of a glass is frozen, and atoms are largely locked into a random, but specific, structure; while a small number of loosely bound atoms may be mobile, most atoms are frozen at room temperature and below. However, recent simulation studies on metallic glasses present a rather different picture. A very significant fraction (~20%) of atoms are not frozen and locally mobile, even down to T = 0 partly because of quantum effect. Here, we describe these observations and discuss their implications on the nature of the glass structure and the origin of the glass transition.

36 MATERIALS SCIENCE↗
NASA NTRS2002

Test and Analysis Capabilities of the Space Environment Effects Team at Marshall Space Flight Center

Marshall Space Flight Center has developed world-class space environmental effects testing facilities to simulate the space environment. The combined environmental effects test system exposes temperature-controlled samples to simultaneous protons, high- and low-energy electrons, vacuum ultraviolet (VUV) radiation, and near-ultraviolet (NUV) radiation. Separate chambers for studying the effects of NUV and VUV at elevated temperatures are also available. The Atomic Oxygen Beam Facility exposes samples to atomic oxygen of 5 eV energy to simulate low-Earth orbit (LEO). The LEO space plasma simulators are used to study current collection to biased spacecraft surfaces, arcing from insulators and electrical conductivity of materials. Plasma propulsion techniques are analyzed using the Marshall magnetic mirror system. The micro light gas gun simulates micrometeoroid and space debris impacts. Candidate materials and hardware for spacecraft can be evaluated for durability in the space environment with a variety of analytical techniques. Mass, solar absorptance, infrared emittance, transmission, reflectance, bidirectional reflectance distribution function, and surface morphology characterization can be performed. The data from the space environmental effects testing facilities, combined with analytical results from flight experiments, enable the Environmental Effects Group to determine optimum materials for use on spacecraft.

Finckenor, M. M.↗
NASA NTRS1983

Computer simulation of surface and film processes

Adequate computer methods, based on interactions between discrete particles, provide information leading to an atomic level understanding of various physical processes. The success of these simulation methods, however, is related to the accuracy of the potential energy function representing the interactions among the particles. The development of a potential energy function for crystalline SiO2 forms that can be employed in lengthy computer modelling procedures was investigated. In many of the simulation methods which deal with discrete particles, semiempirical two body potentials were employed to analyze energy and structure related properties of the system. Many body interactions are required for a proper representation of the total energy for many systems. Many body interactions for simulations based on discrete particles are discussed.

Tiller, W. A.↗
DOE OSTI2025

Accelerated Fuel Qualification of Fast Modular Reactor Fuel in a Thermal Reactor: Modeling and Simulation Paired with Irradiation Testing

The accelerated fuel qualification (AFQ) methodology is applied by simulating accelerated fuel tests of the General Atomics Electromagnetic Systems’ fuel system for its 44-MW(electric) gas-cooled, fast-spectrum fast modular reactor (FMR). This fuel is comprised of UO 2 pellets in SiGA® cladding, a silicon carbide ceramic matrix composite. Fast reactors, like the FMR, offer many benefits, including high fuel utilization and flexibility, but may require a lengthy material design process if tests are performed using fast neutron irradiation alone. A thermal neutron irradiation can instead be used to rapidly test how well key components of the current material models extend to high burnup. Thermal neutrons produce a different radial power distribution within the pin than fast neutrons. However, the temperature and burnup values for the two neutron types are comparable, and the differences between the simulated fuel responses are relatively small, demonstrating the weak sensitivity of the physics-based fuel model calculations on the neutron type and the irradiation rate. Furthermore, the deformation of the SiGA cladding saturates after about 1 displacement per atom for both neutron spectra. In an accelerated fuel test, the irradiation time required to reach the target fuel burnup can be reduced by a factor of 3 by using a small rodlet with a 45% smaller pellet diameter while maintaining the same linear power. Therefore, the time for data collection up to high burnup can be significantly reduced while maintaining the same temperature profile, which largely determines the material response. Tests of fuel rodlets of standard and compact size will be carried out in the Idaho National Laboratory’s Advanced Test Reactor (ATR), including full size and compact rodlets with varying gap sizes. By applying physics-based mechanistic modeling and simulation in accordance with the AFQ methodology, this type of compact rodlet testing in a thermal test reactor captures the necessary phenomena to test fuel material models up to high burnup and to simulate the expected impact of fast neutron radiation on the fuel in FMR operations. Finally, this approach to testing fast reactor fuels in existing thermal test reactors, paired with advanced physics-based mechanistic modeling and simulation, is expected to be applicable to a range of advanced fuels and will decrease the overall fuel qualification timeframe from decades to years.

Advanced test reactor (ATR)↗
DOE OSTI2022

Relation between interfacial shear and friction force in 2D materials

Understanding the interfacial properties between an atomic layer and its substrate is of key interest at both the fundamental and technological levels. From Fermi level pinning to strain engineering and superlubricity, the interaction between a single atomic layer and its substrate governs electronic, mechanical and chemical properties. Here, we measure the hardly accessible interfacial transverse shear modulus of an atomic layer on a substrate. By performing measurements on bulk graphite, and on epitaxial graphene films on SiC with different stacking orders and twisting, as well as in the presence of intercalated hydrogen, we find that the interfacial transverse shear modulus is critically controlled by the stacking order and the atomic layer–substrate interaction. Importantly, we demonstrate that this modulus is a pivotal measurable property to control and predict sliding friction in supported two-dimensional materials. The experiments demonstrate a reciprocal relationship between friction force per unit contact area and interfacial shear modulus. As a result, the same relationship emerges from simulations with simple friction models, where the atomic layer– substrate interaction controls the shear stiffness and therefore the resulting friction dissipation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A nearsighted force-training approach to systematically generate training data for the machine learning of large atomic structures

A challenge of atomistic machine-learning (ML) methods is ensuring that the training data are suitable for the system being simulated, which is particularly challenging for systems with large numbers of atoms. Most atomistic ML approaches rely on the nearsightedness principle (“all chemistry is local”), using information about the position of an atom’s neighbors to predict a per-atom energy. Here, in this work, we develop a framework that exploits the nearsighted nature of ML models to systematically produce an appropriate training set for large structures. We use a per-atom uncertainty estimate to identify the most uncertain atoms and extract chunks centered around these atoms. It is crucial that these small chunks are both large enough to satisfy the ML’s nearsighted principle (that is, filling the cutoff radius) and are large enough to be converged with respect to the electronic structure calculation. We present data indicating when the electronic structure calculations are converged with respect to the structure size, which fundamentally limits the accuracy of any nearsighted ML calculator. These new atomic chunks are calculated in electronic structures, and crucially, only a single force—that of the central atom—is added to the growing training set, preventing the noisy and irrelevant information from the piece’s boundary from interfering with ML training. The resulting ML potentials are robust, despite requiring single-point calculations on only small reference structures and never seeing large training structures. We demonstrated our approach via structure optimization of a 260-atom structure and extended the approach to clusters with up to 1415 atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1979

Hot-atom synthesis of organic compounds on Jupiter

Results of recent laboratory 'simulations' of photochemical processes on Jupiter are combined with available data on mixing rates and exposure times in the Jovian atmosphere to give quantitative predictions of the rate at which hot-atom reactions produce organic molecules. It is shown that abstraction reactions on methane by hot H atoms from solar UV photolysis of H2S will produce no more than 4 times 10 to the -17th power g/sq cm/sec for a steady-state mole fraction of total organics of approximately 10 to the -16th power. This is roughly 10 to the 7th power times less than the limit of detection of the most sensitive gas analysis experiments ever flown on a spacecraft. By far the most common organic molecule produced by this mechanism is CH3SH, methyl mercaptan, which is produced at a rate at least 600 times smaller than the rate of production of ethane by direct photolysis of CH4 at high altitudes.

Lewis, J. S.↗
DOE OSTI2022

Predicted Melt Curve and Liquid Shear Viscosity of RDX up to 30 GPa

Recent grain-scale simulations of HMX and TATB have shown that predictions for hot spot formation in high explosives are particularly sensitive to accurate determinations of the pressure-dependent melt curve and the shear viscosity of the liquid phase. These physics terms are poorly constrained beyond ambient pressure for the explosive RDX. Here we adopt an all-atom modeling approach using molecular dynamics (MD) simulations to predict the melt curve of RDX near to detonation conditions (30 GPa) and determine the shear viscosity of the liquid as a function of temperature and pressure above the melt curve. Phase-coexistence simulations were used to determine the melt curve, which is predicted to vary by almost 1100 K as the pressure increases from 0 GPa to 30 GPa. Equilibrium MD simulations and the Green-Kubo formalism were used to obtain the pressure-temperature-dependent shear viscosity. The shear viscosity of RDX is predicted to be of similar magnitude to the viscosity of TATB at low GPa-range pressures, and to be roughly an order of magnitude lower than the viscosity of HMX. The temperature dependence of the shear viscosity is Arrhenius at a given pressure, and the exponential pre-factor and activation term exhibit a strong, yet complicated, pressure dependence. An empirical pressure-temperature dependent function for RDX shear viscosity is developed that simultaneously captures a wide range of MD predictions while taking an analytic form that extrapolates smoothly beyond the fitted regime. The relative strength of the pressure and temperature dependencies of these two physics terms is found to be of similar magnitude for RDX, HMX, and TATB, which motivates incorporating these results in future RDX grain-scale modeling.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Evaporation kinetics of Mg2SiO4 crystals and melts from molecular dynamics simulations

Computer simulations based on the molecular dynamics (MD) technique were used to study the mechanisms and kinetics of free evaporation from crystalline and molten forsterite (i.e., Mg2SiO4) on an atomic level. The interatomic potential employed for these simulations reproduces the energetics of bonding in forsterite and in gas-phase MgO and SiO2 reasonably accurately. Results of the simulation include predicted evaporation rates, diffusion rates, and reaction mechanisms for Mg2SiO4(s or l) yields 2Mg(g) + 20(g) + SiO2(g).

Kubicki, J. D.↗
DOE OSTI2025

The importance of accounting for the Tolman correction to surface tension for nucleation and growth modeling of Fe clusters

Gibbs free energies of clusters are required for predictive modeling of cluster growth during condensation of a cooling vapor. Here, we present a straightforward method of calculating free energies of cluster formation using the data from molecular dynamics (MD) simulations. We apply this method to iron clusters having from 2 to 100 atoms. The energies obtained are verified by comparing to an MD-simulated equilibrium cluster size distribution in a sub-saturated vapor. We show that these free energies differ significantly from those obtained with a commonly used spherical cluster approximation, which relies on a surface tension coefficient of a flat surface, as it is used in the classical nucleation theory (CNT). We show that the spherical cluster approximation in CNT can be improved by using a cluster-size-dependent Tolman correction for the surface tension. The Tolman length and effective surface tension values were derived for iron clusters, and they significantly differ from the commonly used experimentally measured values. This improved approximation does not account for geometric magic number effects responsible for spikes and troughs in densities of neighbor cluster sizes. Nonetheless, it allows to more accurately model cluster formation from a cooling vapor. It better reproduces the condensation timeline, overall shape of the cluster size distribution, average cluster size, and the distribution width. In contrast, using a constant surface tension coefficient (as done in CNT) resulted in incorrect condensation dynamics and cluster size distributions. The analytical expression for cluster nucleation rate from CNT was updated to account for the size-dependence of cluster surface tension.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

Molecular Dynamics Modeling of Carbon Nanotube Composite Fracture Using ReaxFF

Carbon nanotube (CNT) fiber reinforced composites with specific tensile strengths and moduli approaching those of aerospace grade carbon fiber composites have recently been reported. This achievement was enabled by the emerging availability of high N/tex yarns in kilometer-scale quantities. While the production of this yarn is an impressive advance, its strength is still much lower than that of the individual CNTs comprising the yarn. Closing this gap requires understanding load transfer between CNTs at the nanometer dimensional scale. This work uses reactive molecular dynamics simulations to gain an understanding at the nanometer scale of the key factors that determine CNT nanocomposite mechanical performance, and to place more realistic upper bounds on the target properties. While molecular dynamics simulations using conventional force fields can predict elastic properties, the ReaxFF reactive forcefield can also model fracture behavior because of its ability to accurately describe bond breaking and formation during a simulation. The upper and lower bounds of CNT composite properties are investigated by comparing systems composed of CNTs continuously connected across the periodic boundary with systems composed of finite length CNTs. These lengths, effectively infinite for the continuous tubes and an aspect ratio of 13 for the finite length case, result from practical limitations on the number of atoms that can be included in a simulation. Experimentally measured aspect ratios are typically on the order of 100,000, so the calculated results should represent upper and lower limits on experimental mechanical properties. Finally, the effect of various degrees of covalent crosslinking between the CNTs and amorphous carbon matrix is considered to identify the amount of CNT-matrix covalent bonding that maximizes overall composite properties.

Jensen, Benjamin D.↗