Showing 1–50 of 73 results for author: Millar, T J

A $λ$ 3 mm line survey towards the circumstellar envelope of the carbon-rich AGB star IRC+10216 (CW Leo)

Authors: Juan Tuo, Xiaohu Li, Jixian Sun, Tom J. Millar, Yong Zhang, Jianjie Qiu, Donghui Quan, Jarken Esimbek, Jianjun Zhou, Yu Gao, Qiang Chang, Lin Xiao, Yanan Feng, Zhenzhen Miao, Rong Ma, Ryszard Szczerba, Xuan Fang

Abstract: We present an unbiased $λ$ 3 mm spectral line survey (between 84.5 and 115.8 GHz), conducted by the Purple Mountain Observatory 13.7 meter radio telescope, together with updated modeling results, towards the carbon-rich Asymptotic Giant Branch star, IRC+10216 (CW Leo). A total of 75 spectral lines (96 transitions) are detected, and identified to arise from 19 molecules: C$_2$H, $l$-C$_3$H, C$_4$H,… ▽ More We present an unbiased $λ$ 3 mm spectral line survey (between 84.5 and 115.8 GHz), conducted by the Purple Mountain Observatory 13.7 meter radio telescope, together with updated modeling results, towards the carbon-rich Asymptotic Giant Branch star, IRC+10216 (CW Leo). A total of 75 spectral lines (96 transitions) are detected, and identified to arise from 19 molecules: C$_2$H, $l$-C$_3$H, C$_4$H, CN, C$_3$N, HC$_3$N, HC$_5$N, HCN, HNC, CH$_3$CN, MgNC, CO, $c$-C$_3$H$_2$, SiC$_2$, SiO, SiS, CS, C$_2$S, C$_3$S, and their isotopologues. Among them, one molecular emission line (H$^{13}$CCCN $J=13-12$) is discovered in IRC+10216 for the first time. The excitation temperature, column density, and fractional abundance of the detected species are deduced by assuming they are in local thermodynamic equilibrium. In addition, the isotopic ratios of [$^{12}$C]/[$^{13}$C], [$^{32}$S]/[$^{34}$S], [$^{28}$Si]/[$^{29}$Si], and [$^{12}$C$^{34}$S]/[$^{13}$C$^{32}$S] are obtained and found to be consistent with previous studies. Finally, we summarize all of the 106 species detected in IRC+10216 to date with their observed and modeled column densities for the convenience of future studies. △ Less

Submitted 30 January, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: 71 pages, 39 figures, 10 tables. Accepted for publication in ApJS

Charting Circumstellar Chemistry of Carbon-rich AGB Stars: I. ALMA 3 mm spectral surveys

Authors: R. Unnikrishnan, E. De Beck, L. A. Nyman, H. Olofsson, W. H. T. Vlemmings, D. Tafoya, M. Maercker, S. B. Charnley, M. A. Cordiner, I. de Gregorio, E. Humphreys, T. J. Millar, M. G. Rawlings

Abstract: AGB stars are major contributors to the chemical enrichment of the ISM through nucleosynthesis and extensive mass loss. Most of our current knowledge of AGB atmospheric and circumstellar chemistry, in particular in a C-rich environment, is based on observations of the carbon star IRC+10216. We aim to obtain a more generalised understanding of the chemistry in C-rich AGB CSEs by studying a sample o… ▽ More AGB stars are major contributors to the chemical enrichment of the ISM through nucleosynthesis and extensive mass loss. Most of our current knowledge of AGB atmospheric and circumstellar chemistry, in particular in a C-rich environment, is based on observations of the carbon star IRC+10216. We aim to obtain a more generalised understanding of the chemistry in C-rich AGB CSEs by studying a sample of three carbon stars, IRAS15194-5115, IRAS15082-4808, and IRAS07454-7112, and test the archetypal status often attributed to IRC+10216. We performed spatially resolved, unbiased spectral surveys in ALMA Band 3. We identify a total of 132 rotational transitions from 49 molecular species. There are two main morphologies of the brightness distributions: centrally-peaked (e.g. HCN) and shell-like (e.g. C$_2$H). We estimated the sizes of the molecular emitting regions using azimuthally-averaged radial profiles of the line brightness distributions, and derived abundance estimates. Of the shell distributions, the cyanopolyynes peak at slightly smaller radii than the hydrocarbons, and CN and HNC show the most extended emission. The emitting regions for each species are the smallest for IRAS07454-7112. We find that, within the uncertainties of the analysis, the three stars present similar abundances for most species, also compared to IRC+10216. We find that SiO is more abundant in our three stars compared to IRC+10216. Our estimated isotopic ratios match well the literature values for the sources. The observed circumstellar chemistry appears very similar across our sample and compared to that of IRC+10216, both in terms of the relative location of the emitting regions and molecular abundances. This implies that, to a first approximation, the chemical models tailored to IRC+10216 are able to reproduce the observed chemistry in C-rich envelopes across roughly an order of magnitude in wind density. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 21 pages. 13 figures and 8 tables in the main text. 5 appendices contain additional figures and tables. Appendix tables are available in electronic form at the CDS, along with the reprocessed ALMA cubes and spectra, at http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ or via anonymous ftp to cdsarc.u-strasbg.fr(130.79.128.5). Accepted for publication in Astronomy & Astrophysics

ATOMIUM: Molecular inventory of 17 oxygen-rich evolved stars observed with ALMA

Authors: S. H. J. Wallstrom, T. Danilovich, H. S. P. Muller, C. A. Gottlieb, S. Maes, M. Van de Sande, L. Decin, A. M. S. Richards, A. Baudry, J. Bolte, T. Ceulemans, F. De Ceuster, A. de Koter, I. El Mellah, M. Esseldeurs, S. Etoka, D. Gobrecht, E. Gottlieb, M. Gray, F. Herpin, M. Jeste, D. Kee, P. Kervella, T. Khouri, E. Lagadec , et al. (13 additional authors not shown)

Abstract: The dusty winds of cool evolved stars are a major contributor of the newly synthesised material enriching the Galaxy and future generations of stars. However, the details of the physics and chemistry behind dust formation and wind launching have yet to be pinpointed. Recent spatially resolved observations show the importance of gaining a more comprehensive view of the circumstellar chemistry, but… ▽ More The dusty winds of cool evolved stars are a major contributor of the newly synthesised material enriching the Galaxy and future generations of stars. However, the details of the physics and chemistry behind dust formation and wind launching have yet to be pinpointed. Recent spatially resolved observations show the importance of gaining a more comprehensive view of the circumstellar chemistry, but a comparative study of the intricate interplay between chemistry and physics is still difficult because observational details such as frequencies and angular resolutions are rarely comparable. Aiming to overcome these deficiencies, ATOMIUM is an ALMA Large Programme to study the physics and chemistry of the circumstellar envelopes of a diverse set of oxygen-rich evolved stars under homogeneous observing conditions at three angular resolutions between ~0.02"-1.4". Here we summarize the molecular inventory of these sources, and the correlations between stellar parameters and molecular content. Seventeen oxygen-rich or S-type asymptotic giant branch (AGB) and red supergiant (RSG) stars have been observed in several tunings with ALMA Band 6, targeting a range of molecules to probe the circumstellar envelope and especially the chemistry of dust formation close to the star. We systematically assigned the molecular carriers of the spectral lines and measured their spectroscopic parameters and the angular extent of the emission of each line from integrated intensity maps. Across the ATOMIUM sample, we detect 291 transitions of 24 different molecules and their isotopologues. This includes several first detections in oxygen-rich AGB/RSG stars: PO v=1, SO2 v1=1 and v2=2, and several high energy H2O transitions. We also find several first detections in S-type AGB stars: vibrationally excited HCN v2=2,3 and SiS v=4,5,6, as well as first detections of the molecules SiC, AlCl, and AlF in W Aql... △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: 19 pages plus appendices, forthcoming publication in Astronomy & Astrophysics

The UMIST Database for Astrochemistry 2022

Authors: T J Millar, C Walsh, M Van de Sande, A J Markwick

Abstract: Detailed astrochemical models are a key component to interpret the observations of interstellar and circumstellar molecules since they allow important physical properties of the gas and its evolutionary history to be deduced. We update one of the most widely used astrochemical databases to reflect advances in experimental and theoretical estimates of rate coefficients and to respond to the large i… ▽ More Detailed astrochemical models are a key component to interpret the observations of interstellar and circumstellar molecules since they allow important physical properties of the gas and its evolutionary history to be deduced. We update one of the most widely used astrochemical databases to reflect advances in experimental and theoretical estimates of rate coefficients and to respond to the large increase in the number of molecules detected in space since our last release in 2013. We present the sixth release of the UMIST Database for Astrochemistry (UDfA), a major expansion of the gas-phase chemistry that describes the synthesis of interstellar and circumstellar molecules. Since our last release, we have undertaken a major review of the literature which has increased the number of reactions by over 40% to a total of 8767 and increased the number of species by over 55% to 737. We have made a particular attempt to include many of the new species detected in space over the past decade, including those from the QUIJOTE and GOTHAM surveys, as well as providing references to the original data sources. We use the database to investigate the gas-phase chemistries appropriate to O-rich and C-rich conditions in TMC-1 and to the circumstellar envelope of the C-rich AGB star IRC+10216 and identify successes and failures of gas-phase only models. This update is a significant improvement to the UDfA database. For the dark cloud and C-rich circumstellar envelope models, calculations match around 60% of the abundances of observed species to within an order of magnitude. There are a number of detected species, however, that are not included in the model either because their gas-phase chemistry is unknown or because they are likely formed via surface reactions on icy grains. Future laboratory and theoretical work is needed to include such species in reaction networks. △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: 26 pages, 8 figures, accepted for publication in A&A

Concluding remarks: Faraday Discussion on astrochemistry at high resolution

Authors: T. J. Millar

Abstract: Fifty years on from the first detailed chemical kinetic modelling of astronomical sources, I provide some introductory comments on the history of astrochemistry, summarise some personal views on the topics covered in this discussion meeting, and conclude with some thoughts on its future development. I have left out the jokes. Fifty years on from the first detailed chemical kinetic modelling of astronomical sources, I provide some introductory comments on the history of astrochemistry, summarise some personal views on the topics covered in this discussion meeting, and conclude with some thoughts on its future development. I have left out the jokes. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 13 pages, 2 figures, early on-line publication

Journal ref: Faraday Discussion 'Astrochemistry at high resolution', 2023

