Showing 1–38 of 38 results for author: Willacy, K

Chemistry in the GG Tau A Disk: Constraints from H2D+, N2H+, and DCO+ High Angular Resolution ALMA Observations

Authors: Parashmoni Kashyap, Liton Majumdar, Anne Dutrey, Stéphane Guilloteau, Karen Willacy, Edwige Chapillon, Richard Teague, Dmitry Semenov, Thomas Henning, Neal Turner, Raghvendra Sahai, Ágnes Kóspál, Audrey Coutens, V. Piétu, Pierre Gratier, Maxime Ruaud, N. T. Phuong, E. Di Folco, Chin-Fei Lee, Y. -W. Tang

Abstract: Resolved molecular line observations are essential for gaining insight into the physical and chemical structure of protoplanetary disks, particularly in cold, dense regions where planets form and acquire compositions. However, tracing these regions is challenging because most molecules freeze onto grain surfaces and are not observable in the gas phase. We investigated cold molecular chemistry in t… ▽ More Resolved molecular line observations are essential for gaining insight into the physical and chemical structure of protoplanetary disks, particularly in cold, dense regions where planets form and acquire compositions. However, tracing these regions is challenging because most molecules freeze onto grain surfaces and are not observable in the gas phase. We investigated cold molecular chemistry in the triple stellar T Tauri disk GG Tau A, which harbours a massive gas and dust ring and an outer disk, using ALMA Band 7 observations. We present high angular resolution maps of N2H+ and DCO+ emission, with upper limits reported for H2D+, 13CS, and SO2. The radial intensity profile of N2H+ shows most emission near the ring outer edge, while DCO+ exhibits double peaks, one near the ring inner edge and the other in the outer disk. With complementary observations of lower-lying transitions, we constrained the molecular surface densities and rotation temperatures. We compared the derived quantities with model predictions across different cosmic ray ionization (CRI) rates, carbon-to-oxygen (C/O) ratios, and stellar UV fluxes. Cold molecular chemistry, affecting N2H+, DCO+, and H2D+ abundances, is most sensitive to CRI rates, while stellar UV flux and C/O ratios have minimal impact on these three ions. Our best model requires a low cosmic ray ionization rate of 1e-18 s-1. However, it fails to match the low temperatures derived from N2H+ and DCO+, 12 to 16 K, which are much lower than the CO freezing temperature. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: Resubmitted to ApJ. Thanks to the two anonymous referees for their constructive comments

The Carbon Isotopic Ratio and Planet Formation

Authors: Edwin A. Bergin, Arthur Bosman, Richard Teague, Jenny Calahan, Karen Willacy, L. Ilsedore Cleeves, Kamber Schwarz, Ke Zhang, Simon Bruderer

Abstract: We present the first detection of 13CCH in a protoplanetary disk (TW Hya). Using observations of C2H we measure CCH/13CCH = 65 +/- 20 in gas with a CO isotopic ratio of 12CO/13CO = 21 +/- 5 (Yoshida et al. 2022a). The TW Hya disk exhibits a gas phase C/O that exceeds unity and C2H is the tracer of this excess carbon. We confirm that the TW Hya gaseous disk exhibits two separate carbon isotopic res… ▽ More We present the first detection of 13CCH in a protoplanetary disk (TW Hya). Using observations of C2H we measure CCH/13CCH = 65 +/- 20 in gas with a CO isotopic ratio of 12CO/13CO = 21 +/- 5 (Yoshida et al. 2022a). The TW Hya disk exhibits a gas phase C/O that exceeds unity and C2H is the tracer of this excess carbon. We confirm that the TW Hya gaseous disk exhibits two separate carbon isotopic reservoirs as noted previously (Yoshida et al. 2022a). We explore two theoretical solutions for the development of this dichotomy. One model represents TW Hya today with a protoplanetary disk exposed to a cosmic ray ionization rate that is below interstellar as consistent with current estimates. We find that this model does not have sufficient ionization in cold (T < 40 K) layers to activate carbon isotopic fractionation. The second model investigates a younger TW Hya protostellar disk exposed to an interstellar cosmic ray ionization rate. We find that the younger model has sources of ionization deeper in a colder disk that generates two independent isotopic reservoirs. One reservoir is 12C-enriched carried by methane/hydrocarbon ices and the other is 13C-enriched carried by gaseous CO. The former potentially provides a source of methane/hydrocarbon ices to power the chemistry that generates the anomalously strong C$_2$H emission in this (and other) disk systems in later stages. The latter provides a source of gaseous 13C-rich material to generate isotopic enrichments in forming giant planets as recently detected in the super-Jupiter TYC 8998-760-1 b by Zhang et al. (2021). △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: 16 pages, 8 figures, accepted by the Astrophysical Journal

Retrievals of Protoplanetary Disk Parameters using Thermochemical Models: I. Disk Gas Mass from Hydrogen Deuteride Spectroscopy

Authors: Young Min Seo, Karen Willacy, Geoffrey Bryden, Dariusz C. Lis, Paul F. Goldsmith, Klaus M. Pontoppidan, Wing-Fai Thi

Abstract: We discuss statistical relationships between the mass of protoplanetary disks and the hydrogen deuteride (HD) line emission and the dust spectral energy distribution (SED) determined using 3000 ProDiMo disk models. The models have 15 free parameters describing disk physical properties, the central star, and the local radiation field. The sampling of physical parameters is done using a Monte Carlo… ▽ More We discuss statistical relationships between the mass of protoplanetary disks and the hydrogen deuteride (HD) line emission and the dust spectral energy distribution (SED) determined using 3000 ProDiMo disk models. The models have 15 free parameters describing disk physical properties, the central star, and the local radiation field. The sampling of physical parameters is done using a Monte Carlo approach to evaluate the probability density functions of observables as a function of physical parameters. We find that the HD fractional abundance is almost constant even though the UV flux varies by several orders of magnitude. Probing the statistical relation between the physical quantities and the HD flux, we find that low-mass (optically thin) disks display a tight correlation between the average disk gas temperature and HD line flux, while massive disks show no such correlation. We demonstrate that the central star luminosity, disk size, dust size distribution, and HD flux may be used to determine the disk gas mass to within a factor of three. We also find that the far-IR and sub-mm/mm SEDs and the HD flux may serve as strong constraints for determining the disk gas mass to within a factor of two. If the HD lines are fully spectrally resolved ($R\gtrsim 1.5\times10^6, Δv=0.2~\rm km\,s^{-1}$), the 56 $μ$m and 112 $μ$m HD line profiles alone may constrain the disk gas mass to within a factor of two. △ Less

Submitted 24 April, 2024; v1 submitted 24 October, 2023; originally announced October 2023.

Comments: 31 pages, 16 figures, ApJ accepted, preprint

Statistical Prediction of [CII] Observations by Constructing Probability Density Functions using SOFIA, Herschel, and Spitzer Observations

Authors: Young Min Seo, Karen Willacy, Umaa Rebbapragada

Abstract: We present a statistical algorithm for predicting the [CII] emission from Herschel and Spitzer continuum images using probability density functions between the [CII] emission and continuum emission. The [CII] emission at 158 $μ$m is a critical tracer in studying the life cycle of interstellar medium and galaxy evolution. Unfortunately, its frequency is in the far infrared (FIR), which is opaque th… ▽ More We present a statistical algorithm for predicting the [CII] emission from Herschel and Spitzer continuum images using probability density functions between the [CII] emission and continuum emission. The [CII] emission at 158 $μ$m is a critical tracer in studying the life cycle of interstellar medium and galaxy evolution. Unfortunately, its frequency is in the far infrared (FIR), which is opaque through the troposphere and cannot be observed from the ground except for highly red-shifted sources (z $\gtrsim$ 2). Typically [CII] observations of closer regions have been carried out using suborbital or space observatories. Given the high cost of these facilities and limited time availability, it is important to have highly efficient observations/operations in terms of maximizing science returns. This requires accurate prediction of the strength of emission lines and, therefore, the time required for their observation. However, [CII] emission has been hard to predict due to a lack of strong correlations with other observables. Here we adopt a new approach to making accurate predictions of [CII] emission by relating this emission simultaneously to several tracers of dust emission in the same region. This is done using a statistical methodology utilizing probability density functions (PDFs) among [CII] emission and Spitzer IRAC and Herschel PACS/SPIRE images. Our test result toward a star-forming region, RCW 120, demonstrates that our methodology delivers high-quality predictions with less than 30\% uncertainties over 80\% of the entire observation area, which is more than sufficient to test observation feasibility and maximize science return. The {\it pickle} dump files storing the PDFs and trained neural network module are accessible upon request and will support future far-infrared missions, for example, GUSTO and FIR Probe. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 7 figures

Journal ref: Journal of Astronomical Instrumentation, 2023, 10.1142/S2251171723500071

Vertical distribution of cyclopropenylidene and propadiene in the atmosphere of Titan

Authors: Karen Willacy, SiHe Chen, Danica J. Adams, Yuk L. Yung

