Showing 1–34 of 34 results for author: Van de Sande, M

Chemical tracers of a highly eccentric AGB-main sequence star binary

Authors: T. Danilovich, J. Malfait, M. Van de Sande, M. Montargès, P. Kervella, F. De Ceuster, A. Coenegrachts, T. J. Millar, A. M. S. Richards, L. Decin, C. A. Gottlieb, C. Pinte, E. De Beck, D. J. Price, K. T. Wong, J. Bolte, K. M. Menten, A. Baudry, A. de Koter, S. Etoka, D. Gobrecht, M. Gray, F. Herpin, M. Jeste, E. Lagadec , et al. (10 additional authors not shown)

Abstract: Binary interactions have been proposed to explain a variety of circumstellar structures seen around evolved stars, including asymptotic giant branch (AGB) stars and planetary nebulae. Studies resolving the circumstellar envelopes of AGB stars have revealed spirals, discs and bipolar outflows, with shaping attributed to interactions with a companion. For the first time, we have used a combined chem… ▽ More Binary interactions have been proposed to explain a variety of circumstellar structures seen around evolved stars, including asymptotic giant branch (AGB) stars and planetary nebulae. Studies resolving the circumstellar envelopes of AGB stars have revealed spirals, discs and bipolar outflows, with shaping attributed to interactions with a companion. For the first time, we have used a combined chemical and dynamical analysis to reveal a highly eccentric and long-period orbit for W Aquilae, a binary system containing an AGB star and a main sequence companion. Our results are based on anisotropic SiN emission, the first detections of NS and SiC towards an S-type star, and density structures observed in the CO emission. These features are all interpreted as having formed during periastron interactions. Our astrochemistry-based method can yield stringent constraints on the orbital parameters of long-period binaries containing AGB stars, and will be applicable to other systems. △ Less

Submitted 23 July, 2024; originally announced July 2024.

Comments: This version of the article has been accepted for publication, after peer review (when applicable) and is subject to Springer Nature's AM terms of use, but is not the Version of Record and does not reflect post-acceptance improvements, or any corrections. The Version of Record is available online at: https://doi.org/10.1038/s41550-023-02154-y

Journal ref: Nature Astronomy, 8, pages 308-327 (2024)

Hybrid approach predicts a lower binding energy for benzene on water ice

Authors: Victoria H. J. Clark, David M. Benoit, Marie Van de Sande, Catherine Walsh

Abstract: In this paper we provide a highly accurate value for the binding energy of benzene to proton-ordered crystalline water ice (XIh), as a model for interstellar ices. We compare our computed value to the latest experimental data available from temperature programmed desorption (TPD) experiments and find that our binding energy value agrees well with data obtained from binding to either crystalline or… ▽ More In this paper we provide a highly accurate value for the binding energy of benzene to proton-ordered crystalline water ice (XIh), as a model for interstellar ices. We compare our computed value to the latest experimental data available from temperature programmed desorption (TPD) experiments and find that our binding energy value agrees well with data obtained from binding to either crystalline or amorphous ice. Importantly, our new value is lower than that used in most astrochemical networks by about nearly half its value. We explore the impact of this revised binding energy value for both an AGB outflow and a protoplanetary disk. We find that the lower value of the binding energy predicted here compared with values used in the literature (4050 K versus 7587 K) leads to less depletion of gas-phase benzene in an AGB outflow, and leads to a shift outwards in the benzene snowline in the midplane of a protoplanetary disk. Using this new value, the AGB model predicts lower abundances of benzene in the solid phase throughout the outflow. The disk model also predicts a larger reservoir of gas-phase benzene in the inner disk, which is consistent with the recent detections of benzene for the first time in protoplanetary disks with JWST. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Modelling predicts a molecule-rich disk around the AGB star L2 Puppis

Authors: M. Van de Sande, C. Walsh, T. Danilovich, F. De Ceuster, T. Ceulemans

