Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reaction”

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 235 records · Page 13

Feasibility Study of Two Candidate Reaction Wheel/thruster Hybrid Control Architecture Designs for the Cassini Spacecraft

As the first spacecraft to achieve orbit at Saturn in 2004, Cassini has collected science data throughout its four-year prime mission (2004-08), and has since been approved for a first and second extended mission through 2017. Cassini carries a set of three "fixed" reaction wheels and a backup reaction wheel (reaction wheel #4) is mounted on top of an articulable platform. If necessary, this platform could be articulated to orient the backup reaction wheel with the degraded wheel. The reaction wheels are used primarily for attitude control when precise and stable pointing of a science instrument such as the narrow angle camera is required. In 2001-02, reaction wheel #3 exhibited signs of bearing cage instability. As a result, reaction wheel #4 was articulated to align with reaction wheel #3. Beginning in July 2003, Cassini was controlled using wheel #1, #2, and #4. From their first use in the spring of 2000 until today, reaction wheels #1 and #2 have accumulated more than3.5 billions revolutions each. As such, in spite of very carefully management of the wheel spin rates by the mission operation team, there are some observed increases in the drag torque of the wheels' bearings. Hence, the mission operations team must prepare for the contingency scenario in which the reaction wheel #1 (in addition to wheel #3) had degraded. In this hypothetical fault scenario, the two remaining reaction wheels (#2 and #4) will not be able to provide precise and stable three-axis control of the spacecraft. In this study, we evaluate the feasibility of controlling Cassini using the two remaining reaction wheels and four thrusters to meet the science pointing requirements for two key science operational modes: the Optical Remote Sensing and Downlink, Fields, Particles, & Waves operation modes. The performance (e.g., pointing control error, pointing stability, hydrazine consumption rate, etc.) of the hybrid controllers in both operations scenarios will be compared with those achieved using an all-thruster controller design. Strength and weaknesses of the hybrid control architecture are assessed quantitatively.

Cassini↗

Revealing the role of redox reaction selectivity and mass transfer in current–voltage predictions for ensembles of photocatalysts

Photocatalysts are conceptually simple reaction units where nanoscale semiconductors integrated with catalysts drive a pair of redox reactions on illumination. However, the proximity of reaction sites performing cathodic and anodic reactions poses dire challenges to realize large light-to-fuel conversion efficiencies. In this study, a powerful, yet straightforward, equivalent-circuit detail-balance modeling framework is developed and applied to evaluate the performance of photocatalytic systems featuring multiple light absorbers. Specifically, low bandgap iridium-doped strontium titanate is modeled as a Z-scheme photocatalyst to achieve desirable hydrogen evolution and iron-based redox shuttle oxidation reactions. Our model has unique capabilities to simulate competing redox reactions and address mass-transfer limitations. In a significant departure from state-of-the-art circuit models, our study develops tools to perform load-line analyses by incorporating a net electrochemical load curve that includes both desired and competing redox reactions. Consequently, reaction selectivity is predicted from equivalent circuit models for photocatalytic and photoelectrochemical systems. Our investigation into ensembles comprised of multiple, semi-transparent light absorbers reveals their potential to outperform a single, optically thick light absorber, particularly when operated under mass-transfer-limited conditions. However, this outcome hinges on minimizing mass-transfer rates of select redox species to prevent undesired reactions of hydrogen oxidation and/or redox shuttle reduction. Our findings demonstrate that reaction selectivity can be achieved by tuning asymmetry in redox species mass-transfer even with perfectly symmetric electrocatalytic charge-transfer coefficients. The influences of various kinetic, mass-transfer, and thermodynamic parameters are explored to offer crucial insights for synthesis of the next-generation of photocatalysts and selective coatings, and reactor designs.

25 ENERGY STORAGE↗

A multi-region approach for the analysis of porous materials and gas-solid reactions using USAXS-SAXS-WAXS: application to CaO carbonation

The reaction rate of gas-solid non-catalytic reactions is typically investigated using reactant conversion data over time and ex-situ measurements of the porous solid reactant textural properties; these data enable the experimental estimation of the initial intrinsic reaction rate, whereas the evolution of the textural properties and the reaction rate over time are evaluated theoretically using reaction models. In this work, a different methodology is presented, based on: a) in-situ time-resolved USAXS-SAXS-WAXS measurements of the solid sample, and b) a multi-region modeling approach; this methodology allows for the estimation of the textural properties and the intrinsic reaction rate at any time during the reaction. Referring to a reaction in which a solid product is obtained from a solid reactant and a gaseous reactant, the porous particle is described as consisting of several distinct regions with different microstructural properties, and both a gas-solid reaction model and a SAXS model are derived and applied to the carbonation of porous CaO. The proposed reaction model highlights the role of the inaccessible reactant and accurately predicts the solid reactant conversion versus time profile. The SAXS model accurately predicts the measured linear trends of the Porod invariant and of the pre-factor of the power law scattering profile versus CaO mass fraction; in the absence of chemical reactions, the proposed equations extend classical SAXS theory to porous materials containing macropores and nonporous solid phases in addition to standard nanoscale inhomogeneous regions.

