Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

34 records · Page 2
All sourcesNASA NTRSDOE OSTI
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

A combined experimental and computational study on the reaction dynamics of the 1-propynyl (CH 3 CC, X 2 A 1 ) – propylene (CH 3 CHCH 2 , X 1 A') system: formation of 1,3-dimethylvinylacetylene (CH 3 CCCHCHCH 3 , X 1 A') under single collision conditions

The reaction of the 1-propynyl radical (CH 3 CC; X 2 A 1 ) with propylene (CH 3 CHCH 2 ; X 1 A') was studied in a crossed molecular beam machine at a collision energy of 37 ± 1 kJ mol –1 . Experimental data combined with high-level electronic structure (CCSD(T)-F12/cc-pVTZ-F12//ωB97X-D/6-311G(d,p)) and RRKM calculations reveal the reaction mechanism. Here, the overall barrierless and exoergic reaction involves indirect reaction dynamics and commences preferentially with addition of 1-propynyl with its radical centre to the carbon–carbon double bond at the terminal carbon atom of propylene. This work focuses on molecular mass growth process (hydrogen loss channels) although theory suggests methyl loss as a prevalent channel. In these processes, the C6H9 collision complexes either emit atomic hydrogen or undergo isomerisation followed by atomic hydrogen loss to preferentially yield the cis/trans isomers of 1,3-dimethylvinylacetylene (2-hexen-4-yne) as the primary product. Analysis of reaction dynamics of 1-propynyl and ethynyl radicals with propylene along with their fractional abundance in deep space suggests formation of methyl- and dimethyl derivatives of vinylacetylene in cold molecular clouds. Once formed they may engage in fundamental molecular mass growth processes via the barrierless Hydrogen Abstraction Vinylacetylene Addition mechanism that leads to the formation of methyl- and dimethylnaphthalenes thus providing a versatile route to methyl-substituted PAHs in interstellar medium.

1-propynyl↗
DOE OSTI2020

Revisiting the reactivity between HCO and CH 3 on interstellar grain surfaces

The formation of interstellar complex organic molecules is currently thought to be dominated by the barrierless coupling between radicals on the interstellar icy grain surfaces. Previous standard density functional theory (DFT) results on the reactivity between CH 3 and HCO on amorphous water surfaces showed that the formation of CH 4 + CO by H transfer from HCO to CH 3 assisted by water molecules of the ice was the dominant channel. In contrast, the adopted description of the electronic structure of the biradical (i.e. CH 3 /HCO) system was inadequate [without the broken-symmetry (BS) approach]. In this work, we revisit the original results by means of BS-DFT both in gas phase and with one water molecule simulating the role of the ice. Results indicate that the adoption of BS-DFT is mandatory to describe properly biradical systems. In the presence of the single water molecule, the water-assisted H transfer exhibits a high energy barrier. In contrast, CH 3 CHO formation is found to be barrierless. Yet, direct H transfer from HCO to CH 3 to give CO and CH 4 presents a very low energy barrier, hence being a potential competitive channel to the radical coupling and indicating, moreover, that the physical insights of the original work remain valid.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Inelastic, exchange, and reactive processes in rovibrationally excited collisions of HD with H

ABSTRACT The HD molecule is an important coolant in early universe chemistry models and a tracer of H2 in star-forming regions. Rate coefficients for collisional excitation and de-excitation of HD rotational and vibrational levels form important ingredients in astrophysical models. While collisions with He, H2, and H are the most important, available data for H + HD collisions are largely limited to temperatures less than 1000 K for the vibrational ground state, low-lying rotational levels of the v = 1 HD vibrational level, or computed without reactive contributions. Here, through explicit quantum scattering calculations, we report extensive data for rovibrational transitions in HD induced by H atoms for a range of rotational levels in v = 1 and some v = 0 levels for temperatures up to 1000 K. The significance of the computed results for astrophysical modeling is discussed.

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

Differences in chemical evolution between isolated and embedded prestellar cores

Models of prestellar cores often assume that the cores are isolated from their environment – material outside the core boundary plays no role in the subsequent evolution. This is unlikely to be the case in reality, where cores are located within hierarchically substructured molecular clouds. Here, we investigate the dynamical and chemical evolution of prestellar cores, modelled as Bonnor–Ebert spheres, and show that the density of the ambient medium has a large impact on the resulting chemical properties of the cores. Models embedded in high-density, low-temperature surroundings have greatly enhanced abundances of several molecules, such as CO and CS, compared to models with more diffuse surroundings, corresponding to relatively isolated cores. The predicted intensities and profile shapes of molecular lines are also affected. The density of the ambient medium has a stronger effect on the chemical evolution than whether the cores are initially in or out of equilibrium. This suggests that the impact of environment cannot be neglected when modelling chemistry in prestellar cores; the results of these models are highly sensitive to the assumptions made about the core surroundings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

State-to-state rate coefficients for HCS+ in rotationally inelastic collisions with H2 at low temperatures

ABSTRACT HCS+ ions have been detected in several regions of the interstellar medium (ISM), but an accurate determination of the chemical-physical conditions in the molecular clouds where this molecule is observed requires detailed knowledge of the collisional rate coefficients with the most common colliders in those environments. In this work, we study the dynamics of rotationally inelastic collisions of HCS+ + H2 at low temperature, and report, for the first time, a set of rate coefficients for this system. We used a recently developed potential energy surface for the HCS+–H2 van der Waals complex and computed state-to-state rotational rate coefficients for the lower rotational states of HCS+ in collision with both para- and ortho-H2, analysing the influence of the computed rate coefficients on the determination of critical densities. Additionally, the computed rate coefficients are compared with those obtained by scaling the ones from HCS+ in collision with He (an approximation that is sometimes used when data is lacking), and large differences are found. Furthermore, the approximation of using the rates for the HCO+ + H2 collision as a rough approximation for those of the HCS+ + H2 system is also evaluated. Finally, the complete set of de-excitation rate coefficients for the lowest 30 rotational states of HCS+ by collision with H2 is reported from 5 to 100 K.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Tree-based solvers for adaptive mesh refinement code $\scriptsize{FLASH}$ – IV. An X-ray radiation scheme to couple discrete and diffuse X-ray emission sources to the thermochemistry of the interstellar medium

X-ray radiation, in particular radiation between 0.1 and 10 keV, is evident from both point-like sources, such as compact objects and T-Tauri young stellar objects, and extended emission from hot, cooling gas, such as in supernova remnants. The X-ray radiation is absorbed by nearby gas, providing a source of both heating and ionization. While protoplanetary chemistry models now often include X-ray emission from the central young stellar object, simulations of star-forming regions have yet to include X-ray emission coupled to the chemo-dynamical evolution of the gas. We present an extension of the $\scriptsize{TREERAY}$ reverse ray trace algorithm implemented in the flash magnetohydrodynamic code which enables the inclusion of X-ray radiation from 0.1 keV < E γ < 100 keV, dubbed $\scriptsize{XRAYTHESPOT}$. $\scriptsize{XRAYTHESPOT}$ allows for the use of an arbitrary number of bins, minimum and maximum energies, and both temperature-independent and temperature-dependent user-defined cross-sections, along with the ability to include both point and extended diffuse emission and is coupled to the thermochemical evolution. We demonstrate the method with several multibin benchmarks testing the radiation transfer solution and coupling to the thermochemistry. Finally, we show two example star formation science cases for this module: X-ray emission from protostellar accretion irradiating an accretion disc and simulations of molecular clouds with active chemistry, radiation pressure, and protostellar radiation feedback from infrared to X-ray radiation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The evolution of HCO + in molecular clouds using a novel chemical post-processing algorithm

Modelling the chemistry of molecular clouds is critical to accurately simulating their evolution. To reduce computational cost, 3D simulations generally restrict their chemistry to species with strong heating and cooling effects. Time-dependent information about the evolution of other species is therefore often neglected. We address this gap by post-processing tracer particles in the SILCC-Zoom molecular cloud simulations. Using a chemical network of 39 species and 301 reactions (including freeze-out of CO and H 2 O) and a novel algorithm to reconstruct a density grid from sparse tracer particle data, we produce time-dependent density distributions for various species. We focus upon the evolution of HCO + , which is a critical formation reactant of CO but is not typically modelled on the fly. We find that ∼ 90 per cent of the HCO + content of the cold molecular gas forms in situ around n HCO + ∼ 10 3 –10 4 cm −3 , over a time-scale of approximately 1 Myr. The remaining ∼ 10 per cent forms at high extinction sites, with minimal turbulent mixing out into the less dense gas. We further show that the dominant HCO + formation pathway is dependent on the visual extinction, with the reaction H 3 + + CO contributing 90 per cent of the total HCO + production above A V, 3D = 3. We produce the very first maps of the HCO + column density, N(HCO + ), and show that it reaches values as high as 10 15 cm −2 . We find that 50 per cent of the HCO + mass is located within AV ∼ 10–30 in a density range of 10 3.5 –10 4.5 cm −3 . Our maps of N(HCO + ) are shown to be in good agreement with recent observations of the W49A star-forming region.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2021

Transport of Protostellar Cosmic Rays in Turbulent Dense Cores

We report recent studies have suggested that low-energy cosmic rays (CRs) may be accelerated inside molecular clouds by the shocks associated with star formation. We use a Monte Carlo transport code to model the propagation of CRs accelerated by protostellar accretion shocks through protostellar cores. We calculate the CR attenuation and energy losses and compute the resulting flux and ionization rate as a function of both radial distance from the protostar and angular position. We show that protostellar cores have nonuniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher than the radial average at a given distance. In particular, the CR flux is focused in the direction of the outflow cavity, creating a "flashlight" effect and allowing CRs to leak out of the core. The radially averaged ionization rates are less than the measured value for the Milky Way of ζ ≈ 10 -16 s -1 ; however, within r ≈ 0.03 pc from the protostar, the maximum ionization rates exceed this value. We show that variation in the protostellar parameters, particularly in the accretion rate, may produce ionization rates that are a couple of orders of magnitude higher or lower than our fiducial values. Finally, we use a statistical method to model unresolved subgrid magnetic turbulence in the core. We show that turbulence modifies the CR spectrum and increases the uniformity of the CR distribution but does not significantly affect the resulting ionization rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Laboratory-driven Multiscale Investigation of X-Ray Induced Mass Loss and Photochemical Evolution in Cosmic Carbon and Silicate Dust

We present the results of an integrated laboratory and modeling investigation into the impact of stellar X-rays on cosmic dust. Carbonaceous grains were prepared in a cooled (< 200 K) supersonic expansion from aromatic molecular precursors, and were later irradiated with 970 eV X-rays. Silicate (enstatite) grains were prepared via laser ablation, thermally annealed, and later irradiated with 500 eV X-rays. Infrared spectra of the 3.4 μm band of the carbon sample prepared with benzene revealed 84% ± 5% band area loss for an X-ray dose of 5.2 x 10 23 eV.cm -2 . Infrared spectra of the 8 -12 μm Si-O band of the silicate sample revealed band area loss up to 63% ± 5% for doses of 2.3 x 10 23 eV.cm -2 . A hybrid Monte Carlo particle trajectory approach was used to model the impact of X-rays and ensuing photoelectrons, Auger and collisionally ionized electrons through the bulk. As a result of x-ray ionization and ensuing Coulomb explosions on surface molecules, the calculated mass loss is 60% for the carbonaceous sample and 46% for the silicate sample, within a factor of 2 of the IR band loss, supporting an X-ray induced mass-loss mechanism. We apply the laboratory X-ray destruction rates to estimate the lifetimes of dust grains in protoplanetary disks surrounding 1 M ⊙ and 0.1 M ⊙ G and M stars. In both cases, X-ray destruction timescales are short (a few million years) at the disk surface, but are found to be much longer than typical disk lifetimes (≳10 Myr) over the disk bulk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Collisional Excitation and Non-LTE Modeling of Interstellar Chiral Propylene Oxide

Abstract The first set of theoretical rotational cross sections for propylene oxide (CH 3 CHCH 2 O) colliding with cold He atoms has been obtained at the full quantum level using a high-accuracy potential energy surface. By scaling the collision reduced mass, rotational rate coefficients for collisions with para-H 2 are deduced in the temperature range 5–30 K. These collisional coefficients are combined with radiative data in a non-LTE radiative transfer model in order to reproduce observations of propylene oxide made toward the Sagittarius B2(N) molecular cloud with the Green Bank and Parkes radio telescopes. The three detected absorption lines are found to probe the cold (∼10 K) and translucent ( n H ∼ 2000 cm −3 ) gas in the outer edges of the extended Sgr B2(N) envelope. The derived column density for propylene oxide is N tot ∼ 3 × 10 12 cm −2 , corresponding to a fractional abundance relative to total hydrogen of ∼2.5 × 10 −11 . The present results are expected to help our understanding of the chemistry of propylene oxide, including a potential enantiomeric excess, in the cold interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Rest-frame Submillimeter Spectrum of High-redshift, Dusty, Star-forming Galaxies from the SPT-SZ Survey

We present the average rest-frame spectrum of the final catalog of dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope's SPT-SZ survey and measured with Band 3 of the Atacama Large Millimeter/submillimeter Array. This work builds on the previous average rest-frame spectrum, given in Spilker et al. (2014) for the first 22 sources, and is comprised of a total of 78 sources, normalized by their respective apparent dust masses. The spectrum spans 1.9 < z < 6.9 and covers rest-frame frequencies of 240–800 GHz. Combining this data with low-J CO observations from the Australia Telescope Compact Array, we detect multiple bright line features from 12 CO, [Ci], and H 2 O, as well as fainter molecular transitions from 13 CO, HCN, HCO + , HNC, CN, H 2 O + , and CH. We use these detections, along with limits from other molecules, to characterize the typical properties of the interstellar medium (ISM) for these high-redshift DSFGs. We are able to divide the large sample into subsets in order to explore how the average spectrum changes with various galaxy properties, such as effective dust temperature. We find that systems with hotter dust temperatures exhibit differences in the bright 12 CO emission lines, and contain either warmer and more excited dense gas tracers or larger dense gas reservoirs. These observations will serve as a reference point to studies of the ISM in distant luminous DSFGs (L IR > 10 12 L ⊙ ), and will inform studies of chemical evolution before the peak epoch of star formation at z = 2–3.

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↗