Photochemical origin of SiC$_2$ in the circumstellar envelope of carbon-rich AGB stars revealed by ALMA

Authors: Yanan Feng, Xiaohu Li, Tom J. Millar, Ryszard Szczerba, Ke Wang, Donghui Quan, Shengli Qin, Xuan Fang, Juan Tuo, Zhenzhen Miao, Rong Ma, Fengwei Xu, Jingfei Sun, Biwei Jiang, Qiang Chang, Jianchao Yang, Gao-Lei Hou, Fangfang Li, Yong Zhang

Abstract: Whether SiC$_2$ is a parent species, that is formed in the photosphere or as a by-product of high-temperature dust formation, or a daughter species, formed in a chemistry driven by the photodestruction of parent species in the outer envelope, has been debated for a long time. Here, we analyze the ALMA observations of four SiC$_2$ transitions in the CSEs of three C-rich AGB stars (AI Vol, II Lup, a… ▽ More Whether SiC$_2$ is a parent species, that is formed in the photosphere or as a by-product of high-temperature dust formation, or a daughter species, formed in a chemistry driven by the photodestruction of parent species in the outer envelope, has been debated for a long time. Here, we analyze the ALMA observations of four SiC$_2$ transitions in the CSEs of three C-rich AGB stars (AI Vol, II Lup, and RAFGL 4211), and found that SiC$_2$ exhibits an annular, shell-like distribution in these targets, suggesting that SiC$_2$ can be a daughter species in the CSEs of carbon-rich AGB stars. The results can provide important references for future chemical models. △ Less

Submitted 24 July, 2023; originally announced July 2023.

Comments: Accepted in Frontiers in Astronomy and Space Sciences

Disentangling physics and chemistry in AGB outflows: revealing degeneracies when adding complexity

Authors: Marie Van de Sande, Catherine Walsh, Tom J. Millar

Abstract: Observations of the outflows of asymptotic giant branch (AGB) stars continue to reveal their chemical and dynamical complexity. Spherical asymmetries, such as spirals and disks, are prevalent and thought to be caused by binary interaction with a (sub)stellar companion. Furthermore, high density outflows show evidence of dust-gas interactions. The classical chemical model of these outflows - a gas-… ▽ More Observations of the outflows of asymptotic giant branch (AGB) stars continue to reveal their chemical and dynamical complexity. Spherical asymmetries, such as spirals and disks, are prevalent and thought to be caused by binary interaction with a (sub)stellar companion. Furthermore, high density outflows show evidence of dust-gas interactions. The classical chemical model of these outflows - a gas-phase only, spherically symmetric chemical kinetics model - is hence not appropriate for a majority of observed outflows. We have included several physical and chemical advancements step-by-step: a porous density distribution, dust-gas chemistry, and internal UV photons originating from a close-by stellar companion. Now, we combine these layers of complexity into the most chemically and physically advanced chemical kinetics model of AGB outflows to date. By varying over all model parameters, we obtain a holistic view of the outflow's composition and how it (inter)depends on the different complexities. A stellar companion has the largest influence, especially when combined with a porous outflow. We compile sets of gas-phase molecules that trace the importance of dust-gas chemistry and allow us to infer the presence of a companion and porosity of the outflow. This shows that our new chemical model can be used to infer physical and chemical properties of specific outflows, as long as a suitable range of molecules is observed. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: Faraday Discussions 2023, accepted manuscript. 15 pages, 7 figures, 4 tables. Supplementary information can be found at https://pubs.rsc.org/en/content/articlelanding/2023/fd/d3fd00039g

ATOMIUM: Probing the inner wind of evolved O-rich stars with new, highly excited H$_2$O and OH lines

Authors: A. Baudry, K. T. Wong, S. Etoka, A. M. S. Richards, H. S. P. Muller, F. Herpin, T. Danilovich, M. D. Gray, S. Wallstrom, D. Gobrecht, T. Khouri, L. Decin, C. A. Gottlieb, K. M. Menten, W. Homan, T. J. Millar, M. Montarges, B. Pimpanuwat, J. M. C. Plane, P. Kervella

Abstract: Water and the hydroxyl radical are major constituents of the envelope of O-rich late-type stars. Transitions involving energy levels that are highly excited have been observed in both H$_2$O and OH. These and more recently discovered transitions can now be observed at a high sensitivity and angular resolution with the ALMA Array. Spectra and maps of H$_2$O and OH observed with an angular resolutio… ▽ More Water and the hydroxyl radical are major constituents of the envelope of O-rich late-type stars. Transitions involving energy levels that are highly excited have been observed in both H$_2$O and OH. These and more recently discovered transitions can now be observed at a high sensitivity and angular resolution with the ALMA Array. Spectra and maps of H$_2$O and OH observed with an angular resolution of 20 to $\sim$200 mas were obtained at two epochs with the ALMA array. Observations with the Compact Array were also used to check for time variability of water transitions. Radiative transfer models of water were revisited to characterize masing conditions and up-to-date chemical models were used for comparison with our observations. Ten rotational transitions of H$_2$O with energies up to 9000 K were observed in various vibrational states. All but one are new detections in space, and from these we have derived accurate rest frequencies. Hyperfine split $Λ$-doubling transitions in v = 0, J = 27/2 and 29/2 levels of the $^2Π_{3/2}$ state and, $J = 33/2$ and 35/2 of the $^2Π_{1/2}$ state of OH with excitation energies up to 8900 K were also observed. Four of these transitions are new detections in space. Combining our measurements with earlier observations of OH, the v = 0 and v = 1 $Λ$-doubling frequencies have been improved. Our H$_2$O maps show compact emission and extensions up to twelve stellar radii or more. The 268.149 GHz emission line of water in the v$_2$ = 2 state is time variable, tends to be masing with dominant radiative pumping, and is widely excited. The widespread but weaker 262.898 GHz water line in v$_2$ = 1 also shows signs of maser emission. Emission and absorption of both H$_2$O and OH reveal an infall of matter and complex kinematics influenced by binarity. From our observed column densities, we derived OH/H$_2$O abundance ratios in a few stars. △ Less

Submitted 4 May, 2023; originally announced May 2023.

Comments: 41 pages with references and 25 figures in main text. 4 Tables in Appendix A. 10 figures in Appendix B. 2 figures in Appendix C. 14 figures in Appendix D. 5 figures in Appendix E. 2 figures in Appendix F

Journal ref: A&A 674, A125 (2023)

The VLT/SPHERE view of the ATOMIUM cool evolved star sample. I. Overview: Sample characterization through polarization analysis

Authors: M. Montargès, E. Cannon, A. de Koter, T. Khouri, E. Lagadec, P. Kervella, L. Decin, I. McDonald, W. Homan, L. B. F. M. Waters, R. Sahai, C. A. Gottlieb, J. Malfait, S. Maes, B. Pimpanuwat, M. Jeste, T. Danilovich, F. De Ceuster, M. Van de Sande, D. Gobrecht, S. H. J. Wallström, K. T. Wong, I. El Mellah, J. Bolte, F. Herpin , et al. (10 additional authors not shown)

Abstract: Aims. Through the ATOMIUM project, based on an ALMA large program, we aim to present a consistent view of a sample of 17 nearby cool evolved stars (Aymptotic Giant Branch and red supergiant stars). Methods. Here we present VLT/SPHERE-ZIMPOL polarimetric maps obtained in the visible of 14 out of the 17 ATOMIUM sources. They were obtained contemporaneously with the ALMA high spatial resolution dat… ▽ More Aims. Through the ATOMIUM project, based on an ALMA large program, we aim to present a consistent view of a sample of 17 nearby cool evolved stars (Aymptotic Giant Branch and red supergiant stars). Methods. Here we present VLT/SPHERE-ZIMPOL polarimetric maps obtained in the visible of 14 out of the 17 ATOMIUM sources. They were obtained contemporaneously with the ALMA high spatial resolution data. To help interpret the polarized signal, we produced synthetic maps of light scattering by dust, through 3D radiative transfer simulations with the RADMC3D code. Results. The degree of linear polarization (DoLP) observed by ZIMPOL spreads across several optical filters. We infer that it primarily probes dust located just outside of the point spread function, and in or near the plane of the sky, with a total optical depth close to unity in the line of sight, representing only a fraction of the total circumstellar dust. The maximum DoLP ranges from 0.03-0.38 depending on the source, fractions that can be reproduced by our 3D pilot models for grains composed of common dust species. The spatial structure of the DoLP shows a diverse set of shapes. Only for three sources do we note a correlation between the ALMA CO and SiO lines, which trace the gas density, and the DoLP, which traces the dust. Conclusion. The clumpiness of the DoLP and the lack of a consistent correlation between the gas and the dust location show that, in the inner circumstellar environment (CSE), dust formation occurs at very specific sites. This has potential consequences for the derived mass-loss rates and dust-to-gas ratio in the inner region of the CSE. Except for $π^1$~Gru and perhaps GY Aql, we do not detect interactions between the circumstellar wind and the hypothesized companions that shape the wind at larger scales. This suggests that the orbits of any other companions are tilted out of the plane of the sky. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 22 pages, 15 figures, 5 tables

Journal ref: A&A 671, A96 (2023)

Direct D-atom incorporation in radicals: An overlooked pathway for deuterium fractionation

Authors: Nureshan Dias, Ranil M. Gurusinghe, Bernadette M. Broderick, Tom J Millar, Arthur G. Suits

Abstract: Direct D-H exchange in radicals is investigated in a quasi-uniform flow employing chirped pulse mm-wave spectroscopy. Inspired by the H-atom catalyzed isomerization of C3H2 reported in our previous study, D atom reactions with the propargyl (C3H3) radical and its photoproducts were investigated. We observed very efficient D atom enrichment in the photoproducts through an analogous process of D add… ▽ More Direct D-H exchange in radicals is investigated in a quasi-uniform flow employing chirped pulse mm-wave spectroscopy. Inspired by the H-atom catalyzed isomerization of C3H2 reported in our previous study, D atom reactions with the propargyl (C3H3) radical and its photoproducts were investigated. We observed very efficient D atom enrichment in the photoproducts through an analogous process of D addition/H elimination to C3H2 isomers occurring at 40K or below. Cyclic C3HD is the only deuterated isomer observed, consistent with the expected addition/elimination yielding the lowest energy product. The other expected addition/elimination product, deuterated propargyl, is not directly detected, although its presence is inferred by the observations in the latter part of the flow. There, in the high-density region of the flow, we observed both isotopomers of singly deuterated propyne attributed to stabilization of the H + C3H2D or D + C3H3 adducts. The implications of these observations for the deuterium fractionation of hydrocarbon radicals in astrochemical environments is discussed with the support of a monodeuterated chemical kinetic model. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Investigating Anomalous Photochemistry in the Inner Wind of IRC+10216 Through Interferometric Observations of HC$_3$N

Authors: Mark A. Siebert, Marie Van de Sande, Thomas J. Millar, Anthony J. Remijan

Abstract: In recent years, many questions have arisen regarding the chemistry of photochemical products in the carbon-rich winds of evolved stars. To address them, it is imperative to constrain the distributions of such species through high angular resolution interferometric observations covering multiple rotational transitions. We used archival ALMA observations to map rotational lines involving high energ… ▽ More In recent years, many questions have arisen regarding the chemistry of photochemical products in the carbon-rich winds of evolved stars. To address them, it is imperative to constrain the distributions of such species through high angular resolution interferometric observations covering multiple rotational transitions. We used archival ALMA observations to map rotational lines involving high energy levels of cyanoacetylene (HC$_3$N) toward the inner envelope (radius <8"/1000 AU) of the carbon star IRC+10216. The observed lines include the J=28-27, J=30-29, and J=38-37, transitions of HC$_3$N in its ground vibrational state. In contrast to previous observations of linear carbon chains toward this AGB star which show extended, hollow emission at 15"-20" radii (e.g. C$_4$H, C$_6$H, HC$_5$N), the maps of the HC$_3$N lines here show compact morphologies comprising various arcs and density enhancements, with significant emission from gas clumps at an angular distance of ~3" (350 AU) from the central AGB star. We compared visibility sampled non-LTE radiative transfer models with the observed brightness distributions, and derive a fractional abundance with respect to H$_2$ of $10^{-8}$ for HC$_3$N at the radii probed by these lines. These results are consistent with enhanced photochemistry occurring in warm (~200 K) regions of the circumstellar envelope. After application of a specialized chemical model for IRC+10216, we find evidence that the enhanced HC$_3$N abundances in the inner wind are most likely due to a solar-type binary companion initiating photochemistry in this region. △ Less

Submitted 7 December, 2022; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: 17 pages, 9 figures, 2 tables. Accepted for publication in ApJ

Fevering Interstellar Ices Have More CH3OD

Authors: Beatrice M. Kulterer, Maria N. Drozdovskaya, Stefano Antonellini, Catherine Walsh, Tom J. Millar

Abstract: 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 c… ▽ More 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. △ Less

Submitted 3 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ACS Earth and Space Chemistry; 56 pages, 7 figures

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

Authors: Takashi Tsukagoshi, Hideko Nomura, Takayuki Muto, Ryohei Kawabe, Kazuhiro D. Kanagawa, Satoshi Okuzumi, Shigeru Ida, Catherine Walsh, Tom J. Millar, Sanemichi Z. Takahashi, Jun Hashimoto, Taichi Uyama, Motohide Tamura

Abstract: 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 reconstr… ▽ More 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. △ Less

Submitted 1 February, 2022; originally announced February 2022.

Comments: 17pages, 12 figures, Accepted for publication in The Astrophysical Journal

ATOMIUM: ALMA tracing the origins of molecules in dust forming oxygen rich M-type stars: Motivation, sample, calibration, and initial results

Authors: C. A. Gottlieb, L. Decin, A. M. S. Richards, F. De Ceuster, W. Homan, S. H. J. Wallstrom, T. Danilovich, T. J. Millar, M. Montarges, K. T. Wong, I. McDonald, A. Baudry, J. Bolte, E. Cannon, E. De Beck, A. de Koter, I. El Mellah, S. Etoka, D. Gobrecht, M. Gray, F. Herpin, M. Jeste, P. Kervella, T. Khouri, E. Lagadec , et al. (11 additional authors not shown)

Abstract: This overview paper presents ATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter-submillimeter Array (ALMA). The goal of ATOMIUM is to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical e… ▽ More This overview paper presents ATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter-submillimeter Array (ALMA). The goal of ATOMIUM is to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical environments. 17 oxygen-rich AGB and RSG stars spanning a range in (circum)stellar parameters and evolutionary phases were observed in a homogeneous observing strategy allowing for an unambiguous comparison. Data were obtained between 213.83 and 269.71 GHz at high (0.025-0.050 arcsec), medium (0.13-0.24 arcsec), and low (about 1 arcsec) angular resolution. The sensitivity per 1.3 km/s channel was 1.5-5 mJy/beam. 13 molecules were designated as primary molecules in the survey: CO, SiO, AlO, AlOH, TiO, TiO2, HCN, SO, SO2, SiS, CS, H2O, and NaCl. The scientific motivation, survey design, sample properties, data reduction, and an overview of the data products are described; and we highlight one scientific result - the wind kinematics of the ATOMIUM sources. The ATOMIUM sources often have a slow wind acceleration, and a fraction of the gas reaches a velocity which can be up to a factor of two times larger than previously reported terminal velocities assuming isotropic expansion, and the wind kinematic profiles establish that the radial velocity described by the momentum equation for a spherical wind structure cannot capture the complexity of the velocity field. In 15 sources, some molecular transitions other than 12CO v=0 J=2-1 reach a higher outflow velocity, with a spatial emission zone that is often greater than 30 stellar radii, but much less than the extent of CO. Binary interaction with a (sub)stellar companion might (partly) explain the non-monotonic behaviour of the projected velocity field. △ Less

Submitted 13 December, 2021; v1 submitted 8 December, 2021; originally announced December 2021.

Journal ref: A&A 660, A94 (2022)

The impact of stellar companion UV photons on the chemistry of the circumstellar environments of AGB stars

Authors: M. Van de Sande, T. J. Millar

Abstract: Spherical asymmetries are prevalent within the outflows of AGB stars. Since binary interaction with a stellar or planetary companion is thought to be the underlying mechanism behind large-scale structures, we included the effects of UV radiation originating from a stellar companion in our chemical kinetics model. The one-dimensional model provides a first approximation of its effects on the chemis… ▽ More Spherical asymmetries are prevalent within the outflows of AGB stars. Since binary interaction with a stellar or planetary companion is thought to be the underlying mechanism behind large-scale structures, we included the effects of UV radiation originating from a stellar companion in our chemical kinetics model. The one-dimensional model provides a first approximation of its effects on the chemistry throughout the outflow. The presence of a close-by stellar companion can strongly influence the chemistry within the entire outflow. Its impact depends on the intensity of the radiation (set by the stellar radius and blackbody temperature) and on the extinction the UV radiation experiences (set by the outflow density, density structure, and assumed radius of dust formation). Parent species can be photodissociated by the companion, initiating a rich photon-driven chemistry in the inner parts of the outflow. The outcome depends on the balance between two-body reactions and photoreactions. If two-body reactions dominate, chemical complexity within the outflow increases. This can make the abundance profiles of daughters appear like those of parents, with a larger inner abundance and a gaussian decline. If photoreactions dominate, the outflow can appear molecule-poor. We model three stellar companions. The impact of a red dwarf companion is limited. Solar-like companions show the largest effect, followed by a white dwarf. A stellar companion can also lead to the formation of unexpected species. The outflow's molecular content, especially combined with abundance profiles, can indicate a stellar companion's presence. Our results pave the way for further outflow-specific (three-dimensional) model development. △ Less

Submitted 23 November, 2021; v1 submitted 9 November, 2021; originally announced November 2021.

Comments: Accepted for publication in the Monthly Notices of the Royal Astronomical Society. Main paper: 20 pages, 12 figures, 4 tables. Supplementary material: 42 pages, 66 figures, 2 tables. Revision: corrections to legibility of the label of Fig. A2

ATOMIUM: Halide molecules around the S-type AGB star W Aquilae

Authors: T. Danilovich, M. Van de Sande, J. M. C. Plane, T. J. Millar, P. Royer, M. A. Amor, K. Hammami, L. Decock, C. A. Gottlieb, L. Decin, A. M. S. Richards, E. De Beck, A. Baudry, J. Bolte, E. Cannon, F. De Ceuster, A. de Koter, S. Etoka, D. Gobrecht, M. Gray, F. Herpin, W. Homan, M. Jeste, P. Kervella, T. Khouri , et al. (14 additional authors not shown)

Abstract: S-type asymptotic giant branch (AGB) stars are thought to be intermediates in the evolution of oxygen- to carbon-rich AGB stars. The chemical compositions of their circumstellar envelopes are also intermediate, but have not been studied in as much detail as their carbon- and oxygen-rich counterparts. We aim to determine the abundances of AlCl and AlF from rotational lines, which have been observed… ▽ More S-type asymptotic giant branch (AGB) stars are thought to be intermediates in the evolution of oxygen- to carbon-rich AGB stars. The chemical compositions of their circumstellar envelopes are also intermediate, but have not been studied in as much detail as their carbon- and oxygen-rich counterparts. We aim to determine the abundances of AlCl and AlF from rotational lines, which have been observed for the first time towards an S-type AGB star, W Aql. In combination with models based on PACS observations, we aim to update our chemical kinetics network based on these results. We analyse ALMA observations towards W Aql of AlCl in the ground and first two vibrationally excited states and AlF in the ground vibrational state. Using radiative transfer models, we determine the abundances and spatial abundance distributions of Al$^{35}$Cl, Al$^{37}$Cl, and AlF. We also model HCl and HF emission and compare these models to PACS spectra to constrain the abundances of these species. AlCl is found in clumps very close to the star, with emission confined within 0.1$^{\prime\prime}$ of the star. AlF emission is more extended, with faint emission extending 0.2$^{\prime\prime}$ to 0.6$^{\prime\prime}$ from the continuum peak. We find peak abundances, relative to H$_2$, of $1.7\times 10^{-7}$ for Al$^{35}$Cl, $7\times 10^{-8}$ for Al$^{37}$Cl and $1\times 10^{-7}$ for AlF. From the PACS spectra, we find abundances of $9.7\times 10^{-8}$ and $\leq 10^{-8}$, relative to H$_2$, for HCl and HF, respectively. The AlF abundance exceeds the solar F abundance, indicating that fluorine synthesised in the AGB star has already been dredged up to the surface of the star and ejected into the circumstellar envelope. From our analysis of chemical reactions in the wind, we conclude that AlF may participate in the dust formation process, but we cannot fully explain the rapid depletion of AlCl seen in the wind. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 655, A80 (2021)

High Spatial Resolution Observations of Molecular Lines towards the Protoplanetary Disk around TW Hya with ALMA

Authors: Hideko Nomura, Takashi Tsukagoshi, Ryohei Kawabe, Takayuki Muto, Kazuhiro D. Kanagawa, Yuri Aikawa, Eiji Akiyama, Satoshi Okuzumi, Shigeru Ida, Seokho Lee, Catherine Walsh, T. J. Millar

Abstract: We present molecular line observations of 13CO and C18O J=3-2, CN N = 3 - 2, and CS J=7-6 lines in the protoplanetary disk around TW Hya at a high spatial resolution of ~9 au (angular resolution of 0.15''), using the Atacama Large Millimeter/Submillimeter Array. A possible gas gap is found in the deprojected radial intensity profile of the integrated C18O line around a disk radius of ~58 au, sligh… ▽ More We present molecular line observations of 13CO and C18O J=3-2, CN N = 3 - 2, and CS J=7-6 lines in the protoplanetary disk around TW Hya at a high spatial resolution of ~9 au (angular resolution of 0.15''), using the Atacama Large Millimeter/Submillimeter Array. A possible gas gap is found in the deprojected radial intensity profile of the integrated C18O line around a disk radius of ~58 au, slightly beyond the location of the au-scale dust clump at ~52 au, which resembles predictions from hydrodynamic simulations of planet-disk interaction. In addition, we construct models for the physical and chemical structure of the TW Hya disk, taking account of the dust surface density profile obtained from high spatial resolution dust continuum observations. As a result, the observed flat radial profile of the CN line intensities is reproduced due to a high dust-to-gas surface density ratio inside ~20 au. Meanwhile, the CO isotopologue line intensities trace high temperature gas and increase rapidly inside a disk radius of ~30 au. A model with either CO gas depletion or depletion of gas-phase oxygen elemental abundance is required to reproduce the relatively weak CO isotopologue line intensities observed in the outer disk, consistent with previous atomic and molecular line observations towards the TW Hya disk. {Further observations of line emission of carbon-bearing species, such as atomic carbon and HCN, with high spatial resolution would help to better constrain the distribution of elemental carbon abundance in the disk gas. △ Less

Submitted 3 May, 2021; originally announced May 2021.

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

Atomium: The astounding complexity of the near circumstellar environment of the M-type AGB star R Hydrae. I. Morpho-kinematical interpretation of CO and SiO emission

Authors: Ward Homan, Bannawit Pimpanuwat, Fabrice Herpin, Taissa Danilovich, Iain McDonald, Sofia H. J. Wallström, Anita M. S. Richards, Alain Baudry, Raghvendra Sahai, Tom J. Millar, Alex de Koter, C. A. Gottlieb, Pierre Kervella, Miguel Montargès, Marie Van de Sande, Leen Decin, Albert Zijlstra, Sandra Etoka, Manali Jeste, Holger S. P. Müller, Silke Maes, Jolien Malfait, Karl Menten, John Plane, Kelvin Lee , et al. (14 additional authors not shown)

Abstract: Evolved low- to intermediate-mass stars are known to shed their gaseous envelope into a large, dusty, molecule-rich circumstellar nebula which typically develops a high degree of structural complexity. Most of the large-scale, spatially correlated structures in the nebula are thought to originate from the interaction of the stellar wind with a companion. As part of the Atomium large programme, we… ▽ More Evolved low- to intermediate-mass stars are known to shed their gaseous envelope into a large, dusty, molecule-rich circumstellar nebula which typically develops a high degree of structural complexity. Most of the large-scale, spatially correlated structures in the nebula are thought to originate from the interaction of the stellar wind with a companion. As part of the Atomium large programme, we observed the M-type asymptotic giant branch (AGB) star R Hydrae with ALMA. The morphology of the inner wind of R Hya, which has a known companion at ~3500 au, was determined from maps of CO and SiO obtained at high angular resolution. A map of the CO emission reveals a multi-layered structure consisting of a large elliptical feature at an angular scale of ~10'' that is oriented along the north-south axis. The wind morphology within the elliptical feature is dominated by two hollow bubbles. The bubbles are on opposite sides of the AGB star and lie along an axis with a position angle of ~115 deg. Both bubbles are offset from the central star, and their appearance in the SiO channel maps indicates that they might be shock waves travelling through the AGB wind. An estimate of the dynamical age of the bubbles yields an age of the order of 100 yr, which is in agreement with the previously proposed elapsed time since the star last underwent a thermal pulse. When the CO and SiO emission is examined on subarcsecond angular scales, there is evidence for an inclined, differentially rotating equatorial density enhancement, strongly suggesting the presence of a second nearby companion. The position angle of the major axis of this disc is ~70 deg in the plane of the sky. We tentatively estimate that a lower limit on the mass of the nearby companion is ~0.65 Msol on the basis of the highest measured speeds in the disc and the location of its inner rim at ~6 au from the AGB star. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: 21 pages, 23 figures

Journal ref: A&A 651, A82 (2021)

The Role of Ultraviolet Photons in Circumstellar Astrochemistry

Authors: T. J. Millar

Abstract: Stars with masses between 1 and 8 solar masses (M$_\odot$) lose large amounts of material in the form of gas and dust in the late stages of stellar evolution, during their Asymptotic Giant Branch phase. Such stars supply up to 35% of the dust in the interstellar medium and thus contribute to the material out of which our solar system formed. In addition, the circumstellar envelopes of these stars… ▽ More Stars with masses between 1 and 8 solar masses (M$_\odot$) lose large amounts of material in the form of gas and dust in the late stages of stellar evolution, during their Asymptotic Giant Branch phase. Such stars supply up to 35% of the dust in the interstellar medium and thus contribute to the material out of which our solar system formed. In addition, the circumstellar envelopes of these stars are sites of complex, organic chemistry with over 80 molecules detected in them. We show that internal ultraviolet photons, either emitted by the star itself or from a close-in, orbiting companion, can significantly alter the chemistry that occurs in the envelopes particularly if the envelope is clumpy in nature. At least for the cases explored here, we find that the presence of a stellar companion, such as a white dwarf star, the high flux of UV photons destroys H$_2$O in the inner regions of carbon-rich AGB stars to levels below those observed and produces species such as C$^+$ deep in the envelope in contrast to the expectations of traditional descriptions of circumstellar chemistry. △ Less

Submitted 29 December, 2020; originally announced December 2020.

Comments: To be published in Chinese Journal of Chemical Physics

Chemical modelling of dust-gas chemistry within AGB outflows III. Photoprocessing of the ice and return to the ISM

Authors: M. Van de Sande, C. Walsh, T. J. Millar

Abstract: To explain the properties of dust in the interstellar medium (ISM), the presence of a refractory organic mantle is necessary. The outflows of AGB stars are among the main contributors of stellar dust to the ISM. We present the first study of the refractory organic contribution of AGB stars to the ISM. Based on laboratory experiments, we included a new reaction in our extended chemical kinetics mod… ▽ More To explain the properties of dust in the interstellar medium (ISM), the presence of a refractory organic mantle is necessary. The outflows of AGB stars are among the main contributors of stellar dust to the ISM. We present the first study of the refractory organic contribution of AGB stars to the ISM. Based on laboratory experiments, we included a new reaction in our extended chemical kinetics model: the photoprocessing of volatile complex ices into inert refractory organic material. The refractory organic feedback of AGB outflows to the ISM is estimated using observationally motivated parent species and grids of models of C-rich and O-rich outflows. Refractory organic material is mainly inherited from the gas phase through accretion onto the dust and subsequent photoprocessing. Grain-surface chemistry, initiated by photodissociation of ices, produces only a minor part and takes place in a sub-monolayer regime in almost all outflows. The formation of refractory organic material increases with outflow density and depends on the initial gas-phase composition. While O-rich dust is negligibly covered by refractory organics, C-rich dust has an average coverage of $3-9\%$, but can be as high as $8-22\%$. Although C-rich dust does not enter the ISM bare, its average coverage is too low to influence its evolution in the ISM or significantly contribute to the coverage of interstellar dust. This study opens up questions on the coverage of other dust-producing environments. It highlights the need for an improved understanding of dust formation and for models specific to density structures within the outflow. △ Less

Submitted 23 November, 2020; originally announced November 2020.

Comments: 18 pages, 16 figures, 6 tables, accepted for for publication in the Monthly Notices of the Royal Astronomical Society

Species cycling and the enhancement of ammonia in prestellar cores

Authors: Azrael A. von Procházka, T. J. Millar

Abstract: The quantity of NH3 produced on grain surfaces in the prestellar core is thought to be one of the determining factors regarding the chemical complexity achievable at later stages of stellar birth. In order to investigate how this quantity might be influenced by the gas-grain cycling of molecular material within the cloud, we employ a modified rates gas-grain chemical code and follow the time-depen… ▽ More The quantity of NH3 produced on grain surfaces in the prestellar core is thought to be one of the determining factors regarding the chemical complexity achievable at later stages of stellar birth. In order to investigate how this quantity might be influenced by the gas-grain cycling of molecular material within the cloud, we employ a modified rates gas-grain chemical code and follow the time-dependent chemistry of NH3 as the system evolves. Our models incorporate an updated version of the most recent UDfA network of reaction rate coefficients, desorption from the grains through standard thermal and non-thermal processes, and physisorbed and chemisorbed binding of atomic and molecular hydrogen to a population of carbonaceous and siliceous grains. We find that 1.) observable abundances of NH3 can exist in the gas phase of our models at early times when the N atom is derived from CN via an efficient early-time hydrocarbon chemistry, 2.) a time-dependent gradient exists in the observational agreement between different species classes in our models, consistent with possible physical substructures within the TMC-1 Cyanopolyyne Peak, and 3.) the gaseous and solid-state abundances of NH3 are sensitive to the presence of gas-grain cycling within the system. Our results suggest that the degree of chemical complexity achievable at later stages of the cloud's chemical evolution is indeed influenced by the manner in which the gas-grain cycling occurs. △ Less