CO2 solid sorbents↗

Concepts Relevant for the Kinetic Analysis of Reversible Reaction Systems

The net rate of a reversible chemical reaction is the difference between unidirectional rates of traversal along forward and reverse reaction paths. In a multistep reaction sequence, the forward and reverse trajectories, in general, are not the microscopic reverse of one another; rather, each unidirectional route is comprised of distinct rate-controlling steps, intermediates, and transition states. Consequently, traditional descriptors of rate (e.g., reaction orders) do not reflect intrinsic kinetic information but instead conflate unidirectional contributions determined by (i) the microscopic occurrence of forward/reverse reactions (i.e., unidirectional kinetics) and (ii) the reversibility of reaction (i.e., nonequilibrium thermodynamics). This review aims to provide a comprehensive resource of analytical and conceptual tools which deconvolute the contributions of reaction kinetics and thermodynamics to disambiguate unidirectional reaction trajectories and precisely identify rate- and reversibility-controlling molecular species and steps in reversible reaction systems. The extrication of mechanistic and kinetic information from bidirectional reactions is accomplished through equation-based formalisms (e.g., De Donder relations) grounded in principles of thermodynamics and interpreted in the context of theories of chemical kinetics developed in the past 25 years. Here, the aggregate of mathematical formalisms detailed herein is general to thermochemical and electrochemical reactions and encapsulates a diverse body of scientific literature encompassing chemical physics, thermodynamics, chemical kinetics, catalysis, and kinetic modeling.

42 ENGINEERING↗

Reactants, products, and transition states of elementary chemical reactions based on quantum chemistry

Reaction times, activation energies, branching ratios, yields, and many other quantitative attributes are important for precise organic syntheses and generating detailed reaction mechanisms. Often, it would be useful to be able to classify proposed reactions as fast or slow. However, quantitative chemical reaction data, especially for atom-mapped reactions, are difficult to find in existing databases. Therefore, we used automated potential energy surface exploration to generate 12,000 organic reactions involving H, C, N, and O atoms calculated at the ω B97X-D3/def2-TZVP quantum chemistry level. We report the results of geometry optimizations and frequency calculations for reactants, products, and transition states of all reactions. Additionally, we extracted atom-mapped reaction SMILES, activation energies, and enthalpies of reaction. We believe that this data will accelerate progress in automated methods for organic synthesis and reaction mechanism generation—for example, by enabling the development of novel machine learning models for quantitative reaction prediction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probing the Potential Energy Profile of the I + (H 2 O) 3 → HI + (H 2 O) 2 OH Forward and Reverse Reactions: High Level CCSD(T) Studies with Spin-Orbit Coupling Included