Abstract: Titan's atmosphere is a natural laboratory for exploring the photochemical synthesis of organic molecules. Significant recent advances in the study of the atmosphere of Titan include: (a) detection of C$_3$ molecules: C$_3$H$_6$, CH$_2$CCH$_2$, c-C$_3$H$_2$, and (b) retrieval of C$_6$H$_6$, which is formed primarily via C$_3$ chemistry, from Cassini-UVIS data. The detection of $c$-C$_3$H$_2$ is of… ▽ More Titan's atmosphere is a natural laboratory for exploring the photochemical synthesis of organic molecules. Significant recent advances in the study of the atmosphere of Titan include: (a) detection of C$_3$ molecules: C$_3$H$_6$, CH$_2$CCH$_2$, c-C$_3$H$_2$, and (b) retrieval of C$_6$H$_6$, which is formed primarily via C$_3$ chemistry, from Cassini-UVIS data. The detection of $c$-C$_3$H$_2$ is of particular significance since ring molecules are of great astrobiological importance. Using the Caltech/JPL KINETICS code, along with the best available photochemical rate coefficients and parameterized vertical transport, we are able to account for the recent observations. It is significant that ion chemistry, reminiscent of that in the interstellar medium, plays a major role in the production of c-C$_3$H$_2$ above 1000 km. △ Less

Submitted 27 April, 2022; originally announced April 2022.

Comments: 18 pages (+ 27 page appendix), 16 figures

Comets in context: Comparing comet compositions with protosolar nebula models

Authors: Karen Willacy, Neal Turner, Boncho Bonev, Erika Gibb, Neil Dello Russo, Michael DiSanti, Ronald J. Vervack Jr., Nathan X. Roth

