Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactive molecular dynamics”

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 253 records · Page 14

Self-terminating, heterogeneous solid–electrolyte interphase enables reversible Li–ether cointercalation in graphite anodes

Ether solvents are suitable for formulating solid-electrolyte interphase (SEI)-less ion-solvent cointercalation electrolytes in graphite for Na-ion and K-ion batteries. However, ether-based electrolytes have been historically perceived to cause exfoliation of graphite and cell failure in Li-ion batteries. In this study, we develop strategies to achieve reversible Li–solvent cointercalation in graphite through combining appropriate Li salts and ether solvents. Specifically, we design 1M LiBF 4 1,2-dimethoxyethane (G1), which enables natural graphite to deliver ~91% initial Coulombic efficiency and >88% capacity retention after 400 cycles. We captured the spatial distribution of LiF at various length scales and quantified its heterogeneity. The electrolyte shows self-terminated reactivity on graphite edge planes and results in a grainy, fluorinated pseudo-SEI. The molecular origin of the pseudo-SEI is elucidated by ab initio molecular dynamics (AIMD) simulations. The operando synchrotron analyses further demonstrate the reversible and monotonous phase transformation of cointercalated graphite. Our findings demonstrate the feasibility of Li cointercalation chemistry in graphite for extreme-condition batteries. The work also paves the foundation for understanding and modulating the interphase generated by ether electrolytes in a broad range of electrodes and batteries.

25 ENERGY STORAGE↗

Ensembles of Metastable States Govern Heterogeneous Catalysis on Dynamic Interfaces

Heterogeneous catalysis is at the heart of chemical industry. Being able to tune and design efficient catalysts for processes of interest is of an utmost importance, and for this, the molecular-level understanding of heterogeneous catalysts is the first step, and indeed a prime focus of the modern catalysis research. For a long time, a single most thermodynamically stable structure of the catalytic interface attained in reaction conditions had been envisioned as the reactive phase. However, some catalytic interfaces continue to undergo structural dynamics in the steady state, triggered by high temperatures, pressures, and binding and changing reagents. Among particularly dynamic interfaces are such widely-used catalysts as crystalline and amorphous surfaced supporting (sub-)nano metallic clusters. Recently, it became clear that this dynamic fluxionality causes the supported clusters to populate many distinct structural and stoichiometric states in catalytic conditions. Hence, the catalytic interface should be viewed as an evolving statistical ensemble of many (not one) structures. As such, every member in the ensemble contributes to the properties of the catalyst differently, and in proportion to its probability to be populated. This new notion flips the established paradigm and calls for new theory, modeling approaches, operando measurements, and updated design strategies. The statistical ensemble nature of surface-supported sub-nano cluster catalysts can be exemplified by oxide-supported and adsorbate-covered Pt, Pd, Cu, CuPd clusters, catalytic toward oxidative and non-oxidative dehydrogenation. They have access to a variety of 3D and quasi-2D shapes. The compositions of their thermal ensembles are dependent on the cluster size, leading to size-specific catalytic activities and the famous “every atom counts” phenomenon. The support and adsorbates affect catalyst structures, and state of the reacting species causes the ensemble to change in every reaction intermediate. The most stable member of the ensemble dominates the thermodynamic properties of the corresponding intermediate, whereas the kinetics can be determined by more active but less populated metastable catalyst states, and that suggests that many earlier studies might have overlooked the actual active sites. Both effects depend on the relative time-scales of catalyst restructuring and reaction dynamics. The catalyst may routinely operate off-equilibrium. Ensemble phenomena lead to surprising exceptions from established rules of catalysis, such as scaling relations, and Arrhenius behavior. Catalyst deactivation is also an ensemble property, and its extent of mitigation can be predicted through the new paradigm. These findings were enabled by advances in theory, such as global optimization and subsequent utilization of multiple local minima, and pathways sampling, as well as operando catalyst characterization. The fact that the per-site and per-species resolution is needed for the description and predicting of catalyst properties gives theory the central role in catalysis research, as most experiments provide ensemble-average information and cannot detect the crucial minority species that may be responsible for the catalytic activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Boron nitride nanotube precursor formation during high-temperature synthesis: kinetic and thermodynamic modelling

