SEARCH · Engineering Papers

Results for “chemicals”

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 613 records · Page 34

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Plasma Chemical Conversion and Resource Generation Beyond Low-Earth Orbit

As humanity gears up for its return to the moon after more than half a century, collaborative efforts between NASA, Artemis Accords Partners, and private industry are underway to establish the necessary infrastructure and technologies for lunar habitation and eventual Mars exploration. However, the traditional ISS resupply and waste management model is impractical and economically infeasible for prolonged missions to the Moon and Mars. Advanced chemical conversion technologies are needed to generate vital consumable products from local planetary resources (ISRU) and recycled gasses and waste within semi-closed loop life-support systems. Low-temperature plasma reactors are emerging power-to-gas technologies with the potential to facilitate various chemical synthesis processes with hardware commonality and redundancy. In plasma-based systems, electrical power is used to ionize a feedstock gas, creating a highly reactive environment that leverages electron excitation chemistry to break stable molecular bonds and form value added products. Unlike thermal chemical processes, plasma reactors operate at non-equilibrium conditions, allowing for lower-temperature operation and instantaneous start-up, making them adaptable to intermittent power availability. Moreover, their scalability permits deployment in both portable astronaut systems and large-scale industrial setups for colonies. One promising application of plasmas is for CO 2 conversion. Carbon dioxide comprises 96% of the Martian atmosphere and is a byproduct of human respiration, which typically must be scrubbed and vented from space habitats. A plasma source integrated with membrane separation technology could generate a stream of oxygen for life support and rocket propellant. CO 2 splitting may also be beneficial as a precursor to manufacture carbon-based products and fuels in situ, like methane, methanol, and polyethylene. Plasma-assisted CO 2 conversion is a simpler case to study without the concern of selectivity and is the first step toward complex chemical synthesis. This work presents preliminary experimental case study from a plasma reactor for CO 2 conversion and casts a vision for the potential of plasma technologies in a sub-architecture for resource production to enable the next generation of human spaceflight activities.

Plasma↗

DOE OSTI2022

The Importance of Chemical Research to the U.S. Economy

Chemistry plays a pivotal role in the strength of the U.S. economy and the advancement of humankind. Chemists' achievements include life-saving pharmaceuticals, advanced energy solutions, improved agricultural productivity, and novel materials used in products from clothing to electronic devices. The many sectors reliant on the U.S. chemical economy account for about 25% of the U.S. GDP and support 4.1 million U.S. jobs. However, a new and evolving chemistry landscape requires changes with regard to funding, training, and a focus on integrating sustainability into manufacturing, product usage, and product disposal. This report identifies strategies and options for research investments that will support U.S. leadership while considering environmental sustainability and developing a diverse chemical economy workforce with equitable opportunities for all chemistry talent. The report recommends that funding agencies and philanthropic organizations who support the chemical sciences fund as large a breadth of fundamental research projects as possible. Chemical industry and their partners at universities, scientific research institutions, and national laboratories should align the objectives of fundamental research to directly assist with new practices toward environmental stewardship, sustainability, and clean energy. Additionally, the report recommends that funding agencies make substantial investment toward education research to enable innovative ways of teaching about emerging concepts, tools and technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

ChemComp: Compiling and Computing with Chemical Reaction Networks

The exponential growth in computing demands driven by scientific computing, data analytics, and artificial intelligence is pushing conventional CMOS-based high-performance computing systems to their physical and energy efficiency limits. As we approach the era of post-exascale computing, disruptive approaches are necessary to overcome these barriers and achieve substantial gains in energy efficiency. Analog and hybrid digital-analog computing systems have emerged as promising alternatives, offering the potential for orders-of-magnitude improvements in efficiency. Among these, biochemical computing stands out as a novel paradigm capable of leveraging the natural efficiency of chemical reactions, which have shown promise in solving optimization problems by converging to steady states. By scaling up reaction networks or reaction vessel sizes, biochemical systems present an opportunity to meet the high-performance demands of modern computing tasks. Despite their promise, significant theoretical and practical challenges remain, particularly in formulating and mapping computational problems to chemical reaction networks (CRNs) and designing viable biochemical computing devices. This paper addresses these challenges by introducing new ideas to ChemComp, a compilation and emulation framework for chemical computation. This work describes the mechanisms through which solutions to ordinary differential equations (ODEs) that can be represented as CRN systems can be achieved. Furthermore, we explain the design principles of an ODE dialect implemented as a multi-level intermediate representation (MLIR) compiler extension that will be coupled with existing infrastructure. We demonstrate the potential of our framework through a case study emulating a simplified chemical reservoir computing device. This work establishes foundational tools and methodologies necessary to harness the computational power of chemistry, paving the way for the development of energy-efficient, high-performance computing systems tailored to contemporary and future computational needs.

