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

Catherine Walsh

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

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.GAphysics.ao-phphysics.geo-ph

Catalog footprint

What is connected

40works
5topics
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

40 published item(s)

preprint2022arXiv

A revised lower estimate of ozone columns during Earth's oxygenated history

The history of molecular oxygen (O$_2$) in Earth's atmosphere is still debated; however, geological evidence supports at least two major episodes where O$_2$ increased by an order of magnitude or more: the Great Oxidation Event (GOE) and the Neoproterozoic Oxidation Event. O$_2$ concentrations have likely fluctuated (between $10^{-3}$ and $1.5$ times the present atmospheric level) since the GOE $\sim 2.4$ Gyr ago, resulting in a time-varying ozone (O$_3$) layer. Using a three-dimensional chemistry-climate model, we simulate changes in O$_3$ in Earth's atmosphere since the GOE and consider the implications for surface habitability, and glaciation during the Mesoproterozoic. We find lower O$_3$ columns (reduced by up to $4.68$ times for a given O$_2$ level) compared to previous work; hence, higher fluxes of biologically harmful UV radiation would have reached the surface. Reduced O$_3$ leads to enhanced tropospheric production of the hydroxyl radical (OH) which then substantially reduces the lifetime of methane (CH$_4$). We show that a CH$_4$ supported greenhouse effect during the Mesoproterozoic is highly unlikely. The reduced O$_3$ columns we simulate have important implications for astrobiological and terrestrial habitability, demonstrating the relevance of three-dimensional chemistry-climate simulations when assessing paleoclimates and the habitability of faraway worlds.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

ALMA High-resolution Multiband Analysis for the Protoplanetary Disk around TW Hya

We present a high-resolution (2.5 au) multiband analysis of the protoplanetary disk around TW Hya using ALMA long baseline data at Bands 3, 4, 6, and 7. We aim to reconstruct a high-sensitivity millimeter continuum image and revisit the spectral index distribution. The imaging is performed by combining new ALMA data at Bands 4 and 6 with available archive data. Two methods are employed to reconstruct the images; multi-frequency synthesis (MFS) and the fiducial image-oriented method, where each band is imaged separately and the frequency dependence is fitted pixel by pixel. We find that the MFS imaging with the second order of Taylor expansion can reproduce the frequency dependence of the continuum emission between Bands 3 and 7 in a manner consistent with previous studies and is a reasonable method to reconstruct the spectral index map. The image-oriented method provides a spectral index map consistent with the MFS imaging, but with a two times lower resolution. Mock observations of an intensity model were conducted to validate the images from the two methods. We find that the MFS imaging provides a high-resolution spectral index distribution with an uncertainty of $<10$~\%. Using the submillimeter spectrum reproduced from our MFS images, we directly calculated the optical depth, power-law index of the dust opacity coefficient ($β$), and dust temperature. The derived parameters are consistent with previous works, and the enhancement of $β$ within the intensity gaps is also confirmed, supporting a deficit of millimeter-sized grains within the gaps.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Astrochemical modelling of infrared dark clouds

Infrared dark clouds (IRDCs) are cold, dense regions of the interstellar medium (ISM) that are likely to represent the initial conditions for massive star formation. It is thus important to study the physical and chemical conditions of IRDCs to provide constraints and inputs for theoretical models of these processes. We aim to determine the astrochemical conditions, especially cosmic ray ionisation rate (CRIR) and chemical age, in different regions of the massive IRDC G28.37+00.07 by comparing observed abundances of multiple molecules and molecular ions with the predictions of astrochemical models. We have computed a series of single-zone astrochemical models with a gas-grain network that systematically explores the parameter space of density, temperature, CRIR, and visual extinction. We have also investigated the effects of choices of CO ice binding energy and temperatures achieved in transient heating of grains when struck by cosmic rays. We selected 10 positions across the IRDC that are known to have a variety of star formation activity. We utilised mid-infrared (MIR) extinction maps and sub-mm emission maps to measure the mass surface densities of these regions, needed for abundance and volume density estimates. The sub-mm emission maps were also used to measure temperatures. We then used IRAM-30m observations of various tracers to estimate column densities and thus abundances. Using estimates of the abundances of CO, HCO$^+$ and N$_2$H$^+$ we find consistency with astrochemical models that have relatively low CRIRs of $ζ\sim10^{-18}$ to $\sim10^{-17}\:{\rm s}^{-1}$, with no evidence for systematic variation with the level of star formation activity. Astrochemical ages are found to be < 1 Myr. We discuss potential sources of systematic uncertainties in these results and the overall implications for IRDC evolutionary history and astrochemical models.(abridged for arXiv)

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Different molecular filament widths as tracers of accretion onto filaments

We explore how dense filament widths, when measured using different molecular species, may change as a consequence of gas accretion toward the filament. As a gas parcel falls into the filament, it will experience different density, temperature, and extinction values. The rate at which this environment changes will affect differently the abundance of different molecules. So, a molecule that forms quickly will better reflect the local physical conditions a gas parcel experiences than a slower-forming molecule. Since these differences depend on how the respective timescales compare, the different molecular distributions should reflect how rapidly the environment changes, i.e., the accretion rate toward the filament. We find that the filament widths measured from time-dependent abundances for C2H, CO, CN, CS, and C3H2, are sensitive the most to this effect, being those molecules the ones presenting also the wider filament widths. On the contrary, molecules such as N2H+, NH3, H2CO, HNC or CH3OH are not so sensitive to accretion and present the narrowest filament widths. We propose that ratios of filament widths for different tracers could be a useful tool to estimate the accretion rate onto the filament.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fevering Interstellar Ices Have More CH3OD

Mono-deuterated methanol is thought to form during the prestellar core stage of star formation. Observed variations in the CH2DOH/CH3OD ratio suggest that its formation is strongly dependent on the surrounding cloud conditions. Thus, it is a potential tracer of the physical conditions before the onset of star formation. A single-point physical model representative of a typical prestellar core is coupled to chemical models to investigate potential formation pathways towards deuterated methanol at the prestellar stage. Simple addition reactions of H and D are not able to reproduce observed abundances. The implementation of an experimentally verified abstraction scheme leads to the efficient formation of methyl-deuterated methanol, but lacks sufficient formation of hydroxy-deuterated methanol. CH3OD is most likely formed at a later evolutionarymstage, potentially from H-D exchange reactions in warm ices between HDO (and D2O) and CH3OH. The CH2DOH/CH3OD ratio is not an appropriate tracer of the physical conditions during the prestellar stage, but might be better suited as a tracer of ice heating.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Modeling snowline locations in protostars: The impact of the structure of protostellar cloud cores

