Search NASA⌕ Search

SEARCH · Search NASA

Results for “astrochemistry”

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.

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The reaction NH 2 + CH 2 O: Kinetic measurements and computational studies

The reaction of amidogen with formaldehyde is relevant to astrochemistry and couples nitrogen and carbon chemistry in flames. The first measurements of the temperature dependence of the rate constant k were made, over 409–643 K, and are summarized as k = 7.1 × 10 −12 exp(−18.2 kJ mol −1 /RT) cm 3 molecule −1 s −1 with a confidence interval of ±20%. These results address a controversy over the presence of a significant barrier in the addition/elimination pathway leading to the production of formamide between theoretical models that include the zero-point energy of all modes at the transition state and a model that excludes transitional modes. The latter would lead to a negligible barrier, which is inconsistent with the experiments. Contrary to earlier claims, CCSDT(Q)-corrected energies combined with transition state theory yield quantitative accord with the measurements.

Astrochemistry↗
DOE OSTI2024

Broadband Rotational Spectroscopy in Uniform Supersonic Flows: Chirped Pulse/Uniform Flow for Reaction Dynamics and Low Temperature Kinetics

ConspectusThe study of gas-phase chemical reactions at very low temperatures first became possible with the development and implementation of the CRESU (French acronym for Reaction Kinetics in Uniform Supersonic Flows) technique. CRESU relies on a uniform supersonic flow produced by expansion of a gas through a Laval (convergent-divergent) nozzle to produce a wall-less reactor at temperatures from 10 to 200 K and densities of 1016-1018 cm-3 for the study of low temperature kinetics, with particular application to astrochemistry. In recent years, we have combined uniform flows with revolutionary advances in broadband rotational spectroscopy to yield an instrument that affords near-universal detection for novel applications in photodissociation, reaction dynamics, and kinetics. This combination of uniform supersonic flows with chirped-pulse Fourier-transform microwave spectroscopy (Chirped-Pulse/Uniform Flow, CPUF) permits detection of any species with a modest dipole moment, thermalized to the uniform temperature of the gas flow, with isomer, conformer, and vibrational state specificity. In addition, the use of broadband, high-resolution, and time-dependent (microsecond time scale) micro- and mm-wave spectroscopy makes it an ideal tool for characterizing both transient and stable molecules, as well as studying their spectroscopy and dynamics.In this Account, we review recent advances made using the CPUF technique, including studies of photodissociation, radical-radical reaction dynamics, and low temperature kinetics. These studies highlight both the strength of universal and multiplexed detection and the challenges of coupling it to a high-density collisional environment. Product branching and product evolution as a function of time have been measured for astrochemically relevant systems, relying on the detailed characterization of these flow conditions via experiments and fluid dynamics simulations. In the photodissociation of isoxazole, an unusual heterocyclic molecule with a very low-energy conical intersection, we have identified 7 products in 5 reaction channels and determined the product branching, pointing to both direct and indirect pathways. We have also approached the same system from separated NO and C3H3 reactants to explore a broader range of the potential energy surface, demonstrating the power of multichannel branching measurements for complex radical-radical reactions. We determined the product branching in the C3H2 isomers in the photodissociation of the propargyl radical and identified the importance of a hydrogen atom catalyzed isomerization to the lowest energy cyclic form. This then motivated a study of direct D-H exchange reaction in radicals, in which we demonstrate that it is an important and overlooked pathway for deuterium fractionation in astrochemical environments. Recently, we have shown the measurement of low temperature kinetics inside an extended Laval nozzle, after which a shock-free secondary expansion to low temperature and density affords an ideal environment for detection by rotational spectroscopy. These results highlight the power and potential of the CPUF approach, and future prospects will also be discussed in light of these developments.

Dias, Nureshan↗
DOE OSTI2025

An Unconventional Dark Radical Chemistry in Dense Molecular Clouds: Directed Gas-Phase Formation of Naphthyl Radicals