Bohm Agostini, Nicolas↗

DOE OSTI2026

Chemical and Radiological Compatibility Testing of 3D Printed Materials

To dramatically increase the adaptability, performance, and safety of processes in support of the Defense Waste Processing Facility (DWPF), Savannah River National Lab (SRNL) plans to perform chemical and radiological compatibility testing on a wide variety of 3D printed materials of interest. The 3D printing process provides numerous strategic operational benefits such as rapid prototyping of complex designs and geometry specific to the needs of the nuclear waste disposition process, as well as on-demand rapid prototyping and iteration with materials that aren’t as accessible through traditional manufacturing methods. Reaction chemistry in simulated waste batches can be matched closely to its radioactive counterpart, but glass reactor vessels have limitations. Vessel geometry can play a big factor in mixing transport limitations, process chemistry, and degradation reaction kinetics. In addition, additive manufacturing allows for much more detailed vessel design than traditional alternatives. Waste processing techniques in DWPF also encounter extreme chemical environments including high pH, strong acids, abrasive slurries, and significant irradiation. To meet these challenges, a matrix of various polymer, ceramic, and metal additive manufacturing materials have been exposed to a suite of chemical environments of interest as well as radioactive dose (such as gamma radiation from 60 Co) to properly test their durability under these conditions. Mass change has been monitored over a period of up to a week in these conditions, as well as added characterization for surface modification through Scanning Electron Microscopy/Electron Dispersive X-ray analysis (SEM/EDX). Further chemical characterization has been monitored through Fourier-Transform InfraRed Spectroscopy (FTIR), with planned investigation via thermal and tensile strength degradation. While the direct product of this research is identification of material(s) that can withstand specific hazardous environments encountered by the mercury water wash tank in DWPF process simulation experiments, the reference base of materials will be used for many other nuclear processes in the pursuit of rapidly developed, cost-efficient, and highly specific devices for environmental remediation and much more.