Abstract: The nearby oxygen-rich AGB star L2 Pup hosts a well-studied nearly edge-on disk. To date, disks around AGB stars have not been chemically studied in detail. By combining a parameterisation commonly used for protoplanetary disks and archival ALMA observations, we retrieved an updated density and temperature structure of this disk. This physical model was then used as input to the first chemical mod… ▽ More The nearby oxygen-rich AGB star L2 Pup hosts a well-studied nearly edge-on disk. To date, disks around AGB stars have not been chemically studied in detail. By combining a parameterisation commonly used for protoplanetary disks and archival ALMA observations, we retrieved an updated density and temperature structure of this disk. This physical model was then used as input to the first chemical model of an AGB disk. The model shows that the physical structure of the disk has a large impact on its chemistry, with certain species showing large changes in column density relative to a radial outflow, indicating that chemistry could be used as a tracer of disks that cannot be directly imaged. Despite its oxygen-rich nature, the daughter species formed within the disk are surprisingly carbon-rich. Two chemical regimes can be distinguished: cosmic-ray induced chemistry in the midplane and photochemistry induced by the interstellar radiation field in the outer regions. Certain complex organic molecules are formed in the midplane. This occurs via gas-phase chemistry only, as the disk is too warm for dust-gas chemistry. The photochemistry in the outer regions leads to the efficient formation of (long) carbon-chains. The predictions of the model allow us to tentatively put the disk's age $\lesssim 10^5$ yr. Additional observations are necessary to better constrain the physical structure of L2 Pup's disk and are essential to test the predictions made by the chemical model. Our exploratory work paves the way for a more general study of the chemistry of AGB disks. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: Accepted for publication in MNRAS

MACE: A Machine learning Approach to Chemistry Emulation

Authors: S. Maes, F. De Ceuster, M. Van de Sande, L. Decin

Abstract: The chemistry of an astrophysical environment is closely coupled to its dynamics, the latter often found to be complex. Hence, to properly model these environments a 3D context is necessary. However, solving chemical kinetics within a 3D hydro simulation is computationally infeasible for a even a modest parameter study. In order to develop a feasible 3D hydro-chemical simulation, the classical che… ▽ More The chemistry of an astrophysical environment is closely coupled to its dynamics, the latter often found to be complex. Hence, to properly model these environments a 3D context is necessary. However, solving chemical kinetics within a 3D hydro simulation is computationally infeasible for a even a modest parameter study. In order to develop a feasible 3D hydro-chemical simulation, the classical chemical approach needs to be replaced by a faster alternative. We present mace, a Machine learning Approach to Chemistry Emulation, as a proof-of-concept work on emulating chemistry in a dynamical environment. Using the context of AGB outflows, we have developed an architecture that combines the use of an autoencoder (to reduce the dimensionality of the chemical network) and a set of latent ordinary differential equations (that are solved to perform the temporal evolution of the reduced features). Training this architecture with an integrated scheme makes it possible to successfully reproduce a full chemical pathway in a dynamical environment. mace outperforms its classical analogue on average by a factor 26. Furthermore, its efficient implementation in PyTorch results in a sub-linear scaling with respect to the number of hydrodynamical simulation particles. △ Less

Submitted 6 May, 2024; originally announced May 2024.

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

MELCHIORS: The Mercator Library of High Resolution Stellar Spectroscopy

Authors: P. Royer, T. Merle, K. Dsilva, S. Sekaran, H. Van Winckel, Y. Frémat, M. Van der Swaelmen, S. Gebruers, A. Tkachenko, M. Laverick, M. Dirickx, G. Raskin, H. Hensberge, M. Abdul-Masih, B. Acke, M. L. Alonso, S. Bandhu Mahato, P. G. Beck, N. Behara, S. Bloemen, B. Buysschaert, N. Cox, J. Debosscher, P. De Cat, P. Degroote , et al. (49 additional authors not shown)

Abstract: Over the past decades, libraries of stellar spectra have been used in a large variety of science cases, including as sources of reference spectra for a given object or a given spectral type. Despite the existence of large libraries and the increasing number of projects of large-scale spectral surveys, there is to date only one very high-resolution spectral library offering spectra from a few hundr… ▽ More Over the past decades, libraries of stellar spectra have been used in a large variety of science cases, including as sources of reference spectra for a given object or a given spectral type. Despite the existence of large libraries and the increasing number of projects of large-scale spectral surveys, there is to date only one very high-resolution spectral library offering spectra from a few hundred objects from the southern hemisphere (UVES-POP) . We aim to extend the sample, offering a finer coverage of effective temperatures and surface gravity with a uniform collection of spectra obtained in the northern hemisphere. Between 2010 and 2020, we acquired several thousand echelle spectra of bright stars with the Mercator-HERMES spectrograph located in the Roque de Los Muchachos Observatory in La Palma, whose pipeline offers high-quality data reduction products. We have also developed methods to correct for the instrumental response in order to approach the true shape of the spectral continuum. Additionally, we have devised a normalisation process to provide a homogeneous normalisation of the full spectral range for most of the objects. We present a new spectral library consisting of 3256 spectra covering 2043 stars. It combines high signal-to-noise and high spectral resolution over the entire range of effective temperatures and luminosity classes. The spectra are presented in four versions: raw, corrected from the instrumental response, with and without correction from the atmospheric molecular absorption, and normalised (including the telluric correction). △ Less