Three different pathways for the atomic iodine plus water trimer reaction I + (H 2 O) 3 → HI + (H 2 O) 2 OH were preliminarily examined by the DFT-MPW1K method. Related to previous predictions for the F/Cl/Br + (H 2 O) 3 reactions, three pathways for the I + (H 2 O) 3 reaction are linked in terms of geometry and energetics. To legitimize the results, the “gold standard” CCSD(T) method was employed to investigate the lowest-lying pathway with the correlation-consistent polarized valence basis set up to cc-pVQZ(-PP). According to the CCSD(T)/cc-pVQZ(-PP)//CCSD(T)/cc-pVTZ(-PP) results, the I + (H 2 O) 3 → HI + (H 2 O) 2 OH reaction is predicted to be endothermic by 47.0 kcal mol -1 . The submerged transition state is predicted to lie 43.7 kcal mol -1 above the separated reactants. The I···(H 2 O) 3 entrance complex lies below the separated reactants by 4.1 kcal mol -1 , and spin-orbit coupling has a significant impact on this dissociation energy. The HI···(H 2 O) 2 OH exit complex is bound by 4.3 kcal mol -1 in relation to the separated products. Compared with simpler I + (H 2 O) 2 and I + H 2 O reactions, the I + (H 2 O) 3 reaction is energetically between them in general. It is speculated that the reaction between the iodine atom and the larger water clusters may be energetically analogous to the I + (H 2 O) 3 reaction. The iodine reaction I + (H 2 O) 3 is connected with the analogous valence isoelectronic bromine/chlorine reactions Br/Cl + (H 2 O) 3 but much different from the F + (H 2 O) 3 reaction. Significant difference with other halogen systems, especially for barrier heights, are seen for the iodine systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inertia-induced mixing and reaction maximization in laminar porous media flows

Solute transport and biogeochemical reactions in porous and fractured media flows are controlled by mixing, as are subsurface engineering operations such as contaminant remediation, geothermal energy production, and carbon sequestration. Porous media flows are generally regarded as slow, so the effects of fluid inertia on mixing and reaction are typically ignored. Here, we demonstrate through microfluidic experiments and numerical simulations of mixing-induced reaction that inertial recirculating flows readily emerge in laminar porous media flows and dramatically alter mixing and reaction dynamics. An optimal Reynolds number that maximizes the reaction rate is observed for individual pore throats of different sizes. This reaction maximization is attributed to the effects of recirculation flows on reactant availability, mixing, and reaction completion, which depend on the topology of recirculation relative to the boundary of the reactants or mixing interface. Recirculation enhances mixing and reactant availability, but a further increase in flow velocity reduces the residence time in recirculation, leading to a decrease in reaction rate. The reaction maximization is also confirmed in a flow channel with grain inclusions and randomized porous media. Interestingly, the domain-wide reaction rate shows a dramatic increase with increasing Re in the randomized porous media case. This is because fluid inertia induces complex three-dimensional flows in randomized porous media, which significantly increases transverse spreading and mixing. This study shows how inertial flows control reaction dynamics at the pore scale and beyond, thus having major implications for a wide range of environmental systems.

Chen, Michael A. (ORCID:0000000294614710)↗

Improved Performance of Cu(InGa)(SeS) 2 PV Modules Using the Reaction of Metal Precursors. Final Report