The synthetic pathways to aromatic molecules inside photon-shielded dense molecular clouds remain a fundamental, unsolved enigma in astrochemistry and astrophysics, with low-temperature molecular growth routes involving aromatic radicals, such as prototype bicyclic naphthyl (C 10 H 7 • ), implicated as key sources. Here, exploiting crossed molecular beam experiments augmented by electronic structure calculations, unexpected pathways are exposed leading to the gas-phase formation of 1- and 2-naphthyl via barrierless bimolecular reactions of atomic carbon (C) with indene (C 9 H 8 ) and of dicarbon (C 2 ) with styrene (C 8 H 8 ) accompanied by ring expansion and cyclization together with aromatization. These facile routes challenge conventional wisdom that aromatic radicals are formed in deep space solely via “bright” gas-phase photochemistry of their closed-shell polycyclic aromatic hydrocarbon (PAH) precursors. A hitherto disregarded “dark” aromatic radical chemistry with aromatic radicals synthesized via gas-phase reactions offers new concepts on the chemical evolution of the chemistry of dark molecular clouds eventually culminating in the rapid formation of aromatics, fullerenes, and carbonaceous nanostructures.

Aromatic compounds↗
DOE OSTI2024

Reaction dynamics of S( 3 P) with 1,3-butadiene and isoprene: crossed-beam scattering, low-temperature flow experiments, and high-level electronic structure calculations

Sulfur atoms serve as key players in diverse chemical processes, from astrochemistry at very low temperature to combustion at high temperature. Building upon our prior findings, showing cyclization to thiophenes following the reaction of ground-state sulfur atoms with dienes, we here extend this investigation to include many additional reaction products, guided by detailed theoretical predictions. The outcomes highlight the complex formation of products during intersystem crossing (ISC) to the singlet surfaces. Here, we employed crossed-beam velocity map imaging and high-level ab initio methods to explore the reaction of S( 3 P) with 1,3-butadiene and isoprene under single-collision conditions and in low-temperature flows. For the butadiene reaction, our experimental results show the formation of thiophene via H 2 loss, a 2H-thiophenyl radical through H loss, and thioketene through ethene loss at a slightly higher collision energy compared to previous observations. Complementary Chirped-Pulse Fourier-Transform mmWave spectroscopy (CP-FTmmW) measurements in a uniform flow confirmed the formation of thioketene in the reaction at 20 K. For the isoprene reaction, we observed analogous products along with the 2H-thiophenyl radical arising from methyl loss and C 3 H 4 S (loss of ethene or H 2 + acetylene). CP-FTmmW detected the formation of thioformaldehyde via loss of 1,3-butadiene, again in the 20 K flow. Coupled-cluster calculations on the pathways found by the automated kinetic workflow code KinBot support these findings and indicate ISC to the singlet surface, leading to the generation of various long-lived intermediates, including 5-membered heterocycles.

Lang, Jinxin↗
DOE OSTI2024

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

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

Infrared Spectrum of the Pyrene Anion in the CH Stretching Region

In this work, we report the infrared spectrum of pyrene anion, measured using messenger tagging with up to three Ar atoms. We assign the spectrum using density functional theory and vibrational perturbation theory. Additionally, we discuss our results in the context of computed and experimental spectra from the literature as well as recent observations from astronomical sources, addressing the question whether PAH anions could contribute to the strong infrared emission bands at 3.29 μm from carbon- rich regions of space.

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

Implications of high-pressure oxygen hydrates on radiolytic oxygen in Jovian icy moons

