A Reactive Molecular Dynamics Study of Phenol and Phenolic Polymers in Extreme Environments.
Abstract not provided.
Jones, Keith Andrew↗
SEARCH · Engineering Papers
Results for “molecular reactivity”
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.
Surrogate Fuel Chemistry from ReaxFF Molecular Dynamics
Reactive molecular dynamics with ReaxFF simulations are used to study the chemistry of four component jet fuel mixtures of n-dodecane, iso-octane, n-propylbenzene, and mesitylene. ReaxFF does not require knowledge of the chemistry a priori providing insight into reactions otherwise not seen in established chemical kinetic mechanisms. The results show that a difference in the local chemical environment can change the fuel decomposition pathways across the surrogate mixtures where the same fuel species are used. These simulations highlight the need to consider cross-chemical interactions in chemical kinetic mechanism development and the importance of resonance stabilized polyaromatic hydrocarbon reactions for soot growth and formation. A method for predicting the threshold sooting index from the carbon species formed at the end of the simulations is established. The results of the simulations are consistent with measurements of the threshold sooting index with the predicted values within 20% of the measured data
Surrogate Fuels↗
InteraChem: Virtual Reality Visualizer for Reactive Interactive Molecular Dynamics
Interactive molecular dynamics in virtual reality (IMD-VR) simulations provide a digital molecular playground for students as an alternative or complement to traditional molecular modeling kits or 2D illustrations. Previous IMD-VR studies have used molecular mechanics to enable simulations of macromolecules such as proteins and nanostructures for the classroom setting with considerable success. Here, we present the InteraChem molecular visualizer, intended for reactive IMD-VR simulation using semiempirical and ab initio methods. InteraChem visualizes not only the molecular geometry but also (1) isosurfaces such as molecular orbitals and electrostatic potentials and (2) two-dimensional graphs of time-varying simulation quantities such as kinetic/potential energy, internal coordinates, and user-applied force. Additionally, InteraChem employs speech recognition to facilitate user interaction and introduces a novel “atom happiness” visualization using emojis to indicate the energetic feasibility of a particular bonding arrangement. Furthermore, we include a set of accompanying exercises that we have used to teach chemical reactivity in small molecular systems.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
On the Network Topology of Cross-Linked Acrylate Photopolymers: A Molecular Dynamics Case Study
A reactive molecular dynamics approach is used to simulate cross-linking of acrylate polymer networks. By employing the same force field and reactive scheme and studying three representative multifunctional acrylate monomers, we isolate the importance of the nonreactive moieties within these model monomers. Analyses of reactive trajectories benchmark the estimated gel points, cyclomatic character, and spatially resolved cross-linking tendencies of the acrylates as a function of conversion. The insights into the similarities and differences of the polymerization and resulting networks presented here suggest molecular mechanics as a useful tool in the rational design of photopolymerization resins.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
New insights into the mechanism governing the elasticity of calcium silicate hydrate gels exposed to high temperature: A molecular dynamics study
When exposed to fire, the integrity of cement-based materials is governed by thermally-induced changes in the mechanical properties of their binding phase, i.e., the calcium–silicate–hydrate (C–S–H) gel. However, the effect of temperature on the structure, density, and mechanical properties of C–S–H remains only partially known. Here, based on reactive molecular dynamics simulations, we reveal the nature of thermally-induced damage in C–S–H gels. In general, we show that, at the atomic scale, exposure to high temperature results in partial dehydration, volumetric shrinkage, disordering, and stiffening in the C–S–H grains. However, we show that the thermal response of C–S–H strongly depends on its chemical composition, wherein C–S–H systems associated with lower Ca/Si molar ratios are able to undergo higher temperatures before amorphization. Based on these results, we demonstrate that the stiffness of C–S–H gels (i.e., including porosity—as probed by nanoindentation) is governed by a competition between the stiffening of the grains and the decrease in packing density—wherein the latter eventually become predominant.
36 MATERIALS SCIENCE↗
Unveiling the Mechanism of Mn Dissolution Through a Dynamic Cathode‐Electrolyte Interphase on LiMn2O4
Abstract Understanding the formation and evolution of the cathode‐electrolyte interphase (CEI), which forms at the interface between the cathode and electrolyte, is crucial for revealing degradation mechanisms in cathode materials, especially for developing strategies to stabilize the interphase in the strongly oxidizing conditions that evolve at high operating voltages in next‐generation Li‐ion batteries. However, The present understanding of the CEI is challenged by its complex and dynamic nature. In this work, near‐edge X‐ray absorption fine structure spectroscopy, electrochemical characterization, and reactive molecular dynamics simulations are combined to reveal a mechanism for CEI formation and evolution above model LiMn 2 O 4 (LMO) thin‐film electrodes in contact with conventional carbonate‐based electrolytes. It is found that Mn dissolution from LMO can be understood in terms of repetitive Mn 3 O 4 formation and dissolution behavior during cycling, which is closely connected to electrolyte decomposition and a key aspect of the CEI formation and growth. The behavior of the CEI in this model system offers detailed insight into the dynamic chemistry of the interphase, underscoring the important role of electrolyte composition and cathode surface structure in interphase degradation.
Ou, Wenhan↗
Implementing reactivity in molecular dynamics simulations with harmonic force fields
Not provided.
Science & Technology - Other Topics↗
Prediction of carbon nanostructure mechanical properties and the role of defects using machine learning
Graphene-based nanostructures hold immense potential as strong and lightweight materials, however, their mechanical properties such as modulus and strength are difficult to fully exploit due to challenges in atomic-scale engineering. This study presents a database of over 2,000 pristine and defective nanoscale CNT bundles and other graphitic assemblies, inspired by microscopy, with associated stress–strain curves from reactive molecular dynamics (MD) simulations using the reactive INTERFACE force field (IFF-R). These 3D structures, containing up to 80,000 atoms, enable detailed analyses of structure-stiffness-failure relationships. By leveraging the database and physics- and chemistry-informed machine learning (ML), accurate predictions of elastic moduli and tensile strength are demonstrated at speeds 1,000 to 10,000 times faster than efficient MD simulations. Hierarchical Graph Neural Networks with Spatial Information (HS-GNNs) are introduced, which integrate chemistry knowledge. HS-GNNs as well as extreme gradient boosted trees (XGBoost) achieve forecasts of mechanical properties of arbitrary carbon nanostructures with only 3 to 6% mean relative error. The reliability equals experimental accuracy and is up to 20 times higher than other ML methods. Predictions maintain 8 to 18% accuracy for large CNT bundles, CNT junctions, and carbon fiber cross-sections outside the training distribution. The physics- and chemistry-informed HS-GNN works remarkably well for data outside the training range while XGBoost works well with limited training data inside the training range. The carbon nanostructure database is designed for integration with multimodal experimental and simulation data, scalable beyond 100 nm size, and extendable to chemically similar compounds and broader property ranges. The ML approaches have potential for applications in structural materials, nanoelectronics, and carbon-based catalysts.
Winetrout, Jordan J.↗
Atomistic-Scale Simulations on Graphene Bending Near a Copper Surface
Molecular insights into graphene-catalyst surface interactions can provide useful information for the efficient design of copper current collectors with graphitic anode interfaces. As graphene bending can affect the local electron density, it should reflect its local reactivity as well. Using ReaxFF reactive molecular simulations, we have investigated the possible bending of graphene in vacuum and near copper surfaces. We describe the energy cost for graphene bending and the binding energy with hydrogen and copper with two different ReaxFF parameter sets, demonstrating the relevance of using the more recently developed ReaxFF parameter sets for graphene properties. Moreover, the draping angle at copper step edges obtained from our atomistic simulations is in good agreement with the draping angle determined from experimental measurements, thus validating the ReaxFF results.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Computationally guided synthesis of carbon coated mesoporous silica materials
Mesoporous silica materials (MSMs) are well-suited for biomedical applications due to their unique features, including a large surface area and tunable pore size. To enhance their durability, the small pores in MSMs are filled with carbon precursors and then carbonized to prevent them from interacting with unreacted silicic acid. Here, in this study, we synthesized and healed MSMs using a combination of non-reactive and reactive molecular dynamics (MD) simulations. The non-reactive MD simulation revealed that the self-assembly of Pluronic® L64 polymers in water resulted in nearly 80 % hydrogen bonds between the hydrophilic sections of the micelle and water. In the bond-boosted ReaxFF MD simulations, silicic acid precursors were condensed on the micelle surface, with over 60 % of them leading to the creation of periodic mesoporous silica within the system. Condensation of silicic acid precursors at 300 K with bond-boosting and at 1500 K without it both significantly promoted the polymerization of Si(OH) 4 , with the latter doubling the rate compared to the former. Subsequently, we healed the MSM surface by carbonizing carbon precursors inside an MSM pore. Polyethylene (PE) and high-rank lignite were identified as the most suitable precursors due to their ability to form turbostratic graphene structures. High-rank lignite exhibited the highest carbon conversion ratio to 6-membered rings, closely followed by PE, in the carbon ring formation analysis. Additionally, the production of gases, such as H 2 , increased significantly for PE at both 2200 K and 2600 K, indicating the conversion of a considerable portion of carbon into graphitic or turbostratic structures. The carbonization of PE primarily led to the formation of planar (sp 2 ) structures, while sucrose yielded the least planar structures. Finally, we studied the protective blocking of unreacted silicic acid precursor by considering a PET turbostratic graphene structure in a silica mesopore formed at 2600 K. The trajectory analysis showed that the surface of the silica was effectively coated with PET tar, preventing unreacted silicic acid from interacting with the inner silica pore surface. These findings offer valuable insights into the synthesis and carbonization-based healing processes of MSMs, enhancing their potential for various biomedical applications.
36 MATERIALS SCIENCE↗
Atomistic Mechanics of Torn Back Folded Edges of Triangular Voids in Monolayer WS 2
Triangular nanovoids in 2D materials transition metal dichalcogenides have vertex points that cause stress concentration and lead to sharp crack propagation and failure. Here, the atomistic mechanics of back folding around triangular nanovoids in monolayer WS 2 sheets is examined. Combining atomic-resolution images from annular dark-field scanning transmission electron microscopy with reactive molecular modelling, it is revealed that the folding edge formation has statistical preferences under geometric conditions based on the orientation mismatch. It is further investigated how loading directions and strong interlayer friction, interplay with WS 2 lattice's crack preference, govern the deformation and fracture pattern around folding edges. Finally, these results provide fundamental insights into the combination of fracture and folding in flexible monolayer crystals and the resultant Moiré lattices.
36 MATERIALS SCIENCE↗
Materials genome innovation for computational software (magics) center
Functional layered material (LM) architectures will dominate nanomaterials science in this century. We have developed theory, modeling, simulation, and software and data tools that enhance understanding and AI guide synthesis, enable characterization of complex structures, and improve capabilities in the predictive design and growth of LMs. Research at the Center has focused on: Computational synthesis and characterization: AI guided synthesis and experimental synthesis of stacked LMs with tailored properties via optimized chemical vapor deposition (CVD) growth and liquid-phase exfoliation; study defects, edges, grain boundaries, wrinkling of atomic layers and their effects on chemical, mechanical, electrical, and optical properties. Far-from-equilibrium processes: Joint experimental and simulation based probe of electronic processes with NAQMD and ultrafast X-ray free-electron laser (XFEL) and ultrafast electron diffraction (UED) facilities at Stanford. Experimentally validate NAQMD by ultrafast electron diffraction and X-ray spectroscopy studies of structural and excited state dynamics, shape fluctuations, and phonon dynamics. Scalable software: Simulation engines for desktop-to-exascale platforms using low-overhead, linear-scaling QMD algorithms; divide-conquer-recombine NAQMD with electronic excitations; extended-Lagrangian reactive molecular dynamics (RMD), machine learning (ML) based neural-network quantum molecular dynamics (NNQMD), and super-state accelerated molecular dynamics (AMD) and kinetic Monte Carlo codes; thermal and electrical transport software; and design 3D architectures of LMs with desired functionality using scalable software. Distribution of software and data, and training: Software and simulation-experimental data generated within the Center are distributed to the materials science community via Berkeley Materials Project (MP) framework. We have also organized three workshops for software distribution and training at USC (Nov. 2017, Mar. 2018) and Gaithersburg, MD (Nov. 2018) to train researchers, with the last one in focused on underrepresented groups, in collaboration with Howard University which is one of the largest HBCUs. The Center supported a total of 46 personnel and 6 undergraduate students. These include 14 faculty, 11 postdoctoral research associates, 20 graduate research assistants, and mentored 6 undergraduate students. This resulted in the publications of 63 research papers that include 46 publications on Reactive and Quantum Dynamics Simulations, 13 publications on Machine Learning for Quantum Materials, and 4 publications on Quantum Computing.
2D Materials↗
Theoretical studies of molecular collisions
The following subject areas are covered: (1) total integral reactive cross sections and vibrationally resolved reaction probabilities for F + H2 = HF + H; (2) a theoretical study of inelastic O + N2 collisions; (3) body frame close coupling wave packet approach to gas phase atom-rigit rotor inelastic collisions; (4) wave packet study of gas phase atom-rigit motor scattering; (5) the application of optical potentials for reactive scattering; (6) time dependent, three dimensional body frame quantal wave packet treatment of the H + H2 exchange reaction; (7) a time dependent wave packet approach to atom-diatom reactive collision probabilities; (8) time dependent wave packet for the complete determination of s-matrix elements for reactive molecular collisions in three dimensions; (9) a comparison of three time dependent wave packet methods for calculating electron-atom elastic scattering cross sections; and (10) a numerically exact full wave packet approach to molecule-surface scattering.
Kouri, Donald J.↗
Strain-Induced Growth of Twisted Bilayers during the Coalescence of Monolayer MoS 2 Crystals
Tailoring the grain boundaries (GBs) and twist angles between two-dimensional (2D) crystals are two crucial synthetic challenges to deterministically enable envisioned applications such as moiré excitons, emerging magnetism, or single-photon emission. Here, we reveal how twisted 2D bilayers can be synthesized from the collision and coalescence of two growing monolayer MoS 2 crystals during chemical vapor deposition. The twisted bilayer (TB) moiré angles are found to preserve the misorientation angle (θ) of the colliding crystals. The shapes of the TB regions are rationalized by a kink propagation model that predicts the GB formed by the coalescing crystals. Optical spectroscopy measurements reveal a θ-dependent long-range strain in crystals with stitched grain boundaries and a sharp (θ > 20°) threshold for the appearance of TBs, which relieves this strain. Reactive molecular dynamics simulations explain this strain from the continued growth of the crystals during coalescence due to the insertion of atoms at unsaturated defects along the GB, a process that self-terminates when the defects become saturated. The simulations also reproduce atomic-resolution electron microscopy observations of faceting along the GB, which is shown to arise from the growth-induced long-range strain. These facets align with the axes of the colliding crystals to provide favorable nucleation sites for second-layer growth of a TB with twist angles that preserve the misorientation angle θ. This interplay between strain generation and aligned nucleation provides a synthetic pathway for the growth of TBs with deterministic angles.
2D materials↗
Atomic understanding of structural deformations upon ablation of graphene
We investigate the atomic rearrangement in graphene under femtosecond pulse illumination with reactive molecular dynamics simulations and compare with ultra-fast laser ablation experiments. To model the impact of the laser pulse irradiation, heat is locally applied to a selected area of the graphene layer and the resulting structural deformation is simulated as a function of time, providing a detailed understanding of the bond breaking process under laser illumination and subsequent re-equilibration afterAQ3 the pulse is turned off. Analysis of the atomic dynamics indicates that the types of defects formed depend on the pulse energy and exposure duration. By varying the exposed area, we determine that the shape of the ablated area is not only a function of the pulse energy, but also of the beam spot size and pulse repetition. Furthermore, we apply a machine learning approach to extrapolate our simulated data to experimental length scales and reproduce the trends in ablated area as a function of temperature. Furthermore, our study provides a first step towards understanding the design parameters for graphene nano-patterning.
77 NANOSCIENCE AND NANOTECHNOLOGY↗
From atomistic models to machine learning: Predictive design of nanocarbons under extreme conditions
The formation of technologically valuable nanocarbon structures under extreme conditions, such as those produced during high-explosive detonations, remains poorly understood but holds significant potential for the development of controlled synthesis pathways. While detonation shockwaves provide the high-pressure, high-temperature environment required for nanodiamond formation, subsequent cooling and decompression dictate whether the diamond phase is preserved or transformed into other nanocarbon structures. Here, in this study, we employ GPU-accelerated reactive molecular dynamics (ReaxFF) simulations to investigate the graphitization and structural remodeling of detonation nanodiamond under nonlinear quench and pressure-release trajectories. We further investigate how the initial nanodiamond morphology; cuboctahedral, octahedral, or hexagonal prism influences the resulting transformation products. Evolution of nanostructure, allotrope (via simulated x-ray diffraction), carbon hybridization, and ring statistics are tracked during a two-stage quench from 5000 K to 60 GPa. Rapid cooling combined with slow decompression optimizes cubic diamond retention, whereas slow cooling with rapid pressure release promotes surface-to-core graphitization, producing concentric sp 2 -hybridized layers and hollowed inner shells. Octahedral nanodiamonds evolve into carbon nano-onions, initially forming bucky diamonds that progressively transform into fully sp 2 -hybridized structures, while hexagonal prisms preferentially form parallel-stacked graphite layers resembling carbon dots. Transient hexagonal diamond (lonsdaleite) emerges as an interfacial phase, suggesting potential reversibility in the shock-induced graphite-to-diamond transformation pathway transformation route. To extend predictive capabilities, we trained machine learning (ML) regressors on over 10 5 node-hours of molecular dynamics (MD) trajectories. A multilayer perceptron (MLP) model reliably predicts the number of graphitized layers from temperature–pressure trajectories with a coefficient of determination (R 2 ) exceeding 0.90. This high predictive fidelity enables efficient, high-throughput mapping of the synthesis parameter space for optimized graphitization outcomes. Collectively, morphological control combined with optimized quench–decompression conditions promote the selective synthesis of nanocarbon allotropes. This work establishes a data-driven framework for the rational, a priori design of carbon nanomaterials for applications in energy storage, sensing, and biomedicine.
Detonation nanodiamond remodeling↗
Novel polyelectrolytes
Cationic polyelectrolytes are formed by the polymerization in absence of oxygen of a monomer of the general formula: ##STR1## where x is 3 or more than 6 and Z is I, Br or Cl to form high charge density linear polymers. Segments of the linear polymer may be attached to or formed in the presence of polyfunctional reactive tertiary amines or halogen polymeric substrates or polyfunctional lower molecular reactive polyfunctional substrates to form branched or star polyelectrolytes by a quaternization polymerization reaction.
Rembaum, Alan↗
A model for molecular flow of reactive species through cylindrical reactors
Analytical model for molecular flow of reactive gaseous species through cylindrical reactors with apertures at either end, evaluating surviving flux magnitude
Baker, B. R.↗