Abridged Context: Snowlines during star and disk formation are responsible for a range of effects during the evolution of protostars, such as setting the chemical composition of the envelope and disk. This in turn influences the formation of planets by changing the elemental compositions of solids and affecting the collisional properties and outcomes of dust grains. Snowlines can also reveal accretion bursts, providing insight into the formation process of stars. Methods: A numerical chemical network coupled with a grid of cylindrical-symmetric physical models was used to identify what physical parameters alter the CO and H$_2$O snowline locations. The investigated parameters are the initial molecular abundances, binding energies of CO and H$_2$O, heating source, cloud core density, outflow cavity opening angle, and disk geometry. Simulated molecular line emission maps were used to quantify the change in the snowline location with each parameter. Conclusions: The models presented in this work show that the CO and H$_2$O snowline locations do not occur at a single, well-defined temperature as is commonly assumed. Instead, the snowline position depends on luminosity, cloud core density, and whether a disk is present or not. Inclination and spatial resolution affect the observability and successful measurement of snowline locations. We note that N$_2$H$^+$ and HCO$^+$ emission serve as good observational tracers of CO and H$_2$O snowline locations. However, constraints on whether or not a disk is present, the observation of additional molecular tracers, and estimating envelope density will help in accurately determining the cause of the observed snowline position. Plots of the N$_2$H$^+$ and HCO$^+$ peak emission radius versus luminosity are provided to compare the models with observations of deeply embedded protostars aiming to measure the CO and H$_2$O snowline locations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Molecules with ALMA at Planet-forming Scales (MAPS) III: Characteristics of Radial Chemical Substructures

The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a detailed, high resolution (${\sim}$10-20 au) view of molecular line emission in five protoplanetary disks at spatial scales relevant for planet formation. Here, we present a systematic analysis of chemical substructures in 18 molecular lines toward the MAPS sources: IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. We identify more than 200 chemical substructures, which are found at nearly all radii where line emission is detected. A wide diversity of radial morphologies - including rings, gaps, and plateaus - is observed both within each disk and across the MAPS sample. This diversity in line emission profiles is also present in the innermost 50 au. Overall, this suggests that planets form in varied chemical environments both across disks and at different radii within the same disk. Interior to 150 au, the majority of chemical substructures across the MAPS disks are spatially coincident with substructures in the millimeter continuum, indicative of physical and chemical links between the disk midplane and warm, elevated molecular emission layers. Some chemical substructures in the inner disk and most chemical substructures exterior to 150 au cannot be directly linked to dust substructure, however, which indicates that there are also other causes of chemical substructures, such as snowlines, gradients in UV photon fluxes, ionization, and radially-varying elemental ratios. This implies that chemical substructures could be developed into powerful probes of different disk characteristics, in addition to influencing the environments within which planets assemble. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Molecules with ALMA at Planet-forming Scales (MAPS). A Circumplanetary Disk Candidate in Molecular Line Emission in the AS 209 Disk

