BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Eric Herbst

Eric Herbst appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.SRastro-ph.EPastro-ph.IMastro-ph.COastro-ph.HEastro-phphysics.chem-ph

Catalog footprint

What is connected

45works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

45 published item(s)

preprint2022arXiv

FAUST III. Misaligned rotations of the envelope, outflow, and disks in the multiple protostellar system of VLA 1623$-$2417

We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Ice Age : Chemo-dynamical modeling of Cha-MMS1 to predict new solid-phase species for detection with JWST

Chemical models and experiments indicate that interstellar dust grains and their ice mantles play an important role in the production of complex organic molecules (COMs). To date, the most complex solid-phase molecule detected with certainty in the ISM is methanol, but the James Webb Space Telescope (JWST) may be able to identify still larger organic species. In this study, we use a coupled chemo-dynamical model to predict new candidate species for JWST detection toward the young star-forming core Cha-MMS1, combining the gas-grain chemical kinetic code MAGICKAL with a 1-D radiative hydrodynamics simulation using Athena++. With this model, the relative abundances of the main ice constituents with respect to water toward the core center match well with typical observational values, providing a firm basis to explore the ice chemistry. Six oxygen-bearing COMs (ethanol, dimethyl ether, acetaldehyde, methyl formate, methoxy methanol, and acetic acid), as well as formic acid, show abundances as high as, or exceeding, 0.01% with respect to water ice. Based on the modeled ice composition, the infrared spectrum is synthesized to diagnose the detectability of the new ice species. The contribution of COMs to IR absorption bands is minor compared to the main ice constituents, and the identification of COM ice toward the core center of Cha-MMS1 with the JWST NIRCAM/Wide Field Slitless Spectroscopy (2.4-5.0 micron) may be unlikely. However, MIRI observations (5-28 micron) toward COM-rich environments where solid-phase COM abundances exceed 1% with respect to the water ice column density might reveal the distinctive ice features of COMs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

PDRs4All: A JWST Early Release Science Program on radiative feedback from massive stars

Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the Universe, from the era of vigorous star formation at redshifts of 1-3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation Regions (PDRs) where the far-ultraviolet photons of massive stars create warm regions of gas and dust in the neutral atomic and molecular gas. PDR emission provides a unique tool to study in detail the physical and chemical processes that are relevant for most of the mass in inter- and circumstellar media including diffuse clouds, proto-planetary disks and molecular cloud surfaces, globules, planetary nebulae, and star-forming regions. PDR emission dominates the infrared (IR) spectra of star-forming galaxies. Most of the Galactic and extragalactic observations obtained with the James Webb Space Telescope (JWST) will therefore arise in PDR emission. In this paper we present an Early Release Science program using the MIRI, NIRSpec, and NIRCam instruments dedicated to the observations of an emblematic and nearby PDR: the Orion Bar. These early JWST observations will provide template datasets designed to identify key PDR characteristics in JWST observations. These data will serve to benchmark PDR models and extend them into the JWST era. We also present the Science-Enabling products that we will provide to the community. These template datasets and Science-Enabling products will guide the preparation of future proposals on star-forming regions in our Galaxy and beyond and will facilitate data analysis and interpretation of forthcoming JWST observations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Kiloparsec-scale Neutral Atomic Carbon Outflow in the Nearby Type-2 Seyfert Galaxy NGC 1068: Evidence for Negative AGN Feedback

Active galactic nucleus (AGN) feedback is postulated as a key mechanism for regulating star formation within galaxies. Studying the physical properties of the outflowing gas from AGN is thus crucial for understanding the co-evolution of galaxies and supermassive black holes. Here we report 55 pc resolution ALMA neutral atomic carbon [CI] $^3P_1\text{-}^3P_0$ observations toward the central 1 kpc of the nearby type-2 Seyfert galaxy NGC 1068, supplemented by 55 pc resolution CO($J=1\text{-}0$) observations. We find that [CI] emission within the central kpc is strongly enhanced by a factor of $>$5 compared to the typical [CI]/CO intensity ratio of $\sim$0.2 for nearby starburst galaxies (in units of brightness temperature). The most [CI]-enhanced gas (ratio $>$ 1) exhibits a kpc-scale elongated structure centered at the AGN that matches the known biconical ionized gas outflow entraining molecular gas in the disk. A truncated, decelerating bicone model explains well the kinematics of the elongated structure, indicating that the [CI] enhancement is predominantly driven by the interaction between the ISM in the disk and the highly inclined ionized gas outflow (which is likely driven by the radio jet). Our results strongly favor the "CO dissociation scenario" rather than the "in-situ C formation" one which prefers a perfect bicone geometry. We suggest that the high [CI]/CO intensity ratio gas in NGC 1068 directly traces ISM in the disk that is currently dissociated and entrained by the jet and the outflow, i.e., the "negative" effect of the AGN feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Vibrationally-excited Lines of HC$_{3}$N Associated with the Molecular Disk around the G24.78+0.08 A1 Hyper-compact H$_{\rm {II}}$ Region

We have analyzed Atacama Large Millimeter/submillimeter Array Band 6 data of the hyper-compact H$_{\rm {II}}$ region G24.78+0.08 A1 (G24 HC H$_{\rm {II}}$) and report the detection of vibrationally-excited lines of HC$_{3}$N ($v_{7}=2$, $J=24-23$). The spatial distribution and kinematics of a vibrationally-excited line of HC$_{3}$N ($v_{7}=2$, $J=24-23$, $l=2e$) are found to be similar to the CH$_{3}$CN vibrationally-excited line ($v_{8}=1$), which indicates that the HC$_{3}$N emission is tracing the disk around the G24 HC H$_{\rm {II}}$ region previously identified by the CH$_{3}$CN lines. We derive the $^{13}$CH$_{3}$CN/HC$^{13}$CCN abundance ratios around G24 and compare them to the CH$_{3}$CN/HC$_{3}$N abundance ratios in disks around Herbig Ae and T Tauri stars. The $^{13}$CH$_{3}$CN/HC$^{13}$CCN ratios around G24 ($\sim 3.0-3.5$) are higher than the CH$_{3}$CN/HC$_{3}$N ratios in the other disks ($\sim 0.03-0.11$) by more than one order of magnitude. The higher CH$_{3}$CN/HC$_{3}$N ratios around G24 suggest that the thermal desorption of CH$_{3}$CN in the hot dense gas and efficient destruction of HC$_{3}$N in the region irradiated by the strong UV radiation are occurring. Our results indicate that the vibrationally-excited HC$_{3}$N lines can be used as a disk tracer of massive protostars at the HC H$_{\rm {II}}$ region stage, and the combination of these nitrile species will provide information of not only chemistry but also physical conditions of the disk structures.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Extraterrestrial Photochemistry: Principles and Applications

Energetic processing of interstellar ice mantles and planetary atmospheres via photochemistry is a critical mechanism in the extraterrestrial synthesis of prebiotic molecules. Photochemistry is defined as chemical processes initiated by photon-induced electronic excitation, not involving ionization. In contrast, photons with energies above the ionization threshold initiate radiation chemistry (radiolysis). Vacuum-ultraviolet (6.2-12.4 eV) light may initiate photochemistry and radiation chemistry because the threshold for producing secondary electrons is lower in the condensed phase than in the gas phase. Approximately half of cosmic-ray induced photons incident on interstellar ices in star-forming regions initiate photochemistry while the rest initiate radiation chemistry. While experimental techniques such as velocity map imaging may be used to extract exquisite details about gas-phase photochemistry, such detailed information cannot be obtained for condensed-phase photochemistry, which involves greater complexity, including the production of excitons, excimers, and exciplexes. Because a primary objective of chemistry is to provide molecular-level mechanistic explanations for macroscopic phenomena, our ultimate goal in this book chapter is to critically evaluate our current understanding of the photochemistry that likely leads to the synthesis of extraterrestrial prebiotic molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Formation of complex organic molecules in hot molecular cores through nondiffusive grain-surface and ice-mantle chemistry