In this work, we performed integrated modelling of the chemical pathways of formation for boron nitride nanotube (BNNT) precursors during high-temperature synthesis in a B/N 2 mixture. Integrated modelling includes quantum chemistry, Quantum–classical molecular dynamics, thermodynamic modelling, and kinetic approaches. We demonstrate that BN compounds are formed via the interaction of molecular nitrogen with small boron clusters, rather than through interactions with less reactive liquid boron. (This process can also be described as N 2 molecule fixation.) Liquid boron evaporates to produce these boron clusters (B m with m ≤ 5), which are subsequently converted into B m N n chains. The production of such chains is crucial to the growth of BNNTs because these chains form the building blocks of bigger and longer BN chains and rings, which are in turn the building blocks of fullborenes and BNNTs. Additionally, kinetic modelling revealed that B 4 N 4 and B 5 N 4 species in particular play a major role in the N 2 molecule fixation process. Furthermore, the formation of these species via reactions with B 4 and B 5 clusters is not adequately described under the assumption of thermodynamic equilibrium, as is demonstrated in our kinetic modelling. Thus, the accumulation of both B 4 N 4 and B 5 N 4 depends on the background gas pressure and the gas cooling rate. Long BN chains and rings, which are precursors of the fullborene and BNNT growth, form via self-assembly of components B 4 N 4 and B 5 N 4 . Our modelling results—particularly the increased densities of B 4 N 4 and B 5 N 4 species at higher gas pressures—explain the experimentally observed effect of gas pressure on the yield of high-quality BNNTs. The catalytic role of hydrogen was also studied; it is shown that HBNH molecules can be the main precursor of BNNT synthesis in the presence of hydrogen.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Off-Stoichiometric Restructuring and Sliding Dynamics of Hexagonal Boron Nitride Edges in Conditions of Oxidative Dehydrogenation of Propane

Boron-containing materials, such as hexagonal boron nitride (h-BN), recently shown to be active and selective catalysts for the oxidative dehydrogenation of propane (ODHP), have been shown to undergo significant surface oxyfunctionalization and restructuring. Although experimental ex situ studies have probed the change in chemical environment on the surface, the structural evolution of it under varying reaction conditions has not been established. Herein, we perform global optimization structure search with a grand canonical genetic algorithm to explore the chemical space of off-stoichiometric restructuring of the h-BN surface under ambient as well as ODHP-relevant conditions. A grand canonical ensemble representation of the surface is established, and the predicted 11B solid-state NMR spectra are consistent with previous experimental reports. In addition, we investigated the relative sliding of h-BN sheets and how it influences the surface chemistry with ab initio molecular dynamics simulations. Furthermore, the B–O linkages on the edges are found to be significantly strained during the sliding, causing the metastable sliding configurations to have higher reactivity toward the activation of propane and water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rotation in attosecond vibronic coherence spectroscopy for molecules

Abstract Excitation or ionization of a molecule by ultrafast laser pulses can create a superposition of electronic states, whose dynamics is influenced by the interplay of electronic coherence and nuclear motion, resulting in charge migration and possibly charge transfer. Probing the vibronic coherence is therefore vital to monitoring electronic dynamics and controlling chemical reactivity, as recently demonstrated in molecules via attosecond transient absorption spectroscopy (ATAS). However, theories supporting the interpretation of ATAS experiments neglect the effects of molecular rotation, often leading to inaccurate interpretation of experimental data. Here, we develop a comprehensive theory for ATAS of molecules encompassing the entire pump-probe process. Applying the theory to N 2 , we demonstrate that the emergence of coherent signals critically depends on the consideration of molecular rotation. This work contributes to close the gap between theory and ATAS experiments, paving the way for monitoring electronic motion and controlling chemical reactivity in diverse molecular systems.

Physics↗

Towards the Accurate Thermodynamic Characterization of Enzyme Reaction Mechanisms