We report the discovery of a circumplanetary disk (CPD) candidate embedded in the circumstellar disk of the T Tauri star AS 209 at a radial distance of about 200 au (on-sky separation of 1."4 from the star at a position angle of $161^\circ$), isolated via $^{13}$CO $J=2-1$ emission. This is the first instance of CPD detection via gaseous emission capable of tracing the overall CPD mass. The CPD is spatially unresolved with a $117\times82$ mas beam and manifests as a point source in $^{13}$CO, indicating that its diameter is $\lesssim14$ au. The CPD is embedded within an annular gap in the circumstellar disk previously identified using $^{12}$CO and near-infrared scattered light observations, and is associated with localized velocity perturbations in $^{12}$CO. The coincidence of these features suggests that they have a common origin: an embedded giant planet. We use the $^{13}$CO intensity to constrain the CPD gas temperature and mass. We find that the CPD temperature is $\gtrsim35$ K, higher than the circumstellar disk temperature at the radial location of the CPD, 22 K, suggesting that heating sources localized to the CPD must be present. The CPD gas mass is $\gtrsim 0.095 M_{\rm Jup} \simeq 30 M_{\rm Earth}$ adopting a standard $^{13}$CO abundance. From the non-detection of millimeter continuum emission at the location of the CPD ($3σ$ flux density $\lesssim26.4~μ$Jy), we infer that the CPD dust mass is $\lesssim 0.027 M_{\rm Earth} \simeq 2.2$ lunar masses, indicating a low dust-to-gas mass ratio of $\lesssim9\times10^{-4}$. We discuss the formation mechanism of the CPD-hosting giant planet on a wide orbit in the framework of gravitational instability and pebble accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Sulphur monoxide emission tracing an embedded planet in the HD 100546 protoplanetary disk

Molecular line observations are powerful tracers of the physical and chemical conditions across the different evolutionary stages of star, disk and planet formation. Using the high angular resolution and unprecedented sensitivity of the Atacama Large Millimeter Array (ALMA) there is now a drive to detect small-scale gas structures in protoplanetary disks that can be attributed directly to forming planets. We report high angular resolution ALMA Band 7 observations of sulphur monoxide (SO) in the nearby planet-hosting disk around Herbig star HD 100546. SO is rarely detected in evolved protoplanetary disks but in other environments, it is most often utilised as a tracer of shocks. The SO emission from the HD 100546 disk is primarily originating from gas within the approx. 20 au mm-dust cavity and shows a clear azimuthal brightness asymmetry of a factor of 2. In addition, we see a significant difference in the line profile shape when comparing these new Cycle 7 data to Cycle 0 data of the same SO transitions. We discuss the different physical/chemical mechanisms that could be responsible for this asymmetry and time variability including disk winds, disk warps, and a shock triggered by a (forming) planet. We propose that the SO is enhanced in the cavity due to the presence of a giant planet. The SO asymmetry complements evidence for hot circumplanetary material around the giant planet HD 100546 c traced via CO ro-vibrational emission. This work sets the stage for further observational and modelling efforts to detect and understand the chemical imprint of a forming planet on its parent disk.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Unveiling the outer dust disc of TW Hya with deep ALMA observations

The radial extent of millimetre dust in protoplanetary discs is often far smaller than that of their gas, mostly due to processes such as dust growth and radial drift. However, it has been suggested that current millimetre continuum observations of discs do not trace their full extent due to limited sensitivity. In this Letter, we present deep (19 $μ$Jy beam$^{-1}$) moderate resolution (0.37") ALMA observations at 1 mm of the nearest protoplanetary disc, TW Hya. Using the visibility analysis tool `frank', we reveal a structured millimetre intensity distribution out to 100 au, well beyond previous estimates of 60-70 au. Our analysis suggests the presence of a new millimetre continuum gap at 82 au, coincident with similar features seen in optical/near-infrared scattered light and millimetre molecular line observations. Examination of the fit residuals confirms the presence of the previously reported au-scale continuum excess at 52 au (P.A. = 242.5 degrees). Our results demonstrate the utility of combining deep, moderate resolution observations with super-resolution analysis techniques to probe the faintest regions of protoplanetary discs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An unbiased ALMA spectral survey of the LkCa 15 and MWC 480 protoplanetary disks

The volatile contents of protoplanetary disks both set the potential for planetary chemistry and provide valuable probes of defining disk system characteristics such as stellar mass, gas mass, ionization, and temperature structure. Current disk molecular inventories are fragmented, however, giving an incomplete picture: unbiased spectral line surveys are needed to assess the volatile content. We present here an overview of such a survey of the protoplanetary disks around the Herbig Ae star MWC 480 and the T Tauri star LkCa 15 in ALMA Band 7, spanning $\sim$36 GHz from 275--317 GHz and representing an order of magnitude increase in sensitivity over previous single-dish surveys. We detect 14 molecular species (including isotopologues), with 5 species (C$^{34}$S, $^{13}$CS, H$_{2}$CS, DNC, and C$_2$D) detected for the first time in protoplanetary disks. Significant differences are observed in the molecular inventories of MWC~480 and LkCa~15, and we discuss how these results may be interpreted in light of the different physical conditions of these two disk systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Temperature structures of embedded disks: young disks in Taurus are warm

The chemical composition of gas and ice in disks around young stars set the bulk composition of planets. In contrast to protoplanetary disks (Class II), young disks that are still embedded in their natal envelope (Class 0 and I) are predicted to be too warm for CO to freeze out, as has been confirmed observationally for L1527 IRS. To establish whether young disks are generally warmer than their more evolved counterparts, we observed five young (Class 0/I and Class I) disks in Taurus with the Atacama Large Millimeter/submillimeter Array (ALMA), targeting C$^{17}$O $2-1$, H$_2$CO $3_{1,2}-2_{1,1}$, HDO $3_{1,2}-2_{2,1}$ and CH$_3$OH $5_K-4_K$ transitions at $0.48^{\prime\prime} \times 0.31^{\prime\prime}$ resolution. The different freeze-out temperatures of these species allow us to derive a global temperature structure. C$^{17}$O and H$_2$CO are detected in all disks, with no signs of CO freeze-out in the inner $\sim$100 au, and a CO abundance close to $\sim$10$^{-4}$. H$_2$CO emission originates in the surface layers of the two edge-on disks, as witnessed by the especially beautiful V-shaped emission pattern in IRAS~04302+2247. HDO and CH$_3$OH are not detected, with column density upper limits more than 100 times lower than for hot cores. Young disks are thus found to be warmer than more evolved protoplanetary disks around solar analogues, with no CO freeze-out (or only in the outermost part of $\gtrsim$100 au disks) or CO processing. However, they are not as warm as hot cores or disks around outbursting sources, and therefore do not have a large gas-phase reservoir of complex molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The TW Hya Rosetta Stone Project I: Radial and vertical distributions of DCN and DCO+

Molecular D/H ratios are frequently used to probe the chemical past of Solar System volatiles. Yet it is unclear which parts of the Solar Nebula hosted an active deuterium fractionation chemistry. To address this question, we present 0".2-0".4 ALMA observations of DCO+ and DCN 2-1, 3-2 and 4-3 towards the nearby protoplanetary disk around TW Hya, taken as part of the TW Hya Rosetta Stone project, augmented with archival data. DCO+ is characterized by an excitation temperature of ~40 K across the 70 au radius pebble disk, indicative of emission from a warm, elevated molecular layer. Tentatively, DCN is present at even higher temperatures. Both DCO+ and DCN present substantial emission cavities in the inner disk, while in the outer disk the DCO+ and DCN morphologies diverge: most DCN emission originates from a narrow ring peaking around 30~au, with some additional diffuse DCN emission present at larger radii, while DCO+ is present in a broad structured ring that extends past the pebble disk. Based on parametric disk abundance models, these emission patterns can be explained by a near-constant DCN abundance exterior to the cavity, and an increasing DCO+ abundance with radius. There appears to be an active deuterium fractionation chemistry in multiple disk regions around TW Hya, but not in the cold planetesimal-forming midplane and in the inner disk. More observations are needed to explore whether deuterium fractionation is actually absent in these latter regions, and if its absence is a common feature, or something peculiar to the old TW Hya disk.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The TW Hya Rosetta Stone Project II: Spatially resolved emission of formaldehyde hints at low-temperature gas-phase formation

Formaldehyde (H$_2$CO) is an important precursor to organics like methanol (CH$_3$OH). It is important to understand the conditions that produce H$_2$CO and prebiotic molecules during star and planet formation. H$_2$CO possesses both gas-phase and solid-state formation pathways, involving either UV-produced radical precursors or CO ice and cold ($\lesssim 20$ K) dust grains. To understand which pathway dominates, gaseous H$_2$CO's ortho-to-para ratio (OPR) has been used as a probe, with a value of 3 indicating "warm" conditions and $<3$ linked to cold formation in the solid-state. We present spatially resolved ALMA observations of multiple ortho- and para-H$_2$CO transitions in the TW Hya protoplanetary disk to test H$_2$CO formation theories during planet formation. We find disk-averaged rotational temperatures and column densities of $33\pm2$ K, ($1.1\pm0.1)\times10^{12}$ cm$^{-2}$ and $25\pm2$ K, $(4.4\pm0.3)\times10^{11}$ cm$^{-2}$ for ortho- and para-H$_2$CO, respectively, and an OPR of $2.49\pm0.23$. A radially resolved analysis shows that the observed H$_2$CO emits mostly at rotational temperatures of 30-40 K, corresponding to a layer with $z/R\ge0.25$. The OPR is consistent with 3 within 60 au, the extent of the pebble disk, and decreases beyond 60 au to $2.0\pm0.5$. The latter corresponds to a spin temperature of 12 K, well below the rotational temperature. The combination of relatively uniform emitting conditions, a radial gradient in the OPR, and recent laboratory experiments and theory on OPR ratios after sublimation, lead us to speculate that gas-phase formation is responsible for the observed H$_2$CO across the TW Hya disk.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The TW Hya Rosetta Stone Project III: Resolving the Gaseous Thermal Profile of the Disk

The thermal structure of protoplanetary disks is a fundamental characteristic of the system that has wide reaching effects on disk evolution and planet formation. In this study, we constrain the 2D thermal structure of the protoplanetary disk TW Hya structure utilizing images of seven CO lines. This includes new ALMA observations of 12CO J=2-1 and C18O J=2-1 as well as archival ALMA observations of 12CO J=3-2, 13CO J=3-2, 6-5, C18O J= 3-2, 6-5. Additionally, we reproduce a Herschel observation of the HD J=1-0 line flux, the spectral energy distribution, and utilize a recent quantification of CO radial depletion in TW Hya. These observations were modeled using the thermochemical code RAC2D, and our best fit model reproduces all spatially resolved CO surface brightness profiles. The resulting thermal profile finds a disk mass of 0.025 Msun and a thin upper layer of gas depleted of small dust with a thickness of approx 1.2% of the corresponding radius. Using our final thermal structure, we find that CO alone is not a viable mass tracer as its abundance is degenerate with the total H2 surface density. Different mass models can readily match the spatially resolved CO line profiles with disparate abundance assumptions. Mass determination requires additional knowledge and, in this work, HD provides the additional constraint to derive the gas mass and supports the inference of CO depletion in the TW Hya disk. Our final thermal structure confirms the use of HD as a powerful probe of protoplanetary disk mass. Additionally, the method laid out in this paper is an employable strategy for extraction of disk temperatures and masses in the future.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A primordial origin for molecular oxygen in comets: A chemical kinetics study of the formation and survival of O$_2$ ice from clouds to disks

Molecular oxygen has been confirmed as the fourth most abundant molecule in cometary material O$_2$/H$_2$O $\sim 4$ %) and is thought to have a primordial nature, i.e., coming from the interstellar cloud from which our solar system was formed. However, interstellar O$_2$ gas is notoriously difficult to detect and has only been observed in one potential precursor of a solar-like system. Here, the chemical and physical origin of O$_2$ in comets is investigated using sophisticated astrochemical models. Three origins are considered: i) in dark clouds, ii) during forming protostellar disks, and iii) during luminosity outbursts in disks. The dark cloud models show that reproduction of the observed abundance of O$_2$ and related species in comet 67P/C-G requires a low H/O ratio facilitated by a high total density ($\geq 10^5$ cm$^{-3}$), and a moderate cosmic ray ionisation rate ($\leq 10^{-16}$ s$^{-1}$) while a temperature of 20 K, slightly higher than the typical temperatures found in dark clouds, also enhances the production of O$_2$. Disk models show that O$_2$ can only be formed in the gas phase in intermediate disk layers, and cannot explain the strong correlation between O$_2$ and H$_2$O in comet 67P/C-G together with the weak correlation between other volatiles and H$_2$O. However, primordial O$_2$ ice can survive transport into the comet-forming regions of disks. Taken together, these models favour a dark cloud (or "primordial") origin for O$_2$ in comets, albeit for dark clouds which are warmer and denser than those usually considered as solar system progenitors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