A new, more comprehensive model of gas-grain chemistry in hot molecular cores is presented, in which nondiffusive reaction processes on dust-grain surfaces and in ice mantles are implemented alongside traditional diffusive surface/bulk-ice chemistry. We build on our nondiffusive treatments used for chemistry in cold sources, adopting a standard collapse/warm-up physical model for hot cores. A number of other new chemical model inputs and treatments are also explored in depth, culminating in a final model that demonstrates excellent agreement with gas-phase observational abundances for many molecules, including some (e.g. methoxymethanol) that could not be reproduced by conventional diffusive mechanisms. Observed ratios of structural isomers methyl formate, glycolaldehyde and acetic acid are well reproduced by the models. The main temperature regimes are identified in which various complex organic molecules (COMs) are formed. Nondiffusive chemistry advances the production of many COMs to much earlier times and lower temperatures than in previous model implementations. Those species may form either as by-products of simple-ice production, or via early photochemistry within the ices while external UV photons can still penetrate. Cosmic ray-induced photochemistry is less important than in past models, although it affects some species strongly over long timescales. Another production regime occurs during the high-temperature desorption of solid water, whereby radicals trapped in the ice are released onto the grain/ice surface, where they rapidly react. Several recently-proposed gas-phase COM-production mechanisms are also introduced, but they rarely dominate. New surface/ice reactions involving CH and CH$_2$ are found to contribute substantially to the formation of certain COMs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Carbon Chain Chemistry in Hot-Core Regions around Three Massive Young Stellar Objects Associated with 6.7 GHz Methanol Masers

We have carried out observations of CCH ($N=1-0$), CH$_{3}$CN ($J=5-4$), and three $^{13}$C isotopologues of HC$_{3}$N ($J=10-9$) toward three massive young stellar objects (MYSOs), G12.89+0.49, G16.86--2.16, and G28.28--0.36, with the Nobeyama 45-m radio telescope. Combined with previous results on HC$_{5}$N, the column density ratios of $N$(CCH)/$N$(HC$_{5}$N), hereafter the CCH/HC$_{5}$N ratios, in the MYSOs are derived to be $\sim 15$. This value is lower than that in a low-mass warm carbon chain chemistry (WCCC) source by more than one order of magnitude. We compare the observed CCH/HC$_{5}$N ratios with hot-core model calculations (Taniguchi et al. 2019). The observed ratios in the MYSOs can be best reproduced by models when the gas temperature is $\sim 85$ K, which is higher than in L1527, a low-mass WCCC source ($\sim 35$ K). These results suggest that carbon-chain molecules detected around the MYSOs exist at least partially in higher temperature regions than those in low-mass WCCC sources. There is no significant difference in column density among the three $^{13}$C isotopologues of HC$_{3}$N in G12.89+0.49 and G16.86-2.16, while HCC$^{13}$CN is more abundant than the others in G28.28--0.36. We discuss carbon-chain chemistry around the three MYSOs based on the CCH/HC$_{5}$N ratio and the $^{13}$C isotopic fractionation of HC$_{3}$N.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Detection of Interstellar HC$_4$NC and an Investigation of Isocyanopolyyne Chemistry under TMC-1 Conditions

We report an astronomical detection of HC$_4$NC for the first time in the interstellar medium with the Green Bank Telescope toward the TMC-1 molecular cloud with a minimum significance of $10.5 σ$. The total column density and excitation temperature of HC$_4$NC are determined to be $3.29^{+8.60}_{-1.20}\times 10^{11}$ cm$^{-2}$ and $6.7^{+0.3}_{-0.3}$ K, respectively, using the MCMC analysis. In addition to HC$_4$NC, HCCNC is distinctly detected whereas no clear detection of HC$_6$NC is made. We propose that the dissociative recombination of the protonated cyanopolyyne, HC$_5$NH$^+$, and the protonated isocyanopolyyne, HC$_4$NCH$^+$, are the main formation mechanisms for HC$_4$NC while its destruction is dominated by reactions with simple ions and atomic carbon. With the proposed chemical networks, the observed abundances of HC$_4$NC and HCCNC are reproduced satisfactorily.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Early Science from GOTHAM: Project Overview, Methods, and the Detection of Interstellar Propargyl Cyanide (HCCCH$_2$CN) in TMC-1

We present an overview of the GOTHAM (GBT Observations of TMC-1: Hunting Aromatic Molecules) Large Program on the Green Bank Telescope. This and a related program were launched to explore the depth and breadth of aromatic chemistry in the interstellar medium at the earliest stages of star formation, following our earlier detection of benzonitrile ($c$-C$_6$H$_5$CN) in TMC-1. In this work, details of the observations, use of archival data, and data reduction strategies are provided. Using these observations, the interstellar detection of propargyl cyanide (HCCCH$_2$CN) is described, as well as the accompanying laboratory spectroscopy. We discuss these results, and the survey project as a whole, in the context of investigating a previously unexplored reservoir of complex, gas-phase molecules in pre-stellar sources. A series of companion papers describe other new astronomical detections and analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Efficient Production of S$_8$ in Interstellar Ices: The effects of cosmic ray-driven radiation chemistry and non-diffusive bulk reactions

In this work, we reexamine sulfur chemistry occurring on and in the ice mantles of interstellar dust grains, and report the effects of two new modifications to standard astrochemical models; namely, (a) the incorporation of cosmic ray-driven radiation chemistry and (b) the assumption of fast, non-diffusive reactions for key radicals in the bulk. Results from our models of dense molecular clouds show that these changes can have a profound influence on the abundances of sulfur-bearing species in ice mantles, including a reduction in the abundance of solid-phase H$_2$S and HS, and a significant increase in the abundances of OCS, SO$_2$, as well as pure allotropes of sulfur, especially S$_8$. These pure-sulfur species - though nearly impossible to observe directly - have long been speculated to be potential sulfur reservoirs and our results represent possibly the most accurate estimates yet of their abundances in the dense ISM. Moreover, the results of these updated models are found to be in good agreement with available observational data. Finally, we examine the implications of our findings with regard to the as-yet-unknown sulfur reservoir thought to exist in dense interstellar environments.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Investigation of Chemical Differentiation among the NGC2264 Cluster-Forming Clumps

