SEARCH · Search NASA

Results for “molecular driving forces”

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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Shuffling mode competition leads to directionally anisotropic mobility of faceted Σ11 boundaries in fcc metals

Faceted grain boundaries can migrate in interesting and unexpected ways. For example, faceted Σ11 <110> tilt grain boundaries were observed to exhibit mobility values that could be strongly dependent on the direction of migration. In order to understand whether this directionally-anisotropic mobility is a general phenomenon and to isolate mechanistic explanations for this behavior, molecular dynamics simulations of bicrystals evolved under an artificial driving force are used to study interface migration for a range of boundary plane inclination angles and temperatures in multiple face centered cubic metals (Al, Ni, and Cu). We find that directionally anisotropic mobility is active in a large fraction of these boundaries in Ni and Cu and should therefore impact the coarsening of polycrystalline materials. On the other hand, no such anisotropy is observed in any of the Al boundaries, showing that this behavior is material-dependent. Migration of the faceted boundaries is accomplished through transformation events at facet nodes and incommensurate boundary plane facets, which are termed shuffling modes. Three major shuffling modes have been identified, namely Shockley shuffling, slip plane shuffling, and disordered shuffling. A shift from the first two ordered modes to the third disordered mode is found to be responsible for reducing or removing directionally-anisotropic mobility, especially at the highest temperatures studied.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Neuromorphic heat transport effects in a molecular junction

Understanding energy transport at the nanoscale is an open and fundamental challenge in the molecular sciences with direct implications for the design of new electronics, computing devices, and materials. While nanoscale energy transport under steady-state conditions has been studied extensively, there is much less known about energy transport under time-dependent driving forces, particularly in the far-from-equilibrium regime. In this work, we use nonequilibrium molecular dynamics simulations and stochastic thermodynamics to investigate energy transport in a well-studied nanoscale system—a molecular junction—subjected to a time-periodic temperature gradient. The primary observation is that molecular junctions can exhibit heat transport hysteresis, a phenomenon in which the heat flux through a system depends not only on the instantaneous value of a time-dependent temperature bias but also on the temporal history of that bias. The presented findings illustrate that molecular junctions can exhibit the specific memory effect—heat transport hysteresis—that is essential for the design of thermal neuromorphic computers. This work elucidates a potential pathway toward the realization of such devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Biomimetic Metal-Free Hydride Donor Catalysts for CO 2 Reduction

The catalytic reduction of carbon dioxide to fuels and value-added chemicals is of significance for the development of carbon recycling technologies. One of the main challenges associated with catalytic CO 2 reduction is product selectivity: the formation of carbon monoxide, molecular hydrogen, formate, methanol, and other products occurs with similar thermodynamic driving forces, making it difficult to selectively reduce CO 2 to the target product. Significant scientific effort has been aimed at the development of catalysts that can suppress the undesired hydrogen evolution reaction and direct the reaction toward the selective formation of the desired products, which are easy to handle and store. Inspired by natural photosynthesis, where the CO 2 reduction is achieved using NADPH cofactors in the Calvin cycle, in this study we explore biomimetic metal-free hydride donors as catalysts for the selective reduction of CO 2 to formate. Here, we outline our recent findings on the thermodynamic and kinetic parameters that control the hydride transfer from metal-free hydrides to CO 2 . By experimentally measuring and theoretically calculating the thermodynamic hydricities of a range of metal-free hydride donors, we derive structural and electronic factors that affect their hydride-donating abilities. Two dominant factors that contribute to the stronger hydride donors are identified to be (i) the stabilization of the positive charge formed upon HT via aromatization or by the presence of electron-donating groups and (ii) the destabilization of hydride donors through the anomeric effect or in the presence of significant structural constrains in the hydride molecule. Hydride donors with appropriate thermodynamic hydricities were reacted with CO 2 , and the formation of the formate ion (the first reduction step in CO 2 reduction to methanol) was confirmed experimentally, providing an important proof of principle that organocatalytic CO 2 reduction is feasible. The kinetics of hydride transfer to CO 2 were found to be slow, and the sluggish kinetics were assigned in part to the large self-exchange reorganization energy associated with the organic hydrides in the DMSO solvent. Finally, we outline our approaches to the closure of the catalytic cycle via the electrochemical and photochemical regeneration of the hydride (R–H) from the conjugate hydride acceptors (R + ). We illustrate how proton-coupled electron transfer can be efficiently utilized not only to lower the electrochemical potential at which the hydride regeneration takes place but also to suppress the unwanted dimerization that neutral radical intermediates tend to undergo. Overall, this account provides a summary of important milestones achieved in organocatalytic CO 2 reduction and provides insights into the future research directions needed for the discovery of inexpensive catalysts for carbon recycling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1988

