Engineering Papers⌕ Search

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.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Computer Simulation of Proton Transport in Fuel Cell Membranes (Final Report)

This DOE-supported research grant focused on understanding the nature of proton transport in complex systems such as proton exchange membranes (PEMs). The most unique aspect of the research was the development and implementation of a novel multiscale reactive molecular dynamics (MS-RMD) methodology. In this approach, covalent bonds can dynamically break and form, allowing one to accurately treat the proton hopping process essential to capturing the physics of proton transport. The Voth group applied this method to proton exchange membrane systems, providing insight into their proton transport mechanism. They found that protons can diffuse most rapidly in the water-rich regions, but that protons actually spend so little time in such regions that transport along the hydrophobic – hydrophilic interface controls the membrane performance. In addition, the group worked to increase understanding of acidic solutions, developing methods for simulating and interpreting experimental infrared vibrational spectroscopy for excess protons (acidic solutions). The group also implemented novel tools for developing proton transport reactive MD models using a relative entropy minimization scheme.

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

Plasma-assisted chemical-looping combustion: Mechanistic insights into low temperature methane oxidation with CuO

The low-temperature oxidation of CH 4 by CuO in a coaxial, fixed bed, double dielectric barrier discharge (DBD) reactor was investigated with time-dependent species measurements by an electron-ionization molecular beam mass spectrometer (EI-MBMS). In the experiment, 10% methane carried by noble gasses was flown at 50 sccm through 1 g CuO dispersed in quartz wool both under plasma and non-plasma conditions, while time-dependent gas-phase species profiles were collected. Plasma conditions were explored from 300 to 600 °C while the non-plasma conditions were set from 300 to 900 °C. Mechanistic insights into the oxidation of CH 4 by CuO with plasma discharge at lower temperatures (≤ 600 °C) were obtained by quantifying the fuel oxidation, intermediate species, and CO 2 production in comparison to the non-plasma conditions. Here, we observed significant enhancement of fuel oxidation from the plasma discharge between 400 and 500 °C. The CO 2 production at 500 °C with plasma was greater than that at 700 °C without plasma, reducing fuel oxidation temperature by 200+ °C. During tests, three distinct reaction stages were observed: a gas-phase transport limited stage, a surface reaction limited kinetic stage, and an oxygen ion diffusion limited stage. It was observed that plasma greatly improved the reactivity of the second stage at low temperature. In addition, no carbon deposits were observed on the resultant particles, even under the presence of plasma. M. species such as C 4 H 2 and C 6 H 6 not previously observed or predicted in CuO/CH 4 chemical looping were observed, with some species such as CH 3 OH only becoming detectable as total flowrate was increased from 50 to 1500sccm. A non-plasma reaction pathway for CH 4 based the observed species from the MBMS spectrum and previous predictions from reactive molecular dynamics simulations was created, providing a framework from which more complex plasma CuO mechanisms can be crafted from.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Predicting cost-effective carbon fiber precursors: Unraveling the functionalities of oxygen and nitrogen-containing groups during carbonization from ReaxFF simulations

Blends of polyacrylonitrile (PAN) and poly(p-phenylene-2,6-benzobisoxazole) (PBO) as precursors may offer an opportunity to reduce the cost of carbon fiber (CF) production. The all-carbon ring formations and volatile gas production during the heating and carbonization for 9 PAN/PBO blend precursors with varying mixing ratios are studied via the ReaxFF reactive molecular dynamics simulations. Evolutions of oxygen-containing (O-containing) and nitrogen-containing (N-containing) groups are detailed, in order to reveal the reaction mechanisms of the O and N-based CF precursors. Particularly, the O-containing groups are identified to be more efficient for initiating the carbonization, whereas N-containing groups are far longer retained in the graphitic materials and play a key role in capturing and converting carbon radical species into the graphitic networks. Additionally, the PAN/PBO blend precursor with a mole ratio of 1:1 is compared with the pre-oxidized PAN, PAN, and PBO. As a result, it is found that the PAN/PBO blends could be a promising alternative for the cost-effective PAN-based CF precursors, since they can decrease the cost of the CF production by means of: (a) eliminating the pre-oxidation process, (b) having considerable all-carbon ring formations within a short period, and (c) having a relatively fast conversion rate, reaching 95% 6-membered carbon ring formation.

42 ENGINEERING↗
DOE OSTI2019

Converting PBO fibers into carbon fibers by ultrafast carbonization