We have carried out mapping observations of molecular emission lines of HC$_{3}$N and CH$_{3}$OH toward two massive cluster-forming clumps, NGC2264-C and NGC2264-D, using the Nobeyama 45-m radio telescope. We derive an $I$(HC$_{3}$N)/$I$(CH$_{3}$OH) integrated intensity ratio map, showing a higher value at clumps including 2MASS point sources at the northern part of NGC2264-D. Possible interpretations of the $I$(HC$_{3}$N)/$I$(CH$_{3}$OH) ratio are discussed. We have also observed molecular emission lines from CCS and N$_{2}$H$^{+}$ toward five positions in each clump. We investigate the $N$(N$_{2}$H$^{+}$)/$N$(CCS) and $N$(N$_{2}$H$^{+}$)/$N$(HC$_{3}$N) column density ratios among the ten positions in order to test whether they can be used as chemical evolutionary indicators in these clumps. The $N$(N$_{2}$H$^{+}$)/$N$(CCS) ratio shows a very high value toward a bright embedded IR source (IRS1), whereas the $N$(N$_{2}$H$^{+}$)/$N$(HC$_{3}$N) ratio at IRS1 is comparable with those at the other positions. These results suggest that UV radiation affects the chemistry around IRS1. We find that there are positive correlations between these column density ratios and the excitation temperatures of N$_{2}$H$^{+}$, which implies the chemical evolution of clumps. These chemical evolutionary indicators likely reflect the combination of evolution along the filamentary structure and evolution of each clump.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

CSO and CARMA Observations of L1157. II. Chemical Complexity in the Shocked Outflow