Submitted 14 November, 2020; originally announced November 2020.

Comments: 15 pages, 7 figures, to be published in Monthly Notices of the Royal Astronomical Society

Report number: Report-no: MN-20-3166-MJ.R1

Atomium: A high-resolution view on the highly asymmetric wind of the AGB star Pi1 Gruis. I. First detection of a new companion and its effect on the inner wind

Authors: Ward Homan, Miguel Montarges, Bannawit Pimpanuwat, Anita M. S. Richards, Sofia H. J. Wallstrom, Pierre Kervella, Leen Decin, Albert Zijlstra, Taissa Danilovich, Alex de Koter, Karl Menten, Raghvendra Sahai, John Plane, Kelvin Lee, Rens Waters, Alain Baudry, Ka Tat Wong, Tom J. Millar, Marie Van de Sande, Eric Lagadec, David Gobrecht, Jeremy Yates, Daniel Price, Emily Cannon, Jan Bolte , et al. (13 additional authors not shown)

Abstract: The nebular circumstellar environments of cool evolved stars are known to harbour a rich morphological complexity of gaseous structures on different length scales. A large part of these density structures are thought to be brought about by the interaction of the stellar wind with a close companion. The S-type asymptotic giant branch star Pi1 Gruis, which has a known companion at ~440 au and is tho… ▽ More The nebular circumstellar environments of cool evolved stars are known to harbour a rich morphological complexity of gaseous structures on different length scales. A large part of these density structures are thought to be brought about by the interaction of the stellar wind with a close companion. The S-type asymptotic giant branch star Pi1 Gruis, which has a known companion at ~440 au and is thought to harbour a second, closer-by (<10 au) companion, was observed with the Atacama Large Millimeter/submillimeter Array as part of the ATOMIUM Large programme. In this work, the brightest CO, SiO, and HCN molecular line transitions are analysed. The continuum map shows two maxima, separated by 0.04'' (6 au). The CO data unambiguously reveal that Pi1 Gru's circumstellar environment harbours an inclined, radially outflowing, equatorial density enhancement. It contains a spiral structure at an angle of 38+/-3 deg with the line-of-sight. The HCN emission in the inner wind reveals a clockwise spiral, with a dynamical crossing time of the spiral arms consistent with a companion at a distance of 0.04'' from the asymptotic giant branch star, which is in agreement with the position of the secondary continuum peak. The inner wind dynamics imply a large acceleration region, consistent with a beta-law power of ~6. The CO emission suggests that the spiral is approximately Archimedean within 5'', beyond which this trend breaks down as the succession of the spiral arms becomes less periodic. The SiO emission at scales smaller than 0.5'' exhibits signatures of gas in rotation, which is found to fit the expected behaviour of gas in the wind-companion interaction zone. An investigation of SiO maser emission reveals what could be a stream of gas accelerating from the surface of the AGB star to the companion. Using these dynamics, we have tentatively derived an upper limit on the companion mass to be ~1.1 Msol. △ Less