We employed QM/MM molecular dynamics (MD) simulations to characterize the rate-limiting step of the glycosylation reaction of pancreatic α-amylase with combined DFT/molecular dynamics methods (PBE/def2-SVP : AMBER). Upon careful choice of four starting active site conformations based on thorough reactivity criteria, Gibbs energy profiles were calculated with umbrella sampling simulations within a statistical convergence of 1 -2 kcal · mol -1 . Nevertheless, Gibbs activation barriers and reaction energies still varied from 11.0 to 16.8 kcal · mol -1 and -6.3 to +3.8 kcal · mol -1 depending on the starting conformations, showing that despite significant state-of-the-art QM/MM MD sampling (0.5 ns/profile) the result still depends on the starting structure. In conclusion, the results supported the one step dissociative mechanism of Asp197 glycosylation preceded by an acid-base reaction by the Glu233, which are qualitatively similar to those from multi-PES QM/MM studies, and thus support the use of the latter to determine enzyme reaction mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Formation of Benzene and Naphthalene through Cyclopentadienyl-Mediated Radical–Radical Reactions

Resonantly stabilized free radicals (RSFRs) have been contemplated as fundamental molecular building blocks and reactive intermediates in molecular mass growth processes leading to polycyclic aromatic hydrocarbons (PAHs) and carbonaceous nanoparticles on Earth and in deep space. Here, by combining molecular beams and computational fluid dynamics simulations, we provide compelling evidence on the formation of benzene via the cyclopentadienyl-methyl reaction and of naphthalene through the cyclopentadienyl self-reaction, respectively. These systems offer benchmarks for the conversion of a five-membered ring to the 6π-aromatic (benzene) and the generation of the simplest 10π-PAH (naphthalene) at elevated temperatures. These results uncover molecular mass growth processes from the “bottom up” via RSFRs in high temperature circumstellar environments and combustion systems expanding our fundamental knowledge of the organic, hydrocarbon chemistry in our universe.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomic cluster expansion potential for large scale simulations of hydrocarbons under shock compression

We present an Atomic Cluster Expansion (ACE) machine learned potential developed for high-fidelity atomistic simulations of hydrocarbons, targeting pressures and temperatures near and above supercritical fluid regimes for molecular fluids. A diverse set of stoichiometries were covered in training, including 1:0 (pure carbon), 1:4 (methane), and 1:1 (benzene), and rich bonding environments sampled at supercritical temperatures, hydrogen rich, reactive mixtures where metastable stoichiometries arise, including 1:2 (ethylene) and 1:3 (ethane). A high-fidelity training database was constructed by performing large-scale quantum molecular dynamic simulations [density functional theory (DFT) MD] of diamond, graphite, methane, and benzene. A novel approach to selecting structures from DFT MD is also presented, which allows for the rapid selection of unique DFT MD frames from complex trajectories. Comparisons to DFT and experimental data demonstrate that the presented ACE potential accurately reproduces isotherms, carbon melting curves, radial distribution functions, and shock Hugoniots for carbon and hydrocarbon systems for pressures up to 100 GPa and temperatures up to 6000 K for hydrocarbon systems and up to 9000 K for pure carbon systems. This work delivers a potential that can be used for accurate, large-scale simulations of shocked hydrocarbons and demonstrates a methodology for fitting and validating machine learning interatomic potentials to complex molecular environments, which can be applied to energetic materials in future works.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Level Understanding of Polyethylene Terephthalate (PET) Depolymerization in Base/Alcohol Hybrid Systems