Various icy moons, such as Europa and Ganymede, have thin oxygen atmospheres and exhibit spectral features attributed to oxygen held in their surface ices. The oxygen forms from the radiolysis of water. The interiors of these bodies are subject to high pressures and it is not known how deep into icy moons oxygen-bearing ices can penetrate, or the structures formed by the oxygen–water system at high pressure. Here, we show that oxygen hydrates are stable to 2.6 GPa, allowing them to penetrate deep into icy moons, both above and below proposed sub-surface liquid-water oceans. Similarities between oxygen and hydrogen hydrates indicate potentially enhanced recombination rates, transforming them back into water and offering a resolution to the discrepancy between predicted and measured radiolysis rates. In addition to the low-pressure CS-II clathrate, our results find three high-pressure phases in the oxygen–water system: an ST clathrate, a C 0 hydrate, and a filled ice isomorphous with methane hydrate III. This shows a vast storage potential for molecular oxygen in icy moons and indicates that Europa could still be absorbing oxygen into its crustal ice.

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

Ultraviolet H 2 luminescence in molecular clouds induced by cosmic rays

Context. Galactic cosmic rays (CRs) play a crucial role in ionisation, dissociation, and excitation processes within dense cloud regions where UV radiation is absorbed by dust grains and gas species. CRs regulate the abundance of ions and radicals, leading to the formation of more and more complex molecular species, and determine the charge distribution on dust grains. A quantitative analysis of these effects is essential for understanding the dynamical and chemical evolution of star-forming regions. Aims. The CR-induced photon flux has a significant impact on the evolution of the dense molecular medium in its gas and dust components. This study evaluates the flux of UV photons generated by CRs to calculate the photon-induced dissociation and ionisation rates of a vast number of atomic and molecular species, as well as the integrated UV photon flux. Methods. To achieve these goals, we took advantage of recent developments in the determination of the spectra of secondary electrons, in the calculation of state-resolved excitation cross sections of H 2 by electron impact, and of photodissociation and photoionisation cross sections. Results. We calculated the H 2 level population of each rovibrational level of the X, B, C, B′, D, B″, D′, and a states. We then computed the UV photon spectrum of H 2 in its line and continuum components between 72 and 700 nm, with unprecedented accuracy, as a function of the CR spectrum incident on a molecular cloud, the H 2 column density, the isomeric H 2 composition, and the dust properties. The resulting photodissociation and photoionisation rates are, on average, lower than previous determinations by a factor of about 2, with deviations of up to a factor of 5 for the photodissociation of species such as AlH, C 2 H 2 , C 2 H 3 , C 3 H 3 , LiH, N 2 , NaCl, NaH, O 2 + , S2, SiH, l-C 4 , and l-C 5 H. A special focus is given to the photoionisation rates of H 2 , HF, and N 2 , as well as to the photodissociation of H 2 , which we find to be orders of magnitude higher than previous estimates. We give parameterisations for both the photorates and the integrated UV photon flux as a function of the CR ionisation rate, which implicitly depends on the H 2 column density, as well as the dust properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Modeling the Enceladus dust plume based on in situ measurements performed with the Cassini Cosmic Dust Analyzer

We analyzed data recorded by the Cosmic Dust Analyzer on board the Cassini spacecraft during Enceladus dust plume traversals. Our focus was on profiles of relative abundances of grains of different compositional types derived from mass spectra recorded with the Dust Analyzer subsystem during the Cassini flybys E5 and E17. The E5 profile, corresponding to a steep and fast traversal of the plume, has already been analyzed. In this paper, we included a second profile from the E17 flyby involving a nearly horizontal traversal of the south polar terrain at a significantly lower velocity. Additionally, we incorporated dust detection rates from the High Rate Detector subsystem during flybys E7 and E21. We derived grain size ranges in the different observational data sets and used these data to constrain parameters for a new dust plume model. This model was constructed using a mathematical description of dust ejection implemented in the software package DUDI. Further constraints included published velocities of gas ejection, positions of gas and dust jets, and the mass production rate of the plume. Our model employs two different types of sources: diffuse sources of dust ejected with a lower velocity and jets with a faster and more colimated emission. From our model, we derived dust mass production rates for different compositional grain types, amounting to at least 28 kg s –1 . Previously, salt-rich dust was believed to dominate the plume mass based on E5 data alone. The E17 profile shows a dominance of organic-enriched grains over the south polar terrain, a region not well constrained by E5 data. By including both E5 and E17 profiles, we find the salt-rich dust contribution to be at most 1% by mass. This revision also results from an improved understanding of grain masses of various compositional types that implies smaller sizes for salt-rich grains. Our new model can predict grain numbers and masses for future mission detectors during plume traversals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Chemical templates of the Central Molecular Zone