This project “Improved Performance of Cu(InGa)(SeS) 2 PV Modules using the Reaction of Metal Precursors” was a partnership led by the Institute of Energy Conversion (IEC) at the University of Delaware with Columbia University and the Molecular Foundry at the Lawrence Berkeley National Laboratory. The aim was to develop pathways to improve Cu(InGa)(SeS) 2 (CIGSS) thin film photovoltaic modules using processes compatible with low manufacturing cost. The CIGSS approach investigated was a two-step process including deposition of metal precursor films following by reaction in hydride gases utilizing IEC’s novel reactor. The process was similar to that under commercial development by the project’s industry partner Stion. When Stion went out of business mid-project the focus changed to a rapid thermal process considered more commercially viable. Approaches to improve the performance of solar cells using the reacted films focused on two material innovations. First, the overall Ga content was increased to increase the operating voltage, which is desirable for scale-up to commercial modules. Second, the processing and performance advantages arising from Ag alloying were investigated. Advanced characterization guided process and material development including control of relative composition gradients. Research on the formation of Cu-Ga-In metal precursors utilized sputtering deposition which is normally used in commercial applications. The work resulted in processes for deposition of precursor stacks with increased relative Ga content and effects of deposition parameters on morphology and phase composition were established. It was shown that the metal precursor films have comparable phase composition and morphology so subsequent reaction follows from the same starting point. The addition of Ag to the metal precursors gave more uniform morphology and improved adhesion of reacted films which enable higher reaction temperature for faster processing. A novel outcome was the discovery of a previously undocumented material phase in sputter-deposited and evaporated Ag-Cu-In-Ga thin films. Hydride gas reaction processes including time-temperature-concentration profiles were developed for different precursor compositions. This enables control of composition profiles to engineer through-film gradients for solar cell optimization with characterization and simulations used to correlate measured film composition profiles to measurements of devices. In particular, the gradient of sulfur at the front of the CIGSS film was found to be critical. The simulations guided process development leading to improved reproducibility of devices improved performance with higher Ga content and higher voltage. With Ag-alloyed precursors, the reaction pathways leading were determined. A significant finding was that Ag-alloying increases the reaction rate to completely convert precursor films to the final chalcopyrite which could enable reduced reaction time to benefit manufacturability. To maintain potential commercial viability, the process under investigation was refocused to a rapid thermal process that could potentially be incorporated into an in-line process for manufacturing. Precursors with different composition were capped with an extra selenium layer and reacted in hydrogen sulfide 5-15 minutes, compared to typically 2 hours in the previous multi-step batch process. Critical RTP parameters were identified to control the reaction. Further optimization would be needed for high efficiency solar cells but pathways to high quality devices with further optimization and improved heating uniformity were developed. The project also developed new optoelectronic characterization approaches with a focus on development and application of spatial- and time-resolved photoluminescence and a custom mapping photoluminescence microscope built. It was shown how critical electronic transport properties strongly depend on the chemical composition of the material and that a wide range of samples show inhomogeneity on a length scale larger than the grains in the films. Additionally, two-photon excitation capability was developed to distinguish bulk vs surface losses. The project advances the state-of-the -art for precursor reaction processes in several ways that could impact manufacturing. This includes validation of approaches to increase voltage and establishment of model-guided control to form optimal composition profiles. The application of process control approaches with knowledge of phase formation and reaction pathways can be critically valuable in designing a large-scale process.

14 SOLAR ENERGY↗

Molecular Mass Growth Processes to Polycyclic Aromatic Hydrocarbons through Radical–Radical Reactions Exploiting Photoionization Reflectron Time-of-Flight Mass Spectrometry

Polycyclic aromatic hydrocarbons (PAHs) represent critical building blocks in molecular mass growth processes to carbonaceous nanoparticles, referred to as interstellar and circumstellar grains along with soot particles in astrophysical environments and combustion systems, respectively. Recent advancements on elucidating elementary steps to PAHs have utilized reactions of aromatic radicals, resonantly stabilized free radicals, and aliphatic radicals with closed shell hydrocarbons. However, the role of radical–radical reactions (RRRs) leading to PAHs has remained largely unexplored on the molecular level due to preceding experimental challenges in producing sufficiently high number densities of radical reactants for isomer-selective detection of products from bimolecular and termolecular reactions. This Account offers the latest developments in our knowledge on the mechanisms and pathways to PAHs via RRRs probed in a chemical microreactor at temperatures as high as 1600 K. Product preservation in a molecular beam coupled with synchrotron vacuum ultraviolet photoionization reflectron time-of-flight mass spectrometry and photoelectron photoion coincidence spectroscopy enabled isomer-selective detection of PAHs of up to three rings by their photoionization efficiency curves, which were fit with a linear combination of reference curves for identification. Experiments were combined with computational fluid dynamics modeling of the physicochemical processes in the microreactor, as well as high-level electronic structure calculations to reveal the reaction pathways of each system. Six distinct reaction mechanisms were discovered in this work: propargyl addition─benzannulation (PABA), methyl addition─ring expansion (MARE), cyclopentadienyl addition─naphthylization (CPAN), fulvenallenyl addition─cyclization─aromatization (FACA), benzyl addition─aromatization (BAA), and phenyl addition─pentacyclization (PAP). By systematically varying the number of carbon atoms in the radical reactants, molecular mass growth processes involving reactions between radicals with odd numbers of carbon atoms access aromatics carrying one, two, or three six-membered rings, whereas reactions between even- and odd-carbon-numbered radicals produce aromatics combining five- and six-membered rings. Our investigations reveal unconventional cycloadditions on excited state triplet surfaces, additions of radicals to low spin density carbon-centered radicals, spiroaromatic and fulvene-type intermediates, and highly strained bicyclic reaction intermediates, challenging current perceptions of PAH molecular mass growth processes. All of the listed mechanisms, except for FACA, feature endoergic reactions or barriers which lie above the separated reactants and therefore might be central to circumstellar environments of carbon-rich stars and planetary nebulae as their descendants, but they play no role in the gas phase of cold molecular clouds where temperatures as low as 10 K dominate. Altogether, this work provides detailed reaction mechanisms of PAH growth processes, advancing our knowledge of the chemistry of carbonaceous matter in the universe.