L1157, a molecular dark cloud with an embedded Class 0 protostar possessing a bipolar outflow, is an excellent source for studying shock chemistry, including grain-surface chemistry prior to shocks, and post-shock, gas-phase processing. The L1157-B1 and B2 positions experienced shocks at an estimated ~2000 and 4000 years ago, respectively. Prior to these shock events, temperatures were too low for most complex organic molecules to undergo thermal desorption. Thus, the shocks should have liberated these molecules from the ice grain-surfaces en masse, evidenced by prior observations of SiO and multiple grain mantle species commonly associated with shocks. Grain species, such as OCS, CH3OH, and HNCO, all peak at different positions relative to species that are preferably formed in higher velocity shocks or repeatedly-shocked material, such as SiO and HCN. Here, we present high spatial resolution (~3") maps of CH3OH, HNCO, HCN, and HCO+ in the southern portion of the outflow containing B1 and B2, as observed with CARMA. The HNCO maps are the first interferometric observations of this species in L1157. The maps show distinct differences in the chemistry within the various shocked regions in L1157B. This is further supported through constraints of the molecular abundances using the non-LTE code RADEX (Van der Tak et al. 2007). We find the east/west chemical differentiation in C2 may be explained by the contrast of the shock's interaction with either cold, pristine material or warm, previously-shocked gas, as seen in enhanced HCN abundances. In addition, the enhancement of the HNCO abundance toward the the older shock, B2, suggests the importance of high-temperature O-chemistry in shocked regions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

First detection of gas-phase methanol in a protoplanetary disk

The first detection of gas-phase methanol in a protoplanetary disk (TW Hya) is presented. In addition to being one of the largest molecules detected in disks to date, methanol is also the first disk organic molecule with an unambiguous ice chemistry origin. The stacked methanol emission, as observed with ALMA, is spectrally resolved and detected across six velocity channels ($>3 σ$), reaching a peak signal-to-noise of $5.5σ$, with the kinematic pattern expected for TW~Hya. Using an appropriate disk model, a fractional abundance of $3\times 10^{-12} - 4 \times 10^{-11}$ (with respect to H$_2$) reproduces the stacked line profile and channel maps, with the favoured abundance dependent upon the assumed vertical location (midplane versus molecular layer). The peak emission is offset from the source position suggesting that the methanol emission has a ring-like morphology: the analysis here suggests it peaks at $\approx 30$~AU reaching a column density $\approx 3-6\times10^{12}$~cm$^{-2}$. In the case of TW Hya, the larger (up to mm-sized) grains, residing in the inner 50~AU, may thus host the bulk of the disk ice reservoir. The successful detection of cold gas-phase methanol in a protoplanetary disk implies that the products of ice chemistry can be explored in disks, opening a window to studying complex organic chemistry during planetary system formation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the Inference of the Cosmic Ray Ionization Rate $ζ$ from the HCO$^+$-to-DCO$^+$ Abundance Ratio: The Effect of Nuclear Spin

The chemistry of dense interstellar regions was analyzed using a time-dependent gas-grain astrochemical simulation and a new chemical network that incorporates deuterated chemistry taking into account nuclear spin-states for the hydrogen chemistry and its deuterated isotopologues. With this new network, the utility of the [HCO$^+$]/[DCO$^+$] abundance ratio as a probe of the cosmic ray ionization rate has been reexamined, with special attention paid to the effect of the initial value of the molecular hydrogen ortho-to-para ratio (OPR). After discussing the use of the probe for cold cores, we then compare our results with previous theoretical and observational results for a molecular cloud close to the supernova remnant W51C, which is thought to have an enhanced cosmic ray ionization rate $ζ$ caused by the nearby $γ$-ray source. In addition, we attempt to use our approach to estimate the cosmic ray ionization rate for L1174, a dense core with an embedded star. Beyond the previously known sensitivity of [HCO$^+$]/[DCO$^+$] to $ζ$, we demonstrate its additional dependence on the initial OPR and, secondarily, on the age of the source, its temperature, and its density. We conclude that the usefulness of the [HCO$^+$]/[DCO$^+$] abundance ratio to constrain the cosmic ray ionization rate in dense regions increases with source age and ionization rate as the ratio becomes far less sensitive to the initial value of the OPR.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The ortho-to-para ratio of interstellar NH$_2$: Quasi-classical trajectory calculations and new simulations

Based on recent $Herschel$ results, the ortho-to-para ratio (OPR) of NH$_2$ has been measured towards the following high-mass star-forming regions: W31C (G10.6-0.4), W49N (G43.2-0.1), W51 (G49.5-0.4), and G34.3+0.1. The OPR at thermal equilibrium ranges from the statistical limit of three at high temperatures to infinity as the temperature tends toward zero, unlike the case of H$_{2}$. Depending on the position observed along the lines-of-sight, the OPR was found to lie either slightly below the high temperature limit of three (in the range $2.2-2.9$) or above this limit ($\sim3.5$, $\gtrsim 4.2$, and $\gtrsim 5.0$). In low temperature interstellar gas, where the H$_{2}$ is para-enriched, our nearly pure gas-phase astrochemical models with nuclear-spin chemistry can account for anomalously low observed NH$_2$-OPR values. We have tentatively explained OPR values larger than three by assuming that spin thermalization of NH$_2$ can proceed at least partially by H-atom exchange collisions with atomic hydrogen, thus increasing the OPR with decreasing temperature. In this paper, we present quasi-classical trajectory calculations of the H-exchange reaction NH$_2$ + H, which show the reaction to proceed without a barrier, confirming that the H-exchange will be efficient in the temperature range of interest. With the inclusion of this process, our models suggest both that OPR values below three arise in regions with temperatures $\gtrsim20-25$~K, depending on time, and values above three but lower than the thermal limit arise at still lower temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Theory of Gas Phase Scattering and Reactivity for Astrochemistry

Because of the very peculiar conditions of chemistry in many astrophysical gases (low densities, mostly low temperatures, kinetics-dominated chemical evolution), great efforts have been devoted to study molecular signatures and chemical evolution. While experiments are being performed in many laboratories, it appears that the efforts directed towards theoretical works are not as strong. This report deals with the present status of chemical physics/physical chemistry theory, for the qualitative and quantitative understanding of kinetics of molecular scattering, being it reactive or inelastic. By gathering several types of expertise, from applied mathematics to physical chemistry, dialog is made possible, as a step towards new and more adapted theoretical frameworks, capable of meeting the theoretical, methodological and numerical challenges of kinetics-dominated gas phase chemistry in astrophysical environments. A state of the art panorama is presented, alongside present-day strengths and shortcomings. However, coverage is not complete, being limited in this report to actual attendance of the workshop. Some paths towards relevant progress are proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Unified Microscopic-Macroscopic Monte Carlo Simulations of Complex Organic Molecule Chemistry in Cold Cores

The recent discovery of methyl formate and dimethyl ether in the gas phase of cold cores with temperatures as cold as 10 K challenges our previous astrochemical models concerning the formation of complex organic molecules. The strong correlation between the abundances and distributions of methyl formate and dimethyl ether further shows that current astrochemical models may be missing important chemical processes in cold astronomical sources. We investigate a scenario in which complex organic molecules and the methoxy radical can be formed on dust grains via a so-called "chain reaction" mechanism, in a similar manner to CO$_2$. A unified gas-grain microscopic-macroscopic Monte Carlo approach with both normal and interstitial sites for icy grain mantles is used to perform the chemical simulations. Reactive desorption with varying degrees of efficiency is included to enhance the non-thermal desorption of species formed on cold dust grains. In addition, varying degrees of efficiency for the surface formation of methoxy are also included. The observed abundances of a variety of organic molecules in cold cores can be reproduced in our models. The strong correlation between the abundances of methyl formate and dimethyl ether in cold cores can also be explained. Non-diffusive chemical reactions on dust grain surfaces may play a key role in the formation of some complex organic molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Herschel observations of EXtra-Ordinary Sources: Analysis of the HIFI 1.2 THz Wide Spectral Survey Toward Orion KL II. Chemical Implications

We present chemical implications arising from spectral models fit to the Herschel/HIFI spectral survey toward the Orion Kleinmann-Low nebula (Orion KL). We focus our discussion on the eight complex organics detected within the HIFI survey utilizing a novel technique to identify those molecules emitting in the hottest gas. In particular, we find the complex nitrogen bearing species CH$_{3}$CN, C$_{2}$H$_{3}$CN, C$_{2}$H$_{5}$CN, and NH$_{2}$CHO systematically trace hotter gas than the oxygen bearing organics CH$_{3}$OH, C$_{2}$H$_{5}$OH, CH$_{3}$OCH$_{3}$, and CH$_{3}$OCHO, which do not contain nitrogen. If these complex species form predominantly on grain surfaces, this may indicate N-bearing organics are more difficult to remove from grain surfaces than O-bearing species. Another possibility is that hot (T$_{\rm kin}$$\sim$300 K) gas phase chemistry naturally produces higher complex cyanide abundances while suppressing the formation of O-bearing complex organics. We compare our derived rotation temperatures and molecular abundances to chemical models, which include gas-phase and grain surface pathways. Abundances for a majority of the detected complex organics can be reproduced over timescales $\gtrsim$ 10$^{5}$ years, with several species being under predicted by less than 3$σ$. Derived rotation temperatures for most organics, furthermore, agree reasonably well with the predicted temperatures at peak abundance. We also find that sulfur bearing molecules which also contain oxygen (i.e. SO, SO$_{2}$, and OCS) tend to probe the hottest gas toward Orion KL indicating the formation pathways for these species are most efficient at high temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Herschel observations of interstellar chloronium. II - Detections toward G29.96-0.02, W49N, W51, and W3(OH), and determinations of the ortho-to-para and $^{35}$Cl/$^{37}$Cl isotopic ratios

We report additional detections of the chloronium molecular ion, H$_2$Cl$^+$, toward four bright submillimeter continuum sources: G29.96, W49N, W51, and W3(OH). With the use of the HIFI instrument on the Herschel Space Observatory, we observed the $2_{12}-1_{01}$ transition of ortho-H$_2^{35}$Cl$^+$ at 781.627 GHz in absorption toward all four sources. Much of the detected absorption arises in diffuse foreground clouds that are unassociated with the background continuum sources and in which our best estimates of the $N({\rm H_2Cl^+})/N({\rm H})$ ratio lie in the range $(0.9 - 4.8) \times 10^{-9}$. These chloronium abundances relative to atomic hydrogen can exceed the predictions of current astrochemical models by up to a factor of 5. Toward W49N, we have also detected the $2_{12}-1_{01}$ transition of ortho-H$_2^{37}$Cl$^+$ at 780.053 GHz and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ at 485.418 GHz. These observations imply $\rm H_2^{35}Cl^+/H_2^{37}Cl^+$ column density ratios that are consistent with the solar system $^{35}$Cl/$^{37}$Cl isotopic ratio of 3.1, and chloronium ortho-to-para ratios consistent with 3, the ratio of spin statistical weights.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Water deuteration and ortho-to-para nuclear spin ratio of H2 in molecular clouds formed via the accumulation of HI gas

We investigate the water deuteration ratio and ortho-to-para nuclear spin ratio of H2 (OPR(H2)) during the formation and early evolution of a molecular cloud, following the scenario that accretion flows sweep and accumulate HI gas to form molecular clouds. We follow the physical evolution of post-shock materials using a one-dimensional shock model, with post-processing gas-ice chemistry simulations. This approach allows us to study the evolution of the OPR(H2) and water deuteration ratio without an arbitrary assumption concerning the initial molecular abundances, including the initial OPR(H2). When the conversion of hydrogen into H2 is almost complete, the OPR(H2) is already much smaller than the statistical value of three due to the spin conversion in the gas phase. As the gas accumulates, the OPR(H2) decreases in a non-equilibrium manner. We find that water ice can be deuterium-poor at the end of its main formation stage in the cloud, compared to water vapor observed in the vicinity of low-mass protostars where water ice is likely sublimated. If this is the case, the enrichment of deuterium in water should mostly occur at somewhat later evolutionary stages of star formation, i.e., cold prestellar/protostellar cores. The main mechanism to suppress water ice deuteration in the cloud is the cycle of photodissociation and reformation of water ice, which efficiently removes deuterium from water ice chemistry. The removal efficiency depends on the main formation pathway of water ice. The OPR(H2) plays a minor role in water ice deuteration at the main formation stage of water ice.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Multi-Transition Study of Molecules toward NGC 1068 based on High-Resolution Imaging Observations with ALMA

We present 0.8-mm band molecular images and spectra obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) toward one of the nearest galaxies with an active galactic nucleus (AGN), NGC 1068. Distributions of CO isotopic species ($^{13}$CO and C$^{18}$O) $\it{J}$ = 3--2, CN $\it{N}$ = 3--2 and CS $\it{J}$ = 7--6 are observed toward the circumnuclear disk (CND) and a part of the starburst ring with an angular resolution of $\sim$1.$^{\prime\prime}$3 $\times$ 1.$^{\prime\prime}$2. The physical properties of these molecules and shock-related molecules such as HNCO, CH$_{3}$CN, SO, and CH$_{3}$OH detected in the 3-mm band were estimated using rotation diagrams under the assumption of local thermodynamic equilibrium. The rotational temperatures of the CO isotopic species and the shock-related molecules in the CND are, respectively, 14--22 K and upper limits of 20--40 K. Although the column densities of the CO isotopic species in the CND are only from one-fifth to one-third of that in the starburst ring, those of the shock-related molecules are enhanced by a factor of 3--10 in the CND. We also discuss the chemistry of each species, and compare the fractional abundances in the CND and starburst ring with those of Galactic sources such as cold cores, hot cores, and shocked molecular clouds in order to study the overall characteristics. We find that the abundances of shock-related molecules are more similar to abundances in hot cores and/or shocked clouds than to cold cores. The CND hosts relatively complex molecules, which are often associated with shocked molecular clouds or hot cores. Because a high X-ray flux can dissociate these molecules, they must also reside in regions shielded from X-rays.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A New and Simple Approach to Determine the Abundance of Hydrogen Molecules on Interstellar Ice Mantles