ALMA Observations of a Gap and a Ring in the Protoplanetary Disk around TW Hya

We report the first detection of a gap and a ring in 336 GHz dust continuum emission from the protoplanetary disk around TW Hya, using the Atacama Large Millimeter/Submillimeter Array (ALMA). The gap and ring are located at around 25 and 41 au from the central star, respectively, and are associated with the CO snow line at ~30 au. The gap has a radial width of less than 15 au and a mass deficit of more than 23%, taking into account that the observations are limited to an angular resolution of ~15 au. In addition, the 13CO and C18O J = 3 - 2 lines show a decrement in CO line emission throughout the disk, down to ~10 au, indicating a freeze-out of gas-phase CO onto grain surfaces and possible subsequent surface reactions to form larger molecules. The observed gap could be caused by gravitational interaction between the disk gas and a planet with a mass less than super-Neptune (2M_{Neptune}), or could be the result of the destruction of large dust aggregates due to the sintering of CO ice.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

ALMA reveals the anatomy of the mm-sized dust and molecular gas in the HD 97048 disk

Transitional disks show a lack of excess emission at infrared wavelengths due to a large dust cavity, that is often corroborated by spatially resolved observations at ~ mm wavelengths. We present the first spatially resolved ~ mm-wavelength images of the disk around the Herbig Ae/Be star, HD 97048. Scattered light images show that the disk extends to ~640 au. The ALMA data reveal a circular-symmetric dusty disk extending to ~350 au, and a molecular disk traced in CO J=3-2 emission, extending to ~750 au. The CO emission arises from a flared layer with an opening angle ~ 30 deg - 40 deg. HD 97048 is another source for which the large (~ mm-sized) dust grains are more centrally concentrated than the small (~ μm-sized) grains and molecular gas, likely due to radial drift. The images and visibility data modelling suggests a decrement in continuum emission within ~50 au, consistent with the cavity size determined from mid-infrared imaging (34 +/- 4 au). The extracted continuum intensity profiles show ring-like structures with peaks at ~50, 150, and 300 au, with associated gaps at ~100 and 250 au. This structure should be confirmed in higher-resolution images (FWHM ~ 10 - 20 au). These data confirm the classification of HD 97048 as a transitional disk that also possesses multiple ring-like structures in the dust continuum emission. Additional data are required at multiple and well-separated frequencies to fully characterise the disk structure, and thereby constrain the mechanism(s) responsible for sculpting the HD 97048 disk.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Candidate Water Vapor Lines to Locate the H$_{2}$O Snowline through High-Dispersion Spectroscopic Observations I. The Case of a T Tauri Star