Addition reactions↗

Constraining the destruction rate of K 40 in stellar nucleosynthesis through the study of the Ar 40 ( p , n ) K 40 reaction

Background: We present that K 40 plays a significant role in the radiogenic heating of Earth-like exoplanets, which can affect the development of a habitable environment on their surfaces. The initial amount of K 40 in the interior of these planets depends on the composition of the interstellar clouds from which they formed. Within this context, nuclear reactions that regulate the production of K 40 during stellar evolution can play a critical role. Purpose: In this study, we constrain for the first time the astrophysical reaction rate of K 40 ( n , p ) Ar 40 , which is responsible for the destruction of K 40 during stellar nucleosynthesis. We provide to the nuclear physics community high-resolution data on the cross section and angular distribution of the Ar 40 ( p , n ) K 40 reaction. These are important to various applications involving Ar 40 . The associated reaction rate of the Ar 40 ( p , n ) K 40 process addresses a reaction rate gap in the Joint Institute for Nuclear Astrophysics REACLIB database in the region of intermediate-mass isotopes. Methods: We performed differential cross-section measurements on the Ar 40 ( p , n ) K 40 reaction, for six energies in the center-of-mass system between 3.2 and 4.0 MeV and various angles between 0 ° and 135 ° . The experiment took place at the Edwards Accelerator Laboratory at Ohio University using the beam swinger target location and a standard neutron time-of-flight technique. We extracted total and partial cross sections by integrating the double differential cross sections we measured. Results: The total and partial cross sections varied with energy due to the contribution from isobaric analog states and Ericson type fluctuations. The energy-averaged neutron angular distributions were symmetrical relative to 90 ° . Based on the experimental data, local transmission coefficients were extracted and were used to calculate the astrophysical reaction rates of Ar 40 ( p , n ) K 40 and K 40 ( n , p ) Ar 40 reactions. The new rates were found to vary significantly from the theoretical rates in the REACLIB library. We implemented the new rates in network calculations to study nucleosynthesis via the slow neutron capture process, and we found that the produced abundance of K 40 is reduced by up to 10% compared to calculations with the library rates. At the same time, the above result removes a significant portion of the previous theoretical uncertainty on the K 40 yields from stellar evolution calculations. Conclusions: Our results support a destruction rate of K 40 in massive stars via the K 40 ( n , p ) Ar 40 reaction that is larger compared to previous estimates. The rate of K 40 destruction via the K 40 ( n , p ) Ar 40 reaction now has a dramatically reduced uncertainty based on our measurement. Lastly, this result directly affects the predicted stellar yields of K 40 from nucleosynthesis, which is a critical input parameter for the galactic chemical evolution models that are currently employed for the study of significant properties of exoplanets.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Forward Modeling of Gamma Reaction History Signatures From Anticipated Deuterium-Tritium Filled MagLIF Implosions on Sandia’s Z-Machine

Nuclear reaction history measurements provide a bang time and burn width of Inertial Confinement Fusion (ICF) implosions and are essential for understanding implosion performance to constrain ICF capsule design. When fusion fuel contains Deuterium (D) and Tritium (T) gas, reaction history is informed by measuring the 16.75 MeV gamma rays generated from the D(T,γ) 5 He reaction. Such DT based reaction history measurements have not been made on the Magnetized Laser Inertial Fusion (MagLIF) platform on Sandia’s Z-Machine due to the lack of Tritium being used. The recent development of ICF implosions with tritiated fuel will open the possibility of measuring the gamma reaction history on the Z-Machine. A forward model of the Gamma Reaction History diagnostic on Z (GRH-Z) has been developed using the MCNP6.3 (Monte-Carlo N-Particle) radiation transport code. The model included the Z-Machine geometry of interest to characterize the impact of neutron induced gamma rays on the DT reaction history measurements. In addition, the impulse response functions of the GRH-Z diagnostic to understand the temporal response of the detector, and the minimum yields required to make a reaction history measurement were calculated. This approach also predicted that with T 2 gas doping of MagLIF implosions a reaction history may be made for high performance shots >8e12-2.4e13 depending on the chosen threshold for the detector, with a maximum signal to background ratio of 25%. It was found that for long duration ICF implosions that additional collimation will be needed to prevent the neutron induced gamma rays from modifying the shape of the measured DT reaction history curve.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Chemical-damage MINE: A database of curated and predicted spontaneous metabolic reactions