Context . The Central Molecular Zone (CMZ) of the Milky Way exhibits extreme conditions, including high gas densities, elevated temperatures, enhanced cosmic-ray ionization rates, and large-scale dynamics. This makes it a perfect laboratory for astrochemical studies. With large-scale molecular surveys revealing increasing chemical and physical complexity in the CMZ, it is essential to develop robust methods to decode the chemical information embedded in this extreme region. Aims . A key step to interpreting the molecular richness found in the CMZ is building chemical templates tailored to its diverse conditions. In particular, understanding how CMZ environments affect shock and protostellar chemistry is crucial. The combined impact of high ionization, elevated temperatures, and dense gas remains insufficiently explored for observable tracers. Methods . For this study, we utilized UCLCHEM , a gas-grain time-dependent chemical model, to link physical conditions with their corresponding molecular signatures and identify key tracers of temperature, density, ionization, and shock activity. To achieve this, we ran a grid of models of shocks and protostellar objects representative of typical CMZ conditions, focusing on 24 species, including complex organic molecules. Results . Shocked and protostellar environments show distinct evolutionary timescales (≲10 4 vs. ≳10 4 years); 300 K emerges as a key temperature threshold for chemical differentiation. We find that cosmic-ray ionization and temperature are the main drivers of chemical trends. HCO + , H 2 CO, and CH 3 SH trace ionization, while HCO, HCO + , CH 3 SH, CH 3 NCO, and HCOOCH 3 show consistent abundance contrasts between shocks and protostellar regions over similar temperature ranges. Conclusions . We characterized the behavior of 24 species in protostellar and shock-related environments. While our models underpredict some complex organics in shocks, they reproduce observed trends for most species, supporting scenarios involving a need for recurring shocks in Galactic Center clouds and enhanced ionization toward Sgr B2(N2). Future work should assess the role of shock recurrence and metallicity in shaping chemistry.

Galaxy: center↗
DOE OSTI2026

The stellar origins of 96 Zr excesses in presolar graphites from the Murchison meteorite

Context. Zirconium-96 is a stable isotope that can be synthesized under different neutron-rich nucleosynthetic conditions. Astrophysical models predict its production to occur in various stellar environments: from low-to-intermediate-mass asymptotic giant branch (AGB) stars to massive stars and core-collapse supernovae. Aims. Detections of 96Zr excesses, in combination with other isotopic measurements from presolar grains can provide unique constraints on its stellar origin. Presolar grains are microscopic particles found in primitive Solar System materials, which formed in stellar winds and supernova ejecta. The isotopic composition of each grain can provide us a snapshot of the nucleosynthetic processes that took place during the parent star’s lifetime. Methods. In this study, we measured the stable isotopes of C, N, O, Mo, Zr, and Ru in high-density presolar graphite grains from the Murchison meteorite and found four grains that contain positive isotopic anomalies in 96 Zr carried by their internal subgrains. We analyzed multi-element isotopic datasets from each grain to explore the source of the observed 96 Zr excesses. Results. Comparisons with stellar models indicate that two grains likely condensed in an intermediate-mass AGB star with initial metallicity of Z ≤ 0.014. Their 96 Zr/ 94 Zr ratios also match those predicted for born-again AGB stars undergoing a very late thermal pulse and rapidly accreting white dwarfs. After comparing the relative populations of the aforementioned dust-producing stars, we propose rapidly accreting white dwarfs as a new, and more likely, stellar source for one of the presolar grains. The remaining two grains could have originated in the supernova ejecta of massive stars, due to correlated excesses in the p-nuclides, 92, 94 Mo. Thus, grains with 96 Zr anomalies can have a variety of stellar origins, in agreement with theoretical studies. Conclusions. Our study highlights the importance of multi-element analysis in constraining the types of stars where presolar grains have condensed. These data will help improve our understanding of various nucleosynthesis processes in different stellar phases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Can Third-Body Stabilization of Bimolecular Collision Complexes in Cold Molecular Clouds Happen?