Inside the H$_{2}$O snowline of protoplanetary disks, water evaporates from the dust-grain surface into the gas phase, whereas it is frozen out on to the dust in the cold region beyond the snowline. H$_{2}$O ice enhances the solid material in the cold outer part of a disk, which promotes the formation of gas-giant planet cores. We can regard the H$_{2}$O snowline as the surface that divides the regions between rocky and gaseous giant planet formation. Thus observationally measuring the location of the H$_{2}$O snowline is crucial for understanding the planetesimal and planet formation processes, and the origin of water on Earth. In this paper, we find candidate water lines to locate the H$_{2}$O snowline through future high-dispersion spectroscopic observations. First, we calculate the chemical composition of the disk and investigate the abundance distributions of H$_{2}$O gas and ice, and the position of the H$_{2}$O snowline. We confirm that the abundance of H$_{2}$O gas is high not only in the hot midplane region inside the H$_{2}$O snowline but also in the hot surface layer of the outer disk. Second, we calculate the H$_{2}$O line profiles and identify those H$_{2}$O lines which are promising for locating the H$_{2}$O snowline: the identified lines are those which have small Einstein $A$ coefficients and high upper state energies. The wavelengths of the candidate H$_{2}$O lines range from mid-infrared to sub-millimeter, and they overlap with the regions accessible to ALMA and future mid-infrared high dispersion spectrographs (e.g., TMT/MICHI, SPICA).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Chemical complexity in protoplanetary disks in the era of ALMA and Rosetta

Comets provide a unique insight into the molecular composition and complexity of the material in the primordial solar nebula. Recent results from the Rosetta mission, currently monitoring comet 67P/Churyumov-Gerasimenko in situ, and ALMA (the Atacama Large Millimeter/submillimeter Array), have demonstrated a tantalising link between the chemical complexity now confirmed in disks (via the detection of gas-phase CH3CN; Oberg et al. 2015) and that confirmed on the surface of 67P (Goesmann et al. 2015), raising questions concerning the chemical origin of such species (cloud or inheritance versus disk synthesis). Results from an astrochemical model of a protoplanetary disk are presented in which complex chemistry is included and in which it is assumed that simple ices only are inherited from the parent molecular cloud. The model results show good agreement with the abundances of several COMs observed on the surface of 67P with Philae/COSAC. Cosmic-ray and X-ray-induced photoprocessing of predominantly simple ices inherited by the protoplanetary disk is sufficient to generate a chemical complexity similar to that observed in comets. This indicates that the icy COMs detected on the surface of 67P may have a disk origin. The results also show that gas-phase CH3CN is abundant in the inner warm disk atmosphere where hot gas-phase chemistry dominates and potentially erases the ice chemical signature. Hence, CH3CN may not be an unambiguous tracer of the complex organic ice reservoir. However, a better understanding of the hot gas-phase chemistry of CH3CN is needed to confirm this preliminary conclusion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cometary ices in forming protoplanetary disc midplanes

Low-mass protostars are the extrasolar analogues of the natal Solar System. Sophisticated physicochemical models are used to simulate the formation of two protoplanetary discs from the initial prestellar phase, one dominated by viscous spreading and the other by pure infall. The results show that the volatile prestellar fingerprint is modified by the chemistry en route into the disc. This holds relatively independent of initial abundances and chemical parameters: physical conditions are more important. The amount of CO2 increases via the grain-surface reaction of OH with CO, which is enhanced by photodissociation of H2O ice. Complex organic molecules are produced during transport through the envelope at the expense of CH3OH ice. Their abundances can be comparable to that of methanol ice (few % of water ice) at large disc radii (R > 30 AU). Current Class II disc models may be underestimating the complex organic content. Planet population synthesis models may underestimate the amount of CO2 and overestimate CH3OH ices in planetesimals by disregarding chemical processing between the cloud and disc phases. The overall C/O and C/N ratios differ between the gas and solid phases. The two ice ratios show little variation beyond the inner 10 AU and both are nearly solar in the case of pure infall, but both are sub-solar when viscous spreading dominates. Chemistry in the protostellar envelope en route to the protoplanetary disc sets the initial volatile and prebiotically-significant content of icy planetesimals and cometary bodies. Comets are thus potentially reflecting the provenances of the midplane ices in the Solar Nebula.

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

Grand challenges in protoplanetary disc modelling

The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Setting the volatile composition of (exo)planet-building material. Does chemical evolution in disk midplanes matter?

[Abridged] Chemical evolution in the protoplanetary disk midplane can modify the composition of ices and gases. We have investigated if and how chemical evolution affects the abundances and distributions of key volatile species in the midplane of a protoplanetary disk in the 0.2-30 AU range. A full chemical network including gas-phase, gas-grain interactions and grain-surface chemistry is used to evolve chemistry in time, for 1 Myr. Great diversity is observed in the relative abundance ratios of the main considered species: H2O, CO, CO2, CH4, O2, NH3 and N2. The choice of ionisation level, the choice of initial abundances, as well as the extent of chemical reaction types included are all factors that affect the chemical evolution. The only exception is the inheritance scenario with a low ionisation level, which results in negligible changes compared with the initial abundances, regardless of whether grain-surface chemistry is included. The chemical processing changes the C/O ratios for gas and ice significantly with radius and with model assumptions. For high ionisation levels, chemical processing becomes significant after a few times $10^{5}$~yrs. It is concluded that chemical evolution in the disk midplane needs to be considered in the determination of the volatile composition of planetesimals.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A low-mass protostar's disk-envelope interface: disk-shadowing evidence from ALMA DCO+ observations of VLA1623

Due to instrumental limitations and a lack of disk detections, the structure between the envelope and the rotationally supported disk has been poorly studied. This is now possible with ALMA through observations of CO isotopologs and tracers of freezeout. Class 0 sources are ideal for such studies given their almost intact envelope and young disk. The structure of the disk-envelope interface of the prototypical Class 0 source, VLA1623A which has a confirmed Keplerian disk, is constrained from ALMA observations of DCO+ 3-2 and C18O 2-1. The physical structure of VLA1623 is obtained from the large-scale SED and continuum radiative transfer. An analytic model using a simple network coupled with radial density and temperature profiles is used as input for a 2D line radiative transfer calculation for comparison with the ALMA Cycle 0 12m array and Cycle 2 ACA observations of VLA1623. DCO+ emission shows a clumpy structure bordering VLA1623A's Keplerian disk, suggesting a cold ring-like structure at the disk-envelope interface. The radial position of the observed DCO+ peak is reproduced in our model only if the region's temperature is between 11-16K, lower than expected from models constrained by continuum and SED. Altering the density has little effect on the DCO+ position, but increased density is needed to reproduce the disk traced in C18O. The DCO+ emission around VLA1623A is the product of shadowing of the envelope by the disk. Disk-shadowing causes a drop in the gas temperature outside of the disk on >200AU scales, encouraging deuterated molecule production. This indicates that the physical structure of the disk-envelope interface differs from the rest of the envelope, highlighting the drastic impact that the disk has on the envelope and temperature structure. The results presented here show that DCO+ is an excellent cold temperature tracer.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Photodesorption of H2O, HDO, and D2O ice and its impact on fractionation