Water is usually the main component of ice mantles, which cover the cores of dust grains in cold portions of dense interstellar clouds. When molecular hydrogen is adsorbed onto an icy mantle through physisorption, a common assumption in gas-grain rate equation models is to use an adsorption energy for molecular hydrogen on a pure water substrate. However, at high density and low temperature, when H2 is efficiently adsorbed onto the mantle, its surface abundance can be strongly overestimated if this assumption is still used. Unfortunately, the more detailed microscopic Monte Carlo treatment cannot be used to study the abundance of H2 in ice mantles if a full gas-grain network is utilized. We present a numerical method adapted for rate-equation models that takes into account the possibility that an H2 molecule can, while diffusing on the surface, find itself bound to another hydrogen molecule, with a far weaker bond than the H2-water bond, which can lead to more efficient desorption. We label the ensuing desorption "encounter desorption". The method is implemented first in a simple system consisting only of hydrogen molecules at steady state between gas and dust using the rate-equation approach and comparing the results with the results of a microscopic Monte Carlo calculation. We then discuss the use of the rate-equation approach with encounter desorption embedded in a complete gas-grain chemical network. For both systems, the rate-equation model with encounter desorption reproduces the H2 granular coverage computed by the microscopic Monte Carlo model. The method is especially useful for dense and cold environments, and for time-dependent physical conditions, such as occur in the collapse of dense cores and the formation of protoplanetary disks. It is not significantly CPU time consuming, so can be used for example with complex 3D chemical-hydrodynamical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Complex organic molecules along the accretion flow in isolated and externally irradiated protoplanetary disks

(Abridged) The birth environment of the Sun will have influenced the conditions in the pre-solar nebula, including the attainable chemical complexity, important for prebiotic chemistry. The formation and distribution of complex organic molecules (COMs) in a disk around a T Tauri star is investigated for two scenarios: (i) an isolated disk, and (ii) a disk irradiated externally by a nearby massive star. The chemistry is calculated along the accretion flow from the outer disk inwards using a comprehensive network. Two simulations are performed, one beginning with complex ices and one with simple ices only. For the isolated disk, COMs are transported without major alteration into the inner disk where they thermally desorb into the gas reaching an abundance representative of the initial assumed ice abundance. For simple ices, COMs efficiently form on grain surfaces under the conditions in the outer disk. Gas-phase COMs are released into the molecular layer via photodesorption. For the irradiated disk, complex ices are also transported inwards; however, they undergo thermal processing caused by the warmer conditions in the irradiated disk which tends to reduce their abundance along the accretion flow. For simple ices, grain-surface chemistry cannot synthesise COMs in the outer disk because the necessary grain-surface radicals, which tend to be particularly volatile, are not sufficiently abundant on the grain surfaces. Gas-phase COMs are formed in the inner region of the irradiated disk via gas-phase chemistry induced by the desorption of strongly bound molecules such as methanol; hence, the abundances are not representative of the initial molecular abundances injected into the outer disk. These results suggest that the composition of comets formed in isolated disks may differ from those formed in externally irradiated disks with the latter composed of more simple ices.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Complex organic molecules in protoplanetary disks

(Abridged) Protoplanetary disks are vital objects in star and planet formation, possessing all the material which may form a planetary system orbiting the new star. We investigate the synthesis of complex organic molecules (COMs) in disks to constrain the achievable chemical complexity and predict species and transitions which may be observable with ALMA. We have coupled a 2D model of a protoplanetary disk around a T Tauri star with a gas-grain chemical network including COMs. We compare compare synthesised line intensities and calculated column densities with observations and determine those COMs which may be observable in future. COMs are efficiently formed in the disk midplane via grain-surface chemical reactions, reaching peak grain-surface fractional abundances 1e-6 - 1e-4 that of the H nuclei number density. COMs formed on grain surfaces are returned to the gas phase via non-thermal desorption; however, gas-phase species reach lower fractional abundances than their grain-surface equivalents, 1e-12 - 1e-7. Including the irradiation of grain mantle material helps build further complexity in the ice through the replenishment of grain-surface radicals which take part in further grain-surface reactions. There is reasonable agreement with several line transitions of H2CO observed towards several T Tauri star-disk systems. The synthesised line intensities for CH3OH are consistent with upper limits determined towards all sources. Our models suggest CH3OH should be readily observable in nearby protoplanetary disks with ALMA; however, detection of more complex species may prove challenging. Our grain-surface abundances are consistent with those derived from cometary comae observations providing additional evidence for the hypothesis that comets (and other planetesimals) formed via the coagulation of icy grains in the Sun's natal disk.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Distributions of molecules in the circumnuclear disk and surrounding starburst ring in the Seyfert galaxy NGC 1068 observed with ALMA

Sensitive observations with ALMA allow astronomers to observe the detailed distributions of molecules with relatively weak intensity in nearby galaxies. In particular, we report distributions of several molecular transitions including shock and dust related species ($^{13}$CO $J$ = 1--0, C$^{18}$O $J$ = 1--0, $^{13}$CN $N$ = 1--0, CS $J$ = 2--1, SO $J_N$ = 3$_2$--2$_1$, HNCO $J_{Ka,Kc}$ = 5$_{0,5}$--4$_{0,4}$, HC$_3$N $J$ = 11--10, 12--11, CH$_3$OH $J_K$ = 2$_K$--1$_K$, and CH$_3$CN $J_K$ = 6$_K$--5$_K$) in the nearby Seyfert 2 galaxy NGC 1068 observed with the ALMA early science program. The central $\sim$1 arcmin ($\sim$4.3 kpc) of this galaxy was observed in the 100 GHz region covering $\sim$96--100 GHz and $\sim$108--111 GHz with an angular resolution of $\sim4"\times2"$ (290 pc$\times$140 pc) to study the effects of an active galactic nucleus and its surrounding starburst ring on molecular abundances. Here, we present images and report a classification of molecular distributions into three main categories: (1) Molecules concentrated in the circumnuclear disk (CND) (SO $J_N$ = 3$_2$--2$_1$, HC$_3$N $J$ = 11--10, 12--11, and CH$_3$CN $J_K$ = 6$_K$--5$_K$), (2) Molecules distributed both in the CND and the starburst ring (CS $J$ = 2--1 and CH$_3$OH $J_K$ = 2$_K$--1$_K$), (3) Molecules distributed mainly in the starburst ring ($^{13}$CO $J$ = 1--0 and C$^{18}$O $J$ = 1--0). Since most of the molecules such as HC$_3$N observed in the CND are easily dissociated by UV photons and X-rays, our results indicate that these molecules must be effectively shielded. In the starburst ring, the relative intensity of methanol at each clumpy region is not consistent with those of $^{13}$CO, C$^{18}$O, and CS. This difference is probably caused by the unique formation and destruction mechanisms of CH$_3$OH.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interstellar Simulations Using A Unified Microscopic-Macroscopic Monte Carlo Model with a full Gas-Grain Network including Bulk Diffusion in Ice Mantles