Polyethylene terephthalate (PET) depolymerization in base/alcohol hybrid systems represents a promising low-energy approach for chemically recycling PET waste into valuable monomers. This study investigates the mechanistic pathways of PET depolymerization in NaOH/alcohol solutions, emphasizing the competing roles of hydroxide and alkoxide species. Utilizing a combination of experimental techniques, density functional theory (DFT) calculations, and molecular dynamics (MD) simulations, we explore how factors such as base concentration, alcohol chain length, and pKa values of alcohols influence PET depolymerization efficiency and pathways. Our findings indicate that alkoxide ions (RO⁻) exhibit notably higher reactivity than hydroxide ions (HO⁻), favoring an alcoholysis pathway in the base/alcohol hybrid system. Experimental results across a series of C1 to C5 alcohols show that longer-chain alcohols, particularly 1-butanol, achieve higher PET conversion, although this does not align solely with simple nucleophilicity trends of alkoxides. While DFT calculations reveal comparable activation energies for various alkoxides in PET depolymerization, MD simulations underscore the significant role of alcohol chain length, with longer-chain alcohols forming more stable or frequent interactions with PET. Additionally, the alkoxide concentration, influenced by the alcohol’s pK a , directly impacts PET conversion. These suggest that PET depolymerization is governed by a balance between alkoxide concentration and alkoxide-PET interactions, rather than activation energies or nucleophilicity alone. From a practical perspective, incorporating long-chain alcohols as cosolvents may enhance process efficiency but increases raw material costs by approximately 30%. However, long-chain alcohols present a safer and more sustainable alternative to hazardous cosolvents such as dichloromethane. This work offers a molecular-level understanding of PET depolymerization in base/alcohol systems and provides insights into optimizing these systems for more efficient and sustainable PET recycling processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Concentration Dependent Interfacial Chemistry of the NaOH(aq):Gibbsite Interface

Caustic conditions are often employed for dissolution of a wide variety of minerals, where ion sorption, surface diffusion, and interfacial organization impact surface reactivity. In the case of gibbsite, gamma-Al(OH)3, the chemistry at the NaOH(aq) interface is deeply intertwined with industrial processing of aluminum, including metal production and the disposition of Al-containing wastes. To date, little is known about the structure, speciation, and dynamic behavior of gibbsite interfaces (and that of many other minerals) with NaOH(aq)—particularly as a function of ionic strength. Yet concentration-dependent interfacial organization and dynamics are a critical starting point to develop a fundamental understanding of the factors that influence dissolution. This work reports equilibrium molecular dynamics simulations of the gamma-Al(OH)3:NaOH(aq) interface, revealing the sorption behavior and speciation of ions from 0.5–10 M [NaOH]. As inner-sphere complexes, Na+ primarily coordinates to the side of the gibbsite hexagonal cavities, while OH accepts hydrogen-bonding from the surface-OH groups. The mobility of inner-sphere Na+ and OH ions is significantly reduced due to a strong surface affinity in comparison to previous reports of NaCl, CaCl2, or BaCl2 electrolytes. At high [NaOH], contact ion pairing that is observed in the bulk solution is partially disrupted upon sorption to the gibbsite surface by the individual ion–surface interactions. The molecular-scale changes to surface speciation and competition between ion–surface vs. ion–ion interactions influence surface characterization of gibbsite and potential dissolution processes, providing a valuable baseline for starting conditions needed within future reactive molecular simulations.

alkaline solutions, IDREAM, gibbsite, molecular dy↗

Development and application of hybrid AIMD/cDFT simulations for atomic-to-mesoscale chemistry

Many important chemical processes involve reactivity and dynamics in complex solutions. Gaining a fundamental understanding of these reaction mechanisms is a challenging goal that requires advanced computational and experimental approaches. However, important techniques such as molecular simulation have limitations in terms of scales of time, length, and system complexity. Furthermore, among the currently available solvation models, there are very few designed to describe the interaction between the molecular scale and the mesoscale. To help address this challenge, here, we establish a novel hybrid approach that couples first-principles plane-wave density functional theory with classical density functional theory (cDFT). In this approach, a region of interest described by ab initio molecular dynamics (AIMD) interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. cDFT is a robust but efficient mesoscopic approach to accurate thermodynamics of bulk electrolyte solutions over a wide concentration range (up to 2M concentrations). Benchmarking against commonly used continuum models of solvation, such as SMD, as well as experiments, demonstrates that our hybrid AIMD–cDFT method is able to produce reasonable solvation energies for a variety of molecules and ions. With this model, we also examined the solvent effects on a prototype S N 2 reaction of the nucleophilic attack of a chloride ion on methyl chloride in the solution. The resulting reaction pathway profile and the solution phase barrier agree well with experiment, showing that our AIMD/cDFT hybrid approach can provide insight into the specific role of the solvent on the reaction coordinate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rapid activation of non-oriented mechanophores via shock loading and spallation