Here, we report an ultrafast transformation of poly(p-phenylene-2,6-benzobisoxazole) (PBO) fibers into car- bon fibers (CFs) by direct carbonization that eliminates high-cost oxidation in traditional CF conversion. Ultrafast heating and cooling render PBO-derived CFs superlative mechanical properties. ReaxFF atomistic-scale reactive molecular dynamics simulations unveil the PBO to CF conversion mechanism in which ultrafast heating suppresses the O-contained gas release, facilitating all-carbon rings alignment during the carbonization of PBO. These results offer new guidelines for identifying and evaluating alternative CF precursors and fresh perspectives on smart manufacturing of CFs.

42 ENGINEERING↗
DOE OSTI2024

Evolution of PTCDA-derived seeds prior to graphene nanoribbon growth on Ge(001)

The synthesis of graphene nanoribbons (GNRs) can be realized via CH4 chemical vapor deposition (CVD) on substrates such as Ge(001) that promote highly anisotropic growth kinetics. Small polycyclic aromatic hydrocarbons (PAHs) are first sublimed onto Ge near room temperature to form graphene-like seeds that subsequently initiate GNR growth with CH4 exposure at 1173 K. The behaviors of PAHs in the time in between their sublimation onto Ge and GNR growth are unclear. Here, we study an archetypical PAH, perylene-3,4,9,10-tetracarboxyl acid dianhydride (PTCDA), on Ge(001) using both scanning tunneling microscopy (STM) and density functional theory (DFT) to characterize PAH conformation, surface diffusivity, and clustering, versus temperature. PTCDA becomes mobile above 673 K (consistent with a DFT calculated diffusion barrier of 2.01 eV). The mobile PTCDA molecules meet, cluster, and fuse at Ge step edges. These clusters have an apparent height of 0.4 nm, similar to small graphene islands, and grow laterally with increasing temperature – reaching 1.2 – 2.1 nm in width at 1173 K, consistent with the seed size extrapolated from CVD experiments. These results provide a plausible picture for how PTCDA forms seeds for anisotropic GNR CVD and show that PAHs with reduced surface diffusivity and inter-molecular reactivity are needed to enable more monodisperse PAH-seeded GNR synthesis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Boltzmann machine modeling of layered MoS 2 synthesis on a quantum annealer

Two-dimensional layered materials are usually synthesized by chemical vapor deposition (CVD). Here we use the Restricted Boltzmann Machine (RBM) and its enhancement called Limited Boltzmann Machine (LBM) on a quantum computer to model the chemical vapor deposition (CVD) growth data for a MoS 2 monolayer. The CVD data for Boltzmann machines are generated by reactive molecular dynamics (RMD) simulations. The Boltzmann machines give insight into the semiconducting (2H) and metallic (1T) phases and defects generated during the CVD growth of MoS 2 . We compare the performances of the two types of Boltzmann machine models, i.e. RBM versus augmented topologies with limited intra-layer coupling between the hidden units (LBM). We show that limited topologies have a performance advantage over restricted topologies. We examine connectivity within our topological variants, explore hardware qubit mapping schemes, and discuss what performance differences may imply about locality within the data without prior knowledge.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nanomechanical investigation of the interplay between pore morphology and crack orientation of amorphous silica

Porous amorphous silica (a-SiO 2 ) is of both fundamental and practical interests, as they exhibit a large specific surface area and tunable porous network. However, the brittle nature of a-SiO 2 and the presence of pre-existing cracks at both micro- and nano-scales lead to complex mechanical behavior. Here, we systematically investigate the effects of pre-existing crack and its orientation on the mechanical properties of a-SiO 2 with varying pore shapes using reactive molecular dynamics simulations. We demonstrate that pore shape will primarily influence the Young’s modulus (E) and critical energy release rate (G IC ). We further investigate the impact of pore shape and crack orientation by local characterization of the structural parameters. By defining the high stress and inter-mediate regions, the overall mechanical properties are found to be greatly influenced by the pore shape which can be reflected through the spatial distribution of von Mises stress. Overall, G IC is found to increase with the increase of ligament length (also known as pore wall thickness). Meanwhile, the effect of the pre-existing crack on the crack propagation process is confirmed by analyzing the density distribution evolution. These results highlight the interplay between pore morphology and crack orientation in controlling the fracture behaviors in brittle porous materials.

42 ENGINEERING↗
DOE OSTI2021