Submitted 5 November, 2023; originally announced November 2023.

Comments: 17 pages, 18 figures Preview and access to the library: https://www.royer.se/melchiors.html

Journal ref: A&A 681, A107 (2024)

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

Sensitivity study of chemistry in AGB outflows using chemical kinetics

Authors: S. Maes, M. Van de Sande, T. Danilovich, F. De Ceuster, L. Decin

Abstract: Asymptotic Giant Branch (AGB) stars shed a significant amount of their mass in the form of a stellar wind, creating a vast circumstellar envelope (CSE). Owing to the ideal combination of relatively high densities and cool temperatures, CSEs serve as rich astrochemical laboratories. {While the chemical structure of AGB outflows has been modelled and analysed in detail for specific physical setups,… ▽ More Asymptotic Giant Branch (AGB) stars shed a significant amount of their mass in the form of a stellar wind, creating a vast circumstellar envelope (CSE). Owing to the ideal combination of relatively high densities and cool temperatures, CSEs serve as rich astrochemical laboratories. {While the chemical structure of AGB outflows has been modelled and analysed in detail for specific physical setups, there is a lack of understanding regarding the impact of changes in the physical environment on chemical abundances. A systematic sensitivity study is necessary to comprehend the nuances in the physical parameter space, given the complexity of the chemistry. This is crucial for estimating uncertainties associated with simulations and observations. In this work, we present the first sensitivity study of the impact of varying outflow densities and temperature profiles on the chemistry. With the use of a chemical kinetics model, we report on the uncertainty in abundances, given a specific uncertainty on the physical parameters. }Additionally, we analyse the molecular envelope extent of parent species and compare our findings to observational studies. Mapping the impact of differences in physical parameters throughout the CSE on the chemistry is a strong aid to observational studies. △ Less

Submitted 19 April, 2023; v1 submitted 12 April, 2023; originally announced April 2023.

Comments: 21 pages, 25 figures

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)

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

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)

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)

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

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

Chemical modelling of dust-gas chemistry within AGB outflows -- II. Effect of the dust-grain size distribution

Authors: M. Van de Sande, C. Walsh, T. Danilovich

Abstract: AGB stars are, together with supernovae, the main contributors of stellar dust to the interstellar medium (ISM). Dust grains formed by AGB stars are thought to be large. However, as dust nucleation and growth within their outflows are still not understood, the dust-grain size distribution (GSD) is unknown. This is an important uncertainty regarding our knowledge of the chemical and physical histor… ▽ More AGB stars are, together with supernovae, the main contributors of stellar dust to the interstellar medium (ISM). Dust grains formed by AGB stars are thought to be large. However, as dust nucleation and growth within their outflows are still not understood, the dust-grain size distribution (GSD) is unknown. This is an important uncertainty regarding our knowledge of the chemical and physical history of interstellar dust, as AGB dust forms $\sim$ 70% of the starting point of its evolution. We expand on our chemical kinetics model, which uniquely includes a comprehensive dust-gas chemistry. The GSD is now allowed to deviate from the commonly assumed canonical Mathis et al. (1977) distribution. We find that the specific GSD can significantly influence the dust-gas chemistry within the outflow. Our results show that the level of depletion of gas-phase species depends on the average grain surface area of the GSD. Gas-phase abundance profiles and their possible depletions can be retrieved from observations of molecular emission lines when using a range of transitions. Due to degeneracies within the prescription of GSD, specific parameters cannot be retrieved, only (a lower limit to) the average grain surface area. Nonetheless, this can discriminate between dust composed of predominantly large or small grains. We show that when combined with other observables such as the spectral energy distribution and polarised light, depletion levels from molecular gas-phase abundance profiles can constrain the elusive GSD of the dust delivered to the ISM by AGB outflows. △ Less

Submitted 4 May, 2020; originally announced May 2020.

Comments: 17 pages, 10 figures, 6 tables. Accepted for publication in MNRAS

An ALMA view of SO and SO$_2$ around oxygen-rich AGB stars

Authors: T. Danilovich, A. M. S. Richards, L. Decin, M. Van de Sande, C. A. Gottlieb