Submitted 12 October, 2020; originally announced October 2020.

Comments: 19 pages, 20 figures

Journal ref: A&A 644, A61 (2020)

(Sub)stellar companions shape the winds of evolved stars

Authors: L. Decin, M. Montargès, A. M. S. Richards, C. A. Gottlieb, W. Homan, I. McDonald, I. El Mellah, T. Danilovich, S. H. J. Wallström, A. Zijlstra, A. Baudry, J. Bolte, E. Cannon, E. De Beck, F. De Ceuster, A. de Koter, J. De Ridder, S. Etoka, D. Gobrecht, M. Gray, F. Herpin, M. Jeste, E. Lagadec, P. Kervella, T. Khouri , et al. (10 additional authors not shown)

Abstract: Binary interactions dominate the evolution of massive stars, but their role is less clear for low- and intermediate-mass stars. The evolution of a spherical wind from an asymptotic giant branch (AGB) star into a nonspherical planetary nebula (PN) could be due to binary interactions. We observed a sample of AGB stars with the Atacama Large Millimeter/submillimeter Array (ALMA) and found that their… ▽ More Binary interactions dominate the evolution of massive stars, but their role is less clear for low- and intermediate-mass stars. The evolution of a spherical wind from an asymptotic giant branch (AGB) star into a nonspherical planetary nebula (PN) could be due to binary interactions. We observed a sample of AGB stars with the Atacama Large Millimeter/submillimeter Array (ALMA) and found that their winds exhibit distinct nonspherical geometries with morphological similarities to planetary nebulae (PNe). We infer that the same physics shapes both AGB winds and PNe; additionally, the morphology and AGB mass-loss rate are correlated. These characteristics can be explained by binary interaction. We propose an evolutionary scenario for AGB morphologies that is consistent with observed phenomena in AGB stars and PNe. △ Less