Spontaneous reactions between metabolites are often neglected in favor of emphasizing enzyme-catalyzed chemistry because spontaneous reaction rates are assumed to be insignificant under physiological conditions. However, synthetic biology and engineering efforts can raise natural metabolites' levels or introduce unnatural ones, so that previously innocuous or nonexistent spontaneous reactions become an issue. Problems arise when spontaneous reaction rates exceed the capacity of a platform organism to dispose of toxic or chemically active reaction products. While various reliable sources list competing or toxic enzymatic pathways' side-reactions, no corresponding compilation of spontaneous side-reactions exists, nor is it possible to predict their occurrence. We addressed this deficiency by creating the Chemical Damage (CD)-MINE resource. First, we used literature data to construct a comprehensive database of metabolite reactions that occur spontaneously in physiological conditions. We then leveraged this data to construct 148 reaction rules describing the known spontaneous chemistry in a substrate-generic way. We applied these rules to all compounds in the ModelSEED database, predicting 180,891 spontaneous reactions. The resulting (CD)-MINE is available at https://minedatabase.mcs.anl. gov/cdmine/#/home and through developer tools. We also demonstrate how damage-prone intermediates and end products are widely distributed among metabolic pathways, and how predicting spontaneous chemical damage helps rationalize toxicity and carbon loss using examples from published pathways to commercial products. We explain how analyzing damage-prone areas in metabolism helps design effective engineering strategies. Finally, we use the CD-MINE toolset to predict the formation of the novel damage product N-carbamoyl proline, and present mass spectrometric evidence for its presence in Escherichia coli.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The reaction between the bromine atom and the water trimer: high level theoretical studies

Three different reaction pathways are found for the reaction of bromine atom (Br) with the lowest-energy structure of the water trimer [uud-(H 2 O) 3 ], initially using the MPW1K-DFT method. The three bromine pathways have closely related geometries and energetics, analogous to those found for the fluorine and chlorine reactions. Here, the lowest-energy pathway of the Br + uud-(H 2 O) 3 reaction was further investigated using the “gold standard” CCSD(T) method and the correlation-consistent basis sets up to cc-pVQZ(-PP). Based on the CCSD(T)/cc-pVQZ(-PP)//CCSD(T)/cc-pVTZ(-PP) results, the Br + (H 2 O) 3 reaction is endothermic by 33.3 kcal mol -1 . The classical barrier height is 29.0 kcal mol -1 between the reactants and the exit complex, and there is no barrier for the reverse reaction. The Br···(H 2 O) 3 entrance complex is found to lie 4.7 kcal mol -1 below the separated reactants, and the HBr···(H 2 O) 2 OH exit complex is bound by 6.4 kcal mol -1 relative to the separated products. This potential energy profile is further corrected by the zero point energies and spin–orbit coupling effects. Structurally, the Br + (H 2 O) 3 stationary points can be derived from those of the simpler Br + (H 2 O) 2 reaction by judiciously appending a H 2 O molecule. The Br + (H 2 O) 3 potential energy profile is compared with the Br + (H 2 O) 2 and Br + H 2 O reactions, as well as to the valence isoelectronic Cl + (H 2 O) 3 and F + (H 2 O) 3 systems. It is reasonable to expect that the reactions between the bromine atom and larger water clusters would be similar to the Br + (H 2 O) 3 reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Correlated energy uncertainties in reaction rate calculations