We have designed an improved algorithm that enables us to simulate the chemistry of cold dense interstellar clouds with a full gas-grain reaction network. The chemistry is treated by a unified microscopic-macroscopic Monte Carlo approach that includes photon penetration and bulk diffusion. To determine the significance of these two processes, we simulate the chemistry with three different models. In Model 1, we use an exponential treatment to follow how photons penetrate and photodissociate ice species throughout the grain mantle. Moreover, the products of photodissociation are allowed to diffuse via bulk diffusion and react within the ice mantle. Model 2 is similar to Model 1 but with a slower bulk diffusion rate. A reference Model 0, which only allows photodissociation reactions to occur on the top two layers, is also simulated. Photodesorption is assumed to occur from the top two layers in all three models. We found that the abundances of major stable species in grain mantles do not differ much among these three models, and the results of our simulation for the abundances of these species agree well with observations. Likewise, the abundances of gas-phase species in the three models do not vary. However, the abundances of radicals in grain mantles can differ by up to two orders of magnitude depending upon the degree of photon penetration and the bulk diffusion of photodissociation products. We also found that complex molecules can be formed at temperatures as low as 10 K in all three models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The HNC/HCN Ratio in Star-Forming Regions

HNC and HCN, typically used as dense gas tracers in molecular clouds, are a pair of isomers that have great potential as a temperature probe because of temperature dependent, isomer-specific formation and destruction pathways. Previous observations of the HNC/HCN abundance ratio show that the ratio decreases with increasing temperature, something that standard astrochemical models cannot reproduce. We have undertaken a detailed parameter study on which environmental characteristics and chemical reactions affect the HNC/HCN ratio and can thus contribute to the observed dependence. Using existing gas and gas-grain models updated with new reactions and reaction barriers, we find that in static models the H + HNC gas-phase reaction regulates the HNC/HCN ratio under all conditions, except for very early times. We quantitively constrain the combinations of H abundance and H + HNC reaction barrier that can explain the observed HNC/HCN temperature dependence and discuss the implications in light of new quantum chemical calculations. In warm-up models, gas-grain chemistry contributes significantly to the predicted HNC/HCN ratio and understanding the dynamics of star formation is therefore key to model the HNC/HCN system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Water deuterium fractionation in the high-mass star-forming region G34.26+0.15 based on Herschel/HIFI data

Understanding water deuterium fractionation is important for constraining the mechanisms of water formation in interstellar clouds. Observations of HDO and H$_2^{18}$O transitions were carried out towards the high-mass star-forming region G34.26+0.15 with the HIFI instrument onboard the Herschel Space Observatory, as well as with ground-based single-dish telescopes. Ten HDO lines and three H$_2^{18}$O lines covering a broad range of upper energy levels (22-204 K) were detected. We used a non-LTE 1D analysis to determine the HDO/H$_2$O ratio as a function of radius in the envelope. Models with different water abundance distributions were considered in order to reproduce the observed line profiles. The HDO/H$_2$O ratio is found to be lower in the hot core ($\sim$3.5 $\times$ 10$^{-4}$ - 7.5 $\times$ 10$^{-4}$) than in the colder envelope ($\sim$1.0 $\times$ 10$^{-3}$ - 2.2 $\times$ 10$^{-3}$). This is the first time that a radial variation of the HDO/H$_2$O ratio has been found to occur in a high-mass source. The chemical evolution of this source was modeled as a function of its radius and the observations are relatively well reproduced. The comparison between the chemical model and the observations leads to an age of $\sim$10$^5$ years after the infrared dark cloud stage.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Astrochemistry: Synthesis and Modelling

We discuss models that astrochemists have developed to study the chemical composition of the interstellar medium. These models aim at computing the evolution of the chemical composition of a mixture of gas and dust under as- trophysical conditions. These conditions, as well as the geometry and the physical dynamics, have to be adapted to the objects being studied because different classes of objects have very different characteristics (temperatures, densities, UV radia- tion fields, geometry, history etc); e.g., proto-planetary disks do not have the same characteristics as protostellar envelopes. Chemical models are being improved continually thanks to comparisons with observations but also thanks to laboratory and theoretical work in which the individual processes are studied.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Deep observations of O2 toward a low-mass protostar with Herschel-HIFI

According to traditional gas-phase chemical models, O2 should be abundant in molecular clouds, but until recently, attempts to detect interstellar O2 line emission with ground- and space-based observatories have failed. Following the multi-line detections of O2 with low abundances in the Orion and rho Oph A molecular clouds with Herschel, it is important to investigate other environments, and we here quantify the O2 abundance near a solar-mass protostar. Observations of O2, at 487 GHz toward a deeply embedded low-mass Class 0 protostar, NGC 1333-IRAS 4A, are presented, using the HIFI instrument on the Herschel Space Observatory. Complementary data of the chemically related NO and CO molecules are obtained as well. The high spectral resolution data are analysed using radiative transfer models to infer column densities and abundances, and are tested directly against full gas-grain chemical models. The deep HIFI spectrum fails to show O2 at the velocity of the dense protostellar envelope, implying one of the lowest abundance upper limits of O2/H2 at <6x10^-9 (3 sigma). However, a tentative (4.5 sigma) detection of O2 is seen at the velocity of the surrounding NGC 1333 molecular cloud, shifted by 1 km/s relative to the protostar. For the protostellar envelope, pure gas-phase models and gas-grain chemical models require a long pre-collapse phase (~0.7-1x10^6 years), during which atomic and molecular oxygen are frozen out onto dust grains and fully converted to H2O, to avoid overproduction of O2 in the dense envelope. The same model also reproduces the limits on the chemically related NO molecule. The tentative detection of O2 in the surrounding cloud is consistent with a low-density PDR model with small changes in reaction rates. The low O2 abundance in the collapsing envelope around a low-mass protostar suggests that the gas and ice entering protoplanetary disks is very poor in O2.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Interstellar water chemistry: from laboratory to observations

Water is observed throughout the universe, from diffuse interstellar clouds to protoplanetary disks around young stars, and from comets in our own solar system and exoplanetary atmospheres to galaxies at high redshifts. This review summarizes the spectroscopy and excitation of water in interstellar space as well as the basic chemical processes that form and destroy water under interstellar conditions. Three major routes to water formation are identified: low temperature ion-molecule chemistry, high-temperature neutral-neutral chemistry and gas-ice chemistry. The rate coefficients of several important processes entering the networks are discussed in detail; several of them have been determined only in the last decade through laboratory experiments and theoretical calculations. Astronomical examples of each of the different chemical routes are presented using data from powerful new telescopes, in particular the Herschel Space Observatory. Basic chemical physics studies remain critically important to analyze astronomical data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modeling The Molecular Composition in an AGN Disk