Wilson, Nathan W. [Savannah River National Laborat↗

DOE OSTI2020

Automation of Active Space Selection for Multireference Methods via Machine Learning on Chemical Bond Dissociation

Predicting and understanding the chemical bond is one of the major challenges of computational quantum chemistry. Kohn–Sham density functional theory (KS-DFT) is the most common method, but approximate density functionals may not be able to describe systems where multiple electronic configurations are equally important. Multiconfigurational wave functions, on the other hand, can provide a detailed understanding of the electronic structures and chemical bonds of such systems. In the complete active space self-consistent field (CASSCF) method, one performs a full configuration interaction calculation in an active space consisting of active electrons and active orbitals. However, CASSCF and its variants require the selection of these active spaces. This choice is not black box; it requires significant experience and testing by the user, and thus active space methods are not considered particularly user-friendly and are employed only by a minority of quantum chemists. Our goal is to popularize these methods by making it easier to make good active space choices. Here, we present a machine learning protocol that performs an automated selection of active spaces for chemical bond dissociation calculations of main group diatomic molecules. The protocol shows high prediction performance for a given target system as long as a properly correlated system is chosen for training. Good active spaces are correctly predicted with a considerably better success rate than random guess (larger than 80% precision for most systems studied). Our automated machine learning protocol shows that a “black-box” mode is possible for facilitating and accelerating the large-scale calculations on multireference systems where single-reference methods such as KS-DFT cannot be applied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

The Role of the Solvent in the Condensed-Phase Dynamics and Identity of Chemical Bonds: The Case of the Sodium Dimer Cation in THF

When a solute molecule is placed in solution, is it acceptable to presume that its electronic structure is essentially the same as in the gas phase? In this paper, we address this question from a simulation perspective for the case of the sodium dimer cation (Na 2 + ) molecule in both liquid Ar and liquid tetrahydrofuran (THF). In previous work, we showed that when local specific interactions between a solute and solvent are energetically on the order of a hydrogen bond, the solvent can become part of the chemical identity of the solute. Here, using mixed quantum/classical molecular dynamics simulations, we see that for the Na 2 + molecule, solute--solvent interactions lead to two stable, chemically-distinct coordination states (Na(THF) 4 -Na(THF) 5 + and Na(THF) 5 -Na(THF) 5 + ) that are not only stable themselves as gas-phase molecules but that also have a completely new electronic structure with important implications for the excited-state photodissociation of this molecule in the condensed phase. Furthermore, we show through a set of comparative classical simulations that treating the solute's bonding electron explicitly quantum mechanically is necessary to understand both the ground-state dynamics and chemical identity of this simple diatomic molecule; even use of the quantum-derived potential of mean force is insufficient to describe the behavior of the molecule classically. Lastly, we calculate the results of a proposed transient hole-burning experiment that could be used to spectroscopically disentangle the presence of the different coordination states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Advances in imaging of chemically reacting flows

Many important chemically reacting systems are inherently multi-dimensional with spatial and temporal variations in the thermochemical state, which can be strongly coupled to interactions with transport processes. Fundamental insights into these systems require multi-dimensional measurements of the thermochemical state as well as fluid dynamics quantities. Laser-based imaging diagnostics provide spatially and temporally resolved measurements that help address this need. The state of the art in imaging diagnostics is continually progressing with the goal of attaining simultaneous multi-parameter measurements that capture transient processes, particularly those that lead to stochastic events, such as localized extinction in turbulent combustion. Development efforts in imaging diagnostics benefit from advances in laser and detector technology. This article provides a perspective on the progression of increasing dimensionality of laser-based imaging diagnostics and highlights the evolution from single-point measurements to 1D and 2D multi-parameter imaging and 3D high-speed imaging. This evolution is demonstrated using highlights of laser-based imaging techniques in combustion science research as an exemplar of a complex multi-dimensional chemically reacting system with chemistry–transport coupling. Imaging diagnostics impact basic research in other chemically reacting systems as well, such as measurements of near-surface gases in heterogeneous catalysis. The expanding dimensionality of imaging diagnostics leads to larger and more complex datasets that require increasingly demanding approaches to data analysis and provide opportunities for increased collaboration between experimental and computational researchers in tackling these challenges.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Techno-Economic Optimization of Advanced Energy Plants with Integrated Thermal, Mechanical, and Electro-Chemical Storage (Final Report)

The increasing use of renewable energy sources is leading to increased cycling of fossil-fueled power plants (FFPP) that are designed to operate at base-loaded conditions. Integrating energy storage facilities with the FFPPs can be helpful in reducing load-following operation of FFPPs. Decentralized deployment of energy storage facilities at the FFPP level has considerable advantages due to possibilities of smaller storage capacities, immediate benefits realized by the host power plant due to the increase in the efficiency, cleaner emission and higher plant life to name a few. Most importantly, deployment at the FFPP level can exploit the existing equipment items and facilities at the host power plant, thus reducing the CAPEX and reducing the storage capacity. However, realization of these benefits will critically depend on novel configuration/integration strategies with the least impact on the power plant operation and its configuration. Furthermore, dynamics of the entire integrated system including both the FFPP and the storage technologies must be taken into account to obtain the cost-optimal solution. With these motivations, the objective of this project was to complete a focused evaluation of decentralized deployments of energy storage facilities at the FFPP level. Promising thermal, chemical, mechanical, and electro-chemical storage technologies were evaluated with due consideration of their transient response to obtain various optimal system concepts that can minimize the levelized cost of storage. For thermal storage, cryogenic energy storage along with high temperature heat storage in molten salt as well as phase change material were evaluated. For mechanical storage, cryogenic air storage and pumped hydro storages were evaluated. For chemical storage, hydrogen storage was evaluated. For electrochemical storage, the team evaluated sodium sulfur, vanadium redox flow battery as well as Li-ion batteries. Integration of these storage technologies with the natural gas combined cycle (NGCC) and supercritical pulverized coal (SCPC) plants was considered. Six technologies, namely molten salt, cryogenics, compressed air, pumped hydro, H2 storage and Li-ion battery storage, were downselected based on their levelized cost of storage. It was observed that the ranking of the optimal storage technology can differ based on the host power plant technology even when same demand/supply/price profile for electricity are considered. It was also observed that as the variability in power demand varies, the ranking of the optimal storage technologies vary. However, for the same demand and LMP profile, top six optimal storage technologies for NGCC vs SCPC plants did not differ much even though LCOS for the same technology and optimal size of a given storage technology did differ. Detailed techno-economic assessment of these six technologies was undertaken.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2023

Magnetic and Impedance Analysis of Fe 2 O 3 Nanoparticles for Chemical Warfare Agent Sensing Applications

A dire need for real-time detection of toxic chemical compounds exists in both civilian and military spheres. In this paper, we demonstrate that inexpensive, commercially available Fe 2 O 3 nanoparticles are capable of selective sensing of chemical warfare agents (CWAs) using frequency-dependent impedance spectroscopy, with additional potential as an orthogonal magnetic sensor. X-ray magnetic circular dichroism analysis shows that Fe 2 O 3 nanoparticles possess moderately lowered moment upon exposure to 2-chloroethyl ethyl sulfide (2-CEES) and diisopropyl methylphosphonate (DIMP) and significantly lowered moment upon exposure to dimethyl methylphosphonate (DMMP) and dimethyl chlorophosphate (DMCP). Associated X-ray absorption spectra confirm a redox reaction in the Fe 2 O 3 nanoparticles due to CWA structural analog exposure, with differentiable energy-dependent features that suggest selective sensing is possible, given the correct method. Impedance spectroscopy performed on samples dosed with DMMP, DMCP, and tabun (GA, chemical warfare nerve agent) showed strong, differentiable, frequency-dependent responses. The frequency profiles provide unique “shift fingerprints” with which high specificity can be determined, even amongst similar analytes. The results suggest that frequency-dependent impedance fingerprinting using commercially available Fe 2 O 3 nanoparticles as a sensor material is a feasible route to selective detection.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Data for "Decompartmentalization of the yeast mitochondrial metabolism to improve chemical production in Issatchenkia orientalis "

Microbial production of chemicals may suffer from inadequate cofactor provision, a challenge further exacerbated in yeasts due to compartmentalized cofactor metabolism. Here, we perform cofactor engineering through the decompartmentalization of mitochondrial metabolism to improve succinic acid (SA) production in Issatchenkia orientalis . We localize the reducing equivalents of mitochondrial NADH to the cytosol through cytosolic expression of its pyruvate dehydrogenase (PDH) complex and couple a reductive tricarboxylic acid pathway with a glyoxylate shunt, partially bypassing an NADH-dependent malate dehydrogenase to conserve NADH. Cytosolic SA production reaches a titer of 104 g/L and a yield of 0.85 g/g glucose, surpassing the yield of 0.66 g/g glucose constrained by cytosolic NADH availability. Additionally, expressing cytosolic PDH, we expand our I. orientalis platform to enhance acetyl-CoA-derived citramalic acid and triacetic acid lactone production by 1.22- and 4.35-fold, respectively. Our work establishes I. orientalis as a versatile platform to produce markedly reduced and acetyl-CoA-derived chemicals.

bioproducts↗

DOE OSTI2020

Autonomously Responsive Membranes for Chemical Warfare Protection

Stimuli-responsive materials offer new opportunities to resolve long-standing material challenges and are rapidly gaining pivotal roles in diverse applications. For example, smart protective garments that rapidly transport water vapor and autonomously block chemical threats are expected to enable an effective new paradigm of adaptive personal protection. However, the incorporation of these seemingly incompatible properties into a single responsive system remains elusive. Herein, a bistable membrane that can rapidly, selectively, and reversibly transition from a highly breathable state in a safe environment to a chemically protective state when exposed to organophosphate threats such as sarin is demonstrated. Dynamic response to chemical stimuli is achieved through the physical collapse of an ultrathin copolymer layer on the membrane surface, which efficiently gates transport through membrane pores composed of single-walled carbon nanotubes (SWNTs). The adoption of nanometer-wide SWNTs for ultrafast moisture conduction enables a simultaneous boost in size-sieving selectivity and water-vapor permeability by decreasing nanotube diameter, thereby overcoming the breathability/protection trade-off that limits conventional membrane materials. Adaptive multifunctional membranes based on this platform greatly extend the active use of a protective garment and present exciting opportunities in many other areas including separation processes, sensing, and smart delivery.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Rational Design of A Chemical Bath Deposition Based Tin Oxide Electron Transport Layer for Perovskite Photovoltaics

Chemical bath deposition is widely used to deposit SnO x as an electron transport layer in perovskite solar cells (PSCs). The conventional recipe uses thioglycolic acid (TGA) to facilitate attachments of SnO x particles onto the substrate. However, nonvolatile TGA has been reported to harm the operational stability of PSCs. In this work, we introduced a volatile oxalic acid (OA) as an alternative to TGA. OA, a dicarboxylic acid, functions as a chemical linker for the nucleation and attachment of particles to the substrate in the chemical bath. Moreover, OA can be readily removed through thermal annealing followed by a mild H 2 O 2 treatment, as shown by FTIR measurements. Synergistically, the mild H 2 O 2 treatment selectively oxidizes the surface of the SnOx layer, minimizing nonradiative interface carrier recombination. EELS (electron-energy-loss-spectroscopy) confirms that the SnO x surface is dominated by Sn 4+ , while the bulk is a mixture of Sn 2+ and Sn 4+ . This rational design of a CBD SnO x layer leads to devices with T85~1,500h, a significant improvement over the TGA-based device with T80~250h. Our champion device reached a power conversion efficiency of 24.6%. This work offers a rationale for optimizing the complex parameter space of CBD SnO x to achieve efficient and stable PSCs.

14 SOLAR ENERGY↗

DOE OSTI2023

New Chemical Dopant and Counterion Mechanism for Organic Electrochemical Transistors and Organic Mixed Ionic–Electronic Conductors

Abstract Organic mixed ionic–electronic conductors (OMIECs) have varied performance requirements across a diverse application space. Chemically doping the OMIEC can be a simple, low‐cost approach for adapting performance metrics. However, complex challenges, such as identifying new dopant materials and elucidating design rules, inhibit its realization. Here, these challenges are approached by introducing a new n‐dopant, tetrabutylammonium hydroxide (TBA‐OH), and identifying a new design consideration underpinning its success. TBA‐OH behaves as both a chemical n‐dopant and morphology additive in donor acceptor co‐polymer naphthodithiophene diimide‐based polymer, which serves as an electron transporting material in organic electrochemical transistors (OECTs). The combined effects enhance OECT transconductance, charge carrier mobility, and volumetric capacitance, representative of the key metrics underpinning all OMIEC applications. Additionally, when the TBA + counterion adopts an “edge‐on” location relative to the polymer backbone, Coulombic interaction between the counterion and polaron is reduced, and polaron delocalization increases. This is the first time such mechanisms are identified in doped‐OECTs and doped‐OMIECs. The work herein therefore takes the first steps toward developing the design guidelines needed to realize chemical doping as a generic strategy for tailoring performance metrics in OECTs and OMIECs.

36 MATERIALS SCIENCE↗

DOE OSTI2020

CheKiPEUQ Intro 2: Harnessing Uncertainties from Data Sets, Bayesian Design of Experiments in Chemical Kinetics**

When choosing experimental conditions, Bayesian statistical tools can predict the experimental choices which will yield the highest information gain. Experimental choices could be temperature, pressure, reaction time, number of measurements, reactor volume, etc.. Three example analyses are presented here, each using the software Chemical Kinetics Parameter Estimation and Uncertainty Quantification (CheKiPEUQ). Information gain is a measure of reduction of uncertainty in a model's parameters. The three chemical system examples presented each illustrate Bayesian Design of Experiments using information gain. In the first chemical example, temperature selection impacts the information gain for the free energy of reaction in a two-component equilibrium reaction. In the second example, temperature and pressure are explored for a competitive adsorption Langmuir replacement reaction system. Finally, the third example is a catalytic membrane reactor which is a culmination of the previous examples. The catalytic membrane reactor has a complex and nonlinear response in the observables which is solved by numerical evaluation. In the three examples, the experimental conditions are treated as design variables for maximizing information gain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

On the Structural Transformation of Ni/BaH 2 During a N 2 -H 2 Chemical Looping Process for Ammonia Synthesis: A Joint In Situ Inelastic Neutron Scattering and First-Principles Simulation Study

The demand for decarbonizing the ammonia industry by using renewable energy has invoked increasing research interests into catalyst development for effective N 2 reduction under mild conditions. Hydride-based materials are among some of the emerging catalysts for ammonia synthesis at ambient pressure and low temperatures (< 673 K). A recent chemical looping process based on Ni/BaH 2 showed the most promise as it can realize ammonia production at a temperature as low as 373 K and under ambient pressure. However, the chemical transformation of the hydride catalyst at the molecular level remains unclear in this process. Here, we report detailed in situ neutron spectroscopy and diffraction investigations along with first-principles simulations on the structural transformation of Ni/BaH 2 during the nitridation and hydrogenation steps in the chemical looping process for ammonia synthesis. It was shown that a ball-milling process of the starting Ni/BaH 2 could significantly decrease the size of BaH 2 and increase the density of defects, thus potentially enhancing the reactivity of the hydride. The evolution from BaH 2 to barium imide (BaNH) was evidenced in the inelastic neutron scattering (INS) and neutron diffraction results during the N 2 reaction step. During the hydrogenation study, in addition to the recovery of BaH 2 , a possible intermediate species, N-deficient barium imide, was also detected. In comparing the N 2 and H 2 reaction steps, the neutron results indicate that the hydrogenation step appears more difficult than the nitridation step, confirming the facile N 2 fixation property of Ni/BaH 2 catalyst in ammonia synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Chemical redox of lithium-ion solid electroactive material in a packed bed flow reactor

Chemical redox reactions between solid electroactive materials and dissolved electroactive compounds are a necessary component of emerging technologies relevant to renewable energy including high energy density flow batteries and battery material recycling. In this work, an initial investigation of heterogeneous chemical redox for a packed bed reactor configuration between solid electroactive material and dissolved redox shuttles will be described. Experimental conditions including the height of the packed bed, the redox shuttle solution concentration, solution flow rate, and operating temperature were varied and their impact on the molar conversion of the solid electroactive material in the packed bed was quantified using electroanalytical techniques on the reactor effluent. The progression of the reaction and its dependence on the different variables explored will be discussed in the context of the limiting processes for porous packed aggregates of the solid electroactive material undergoing chemical oxidation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

A projection-based analytical Jacobian framework for chemical kinetics applications

A major challenge in simulating complex combustion systems with detailed chemical kinetic models is the cost of integrating the chemical source terms, often done using stiff ODE solvers that require frequent Jacobian evaluations. Using analytically derived Jacobian matrices instead of divided-difference-based numerical Jacobian approximations can significantly reduce the associated computational cost. However, ambiguities arise in the formulation of analytical Jacobians because the chemical state of the system, or state vector, can be expressed in multiple ways, involving variables that are typically not independent from one another. Here, in this work, the consequences of those ambiguities on practical calculations are characterized in detail, and a generalized, projection-based framework is proposed as a mitigation strategy. Performances are assessed in a series of test cases involving a variety of configurations and numerical solution approaches. Results show that with proper treatment, commonly used analytical Jacobian formulations can be considered as equivalent for practical purposes, thereby alleviating concerns that the state vector chosen to express the governing equations and corresponding analytical Jacobian may significantly impact the accuracy of the simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