Abstract: We present and analyse SO and SO$_2$, recently observed with high angular resolution and sensitivity in a spectral line survey with ALMA, for two oxygen-rich AGB stars: the low mass-loss rate R Dor and high mass-loss rate IK Tau. We analyse 8 lines of SO detected towards both stars, 78 lines of SO$_2$detected towards R Dor and 52 lines of SO$_2$ detected towards IK Tau. We detect several lines of… ▽ More We present and analyse SO and SO$_2$, recently observed with high angular resolution and sensitivity in a spectral line survey with ALMA, for two oxygen-rich AGB stars: the low mass-loss rate R Dor and high mass-loss rate IK Tau. We analyse 8 lines of SO detected towards both stars, 78 lines of SO$_2$detected towards R Dor and 52 lines of SO$_2$ detected towards IK Tau. We detect several lines of $^{34}$SO, $^{33}$SO and $^{34}$SO$_2$ towards both stars, and tentatively S$^{18}$O towards R Dor, and hence derive isotopic ratios for these species. The spatially resolved observations show us that the two sulphur oxides are co-located towards R Dor and trace out the same wind structures in the circumstellar envelope (CSE). Much of the emission is well reproduced with a Gaussian abundance distribution spatially centred on the star. Emission from the higher energy levels of SO and SO$_2$ towards R Dor provide evidence in support of a rotating inner region of gas identified in earlier work. The new observations allow us to refine the abundance distribution of SO in IK Tau derived from prior observations with single antennas, and confirm the distribution is shell-like with the peak in the fractional abundance not centred on the star. The confirmation of different types of SO abundance distributions will help fine-tune chemical models and allows for an additional method to discriminate between low and high mass-loss rates for oxygen-rich AGB stars. △ Less

Submitted 9 March, 2020; originally announced March 2020.

Comments: Accepted for publication in MNRAS

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

Dust-gas chemistry in AGB outflows: Chemical modelling of dust-gas chemistry within AGB outflows I. Effect on the gas-phase chemistry

Authors: M. Van de Sande, C. Walsh, T. P. Mangan, L. Decin

Abstract: Chemical modelling of AGB outflows is typically focused on either non-thermodynamic equilibrium chemistry in the inner region or photon-driven chemistry in the outer region. We include, for the first time, a comprehensive dust-gas chemistry in our AGB outflow chemical kinetics model, including both dust-gas interactions and grain-surface chemistry. The dust is assumed to have formed in the inner r… ▽ More Chemical modelling of AGB outflows is typically focused on either non-thermodynamic equilibrium chemistry in the inner region or photon-driven chemistry in the outer region. We include, for the first time, a comprehensive dust-gas chemistry in our AGB outflow chemical kinetics model, including both dust-gas interactions and grain-surface chemistry. The dust is assumed to have formed in the inner region, and follows an interstellar-like dust-size distribution. Using radiative transfer modelling, we obtain dust temperature profiles for different dust types in an O-rich and a C-rich outflow. We calculate a grid of models, sampling different outflow densities, drift velocities between the dust and gas, and dust types. Dust-gas chemistry can significantly affect the gas-phase composition, depleting parent and daughter species and increasing the abundance of certain daughter species via grain-surface formation followed by desorption/sputtering. Its influence depends on four factors: outflow density, dust temperature, initial composition, and drift velocity. The largest effects are for higher density outflows with cold dust and O-rich parent species, as these species generally have a larger binding energy. At drift velocities larger than $\sim 10$ km s$^{-1}$, ice mantles undergo sputtering; however, they are not fully destroyed. Models with dust-gas chemistry can better reproduce the observed depletion of species in O-rich outflows. When including colder dust in the C-rich outflows and adjusting the binding energy of CS, the depletion in C-rich outflows is also better reproduced. To best interpret high-resolution molecular line observations from AGB outflows, dust-gas interactions are needed in chemical kinetics models. △ Less

Submitted 24 September, 2019; v1 submitted 23 September, 2019; originally announced September 2019.

Comments: 20 pages, 14 figures, accepted for for publication in the Monthly Notices of the Royal Astronomical Society

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

An ALMA view of CS and SiS around oxygen-rich AGB stars

Authors: T. Danilovich, A. M. S. Richards, A. I. Karakas, M. Van de Sande, L. Decin, F. De Ceuster