For more than half a century, networks of radiative association, dissociative recombination, and bimolecular reactions have been postulated to drive the low-temperature chemistry of cold molecular clouds. Third-body stabilizations of collision complexes have been assumed to be 'irrelevant' due to short lifetime of such complexes. Here, we conduct crossed molecular beam studies of ground state atomic silicon with diacetylene in combination with electronic structure calculations and microcanonical kinetics models operating under cold molecular cloud conditions. Our combined experimental, electronic structure, and microcanonical kinetics modelling investigations provide compelling evidence that three-body collisions of molecular hydrogen with long-lived reaction intermediates accessed through intersystem crossing are prevalent deep inside molecular clouds. This concept might be exportable to reactions involving polycyclic aromatic hydrocarbons thus affording a versatile machinery to complex organics via third-body stabilizations of bimolecular collision complexes deep inside cold molecular clouds.

astrochemistry↗
DOE OSTI2023

CO and [C ii] line emission of molecular clouds: the impact of stellar feedback and non-equilibrium chemistry

We analyse synthetic 12 CO, 13 CO, and [C ii] emission maps of molecular cloud (MC) simulations from the SILCC-Zoom project. We present radiation, magnetohydrodynamic zoom-in simulations of individual clouds, both with and without radiative stellar feedback, forming in a turbulent multiphase interstellar medium following on-the-fly the evolution of e.g. H 2 , CO, and C + . We introduce a novel post-processing routine based on cloudy which accounts for higher ionization states of carbon due to stellar radiation in H ii regions. Synthetic emission maps of [C ii] in and around feedback bubbles show that the bubbles are largely devoid of [C ii], as recently found in observations, which we attribute to the further ionization of C+ into C 2+ . For both 12 CO and 13 CO, the cloud-averaged luminosity ratio, $L_\rm {CO}/L_\rm {[C\, \small {II}]}$, can neither be used as a reliable measure of the H 2 mass fraction nor of the evolutionary stage of the clouds. We note a relation between the $I_\rm {CO}/I_\rm {[C\, \small {II}]}$ intensity ratio and the H 2 mass fraction for individual pixels of our synthetic maps. The scatter, however, is too large to reliably infer the H 2 mass fraction. Finally, the assumption of chemical equilibrium overestimates H2 and CO masses by up to 150 and 50 per cent, respectively, and $L_\rm {CO}$ by up to 60 per cent. The masses of H and C + would be underestimated by 65 and 30 per cent, respectively, and $L_\rm {[C\, \small {II}]}$ by up to 35 per cent. Hence, the assumption of chemical equilibrium in MC simulations introduces intrinsic errors of a factor of 2 in chemical abundances, luminosities, and luminosity ratios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A new way to unravel the $^{12}$C($\alpha $,$\gamma $)$^{16}$O cross section components using neural networks