Topological origin of phase separation in hydrated gels

Depending on their composition, hydrated gels can be homogeneous or phase-separated, which, in turn, affects their dynamical and mechanical properties. However, the nature of the structural features, if any, that govern the propensity for a given gel to phase-separate remains largely unknown. Here, we argue that the propensity for hydrated gels to phase-separate is topological in nature. We employ reactive molecular dynamics simulations to model the early-age precipitation of calcium–alumino–silicate–hydrate (Csingle bondAsingle bondSsingle bondH) gels with varying compositions, i.e., (CaO) 1.7 (Al 2 O 3 ) x (SiO 2 ) 1 – x (H 2 O) 3.7 + x . By adopting topological constraint theory, we investigate the structural origin of phase separation in hydrated gels. We report the existence of a homogeneous-to-phase-separated transition, wherein Si-rich (x ≤ 0.10) Csingle bondAsingle bondSsingle bondH gels are homogeneous, whereas Al-rich (x > 0.10) Csingle bondAsingle bondSsingle bondH gels tend to phase-separate. Furthermore, we demonstrate that this transition is correlated to a topological flexible-to-rigid transition within the atomic network. Finally, we reveal that the propensity for topologically-overconstrained gels to phase-separate arises from the existence of some internal stress within their atomic network, which acts as an energy penalty that drives phase separation.

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

Intramolecular Electrostatic Effects on O 2 , CO 2 , and Acetate Binding to a Cationic Iron Porphyrin

Noncovalent electrostatic interactions are important in many biological and chemical reactions, especially those that involve charged intermediates. There has been a growing interest in using electrostatic ligand designs—placing charges in the second coordination sphere—to improve molecular reactivity, catalysis, and electrocatalysis. For instance, an iron porphyrin bearing four cationic ortho -trimethylanilinium groups, Fe( o -TMA), has been reported to be an exceptional electrocatalyst for both the carbon dioxide reduction reaction (CO 2 RR) and the oxygen reduction reaction (ORR). These reactions involve many different steps, and it is not evident which steps are affected by the four positive charges, or why. By comparing Fe( o -TMA) with the related iron-tetraphenylporphyrin, this work examines how covalently positioned charged groups affect substrate binding and other key pre-equilibria of both the ORR and CO 2 RR, specifically acetate, dioxygen, and carbon dioxide binding. This study is among the first to directly measure the effects of electrostatics on ligand-binding. The results show that adding electrostatic groups to a catalyst design often results in a complex interplay of multiple effects, including changes in pre-equilibria prior to substrate binding, combinations of through-space and inductive contributions, and effects of ionic strength and solution dielectric. The inverse half-order dependence of binding constant on ionic strength is proposed as a clear marker for an electrostatic effect. Overall, the conclusions provide guidance for the increasingly popular electrostatic ligand designs in catalysis and other reactivity.

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

Improving Bond Dissociations of Reactive Machine Learning Potentials through Physics-Constrained Data Augmentation

In the field of computational chemistry, predicting bond dissociation energies (BDEs) presents well-known challenges, particularly due to the multireference character of reactive systems. Many chemical reactions involve configurations where single-reference methods fall short, as the electronic structure can significantly change during bond breaking. As generating training data for partially broken bonds is a challenging task, even state-of-the-art reactive machine learning interatomic potentials (MLIPs) often fail to predict reliable BDEs and smooth dissociation curves. By contrast, simple and inexpensive physics-based models, such as the well-established Morse potential, do not suffer from any such limitations. This work leverages the Morse potential to improve reactive MLIPs by augmenting the training data set with inexpensive Morse data along the dissociation pathways. Further, this physics-constrained data augmentation (PCDA) approach results in MLIPs with smooth bond dissociation curves as well as near coupled-cluster level BDEs, all without requiring any expensive multireference quantum mechanical calculations. A case study for methane combustion demonstrates how the PCDA approach can improve an existing reactive MLIP, namely, ANI-1xnr. In conclusion, not only are the BDEs and bond dissociation curves for all radicals and molecules significantly improved compared to ANI-1xnr but the PCDA-trained MLIP retains the reliability of ANI-1xnr when performing reactive molecular dynamics simulations.

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

Semi-Automated Creation of Density Functional Tight Binding Models through Leveraging Chebyshev Polynomial-Based Force Fields

