Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic catalysis”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Structural Basis for Flip-Flop Action of Thiamin Pyrophosphate-dependent Enzymes Revealed by Human Pyruvate Dehydrogenase

The derivative of vitamin B1, thiamin pyrophosphate, is a cofactor of enzymes performing catalysis in pathways of energy production. In alpha (sub 2) beta (sub 2)-heterotetrameric human pyruvate dehydrogenase, this cofactor is used to cleave the C(sup alpha) -C(=O) bond of pyruvate followed by reductive acetyl transfer to lipoyl-dihydrolipoamide acetyltransferase. The dynamic nonequivalence of two, otherwise chemically equivalent, catalytic sites has not yet been understood. To understand the mechanism of action of this enzyme, we determined the crystal structure of the holo-form of human pyruvate dehydrogenase at 1.95-Angstrom resolution. We propose a model for the flip-flop action of this enzyme through a concerted approximately 2-Angstrom shuttle-like motion of its heterodimers. Similarity of thiamin pyrophosphate binding in human pyruvate dehydrogenase with functionally related enzymes suggests that this newly defined shuttle-like motion of domains is common to the family of thiamin pyrophosphate-dependent enzymes.

Ciszak, Ewa M.↗
DOE OSTI2024

Conducting polymer-based electrochemical sensors: Progress, challenges, and future perspectives

Conducting polymers are promising due to their unique properties, such as excellent electrical and optical properties, physical and chemical stability, high conductivity, and effective redox properties with high-temperature stability and biocompatibility. Due to these properties, conducting polymers are useful in diverse applications like sensors, batteries, oil industries, biosensors, biomedicines, catalysis, cancer treatment, etc. This review article aims to discuss the recent trends and analysis of conducting polymer-based electrochemical sensors in diverse areas with all required sensor characteristics, such as the derived limit of detection, utilized techniques for the sensing analysis and derived linear dynamic range with the stability of the sensors. Conducting polymers and their nanocomposites-based electrochemical sensors have demonstrated exceptional capabilities towards detecting various biomolecules, heavy metals, pesticides, and viruses like SARS-COV-2. Incorporation of redox mediators, use of conducting hydrogels, and molecular imprinting are promising strategies for better performance of the derived sensor. The article has demonstrated the existing challenges and limitations and provided solutions in the field. In the future, conducting polymers-based electrochemical sensors can be utilized in wearable sensors and integrated with IoT devices for better reach in real-time applications. They can also be made more accessible with precise control and data output by following specific methodologies. Utilizing green and sustainable conducting polymers can be crucial in advancing eco-friendly practices in the future. Conducting polymer-based electrochemical sensors has affectivity in neurochemical and pathogen sensing, which is essential for brain function and mental health.

42 ENGINEERING↗
DOE OSTI2025

Atomic-scale secondary electron imaging for heterogeneous catalysis research

Surface characterization at the atomic scale is essential for understanding the catalytic properties of supported metal nanoparticles. Secondary electron (SE) imaging in scanning transmission electron microscopy (STEM) provides three-dimensional surface topographic information, enabling the characterization of the size, morphology, and distribution of supported nanoparticles. Furthermore, real-time observation of catalyst materials in a gaseous environment would enhance the understanding of catalyst dynamics under operational conditions. Ongoing technical developments in SE-STEM and advancements in computational methods are expected to facilitate atomic-scale surface observations and enable more quantitative and statistical analyses. This progress will not only elucidate fundamental mechanisms at the atomic level but also provide comprehensive and universal insights into catalyst performances. Here, this minireview showcases the recent advancements and research findings in surface-sensitive SE imaging in STEM for the characterization of active catalyst materials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

In Situ Structural Determination of a Homogeneous Ruthenium Racemization Catalyst and Its Activated Intermediates Using X-Ray Absorption Spectroscopy

The activation process of a known Ru-catalyst, dicarbonyl(pentaphenylcyclopentadienyl)ruthenium chloride, has been studied in detail using time resolved in situ X-ray absorption spectroscopy. The data provide bond lengths of the species involved in the process as well as information about bond formation and bond breaking. On addition of potassium tert-butoxide, the catalyst is activated and an alkoxide complex is formed. The catalyst activation proceeds via a key acyl intermediate, which gives rise to a complete structural change in the coordination environment around the Ru atom. The rate of activation for the different catalysts was found to be highly dependent on the electronic properties of the cyclopentadienyl ligand. During catalytic racemization of 1-phenylethanol a fast-dynamic equilibrium was observed.

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

Second Coordination Sphere Effects in an Evolved Ru Complex Based on Highly Adaptable Ligand Results in Rapid Water Oxidation Catalysis