We use a high-temperature chemical network to derive the molecular abundances in axisymmetric accretion disk models around active galactic nuclei (AGNs) within 100 pc using simple radial and vertical density and temperature distributions motivated by more detailed physical models. We explore the effects of X-ray irradiation and cosmic ray ionization on the spatial distribution of the molecular abundances of CO, CN, CS, HCN, HCO+, HC3N, C2H, and c-C3H2 using a variety of plausible disk structures. These simple models have molecular regions with a layer of X-ray dominated regions, a midplane without the strong influence of X-rays, and a high-temperature region in the inner portion with moderate X-ray flux where families of polyynes (C$_{\rm n}$H$_{2}$) and cyanopolyynes can be enhanced. For the high midplane density disks we explore, we find that cosmic rays produced by supernovae do not significantly affect the regions unless the star formation efficiency significantly exceeds that of the Milky Way. We highlight molecular abundance observations and ratios that may distinguish among theoretical models of the density distribution in AGN disks. Finally, we assess the importance of the shock crossing time and the accretion time relative to the formation time for various chemical species. Vertical column densities are tabulated for a number of molecular species at both the characteristic shock crossing time and steady state. Although we do not attempt to fit any particular system or set of observations, we discuss our models and results in the context of the nearby AGN NGC 1068.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Reactive Desorption and Radiative Association as Possible Drivers of Complex Molecule Formation in the Cold Interstellar Medium

The recent discovery of terrestrial-type organic species such as methyl formate and dimethyl ether in the cold interstellar gas has proved that the formation of organic matter in the Galaxy begins at a much earlier stage of star formation than was thought before. This discovery represents a challenge for astrochemical modelers. The abundances of these molecules cannot be explained by the previously developed "warm-up" scenario, in which organic molecules are formed via diffusive chemistry on surfaces of interstellar grains starting at 30 K, and then released to the gas at higher temperatures during later stages of star formation. In this article, we investigate an alternative scenario in which complex organic species are formed via a sequence of gas-phase reactions between precursor species formed on grain surfaces and then ejected into the gas via efficient reactive desorption, a process in which non-thermal desorption occurs as a result of conversion of the exothermicity of chemical reactions into the ejection of products from the surface. The proposed scenario leads to reasonable if somewhat mixed results at temperatures as low as 10 K and may be considered as a step towards the explanation of abundances of terrestrial-like organic species observed during the earliest stages of star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Unified Microscopic-Macroscopic Monte Carlo Simulation of Gas-Grain Chemistry in Cold Dense Interstellar Clouds

For the first time, we report a unified microscopic-macroscopic Monte Carlo simulation of gas-grain chemistry in cold interstellar clouds in which both the gas-phase and the grain surface chemistry are simulated by a stochastic technique. The surface chemistry is simulated with a microscopic Monte Carlo method in which the chemistry occurs on an initially flat surface. The surface chemical network consists of 29 reactions initiated by the accreting species H, O, C, and CO. Four different models are run with diverse but homogeneous physical conditions including temperature, gas density, and diffusion-barrier-to-desorption energy ratio. As time increases, icy interstellar mantles begin to grow. Our approach allows us to determine the morphology of the ice, layer by layer, as a function of time, and to ascertain the environment or environments for individual molecules. Our calculated abundances can be compared with observations of ices and gas-phase species, as well as the results of other models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

From Prestellar to Protostellar Cores II. Time Dependence and Deuterium Fractionation

We investigate the molecular evolution and D/H abundance ratios that develop as star formation proceeds from a dense-cloud core to a protostellar core, by solving a gas-grain reaction network applied to a 1-D radiative hydrodynamic model with infalling fluid parcels. Spatial distributions of gas and ice-mantle species are calculated at the first-core stage, and at times after the birth of a protostar. Gas-phase methanol and methane are more abundant than CO at radii $r\lesssim 100$ AU in the first-core stage, but gradually decrease with time, while abundances of larger organic species increase. The warm-up phase, when complex organic molecules are efficiently formed, is longer-lived for those fluid parcels in-falling at later stages. The formation of unsaturated carbon chains (warm carbon-chain chemistry) is also more effective in later stages; C$^+$, which reacts with CH$_4$ to form carbon chains, increases in abundance as the envelope density decreases. The large organic molecules and carbon chains are strongly deuterated, mainly due to high D/H ratios in the parent molecules, determined in the cold phase. We also extend our model to simulate simply the chemistry in circumstellar disks, by suspending the 1-D infall of a fluid parcel at constant disk radii. The species CH$_3$OCH$_3$ and HCOOCH$_3$ increase in abundance in $10^4-10^5$ yr at the fixed warm temperature; both also have high D/H ratios.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Carbon-Chain Species in Warm-up Models

In previous warm-up chemical models of the low-mass star-forming region L1527, we investigated the evolution of carbon-chain unsaturated hydrocarbon species when the envelope temperature is slightly elevated to $T\approx 30$ K. These models demonstrated that enhanced abundances of such species can be explained by gas-phase ion-molecule chemistry following the partial sublimation of methane from grain surfaces. We also concluded that the abundances of hydrocarbon radicals such as the C$_{\rm n}$H family should be further enhanced as the temperatures increase to higher values, but this conclusion stood in contrast with the lack of unambiguous detection of these species toward hot core and corino sources. Meanwhile, observational surveys have identified C$_2$H, C$_4$H, CH$_3$CCH, and CH$_3$OH toward hot corinos (especially IRAS 16293-2422) as well as towards L1527, with lower abundances for the carbon chain radicals and higher abundances for the other two species toward the hot corinos. In addition, the {\it Herschel Space Telescope} has detected the bare linear chain C$_3$ in 50 K material surrounding young high-mass stellar objects. To understand these new results, we revisit previous warm-up models with an augmented gas-grain network that incorporated reactions from a gas-phase network that was constructed for use with increased temperature up to 800 K. Some of the newly adopted reactions between carbon-chain species and abundant H$_2$ possess chemical activation energy barriers. The revised model results now better reproduce the observed abundances of unsaturated carbon chains under hot-corino (100 K) conditions and make predictions for the abundances of bare carbon chains in the 50 K regions observed by Herschel HIFI.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Propagation of Low-Energy Cosmic Rays in Molecular Clouds: Calculations in Two Dimensions

We calculate cosmic ray transport with a collisional Boltzmann Transport Equation, including E-M forces. We apply the Crank-Nicholson Method to solve this equation. At each time step, the spatial distribution of cosmic rays is applied to the ZEUS 2D magnetohydrodynamics model, which is then utilized to calculate the resulting E-M field. Finally, the field is applied to the initial equation. This sequence is repeated over many time steps until a steady state is reached. We consider results from t = 0 until steady state for an isotropic low energy cosmic ray flux, and also for an enhanced cosmic ray flux impinging on one side of a molecular cloud. This cosmic ray flux is used to determine an ionization rate of interstellar hydrogen by cosmic rays, zeta. Astrochemical implications are briefly mentioned.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Interstellar OH+, H2O+ and H3O+ along the sight-line to G10.6-0.4

We report the detection of absorption lines by the reactive ions OH+, H2O+ and H3O+ along the line of sight to the submillimeter continuum source G10.6$-$0.4 (W31C). We used the Herschel HIFI instrument in dual beam switch mode to observe the ground state rotational transitions of OH+ at 971 GHz, H2O+ at 1115 and 607 GHz, and H3O+ at 984 GHz. The resultant spectra show deep absorption over a broad velocity range that originates in the interstellar matter along the line of sight to G10.6$-$0.4 as well as in the molecular gas directly associated with that source. The OH+ spectrum reaches saturation over most velocities corresponding to the foreground gas, while the opacity of the H2O+ lines remains lower than 1 in the same velocity range, and the H3O+ line shows only weak absorption. For LSR velocities between 7 and 50 kms$^{-1}$ we estimate total column densities of $N$(OH+) $> 2.5 \times 10^{14}$ cm$^{-2}$, $N$(H2O+) $\sim 6 \times 10^{13}$ cm$^{-2}$ and $N$(H3O+) $\sim 4.0 \times 10^{13}$ cm$^{-2}$. These detections confirm the role of O$^+$ and OH$^+$ in initiating the oxygen chemistry in diffuse molecular gas and strengthen our understanding of the gas phase production of water. The high ratio of the OH+ by the H2O+ column density implies that these species predominantly trace low-density gas with a small fraction of hydrogen in molecular form.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Laboratory Astrophysics and the State of Astronomy and Astrophysics