The HDO/H2O ratio in interstellar gas is often used to draw conclusions on the origin of water in star-forming regions and on Earth. In cold cores and in the outer regions of protoplanetary disks, gas-phase water comes from photodesorption of water ice. We present fitting formulae for implementation in astrochemical models using photodesorption efficiencies for all water ice isotopologues obtained using classical molecular dynamics simulations. We investigate if the gas-phase HDO/H2O ratio reflects that present in the ice or whether fractionation can occur during photodesorption. Probabilities for the top four monolayers are presented for photodesorption of X (X=H,D) atoms, OX radicals, and X2O and HDO molecules following photodissociation of H2O, D2O, and HDO in H2O amorphous ice at temperatures from 10-100 K. Isotope effects are found for all products: (1) H atom photodesorption probabilities from H2O ice are larger than those for D atom photodesorption from D2O ice by a factor of 1.1; the ratio of H and D photodesorbed upon HDO photodissociation is a factor of 2. This process will enrich the ice in deuterium atoms over time; (2) the OD/OH photodesorption ratio upon D2O and H2O photodissociation is on average a factor of 2, but the ratio upon HDO photodissociation is almost constant at unity for all temperatures; (3) D atoms are more effective in kicking out neighbouring water molecules than H atoms. However, the ratio of the photodesorbed HDO and H2O molecules is equal to the HDO/H2O ratio in the ice, therefore, there is no isotope fractionation upon HDO and H2O photodesorption. Nevertheless, the enrichment of the ice in D atoms due to photodesorption can over time lead to an enhanced HDO/H2O ratio in the ice, and, when photodesorbed, also in the gas. The extent to which the ortho/para ratio of H2O can be modified by the photodesorption process is also discussed. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resolving structure of the disk around HD100546 at 7 mm with ATCA

There is much evidence that planet formation is occurring in the disk around the Herbig Be star HD100546. To learn more about the processes occurring in this disk we conducted high resolution imaging at 43/45 GHz with the Australia Telescope Compact Array (ATCA). Multiple array configurations were used, providing a best spatial resolution of $\sim$ 0.15 arcsec, or 15 AU at HD100546's distance of $\sim$ 100 pc. Significant structure is revealed, but its precise form is dependent on the $u-v$ plane sampling used for the image reconstruction. At a resolution of $\leq$ 30 AU we detected an inner gap in the disk with a radius of $\sim$ 25 AU and a position angle approximately along the known disk major axis. With different weighting, and an achieved resolution of $\sim$ 15 AU, emission appears at the centre and the disk takes on the shape of an incomplete ring, much like a horseshoe, again with a gap radius of $\sim$ 25 AU. The position angle of the disk major axis and its inclination from face-on are determined to be $140^{\circ}\pm5^{\circ}$ and $40^{\circ}\pm5^{\circ}$ respectively. The $\sim$ 25 AU gap radius is confirmed by a null in the real part of the binned visibilities at 320$\pm$10 k$λ$, whilst the non-axisymmetric nature is also confirmed through significant structure in the imaginary component. The emission mechanism at the central peak is most likely to be free-free emission from a stellar or disk wind. Overall our data support the picture of at least one, but probably several, giant planets orbiting HD100546 within 25 AU.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The complex chemistry of outflow cavity walls exposed: the case of low-mass protostars

Complex organic molecules are ubiquitous companions of young low-mass protostars. Recent observations suggest that their emission stems, not only from the traditional hot corino, but also from offset positions. In this work, 2D physicochemical modelling of an envelope-cavity system is carried out. Wavelength-dependent radiative transfer calculations are performed and a comprehensive gas-grain chemical network is used to simulate the physical and chemical structure. The morphology of the system delineates three distinct regions: the cavity wall layer with time-dependent and species-variant enhancements; a torus rich in complex organic ices, but not reflected in gas-phase abundances; and the remaining outer envelope abundant in simpler solid and gaseous molecules. Strongly irradiated regions, such as the cavity wall layer, are subject to frequent photodissociation in the solid phase. Subsequent recombination of the photoproducts leads to frequent reactive desorption, causing gas-phase enhancements of several orders of magnitude. This mechanism remains to be quantified with laboratory experiments. Direct photodesorption is found to be relatively inefficient. If radicals are not produced directly in the icy mantle, the formation of complex organics is impeded. For efficiency, a sufficient number of FUV photons needs to penetrate the envelope; and elevated cool dust temperatures need to enable grain-surface radical mobility. As a result, a high stellar luminosity and a sufficiently wide cavity favor chemical complexity. Furthermore within this paradigm, complex organics are demonstrated to have unique lifetimes and be grouped into early (formaldehyde, ketene, methanol, formic acid, methyl formate, acetic acid, glycolaldehyde) and late (acetaldehyde, dimethyl ether, ethanol) species.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The molecular composition of the planet-forming regions of protoplanetary disks across the luminosity regime

(Abridged) Near- to mid-IR observations of protoplanetary disks show that the inner regions (<10AU) are rich in small organic volatiles (e.g., C2H2 and HCN). Trends in the data suggest that disks around cooler stars (~3000K) are potentially more carbon- and molecule-rich than their hotter counterparts. Our aims are to explore the composition of the planet-forming region of disks around stars from M dwarf to Herbig Ae and compare with the observed trends. Models of the disk physical structure are coupled with a gas-grain chemical network to map the abundances in the planet-forming zone. N2 self shielding, X-ray-induced chemistry, and initial abundances, are investigated. The composition in the 'observable' atmosphere is compared with that in the midplane where the planet-building reservoir resides. M dwarf disk atmospheres are relatively more molecule rich than those for T Tauri or Herbig Ae disks. The weak far-UV flux helps retain this complexity which is enhanced by X-ray-induced ion-molecule chemistry. N2 self shielding has only a small effect and does not explain the higher C2H2/HCN ratios observed towards cooler stars. The models underproduce the OH/H2O column density ratios constrained in Herbig Ae disks, despite reproducing the absolute value for OH: H2O self shielding only increases this discrepency. The disk midplane content is sensitive to the initial main elemental reservoirs. The gas in the inner disk is generally more carbon rich than the midplane ices and is most significant for disks around cooler stars. The atmospheric C/O ratio appears larger than it actually is when calculated using observable tracers only because gas-phase O2 is predicted to be a significant oxygen reservoir. The models suggest that the gas in the inner regions of disks around cooler stars is more carbon rich; however, calculations of the molecular emission are necessary to confirm the observed trends.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