Abstract: Comets provide a valuable window into the chemical and physical conditions at the time of their formation in the young solar system. We seek insights into where and when these objects formed by comparing the range of abundances observed for nine molecules and their average values across a sample of 29 comets to the predicted midplane ice abundances from models of the protosolar nebula. Our fiducia… ▽ More Comets provide a valuable window into the chemical and physical conditions at the time of their formation in the young solar system. We seek insights into where and when these objects formed by comparing the range of abundances observed for nine molecules and their average values across a sample of 29 comets to the predicted midplane ice abundances from models of the protosolar nebula. Our fiducial model, where ices are inherited from the interstellar medium, can account for the observed mixing ratio ranges of each molecule considered, but no single location or time reproduces the abundances of all molecules simultaneously. This suggests that each comet consists of material processed under a range of conditions. In contrast, a model where the initial composition of disk material is `reset', wiping out any previous chemical history, cannot account for the complete range of abundances observed in comets. Using toy models that combine material processed under different thermal conditions we find that a combination of warm (CO-poor) and cold (CO-rich) material is required to account for both the average properties of the Jupiter-family and Oort cloud comets, and the individual comets we consider. This could occur by the transport (either radial or vertical) of ice-coated dust grains in the early solar system. Comparison of the models to the average Jupiter-family and Oort cloud comet compositions suggest the two families formed in overlapping regions of the disk, in agreement with the findings of A'Hearn et al. (2012) and with the predictions of the Nice model (Gomes et al. 2005, Tsiganis et al. 2005). △ Less

Submitted 15 April, 2022; originally announced April 2022.

Comments: 32 pages, 8 figures

Linking atmospheric chemistry of the hot Jupiter HD 209458b to its formation location through infrared transmission and emission spectra

Authors: Spandan Dash, Liton Majumdar, Karen Willacy, Shang-Min Tsai, Neal Turner, P. B. Rimmer, Murthy S. Gudipati, Wladimir Lyra, Anil Bhardwaj

Abstract: The elemental ratios of carbon, nitrogen, and oxygen in the atmospheres of hot Jupiters may hold clues to their formation locations in the protostellar disc. In this work, we adopt gas phase chemical abundances of C, N and O from several locations in a disc chemical kinetics model as sources for the envelope of the hot Jupiter HD 209458b and evolve the planet's atmospheric composition using a 1D c… ▽ More The elemental ratios of carbon, nitrogen, and oxygen in the atmospheres of hot Jupiters may hold clues to their formation locations in the protostellar disc. In this work, we adopt gas phase chemical abundances of C, N and O from several locations in a disc chemical kinetics model as sources for the envelope of the hot Jupiter HD 209458b and evolve the planet's atmospheric composition using a 1D chemical kinetics model, treating both vertical mixing and photochemistry. We consider two atmospheric pressure-temperature profiles, one with and one without a thermal inversion. From each of the resulting 32 atmospheric composition profiles, we find that the molecules CH4, NH3, HCN, and C2H2 are more prominent in the atmospheres computed using a realistic non-inverted P-T profile in comparison to a prior equilibrium chemistry based work which used an analytical P-T profile. We also compute the synthetic transmission and emission spectra for these atmospheres and find that many spectral features vary with the location in the disc where the planet's envelope was accreted. By comparing with the species detected using the latest high-resolution ground-based observations, our model suggests HD 209458b could have accreted most of its gas between the CO2 and CH4 icelines with a super solar C/O ratio from its protostellar disc, which in turn directly inherited its chemical abundances from the protostellar cloud. Finally, we simulate observing the planet with the James Webb Space Telescope (JWST) and show that differences in spectral signatures of key species can be recognized. Our study demonstrates the enormous importance of JWST in providing new insights into hot Jupiter's formation environments. △ Less

Submitted 11 May, 2022; v1 submitted 8 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ApJ, 28 Pages, 11 figures, 2 Tables (Appendix: 5 Figures, 1 Table)

Non-ideal MHD simulations of subcritical prestellar cores with non-equilibrium chemistry

Authors: Aris Tritsis, Christoph Federrath, Karen Willacy, Konstantinos Tassis

Abstract: Non-ideal magnetohydrodynamic (MHD) effects are thought to be gravity's closest ally in overcoming the support of magnetic fields and in forming stars. Here, we modify the publicly available version of the adaptive mesh refinement code FLASH (Fryxell et al. 2000; Dubey et al. 2008) to include a detailed treatment of non-ideal MHD and study such effects in collapsing prestellar cores. We implement… ▽ More Non-ideal magnetohydrodynamic (MHD) effects are thought to be gravity's closest ally in overcoming the support of magnetic fields and in forming stars. Here, we modify the publicly available version of the adaptive mesh refinement code FLASH (Fryxell et al. 2000; Dubey et al. 2008) to include a detailed treatment of non-ideal MHD and study such effects in collapsing prestellar cores. We implement two very extended non-equilibrium chemical networks, the largest of which is comprised of $\sim$ 300 species and includes a detailed description of deuterium chemistry. The ambipolar-diffusion, Ohmic and Hall resistivities are then self-consistently calculated from the abundances of charged species. We present a series of 2-dimensional axisymmetric simulations where we vary the chemical model, cosmic-ray ionization rate, and grain distribution. We benchmark our implementation against ideal MHD simulations and previously-published results. We show that, at high densities ($n_{\rm{H_2}}>~10^6~\rm{cm^{-3}}$), the ion that carries most of the perpendicular and parallel conductivities is not $\rm{H_3^+}$ as was previously thought, but is instead $\rm{D_3^+}$. △ Less

Submitted 21 December, 2021; originally announced December 2021.

Comments: 18 pages, 14 figures. Accepted for publication in MNRAS

Physical and Chemical Structure of the Disk and Envelope of the Class 0/I Protostar L1527

Authors: Lizxandra Flores-Rivera, Susan Terebey, Karen Willacy, Andrea Isella, Neal Turner, Mario Flock

Abstract: Sub-millimeter spectral line and continuum emission from the protoplanetary disks and envelopes of protostars are powerful probes of their structure, chemistry, and dynamics. Here we present a benchmark study of our modeling code, RadChemT, that for the first time uses a chemical model to reproduce ALMA C$^{18}$O (2-1) and CARMA $^{12}$CO (1-0) and N$_{2}$H$^{+}$ (1-0) observations of L1527, that… ▽ More Sub-millimeter spectral line and continuum emission from the protoplanetary disks and envelopes of protostars are powerful probes of their structure, chemistry, and dynamics. Here we present a benchmark study of our modeling code, RadChemT, that for the first time uses a chemical model to reproduce ALMA C$^{18}$O (2-1) and CARMA $^{12}$CO (1-0) and N$_{2}$H$^{+}$ (1-0) observations of L1527, that allow us to distinguish the disk, the infalling envelope and outflow of this Class 0/I protostar. RadChemT combines dynamics, radiative transfer, gas chemistry and gas-grain reactions to generate models which can be directly compared with observations for individual protostars. Rather than individually fit abundances to a large number of free parameters, we aim to best match the spectral line maps by (i) adopting a physical model based on density structure and luminosity derived primarily from previous work that fit SED and 2D imaging data, updating it to include a narrow jet detected in CARMA and ALMA data near ($\leq 75$au) the protostar, and then (ii) computing the resulting astrochemical abundances for 292 chemical species. Our model reproduces the C$^{18}$O and N$_{2}$H$^{+}$ line strengths within a factor of 3.0; this is encouraging considering the pronounced abundance variation (factor $> 10^3$) between the outflow shell and CO snowline region near the midplane. Further, our modeling confirms suggestions regarding the anti-correlation between N$_{2}$H$^{+}$ and the CO snowline between 400 au to 2,000 au from the central star. Our modeling tools represent a new and powerful capability with which to exploit the richness of spectral line imaging provided by modern submillimeter interferometers. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Comments: 25 pages, 9 figures, and 7 tables. Accepted by ApJ

An Ammonia Spectral Map of the L1495-B218 Filaments in the Taurus Molecular Cloud: II CCS & HC$_7$N Chemistry and Three Modes of Star Formation in the Filaments

Authors: Young Min Seo, Liton Majumdar, Paul F. Goldsmith, Yancy L. Shirley, Karen Willacy, Derek Ward-Thompson, Rachel Friesen, David Frayer, Sarah E. Church, Dongwoo Chung, Kieran Cleary, Nichol Cunningham, Kiruthika Devaraj, Dennis Egan, Todd Gaier, Rohit Gawande, Joshua O. Gundersen, Andrew I. Harris, Pekka Kangaslahti, Anthony C. S. Readhead, Lorene Samoska, Matthew Sieth, Michael Stennes, Patricia Voll, Steve White

Abstract: We present deep CCS and HC$_7$N observations of the L1495-B218 filaments in the Taurus molecular cloud obtained using the K-band focal plane array on the 100m Green Bank Telescope. We observed the L1495-B218 filaments in CCS $J_N$ = 2$_1$$-$1$_0$ and HC$_7$N $J$ = 21$-$20 with a spectral resolution of 0.038 km s$^{-1}$ and an angular resolution of 31$''$. We observed strong CCS emission in both ev… ▽ More We present deep CCS and HC$_7$N observations of the L1495-B218 filaments in the Taurus molecular cloud obtained using the K-band focal plane array on the 100m Green Bank Telescope. We observed the L1495-B218 filaments in CCS $J_N$ = 2$_1$$-$1$_0$ and HC$_7$N $J$ = 21$-$20 with a spectral resolution of 0.038 km s$^{-1}$ and an angular resolution of 31$''$. We observed strong CCS emission in both evolved and young regions and weak emission in two evolved regions. HC$_7$N emission is observed only in L1495A-N and L1521D. We find that CCS and HC$_7$N intensity peaks do not coincide with NH$_3$ or dust continuum intensity peaks. We also find that the fractional abundance of CCS does not show a clear correlation with the dynamical evolutionary stage of dense cores. Our findings and chemical modeling indicate that the fractional abundances of CCS and HC$_7$N are sensitive to the initial gas-phase C/O ratio, and they are good tracers of young condensed gas only when the initial C/O is close to solar value. Kinematic analysis using multiple lines including NH$_3$, HC$_7$N, CCS, CO, HCN, \& HCO$^+$ suggests that there may be three different star formation modes in the L1495-B218 filaments. At the hub of the filaments, L1495A/B7N has formed a stellar cluster with large-scale inward flows (fast mode), while L1521D, a core embedded in a filament, is slowly contracting due to its self-gravity (slow mode). There is also one isolated core that appears to be marginally stable and may undergo quasi-static evolution (isolated mode). △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: ApJ accepted

Ionization-Driven Depletion and Redistribution of CO in Protoplanetary Disks

Authors: Sarah E. Dodson-Robinson, Neal J. Evans II, Alyssa Ramos, Mo Yu, Karen Willacy

Abstract: Based on the interstellar CO/H2 ratio, carbon monoxide-based censuses of protoplanetary disks in Lupus, sigma Orionis, and Chamaeleon I found no disks more massive than the minimum-mass solar nebula, which is inconsistent with the existence of exoplanets more massive than Jupiter. Observations and models are converging on the idea that ionization-driven chemistry depletes carbon monoxide in T-Taur… ▽ More Based on the interstellar CO/H2 ratio, carbon monoxide-based censuses of protoplanetary disks in Lupus, sigma Orionis, and Chamaeleon I found no disks more massive than the minimum-mass solar nebula, which is inconsistent with the existence of exoplanets more massive than Jupiter. Observations and models are converging on the idea that ionization-driven chemistry depletes carbon monoxide in T-Tauri disks. Yet the extent of CO depletion depends on the incident flux of ionizing radiation, and some T-Tauri stars may have winds strong enough to shield their disks from cosmic rays. There is also a range of X-ray luminosities possible for a given stellar mass. Here we use a suite of chemical models, each with a different incident X-ray or cosmic-ray flux, to assess whether CO depletion is a typical outcome for T-Tauri disks. We find that CO dissociation in the outer disk is a robust result for realistic ionization rates, with abundance reductions between 70% and 99.99% over 2 Myr of evolution. Furthermore, after the initial dissociation epoch, the inner disk shows some recovery of the CO abundance from CO2 dissociation. In highly ionized disks, CO recovery in the inner disk combined with depletion in the outer disk creates a centrally peaked CO abundance distribution. The emitting area in rare CO isotopologues may be an indirect ionization indicator: in a cluster of disks with similar ages, those with the most compact CO isotopologue emission see the highest ionization rates. △ Less

Submitted 14 November, 2018; originally announced November 2018.

Comments: 8 pages, 4 figures; accepted for publication in ApJ Letters

Buildup of Abiotic Oxygen and Ozone in Moist Atmospheres of Temperate Terrestrial Exoplanets and its Impact on the Spectral Fingerprint in Transit Observations

Authors: Armin Kleinböhl, Karen Willacy, A. James Friedson, Pin Chen, Mark R. Swain

Abstract: We investigate the abiotic production of oxygen and its photochemical byproduct ozone through water vapor photolysis in moist atmospheres of temperate terrestrial exoplanets. The amount of water vapor available for photolysis in the middle atmosphere of a planet can be limited by an atmospheric cold-trap, the formation of which largely depends on the amount of non-condensable gases. We study this… ▽ More We investigate the abiotic production of oxygen and its photochemical byproduct ozone through water vapor photolysis in moist atmospheres of temperate terrestrial exoplanets. The amount of water vapor available for photolysis in the middle atmosphere of a planet can be limited by an atmospheric cold-trap, the formation of which largely depends on the amount of non-condensable gases. We study this effect using a photochemical model coupled to a 1D radiative-convective equilibrium model in atmospheres with N2, CO2 and H2O as the main constituents. We find that in atmospheres with a low N2 inventory, water vapor mixing ratios in the middle atmosphere can be over two orders of magnitude higher compared to atmospheres with an Earth-like N2 inventory. Without a strong surface sink, the non-condensable oxygen can build up rapidly, drying out the upper atmosphere. With a moderate surface sink, the planet can approach a steady state with significant oxygen mixing ratios in which oxygen production is balanced by surface uptake. We use a radiative transfer model to study the spectroscopic fingerprint of these atmospheres in transit observations. Spectral signatures of abiotic oxygen and ozone can be of comparable magnitude as in spectra of Earth seen as an exoplanet. Middle atmospheric water vapor is unlikely to be a usable indicator of the abiotic origin of oxygen because of the influence of oxygen on the water vapor distribution. This suggests that atmospheric oxygen and ozone cannot be used as binary bioindicators and their interpretation will likely require atmospheric and planetary models. △ Less

Submitted 25 May, 2018; originally announced May 2018.

Comments: 16 pages, 10 figures

Detection of HOCO+ in the protostar IRAS 16293-2422

Authors: Liton Majumdar, Pierre Gratier, Valentine Wakelam, Emmanuel Caux, Karen Willacy, Michael E. Ressler

Abstract: The protonated form of CO2, HOCO+, is assumed to be an indirect tracer of CO2 in the millimeter/submillimeter regime since CO2 lacks a permanent dipole moment. Here, we report the detection of two rotational emission lines (4 0,4-3 0,3) and (5 0,5-4 0,4) of HOCO+ in IRAS 16293-2422. For our observations, we have used EMIR heterodyne 3 mm receiver of the IRAM 30m telescope. The observed abundance o… ▽ More The protonated form of CO2, HOCO+, is assumed to be an indirect tracer of CO2 in the millimeter/submillimeter regime since CO2 lacks a permanent dipole moment. Here, we report the detection of two rotational emission lines (4 0,4-3 0,3) and (5 0,5-4 0,4) of HOCO+ in IRAS 16293-2422. For our observations, we have used EMIR heterodyne 3 mm receiver of the IRAM 30m telescope. The observed abundance of HOCO+ is compared with the simulations using the 3-phase NAUTILUS chemical model. Implications of the measured abundances of HOCO+ to study the chemistry of CO2 ices using JWST-MIRI and NIRSpec are discussed as well. △ Less

Submitted 14 March, 2018; originally announced March 2018.

Comments: Accepted for publication in MNRAS on March 13th, 7 pages, 3 figures, 2 tables

The Effects of Protostellar Disk Turbulence on CO Emission Lines: A Comparison Study of Disks with Constant CO Abundance vs. Chemically Evolving Disks

Authors: Mo Yu, Neal J. Evans, Sarah E. Dodson-Robinson, Karen Willacy, Neal J. Turner

Abstract: Turbulence is the leading candidate for angular momentum transport in protoplanetary disks and therefore influences disk lifetimes and planet formation timescales. However, the turbulent properties of protoplanetary disks are poorly constrained observationally. Recent studies have found turbulent speeds smaller than what fully-developed MRI would produce (Flaherty et al. 2015, 2017). However, exis… ▽ More Turbulence is the leading candidate for angular momentum transport in protoplanetary disks and therefore influences disk lifetimes and planet formation timescales. However, the turbulent properties of protoplanetary disks are poorly constrained observationally. Recent studies have found turbulent speeds smaller than what fully-developed MRI would produce (Flaherty et al. 2015, 2017). However, existing studies assumed a constant CO/H2 ratio of 0.0001 in locations where CO is not frozen-out or photo-dissociated. Our previous studies of evolving disk chemistry indicate that CO is depleted by incorporation into complex organic molecules well inside the freeze-out radius of CO. We consider the effects of this chemical depletion on measurements of turbulence. Simon et al. (2015) suggested that the ratio of the peak line flux to the flux at line center of the CO J=3-2 transition is a reasonable diagnostic of turbulence, so we focus on that metric, while adding some analysis of the more complex effects on spatial distribution. We simulate the emission lines of CO based on chemical evolution models presented in Yu et al. (2016), and find that the peak-to-trough ratio changes as a function of time as CO is destroyed. Specifically, a CO-depleted disk with high turbulent velocity mimics the peak-to-trough ratios of a non-CO-depleted disk with lower turbulent velocity. We suggest that disk observers and modelers take into account the possibility of CO depletion when using line peak-to-trough ratios to constrain the degree of turbulence in disks. Assuming that CO/H2 = 0.0001 at all disk radii can lead to underestimates of turbulent speeds in the disk by at least 0.2 km/s. △ Less

Submitted 7 November, 2017; v1 submitted 10 October, 2017; originally announced October 2017.

Comments: 12 pages, including 10 figures. Accepted for publication in Astrophysical Journal. Revised following comments from community

Quantifying the Impact of Spectral Coverage on the Retrieval of Molecular Abundances from Exoplanet Transmission Spectra

Authors: John W. Chapman, Robert T. Zellem, Michael R. Line, Geoff Bryden, Karen Willacy, Aishwarya R. Iyer, Gautam Vasisht, Jacob Bean, Nicolas B. Cowan, Jonathan J. Fortney, Caitlin A. Griffith, Tiffany Kataria, Eliza M. -R. Kempton, Laura Kreidberg, Julianne I. Moses, Kevin B. Stevenson, Mark R. Swain

Abstract: Using forward models for representative exoplanet atmospheres and a radiometric instrument model, we have generated synthetic observational data to explore how well the major C- and O-bearing chemical species (CO, CO2, CH4, and H2O), important for determining atmospheric opacity and radiation balance, can be constrained by transit measurements as a function of spectral wavelength coverage. This wo… ▽ More Using forward models for representative exoplanet atmospheres and a radiometric instrument model, we have generated synthetic observational data to explore how well the major C- and O-bearing chemical species (CO, CO2, CH4, and H2O), important for determining atmospheric opacity and radiation balance, can be constrained by transit measurements as a function of spectral wavelength coverage. This work features simulations for a notional transit spectroscopy mission and compares two cases for instrument spectral coverage (wavelength coverage from 0.5-2.5 $μ$m and 0.5-5 $μ$m. The simulation is conducted on a grid with a range of stellar magnitudes and incorporates a full retrieval of atmospheric model parameters. We consider a range of planets from sub-Neptunes to hot Jupiters and include both low and high mean molecular weight atmospheres. We find that including the 2.5-5 $μ$m wavelength range provides a significant improvement in the degree of constraint on the retrieved molecular abundances: up to ~3 orders of magnitude for a low mean molecular weight atmosphere ($μ$=2.3) and up to a factor of ~6 for a high mean molecular weight atmosphere ($μ$=28). These decreased uncertainties imply that broad spectral coverage between the visible and the mid-infrared is an important tool for understanding the chemistry and composition of exoplanet atmospheres. This analysis suggests that the JWST/NIRSpec 0.6-5 $μ$m prism spectroscopy mode, or similar wavelength coverage with possible future missions, will be an important resource for exoplanet atmospheric characterization. △ Less

Submitted 15 August, 2017; v1 submitted 15 May, 2017; originally announced May 2017.

Comments: Accepted by PASP; 18 pages, 7 figures

Disk masses around solar-mass stars are underestimated by CO observations

Authors: Mo Yu, Neal J. Evans, Sarah E. Dodson-Robinson, Karen Willacy, Neal J. Turner

Abstract: Gas in protostellar disks provides the raw material for giant planet formation and controls the dynamics of the planetesimal-building dust grains. Accurate gas mass measurements help map the observed properties of planet-forming disks onto the formation environments of known exoplanets. Rare isotopologues of carbon monoxide (CO) have been used as gas mass tracers for disks in the Lupus star-formin… ▽ More Gas in protostellar disks provides the raw material for giant planet formation and controls the dynamics of the planetesimal-building dust grains. Accurate gas mass measurements help map the observed properties of planet-forming disks onto the formation environments of known exoplanets. Rare isotopologues of carbon monoxide (CO) have been used as gas mass tracers for disks in the Lupus star-forming region, with an assumed interstellar CO/H$_2$ abundance ratio. Unfortunately, observations of T-Tauri disks show that CO abundance is not interstellar---a finding reproduced by models that show CO abundance decreasing both with distance from the star and as a function of time. Here we present radiative transfer simulations that assess the accuracy of CO-based disk mass measurements. We find that the combination of CO chemical depletion in the outer disk and optically thick emission from the inner disk leads observers to underestimate gas mass by more than an order of magnitude if they use the standard assumptions of interstellar CO/H$_2$ ratio and optically thin emission. Furthermore, CO abundance changes on million-year timescales, introducing an age/mass degeneracy into observations. To reach factor of a few accuracy for CO-based disk mass measurements, we suggest that observers and modelers adopt the following strategies: (1) select the low-$J$ transitions; (2) observe multiple CO isotopologues and use either intensity ratios or normalized line profiles to diagnose CO chemical depletion; and (3) use spatially resolved observations to measure the CO abundance distribution. △ Less

Submitted 18 April, 2017; originally announced April 2017.

Comments: 22 pages, 16 figures, Accepted for publication in ApJ

A new astrobiological model of the atmosphere of Titan

Authors: Karen Willacy, Mark Allen, Yuk Yung

Abstract: We present results of an investigation into the formation of nitrogen-bearing molecules in the atmosphere of Titan. We extend a previous model (Li et al. 2015, 2016) to cover the region below the tropopause, so the new model treats the atmosphere from Titan's surface to an altitude of 1500 km. We consider the effects of condensation and sublimation using a continuous, numerically stable method. Th… ▽ More We present results of an investigation into the formation of nitrogen-bearing molecules in the atmosphere of Titan. We extend a previous model (Li et al. 2015, 2016) to cover the region below the tropopause, so the new model treats the atmosphere from Titan's surface to an altitude of 1500 km. We consider the effects of condensation and sublimation using a continuous, numerically stable method. This is coupled with parameterized treatments of the sedimentation of the aerosols and their condensates, and the formation of haze particles. These processes affect the abundances of heavier species such as the nitrogen-bearing molecules, but have less effect on the abundances of lighter molecules. Removal of molecules to form aerosols also plays a role in determining the mixing ratios, in particular of HNC, HC3N and HCN. We find good agreement with the recently detected mixing ratios of C2H5CN, with condensation playing an important role in determining the abundance of this molecule below 500 km. Of particular interest is the chemistry of acrylonitrile (C2H3CN) which has been suggested by Stevenson et al. (2015) as a molecule that could form biological membranes in an oxygen-deficient environment. With the inclusion of haze formation we find good agreement of our model predictions of acrylonitrile with the available observations. △ Less

Submitted 7 September, 2016; v1 submitted 15 July, 2016; originally announced July 2016.

Comments: 17 pages, 6 figures, Accepted by ApJ

Probing Planet Forming Zones with Rare CO Isotopologues

Authors: Mo Yu, Karen Willacy, Sarah E. Dodson-Robinson, Neal J. Turner, Neal J. Evans II

Abstract: The gas near the midplanes of planet-forming protostellar disks remains largely unprobed by observations due to the high optical depth of commonly observed molecules such as CO and H$_2$O. However, rotational emission lines from rare molecules may have optical depths near unity in the vertical direction, so that the lines are strong enough to be detected, yet remain transparent enough to trace the… ▽ More The gas near the midplanes of planet-forming protostellar disks remains largely unprobed by observations due to the high optical depth of commonly observed molecules such as CO and H$_2$O. However, rotational emission lines from rare molecules may have optical depths near unity in the vertical direction, so that the lines are strong enough to be detected, yet remain transparent enough to trace the disk midplane. Here we present a chemical model of an evolving T-Tauri disk and predict the optical depths of rotational transitions of $^{12}$C$^{16}$O, $^{13}$C$^{16}$O, $^{12}$C$^{17}$O and $^{12}$C$^{18}$O. The MRI-active disk is primarily heated by the central star due to the formation of the dead zone. CO does not freeze out in our modeled region within $70$AU around a sunlike star. However, the abundance of CO decreases because of the formation of complex organic molecules (COM), producing an effect that can be misinterpreted as the "snow line". These results are robust to variations in our assumptions about the evolution of the gas to dust ratio. The optical depths of low-order rotational lines of C$^{17}$O are around unity, making it possible to see into the disk midplane using C$^{17}$O. Combining observations with modeled C$^{17}$O$/$H$_2$ ratios, like those we provide, can yield estimates of protoplanetary disks' gas masses. △ Less

Submitted 29 March, 2016; originally announced March 2016.

Comments: Accepted for publication in ApJ

Chemistry as a diagnostic of prestellar core geometry

Authors: A. Tritsis, K. Tassis, K. Willacy

Abstract: We present a new method for assessing the intrinsic 3D shape of prestellar cores from molecular column densities. We have employed hydrodynamic simulations of contracting, isothermal cores considering three intrinsic geometries: spherical, cylindrical/filamentary and disk-like. We have coupled our hydrodynamic simulations with non-equilibrium chemistry. We find that a) when cores are observed very… ▽ More We present a new method for assessing the intrinsic 3D shape of prestellar cores from molecular column densities. We have employed hydrodynamic simulations of contracting, isothermal cores considering three intrinsic geometries: spherical, cylindrical/filamentary and disk-like. We have coupled our hydrodynamic simulations with non-equilibrium chemistry. We find that a) when cores are observed very elongated (i.e. for aspect ratios $\le$ 0.15) the intrinsic 3D geometry can be probed by their 2D molecular emission maps, since these exhibit significant qualitative morphological differences between cylindrical and disk-like cores. Specifically, if a disk-like core is observed as a filamentary object in dust emission, then it will be observed as two parallel filaments in $\rm{N_2H^{+}}$; b) for cores with higher aspect ratios (i.e. 0.15 $\sim$ 0.9) we define a metric $Δ$ that quantifies whether a molecular column density profile is centrally peaked, depressed or flat. We have identified one molecule ($\rm{CN}$) for which $Δ$ as a function of the aspect ratio probes the 3D geometry of the core; and c) for cores with almost circular projections (i.e. for aspect ratios $\sim$ 1), we have identified three molecules ($\rm{OH}$, $\rm{CO}$ and $\rm{H_2CO}$) that can be used to probe the intrinsic 3D shape by close inspection of their molecular column density radial profiles. We alter the temperature and the cosmic-ray ionization rate and demonstrate that our method is robust against the choice of parameters. △ Less

Submitted 9 February, 2016; originally announced February 2016.

Comments: 14 pages, 14 figures, 2 tables. Accepted for publication by MNRAS

The composition of the protosolar disk and the formation conditions for comets

Authors: K. Willacy, C. Alexander, M. Ali-Dib, C. Ceccarelli, S. B. Charnley, M. Doronin, Y. Ellinger, P. Gast, E. Gibb, S. N. Milam, O. Mousis, E. Pauzat, C. Tornow, E. S. Wirstrom, E Zicler

Abstract: Conditions in the protosolar nebula have left their mark in the composition of cometary volatiles, thought to be some of the most pristine material in the solar system. Cometary compositions represent the end point of processing that began in the parent molecular cloud core and continued through the collapse of that core to form the protosun and the solar nebula, and finally during the evolution o… ▽ More Conditions in the protosolar nebula have left their mark in the composition of cometary volatiles, thought to be some of the most pristine material in the solar system. Cometary compositions represent the end point of processing that began in the parent molecular cloud core and continued through the collapse of that core to form the protosun and the solar nebula, and finally during the evolution of the solar nebula itself as the cometary bodies were accreting. Disentangling the effects of the various epochs on the final composition of a comet is complicated. But comets are not the only source of information about the solar nebula. Protostellar disks around young stars similar to the protosun provide a way of investigating the evolution of disks similar to the solar nebula while they are in the process of evolving to form their own solar systems. In this way we can learn about the physical and chemical conditions under which comets formed, and about the types of dynamical processing that shaped the solar system we see today. This paper summarizes some recent contributions to our understanding of both cometary volatiles and the composition, structure and evolution of protostellar disks. △ Less

Submitted 8 July, 2015; originally announced July 2015.

Comments: To appear in Space Science Reviews. The final publication is available at Springer via http://dx.doi.org/10.1007/s11214-015-0167-6

Effect of OH depletion on measurements of the mass-to-flux ratio in molecular cloud cores

Authors: K. Tassis, K. Willacy, H. W. Yorke, N. J. Turner

Abstract: The ratio of mass and magnetic flux determines the relative importance of magnetic and gravitational forces in the evolution of molecular clouds and their cores. Its measurement is thus central in discriminating between different theories of core formation and evolution. Here we discuss the effect of chemical depletion on measurements of the mass-to-flux ratio using the same molecule (OH) both for… ▽ More The ratio of mass and magnetic flux determines the relative importance of magnetic and gravitational forces in the evolution of molecular clouds and their cores. Its measurement is thus central in discriminating between different theories of core formation and evolution. Here we discuss the effect of chemical depletion on measurements of the mass-to-flux ratio using the same molecule (OH) both for Zeeman measurements of the magnetic field and the determination of the mass of the region. The uncertainties entering through the OH abundance in determining separately the magnetic field and the mass of a region have been recognized in the literature. It has been proposed however that, when comparing two regions of the same cloud, the abundance will in both cases be the same. We show that this assumption is invalid. We demonstrate that when comparing regions with different densities, the effect of OH depletion in measuring changes of the mass-to-flux ratio between different parts of the same cloud can even reverse the direction of the underlying trends (for example, the mass-to-flux ratio may appear to decrease as we move to higher density regions). The systematic errors enter primarily through the inadequate estimation of the mass of the region. △ Less

Submitted 25 August, 2014; originally announced August 2014.

Comments: 5 pages, 3 figures, accepted for publication in MNRAS

A Herschel [C II] Galactic plane survey II: CO-dark H2 in clouds

Authors: W. D. Langer, T. Velusamy, J. L. Pineda, K. Willacy, P. F. Goldsmith

Abstract: ABRIDGED: Context: HI and CO large scale surveys of the Milky Way trace the diffuse atomic clouds and the dense shielded regions of molecular hydrogen clouds. However, until recently, we have not had spectrally resolved C+ surveys to characterize the photon dominated interstellar medium, including, the H2 gas without C, the CO-dark H2, in a large sample of clouds. Aims: To use a sparse Galactic pl… ▽ More ABRIDGED: Context: HI and CO large scale surveys of the Milky Way trace the diffuse atomic clouds and the dense shielded regions of molecular hydrogen clouds. However, until recently, we have not had spectrally resolved C+ surveys to characterize the photon dominated interstellar medium, including, the H2 gas without C, the CO-dark H2, in a large sample of clouds. Aims: To use a sparse Galactic plane survey of the 1.9 THz [C II] spectral line from the Herschel Open Time Key Programme, Galactic Observations of Terahertz C+ (GOT C+), to characterize the H2 gas without CO in a statistically significant sample of clouds. Methods: We identify individual clouds in the inner Galaxy by fitting [CII] and CO isotopologue spectra along each line of sight. We combine these with HI spectra, along with excitation models and cloud models of C+, to determine the column densities and fractional mass of CO-dark H2 clouds. Results: We identify 1804 narrow velocity [CII] interstellar cloud components in different categories. About 840 are diffuse molecular clouds with no CO, 510 are transition clouds containing [CII] and 12CO, but no 13CO, and the remainder are dense molecular clouds containing 13CO emission. The CO-dark H2 clouds are concentrated between Galactic radii 3.5 to 7.5 kpc and the column density of the CO-dark H2 layer varies significantly from cloud-to-cloud with an average 9X10^(20) cm-2. These clouds contain a significant fraction of CO-dark H2 mass, varying from ~75% for diffuse molecular clouds to ~20% for dense molecular clouds. Conclusions: We find a significant fraction of the warm molecular ISM gas is invisible in HI and CO, but is detected in [CII]. The fraction of CO-dark H2 is greatest in the diffuse clouds and decreases with increasing total column density, and is lowest in the massive clouds. △ Less

Submitted 7 January, 2014; v1 submitted 11 December, 2013; originally announced December 2013.

Comments: 21 pages, 19 figures, accepted for publication in A&A (2014)

A search for co-evolving ion and neutral gas species in prestellar molecular cloud cores

Authors: Konstantinos Tassis, Talayeh Hezareh, Karen Willacy

Abstract: Comparison of linewidths of spectral line profiles of ions and neutral molecules have been recently used to estimate the strength of the magnetic field in turbulent star-forming regions. However, the ion (HCO+) and neutral (HCN) species used in such studies may not be necessarily co-evolving at every scale and density and may thus not trace the same regions. Here, we use coupled chemical/dynamical… ▽ More Comparison of linewidths of spectral line profiles of ions and neutral molecules have been recently used to estimate the strength of the magnetic field in turbulent star-forming regions. However, the ion (HCO+) and neutral (HCN) species used in such studies may not be necessarily co-evolving at every scale and density and may thus not trace the same regions. Here, we use coupled chemical/dynamical models of evolving prestellar molecular cloud cores including non-equilibrium chemistry, with and without magnetic fields, to study the spatial distribution of HCO+ and HCN, which have been used in observations of spectral linewidth differences to date. In addition, we seek new ion-neutral pairs that are good candidates for such observations because they have similar evolution and are approximately co-spatial in our models. We identify three such good candidate pairs: HCO+/NO, HCO+/CO, and NO+/NO. △ Less

Submitted 25 September, 2012; originally announced September 2012.

Comments: 9 pages, 5 figures, accepted for publication in ApJ

Non-Equilibrium Chemistry of Dynamically Evolving Prestellar Cores: II. Ionization and Magnetic Field

Authors: Konstantinos Tassis, Karen Willacy, Harold W. Yorke, Neal Turner

Abstract: We study the effect that non-equilibrium chemistry in dynamical models of collapsing molecular cloud cores has on measurements of the magnetic field in these cores, the degree of ionization, and the mean molecular weight of ions. We find that OH and CN, usually used in Zeeman observations of the line-of-sight magnetic field, have an abundance that decreases toward the center of the core much faste… ▽ More We study the effect that non-equilibrium chemistry in dynamical models of collapsing molecular cloud cores has on measurements of the magnetic field in these cores, the degree of ionization, and the mean molecular weight of ions. We find that OH and CN, usually used in Zeeman observations of the line-of-sight magnetic field, have an abundance that decreases toward the center of the core much faster than the density increases. As a result, Zeeman observations tend to sample the outer layers of the core and consistently underestimate the core magnetic field. The degree of ionization follows a complicated dependence on the number density at central densities up to 10^5 cm^{-3} for magnetic models and 10^6 cm^{-3} in non-magnetic models. At higher central densities the scaling approaches a power-law with a slope of -0.6 and a normalization which depends on the cosmic-ray ionization rate ζ and the temperature T as (ζT)^1/2. The mean molecular weight of ions is systematically lower than the usually assumed value of 20 - 30, and, at high densities, approaches a value of 3 due to the asymptotic dominance of the H3+ ion. This significantly lower value implies that ambipolar diffusion operates faster. △ Less

Submitted 22 May, 2012; v1 submitted 17 November, 2011; originally announced November 2011.

Comments: 7 pages, 5 figures, accepted for publication in ApJ

Non-Equilibrium Chemistry of Dynamically Evolving Prestellar Cores: I. Basic Magnetic and Non-Magnetic Models and Parameter Studies

Authors: Konstantinos Tassis, Karen Willacy, Harold W. Yorke, Neal J. Turner

Abstract: We combine dynamical and non-equilibrium chemical modeling of evolving prestellar molecular cloud cores, and explore the evolution of molecular abundances in the contracting core. We model both magnetic cores, with varying degrees of initial magnetic support, and non-magnetic cores, with varying collapse delay times. We explore, through a parameter study, the competing effects of various model par… ▽ More We combine dynamical and non-equilibrium chemical modeling of evolving prestellar molecular cloud cores, and explore the evolution of molecular abundances in the contracting core. We model both magnetic cores, with varying degrees of initial magnetic support, and non-magnetic cores, with varying collapse delay times. We explore, through a parameter study, the competing effects of various model parameters in the evolving molecular abundances, including the elemental C/O ratio, the temperature, and the cosmic-ray ionization rate. We find that different models show their largest quantitative differences at the center of the core, whereas the outer layers, which evolve slower, have abundances which are severely degenerate among different dynamical models. There is a large range of possible abundance values for different models at a fixed evolutionary stage (central density), which demonstrates the large potential of chemical differentiation in prestellar cores. However, degeneracies among different models, compounded with uncertainties induced by other model parameters, make it difficult to discriminate among dynamical models. To address these difficulties, we identify abundance ratios between particular molecules, the measurement of which would have maximal potential for discrimination among the different models examined here. In particular, we find that the ratios between NH3 and CO; NH2 and CO; NH3 and HCO+ are sensitive to the evolutionary timescale, and that the ratio between HCN and OH is sensitive to the C/O ratio. Finally, we demonstrate that measurements of the central deviation (central depletion or enhancement) of abundances of certain molecules are good indicators of the dynamics of the core. △ Less

Submitted 26 April, 2012; v1 submitted 16 November, 2011; originally announced November 2011.

Comments: 20 pages, 15 figures, accepted for publication in ApJ

The Influence of Deuteration and Turbulent Diffusion on the Observed D/H Ratio

Authors: T. A. Bell, K. Willacy, T. G. Phillips, M. Allen, D. C. Lis

Abstract: The influence of turbulent mixing on the chemistry of the interstellar medium has so far received little attention. Previous studies of this effect have suggested that it might play an important role in mixing the various phases of the interstellar medium. In this paper we examine the potential effects of turbulent diffusion on the deuterium chemistry within molecular clouds. We find that such mix… ▽ More The influence of turbulent mixing on the chemistry of the interstellar medium has so far received little attention. Previous studies of this effect have suggested that it might play an important role in mixing the various phases of the interstellar medium. In this paper we examine the potential effects of turbulent diffusion on the deuterium chemistry within molecular clouds. We find that such mixing acts to reduce the efficiency of deuteration in these clouds by increasing the ionization fraction and reducing freeze-out of heavy molecules. This leads to lower abundances for many deuterated species. We also examine the influence of turbulent mixing on the transition from atomic hydrogen to H2 and from atomic deuterium to HD near the cloud edge. We find that including turbulent diffusion in our models serves to push these transitions deeper into the cloud and helps maintain a higher atomic fraction throughout the cloud envelope. Based on these findings, we propose a new process to account for the significant scatter in the observed atomic D/H ratio for galactic sightlines extending beyond the Local Bubble. Although several mechanisms have been put forward to explain this scatter, they are unable to fully account for the range in D/H values. We suggest a scenario in which turbulent mixing of atomic and molecular gas at the edges of molecular clouds causes the observed atomic D/H ratio to vary by a factor of ~2. △ Less

Submitted 14 February, 2011; originally announced February 2011.

Comments: 14 pages, 14 figures, accepted for publication in ApJ

Circumstellar molecular composition of the oxygen-rich AGB star IK~Tau: II. In-depth non-LTE chemical abundance analysis

Authors: L. Decin, E. De Beck, S. Brunken, H. S. P. Muller, K. M. Menten, H. Kim, K. Willacy, A. de Koter, F. Wyrowski

Abstract: Aims: Little information exists on the circumstellar molecular abundance stratifications of many molecules. The aim is to study the circumstellar chemical abundance pattern of 11 molecules and isotopologs ($^{12}$CO, $^{13}$CO, SiS, $^{28}$SiO, $^{29}$SiO, $^{30}$SiO, HCN, CN, CS, SO, SO$_2$) in the oxygen-rich evolved star IK~Tau. Methods: We have performed an in-depth analysis of a large number… ▽ More Aims: Little information exists on the circumstellar molecular abundance stratifications of many molecules. The aim is to study the circumstellar chemical abundance pattern of 11 molecules and isotopologs ($^{12}$CO, $^{13}$CO, SiS, $^{28}$SiO, $^{29}$SiO, $^{30}$SiO, HCN, CN, CS, SO, SO$_2$) in the oxygen-rich evolved star IK~Tau. Methods: We have performed an in-depth analysis of a large number of molecular emission lines excited in the circumstellar envelope around IK~Tau. The analysis is done based on a non-local thermodynamic equilibrium (non-LTE) radiative transfer analysis, which calculates the temperature and velocity structure in a self-consistent way. The chemical abundance pattern is coupled to theoretical outer wind model predictions including photodestruction and cosmic ray ionization. Not only the integrated line intensities, but also the line shapes, are used as diagnostic tool to study the envelope structure. Results: The deduced wind acceleration is much slower than predicted from classical theories. SiO and SiS are depleted in the envelope, possibly due to the adsorption onto dust grains. For HCN and CS a clear difference with respect to inner wind non-equilibrium predictions is found, either indicating uncertainties in the inner wind theoretical modeling or the possibility that HCN and CS (or the radical CN) participate in the dust formation. The low signal-to-noise profiles of SO and CN prohibit an accurate abundance determination; the modeling of high-excitation SO$_2$ lines is cumbersome, possibly related to line misidentifications or problems with the collisional rates. The SiO isotopic ratios ($^{29}$SiO/$^{28}$SiO and $^{30}$SiO/$^{28}$SiO) point toward an enhancement in $^{28}$SiO compared to results of classical stellar evolution codes. Predictions for H$_2$O lines in the spectral range of the Herschel/HIFI mission are performed. [abbreviated] △ Less

Submitted 12 April, 2010; originally announced April 2010.

Comments: 24 pagees, accepted for publication in Astronomy & Astrophysics

Deuterium chemistry in protoplanetary disks II The inner 30 AU

Authors: K. Willacy, P. Woods

Abstract: We present the results of models of the chemistry, including deuterium, in the inner regions of protostellar disks. We find good agreement with recent gas phase observations of several (non--deuterated) species. We also compare our results with observations of comets and find that in the absence of other processing e.g. in the accretion shock at the surface of the disk, or by mixing in the disk,… ▽ More We present the results of models of the chemistry, including deuterium, in the inner regions of protostellar disks. We find good agreement with recent gas phase observations of several (non--deuterated) species. We also compare our results with observations of comets and find that in the absence of other processing e.g. in the accretion shock at the surface of the disk, or by mixing in the disk, the calculated D/H ratios in ices are higher than measured and reflect the D/H ratio set in the molecular cloud phase. Our models give quite different abundances and molecular distributions to other inner disk models because of the differences in physical conditions in the model disk. This emphasizes how changes in the assumptions about the density and temperature distribution can radically affect the results of chemical models. △ Less

Submitted 7 August, 2009; originally announced August 2009.

Comments: Accepted by Astrophysical Journal

Spitzer mapping of molecular hydrogen pure rotational lines in NGC 1333: A detailed study of feedback in star formation

Authors: Sebastien Maret, Edwin A. Bergin, David A. Neufeld, Joel D. Green, Dan M. Watson, Martin O. Harwit, Lars E. Kristensen, Gary J. Melnick, Paule Sonnentrucker, Volker Tolls, Michael W. Werner, Karen Willacy, Yuan Yuan

Abstract: We present mid-infrared spectral maps of the NGC 1333 star forming region, obtained with the the Infrared Spectrometer on board the Spitzer Space Telescope. Eight pure H2 rotational lines, from S (0) to S (7), are detected and mapped. The H2 emission appears to be associated with the warm gas shocked by the multiple outflows present in the region. A comparison between the observed intensities an… ▽ More We present mid-infrared spectral maps of the NGC 1333 star forming region, obtained with the the Infrared Spectrometer on board the Spitzer Space Telescope. Eight pure H2 rotational lines, from S (0) to S (7), are detected and mapped. The H2 emission appears to be associated with the warm gas shocked by the multiple outflows present in the region. A comparison between the observed intensities and the predictions of detailed shock models indicates that the emission arises in both slow (12 - 24 km/s) and fast (36 - 53 km/s) C-type shocks with an initial ortho-to-para ratio of ~ 1. The present H2 ortho-to-para ratio exhibits a large degree of spatial variations. In the post-shocked gas, it is usually about 2, i.e. close to the equilibrium value (~ 3). However, around at least two outflows, we observe a region with a much lower (~ 0.5) ortho-to-para ratio. This region probably corresponds to gas which has been heated-up recently by the passage of a shock front, but whose ortho-to-para has not reached equilibrium yet. This, together with the low initial ortho-to-para ratio needed to reproduce the observed emission, provide strong evidence that H2 is mostly in para form in cold molecular clouds. The H2 lines are found to contribute to 25 - 50% of the total outflow luminosity, and thus can be used to ascertain the importance of star formation feedback on the natal cloud. From these lines, we determine the outflow mass loss rate and, indirectly, the stellar infall rate, the outflow momentum and the kinetic energy injected into the cloud over the embedded phase. The latter is found to exceed the binding energy of individual cores, suggesting that outflows could be the main mechanism for core disruption. △ Less

Submitted 3 April, 2009; originally announced April 2009.

Comments: Accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J.698:1244-1260,2009

Carbon isotope fractionation in protoplanetary disks

Authors: Paul M. Woods, Karen Willacy

Abstract: We investigate the gas-phase and grain-surface chemistry in the inner 30 AU of a typical protoplanetary disk using a new model which calculates the gas temperature by solving the gas heating and cooling balance and which has an improved treatment of the UV radiation field. We discuss inner-disk chemistry in general, obtaining excellent agreement with recent observations which have probed the mat… ▽ More We investigate the gas-phase and grain-surface chemistry in the inner 30 AU of a typical protoplanetary disk using a new model which calculates the gas temperature by solving the gas heating and cooling balance and which has an improved treatment of the UV radiation field. We discuss inner-disk chemistry in general, obtaining excellent agreement with recent observations which have probed the material in the inner regions of protoplanetary disks. We also apply our model to study the isotopic fractionation of carbon. Results show that the fractionation ratio, 12C/13C, of the system varies with radius and height in the disk. Different behaviour is seen in the fractionation of different species. We compare our results with 12C/13C ratios in the Solar System comets, and find a stark contrast, indicative of reprocessing. △ Less

Submitted 1 December, 2008; originally announced December 2008.

Comments: 50 pages, 10 figures, accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J.693:1360-1378,2009

Saturn Forms by Core Accretion in 3.4 Myr

Authors: Sarah E. Dodson-Robinson, Peter Bodenheimer, Gregory Laughlin, Karen Willacy, Neal J. Turner, C. A. Beichman

Abstract: We present two new in situ core accretion simulations of Saturn with planet formation timescales of 3.37 Myr (model S0) and 3.48 Myr (model S1), consistent with observed protostellar disk lifetimes. In model S0, we assume rapid grain settling reduces opacity due to grains from full interstellar values (Podolak 2003). In model S1, we do not invoke grain settling, instead assigning full interstell… ▽ More We present two new in situ core accretion simulations of Saturn with planet formation timescales of 3.37 Myr (model S0) and 3.48 Myr (model S1), consistent with observed protostellar disk lifetimes. In model S0, we assume rapid grain settling reduces opacity due to grains from full interstellar values (Podolak 2003). In model S1, we do not invoke grain settling, instead assigning full interstellar opacities to grains in the envelope. Surprisingly, the two models produce nearly identical formation timescales and core/atmosphere mass ratios. We therefore observe a new manifestation of core accretion theory: at large heliocentric distances, the solid core growth rate (limited by Keplerian orbital velocity) controls the planet formation timescale. We argue that this paradigm should apply to Uranus and Neptune as well. △ Less

Submitted 1 October, 2008; originally announced October 2008.

Comments: 4 pages, including 1 figure, submitted to ApJ Letters

Ice Lines, Planetesimal Composition and Solid Surface Density in the Solar Nebula

Authors: Sarah E. Dodson-Robinson, Karen Willacy, Peter Bodenheimer, Neal J. Turner, C. A. Beichman

Abstract: To date, there is no core accretion simulation that can successfully account for the formation of Uranus or Neptune within the observed 2-3 Myr lifetimes of protoplanetary disks. Since solid accretion rate is directly proportional to the available planetesimal surface density, one way to speed up planet formation is to take a full accounting of all the planetesimal-forming solids present in the… ▽ More To date, there is no core accretion simulation that can successfully account for the formation of Uranus or Neptune within the observed 2-3 Myr lifetimes of protoplanetary disks. Since solid accretion rate is directly proportional to the available planetesimal surface density, one way to speed up planet formation is to take a full accounting of all the planetesimal-forming solids present in the solar nebula. By combining a viscously evolving protostellar disk with a kinetic model of ice formation, we calculate the solid surface density in the solar nebula as a function of heliocentric distance and time. We find three effects that strongly favor giant planet formation: (1) a decretion flow that brings mass from the inner solar nebula to the giant planet-forming region, (2) recent lab results (Collings et al. 2004) showing that the ammonia and water ice lines should coincide, and (3) the presence of a substantial amount of methane ice in the trans-Saturnian region. Our results show higher solid surface densities than assumed in the core accretion models of Pollack et al. (1996) by a factor of 3 to 4 throughout the trans-Saturnian region. We also discuss the location of ice lines and their movement through the solar nebula, and provide new constraints on the possible initial disk configurations from gravitational stability arguments. △ Less

Submitted 4 December, 2008; v1 submitted 23 June, 2008; originally announced June 2008.

Comments: Version 2: reflects lead author's name and affiliation change, contains minor changes to text from version 1. 12 figures, 7 tables, accepted for publication in Icarus

The chemistry of multiply deuterated molecules in protoplanetary disks. I. The outer disk

Authors: K. Willacy

Abstract: We present new models of the deuterium chemistry in protoplanetary disks, including, for the first time, multiply deuterated species. We use these models to explore whether observations in combination with models can give us clues as to which desorption processes occur in disks. We find, in common with other authors, that photodesorption can allow strongly bound molecules such as HDO to exist in… ▽ More We present new models of the deuterium chemistry in protoplanetary disks, including, for the first time, multiply deuterated species. We use these models to explore whether observations in combination with models can give us clues as to which desorption processes occur in disks. We find, in common with other authors, that photodesorption can allow strongly bound molecules such as HDO to exist in the gas phase in a layer above the midplane. Models including this process give the best agreement with the observations. In the midplane, cosmic ray heating can desorb weakly bound molecules such as CO and N$_2$. We find the observations suggest that N$_2$ is gaseous in this region, but that CO must be retained on the grains to account for the observed DCO$^+$/HCO$^+$. This could be achieved by CO having a higher binding energy than N$_2$ (as may be the case when these molecules are accreted onto water ice) or by a smaller cosmic ray desorption rate for CO than assumed here, as suggested by recent theoretical work. For gaseous molecules the calculated deuteration can be greatly changed by chemical processing in the disk from the input molecular cloud values. On the grains singly deuterated species tend to retain the D/H ratio set in the molecular cloud, whereas multiply deuterated species are more affected by the disk chemistry. Consequently the D/H ratios observed in comets may be partly set in the parent cloud and partly in the disk, depending on the molecule. △ Less

Submitted 16 January, 2007; originally announced January 2007.

Comments: Accepted for publication in ApJ. 48 pages, 8 figures

Journal ref: Astrophys.J.660:441-460,2007

Benzene formation in the inner regions of protostellar disks

Authors: Paul M. Woods, Karen Willacy

Abstract: Benzene (c-C6H6) formation in the inner 3 AU of a protostellar disk can be efficient, resulting in high abundances of benzene in the midplane region. The formation mechanism is different to that found in interstellar clouds and in protoplanetary nebulae, and proceeds mainly through the reaction between allene (C3H4) and its ion. This has implications for PAH formation, in that some fraction of P… ▽ More Benzene (c-C6H6) formation in the inner 3 AU of a protostellar disk can be efficient, resulting in high abundances of benzene in the midplane region. The formation mechanism is different to that found in interstellar clouds and in protoplanetary nebulae, and proceeds mainly through the reaction between allene (C3H4) and its ion. This has implications for PAH formation, in that some fraction of PAHs seen in the solar system could be native rather than inherited from the interstellar medium. △ Less

Submitted 8 December, 2006; originally announced December 2006.

Comments: 9 pages, 2 colour figures, to be published in the Astrophysical Journal Letters

Journal ref: Astrophys.J.655:L49-L52,2007

Turbulence driven diffusion in protoplanetary disks - chemical effects in the outer disk

Authors: K. Willacy, W. D. Langer, M. Allen, G. Bryden

Abstract: The dynamics and chemistry of protostellar disks are likely to be intricately linked, with dynamical processes altering the chemical composition, and chemistry, in turn, controlling the ionization structure and hence the ability of the magneto-rotational instability to drive the disk turbulence. Here we present the results from the first chemical models of the outer regions (R > 100 AU) of proto… ▽ More The dynamics and chemistry of protostellar disks are likely to be intricately linked, with dynamical processes altering the chemical composition, and chemistry, in turn, controlling the ionization structure and hence the ability of the magneto-rotational instability to drive the disk turbulence. Here we present the results from the first chemical models of the outer regions (R > 100 AU) of protoplanetary disks to consider the effects of turbulence driven diffusive mixing in the vertical direction. We show that vertical diffusion can greatly affect the column densities of many species, increasing them by factors of up to two orders of magnitude. Previous disk models have shown that disks can be divided into three chemically distinct layers, with the bulk of the observed molecular emission coming from a region between an atomic/ionic layer on the surface of the disk and the midplane regoin where the bulk of molecules are frozen onto grains. Diffusion retains this three layer structure, but increases the depth of the molecular layer by bringing atoms and atomic ions form by photodissociation in the surface layers into the shielded molecular layer where molecules can reform. For other species, notably NH3 and N2H+, the column densities are relatively unaffected by diffusion. These species peak in abundance near the midplane where most other molecules are heavily depleted, rather than in the molecular layer above. Diffusion only affects the abundances of those molecules with peak abundances in the molecular layer. We find that diffusion does not affect the ionization fraction of the disk. We compare the calculated column densities to observations of DM Tau, LkCa 15 and TW Hya and find good agreement for many molecules with a diffusion coefficient of 1e18 cm^2 s^-1. △ Less

Submitted 3 March, 2006; originally announced March 2006.

Comments: Accepted in ApJ

Journal ref: Astrophys.J. 644 (2006) 1202-1213

Turbulent Mixing in the Outer Solar Nebula

Authors: N. J. Turner, K. Willacy, G. Bryden, H. W. Yorke

Abstract: The effects of turbulence on the mixing of gases and dust in the outer Solar nebula are examined using 3-D MHD calculations in the shearing-box approximation with vertical stratification. The turbulence is driven by the magneto-rotational instability. The magnetic and hydrodynamic stresses in the turbulence correspond to an accretion time at the midplane about equal to the lifetimes of T Tauri d… ▽ More The effects of turbulence on the mixing of gases and dust in the outer Solar nebula are examined using 3-D MHD calculations in the shearing-box approximation with vertical stratification. The turbulence is driven by the magneto-rotational instability. The magnetic and hydrodynamic stresses in the turbulence correspond to an accretion time at the midplane about equal to the lifetimes of T Tauri disks, while accretion in the surface layers is thirty times faster. The mixing resulting from the turbulence is also fastest in the surface layers. The mixing rate is similar to the rate of radial exchange of orbital angular momentum, so that the Schmidt number is near unity. The vertical spreading of a trace species is well-matched by solutions of a damped wave equation when the flow is horizontally-averaged. The damped wave description can be used to inexpensively treat mixing in 1-D chemical models. However, even in calculations reaching a statistical steady state, the concentration at any given time varies substantially over horizontal planes, due to fluctuations in the rate and direction of the transport. In addition to mixing species that are formed under widely varying conditions, the turbulence intermittently forces the nebula away from local chemical equilibrium. The different transport rates in the surface layers and interior may affect estimates of the grain evolution and molecular abundances during the formation of the Solar system. △ Less

Submitted 10 November, 2005; originally announced November 2005.

Comments: To appear in the Astrophysical Journal; 20 pages, 9 figures

Journal ref: Astrophys.J.639:1218-1226,2006

Electron Abundance in Protostellar Cores

Authors: Paolo Padoan, Karen Willacy, William Langer, Mika Juvela

Abstract: The determination of the fractional electron abundance, Xe, in protostellar cores relies on observations of molecules, such as DCO+, H13CO+ and CO, and on chemical models to interpret their abundance. Studies of protostellar cores have revealed significant variations of Xe from core to core within a range 10^-8<Xe<10^-6. The physical origin of these large variations in Xe is not well understood,… ▽ More The determination of the fractional electron abundance, Xe, in protostellar cores relies on observations of molecules, such as DCO+, H13CO+ and CO, and on chemical models to interpret their abundance. Studies of protostellar cores have revealed significant variations of Xe from core to core within a range 10^-8<Xe<10^-6. The physical origin of these large variations in Xe is not well understood, unless unlikely variations in the cosmic ray ionization rate or ad hoc values of metal depletion are assumed. In this work we explore other potential causes of these variations in Xe, including core age, extinction and density. We compute numerically the intensity of the radiation field within a density distribution generated by supersonic turbulence. Taking into account the lines of sight in all directions, the effective visual extinction in dense regions is found to be always much lower than the extinction derived from the column density along a fixed line of sight. Dense cores with volume and column densities comparable to observed protostellar cores have relatively low mass-averaged visual extinction, 2mag <= A_V <= 5mag, such that photo-ionization can sometimes be as important as cosmic ray ionization. Chemical models, including gas-grain chemistry and time dependent gas depletion and desorption, are computed for values of visual extinction in the range 2mag <= A_V <= 6mag, and for a hydrogen gas density of 10^4cm^-3$, typical of protostellar cores. The models presented here can reproduce some of the observed variations of ion abundance from core to core as the combined effect of visual extinction and age variations. The range of electron abundances predicted by the models is relatively insenstive to density over 10^4 to 10^6 cm^{-3}. △ Less

Submitted 23 June, 2004; originally announced June 2004.

Comments: 25 pages, 9 figures. Astrophysical Journal, in press

Journal ref: Astrophys.J. 614 (2004) 203-210

Carbon molecules in the inner wind of the oxygen-rich Mira IK Tauri

Authors: Debiprosad Duari, Isabelle Cherchneff, Karen Willacy

Abstract: The gas-phase non-equilibrium chemistry of the inner wind of the oxygen-rich Mira variable IK Tau (NML Tau) is investigated using a physio-chemical model describing the periodically shocked gas close to the stellar photosphere. We predict the formation in large amounts of a few carbon species in the inner envelope, in particular CO$_2$, which has been recently detected with ISO in the spectra of… ▽ More The gas-phase non-equilibrium chemistry of the inner wind of the oxygen-rich Mira variable IK Tau (NML Tau) is investigated using a physio-chemical model describing the periodically shocked gas close to the stellar photosphere. We predict the formation in large amounts of a few carbon species in the inner envelope, in particular CO$_2$, which has been recently detected with ISO in the spectra of several oxygen-rich semi-regular and Mira stars. The theoretical abundances are also in excellent agreement with values derived from millimeter and sub-millimeter observations, pointing to the fact that some carbon species in oxygen stars do form from shock chemistry in the inner layers and travel the envelope as ``parent'' species. △ Less

Submitted 17 November, 1998; originally announced November 1998.

Comments: 4 pages (accepted for publication in A&A Letters)

Journal ref: Astron.Astrophys.341:L47-L50,1999