Density functional tight binding (DFTB) is an attractive method for accelerated quantum simulations of condensed matter due to its enhanced computational efficiency over standard density functional theory (DFT) approaches. However, DFTB models can be challenging to determine for individual systems of interest, especially for metallic and interfacial systems where different bonding arrangements can lead to significant changes in electronic states. In this regard, we have created a rapid-screening approach for determining systematically improvable DFTB interaction potentials that can yield transferable models for a variety of conditions. Our method leverages a recent reactive molecular dynamics force field where many-body interactions are represented by linear combinations of Chebyshev polynomials. This allows for the efficient creation of multi-center representations with relative ease, requiring only a small investment in initial DFT calculations. Here, we have focused our workflow on TiH 2 as a model system and show that a relatively small training set based on unit-cell-sized calculations yields a model accurate for both bulk and surface properties. Our approach is easy to implement and can yield reliable DFTB models over a broad range of thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

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

Insight into the Chemistry of PETN Under Shock Compression Through Ultrafast Broadband Mid-Infrared Absorption Spectroscopy

In this study, thin films of pentaerythritol tetranitrate (PETN) were shock compressed using the laser driven shock apparatus at Los Alamos National Laboratory (LANL). Two spectroscopic probes were available to this apparatus: visible white light transient absorption spectroscopy (VIS) from 400 to 700 nm and mid-infrared transient absorption spectroscopy (MIR) from 1150 to 3800 cm –1 . Important PETN vibrational modes are the symmetric and antisymmetric NO 2 stretches at 1280 and 1650 cm –1 , respectively, as well as CH stretches at ~2900 cm –1 . Shock strength was varied from approximately 3 to 55 GPa to span from the chemically unreactive regime to the regime in which fast chemical reaction took place on the 250 ps time scale of the measurements. VIS and MIR results suggest irreversible chemistry was induced in PETN at pressures above 30 GPa. At lower shock pressures, the spectroscopy showed minimal changes attributable to pressure induced effects. Under the higher-pressure reactive conditions, the frequency region at the antisymmetric NO 2 stretch mode had a significantly increased absorption while the region around the symmetric NO 2 stretch did not. No observable increased absorption occurred in the higher frequency regions where CH-, NH-, and OH- bond absorptions would be observed. A broad absorption appeared on the shoulder at the red-edge of the CO 2 vibrational band around 2200 cm –1 . In addition to the experiments, reactive molecular dynamics were carried out under equivalent shock conditions to correlate the evolution of the infrared spectrum to molecular processes. The simulations show results consistent to experiments up to 30 GPa but suggest that NO and NO 2 related features provided the strongest contributions to the shocked infrared changes. Proposed mechanisms for shocked PETN chemistry are analyzed as consistent or inconsistent with the data presented here. Our experimental data suggests C≡O or N 2 O bond formation, nitrite formation, and absence of significant hydroxyl or amine concentrations in the initial chemistry steps in PETN shocked above 30 GPa.

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

Multiscale Reactive Model for 1,3,5-Triamino-2,4,6-trinitrobenzene Inferred by Reactive MD Simulations and Unsupervised Learning

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

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

Surface Segregation of Liquid Metal Plasma-Facing Component Alloys: A ReaxFF Investigation

Engineering liquid metal alloys offers a transformative pathway for plasma-facing components by enabling chemically tailored surfaces that can simultaneously optimize plasma-material interactions, reduce divertor heat flux, and enhance core plasma confinement, thereby advancing the commercial viability of nuclear fusion power plants. This study, employing an atomistic simulation approach, provides direct evidence that incorporating nonmetal surface-active agents (such as O and H, or their combination) enables strong surface segregation. This capability makes tin−aluminum (Sn−Al) and tin−lithium (Sn−Li) alloys, with suitable compositions, good candidates for PFC applications. Specifically, the presence of low-Z solutes (Li, Al) leads to preferential surface enrichment, which imparts low-Z sputtering characteristics, while the Sn solvent maintains thermophysical stability. To systematically examine this behavior, we developed a ReaxFF force field spanning the full Sn/Al/Li/O/H chemistry, validated it against formation energies and elastic constants, and applied it in reactive molecular dynamics simulations at fusion-relevant temperatures. We also introduced an overlapbased segregation index that captures interfacial compositional separation directly from atomistic density distributions. Here, this metric reveals a clear hierarchy of segregation regimes and provides a unified view across all systems studied. Together, these findings establish a mechanistic link between nonmetal chemistry and interfacial structure, providing a predictive framework for designing self-adaptive, low-sputtering liquid metal alloys for fusion applications.