Submitted 28 September, 2020; v1 submitted 24 September, 2020; originally announced September 2020.

Comments: 19 pages main journal, 97 pages Supplementary Information

Journal ref: Science, 2000, Vol 369, Issue 6510, p. 1497

A detailed view on the circumstellar environment of the M-type AGB star EP Aquarii I. High-resolution ALMA and SPHERE observations

Authors: Ward Homan, Emily Cannon, Miguel Montargès, Anita M. S. Richards, Tom J. Millar, Leen Decin

Abstract: Cool evolved stars are known to be significant contributors to the enrichment of the interstellar medium through their dense and dusty stellar winds. High resolution observations of these outflows have shown them to possess high degrees of morphological complexity. We observed the asymptotic giant branch (AGB) star EP Aquarii with ALMA in band 6 and VLT/SPHERE/ZIMPOL in four filters the visible. B… ▽ More Cool evolved stars are known to be significant contributors to the enrichment of the interstellar medium through their dense and dusty stellar winds. High resolution observations of these outflows have shown them to possess high degrees of morphological complexity. We observed the asymptotic giant branch (AGB) star EP Aquarii with ALMA in band 6 and VLT/SPHERE/ZIMPOL in four filters the visible. Both instruments had an angular resolution of 0.025''. These are follow-up observations to the lower-resolution 2016 ALMA analysis of EP Aquarii, which revealed that its wind possesses a nearly face-on, spiral-harbouring equatorial density enhancement, with a nearly pole-on bi-conical outflow. At the base of the spiral, the SiO emission revealed a distinct emission void approximately 0.4'' to the west of the continuum brightness peak, which was proposed to be linked to the presence of a companion. The new ALMA data better resolve the inner wind and reveal that its morphology as observed in CO is consistent with hydrodynamical companion-induced perturbations. Assuming that photodissociation by the UV-field of the companion is responsible for the emission void in SiO, we deduced the spectral properties of the tentative companion from the size of the hole. We conclude that the most probable companion candidate is a white dwarf with a mass between 0.65 and 0.8 Msol , though a solar-like companion could not be definitively excluded. The radial SiO emission shows periodic, low-amplitude perturbations. We tentatively propose that they could be the consequence of the interaction of the AGB wind with another much closer low-mass companion. The polarised SPHERE/ZIMPOL data show a circular signal surrounding the AGB star with a radius of 0.1''... △ Less

Submitted 19 August, 2020; originally announced August 2020.

Comments: 24 pages, 25 Figures

Journal ref: A&A 642, A93 (2020)

Measurements of Low Temperature Rate Coefficients for the Reaction of CH with CH$_{2}$O and Application to Dark Cloud and AGB Stellar Wind Models

Authors: Niclas A. West, Tom J. Millar, Marie Van de Sande, Edward Rutter, Mark A. Blitz, Leen Decin, Dwayne E. Heard

Abstract: Rate coefficients have been measured for the reaction of CH radicals with formaldehyde, CH$_{2}$O, over the temperature range 31 - 133 K using a pulsed Laval nozzle apparatus combined with pulsed laser photolysis and laser induced fluorescence spectroscopy. The rate coefficients are very large and display a distinct decrease with decreasing temperature below 70 K, although classical collision rate… ▽ More Rate coefficients have been measured for the reaction of CH radicals with formaldehyde, CH$_{2}$O, over the temperature range 31 - 133 K using a pulsed Laval nozzle apparatus combined with pulsed laser photolysis and laser induced fluorescence spectroscopy. The rate coefficients are very large and display a distinct decrease with decreasing temperature below 70 K, although classical collision rate theory fails to reproduce this temperature dependence. The measured rate coefficients have been parameterized and used as input for astrochemical models for both dark cloud and AGB stellar outflow scenarios. The models predict a distinct change (up to a factor of two) in the abundance of ketene, H$_{2}$CCO, which is the major expected molecular product of the CH + CH$_{2}$O reaction. △ Less

Submitted 17 October, 2019; originally announced October 2019.

Comments: 17 pages, 9 Figures, 5 Tables, Accepted for publication in The Astrophysical Journal

Discovery of an au-scale excess in millimeter emission from the protoplanetary disk around TW Hya

Authors: Takashi Tsukagoshi, Takayuki Muto, Hideko Nomura, Ryohei Kawabe, Kazuhiro D. Kanagawa, Satoshi Okuzumi, Shigeru Ida, Catherine Walsh, Tom J. Millar, Sanemichi Z. Takahashi, Jun Hashimoto, Taichi Uyama, Motohide Tamura

Abstract: We report the detection of an excess in dust continuum emission at 233~GHz (1.3~mm in wavelength) in the protoplanetary disk around TW~Hya revealed through high-sensitivity observations at $\sim$3~au resolution with the Atacama Large Millimeter/submillimeter Array (ALMA). The sensitivity of the 233~GHz image has been improved by a factor of 3 with regard to that of our previous cycle 3 observation… ▽ More We report the detection of an excess in dust continuum emission at 233~GHz (1.3~mm in wavelength) in the protoplanetary disk around TW~Hya revealed through high-sensitivity observations at $\sim$3~au resolution with the Atacama Large Millimeter/submillimeter Array (ALMA). The sensitivity of the 233~GHz image has been improved by a factor of 3 with regard to that of our previous cycle 3 observations. The overall structure is mostly axisymmetric, and there are apparent gaps at 25 and 41 au as previously reported. The most remarkable new finding is a few au-scale excess emission in the south-west part of the protoplanetary disk. The excess emission is located at 52 au from the disk center and is 1.5 times brighter than the surrounding protoplanetary disk at a significance of 12$σ$. We performed a visibility fitting to the extracted emission after subtracting the axisymmetric protoplanetary disk emission and found that the inferred size and the total flux density of the excess emission are 4.4$\times$1.0~au and 250~$μ$Jy, respectively. The dust mass of the excess emission corresponds to 0.03~$M_\oplus$ if a dust temperature of 18~K is assumed. Since the excess emission can also be marginally identified in the Band 7 image at almost the same position, the feature is unlikely to be a background source. The excess emission can be explained by a dust clump accumulated in a small elongated vortex or a massive circumplanetary disk around a Neptune mass forming-planet. △ Less

Submitted 11 September, 2019; v1 submitted 20 May, 2019; originally announced May 2019.