Laboratory astrophysics and complementary theoretical calculations are the foundations of astronomy and astrophysics and will remain so into the foreseeable future. The impact of laboratory astrophysics ranges from the scientific conception stage for ground-based, airborne, and space-based observatories, all the way through to the scientific return of these projects and missions. It is our understanding of the under-lying physical processes and the measurements of critical physical parameters that allows us to address fundamental questions in astronomy and astrophysics. In this regard, laboratory astrophysics is much like detector and instrument development at NASA, NSF, and DOE. These efforts are necessary for the success of astronomical research being funded by the agencies. Without concomitant efforts in all three directions (observational facilities, detector/instrument development, and laboratory astrophysics) the future progress of astronomy and astrophysics is imperiled. In addition, new developments in experimental technologies have allowed laboratory studies to take on a new role as some questions which previously could only be studied theoretically can now be addressed directly in the lab. With this in mind we, the members of the AAS Working Group on Laboratory Astrophysics, have prepared this State of the Profession Position Paper on the laboratory astrophysics infrastructure needed to ensure the advancement of astronomy and astrophysics in the next decade.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

New Discoveries in Cosmology and Fundamental Physics through Advances in Laboratory Astrophysics

As the Cosmology and Fundamental Physics (CFP) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the early universe, the microwave background, the reionization and galaxy formation up to virialization of protogalaxies, large scale structure, the intergalactic medium, the determination of cosmological parameters, dark matter, dark energy, tests of gravity, astronomically determined physical constants, and high energy physics using astronomical messengers. Central to the progress in these areas are the corresponding advances in laboratory astrophysics which are required for fully realizing the CFP scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics which produce the observed astrophysical processes. The 5 areas of laboratory astrophysics which we have identified as relevant to the CFP panel are atomic, molecular, plasma, nuclear, and particle physics. Here, Section 2 describes some of the new scientific opportunities and compelling scientific themes which will be enabled by advances in laboratory astrophysics. In Section 3, we provide the scientific context for these opportunities. Section 4 briefly discusses some of the experimental and theoretical advances in laboratory astrophysics required to realize the CFP scientific opportunities of the next decade. As requested in the Call for White Papers, Section 5 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 6.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

New Discoveries in Galaxies across Cosmic Time through Advances in Laboratory Astrophysics

As the Galaxies across Cosmic Time (GCT) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the formation, evolution, and global properties of galaxies and galaxy clusters, as well as active galactic nuclei and QSOs, mergers, star formation rate, gas accretion, and supermassive black holes. Central to the progress in these areas are the corresponding advances in laboratory astrophysics that are required for fully realizing the GCT scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics that produce the observed astrophysical processes. The 5 areas of laboratory astrophysics that we have identified as relevant to the CFP panel are atomic, molecular, solid matter, plasma, nuclear, and particle physics. In this white paper, we describe in Section 2 some of the new scientific opportunities and compelling scientific themes that will be enabled by advances in laboratory astrophysics. In Section 3, we provide the scientific context for these opportunities. Section 4 briefly discusses some of the experimental and theoretical advances in laboratory astrophysics required to realize the GCT scientific opportunities of the next decade. As requested in the Call for White Papers, Section 5 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 6.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

New Discoveries in Stars and Stellar Evolution through Advances in Laboratory Astrophysics

As the Stars and Stellar Evolution (SSE) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the Sun as a star, stellar astrophysics, the structure and evolution of single and multiple stars, compact objects, SNe, gamma-ray bursts, solar neutrinos, and extreme physics on stellar scales. Central to the progress in these areas are the corresponding advances in laboratory astrophysics, required to fully realize the SSE scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics that produces the observed astrophysical processes. The 6 areas of laboratory astrophysics, which we have identified as relevant to the CFP panel, are atomic, molecular, solid matter, plasma, nuclear physics, and particle physics. In this white paper, we describe in Section 2 the scientific context and some of the new scientific opportunities and compelling scientific themes which will be enabled by advances in laboratory astrophysics. In Section 3, we discuss some of the experimental and theoretical advances in laboratory astrophysics required to realize the SSE scientific opportunities of the next decade. As requested in the Call for White Papers, Section 4 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 5.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Modeling Carbon Chain Anions in L1527

The low-mass protostellar region L1527 is unusual because it contains observable abundances of unsaturated carbon-chain molecules including CnH radicals, H2Cn carbenes, cyanopolyynes, and the negative ions C4H- and C6H-, all of which are more associated with cold cores than with protostellar regions. Sakai et al. suggested that these molecules are formed in L1527 from the chemical precursor methane, which evaporates from the grains during the heat-up of the region. With the gas-phase osu.03.2008 network extended to include negative ions of the families Cn-, and CnH-, as well as the newly detected C3N-, we modeled the chemistry that occurs following methane evaporation at T~ 25-30 K. We are able to reproduce most of the observed molecular abundances in L1527 at a time of ~5000 yr. At later times, the overall abundance of anions become greater than that of electrons, which has an impact on many organic species and ions. The anion-to-neutral ratio in our calculation is in good agreement with observation for C6H- but exceeds the observed ratio by more than three orders of magnitude for C4H-. In order to explain this difference, further investigation is needed on the rate coefficients for electron attachment and other reactions regarding anions.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Chemical sensitivity to the ratio of the cosmic-ray ionization rates of He and H2 in dense clouds

Aim: To determine whether or not gas-phase chemical models with homogeneous and time-independent physical conditions explain the many observed molecular abundances in astrophysical sources, it is crucial to estimate the uncertainties in the calculated abundances and compare them with the observed abundances and their uncertainties. Non linear amplification of the error and bifurcation may limit the applicability of chemical models. Here we study such effects on dense cloud chemistry. Method: Using a previously studied approach to uncertainties based on the representation of rate coefficient errors as log normal distributions, we attempted to apply our approach using as input a variety of different elemental abundances from those studied previously. In this approach, all rate coefficients are varied randomly within their log normal (Gaussian) distribution, and the time-dependent chemistry calculated anew many times so as to obtain good statistics for the uncertainties in the calculated abundances. Results: Starting with so-called ``high-metal'' elemental abundances, we found bimodal rather than Gaussian like distributions for the abundances of many species and traced these strange distributions to an extreme sensitivity of the system to changes in the ratio of the cosmic ray ionization rate zeta\_He for He and that for molecular hydrogen zeta\_H2. The sensitivity can be so extreme as to cause a region of bistability, which was subsequently found to be more extensive for another choice of elemental abundances. To the best of our knowledge, the bistable solutions found in this way are the same as found previously by other authors, but it is best to think of the ratio zeta\_He/zeta\_H2 as a control parameter perpendicular to the ''standard'' control parameter zeta/n\_H.

0 citations0 reviewsNo code linkedOpen access