Alloys↗
DOE OSTI2022

Extemporaneous Mechanochemistry: Shock-Wave-Induced Ultrafast Chemical Reactions Due to Intramolecular Strain Energy

We report regions of energy localization referred to as hotspots are known to govern shock initiation and the run-to-detonation in energetic materials. Mounting computational evidence points to accelerated chemistry in hotspots from large intramolecular strains induced via the interactions between the shock wave and microstructure. However, definite evidence mapping intramolecular strain to accelerated or altered chemical reactions has so far been elusive. From a large-scale reactive molecular dynamics simulation of the energetic material 1,3,5-triamino-2,4,6-trinitrobenzene, we map decomposition kinetics to molecular temperature and intramolecular strain energy prior to reaction. Both temperature and intramolecular strain are shown to accelerate chemical kinetics. A detailed analysis of the atomistic trajectory shows that intramolecular strain can induce a mechanochemical alteration of decomposition mechanisms. The results in this paper could inform continuum-level chemistry models to account for a wide range of mechanochemical effects.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Surface Transfer Doping in MoO 3– x /Hydrogenated Diamond Heterostructure

Surface transfer doping is proposed to be a potential solution for doping diamond, which is hard to dope for applications in high-power electronics. While MoO 3 is found to be an effective surface electron acceptor for hydrogen-terminated diamond with a negative electron affinity, the effects of commonly existing oxygen vacancies remain elusive. We have performed reactive molecular dynamics simulations to study the deposition of MoO 3–x on a hydrogenated diamond (111) surface and used first-principles calculations based on density functional theory to investigate the electronic structures and charge transfer mechanisms. We find that MoO 3–x is an effective surface electron acceptor and the spatial extent of doped holes in hydrogenated diamond is extended, promoting excellent transport properties. Charge transfer is found to monotonically decrease with the level of oxygen vacancy, providing guidance for engineering of the surface transfer doping process.

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

ReaxFF Simulations of Laser-Induced Graphene (LIG) Formation for Multifunctional Polymer Nanocomposites

Irradiation of polymer films by a CO 2 infrared laser under ambient conditions converts the polymer into porous graphene or laser-induced graphene (LIG). In this work, we simulate the formation of LIG from five different commercially available polymers using reactive molecular dynamics. We determined that the molecular structure of the parent polymer has a significant effect on the final graphitic structure. CO is liberated during the initial part of the LIG formation process when the polymer is converted into an amorphous structure, while H 2 is evolved steadily as the amorphous structure is converted to an ordered graphitic structure. The LIG structure has out-of-plane undulations and bends due to a significant number of 5- and 7-member carbon rings present throughout the structure. We find that the simulated molecular structure compares well with recent experimental observations from the literature. We also demonstrate that the yield of LIG is higher in inert conditions, compared to environments with oxygen. Polybenzimidazole-derived LIG has the highest surface area and yield among the five polymers examined. These findings provide knowledge of LIG formation mechanisms that can be leveraged for bulk LIG applications such as sensors, electrocatalysts, microfluidics, and targeted heating for welding polymers.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Atomistic Mechanisms of Stress-Dependent Molten Salt Corrosion in NiCr Alloys

Ni-based structural alloys in molten salt environments often experience simultaneous mechanical loading and corrosive attack, yet the mechanisms governing stress-corrosion interactions remain unclear. Prior studies largely emphasize tensile stress, while the role of compressive stress has received limited attention. Here, reactive molecular dynamics simulations are used to investigate the coupled effects of applied strain and corrosion in Ni 0.75 Cr 0.25 exposed to molten FLiNaK at 800 °C. A Σ5(210) grain boundary model is subjected to tensile (+4%) to compressive (−4%) uniaxial strains, and corrosion behavior is evaluated through fluorine adsorption, charge redistribution, and grain boundary evolution. Tensile strain accelerates intergranular corrosion susceptibility by reducing local atomic packing through elastic dilation and increasing excess free volume at the grain boundary, which enhances atomic mobility and salt infiltration. In contrast, compressive strain can suppress corrosion by promoting the formation of a ridge-like surface layer along the grain boundary, limiting salt access to the underlying alloy. These results provide atomistic insight into how stress states influence grain boundary corrosion in molten salts.

36 - MATERIALS SCIENCE↗