A new Ru complex containing the deprotonated 2,2':6',2"-terpyridine-6,6"-diphosphonic acid (H 4 tPa) and pyridine (py) of general formula [Ru II (H 3 tPa-κ-N 3 O)(py) 2 ] + , 2 + , has been prepared and thoroughly characterized by means of spectroscopic and electrochemical techniques, X-ray diffraction analysis, and density functional theory (DFT) calculations. Complex 2 + presents a dynamic behavior in the solution that involves the synchronous coordination and the decoordination of the dangling phosphonic groups of the tPa 4– ligand. Yet, at oxidation state IV, complex 2 + becomes seven coordinated with the two phosphonic groups now bonded to the metal center. Further, at this oxidation state at neutral and basic pH, the Ru complex undergoes the coordination of an exogenous OH – group from the solvent that leads to an intramolecular aromatic O atom insertion into the CH bond of one of the pyridyl groups, forming the corresponding phenoxo-phosphonate Ru complex [Ru III (tPaO-κ-N 2 O P O C )(py) 2 ] 2– , 4 2– , where tPaO 5– is the 3-(hydroxo-[2,2':6',2"-terpyridine]-6,6"-diyl)bis(phosphonate) ligand. This new in situ generated Ru complex, 4 2– , has been isolated and spectroscopically and electrochemically characterized. In addition, a crystal structure has been also obtained using single-crystal X-ray diffraction techniques. Complex 4 2– turns out to be an exceptional water oxidation catalyst achieving record maximum turnover frequencies (TOF max ) on the order of 16 000 s –1 . A mechanistic analysis complemented with DFT calculations has further been carried out, showing the critical role of intramolecular second coordination sphere effects exerted by the phosphonate groups in lowering the activation energy at the rate-determining step.

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

Nanoparticle Tracking to Probe Transport in Porous Media

From the granular and fractured subsurface environment to highly engineered polymer membranes used in pharmaceutical purification, porous materials are ubiquitous in nature and industrial applications. In particular, porous media are used extensively in processes including water treatment, pharmaceutical sterilization, food/beverage processing and heterogeneous catalysis, where hindered mass transport is either essential to the process or a necessary, but undesirable limitation. Unfortunately, there are currently no universal models capable of predicting mass transport based on a description of the porous material because real porous materials are complex, and because many coupled dynamic mechanisms (e.g. adsorption, steric effects, hydrodynamic effects, electrostatic interactions, etc.) give rise to the observed macroscopic transport phenomena. While classical techniques, like nuclear magnetic resonance and dynamic light scattering, provide useful information about mass transport in porous media at the ensemble level, they provide limited insight into the microscopic mechanisms that give rise to complex phenomena such as anomalous diffusion, hindered pore-space accessibility, and unexpected retention under flow, among many others. To address this issue, we have developed refractive index matching imaging systems, combined with single-particle tracking methods, allowing the direct visualization of single particle motion within a variety of porous materials. Here, we summarize our recent efforts to advance the understanding of nanoparticle transport in porous media using single-particle tracking methods in both fundamental and applied scenarios. First, we describe the basic principles for two-dimensional and three-dimensional single-particle tracking in porous materials. Then, we describe concrete examples of nanoparticle transport in porous materials from two perspectives: (1) understanding fundamental elementary particle transport processes in porous media, including pore accessibility and cavity escape, which limit transport in porous media; (2) facilitating applications in industrial processes, e.g. by understanding the mechanisms of particle fouling and remobilization in filtration membranes. Finally, we provide an outlook of opportunities associated with investigating other types of mass transport in confined environments using single-particle tracking methods, including electrophoretic and self-propelled motion.

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

Spin-State and Reorganization Energy Considerations for Metal-Centered Photoredox Catalysis

Transition-metal complexes featuring metal-centered excited states have recently emerged as mechanistically distinct platforms for selective photochemistry, including photoredox catalysis. Among these, Co(III) complexes have demonstrated productive photoinduced electron transfer via the 3 T 1 metal-centered state. In contrast, photoreactivity from the 5 T 2 metal-centered state in Fe(II) polypyridyl complexes remains limited. Building on our prior report concerning reactivity associated with the 5 T 2 state in [Fe(tren(py) 3 )] 2+ (tren(py) 3 = tris(2-pyridylmethyliminoethyl)-amine), we introduced stronger-field ligands in an effort to increase excited-state energies of Fe(II) polypyridyl complexes and enhance reactivity. Despite achieving nanosecond-scale excited-state lifetimes and favorable thermodynamic driving forces, no photoreactivity was observed. Reinvestigation of the observations previously reported for [Fe(tren(py) 3 )] 2+ revealed interactions between the metal complex and the substrate in their respective ground states that mimicked dynamic quenching of the chromophore, prompting a reassessment of mechanistic considerations inherent in leveraging reductive chemistry from the 5 T 2 excited state of Fe(II). Our analysis indicates that electron transfer from the 5 T 2 excited state of a low-spin d6 metal is subject to significant barriers both in terms of reorganization energies and spin conservation that undermines its ability to act as an electron donor for photoredox catalysis. In contrast, ligand fields that are sufficient to stabilize the 3 T 1 excited state have available to them numerous spin-allowed and, in certain cases, near-barrierless pathways to engage in excited-state electron transfer (both oxidative and reductive depending on the identity of the metal). These results highlight the critical role of spin-state changes and their associated reorganization energy requirements in metal-centered photoredox catalysis.

charge transfer↗
DOE OSTI2025

Photoredox Catalysis with Spin Magnetic Field Effects: A Scheme for Chiral Resolution

Quantities that break both mirror symmetry and time-reversal symmetry, such as the orbital angular momentum, are known to connect molecular chirality with an applied magnetic field. This concept has led to observations such as magneto-chiral dichroism and chirality-induced spin selectivity (CISS). However, being small, these effects often require additional amplification procedures such as flow chemistry to achieve bulk enantioseparation. In this work, we demonstrate how the magnetic field effect on photogenerated radical pairs, which also breaks time-reversal symmetry, can be harnessed for enantiopurification. Fundamental to this process is the collective decay of the singlet and triplet radical-pair states made possible by an applied magnetic field. Because opposite enantiomers exhibit spin-orbit coupling matrix elements of opposite signs, the singlet and triplet decay channels interfere constructively in one enantiomer. Meanwhile, molecules of the other enantiomer are funnelled into the first enantiomer through excited-state chirality inversion, achieving (dynamic kinetic) chiral resolution. Using an axially chiral binaphthyl derivative and a borane photosensitizer as prototype, we predict an appreciable enantiomeric excess (e.e.) of 90% to be possible at steady state, attained within hundreds of milliseconds when irradiated by a laser. Importantly, our analytical results showcase regimes of perfect enantioselectivity (100% e.e.), accessible by further chemical optimization of the photosensitizer for which general strategies are discussed. Altogether, this work illustrates a so-far untapped but powerful control knob for photoredox catalysis based on spin chemistry principles.

Magnetic properties↗
DOE OSTI2022

Insights into Sorption and Molecular Transport of Aqueous Glucose into Zeolite Nanopores

Liquid-phase heterogeneous catalysis using zeolites is important for biomass conversion to fuels and chemicals. There is a substantial body of work on gas-phase sorption in zeolites with different topologies; however, studies investigating the diffusion of complex molecules in liquid medium into zeolitic nanopores are scarce. Here, we present a molecular dynamics study to understand the sorption and diffusion of aqueous β-d-glucose into β-zeolite silicate at T = 395 K and P = 1 bar. Through 2-μs-long molecular dynamics trajectories, we reveal the role of the solvent, the kinetics of the pore filling, and the effect of the water model on these properties. We find that the glucose and water loading is a function of the initial glucose concentration. Although the glucose concentration increases monotonically with the initial glucose concentration, the water loading exhibits a nonmonotonic behavior. At the highest initial concentration (~20 wt %), we find that the equilibrium loading of glucose is approximately five molecules per unit cell and displays a weak dependence on the water model. Glucose molecules follow a single-file diffusion in the nanopores due to confinement. Here, the dynamics of glucose and water molecules slows significantly at the interface. The average residence time for glucose molecules is an order of magnitude larger than that in the bulk solution, while it is about twice as large for the water molecules. Our simulations reveal critical molecular details of the glucose molecule’s local environment inside the zeolite pore relevant to catalytic conversion of biomass to valuable chemicals.

09 BIOMASS FUELS↗
DOE OSTI2024

Active and Transfer Learning of High-Dimensional Neural Network Potentials for Transition Metals

Classical molecular dynamics (MD) simulations represent a very popular and powerful tool for materials modeling and design. The predictive power of MD hinges on the ability of the interatomic potential to capture the underlying physics and chemistry. There have been decades of seminal work on developing interatomic potentials, albeit with a focus predominantly on capturing the properties of bulk materials. Such physics-based models, while extensively deployed for predicting the dynamics and properties of nanoscale systems over the past two decades, tend to perform poorly in predicting nanoscale potential energy surfaces (PESs) when compared to high-fidelity first-principles calculations. These limitations stem from the lack of flexibility in such models, which rely on a predefined functional form. Machine learning (ML) models and approaches have emerged as a viable alternative to capture the diverse size-dependent cluster geometries, nanoscale dynamics, and the complex nanoscale PESs, without sacrificing the bulk properties. Here, in this study, we introduce an ML workflow that combines transfer and active learning strategies to develop high-dimensional neural networks (NNs) for capturing the cluster and bulk properties for several different transition metals with applications in catalysis, microelectronics, and energy storage, to name a few. Our NN first learns the bulk PES from the high-quality physics-based models in literature and subsequently augments this learning via retraining with a higher-fidelity first-principles training data set to concurrently capture both the nanoscale and bulk PES. Our workflow departs from status-quo in its ability to learn from a sparsely sampled data set that nonetheless covers a diverse range of cluster configurations from near-equilibrium to highly nonequilibrium as well as learning strategies that iteratively improve the fingerprinting depending on model fidelity. All the developed models are rigorously tested against an extensive first-principles data set of energies and forces of cluster configurations as well as several properties of bulk configurations for 10 different transition metals. Our approach is material agnostic and provides a methodology to transfer and build upon the learnings from decades of seminal work in molecular simulations on to a new generation of ML-trained potentials to accelerate materials discovery and design.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Deciphering electrocatalysts with multimodal operando approaches

Here, the electrocatalytic processes of a copper catalyst during nitrate electroreduction are unveiled by correlated operando microscopy and spectroscopy. Catalysts are vital to the modern chemical industry, yet their development has largely relied on trial-and-error approaches. Optimizing catalysts requires a fundamental understanding of their structure–chemical property behaviour under operational conditions. However, operando characterization remains challenging, especially for reactions occurring in the complex and dynamic liquid environments of electrochemical systems. Over the past decade, the rapid development of in situ and operando environmental transmission electron microscopy (ETEM) and microelectromechanical system (MEMS)-based closed-cell holders, enabling environmental studies within the vacuum environment of transmission electron microscopy (TEM), has substantively advanced understanding of heterogeneous catalysis, notably for gas-phase reactions. In situ ETEM enables the direct observation of the catalyst evolution in terms of structure, morphology and chemical state at the nano-to-atomic scale, providing insights into their correlation with catalytic performance.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Plant-wide modeling and techno-economic analysis of a direct non-oxidative methane dehydroaromatization process via conventional and microwave-assisted catalysis

Direct non-oxidative methane dehydroaromatization (DHA) process via conventional and microwave (MW)-assisted thermo-catalytic catalysis is studied. Rate models for methane DHA reactions, including the effect of catalyst deactivation, are developed by using the in-house experimental data. Model results for gas concentration profile and catalyst deactivation are in good agreement with the experimental data. This rate model is then used for the development of dynamic multi-scale, multi-physics commercial-scale reactor models. Total number of fixed bed reactors desired for a cyclic steady state process is estimated. Plant-wide models are then developed for conventional and MW-assisted processes for producing products of desired specifications. Techno-economic analysis of the methane DHA process is undertaken. Economics of these methane DHA processes are compared with the typical multi-step natural gas to aromatics production process via methanol synthesis. Sensitivity of internal rate of return (IRR) and net present value (NPV) to various economic and process parameters such as plant scale, desired rate of return, reactor cost, feedstock and utility cost, catalyst variable cost, and MW reactor cost is studied. Here, electric equivalent efficiency of the conventional methane DHA process is found to be 69.2 % and 67.3 % at 750 °C and 800 °C, respectively, while the MW-assisted methane DHA process has the electric equivalent efficiency of 48.9 % at 800 °C. IRRs of the conventional methane DHA process at 750 °C and 800 °C, and MW-assisted process are 15.2 %, 17.5 %, and 18.8 %, respectively for a methane feed flowrate of 19,782 kg/h, while the IRR of the multi-step natural gas to aromatics production process is estimated to be 0 % for the same plant scale. Impact of change in the methane price, electricity price, and catalyst cost is found to be considerable on the process economics, while the cost of the MW reactor is found to have negligible impact.

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

Dynamics of McMillan mappings III. Symmetric map with mixed nonlinearity

This article extends the study of the dynamical properties of the symmetric McMillan map, emphasizing its utility in understanding and modeling complex nonlinear systems. Although the map features six parameters, we demonstrate that only two are irreducible: the linearized rotation number at the fixed point and a nonlinear parameter representing the ratio of terms in the biquadratic invariant. Through a detailed analysis, we classify regimes of stable motion, provide exact solutions to the mapping equations, and derive a canonical set of action-angle variables, offering analytical expressions for the rotation number and nonlinear tune shift. We further establish connections between general standard-form mappings and the symmetric McMillan map, using the area-preserving Hénon map and accelerator lattices with thin sextupole magnet as representative case studies. Our results show that, despite being a second-order approximation, the symmetric McMillan map provides a highly accurate depiction of dynamics across a wide range of system parameters, demonstrating its practical relevance in both theoretical and applied contexts.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Enhanced Menshutkin S N 2 Reactivity in Mesoporous Silica: The Influence of Surface Catalysis and Confinement

A significant enhancement in the Menshutkin S N 2 reaction between 1-methylimidazole (MeIm) and methyl thiocyanate (MeSCN) is observed when the reaction is confined in the nanoscale silica pores of MCM41 and SBA15. The experiments in the silica pores are conducted without surrounding bulk reaction mixture. The influences of temperature, pore radius, and surface chemistry on the kinetics of the confined reaction are analyzed with time-dependent infrared spectroscopy, molecular dynamics simulations, and ab initio calculations. The rate constant of the pseudo-first order reaction increases with decreasing pore size, and the activation energy is found to decrease by 5.6 kJ/mol in the smallest pore studied (2.8 nm) relative to the bulk reaction. The rate constant dependence on pore size is accurately described by a two-state model in which molecules within the 4.6 Å interfacial layer experience a 2.4-fold rate constant increase relative to those reacting at the bulk rate further away from the interface. The removal of polar silanol groups from the silica surface via passivation with trimethylsilyl chloride results in bulk-like kinetics despite a reduction in the pore diameter, demonstrating the role of silanols as catalytic sites. Here, electronic structure calculations of the energy profile on a model silica surface confirm that silanol groups, particularly those of the vicinal type, can reduce the activation energy and reaction endothermicity through the donation of hydrogen bonds to the reactant, transition state, and product complexes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

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

Assessing entropy for catalytic processes at complex reactive interfaces

When chemical reactions are accelerated by a catalyst, entropy differences between reactants and their transient intermediates can be the driving force behind the promotion or inhibition of desired and parasitic chemical pathways. Understanding and controlling catalytic processes therefore requires both a fundamental and practicable understanding of entropy in addition to enthalpy. In unstructured media such as the vapor phase equilibrated with sparsely covered surfaces, entropy can be adequately accounted for by well-established approaches based on translational, rotational, and harmonic vibrational partition functions. However, these approximations become inadequate in more complex condensed phase environments, e.g., solid liquid interfaces of confined reaction spaces. In this chapter, we provide an overview of the state-of-art in the computational quantification of entropy and its known ramifications on catalysis. The fundamental roles of thermodynamics and kinetics in catalysis are covered in enough detail to appreciate and contextualize the computational methods employed to compute chemically accurate estimates of entropy. These methods are discussed in appropriate detail and range from the ubiquitous harmonic oscillator approximation where entropy unrelated to high frequency oscillations is typically underestimated, to enhanced free energy sampling with molecular dynamics where the desired accuracy must be weighed against the associated computational cost of obtaining it. The rising importance of machine learning and artificial intelligence in accelerating methodological progress in this field is touched upon, as well. Finally, applications, successes, and pitfalls of using these methods are provided to showcase past and present accomplishments while clarifying where improvements in both understanding and methodology are still needed.

Kollias, Loukas↗
DOE OSTI2024

THD-C Sheet: A Novel Nonbenzenoid Carbon Allotrope with Tetra-, Hexa-, and Dodeca-Membered Rings

Here, we propose a novel two-dimensional carbon-based structure with tetra-, hexa-, and dodeca-membered rings, which we refer to by the abbreviated name, THD-C. The structure presents a mixture of sp–sp 2 hybridization and can potentially be synthesized by the topological assembly of 4-ethynyldiphenylacetylene molecules. By employing first-principles calculations, the stability and ease of synthesis of the sheet are investigated and compared with various C-allotropes. We predict its metallic behavior and excellent kinetic and dynamic stability. Due to the crystal structure of the sheet, a strong mechanical anisotropy is observed. The effects of functionalization on the electronic properties of the material are also studied, and different semiconducting systems are obtained. The potential of THD-C for energy storage in metal-based batteries, hydrogen storage, and catalysis is also investigated, and we find a superior performance in comparison to graphite and other allotropes. The quantum confinement effect is investigated by constructing nanoribbons and nanotubes of various sizes. For ribbons, we find that tailor-made electronic and magnetic properties can be obtained and explored in potential spintronic devices. Additionally, we observe that nanotubes are conducting irrespective of their chirality and can potentially be used for capture, storage, and separation of industrially relevant small gas molecules.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