The $^{12}$C($\alpha $,$\gamma $)$^{16}$O reaction rate is crucial in determining the carbon-to-oxygen abundance ratio in stellar nucleosynthesis. Measuring this reaction’s cross section at stellar energies is challenging due to its extremely small value, approximately 10 -17 barn at E c.m. = 300 keV. To address this, R-matrix calculations are employed to extrapolate data to lower energies, requiring a comprehensive understanding of each contribution to the cross section. The dominant contributions to the cross section at stellar energies arise from electric dipole (E1) and electric quadrupole (E2) transitions to the ground state of 16 O, along with a significant cascade contribution. Traditionally, these contributions have been separated using the γ-ray angular distribution. In this work, we propose a novel technique using the energy distribution of the 16 O recoils at the focal plane. This method involves a neural network trained on detailed Monte Carlo simulations of the energy distribution of recoils transported through the recoil mass separator ERNA. This approach enables the simultaneous determination of all three contributions with errors around 10% in the energy range E c.m. = 1.0–2.2 MeV. In conclusion, by employing this new technique, we aim to significantly improve the accuracy of determining the cross section of the $^{12}$C($\alpha $,$\gamma $)$^{16}$O reaction at astrophysical energies.

Astrochemistry↗
DOE OSTI2024

Ultraviolet H 2 luminescence in molecular clouds induced by cosmic rays

Context. Galactic cosmic rays (CRs) play a crucial role in ionisation, dissociation, and excitation processes within dense cloud regions where UV radiation is absorbed by dust grains and gas species. CRs regulate the abundance of ions and radicals, leading to the formation of more and more complex molecular species, and determine the charge distribution on dust grains. A quantitative analysis of these effects is essential for understanding the dynamical and chemical evolution of star-forming regions. Aims. The CR-induced photon flux has a significant impact on the evolution of the dense molecular medium in its gas and dust components. This study evaluates the flux of UV photons generated by CRs to calculate the photon-induced dissociation and ionisation rates of a vast number of atomic and molecular species, as well as the integrated UV photon flux. Methods. To achieve these goals, we took advantage of recent developments in the determination of the spectra of secondary electrons, in the calculation of state-resolved excitation cross sections of H 2 by electron impact, and of photodissociation and photoionisation cross sections. Results. We calculated the H 2 level population of each rovibrational level of the X, B, C, B', D, B", D', and a states. We then computed the UV photon spectrum of H 2 in its line and continuum components between 72 and 700 nm, with unprecedented accuracy, as a function of the CR spectrum incident on a molecular cloud, the H 2 column density, the isomeric H 2 composition, and the dust properties. The resulting photodissociation and photoionisation rates are, on average, lower than previous determinations by a factor of about 2, with deviations of up to a factor of 5 for the photodissociation of species such as AlH, C 2 H 2 , C 2 H 3 , C 3 H 3 , LiH, N 2 , NaCl, NaH, O$^+_2$, S 2 , SiH, l-C 4 , and l-C 5 H. A special focus is given to the photoionisation rates of H 2 , HF, and N 2 , as well as to the photodissociation of H 2 , which we find to be orders of magnitude higher than previous estimates. We give parameterisations for both the photorates and the integrated UV photon flux as a function of the CR ionisation rate, which implicitly depends on the H 2 column density, as well as the dust properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Implementation of Chemistry in the Athena++ Code

Chemistry plays a key role in many aspects of astrophysical fluids. Atoms and molecules are agents for heating and cooling, determine the ionization fraction, serve as observational tracers, and build the molecular foundation of life. We present the implementation of a chemistry module in the publicly available magnetohydrodynamic code Athena++ . We implement several chemical networks and heating and cooling processes suitable for simulating the interstellar medium (ISM). A general chemical network framework in the KIDA format is also included, allowing users to easily implement their own chemistry. Radiation transfer and cosmic-ray ionization are coupled with chemistry and solved with the simple six-ray approximation. The chemical and thermal processes are evolved as a system of coupled ordinary differential equations with an implicit solver from the CVODE library. We perform and present a series of tests to ensure the numerical accuracy and convergence of the code. Many tests combine chemistry with gas dynamics, including comparisons with analytic solutions, 1D problems of the photodissociation regions and shocks, and realistic 3D simulations of the turbulent ISM. We release the code with the new public version of Athena++ , aiming to provide a robust and flexible code for the astrochemical simulation community.

79 ASTRONOMY AND ASTROPHYSICS↗