Electrohydrodynamic migration of charged droplets in an insulating fluid

The motion of charged, conducting droplets present in an insulating fluid medium is analyzed under the action of an electric field, in microgravity. Previous analyses of this problem have considered the Maxwell stresses as the only driving force. In the present study, arguments from macroscopic thermodynamics and the molecular theory of surface tension are used to show that the surface tension gradients can be induced due to the variation of the electric potential on the interface. In the limit of Reynolds numbers small compared to unity, the terminal velocity of migration of the droplet is calculated under the combined action of the Maxwell stresses and the surface tension gradients. The results show that there are no surface tension gradients (i.e., no electric potential variation at the interface) in a case that is due to the convection of the surface charges, surface tension gradients do exist and tend to reduce the terminal velocity of the droplet. The shape of the droplet altered by the motion was also calculated, when the deformations from the spherical shape are small.

Balasubramaniam, R.↗

NASA NTRS1988

Electrohydrodynamic migration of charged droplets in an insulating fluid

Balasubramaniam, R.↗

DOE OSTI2025

Microscopic insights into the solvation of polyethylene glycol chains in water: A machine learning potential approach

Polyethylene glycol (PEG) is a structurally simple, nontoxic, and water-soluble polymer widely utilized in medical and pharmaceutical applications. Notably, when a PEG chain is immersed in water, the surrounding water molecules play a key role in driving conformational changes of this macromolecule. In this study, we explore the solvation behavior of PEG under mechanical strain using molecular dynamics simulations, with an interatomic potential obtained from machine learning. Our focus is on the transition from the favored coil-like conformation to an extended one under external force. Through analyses of radial distribution functions, hydrogen bonding, and solvation dynamics, we uncover how mechanical stretching influences the local hydration environment. Furthermore, we disentangle the enthalpic and entropic contributions to the conformational stability of PEG in water. Surprisingly, our neural network potential model identifies dewetting of PEG C-atoms, and not water H-bonding with PEG O-atoms, as the main enthalpic driving force for the coiling of PEG in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Ion-dependent protein–surface interactions from intrinsic solvent response

Significance Hard–soft interfaces between inorganic surfaces and biomolecules promote self-assembly processes with broad implications in biogeochemistry, energy sciences, nanomedicine, and origins of life. Yet, detailed molecular-scale understanding of inorganic–biomolecule interactions and their dependence on solution conditions is missing. We present a theory for the initial stages of inorganic–biomolecule assembly based on the far-field response of water, using experimentally characterized interactions between muscovite surfaces and mica-binding proteins as model systems. Our work connects molecular details of the solution to assembly outcomes and suggests the initial driving forces for assembly are dominated by long-range, ion-specific interactions. The connections made between interfacial structure and long-range surface–biomolecule interactions provide insights toward a predictive understanding of biomolecular self-assembly on mineral surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2001

A Structural and Molecular Approach for the Study Biomarkers