Mechanophores (MPs), stimuli-responsive molecules that respond chromatically to mechanochemical reactions, are important for understanding the coupling between mechanics and chemistry as well as in engineering applications. However, the atomic-level understanding of their activation originates from gas phase studies or under simple linear elongation forces directly on molecules or polymer chains containing MPs. The effect of many-body distortions, pervasive in condensed-phase applications, is not understood. Therefore, we performed large-scale molecular dynamics (MD) simulations of a poly(methyl methacrylate)-spiropyran copolymer under dynamic mechanical loading and studied the activation of the MP under various conditions from dynamical compression to tension during unloading. Detailed analysis of the all-atom MD trajectories shows that the MP blocks experience significant many-body intramolecular distortion that can significantly decrease the activation barrier as compared with when deformation rates are slow relative to molecular relaxation time scales. Here we find that the reactivity of MPs under material compression states is governed by many-body effects of intramolecular torsions, whereas under tension, the reactions are governed by tensile stresses.

36 MATERIALS SCIENCE↗

Chlorine gas and anion radical reactivity in molten salts and the link to chlorobasicity

Next-generation nuclear power plants may include exciting novel designs in which molten salts are the coolant or a combination of the coolant and fuel. Whereas it is straightforward to see why having a low volatility coolant can be advantageous for safety, much is not understood about the production of volatile halogen gases as a result of radiation and even less is known about the distribution of these species at and away from interfaces. Using first principles molecular dynamics simulations, we investigate the product of the disproportionation reaction between chlorine anion radicals (nominally Cl 2 ˙ − ) in the bulk and slab configurations. We find that the product depends on the chlorobasicity of the medium. For example, in ZnCl 2 , Cl 2 forms, but in a eutectic mixture of LiCl and KCl, Cl 3 − is formed as a product. We also find that Cl 3 − prefers to form at the vapor interface and this may have implications for corrosion and reactivity. Furthermore, the mechanisms of the mobility of Cl 2 and Cl 3 − are radically different, the first one being vehicular and the second Grotthus-like. Chlorobasicity is linked to the electronic structure of the host melt; ZnCl 2 forms extended networks along which metal ions and anionic counterions have significant electronic orbital overlap forming long, linear, molecular-like constructs; the opposite is true for the alkali metal eutectic salt.