ALMA hints at the presence of two companions in the disk around HD 100546

HD 100546 is a well-studied Herbig Be star-disk system that likely hosts a close-in companion with compelling observational evidence for an embedded protoplanet at 68 AU. We present ALMA observations of the HD 100546 disk which resolve the gas and dust structure at (sub)mm wavelengths. The CO emission (at 345.795 GHz) originates from an extensive molecular disk (390+/-20 AU in radius) whereas the continuum emission is more compact (230+/-20 AU in radius) suggesting radial drift of the mm-sized grains. The CO emission is similar in extent to scattered light images indicating well-mixed gas and um-sized grains in the disk atmosphere. Assuming azimuthal symmetry, a single-component power-law model cannot reproduce the continuum visibilities. The visibilities and images are better reproduced by a double-component model: a compact ring with a width of 21 AU centered at 26 AU and an outer ring with a width of 75+/-3 AU centered at 190+/-3 AU. The influence of a companion and protoplanet on the dust evolution is investigated. The companion at 10 AU facilitates the accumulation of mm-sized grains within a compact ring, ~ 20 - 30 AU, by ~ 10 Myr. The injection of a protoplanet at 1 Myr hastens the ring formation (~ 1.2 Myr) and also triggers the development of an outer ring (~ 100 - 200 AU). These observations provide additional evidence for the presence of a close-in companion and hint at dynamical clearing by a protoplanet in the outer disk.

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

Methanol Along the Path from Envelope to Protoplanetary Disc

Interstellar methanol is considered to be a parent species of larger, more complex organic molecules. A physicochemical simulation of infalling parcels of matter is performed for a low-mass star-forming system to trace the chemical evolution from cloud to disc. An axisymmetric 2D semi-analytic model generates the time-dependent density and velocity distributions, and full continuum radiative transfer is performed to calculate the dust temperature and the UV radiation field at each position as a function of time. A comprehensive gas-grain chemical network is employed to compute the chemical abundances along infall trajectories. Two physical scenarios are studied, one in which the dominant disc growth mechanism is viscous spreading, and another in which continuous infall of matter prevails. The results show that the infall path influences the abundance of methanol entering each type of disc, ranging from complete loss of methanol to an enhancement by a factor of > 1 relative to the prestellar phase. Critical chemical processes and parameters for the methanol chemistry under different physical conditions are identified. The exact abundance and distribution of methanol is important for the budget of complex organic molecules in discs, which will be incorporated into forming planetary system objects such as protoplanets and comets. These simulations show that the comet-forming zone contains less methanol than in the precollapse phase, which is dominantly of prestellar origin, but also with additional layers built up in the envelope during infall. Such intriguing links will soon be tested by upcoming data from the Rosetta mission.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Photodissociation and chemistry of N$_2$ in the circumstellar envelope of carbon-rich AGB stars

The envelopes of AGB stars are irradiated externally by ultraviolet photons; hence, the chemistry is sensitive to the photodissociation of N$_2$ and CO, which are major reservoirs of nitrogen and carbon, respectively. The photodissociation of N$_2$ has recently been quantified by laboratory and theoretical studies. Improvements have also been made for CO photodissociation. For the first time, we use accurate N$_2$ and CO photodissociation rates and shielding functions in a model of the circumstellar envelope of the carbon-rich AGB star, IRC +10216. We use a state-of-the-art chemical model of an AGB envelope, the latest CO and N$_2$ photodissociation data, and a new method for implementing molecular shielding functions in full spherical geometry with isotropic incident radiation. We compare computed column densities and radial distributions of molecules with observations. The transition of N$_2$ $\to$ N (also, CO $\to$ C $\to$ C$^+$) is shifted towards the outer envelope relative to previous models. This leads to different column densities and radial distributions of N-bearing species, especially those species whose formation/destruction processes largely depend on the availability of atomic or molecular nitrogen, for example, C$_n$N ($n$=1, 3, 5), C$_n$N$^-$ ($n$=1, 3, 5), HC$_n$N ($n$=1, 3, 5, 7, 9), H$_2$CN and CH$_2$CN. The chemistry of many species is directly or indirectly affected by the photodissociation of N$_2$ and CO, especially in the outer shell of AGB stars where photodissociation is important. Thus, it is important to include N$_2$ and CO shielding in astrochemical models of AGB envelopes and other irradiated environments. In general, while differences remain between our model of IRC +10216 and the observed molecular column densities, better agreement is found between the calculated and observed radii of peak abundance.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Molecular line emission from a protoplanetary disk irradiated externally by a nearby massive star

Star formation often occurs within or nearby stellar clusters. Irradiation by nearby massive stars can photoevaporate protoplanetary disks around young stars (so-called proplyds) which raises questions regarding the ability of planet formation to take place in these environments. We investigate the two-dimensional physical and chemical structure of a protoplanetary disk surrounding a low-mass (T Tauri) star which is irradiated by a nearby massive O-type star to determine the survivability and observability of molecules in proplyds. Compared with an isolated star-disk system, the gas temperature ranges from a factor of a few (in the disk midplane) to around two orders of magnitude (in the disk surface) higher in the irradiated disk. Although the UV flux in the outer disk, in particular, is several orders of magnitude higher, the surface density of the disk is sufficient for effective shielding of the disk midplane so that the disk remains predominantly molecular in nature. We also find that non-volatile molecules, such as HCN and H2O, are able to freeze out onto dust grains in the disk midplane so that the formation of icy planetesimals, e.g., comets, may also be possible in proplyds. We have calculated the molecular line emission from the disk assuming LTE and determined that multiple transitions of atomic carbon, CO (and isotopologues, 13CO and C18O), HCO+, CN, and HCN may be observable with ALMA, allowing characterization of the gas column density, temperature, and optical depth in proplyds at the distance of Orion (~400 pc).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Chemical Processes in Protoplanetary Disks II. On the importance of photochemistry and X-ray ionization