Investigation of the nucleation and growth of crystals in both abiotic and biotic systems is critical to seemingly diverse disciplines of geology, biology, environmental science, and astrobiology. While there are abundant studies devoted to the determination of the structure and composition of inorganic crystals, as well as to the development of thermodynamic and kinetic models, it is only recently that research efforts have been directed towards understanding mineralization in biological systems (i.e., biomineralization). Biomineralization refers to the processes by which living organisms form inorganic solids. Studies of the processes of biomineralization under low temperature aqueous conditions have focused primarily on magnetite forming bacteria and shell forming marine organisms. Many of the biological building materials consist of inorganic minerals (calcium carbonate, calcium phosphate, silica or iron oxide) intricately combined with organic polymers (like proteins). More recently, efforts have been undertaken to explore the nature of biological activities in ancient rocks. In the absence of well-preserved microorganisms or genetic material required for the polmerase chain reaction (PCR) method in molecular phylogenetic studies, using biominerals as biomarkers offers an alternative approach for the recognition of biogenic activity in both terrestrial and extraterrestrial environments. The primary driving force in biomineralization is the interaction between organic and inorganic phases. Thus, the investigation of the ultrastructure and the nature of reactions at the molecular level occurring at the interface between inorganic and organic phases is essential to understanding the processes leading to the nucleation and growth of crystals. It is recognized that crystal surfaces can serve as the substrate for the organization of organic molecules that lead to the formation of polymers and other complex organic molecules, and in discussions of the origins of life, is referred to as organic synthesis on mineral surfaces. Furthermore, it is suggested that the interaction between mineral surfaces and simple organic molecules resulted in the formation of amino acids, RNA, and perhaps other more complex molecules such as proteins. On the other hand, in natural systems, it is recognized that functional groups on cell walls or membranes of microorganisms serve as sites of nucleation and crystallization. The precise replication of biominerals with controlled structure, morphology, size and texture is not confined to higher organisms as it also occurs in primitive prokaryotic cells such as magnetotactic bacteria and cyanobacteria. This suggests that the principal strategies of biomineralization were established early on in the evolutionary history of organisms. It is critical, therefore, to search for common mechanisms within diverse biological systems. One such common factor is the capability for organization and self-assembly. Organic macromolecules such as proteins and lipids can aggregate and polymerize forming membranes or extracellular matrix. At the organic-inorganic interface, several factors such as lattice geometry, polarity, stereochemistry and topography may act in concert to control nucleation and growth of crystals. Although several models have been proposed that discuss the significance of these factors for biomineralization, no comprehensive experimental data are available. In contrast to crystallization in exclusively inorganic systems, the kinetics of reaction and structural relationships between organic and inorganic phases in biominerals or biomimetic material is poorly understood. For example, it is not clear if the concept of epitactic growth (geometrical matching of unit cells at the interface of a secondary crystal growing on a primary crystal) applies to organic-inorganic systems. In contrast to inorganic templates that often have a smooth and rigid surface that promotes epitactic growth, biological substrates are usually rough and result in a large degree of mismatch. It is apparent that factors controlling the reaction at the crystal-matrix interface are strongly dependent upon the nature of the substrate. Therefore, characterization of the assembled organic surface and surface structure of the inorganic phase is crucial to understanding the processes of biomineralization. The focus of our research is the investigation of the processes leading to the nucleation and growth of crystals on both natural and synthetic systems through an interdisciplinary approach that integrates molecular biology, morphology and mineralogy using advanced preparation and analytical techniques. We have studied run-products, particularly magnetite, siderite and other carbonates, that resulted from extracellular biomineralization by extremophiles isolated from a variety of extreme environments ranging from permafrost to hydrothermal vent systems. The results of this study are critical to recognizing biomarkers in terrestrial and extraterrestrial environments.

Thomas-Keprta, Kathie↗

DOE OSTI2022

The molecular mechanism of load adaptation by branched actin networks

Branched actin networks are self-assembling molecular motors that move biological membranes and drive many important cellular processes, including phagocytosis, endocytosis, and pseudopod protrusion. When confronted with opposing forces, the growth rate of these networks slows and their density increases, but the stoichiometry of key components does not change. The molecular mechanisms governing this force response are not well understood, so we used single-molecule imaging and AFM cantilever deflection to measure how applied forces affect each step in branched actin network assembly. Although load forces are observed to increase the density of growing filaments, we find that they actually decrease the rate of filament nucleation due to inhibitory interactions between actin filament ends and nucleation promoting factors. The force-induced increase in network density turns out to result from an exponential drop in the rate constant that governs filament capping. The force dependence of filament capping matches that of filament elongation and can be explained by expanding Brownian Ratchet theory to cover both processes. We tested a key prediction of this expanded theory by measuring the force-dependent activity of engineered capping protein variants and found that increasing the size of the capping protein increases its sensitivity to applied forces. In summary, we find that Brownian Ratchets underlie not only the ability of growing actin filaments to generate force but also the ability of branched actin networks to adapt their architecture to changing loads.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Clustering Effects on the Structure of Ionomer Solutions: A Combined SANS and Simulations Study