Abstract: We aim to determine the distributions of molecular SiS and CS in the circumstellar envelopes of oxygen-rich asymptotic giant branch stars and how these distributions differ between stars that lose mass at different rates. In this study we analyse ALMA observations of SiS and CS emission lines for three oxygen-rich galactic AGB stars: IK Tau, with a moderately high mass-loss rate of… ▽ More We aim to determine the distributions of molecular SiS and CS in the circumstellar envelopes of oxygen-rich asymptotic giant branch stars and how these distributions differ between stars that lose mass at different rates. In this study we analyse ALMA observations of SiS and CS emission lines for three oxygen-rich galactic AGB stars: IK Tau, with a moderately high mass-loss rate of $5\times10^{-6}$M$_\odot$ yr$^{-1}$, and W Hya and R Dor with low mass loss rates of $\sim1\times10^{-7}$M$_\odot$ yr$^{-1}$. These molecules are usually more abundant in carbon stars but the high sensitivity of ALMA allows us to detect their faint emission in the low mass-loss rate AGB stars. The high spatial resolution of ALMA also allows us to precisely determine the spatial distribution of these molecules in the circumstellar envelopes. We run radiative transfer models to calculate the molecular abundances and abundance distributions for each star. We find a spread of peak SiS abundances with $\sim10^{-8}$ for R Dor, $\sim10^{-7}$ for W Hya, and $\sim3\times10^{-6}$ for IK Tau relative to H$_2$. We find lower peak CS abundances of $\sim7\times10^{-9}$ for R Dor, $\sim7\times10^{-8}$ for W Hya and $\sim4\times10^{-7}$ for IK Tau, with some stratifications in the abundance distributions. For IK Tau we also calculate abundances for the detected isotopologues: C$^{34}$S, $^{29}$SiS, $^{30}$SiS, Si$^{33}$S, Si$^{34}$S, $^{29}$Si$^{34}$S, and $^{30}$Si$^{34}$S. Overall the isotopic ratios we derive for IK Tau suggest a lower metallicity than solar. △ Less

Submitted 31 December, 2018; originally announced January 2019.

Comments: 16 pages

ALMA detection of a tentative nearly edge-on rotating disk around the nearby AGB star R~Doradus

Authors: Ward Homan, Taissa Danilovich, Leen Decin, Alex de Koter, Joseph Nuth, Marie Van de Sande