Nguyen, Hung H. [Univ. of Iowa, Iowa City, IA (Uni↗

VUV Photoionization Dynamics and Reactivity of Heterogeneous Water Clusters

This feature focuses on bridging isolated water and bulk water studies. It assembles experiments and theory on water clusters, mainly probed by vacuum ultraviolet (VUV) radiation, and summarizes what the effects are of intermolecular interactions on both the spectroscopy and the VUV-induced processes in water. In particular, it highlights studies of heterogeneous water clusters─those incorporating other molecular species such as naphthalene, methane, formic acid, and glycerol─which serve as model systems to investigate the fundamental roles of water in hydrogen bonding networks, proton transfer, and astrochemical processes. These mixed clusters provide a platform to revisit the stability of protonated water clusters and compile observations of ionization-induced structural rearrangements and fragmentation, especially in systems involving hydroxyl-rich cosolvents. This feature then explores energy transfer mechanisms in molecular clusters, following VUV photoexcitation. It concludes with potential future directions. First, we investigate excited-state dynamics in molecular clusters through direct probing, complementing the synchrotron studies discussed here. Second, exploring gas-phase molecule evaporation from confined spaces and interfaces using advanced spectroscopic techniques sheds light on these ubiquitous, yet currently debated, molecular processes.

Cluster chemistry↗

Molecular Design of Functional Polymers for Silica Scale Inhibition

Silica polymerization, which involves the condensation reaction of silicic acid, is a fundamental process with wide-ranging implications in biological systems, material synthesis, and scale formation. The formation of a silica-based scale poses significant technological challenges to energy-efficient operations in various industrial processes, including heat exchangers and water treatment membranes. Despite the common strategy of applying functional polymers for inhibiting silica polymerization, the underlying mechanisms of inhibition remain elusive. In this study, we synthesized a series of nitrogen-containing polymers as silica inhibitors and elucidated the role of their molecular structures in stabilizing silicic acids. Polymers with both charged amine and uncharged amide groups in their backbones exhibit superior inhibition performance, retaining up to 430 ppm of reactive silica intact for 8 h under neutral pH conditions. In contrast, monomers of these amine/amide-containing polymers as well as polymers containing only amine or amide functionalities present insignificant inhibition. Molecular dynamics simulations reveal strong binding between the deprotonated silicic acid and a polymer when the amine groups in the polymer are protonated. Notably, an extended chain conformation of the polymer is crucial to prevent proximity between the interacting monomeric silica species, thereby facilitating effective silica inhibition. Furthermore, the hydrophobic nature of alkyl segments in polymer chains disrupts the hydration shell around the polymer, resulting in enhanced binding with ionized silicic acid precursors compared to monomers. Finally, our findings provide novel mechanistic insights into the stabilization of silicic acids with functional polymers, highlighting the molecular design principles of effective inhibitors for silica polymerization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Anisotropic strength behavior of single-crystal TATB

High-rate strength behavior plays an important role in the shock initiation of high explosives, with plastic deformation serving to localize heat into hot spots and as a mechanochemical means to enhance reactivity. Recent simulations predict that detonation-like shocks produce highly reactive nanoscale shear bands in the layered crystalline explosive TATB (1,3,5-triamino-2,4,6-trinitrobenzene), but the thresholds leading to this response are poorly understood. We utilize molecular dynamics to simulate the high-rate compressive stress–strain response of TATB, with a focus on understanding flow behavior. Here, the dependence of strength on pressure and loading axis (crystal orientation) is explored. The deformation mechanisms fall broadly into two categories, with compression along crystal layers activating a buckling/twinning mode and compression normal to the layers producing nanoscale shear bands. Despite the complexity of the underlying mechanisms, the crystal exhibits relatively straightforward stress–strain curves. Most of the crystal orientations studied show rapid strain softening following the onset of yielding, which settles to a steady flow state. Trajectories are analyzed using five metrics for local states and structural order, but most of these metrics yield similar distributions for these deformation mechanisms. On the other hand, a recently proposed measure of intramolecular strain energy is found to most cleanly distinguish between these mechanisms, while also providing a plausible connection with mechanochemically accelerated decomposition kinetics. Localization of intramolecular strain energy is found to depend strongly on crystal orientation and pressure.

36 MATERIALS SCIENCE↗

Lunar and Planetary Science XXXVI, Part 9

The following topics were discussed: Monitoring floods with NASA's ST6 autonomous spacecraft experiment; Dynamical cloud models constrained by high resolution spectroscopy of zodiacal light; The oxygen isotopic composition of the sun and implications for oxygen processing in molecular clouds; A nochian/hisperian hiatus and erosive reactivation of martian valley networks; Hard x-ray spectro-microscopy techniques; Thermoluminescence studies of carbonaceous chondrites, etc.

Source record↗

First Principles Based Reactive Atomistic Simulations to Understand the Effects of Molecular Hypervelocity Impact on Cassini's Ion and Neutral Mass Spectrometer

We report here on the predicted impact of species such as ice-water, CO2, CH4, and NH3, on oxidized titanium, as well as HC species on diamond surfaces. These simulations provide the dynamics of product distributions during and after a hypervelocity impact event, ionization fractions, and dissociation probabilities for the various species of interest as a function of impact velocity (energy). We are using these results to determine the relevance of the fragmentation process to Cassini INMS results, and to quantify its effects on the observed spectra.

hypervelocity collisions↗