Monte Carlo methods can be used to evaluate the uncertainty of a reaction rate that arises from many uncertain nuclear inputs. However, until now no attempt has been made to find the effect of correlated energy uncertainties in input resonance parameters. Our goal is to investigate the impact of correlated resonance energy uncertainties on reaction rates. Using a combination of numerical and Monte Carlo variation of resonance energies, the effect of correlations are investigated. Five reactions are considered: two fictional, illustrative cases and three reactions whose rates are of current interest. The effect of correlations in resonance energies depends on the specific reaction cross section and temperatures considered. When several resonances contribute equally to a reaction rate, and when they are located on either side of the Gamow peak, correlations between their energies dilute their effect on reaction rate uncertainties. If they are both located above or below the maximum of the Gamow peak, however, correlations between their resonance energies can increase the reaction rate uncertainties. This effect can be hard to predict for complex reactions with wide and narrow resonances contributing to the reaction rate.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurement of the 230 Th( p ,2n)Pa229 and 230 Th( p ,3n)Pa228 reaction cross sections from 14.1 to 16.9 MeV

Actinium-225 is of interest for medical isotope production and there is on-going research into methods of producing Ac 225 , either directly or via the decay of its parent isotopes ( Th 229 , Pa 229 , and Ra 225 ). One method that has been suggested is the Th 230 ( p , 2 n ) Pa 229 reaction. However, there is no available cross-section data for this reaction in the literature. Purpose: Measure the Th 230 ( p , 2 n ) and Th 230 ( p , 3 n ) reaction cross sections in the energy range where the ( p , 2 n ) reaction is predicted to peak to determine the feasibility of Ac 225 production via the Th 230 ( p , 2 n ) reaction. Methods: Targets naturally enriched in Th 230 were irradiated at the Center for Accelerator Mass Spectrometry at Lawrence Livermore National Laboratory with energies ranging from 14.1 to 16.9 MeV. Furthermore, chemical processing was used to separate the protactinium activation products, followed by γ -ray spectroscopy to measure the activities of Pa 228 , 229 , 230 , 232 produced in the irradiation. Results: We find that excitation functions are reported for the first time in the literature for the Th 230 ( p , 2 n ) and Th 230 ( p , 3 n ) reactions in this energy range. The peak measured value of the Th 230 ( p , 2 n ) reaction was found to be 182 ± 12 mb at 14.4 ± 0.1 MeV. The Th 232 ( p , n ) Pa 232 reaction was used to verify the experimental conditions, the measured values are reported and are comparable to the existing literature values. From the γ -ray spectrometry data, the half-life of Pa 229 was measured as 1.5 ± 0.1 days, which is within the error of the half-life reported in the evaluated nuclear data as well as in the recent measurements, and the half-life of Pa 228 was measured as 19.5 ± 0.4 hours. Conclusions: Overall, the Th 230 ( p , 2 n ) Pa 229 reaction could reasonably be used for Ac 225 isotope production, although significant amounts of relatively isotopically pure Th 230 would be needed for significant production because the low alpha-decay branching ratio of Pa 229 and long half-life of Th 229 inhibit the in-growth of significant amounts of Ac 225 .

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

First principles reaction discovery: from the Schrodinger equation to experimental prediction for methane pyrolysis

Our recent success in exploiting graphical processing units (GPUs) to accelerate quantum chemistry computations led to the development of the ab initio nanoreactor, a computational framework for automatic reaction discovery and kinetic model construction. In this work, we apply the ab initio nanoreactor to methane pyrolysis, from automatic reaction discovery to path refinement and kinetic modeling. Elementary reactions occurring during methane pyrolysis are revealed using GPU-accelerated ab initio molecular dynamics simulations. Subsequently, these reaction paths are refined at a higher level of theory with optimized reactant, product, and transition state geometries. Reaction rate coefficients are calculated by transition state theory based on the optimized reaction paths. The discovered reactions lead to a kinetic model with 53 species and 134 reactions, which is validated against experimental data and simulations using literature kinetic models. We highlight the advantage of leveraging local brute force and Monte Carlo sensitivity analysis approaches for efficient identification of important reactions. Both sensitivity approaches can further improve the accuracy of the methane pyrolysis kinetic model. The results in this work demonstrate the power of the ab initio nanoreactor framework for computationally affordable systematic reaction discovery and accurate kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Heterogeneous chemistry related to Antarctic ozone depletion: Reaction of ClONO2 and N2O5 on ice surfaces