We investigate the impact of photochemistry and X-ray ionization on the molecular composition of, and ionization fraction in, a protoplanetary disk surrounding a typical T Tauri star. We use a sophisticated physical model, which includes a robust treatment of the radiative transfer of UV and X-ray radiation, and calculate the time-dependent chemical structure using a comprehensive chemical network. In previous work, we approximated the photochemistry and X-ray ionization, here, we recalculate the photoreaction rates using the explicit UV wavelength spectrum and wavelength-dependent reaction cross sections. We recalculate the X-ray ionization rate using our explicit elemental composition and X-ray energy spectrum. We find photochemistry has a larger influence on the molecular composition than X-ray ionization. Observable molecules sensitive to the photorates include OH, HCO+, N2H+, H2O, CO2 and CH3OH. The only molecule significantly affected by the X-ray ionization is N2H+ indicating it is safe to adopt existing approximations of the X-ray ionization rate in typical T Tauri star-disk systems. The recalculation of the photorates increases the abundances of neutral molecules in the outer disk, highlighting the importance of taking into account the shape of the UV spectrum in protoplanetary disks. A recalculation of the photoreaction rates also affects the gas-phase chemistry due to the adjustment of the H/H2 and C+/C ratios. The disk ionization fraction is not significantly affected by the methods adopted to calculate the photochemistry and X-ray ionization. We determine there is a probable 'dead zone' where accretion is suppressed, present in a layer, Z/R <~ 0.1 - 0.2, in the disk midplane, within R \approx 200 AU.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The chemistry of extragalactic carbon stars

Prompted by the ongoing interest in Spitzer Infrared Spectrometer spectra of carbon stars in the Large Magellanic Cloud, we have investigated the circumstellar chemistry of carbon stars in low-metallicity environments. Consistent with observations, our models show that acetylene is particularly abundant in the inner regions of low metallicity carbon-rich AGB stars -- more abundant than carbon monoxide. As a consequence, larger hydrocarbons have higher abundances at the metallicities of the Magellanic Clouds than in stars with solar metallicity. We also find the oxygen and nitrogen chemistry is suppressed at lower metallicity, as expected. Finally, we calculate molecular line emission from carbon stars in the Large and Small Magellanic Cloud and find that several molecules should be readily detectable with the Atacama Large Millimeter Array at Full Science operations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

CH abundance gradient in TMC-1

We observed the 9-cm Lambda-doubling lines of CH along the dense filament of TMC-1. The CH column densities were compared with the total H2 column densities derived using the 2MASS NIR data and previously published SCUBA maps and with OH column densities derived using previous observations with Effelsberg. We also modelled the chemical evolution of TMC-1 adopting physical conditions typical of dark clouds using the UMIST Database for Astrochemistry gas-phase reaction network to aid the interpretation of the observed OH/CH abundance ratios. The CH column density has a clear peak in the vicinity of the cyanopolyyne maximum of TMC-1. The fractional CH abundance relative to H2 increases steadily from the northwestern end of the filament where it lies around 1.0e-8, to the southeast where it reaches a value of 2.0e-8. The OH and CH column densities are well correlated, and we obtained OH/CH abundance ratios of ~ 16 - 20. These values are clearly larger than what has been measured recently in diffuse interstellar gas and is likely to be related to C to CO conversion at higher densities. The good correlation between CH and OH can be explained by similar production and destruction pathways. We suggest that the observed CH and OH abundance gradients are mainly due to enhanced abundances in a low-density envelope which becomes more prominent in the southeastern part and seems to continue beyond the dense filament. An extensive envelope probably signifies an early stage of dynamical evolution, and conforms with the detection of a large CH abundance in the southeastern part of the cloud. The implied presence of other simple forms of carbon in the gas phase provides a natural explanation for the observation of "early-type" molecules in this region.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Chemical evolution of protoplanetary disks - the effects of viscous accretion, turbulent mixing and disk winds

We calculate the chemical evolution of protoplanetary disks considering radial viscous accretion, vertical turbulent mixing and vertical disk winds. We study the effects on the disk chemical structure when different models for the formation of molecular hydrogen on dust grains are adopted. Our gas-phase chemistry is extracted from the UMIST Database for Astrochemistry (Rate06) to which we have added detailed gas-grain interactions. We use our chemical model results to generate synthetic near- and mid-infrared LTE line emission spectra and compare these with recent Spitzer observations. Our results show that if H2 formation on warm grains is taken into consideration, the H2O and OH abundances in the disk surface increase significantly. We find the radial accretion flow strongly influences the molecular abundances, with those in the cold midplane layers particularly affected. On the other hand, we show that diffusive turbulent mixing affects the disk chemistry in the warm molecular layers, influencing the line emission from the disk and subsequently improving agreement with observations. We find that NH3, CH3OH, C2H2 and sulphur-containing species are greatly enhanced by the inclusion of turbulent mixing. We demonstrate that disk winds potentially affect the disk chemistry and the resulting molecular line emission in a similar manner to that found when mixing is included.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Chemical Processes in Protoplanetary Disks

We have developed a high resolution combined physical and chemical model of a protoplanetary disk surrounding a typical T Tauri star. Our aims were to use our model to calculate the chemical structure of disks on small scales (sub-milli-arcsecond in the inner disk for objects at the distance of Taurus, ~ 140 pc) to investigate the various chemical processes thought to be important in disks and to determine potential molecular tracers of each process. Our gas-phase network was extracted from the UMIST Database for Astrochemistry to which we added gas-grain interactions including freeze out and thermal and non-thermal desorption (cosmic-ray induced desorption, photodesorption and X-ray desorption) and a grain-surface network. We find that cosmic-ray induced desorption has the least effect on our disk chemical structure while photodesorption has a significant effect, enhancing the abundances of most gas-phase molecules throughout the disk and affecting the abundances and distribution of HCN, CN and CS, in particular. In the outer disk, we also see enhancements in the abundances of H2O and CO2. X-ray desorption is a potentially powerful mechanism in disks, acting to homogenise the fractional abundances of gas-phase species across the depth and increasing the column densities of most molecules although there remain significant uncertainties in the rates adopted for this process. The addition of grain-surface chemistry enhances the fractional abundances of several small complex organic molecules including CH3OH, HCOOCH3 and CH3OCH3 to potentially observable values (i.e. a fractional abundance of >~ 1.0E-11).

0 citations0 reviewsNo code linkedOpen access