Comments: 9pages, 5 figures, Accepted for publication in The Astrophysical Journal Letters

Dust continuum emission and the upper limit fluxes of sub-millimeter water lines of the protoplanetary disk around HD 163296 observed by ALMA

Authors: Shota Notsu, Eiji Akiyama, Alice Booth, Hideko Nomura, Catherine Walsh, Tomoya Hirota, Mitsuhiko Honda, Takashi Tsukagoshi, T. J. Millar

Abstract: In this paper, we analyse the upper limit fluxes of sub-millimeter ortho-H$_{2}$$^{16}$O 321 GHz, para-H$_{2}$$^{18}$O 322 GHz, and HDO 335 GHz lines from the protoplanetary disk around the Herbig Ae star HD 163296, using the Atacama Large Millimeter/Submillimeter Array (ALMA). These water lines are considered to be the best candidate sub-millimeter lines to locate the position of the H$_{2}$O sno… ▽ More In this paper, we analyse the upper limit fluxes of sub-millimeter ortho-H$_{2}$$^{16}$O 321 GHz, para-H$_{2}$$^{18}$O 322 GHz, and HDO 335 GHz lines from the protoplanetary disk around the Herbig Ae star HD 163296, using the Atacama Large Millimeter/Submillimeter Array (ALMA). These water lines are considered to be the best candidate sub-millimeter lines to locate the position of the H$_{2}$O snowline, on the basis of our previous model calculations. We compare the upper limit fluxes with the values calculated by our models with dust emission included, and we constrain the line emitting region and the dust opacity from the observations. We conclude that, if the outer edge of the region with high water vapor abundance and if the position of the water snowline are beyond 8 au also, the mm dust opacity will have a value larger than 2.0 cm$^{2}$ g$^{-1}$. In addition, the position of the water snowline will be inside 20 au, if the mm dust opacity is 2.0 cm$^{2}$ g$^{-1}$. Future observations of the dust continuum emission at higher angular resolution and sub-millimeter water lines with longer observation time are required to clarify the detailed structures and the position of the H$_{2}$O snowline in the disk midplane. △ Less

Submitted 26 February, 2019; originally announced February 2019.

Comments: 17 pages, 8 figures, and 3 tables are contained in this paper. It was accepted for publication in The Astrophysical Journal (ApJ) on February 26th, 2019

The Role of Internal Photons on the Chemistry of the Circumstellar Envelopes of AGB Stars

Authors: M Van de Sande, T J Millar

Abstract: Recent high spatial resolution observations of gas and dust in the circumstellar envelopes (CSEs) of AGB stars indicate morphologies much more complex than the smooth density distributions generated by spherically symmetric, constant mass loss rates. In particular, the observation of spiral arcs and disks indicate the likely presence of a binary companion which in some cases give rise to the UV ph… ▽ More Recent high spatial resolution observations of gas and dust in the circumstellar envelopes (CSEs) of AGB stars indicate morphologies much more complex than the smooth density distributions generated by spherically symmetric, constant mass loss rates. In particular, the observation of spiral arcs and disks indicate the likely presence of a binary companion which in some cases give rise to the UV photons detected by GALEX. In this Article, we extend our recent model of the chemistry in a clumpy, porous CSE around an AGB star to include the influence of stellar blackbody photons on the CSE chemistry. Our results indicate that internal photons, in a clumpy, porous CSE, can alter chemistry within a few stellar radii and, for some molecules, alter abundances out to several hundred stellar radii. They further suggest that harder radiation from companion stars or accretion disks will have a substantial impact on chemistry in the dust formation zones and inner CSEs of AGB stars. △ Less

Submitted 28 February, 2019; v1 submitted 1 February, 2019; originally announced February 2019.

Comments: 10 pages, 4 figures, 2 tables, references updated

The effect of carbon grain destruction on the chemical structure of protoplanetary disks

Authors: Chen-En Wei, Hideko Nomura, Jeong-Eun Lee, Wing-Huen Ip, Catherine Walsh, T. J. Millar

Abstract: The bulk composition of the Earth is dramatically carbon poor compared to that of the interstellar medium, and this phenomenon extends to the asteroid belt. To interpret this carbon deficit problem, the carbonaceous component in grains must have been converted into the gas-phase in the inner regions of protoplanetary disks prior to planetary formation. We examine the effect of carbon grain destruc… ▽ More The bulk composition of the Earth is dramatically carbon poor compared to that of the interstellar medium, and this phenomenon extends to the asteroid belt. To interpret this carbon deficit problem, the carbonaceous component in grains must have been converted into the gas-phase in the inner regions of protoplanetary disks prior to planetary formation. We examine the effect of carbon grain destruction on the chemical structure of disks by calculating the molecular abundances and distributions using a comprehensive chemical reaction network. When carbon grains are destroyed and the elemental abundance of the gas becomes carbon-rich, the abundances of carbon-bearing molecules, such as HCN and carbon-chain molecules, increase dramatically near the midplane, while oxygen-bearing molecules, such as H$_2$O and CO$_2$, are depleted. We compare the results of these model calculations with the solid carbon fraction in the solar system. Although we find a carbon depletion gradient, there are some quantitative discrepancies: the model shows a higher value at the position of asteroid belt and a lower value at the location of the Earth. In addition, using the obtained molecular abundances distributions, coupled with line radiative transfer calculations, we make predictions for ALMA to potentially observe the effect of carbon grain destruction in nearby protoplanetary disks. The results indicate that HCN, H$^{13}$CN and c-C$_3$H$_2$ may be good tracers. △ Less

Submitted 11 December, 2018; v1 submitted 26 November, 2018; originally announced November 2018.

Comments: 31 pages, 21 figures, Accepted for publication in The Astrophysical Journal

Determining the effects of clumping and porosity on the chemistry in a non-uniform AGB outflow

Authors: M. Van de Sande, J. O. Sundqvist, T. J. Millar, D. Keller, W. Homan, A. de Koter, L. Decin, F. De Ceuster

Abstract: (abridged) In the inner regions of AGB outflows, several molecules have been detected with abundances much higher than those predicted from thermodynamic equilibrium (TE) chemical models. The presence of the majority of these species can be explained by shock-induced non-TE chemical models, where shocks caused by the pulsating star take the chemistry out of TE in the inner region. Moreover, a non-… ▽ More (abridged) In the inner regions of AGB outflows, several molecules have been detected with abundances much higher than those predicted from thermodynamic equilibrium (TE) chemical models. The presence of the majority of these species can be explained by shock-induced non-TE chemical models, where shocks caused by the pulsating star take the chemistry out of TE in the inner region. Moreover, a non-uniform density structure has been detected in several AGB outflows. A detailed parameter study on the quantitative effects of a non-homogeneous outflow has so far not been performed. We implement a porosity formalism for treating the increased leakage of light associated with radiation transport through a clumpy, porous medium. The effects from the altered UV radiation field penetration on the chemistry, accounting also for the increased reaction rates of two-body processes in the overdense clumps, are examined. We present a parameter study of the effect of clumping and porosity on the chemistry throughout the outflow. Both the higher density within the clumps and the increased UV radiation field penetration have an important impact on the chemistry, as they both alter the chemical pathways. The increased amount of UV radiation in the inner region leads to photodissociation of parent species, releasing the otherwise deficient elements. We find an increased abundance in the inner region of all species not expected to be present assuming TE chemistry, such as HCN in O-rich outflows, H$_2$O in C-rich outflows, and NH$_3$ in both. Outflows whose clumps have a large overdensity and that are very porous to the interstellar UV radiation field yield abundances comparable to those observed in O- and C-rich outflows for most of the unexpected species investigated. The inner wind abundances of H$_2$O in C-rich outflows and of NH$_3$ in O- and C-rich outflows are however underpredicted. △ Less

Submitted 6 March, 2018; v1 submitted 5 March, 2018; originally announced March 2018.

Comments: 33 pages, 20 figures, 15 tables, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 616, A106 (2018)

Candidate Water Vapor Lines to Locate the $\mathrm{H_2O}$ Snowline through High-dispersion Spectroscopic Observations. III. Submillimeter $\mathrm{H_2}$$^{16}\mathrm{O}$ and $\mathrm{H_2}$$^{18}\mathrm{O}$ Lines

Authors: Shota Notsu, Hideko Nomura, Catherine Walsh, Mitsuhiko Honda, Tomoya Hirota, Eiji Akiyama, T. J. Millar

Abstract: In this paper, we extend the results presented in our former papers (Notsu et al. 2016, 2017) on using ortho-$\mathrm{H_2}$$^{16}\mathrm{O}$ line profiles to constrain the location of the $\mathrm{H_2O}$ snowline in T Tauri and Herbig Ae disks, to include sub-millimeter para-$\mathrm{H_2}$$^{16}\mathrm{O}$ and ortho- and para-$\mathrm{H_2}$$^{18}\mathrm{O}$ lines. Since the number densities of the… ▽ More In this paper, we extend the results presented in our former papers (Notsu et al. 2016, 2017) on using ortho-$\mathrm{H_2}$$^{16}\mathrm{O}$ line profiles to constrain the location of the $\mathrm{H_2O}$ snowline in T Tauri and Herbig Ae disks, to include sub-millimeter para-$\mathrm{H_2}$$^{16}\mathrm{O}$ and ortho- and para-$\mathrm{H_2}$$^{18}\mathrm{O}$ lines. Since the number densities of the ortho- and para-H$_{2}$$^{18}$O molecules are about 560 times smaller than their $^{16}$O analogues, they trace deeper into the disk than the ortho-H$_{2}$$^{16}$O lines (down to $z=0$, i.e., the midplane). Thus these H$_{2}$$^{18}$O lines are potentially better probes of the position of the H$_{2}$O snowline at the disk midplane, depending on the dust optical depth. The values of the Einstein $A$ coefficients of sub-millimeter candidate water lines tend to be lower (typically $<$$10^{-4}$ s$^{-1}$) than infrared candidate water lines (Notsu et al. 2017). Thus in the sub-millimeter candidate water line cases, the local intensity from the outer optically thin region in the disk is around $10^{4}$ times smaller than that in the infrared candidate water line cases. Therefore, in the sub-millimeter lines, especially H$_{2}$$^{18}$O and para-H$_{2}$$^{16}$O lines with relatively lower upper state energies ($\sim$ a few 100K) can also locate the position of the $\mathrm{H_2O}$ snowline. We also investigate the possibility of future observations with ALMA to identify the position of the water snowline. There are several candidate water lines that trace the hot water vapor inside the $\mathrm{H_2O}$ snowline in ALMA Bands $5-10$. △ Less