Laboratory studies of heterogeneous reactions of possible importance for Antarctic ozone depletion were performed. In particular, the reactions of chlorine nitrate (ClONO2) and dinitrogen pentoxide (N2O5) were investigated on ice and HCl/ice surfaces. These reactions occur on the surfaces of polar stratospheric clouds (PSCs) over Antarctica. One reaction transforms the stable chlorine reservoir species (ClONO2 and HCl) into photochemically active chlorine in the form of HOCl and Cl2. Condensation of HNO3 in the reactions removes odd nitrogen from the stratosphere, a requirement in nearly all models of Antarctic ozone depletion. Other reactions may also be important for Antarctic ozone depletion. Like the reactions of chlorine nitrate, these reactions deplete odd nitrogen through HNO3 condensation. In addition, one reaction converts a stable chlorine reservior species (HCl) into photochemically active chlorine (ClNO2). These reactions were studied with a modified version of a Knudsen cell flow reactor.

Tolbert, Margaret A.↗

Temperature-Dependent Kinetics Studies of the Reactions Br((sup 2)P3/2) + H2S yields SH + HBr and Br((sup 2)P3/2) + CH3SH yields CH3S + HBr. Heats of Formation of SH and CH3S Radicals

Time resolved resonance fluorescence detection of Br(sup 2)P3/2) atom disappearance or appearance following 266-nm laser flash photolysis of CF2Br2/H2S/H2/N2, CF2Br2/CH3SH/H2/N2, Cl2CO/H2S/HBr/N2, and CH3SSCH3/HBr/H2/N2 mixtures has been employed to study the kinetics of the reactions Br((sup 2)P3/2) + H2S = SH + HBr (1,-1) and Br((sup2)P3/2) + CH3SH = CH3S + HBr (2, -2) as a function of temperature over the range 273-431K. Arrhenius expressions in units of 10(exp -12) cu cm/molecule/s which describe the results are k1 = (14.2 +/- 3.4) exp[(-2752 +/- 90)/T],(k-1) = (4.40 +/- 0.92) exp[(-971 +/- 73)/T],k(2) = (9.24 +/- 1.15) exp[(-386 +/- 41)/T], and k(-2) = (1.46 +/-0.21) exp[(-399 +/-41)/T; errors are 2 sigma and represent precision only. By examining Br((sup 2)P3/2) equilibrium kinetics following 355nm laser flash photolysis of Br2/CH3SH/H2/N2 mixtures, a 298 K rate coefficient of (1.7 +/- 0.5) x 10(exp -10) cu cm/molecule/s has been obtained for the reaction CH3S + Br2 yields CH3SBr + Br. To our knowledge, these are the first kinetic data reported for each of the reactions studied. Measured rate coefficients, along with known rate coefficients for similar radical + H2S, CH3SH, HBr,Br2 reactions are considered in terms of possible correlations of reactivity with reaction thermochemistry and with IP - EA, the difference between the ionization potential of the electron donor and the electron affinity of the electron acceptor. Both thermochemical and charge-transfer effects appear to be important in controlling observed reactivities. Second and third law analyses of the equilibrium data for reactions 1 and 2 have been employed to obtain the following enthalpies of reaction in units of kcal/mol: for reaction 1, Delta-H(298) = 3.64 +/- 0.43 and Delta-H(0) = 3.26 +/-0.45; for reaction 2, Delta-H(298) = -0.14 +/- 0.28 and Delta-H(0) = -0.65 +/- 0.36. Combining the above enthalpies of reaction with the well-known heats of formation of Br, HBr, H2S, and CH3SH gives the following heats of formation for the RS radicals in units of kcal/mol: Delta-H(sub f)(sub 0)(SH) = 34.07 +/- 0.72, Delta-H(sub f)(sub 298)(SH) = 34.18 +/- 0.68, Delta-H(sub f)(sub 0)(CH3S) = 31.44 +/- 0.54, Delta-H(sub f)(sub 298)(CH3S) = 29.78 +/- 0.44; errors are 2 sigma and represent estimates of absolute accuracy. The SH heat of formation determined from our data agrees well with literature values but has reduced error limits compared to other available values. The CH3S heat of formation determined from our date is near the low end of the range of previous estimates and is 3-4 kcal/mol lower than values derived from recent molecular beam photofragmentation studies.

Nicovich, J. M.↗