Abstract: A spectral scan of the circumstellar environment of the asymptotic giant branch (AGB) star R~Doradus was taken with ALMA in cycle 2 at frequencies between 335 and 362 GHz and with a spatial resolution of $\sim$150 milliarcseconds. Many molecular lines show a spatial offset between the blue and red shifted emission in the innermost regions of the wind. The position-velocity diagrams of this feature… ▽ More A spectral scan of the circumstellar environment of the asymptotic giant branch (AGB) star R~Doradus was taken with ALMA in cycle 2 at frequencies between 335 and 362 GHz and with a spatial resolution of $\sim$150 milliarcseconds. Many molecular lines show a spatial offset between the blue and red shifted emission in the innermost regions of the wind. The position-velocity diagrams of this feature, in combination with previous SPHERE data and theoretical work point towards the presence of a compact differentially rotating disk, orientated nearly edge-on. We model the $^{\rm 28}$SiO ($v=1,~J=8\to7$) emission with a disk model. We estimate the disk mass and angular momentum to be $3 \times 10^{-6}$ Solar masses and $5 \times 10^{40}\ {\rm m^2 kg/s}$. The latter presents an `angular momentum problem' that may be solved by assuming that the disk is the result of wind-companion interactions with a companion of at least 2.5 earth masses, located at 6 AU, the tentatively determined location of the disk's inner rim. An isolated clump of emission is detected to the south-east with a velocity that is high compared to the previously determined terminal velocity of the wind. Its position and mean velocity suggest that it may be associated with a companion planet, located at the disk's inner rim. △ Less

Submitted 16 March, 2018; originally announced March 2018.

Comments: 11 pages, 7 Figures

Journal ref: A&A 614, A113 (2018)

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)

Chemical content of the circumstellar envelope of the oxygen-rich AGB star R Dor: Non-LTE abundance analysis of CO, SiO, and HCN

Authors: M. Van de Sande, L. Decin, R. Lombaert, T. Khouri, A. de Koter, F. Wyrowski, R. De Nutte, W. Homan

Abstract: (abridged) Our aim is to determine the radial abundance profile of SiO and HCN throughout the stellar outflow of R Dor, an oxygen-rich AGB star with a low mass-loss rate. We have analysed molecular transitions of CO, SiO, and HCN measured with the APEX telescope and all three instruments on the Herschel Space Observatory, together with literature data. Photometric data and the infrared spectrum me… ▽ More (abridged) Our aim is to determine the radial abundance profile of SiO and HCN throughout the stellar outflow of R Dor, an oxygen-rich AGB star with a low mass-loss rate. We have analysed molecular transitions of CO, SiO, and HCN measured with the APEX telescope and all three instruments on the Herschel Space Observatory, together with literature data. Photometric data and the infrared spectrum measured by ISO-SWS were used to constrain the dust component of the outflow. Using both continuum and line radiative transfer methods, a physical envelope model of both gas and dust was established. We have performed an analysis of the SiO and HCN molecular transitions in order to calculate their abundances. We have obtained an envelope model that describes the dust and the gas in the outflow, and determined the abundance of SiO and HCN throughout the region of the outflow probed by our molecular data. For SiO, we find that the initial abundance lies between $5.5 \times 10^{-5}$ and $6.0 \times 10^{-5}$ w.r.t. H$_2$. The abundance profile is constant up to $60\ \pm 10\ R_*$, after which it declines following a Gaussian profile with an $e$-folding radius of $3.5 \pm 0.5 \times 10^{13}$ cm. For HCN, we find an initial abundance of $5.0 \times 10^{-7}$ w.r.t. H$_2$. The Gaussian profile that describes the decline starts at the stellar surface and has an $e$-folding radius $r_e$ of $1.85 \pm 0.05 \times 10^{15}$ cm. We cannot to unambiguously identify the mechanism by which SiO is destroyed at $60\ \pm 10\ R_*$. The initial abundances found are larger than previously determined (except for one previous study on SiO), which might be due to the inclusion of higher-$J$ transitions. The difference in abundance for SiO and HCN compared to high mass-loss rate Mira star IK Tau might be due to different pulsation characteristics of the central star and/or a difference in dust condensation physics. △ Less

Submitted 30 August, 2017; originally announced August 2017.

Comments: 16 pages, 13 figures, 7 tables, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 609, A63 (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)

Tracing the phase transition of Al-bearing species from molecules to dust in AGB winds. Constraining the presence of gas-phase (Al2O3)n clusters

Authors: L. Decin, A. M. S. Richards, L. B. F. M. Waters, T. Danilovich, D. Gobrecht, T. Khouri, W. Homan, J. Bakker, M. Van de Sande, J. A. Nuth, E. De Beck

Abstract: (abbreviated) We aim to constrain the dust formation histories in the winds of oxygen-rich AGB stars. We have obtained ALMA observations with a spatial resolution of 120x150 mas tracing the dust formation region of a low mass-loss rate and a high mass-loss rate AGB star, R Dor and IK Tau. Emission line profiles of AlO, AlOH and AlCl are detected and are used to derive a lower limit of atomic alumi… ▽ More (abbreviated) We aim to constrain the dust formation histories in the winds of oxygen-rich AGB stars. We have obtained ALMA observations with a spatial resolution of 120x150 mas tracing the dust formation region of a low mass-loss rate and a high mass-loss rate AGB star, R Dor and IK Tau. Emission line profiles of AlO, AlOH and AlCl are detected and are used to derive a lower limit of atomic aluminium incorporated in molecules. We show that the gas-phase aluminium chemistry is completely different in both stars, with a remarkable difference in the AlO and AlOH abundance stratification. The amount of aluminium locked up in these 3 molecules is small, <=1.1e-7, for both stars, i.e. only <=2% of the total aluminium budget. This leaves ample of room for aluminium to be incorporated in grains. A fundamental result is that AlO and AlOH, being the direct precursors of alumina grains, are detected well beyond the onset of the dust condensation proving that the aluminium oxide condensation cycle is not fully efficient. The ALMA observations allow us to quantitatively assess the current generation of theoretical dynamical-chemical models for AGB winds. We discuss how the current proposed scenario of aluminium dust condensation for low mass-loss rate AGB stars at a distance of ~1.5 Rstar, in particular for the stars R Dor and W Hya, poses a challenge if one wishes to explain both the dust spectral features in the spectral energy distribution (SED), in interferometric data, and in polarized light signal. In particular, the estimated grain temperature of Al2O3 is too high for the grains to retain their amorphous structure. We propose that large gas-phase (Al2O3)n-clusters (n>34) can be the potential agents of the broad 11 micron feature in the SED and in the interferometric data and we explain how these large clusters can be formed. △ Less

Submitted 23 April, 2017; v1 submitted 18 April, 2017; originally announced April 2017.

Comments: 21 pages, submitted to Astronomy & Astrophysics (A&A), figures at reduced resolution in arxiv paper; full resolution in Astronomy & Astrophysics paper

Journal ref: A&A 608, A55 (2017)

Search for exoplanets around pulsating stars of A--F type in Kepler Short Cadence data and the case of KIC 8197761

Authors: Paulina Sowicka, Gerald Handler, Bartłomiej Dębski, David Jones, Marie Van de Sande, Péter I. Pápics

Abstract: We searched for extrasolar planets around pulsating stars by examining $\textit{Kepler}$ data for transit-like events hidden in the intrinsic variability. All Short Cadence observations for targets with 6000 K $< T_{\rm eff} <$ 8500 K were visually inspected for transit-like events following the removal of pulsational signals by sinusoidal fits. Clear transit-like events were detected in KIC 56133… ▽ More We searched for extrasolar planets around pulsating stars by examining $\textit{Kepler}$ data for transit-like events hidden in the intrinsic variability. All Short Cadence observations for targets with 6000 K $< T_{\rm eff} <$ 8500 K were visually inspected for transit-like events following the removal of pulsational signals by sinusoidal fits. Clear transit-like events were detected in KIC 5613330 and KIC 8197761. KIC 5613330 is a confirmed exoplanet host (Kepler-635b), where the transit period determined here is consistent with the literature value. KIC 8197761 is a $γ$ Doradus - $δ$ Scuti star exhibiting eclipses/transits occurring every 9.8686667(27) d, having durations of 8.37 h, and causing brightness drops $\frac{ΔF}{F} = 0.00629(29)$. The star's pulsation spectrum contains several mode doublets and triplets, identified as $l = 1$, with a mean spacing of 0.001659(15) c/d, implying an internal rotation period of $301\pm3$ d. Trials to calculate the size of the light travel time effect (LTTE) from the pulsations to constrain the companion's mass ended inconclusive. Finding planets around $γ$ Doradus stars from the pulsational LTTE, therefore, is concluded to be unrealistic. Spectroscopic monitoring of KIC 8197761 revealed sinusoidal radial velocity variations with a semi-amplitude of $19.75 \pm 0.32$ km/s, while individual spectra present rotational broadening consistent with $v \sin i = 9\pm1$ km/s. This suggests that the stellar surface rotation is synchronized with the orbit, whereas the stellar core rotates $\sim$30 times slower. Combining the observed radial velocity variability with the transit photometry, constrains the companion's mass to be $\approx 0.28$ M$_{\odot}$, ruling out an exoplanet hypothesis. △ Less

Submitted 16 February, 2017; originally announced February 2017.

Comments: 12 pages, 9 figures, accepted for publication in MNRAS

Journal ref: MNRAS (2017) 467 (4): 4663-4673

Signatures of internal rotation discovered in the Kepler data of five slowly pulsating B stars

Authors: P. I. Pápics, A. Tkachenko, T. Van Reeth, C. Aerts, E. Moravveji, M. Van de Sande, K. De Smedt, S. Bloemen, J. Southworth, J. Debosscher, E. Niemczura, J. F. Gameiro

Abstract: Massive stars are important building blocks of the Universe, and their stellar structure and evolution models are fundamental cornerstones of various fields in modern astrophysics. The precision of these models is limited by our lack of understanding of various internal mixing processes that significantly influence the lifetime of these objects (e.g. core overshoot, chemical mixing, or the interna… ▽ More Massive stars are important building blocks of the Universe, and their stellar structure and evolution models are fundamental cornerstones of various fields in modern astrophysics. The precision of these models is limited by our lack of understanding of various internal mixing processes that significantly influence the lifetime of these objects (e.g. core overshoot, chemical mixing, or the internal differential rotation). Our goal is to calibrate models by extending the sample of available seismic studies of slowly pulsating B (SPB) stars, providing input for theoretical modelling efforts that will deliver precise constraints on the parameters describing the internal mixing processes in these objects. We used spectral synthesis and disentangling techniques to derive fundamental parameters and to determine precise orbital parameters from high-resolution spectra. We employed custom masks to construct light curves from the Kepler pixel data and used standard time-series analysis tools to construct a set of significant frequencies for each target. These sets were first filtered from combination frequencies, and then screened for period spacing patterns. We detect gravity mode period series of modes of the same degree l with consecutive radial order n in four new and one revisited SPB star. These predominantly prograde dipole series (covering ~10-40 radial orders) are influenced by rotation and carry signatures of chemical mixing processes. Our spectroscopic analysis, in addition to placing each object inside the SPB instability strip and identifying KIC 4930889 as an SB2 binary, reveals that KIC 11971405 is a fast rotator that shows very weak Be signatures. Together with the observed photometric outbursts this illustrates that this Be star is a fast rotating SPB star. We hypothesise that the outbursts might be connected to its very densely compressed oscillation spectrum. △ Less

Submitted 21 November, 2016; originally announced November 2016.

Comments: Accepted for publication in Astronomy and Astrophysics; 23 pages containing 29 figures and 14 tables

Journal ref: A&A 598, A74 (2017)

Nucleosynthesis in AGB stars traced by oxygen isotopic ratios. I - Determining the stellar initial mass by means of the $^{17}$O/$^{18}$O ratio

Authors: R. De Nutte, L. Decin, H. Olofsson, R. Lombaert, A. de Koter, A. Karakas, S. Milam, S. Ramstedt, R. J. Stancliffe, W. Homan, M. Van de Sande

Abstract: The aim of this paper is to investigate the $^{17}$O/$^{18}$O ratio for a sample of AGB stars, containing M-, S- and C-type stars. These ratios are evaluated in relation to fundamental stellar evolution parameters: the stellar initial mass and pulsation period. Circumstellar $^{13}$C$^{16}$O, $^{12}$C$^{17}$O and $^{12}$C$^{18}$O line observations were obtained for a sample of nine stars with vari… ▽ More The aim of this paper is to investigate the $^{17}$O/$^{18}$O ratio for a sample of AGB stars, containing M-, S- and C-type stars. These ratios are evaluated in relation to fundamental stellar evolution parameters: the stellar initial mass and pulsation period. Circumstellar $^{13}$C$^{16}$O, $^{12}$C$^{17}$O and $^{12}$C$^{18}$O line observations were obtained for a sample of nine stars with various single-dish long-wavelength facilities. Line intensity ratios are shown to relate directly to the surface $^{17}$O/$^{18}$O abundance ratio. Stellar evolution models predict the $^{17}$O/$^{18}$O ratio to be a sensitive function of initial mass and to remain constant throughout the entire TP-AGB phase for stars initially less massive than 5\,$M_{\odot}$. This makes the measured ratio a probe of the initial stellar mass. Observed $^{17}$O/$^{18}$O ratios are found to be well in the range predicted by stellar evolution models that do not consider convective overshooting. From this, accurate initial mass estimates are calculated for seven sources. For the remaining two sources two mass solutions result, though with a larger probability that the low-mass solution is the correct one. Finally, hints at a possible separation between M/S- and C-type stars when comparing the $^{17}$O/$^{18}$O ratio to the stellar pulsation period are presented. △ Less

Submitted 14 December, 2016; v1 submitted 23 June, 2016; originally announced June 2016.

Comments: Accepted for publication in Astronomy and Astrophysics. Revised: implemented comments from referee

Journal ref: A&A 600, A71 (2017)

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)

The rms-flux relation in accreting white dwarfs: another nova-like variable and the first dwarf nova

Authors: M. Van de Sande, S. Scaringi, C. Knigge

Abstract: We report on the detection of the linear rms-flux relation in two accreting white dwarf binary systems: V1504 Cyg and KIC 8751494. The rms-flux relation relates the absolute root-mean-square (rms) variability of the light curve to its mean flux. The light curves analysed were obtained with the Kepler satellite at a 58.8 s cadence. The rms-flux relation was previously detected in only one other cat… ▽ More We report on the detection of the linear rms-flux relation in two accreting white dwarf binary systems: V1504 Cyg and KIC 8751494. The rms-flux relation relates the absolute root-mean-square (rms) variability of the light curve to its mean flux. The light curves analysed were obtained with the Kepler satellite at a 58.8 s cadence. The rms-flux relation was previously detected in only one other cataclysmic variable, MV Lyr. This result reenforces the ubiquity of the linear rms-flux relation as a characteristic property of accretion-induced variability, since it has been observed in several black hole binaries, neutron star binaries and active galactic nuclei. Moreover, its detection in V1504 Cyg is the first time the rms-flux relation has been detected in a dwarf nova-type CV during quiescence. This result, together with previous studies, hence points towards a common physical origin of accretion-induced variability, independent of the size, mass, or type of the central accreting compact object. △ Less

Submitted 21 January, 2015; originally announced January 2015.

Comments: 8 pages, 6 figures. Accepted to MNRAS