-
Atacama Large Aperture Submillimeter Telescope (AtLAST) Science: Our Galaxy
Authors:
Pamela Klaassen,
Alessio Traficante,
Maria T. Beltrán,
Kate Pattle,
Mark Booth,
Joshua B. Lovell,
Jonathan P. Marshall,
Alvaro Hacar,
Brandt A. L. Gaches,
Caroline Bot,
Nicolas Peretto,
Thomas Stanke,
Doris Arzoumanian,
Ana Duarte Cabral,
Gaspard Duchêne,
David J. Eden,
Antonio Hales,
Jens Kauffmann,
Patricia Luppe,
Sebastian Marino,
Elena Redaelli,
Andrew J. Rigby,
Álvaro Sánchez-Monge,
Eugenio Schisano,
Dmitry A. Semenov
, et al. (16 additional authors not shown)
Abstract:
As we learn more about the multi-scale interstellar medium (ISM) of our Galaxy, we develop a greater understanding for the complex relationships between the large-scale diffuse gas and dust in Giant Molecular Clouds (GMCs), how it moves, how it is affected by the nearby massive stars, and which portions of those GMCs eventually collapse into star forming regions. The complex interactions of those…
▽ More
As we learn more about the multi-scale interstellar medium (ISM) of our Galaxy, we develop a greater understanding for the complex relationships between the large-scale diffuse gas and dust in Giant Molecular Clouds (GMCs), how it moves, how it is affected by the nearby massive stars, and which portions of those GMCs eventually collapse into star forming regions. The complex interactions of those gas, dust and stellar populations form what has come to be known as the ecology of our Galaxy. Because we are deeply embedded in the plane of our Galaxy, it takes up a significant fraction of the sky, with complex dust lanes scattered throughout the optically recognisable bands of the Milky Way. These bands become bright at (sub-)millimetre wavelengths, where we can study dust thermal emission and the chemical and kinematic signatures of the gas. To properly study such large-scale environments, requires deep, large area surveys that are not possible with current facilities. Moreover, where stars form, so too do planetary systems, growing from the dust and gas in circumstellar discs, to planets and planetesimal belts. Understanding the evolution of these belts requires deep imaging capable of studying belts around young stellar objects to Kuiper belt analogues around the nearest stars. Here we present a plan for observing the Galactic Plane and circumstellar environments to quantify the physical structure, the magnetic fields, the dynamics, chemistry, star formation, and planetary system evolution of the galaxy in which we live with AtLAST; a concept for a new, 50m single-dish sub-mm telescope with a large field of view which is the only type of facility that will allow us to observe our Galaxy deeply and widely enough to make a leap forward in our understanding of our local ecology.
△ Less
Submitted 1 March, 2024;
originally announced March 2024.
-
The Astrochemistry Low-energy Electron Cross-Section (ALeCS) database I. Semi-empirical electron-impact ionization cross-section calculations and ionization rates
Authors:
Brandt A. L. Gaches,
Tommaso Grassi,
Stefan Vogt-Geisse,
Giulia M. Bovolenta,
Claire Vallance,
David Heathcote,
Marco Padovani,
Stefano Bovino,
Prasanta Gorai
Abstract:
(Abridged) Electron-molecule interaction is a fundamental process in radiation-driven chemistry in space, from the interstellar medium to comets. Therefore, knowledge of interaction cross-sections is key. While there has been a plethora of studies of total ionization cross-sections, data is often spread over many sources, or not public or readily available. We introduce the Astrochemistry Low-ener…
▽ More
(Abridged) Electron-molecule interaction is a fundamental process in radiation-driven chemistry in space, from the interstellar medium to comets. Therefore, knowledge of interaction cross-sections is key. While there has been a plethora of studies of total ionization cross-sections, data is often spread over many sources, or not public or readily available. We introduce the Astrochemistry Low-energy Electron Cross-Section (ALeCS) database, a public database for electron interaction cross-sections and ionization rates for molecules of astrochemical interest. In this work, we present the first data release comprising total ionization cross-sections and ionization rates for over 200 neutral molecules. We include optimized geometries and molecular orbital energies at various levels of theory, and for a subset of the molecules, the ionization potentials. We compute total ionization cross-sections using the binary-encounter Bethe model and screening-corrected additivity rule, and ionization rates and reaction network coefficients for molecular cloud environments for $>$200 neutral molecules ranging from diatomics to complex organics. We demonstrate that our binary-encounter Bethe cross-sections agree well with experimental data. We show that the ionization rates scale roughly linearly with the number of constituent atoms in the molecule. We introduce and describe the public ALeCS database. For the initial release, we include total ionization cross-sections for $>$200 neutral molecules and several cations and anions calculated with different levels of quantum chemistry theory, the chemical reaction rates for the ionization, and network files in the formats of the two most popular astrochemical networks, the KIDA and UMIST. The database will be continuously updated for more molecules and interactions.
△ Less
Submitted 13 February, 2024; v1 submitted 16 October, 2023;
originally announced October 2023.
-
Tree-based solvers for adaptive mesh refinement code FLASH -- IV: An X-ray radiation scheme to couple discrete and diffuse X-ray emission sources to the thermochemistry of the interstellar medium
Authors:
Brandt A. L. Gaches,
Stefanie Walch,
Richard Wünsch,
Jonathan Mackey
Abstract:
X-ray radiation, in particular radiation between 0.1 keV 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…
▽ More
X-ray radiation, in particular radiation between 0.1 keV 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 {\sc TreeRay} reverse raytrace algorithm implemented in the {\sc Flash} magneto-hydrodynamic code which enables the inclusion of X-ray radiation from 0.1 keV $< E_γ <$ 100 keV, dubbed {\rm XrayTheSpot}. {\sc 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 multi-bin 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 disk and simulations of molecular clouds with active chemistry, radiation pressure, protostellar radiation feedback from infrared to X-ray radiation.
△ Less
Submitted 29 April, 2023; v1 submitted 30 January, 2023;
originally announced January 2023.
-
Cosmic-ray-induced H$_2$ line emission: Astrochemical modeling and implications for JWST observations
Authors:
Brandt A. L. Gaches,
Shmuel Bialy,
Thomas G. Bisbas,
Marco Padovani,
Daniel Seifried,
Stefanie Walch
Abstract:
Context: It has been proposed that H$_2$ near-infrared lines may be excited by cosmic rays and allow for a determination of the cosmic-ray ionization rate in dense gas. One-dimensional models show that measuring both the H$_2$ gas column density and H$_2$ line intensity enables a constraint on the cosmic-ray ionization rate as well as the spectral slope of low-energy cosmic-ray protons in the inte…
▽ More
Context: It has been proposed that H$_2$ near-infrared lines may be excited by cosmic rays and allow for a determination of the cosmic-ray ionization rate in dense gas. One-dimensional models show that measuring both the H$_2$ gas column density and H$_2$ line intensity enables a constraint on the cosmic-ray ionization rate as well as the spectral slope of low-energy cosmic-ray protons in the interstellar medium (ISM). Aims: We aim to investigate the impact of certain assumptions regarding the H$_2$ chemical models and ISM density distributions on the emission of cosmic-ray induced H$_2$ emission lines. This is of particular importance for utilizing observations of these lines with the James Webb Space Telescope to constrain the cosmic-ray ionization rate. Methods: We compare the predicted emission from cosmic-ray induced, ro-vibrationally excited H$_2$ emission lines for different one- and three-dimensional models with varying assumptions on the gas chemistry and density distribution. Results: We find that the model predictions of the H$_2$ line intensities for the (1-0)S(0), (1-0)Q(2), (1-0)O(2) and (1-0)O(4) transitions at 2.22, 2.41, 2.63 and 3.00 $μ$m, respectively, are relatively independent of the astro-chemical model and the gas density distribution when compared against the H$_2$ column density, making them robust tracer of the cosmic-ray ionization rate. Conclusions: We recommend the use of ro-vibrational H$_2$ line emission in combination with estimation of the cloud's H$_2$ column density, to constrain the ionization rate and the spectrum of low energy cosmic-rays.
△ Less
Submitted 6 July, 2022; v1 submitted 20 June, 2022;
originally announced June 2022.
-
The impact of cosmic-ray attenuation on the carbon cycle emission in molecular clouds
Authors:
Brandt A. L. Gaches,
Thomas G. Bisbas,
Shmuel Bialy
Abstract:
Observations of the emission of the carbon cycle species (C, C+ CO) are commonly used to diagnose gas properties in the interstellar medium but are significantly sensitive to the cosmic-ray ionization rate. The carbon-cycle chemistry is known to be quite sensitive to the cosmic-ray ionization rate, $ζ$, controlled by the flux of low-energy cosmic rays which get attenuated through molecular clouds.…
▽ More
Observations of the emission of the carbon cycle species (C, C+ CO) are commonly used to diagnose gas properties in the interstellar medium but are significantly sensitive to the cosmic-ray ionization rate. The carbon-cycle chemistry is known to be quite sensitive to the cosmic-ray ionization rate, $ζ$, controlled by the flux of low-energy cosmic rays which get attenuated through molecular clouds. However, astrochemical models commonly assume a constant cosmic-ray ionization rate in the clouds. We investigate the effect of cosmic-ray attenuation on the emission of carbon cycle species from molecular clouds, in particular the [CII] 158 $μ$m, [CI] 609 $μ$m and CO (J = 1 - 0) 115.27 GHz lines. We use a post-processed chemical model of diffuse and dense simulated molecular clouds and quantify the variation in both column densities and velocity integrated line emission of the carbon cycle with different cosmic-ray ionization rate models. We find that the abundances and column densities of carbon cycle species is significantly impacted by the chosen cosmic-ray ionization rate model: no single constant ionization rate can reproduce the abundances modelled with an attenuated cosmic-ray model. Further, we show that constant ionization rate models fail to simultaneously reproduce the integrated emission of the lines we consider, and their deviations from a physically derived cosmic-ray attenuation model is too complex to be simply corrected. We demonstrate that the two clouds we model exhibit a similar average $A_{\rm V, eff}$ -- $n_{\rm H}$ relationship, resulting in an average relation between the cosmic-ray ionization rate and density $ζ(n_{\rm H})$. We conclude by providing a number of implementation recommendations for CRs in astrochemical models, but emphasize the necessity for column-dependent cosmic-ray ionization rate prescriptions.
△ Less
Submitted 11 January, 2022; v1 submitted 13 October, 2021;
originally announced October 2021.
-
CRAFT (Cosmic Ray Acceleration From Turbulence) in Molecular Clouds
Authors:
Brandt A. L. Gaches,
Stefanie Walch,
Alex Lazarian
Abstract:
Low-energy cosmic rays, in particular protons with energies below 1 GeV, are significant drivers of the thermochemistry of molecular clouds. However, these cosmic rays are also greatly impacted by energy losses and magnetic field transport effects in molecular gas. Explaining cosmic ray ionization rates of $10^{-16}$ s$^{-1}$ or greater in dense gas requires either a high external cosmic ray flux,…
▽ More
Low-energy cosmic rays, in particular protons with energies below 1 GeV, are significant drivers of the thermochemistry of molecular clouds. However, these cosmic rays are also greatly impacted by energy losses and magnetic field transport effects in molecular gas. Explaining cosmic ray ionization rates of $10^{-16}$ s$^{-1}$ or greater in dense gas requires either a high external cosmic ray flux, or local sources of MeV-GeV cosmic ray protons. We present a new local source of low-energy cosmic rays in molecular clouds: first order Fermi-acceleration of protons in regions undergoing turbulent reconnection in molecular clouds. We show from energetic-based arguments there is sufficient energy within the magneto-hydrodynamic turbulent cascade to produce ionization rates compatible with inferred ionization rates in molecular clouds. As turbulent reconnection is a volume-filling process, the proposed mechanism can produce a near-homogeneous distribution of low-energy cosmic rays within molecular clouds.
△ Less
Submitted 6 August, 2021;
originally announced August 2021.
-
TIMES II: Investigating the Relation Between Turbulence and Star-forming Environments in Molecular Clouds
Authors:
Hyeong-Sik Yun,
Jeong-Eun Lee,
Neal J. Evans II,
Stella S. R. Offner,
Mark H. Heyer,
Jungyeon Cho,
Brandt A. L. Gaches,
Yao-Lun Yang,
How-Huan Chen,
Yunhee Choi,
Yong-Hee Lee,
Giseon Baek,
Minho Choi,
Jongsoo Kim,
Hyunwoo Kang,
Seokho Lee,
Ken'ichi Tatematsu
Abstract:
We investigate the effect of star formation on turbulence in the Orion A and Ophiuchus clouds using principal component analysis (PCA). We measure the properties of turbulence by applying PCA on the spectral maps in $^{13}$CO, C$^{18}$O, HCO$^+$ $J=$1$-$0, and CS $J=$2$-$1. First, the scaling relations derived from PCA of the $^{13}$CO maps show that the velocity difference ($δv$) for a given spat…
▽ More
We investigate the effect of star formation on turbulence in the Orion A and Ophiuchus clouds using principal component analysis (PCA). We measure the properties of turbulence by applying PCA on the spectral maps in $^{13}$CO, C$^{18}$O, HCO$^+$ $J=$1$-$0, and CS $J=$2$-$1. First, the scaling relations derived from PCA of the $^{13}$CO maps show that the velocity difference ($δv$) for a given spatial scale ($L$) is the highest in the integral shaped filament (ISF) and L1688, where the most active star formation occurs in the two clouds. The $δv$ increases with the number density and total bolometric luminosity of the protostars in the sub-regions. Second, in the ISF and L1688 regions, the $δv$ of C$^{18}$O, HCO$^+$, and CS are generally higher than that of $^{13}$CO, which implies that the dense gas is more turbulent than the diffuse gas in the star-forming regions; stars form in dense gas, and dynamical activities associated with star formation, such as jets and outflows, can provide energy into the surrounding gas to enhance turbulent motions.
△ Less
Submitted 28 July, 2021;
originally announced July 2021.
-
TIMES I: a Systematic Observation in Multiple Molecular Lines Toward the Orion A and Ophiuchus Clouds
Authors:
Hyeong-Sik Yun,
Jeong-Eun Lee,
Yunhee Choi,
Neal J. Evans II,
Stella S. R. Offner,
Mark H. Heyer,
Brandt A. L. Gaches,
Yong-Hee Lee,
Giseon Baek,
Minho Choi,
Hyunwoo Kang,
Seokho Lee,
Ken'ichi Tatematsu,
Yao-Lun Yang,
How-Huan Chen,
Youngung Lee,
Jae Hoon Jung,
Changhoon Lee,
Jungyeon Cho
Abstract:
We have used the Taeduk Radio Astronomy Observatory to observe the Orion A and Ophiuchus clouds in the $J=$1$-$0 lines of $^{13}$CO, C$^{18}$O, HCN, HCO$^+$, and N$_2$H$^+$ and the $J=$2$-$1 line of CS. The fully sampled maps with uniform noise levels are used to create moment maps. The variations of the line intensity and velocity dispersion with total column density, derived from dust emission m…
▽ More
We have used the Taeduk Radio Astronomy Observatory to observe the Orion A and Ophiuchus clouds in the $J=$1$-$0 lines of $^{13}$CO, C$^{18}$O, HCN, HCO$^+$, and N$_2$H$^+$ and the $J=$2$-$1 line of CS. The fully sampled maps with uniform noise levels are used to create moment maps. The variations of the line intensity and velocity dispersion with total column density, derived from dust emission maps, are presented and compared to previous work. The CS line traces dust column density over more than one order of magnitude, and the N$_2$H$^+$ line best traces the highest column density regime ($\log(N_\mathrm{H_2}$) $>$ 22.8). Line luminosities, integrated over the cloud, are compared to those seen in other galaxies. The HCO$^+$-to-HCN luminosity ratio in the Orion A cloud is similar to that of starburst galaxies, while that in the Ophiuchus cloud is in between those of active galactic nuclei and starburst galaxies.
△ Less
Submitted 8 June, 2021;
originally announced June 2021.
-
Transport of Protostellar Cosmic Rays in Turbulent Dense Cores
Authors:
Margot Fitz Axen,
Stella S. S. Offner,
Brandt A. L. Gaches,
Chris L. Fryer,
Aimee Hungerford,
Kedron Silsbee
Abstract:
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…
▽ More
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 non-uniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher then 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 $ζ\approx 10^{-16} \rm s^{-1}$; however, within $r \approx 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 sub-grid 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.
△ Less
Submitted 30 April, 2021;
originally announced May 2021.
-
Aluminum-26 Enrichment in the Surface of Protostellar Disks Due to Protostellar Cosmic Rays
Authors:
Brandt A. L. Gaches,
Stefanie Walch,
Stella S. R. Walch,
Carsten Münker
Abstract:
The radioactive decay of aluminum-26 ($^{26}$Al) is an important heating source in early planet formation. Since its discovery, there have been several mechanisms proposed to introduce $^{26}$Al into protoplanetary disks, primarily through contamination by external sources. We propose a local mechanism to enrich protostellar disks with $^{26}$Al through irradiation of the protostellar disk surface…
▽ More
The radioactive decay of aluminum-26 ($^{26}$Al) is an important heating source in early planet formation. Since its discovery, there have been several mechanisms proposed to introduce $^{26}$Al into protoplanetary disks, primarily through contamination by external sources. We propose a local mechanism to enrich protostellar disks with $^{26}$Al through irradiation of the protostellar disk surface by cosmic rays accelerated in the protostellar accretion shock. We calculate the $^{26}$Al enrichment, [$^{26}$Al/$^{27}$Al], at the surface of the protostellar disk in the inner AU throughout the evolution of low-mass stars, from M-dwarfs to proto-Suns. Assuming constant mass accretion rates, $\dot{m}$, we find that irradiation by MeV cosmic rays can provide significant enrichment on the disk surface if the cosmic rays are not completely coupled to the gas in the accretion flow. Importantly, we find that low accretion rates, $\dot{m} < 10^{-7}$ M$_{\odot}$ yr$^{-1}$, are able to produce canonical amounts of $^{26}$Al, $[^{26}{\rm Al}/^{27}{\rm Al}] \approx 5\times10^{-5}$. These accretion rates are experienced at the transition from Class I- to Class II-type protostars, when it is assumed that calcium-aluminum-rich inclusions condense in the inner disk. We conclude that irradiation of the inner disk surface by cosmic ray protons accelerated in accretion shocks at the protostellar surface may be an important mechanism to produce $^{26}$Al. Our models show protostellar cosmic rays may be a viable model to explain the enrichment of $^{26}$Al found in the Solar System.
△ Less
Submitted 24 July, 2020;
originally announced July 2020.
-
Impact of Cosmic-Ray Feedback on Accretion and Chemistry in Circumstellar Disks
Authors:
S. S. R. Offner,
B. A. L. Gaches,
J. Holdship
Abstract:
We use the gas-grain chemistry code UCLCHEM to explore the impact of cosmic-ray feedback on the chemistry of circumstellar disks. We model the attenuation and energy losses of the cosmic-rays as they propagate outwards from the star and also consider ionization due to stellar radiation and radionuclides. For accretion rates typical of young stars, $\dot M_* \sim 10^{-9}-10^{-6}$ M_\odot yr$^{-1}$,…
▽ More
We use the gas-grain chemistry code UCLCHEM to explore the impact of cosmic-ray feedback on the chemistry of circumstellar disks. We model the attenuation and energy losses of the cosmic-rays as they propagate outwards from the star and also consider ionization due to stellar radiation and radionuclides. For accretion rates typical of young stars, $\dot M_* \sim 10^{-9}-10^{-6}$ M_\odot yr$^{-1}$, we show that cosmic rays accelerated by the stellar accretion shock produce a cosmic-ray ionization rate at the disk surface $ζ\gtrsim 10^{-15}$ s$^{-1}$, at least an order of magnitude higher than the ionization rate associated with the Galactic cosmic-ray background. The incident cosmic-ray flux enhances the disk ionization at intermediate to high surface densities ($Σ> 10$ g cm$^{-2}$) particularly within 10 au of the star. We find the dominant ions are C$^+$, S$^+$ and Mg$^+$ in the disk surface layers, while the H$_3^+$ ion dominates at surface densities above 1.0 g cm$^{-2}$. We predict the radii and column densities at which the magneto-rotational instability (MRI) is active in T Tauri disks and show that ionization by cosmic-ray feedback extends the MRI-active region towards the disk mid-plane. However, the MRI is only active at the mid-plane of a minimum mass solar nebula disk if cosmic-rays propagate diffusively ($ζ\propto r^{-1}$) away from the star. The relationship between accretion, which accelerates cosmic rays, the dense accretion columns, which attenuate cosmic rays, and the MRI, which facilitates accretion, create a cosmic-ray feedback loop that mediates accretion and may produce luminosity variability.
△ Less
Submitted 10 September, 2019; v1 submitted 21 August, 2019;
originally announced August 2019.
-
The Astrochemical Impact of Cosmic Rays in Protoclusters II: CI-to-H$_2$ and CO-to-H$_2$ Conversion Factors
Authors:
Brandt A. L. Gaches,
Stella S. R. Offner,
Thomas G. Bisbas
Abstract:
We utilize a modified astrochemistry code which includes cosmic ray attenuation in-situ to quantify the impact of different cosmic-ray models on the CO-to-H$_2$ and CI-to-H$_2$ conversion factors, $X_{\rm CO}$ and $X_{\rm CI}$, respectively. We consider the impact of cosmic rays accelerated by accretion shocks, and show that clouds with star formation efficiencies greater than 2\% have…
▽ More
We utilize a modified astrochemistry code which includes cosmic ray attenuation in-situ to quantify the impact of different cosmic-ray models on the CO-to-H$_2$ and CI-to-H$_2$ conversion factors, $X_{\rm CO}$ and $X_{\rm CI}$, respectively. We consider the impact of cosmic rays accelerated by accretion shocks, and show that clouds with star formation efficiencies greater than 2\% have $X_{\rm CO} = (2.5 \pm 1)\times10^{20}$ cm$^{-2}$(K km s$^{-1}$)$^{-1}$, consistent with Milky Way observations. We find that changing the cosmic ray ionization rate from external sources from the canonical $ζ\approx 10^{-17}$ to $ζ\approx 10^{-16}$ s$^{-1}$, which better represents observations in diffuse gas, reduces $X_{\rm CO}$ by 0.2 dex for clusters with surface densities below 3 g cm$^{-2}$. We show that embedded sources regulate $X_{\rm CO}$ and decrease its variance across a wide range of surface densities and star formation efficiencies. Our models reproduce the trends of a decreased $X_{\rm CO}$ in extreme cosmic ray environments. $X_{\rm CI}$ has been proposed as an alternative to $X_{\rm CO}$ due to its brightness at high redshifts. The inclusion of internal cosmic ray sources leads to 1.2 dex dispersion in $X_{\rm CI}$ ranging from $2\times10^{20} < X_{\rm CI} < 4\times10^{21}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$. We show that $X_{\rm CI}$ is highly sensitive to the underlying cosmic ray model.
△ Less
Submitted 19 August, 2019;
originally announced August 2019.
-
The Astrochemical Impact of Cosmic Rays in Protoclusters I: Molecular Cloud Chemistry
Authors:
Brandt A. L. Gaches,
Stella S. R. Offner,
Thomas G. Bisbas
Abstract:
We present astrochemical photo-dissociation region models in which cosmic ray attenuation has been fully coupled to the chemical evolution of the gas. We model the astrochemical impact of cosmic rays, including those accelerated by protostellar accretion shocks, on molecular clouds hosting protoclusters. Our models with embedded protostars reproduce observed ionization rates. We study the imprint…
▽ More
We present astrochemical photo-dissociation region models in which cosmic ray attenuation has been fully coupled to the chemical evolution of the gas. We model the astrochemical impact of cosmic rays, including those accelerated by protostellar accretion shocks, on molecular clouds hosting protoclusters. Our models with embedded protostars reproduce observed ionization rates. We study the imprint of cosmic ray attenuation on ions for models with different surface cosmic ray spectra and different star formation efficiencies. We find that abundances, particularly ions, are sensitive to the treatment of cosmic rays. We show the column densities of ions are under predicted by the `classic' treatment of cosmic rays by an order of magnitude. We also test two common chemistry approximations used to infer ionization rates. We conclude that the approximation based on the H$_3^+$ abundance under predicts the ionization rate except in regions where the cosmic rays dominate the chemistry. Our models suggest the chemistry in dense gas will be significantly impacted by the increased ionization rates, leading to a reduction in molecules such as NH$_3$ and causing H$_2$-rich gas to become [C II] bright.
△ Less
Submitted 9 May, 2019; v1 submitted 6 May, 2019;
originally announced May 2019.
-
Exploration of Cosmic Ray Acceleration in Protostellar Accretion Shocks and A Model for Ionization Rates in Embedded Protoclusters
Authors:
Brandt A. L. Gaches,
Stella S. R. Offner
Abstract:
We construct a model for cosmic ray acceleration from protostellar accretion shocks and calculate the resulting cosmic ray ionization rate within star-forming molecular clouds. We couple a protostar cluster model with an analytic accretion shock model to calculate the cosmic ray acceleration from protostellar surfaces. We present the cosmic ray flux spectrum from keV to GeV energies for a typical…
▽ More
We construct a model for cosmic ray acceleration from protostellar accretion shocks and calculate the resulting cosmic ray ionization rate within star-forming molecular clouds. We couple a protostar cluster model with an analytic accretion shock model to calculate the cosmic ray acceleration from protostellar surfaces. We present the cosmic ray flux spectrum from keV to GeV energies for a typical low-mass protostar. We find that at the shock surface the spectrum follows a power-law trend across 6 orders of magnitude in energy. After attenuation, the spectrum at high energies steepens, while at low energies it is relatively flat. We calculate the cosmic ray pressure and cosmic ray ionization rate from relativistic protons at the protostellar surface and at the edge of the core. We present the cosmic ray ionization rate for individual protostars as a function of their instantaneous mass and final mass. The protostellar cosmic ray ionization rate is $ζ\approx 0.01 - 1$ s$^{-1}$ at the accretion shock surface. However, at the edge of the core, the cosmic ray ionization rate drops substantially to between $ζ\approx 10^{-20}$ to $10^{-17}$ s$^{-1}$. There is a large spatial gradient in the cosmic ray ionization rate, such that inner regions may experience cosmic ray ionization rates larger than the often assumed fiducial rate, $ζ= 3\times10^{-17}$ s$^{-1}$. Finally, we calculate the cosmic ray ionization rate for protostellar clusters over 5 orders of magnitude of cluster size. We find that clusters with more than approximately 200 protostars produce a higher cosmic ray ionization rate within their natal cloud than the fiducial galactic value.
△ Less
Submitted 25 May, 2018; v1 submitted 8 May, 2018;
originally announced May 2018.
-
A Model for Protostellar Cluster Luminosities and the Impact on the CO-H$_2$ Conversion Factor
Authors:
Brandt A. L. Gaches,
Stella S. R. Offner
Abstract:
We construct a semi-analytic model to study the effect of far-ultraviolet (FUV) radiation on gas chemistry from embedded protostars. We use the Protostellar Luminosity Function (PLF) formalism of Offner & McKee (2011) to calculate the total, FUV, and ionizing cluster luminosity for various protostellar accretion histories and cluster sizes. We compare the model predictions with surveys of Gould Be…
▽ More
We construct a semi-analytic model to study the effect of far-ultraviolet (FUV) radiation on gas chemistry from embedded protostars. We use the Protostellar Luminosity Function (PLF) formalism of Offner & McKee (2011) to calculate the total, FUV, and ionizing cluster luminosity for various protostellar accretion histories and cluster sizes. We compare the model predictions with surveys of Gould Belt star-forming regions and find the Tapered Turbulent Core model matches best the mean luminosities and the spread in the data. We combine the cluster model with the photo-dissociation region astrochemistry code, {\sc 3d-pdr}, to compute the impact of the FUV luminosity from embedded protostars on the CO to H$_2$ conversion factor, $X_{\rm CO}$, as a function of cluster size, gas mass and star formation efficiency. We find that $X_{\rm CO}$ has a weak dependence on the FUV radiation from embedded sources for large clusters due to high cloud optical depths. In smaller and more efficient clusters the embedded FUV increases X$_{\rm CO}$ to levels consistent with the average Milky Way values. The internal physical and chemical structure of the cloud are significantly altered, and $X_{\rm CO}$ depends strongly on the protostellar cluster mass for small efficient clouds.
△ Less
Submitted 25 January, 2018;
originally announced January 2018.
-
Astrochemical Correlations in Molecular Clouds
Authors:
Brandt A. L. Gaches,
Stella S. R. Offner,
Erik W. Rosolowsky,
Thomas G. Bisbas
Abstract:
We investigate the spectral correlations between different species used to observe molecular clouds. We use hydrodynamic simulations and a full chemical network to study the abundances of over 150 species in typical Milky Way molecular clouds. We perform synthetic observations in order to produce emission maps of a subset of these tracers. We study the effects of different lines of sight and spati…
▽ More
We investigate the spectral correlations between different species used to observe molecular clouds. We use hydrodynamic simulations and a full chemical network to study the abundances of over 150 species in typical Milky Way molecular clouds. We perform synthetic observations in order to produce emission maps of a subset of these tracers. We study the effects of different lines of sight and spatial resolution on the emission distribution and perform a robust quantitative comparison of the species to each other. We use the Spectral Correlation Function (SCF), which quantifies the root mean squared difference between spectra separated by some length scale, to characterize the structure of the simulated cloud in position-position-velocity (PPV) space. We predict the observed SCF for a broad range of observational tracers, and thus, identify homologous species. In particular, we show that the pairs C and CO, C$^{+}$ and CN, NH$_3$ and H$_2$CS have very similar SCFs. We measure the SCF slope variation as a function of beam size for all species and demonstrate that the beam size has a distinct effect on different species emission. However, for beams of up to 10", placing the cloud at 1 kpc, the change is not large enough to move the SCF slopes into different regions of parameter space. The results from this study provide observational guidance for choosing the best tracer to probe various cloud length scales.
△ Less
Submitted 8 December, 2014;
originally announced December 2014.