Ionic assemblies, or clusters, determine the structure and dynamics of ionizable polymers and enable their many applications. Fundamental to attaining well-defined materials is controlling the balance between van der Waals interactions that govern the backbone behavior and the forces that drive the formation of ionic clusters. Here, using small angle neutron scattering and fully atomistic molecular dynamics simulations, the structure of a model ionomer, sulfonated polystyrene in toluene solutions, was investigated as the cluster cohesion was tweaked by the addition of ethanol. The static structure factor was measured by both techniques and correlated with the size of the ionic clusters as the polymer concentration was varied. The conjunction of SANS results and molecular insight from MD simulations enabled the determination of the structure in these inhomogeneous networks on multiple length scales. Additionally, we find that across the entire concentration range studied, a network driven by the formation of ionic clusters was formed, where the size of the clusters drives the inhomogeneity of these systems. Tweaking the ionic clusters through the addition of ethanol impacts the packing of the sulfonated groups, their shape, and their size distribution, which in turn, affects the structure of these networks.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Tailoring Interfaces for Enhanced Methanol Production from Photoelectrochemical CO 2 Reduction

Efficient and stable photoelectrochemical reduction of CO 2 into highly reduced liquid fuels remains a formidable challenge, which requires an innovative semiconductor/catalyst interface to tackle. In this study, we introduce a strategy involving the fabrication of a silicon micropillar array structure coated with a superhydrophobic fluorinated carbon layer for the photoelectrochemical conversion of CO 2 into methanol. The pillars increase the electrode surface area, improve catalyst loading and adhesion without compromising light absorption, and help confine gaseous intermediates near the catalyst surface. The superhydrophobic coating passivates parasitic side reactions and further enhances local accumulation of reaction intermediates. Upon one-electron reduction of the molecular catalyst, the semiconductor–catalyst interface changes from adaptive to buried junctions, providing a sufficient thermodynamic driving force for CO 2 reduction. These structures together create a unique microenvironment for effective reduction of CO 2 to methanol, leading to a remarkable Faradaic efficiency reaching 20% together with a partial current density of 3.4 mA cm –2 , surpassing the previous record based on planar silicon photoelectrodes by a notable factor of 17. Furthermore, this work demonstrates a new pathway for enhancing photoelectrocatalytic CO 2 reduction through meticulous interface and microenvironment tailoring and sets a benchmark for both Faradaic efficiency and current density in solar liquid fuel production.

14 SOLAR ENERGY↗

DOE OSTI2024

Thermodynamic driving forces in contact electrification between polymeric materials

Abstract Contact electrification, or contact charging, refers to the process of static charge accumulation after rubbing, or even simple touching, of two materials. Despite its relevance in static electricity, various natural phenomena, and numerous technologies, contact charging remains poorly understood. For insulating materials, even the species of charge carrier may be unknown, and the direction of charge-transfer lacks firm molecular-level explanation. Here, we use all-atom molecular dynamics simulations to investigate whether thermodynamics can explain contact charging between insulating polymers. Based on prior work suggesting that water-ions, such as hydronium and hydroxide ions, are potential charge carriers, we predict preferred directions of charge-transfer between polymer surfaces according to the free energy of water-ions within water droplets on such surfaces. Broad agreement between our predictions and experimental triboelectric series indicate that thermodynamically driven ion-transfer likely influences contact charging of polymers. Furthermore, simulation analyses reveal how specific interactions of water and water-ions proximate to the polymer-water interface explain observed trends. This study establishes relevance of thermodynamic driving forces in contact charging of insulators with new evidence informed by molecular-level interactions. These insights have direct implications for future mechanistic studies and applications of contact charging involving polymeric materials.

36 MATERIALS SCIENCE↗

NASA NTRS1990

Titan

Saturn's satellite Titan is the second-largest in the solar system. Its dense atmosphere is mostly molecular nitrogen with an admixture of methane, a surface pressure of 1.5 bars and a surface temperature of 94K. The fundamental driving force in the long-term evolution of Titan's atmosphere is the photolysis of methane in the stratosphere to form higher hydrocarbons and aerosols. The current rate of photolysis and undersaturation of methane in the lower troposphere suggests the presence of a massive ethane-methane-nitrogen ocean. The ocean evolves to a more ethane-rich state over geologic time, driving changes in the atmospheric thermal structure. An outstanding issue concerning Titan's earliest history is the origin of atmospheric nitrogen: was it introduced into Titan as molecular nitrogen or ammonia? Measurement of the argon-to-nitrogen ratio in the present atmosphere provides a diagnostic test of these competing hypotheses. Many of the questions raised by the Voyager encounters about Titan and its atmosphere can be adequately addressed only by an entry probe, such as that planned for the Cassini mission.