Submitted 30 May, 2018; v1 submitted 23 January, 2018; originally announced January 2018.

Comments: 24 pages, 10 figures, and 4 tables are contained in this paper. It was accepted for publication in The Astrophysical Journal (ApJ) on January 22th, 2018

Sulphur-bearing molecules in AGB stars I: The occurrence of hydrogen sulfide

Authors: T. Danilovich, M. Van de Sande, E. De Beck, L. Decin, H. Olofsson, S. Ramstedt, T. J. Millar

Abstract: Through a survey of (sub-)millimetre emission lines of various sulphur-bearing molecules, we aim to determine which molecules are the primary carriers of sulphur in different types of AGB stars. In this paper, the first in a series, we investigate the occurrence of H$_2$S in AGB circumstellar envelopes and determine its abundance, where possible. We have surveyed 20 AGB stars with a range of mass-… ▽ More Through a survey of (sub-)millimetre emission lines of various sulphur-bearing molecules, we aim to determine which molecules are the primary carriers of sulphur in different types of AGB stars. In this paper, the first in a series, we investigate the occurrence of H$_2$S in AGB circumstellar envelopes and determine its abundance, where possible. We have surveyed 20 AGB stars with a range of mass-loss rates and of different chemical types using the APEX telescope to search for rotational transition lines of five key sulphur-bearing molecules: CS, SiS, SO, SO$_2$ and H$_2$S. Here we present our results for H$_2$S, including detections, non-detections and detailed radiative transfer modelling of the detected lines. We compare results based on different descriptions of the molecular excitation of H$_2$S and different abundance distributions, including those derived from chemical modelling results. We detected H$_2$S towards five AGB stars, all of which have high mass-loss rates of $\dot{M}\geq 5\times 10^{-6}M_\odot$ yr$^{-1}$ and are oxygen-rich. H$_2$S was not detected towards the carbon or S-type stars that fall in a similar mass-loss range. For the stars in our sample with detections, we find peak o-H$_2$S abundances relative to H$_2$ between $4\times 10^{-7}$ and $2.5\times 10^{-5}$. Overall, we conclude that H$_2$S can play a significant role in oxygen-rich AGB stars with higher mass-loss rates, but is unlikely to play a key role in stars of other chemical types or the lower mass-loss rate oxygen-rich stars. For two sources, V1300 Aql and GX Mon, H$_2$S is most likely the dominant sulphur-bearing molecule in the circumstellar envelope. △ Less

Submitted 19 July, 2017; originally announced July 2017.

Comments: 14 pages, 7 figures, accepted in A&A

Journal ref: A&A 606, A124 (2017)

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

Authors: Shota Notsu, Hideko Nomura, Daiki Ishimoto, Catherine Walsh, Mitsuhiko Honda, Tomoya Hirota, T. J. Millar

Abstract: 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 disks around Herbig Ae stars ($T_{\mathrm{*}}\sim$ 10,000K, $M_{\mathrm{*}}\gtrsim$ 2.5$M_{\bigodot}$), the position of the H$_{2}$O snowline is further from the central star compared with that around cooler, and less mas… ▽ More 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 disks around Herbig Ae stars ($T_{\mathrm{*}}\sim$ 10,000K, $M_{\mathrm{*}}\gtrsim$ 2.5$M_{\bigodot}$), the position of the H$_{2}$O snowline is further from the central star compared with that around cooler, and less massive T Tauri stars. Thus, the H$_{2}$O emission line fluxes from the region within the H$_{2}$O snowline are expected to be stronger. In this paper, we calculate the chemical composition of a Herbig Ae disk using chemical kinetics. Next, we calculate the H$_{2}$O emission line profiles, and investigate the properties of candidate water lines across a wide range of wavelengths (from mid-infrared to sub-millimeter) that can locate the position of the H$_{2}$O snowline. Those line identified have small Einstein $A$ coefficients ($\sim 10^{-6} -10^{-3}$ s$^{-1}$) and relatively high upper state energies ($\sim$ 1000K). The total fluxes tend to increase with decreasing wavelengths. We investigate the possibility of future observations (e.g., ALMA, SPICA/SMI-HRS) to locate the position of the H$_{2}$O snowline. Since the fluxes of those identified lines from Herbig Ae disks are stronger than those from T Tauri disks, the possibility of a successful detection is expected to increase for a Herbig Ae disk. △ Less

Submitted 16 January, 2017; originally announced January 2017.

Comments: 39 pages, 14 figures, and 2 tables are contained in this paper. It was received by The Astrophysical Journal (ApJ) on October 27th, 2016, and was accepted on January 13th, 2017. arXiv admin note: text overlap with arXiv:1606.05828

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

Authors: Catherine Walsh, Attila Juhász, Gwendolyn Meeus, William R. F. Dent, Luke Maud, Yuri Aikawa, Tom J. Millar, Hideko Nomura

Abstract: 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-symmet… ▽ More 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. △ Less

Submitted 10 November, 2016; v1 submitted 7 September, 2016; originally announced September 2016.

Comments: 44 pages, 13 figures, published in ApJ, minor corrections applied to match proofs

Journal ref: 2016, ApJ, 831, 200

First detection of gas-phase methanol in a protoplanetary disk

Authors: Catherine Walsh, Ryan A. Loomis, Karin I. Oberg, Mihkel Kama, Merel L. R. van 't Hoff, Tom J. Millar, Yuri Aikawa, Eric Herbst, Susanna L. Widicus Weaver, Hideko Nomura

Abstract: 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 pe… ▽ More 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. △ Less

Submitted 21 June, 2016; originally announced June 2016.

Comments: 14 pages, 4 figures, 1 table, published online in ApJL on 13th May 2016

Journal ref: ApJL, 823, L10 (2016)

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

Authors: Shota Notsu, Hideko Nomura, Daiki Ishimoto, Catherine Walsh, Mitsuhiko Honda, Tomoya Hirota, T. J. Millar

Abstract: 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 re… ▽ More 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). △ Less

Submitted 18 June, 2016; originally announced June 2016.

Comments: 33 pages, 9 figures, and 2 tables are contained in this paper. It was received by The Astrophysical Journal (ApJ) on December 30th, 2015, and was accepted on June 16th, 2016

Chemistry in AGB stars: successes and challenges

Authors: T. J. Millar

Abstract: Emission and absorption line observations of molecules in late-type stars are a vital component in our understanding of stellar evolution, dust formation and mass loss in these objects. The molecular composition of the gas in the circumstellar envelopes of AGB stars reflects chemical processes in gas whose properties are strong functions of radius with density and temperature varying by more than… ▽ More Emission and absorption line observations of molecules in late-type stars are a vital component in our understanding of stellar evolution, dust formation and mass loss in these objects. The molecular composition of the gas in the circumstellar envelopes of AGB stars reflects chemical processes in gas whose properties are strong functions of radius with density and temperature varying by more than ten and two orders of magnitude, respectively. In addition, the interstellar UV field plays a critical role in determining not only molecular abundances but also their radial distributions. In this article, I shall briefly review some recent successful approaches to describing chemistry in both the inner and outer envelopes and outline areas of challenge for the future. △ Less

Submitted 11 May, 2016; originally announced May 2016.

Comments: To appear in the refereed proceedings of the 11th Pacific Rim Conference on Stellar Astrophysics held in Hong Kong, Dec.14-17, 2015, Sun Kwok (editor). 10 pages, 2 figures

ALMA-resolved salt emission traces the chemical footprint and inner wind morphology of VY CMa

Authors: L. Decin, A. M. S. Richards, T. J. Millar, A. Baudry, E. De Beck, W. Homan, N. Smith, M. Van de Sande, C. Walsh

Abstract: (abreviated) We aim to study the inner-wind structure (R<250 Rstar) of the well-known red supergiant VY CMa. We analyse high spatial resolution (~0".24x0".13) ALMA Science Verification (SV) data in band 7 in which four thermal emission lines of gaseous sodium chloride (NaCl) are present at high signal-to-noise ratio. For the first time, the NaCl emission in the inner wind region of VY CMa is spati… ▽ More (abreviated) We aim to study the inner-wind structure (R<250 Rstar) of the well-known red supergiant VY CMa. We analyse high spatial resolution (~0".24x0".13) ALMA Science Verification (SV) data in band 7 in which four thermal emission lines of gaseous sodium chloride (NaCl) are present at high signal-to-noise ratio. For the first time, the NaCl emission in the inner wind region of VY CMa is spatially resolved. The ALMA observations reveal the contribution of up to four different spatial regions. The NaCl emission pattern is different compared to the dust continuum and TiO2 emission already analysed from the ALMA SV data. The emission can be reconciled with an axisymmetric geometry, where the lower density polar/rotation axis has a position angle of ~50 degrees measured from north to east. However, this picture can not capture the full morphological diversity, and discrete mass ejection events need to be invoked to explain localized higher-density regions. The velocity traced by the gaseous NaCl line profiles is significantly lower than the average wind terminal velocity, and much slower than some of the fastest mass ejections, signalling a wide range of characteristic speeds for the mass loss. Gaseous NaCl is detected far beyond the main dust condensation region. Realising the refractory nature of this metal halide, this hints at a chemical process preventing all NaCl from condensing onto dust grains. We show that in the case of the ratio of the surface binding temperature to the grain temperature being ~50, only some 10% of NaCl remains in gaseous form, while for lower values of this ratio thermal desorption efficiently evaporates NaCl. Photodesorption by stellar photons seems not to be a viable explanation for the detection of gaseous NaCl at 220 Rstar from the central star, and instead, we propose shock-induced sputtering driven by localized mass ejection events as alternative. △ Less

Submitted 10 May, 2016; originally announced May 2016.

Comments: 17 pages, Astronomy & Astrophysics

Journal ref: A&A 592, A76 (2016)

A Gap with a Deficit of Large Grains in the protoplanetary disk around TW Hya

Authors: Takashi Tsukagoshi, Hideko Nomura, Takayuki Muto, Ryohei Kawabe, Daiki Ishimoto, Kazuhiro D. Kanagawa, Satoshi Okuzumi, Shigeru Ida, Catherine Walsh, Tom J. Millar

Abstract: We report $\sim$3 au resolution imaging observations of the protoplanetary disk around TW Hya at 145 and 233 GHz with the Atacama Large Millimeter/Submillimeter Array. Our observations revealed two deep gaps ($\sim$25--50 \%) at 22 and 37~au and shallower gaps (a few \%) at 6, 28, and 44~au, as recently reported by Andrews et al. (2016). The central hole with a radius of $\sim3$~au was also margin… ▽ More We report $\sim$3 au resolution imaging observations of the protoplanetary disk around TW Hya at 145 and 233 GHz with the Atacama Large Millimeter/Submillimeter Array. Our observations revealed two deep gaps ($\sim$25--50 \%) at 22 and 37~au and shallower gaps (a few \%) at 6, 28, and 44~au, as recently reported by Andrews et al. (2016). The central hole with a radius of $\sim3$~au was also marginally resolved. The most remarkable finding is that the spectral index $α(R)$ between bands 4 and 6 peaks at the 22~au gap. The derived power-law index of the dust opacity $β(R)$ is $\sim1.7$ at the 22~au gap and decreases toward the disk center to $\sim0$. The most prominent gap at 22~au could be caused by the gravitational interaction between the disk and an unseen planet with a mass of $\lesssim$1.5 $M_\mathrm{Neptune}$ although other origins may be possible. The planet-induced gap is supported by the fact that $β(R)$ is enhanced at the 22~au gap, indicating a deficit of $\sim$mm-sized grains within the gap due to dust filtration by a planet. △ Less

Submitted 31 August, 2016; v1 submitted 1 May, 2016; originally announced May 2016.

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

Authors: Hideko Nomura, Takashi Tsukagoshi, Ryohei Kawabe, Daiki Ishimoto, Satoshi Okuzumi, Takayuki Muto, Kazuhiro D. Kanagawa, Shigeru Ida, Catherine Walsh, T. J. Millar, Xue-Ning Bai

Abstract: 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 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. △ Less

Submitted 26 February, 2016; v1 submitted 16 December, 2015; originally announced December 2015.

Comments: 6 pages, 5 figures

Journal ref: ApJL, 819, L7 (2016)

The irradiation of water ice by C$^+$ ions in the cosmic environment

Authors: E. J. McBride, T. J. Millar, J. J. Kohanoff

Abstract: We present a first principles molecular dynamics (FPMD) study of the interaction of low energy, positively charged, carbon (C+) projectiles with amorphous solid water clusters at 30 K. Reactions involving the carbon ion at an initial energy of 11 eV and 1.7 eV with 30-molecule clusters have been investigated. Simulations indicate that the neutral isoformyl radical, COH, and carbon monoxide, CO, ar… ▽ More We present a first principles molecular dynamics (FPMD) study of the interaction of low energy, positively charged, carbon (C+) projectiles with amorphous solid water clusters at 30 K. Reactions involving the carbon ion at an initial energy of 11 eV and 1.7 eV with 30-molecule clusters have been investigated. Simulations indicate that the neutral isoformyl radical, COH, and carbon monoxide, CO, are the dominant products of these reactions. All these reactions are accompanied by the transfer of a proton from the reacting water molecule to the ice, where it forms a hydronium ion. We find that COH is formed either via a direct, "knock-out", mechanism following the impact of the C+ projectile upon a water molecule or by creation of a COH_2^+ intermediate. The direct mechanism is more prominent at higher energies. CO is generally produced following the dissociation of COH. More frequent production of the formyl radical, HCO, is observed here than in gas phase calculations. A less commonly occurring product is the dihydroxymethyl, CH(OH)_2, radical. Although a minor result, its existence gives an indication of the increasing chemical complexity which is possible in such heterogeneous environments. △ Less

Submitted 6 August, 2014; originally announced August 2014.

Comments: J Phys Chem A, in press, 24 pages, 6 figures

New detections of HC$_{5}$N toward hot cores associated with 6.7 GHz methanol masers

Authors: Claire-Elise Green, James A. Green, Michael G. Burton, Shinji Horiuchi, Nicholas F. H. Tothill, Andrew J. Walsh, Cormac R. Purcell, James E. J. Lovell, Thomas J. Millar

Abstract: We present new detections of cyanodiacetylene (HC$_{5}$N) toward hot molecular cores, observed with the Tidbinbilla 34 m radio telescope (DSS-34). In a sample of 79 hot molecular cores, HC$_{5}$N was detected towards 35. These results are counter to the expectation that long chain cyanopolyynes, such as HC$_{5}$N, are not typically found in hot molecular cores, unlike their shorter chain counterpa… ▽ More We present new detections of cyanodiacetylene (HC$_{5}$N) toward hot molecular cores, observed with the Tidbinbilla 34 m radio telescope (DSS-34). In a sample of 79 hot molecular cores, HC$_{5}$N was detected towards 35. These results are counter to the expectation that long chain cyanopolyynes, such as HC$_{5}$N, are not typically found in hot molecular cores, unlike their shorter chain counterpart HC$_{3}$N. However it is consistent with recent models which suggest HC$_{5}$N may exist for a limited period during the evolution of hot molecular cores. △ Less

Submitted 7 July, 2014; originally announced July 2014.

Comments: 13 pages, 7 figures. Accepted for publication in MNRAS, July 3 2014

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

Authors: Xiaohu Li, T. J. Millar, Catherine Walsh, Alan N. Heays, Ewine F. van Dishoeck

Abstract: 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 u… ▽ More 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. △ Less

Submitted 28 June, 2014; originally announced June 2014.

Comments: 13 pages, 17 figures, 3 tables. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 568, A111 (2014)

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

Authors: Catherine Walsh, Eric Herbst, Hideko Nomura, T. J. Millar, Susanna Widicus Weaver

Abstract: (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… ▽ More (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. △ Less

Submitted 3 June, 2014; originally announced June 2014.

Comments: Accepted for publication in Faraday Discussion, 168 "Astrochemistry of Dust, Gas, and Ice"

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

Authors: Catherine Walsh, Attila Juhasz, Paola Pinilla, Daniel Harsono, Geoffrey S. Mathews, William R. F. Dent, Michiel R. Hogerheijde, T. Birnstiel, Gwendolyn Meeus, Hideko Nomura, Yuri Aikawa, T. J. Millar, Goran Sandell

Abstract: 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… ▽ More 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. △ Less

Submitted 10 July, 2014; v1 submitted 26 May, 2014; originally announced May 2014.

Comments: Accepted for publication in the Astrophysical Journal Letters

Complex organic molecules in protoplanetary disks

Authors: Catherine Walsh, T. J. Millar, Hideko Nomura, Eric Herbst, Susanna L. Widicus Weaver, Yuri Aikawa, Jake C. Laas, Anton I. Vasyunin

Abstract: (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 protoplane… ▽ More (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. △ Less

Submitted 3 March, 2014; originally announced March 2014.

Comments: 34 pages, 12 figures, 7 tables

Journal ref: Walsh et al. 2014, A&A, 536, A33 (available here: http://www.aanda.org/articles/aa/abs/2014/03/aa22446-13/aa22446-13.html)

Photodetachment as destruction mechanism for CN- and C3N- anions in circumstellar envelopes

Authors: S. S. Kumar, D. Hauser, R. Jindra, T. Best, S. Roucka, W. D. Geppert, T. J. Millar, R. Wester

Abstract: Absolute photodetachment cross sections of two anions of astrophysical importance CN- and C3N- were measured to be (1.18 +- (0.03)_stat (0.17)_sys) * 10^-17 cm^2 and (1.43 +- (0.14)_stat (0.37)_sys) * 10^-17 cm^2 respectively at the ultraviolet wavelength of 266 nm (4.66 eV). These relatively large values of the cross sections imply that photodetachment can play a major role in the destruction mec… ▽ More Absolute photodetachment cross sections of two anions of astrophysical importance CN- and C3N- were measured to be (1.18 +- (0.03)_stat (0.17)_sys) * 10^-17 cm^2 and (1.43 +- (0.14)_stat (0.37)_sys) * 10^-17 cm^2 respectively at the ultraviolet wavelength of 266 nm (4.66 eV). These relatively large values of the cross sections imply that photodetachment can play a major role in the destruction mechanisms of these anions particularly in photon-dominated regions. We have therefore carried out model calculations using the newly measured cross sections to investigate the abundance of these molecular anions in the cirumstellar envelope of the carbon-rich star IRC+10216. The model predicts the relative importance of the various mechanisms of formation and destruction of these species in different regions of the envelope. UV photodetachment was found to be the major destruction mechanism for both CN- and C3N- anions in those regions of the envelope, where they occur in peak abundance. It was also found that photodetachment plays a crucial role in the degradation of these anions throughout the circumstellar envelope. △ Less

Submitted 10 July, 2013; originally announced July 2013.

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

Authors: Catherine Walsh, T. J. Millar, Hideko Nomura

Abstract: 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 T… ▽ More 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). △ Less

Submitted 20 March, 2013; originally announced March 2013.

Comments: 24 pages, 4 figures, 1 table, journal-formatted article available at http://iopscience.iop.org/2041-8205/766/2/L23/

Journal ref: Walsh, Millar, & Nomura 2013, ApJL, 766, L23

The UMIST database for astrochemistry 2012

Authors: D. McElroy, C. Walsh, A. J. Markwick, M. A. Cordiner, K. Smith, T. J. Millar

Abstract: We present the fifth release of the UMIST Database for Astrochemistry (UDfA). The new reaction network contains 6173 gas-phase reactions, involving 467 species, 47 of which are new to this release. We have updated rate coefficients across all reaction types. We have included 1171 new anion reactions and updated and reviewed all photorates. In addition to the usual reaction network, we also now inc… ▽ More We present the fifth release of the UMIST Database for Astrochemistry (UDfA). The new reaction network contains 6173 gas-phase reactions, involving 467 species, 47 of which are new to this release. We have updated rate coefficients across all reaction types. We have included 1171 new anion reactions and updated and reviewed all photorates. In addition to the usual reaction network, we also now include, for download, state-specific deuterated rate coefficients, deuterium exchange reactions and a list of surface binding energies for many neutral species. Where possible, we have referenced the original source of all new and existing data. We have tested the main reaction network using a dark cloud model and a carbon-rich circumstellar envelope model. We present and briefly discuss the results of these models. △ Less

Submitted 27 December, 2012; originally announced December 2012.

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

Authors: Catherine Walsh, Hideko Nomura, T. J. Millar, Yuri Aikawa

Abstract: 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… ▽ More 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. △ Less

Submitted 12 January, 2012; originally announced January 2012.

Comments: 21 pages, 11 figures, 6 tables, accepted for publication in ApJ