Lunine, Jonathan I.↗

DOE OSTI2020

Effect of Molecular Structure of Quinones and Carbon Electrode Surfaces on the Interfacial Electron Transfer Process

Quinones can undergo thermodynamically reversible proton-coupled electron transfer reactions and are being applied as electroactive compounds in aqueous organic batteries. However, the electrochemical reversibility of these compounds is affected not only by their molecular structure but also by the properties of a carbon-based electrode surface. This study combines experimental and theoretical approaches to understand this dependence. Specifically, we study the electron transfer kinetics of two synthesized quinone derivatives and two commercially available ones with a glassy carbon, a highly ordered pyrolytic graphite, and a high-edge-density graphite electrode (HEDGE). The electrochemical reversibility is notably improved on the HEDGE, which shows a higher density of defects and presents oxygenated functional groups at its surface. The electron transfer kinetics are controlled by adsorbed species onto the HEDGE. Molecular dynamics simulation and quantum mechanics calculations suggest defects with oxygen-containing functional groups, such as C-O and C$=$O, on HEDGE surfaces drive the interaction with the functional groups of the molecules, during physisorption from van der Waals forces. The presence of sulfonic acid side groups and a greater number of aromatic rings in the molecular structure may contribute to a higher stabilization of quinone derivatives on HEDGEs. We propose that high-performance carbon-based electrodes can be obtained without catalysts for organic batteries, by the engineering of carbon-based surfaces with edge-like defects and oxygenated functional groups.

25 ENERGY STORAGE↗

DOE OSTI2020

In Situ Oxidation Studies of High-Entropy Alloy Nanoparticles

Although high-entropy alloys (HEAs) have shown tremendous potential for elevated temperature, anticorrosion, and catalysis applications, little is known on how HEA materials behave under complex service environments. Herein, we studied the high-temperature oxidation behavior of Fe 0.28 Co 0.21 Ni 0.20 Cu 0.08 Pt 0.23 HEA nanoparticles (NPs) in an atmospheric pressure dry air environment by in situ gas-cell transmission electron microscopy. It is found that the oxidation of HEA NPs is governed by Kirkendall effects with logarithmic oxidation rates rather than parabolic as predicted by Wagner's theory. Further, the HEA NPs are found to oxidize at a significantly slower rate compared to monometallic NPs. The outward diffusion of transition metals and formation of disordered oxide layer are observed in real time and confirmed through analytical energy dispersive spectroscopy, and electron energy loss spectroscopy characterizations. Localized ordered lattices are identified in the oxide, suggesting the formation of Fe 2 O 3 , CoO, NiO, and CuO crystallites in an overall disordered matrix. Hybrid Monte Carlo and molecular dynamics simulations based on first-principles energies and forces support these findings and show that the oxidation drives surface segregation of Fe, Co, Ni, and Cu, while Pt stays in the core region. The present work offers key insights into how HEA NPs behave under high-temperature oxidizing environment and sheds light on future design of highly stable alloys under complex service conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Vacuum-assisted carbon molecular sieve membrane reactor for non-oxidative ethane dehydrogenation

Non-oxidative ethane dehydrogenation (EDH) is equilibrium-limited and endothermic. Selective hydrogen removal using a gas-permeable membrane within the EDH reaction zone can overcome the thermodynamic equilibrium, enabling higher ethane conversions. Employing vacuum as the permeation driving force, rather than a sweep gas, enhances the industrial viability of membrane reactors by eliminating additional post-reaction separation units. This study presents a membrane reactor that integrates H 2 -permeable carbon molecular sieve (CMS) hollow fiber membranes embedded in a fixed bed of cobalt in a dealuminated beta zeolite (Co@DeAl-BEA) catalyst, utilizing a vacuum to remove hydrogen efficiently. The CMS membrane exhibits high hydrogen permeance and an excellent H 2 /C 2 H 6 separation factor. The membrane reactor significantly enhanced the ethane conversion under reaction conditions comparable to those reported in the literature. A Langmuir-Hinshelwood kinetic rate expression was developed and incorporated into a one-dimensional steady-state reactor model. The experimentally validated model indicates that increasing the number of hollow fibers improves ethane conversion, although ethane loss to the permeate limits the benefit. The contact area between the catalyst and the membrane limits the reactor performance more than the catalytic throughput. Furthermore, we find that the location of the catalyst packing relative to the hollow fiber membranes influences ethane loss and conversion. Higher reactor pressures and inlet ethane flow rates improve space-time yield at the expense of lower ethane conversion. Increasing reactor temperature or packing length promotes both performance metrics. The EDH membrane reactor demonstrated durability over 200 h of continuous operation, maintaining record-low deactivation rates and high ethylene selectivity. Protocols for catalyst regeneration were developed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Tailoring multi-dimensional hierarchical self-assembly of metallacages through balancing non-covalent interactions

Despite advances in non-covalent interactions and complex self-assembly, precise control of multi-dimensional hierarchical self-assembly (HSA) from the molecular level to larger scales remains challenging. Herein, we designed and synthesized two rigid tetratopic ligands with multiple preset driving forces, and assembled them with Zn(II) ions to obtain metallacages SC6 and SC12 . Through balancing multiple non-covalent interactions, including metal–organic coordination, hydrophobic, and π–π interactions, SC6 can hierarchically self-assemble driven by a poor solvent, from one-dimensional nanowires to super-helical nanostructures via non-equilibrium self-assembly, progressing continuously with time and the increasing proportion of the poor solvent. However, for SC12 with enhanced hydrophobic interactions, two-dimensional monolayer nanogrids were formed by hierarchical self-assembly. Notably, these structures can be recycled back to primary metallacages through simple dissolution, highlighting their potential for efficient recycling and reuse. These results demonstrate that multi-dimensional hierarchical structures enable precise construction by balancing non-covalent interactions through a bottom-up self-assembly approach. This study provides deeper insight into the mechanisms of HSA and a promising strategy for the tailored creation of complex structures and sustainable porous materials.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Third-body stabilization of supercritical CO 2 in CO oxidation: development and application of a ReaxFF force field for the CO/O/CO 2 system

Supercritical CO 2 (scCO 2 ) plays a crucial role as a solvent in separation processes, advanced power cycles, and materials processing. Nonetheless, the atomistic comprehension of how the dense scCO 2 matrix influences the fundamental reaction of carbon monoxide (CO) is still insufficiently explored. Experimental studies and molecular dynamics (MD) simulations frequently fail to detect the highly reactive, transient intermediates, such as atomic oxygen (O), that drive these reactions. Here, to address this issue, we have developed a novel ReaxFF reactive force field for the CO 2 /CO/O system. The force field parameters were calibrated using density functional theory and second-order Møller-Plesset calculations to model CO 2 crystal properties, intermolecular interactions, bond dissociation curves, and reaction energy barriers. The force field reproduces the cohesive energy of the CO 2 crystal, the pressure characteristics of bulk scCO 2 , the equation-of-state behavior over a wide pressure–density range, the pressure dependence of the C–O bond length under compression, and the structural properties of liquid and scCO 2 , as documented by experiments, ab-initio MD, and prominent non-reactive models. The force field was subsequently applied to study the CO + O → CO 2 reaction. In a dilute environment, the reaction is inefficient as the newly formed CO 2 rapidly dissociates due to excess kinetic and potential energy acquired from the exothermic reaction. Conversely, in a dense scCO 2 environment, the surrounding matrix acts as an efficient third body, stabilizing the emerging CO2 product via molecular collisions. Statistical analysis confirms an average excess energy dissipation of 133.9 ± 3.6 kcal/mol over 112.4 ± 17.9 ps. Kinetic energy decomposition reveals that ∼ 92% of the excess kinetic energy is stored in internal (rotational and vibrational) degrees of freedom. This ReaxFF force field establishes a mechanistic foundation for third-body stabilization in dense reactive environments.

Chowdhury, Emdadul Haque [Pennsylvania State Univ.↗