-
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.
-
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.
-
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.
-
Bayesian model reconstruction based on spectral line observations
Authors:
Frederik De Ceuster,
Thomas Ceulemans,
Leen Decin,
Taïssa Danilovich,
Jeremy Yates
Abstract:
Spectral line observations encode a wealth of information. A key challenge, therefore, lies in the interpretation of these observations in terms of models to derive the physical and chemical properties of the astronomical environments from which they arise. In this paper, we present pomme: an open-source Python package that allows users to retrieve 1D or 3D models of physical properties, such as c…
▽ More
Spectral line observations encode a wealth of information. A key challenge, therefore, lies in the interpretation of these observations in terms of models to derive the physical and chemical properties of the astronomical environments from which they arise. In this paper, we present pomme: an open-source Python package that allows users to retrieve 1D or 3D models of physical properties, such as chemical abundance, velocity, and temperature distributions of (optically thin) astrophysical media, based on spectral line observations. We discuss how prior knowledge, for instance, in the form of a steady-state hydrodynamics model, can be used to guide the retrieval process, and demonstrate our methods both on synthetic and real observations of cool stellar winds.
△ Less
Submitted 28 February, 2024;
originally announced February 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.
-
Self-consistent dust and non-LTE line radiative transfer with SKIRT
Authors:
Kosei Matsumoto,
Peter Camps,
Maarten Baes,
Frederik De Ceuster,
Keiichi Wada,
Takao Nakagawa,
Kentaro Nagamine
Abstract:
We introduce Monte Carlo-based non-LTE line radiative transfer calculations in the 3D dust radiative transfer code SKIRT, which was originally set up as a dust radiative transfer code. By doing so, we develop a generic and powerful 3D radiative transfer code that can self-consistently generate spectra with molecular and atomic lines against the underlying continuum. We test the accuracy of the non…
▽ More
We introduce Monte Carlo-based non-LTE line radiative transfer calculations in the 3D dust radiative transfer code SKIRT, which was originally set up as a dust radiative transfer code. By doing so, we develop a generic and powerful 3D radiative transfer code that can self-consistently generate spectra with molecular and atomic lines against the underlying continuum. We test the accuracy of the non-LTE line radiative transfer module in the extended SKIRT code using standard benchmarks. We find excellent agreement between the SKIRT results, the published benchmark results, and results obtained using the ray-tracing non-LTE line radiative transfer code MAGRITTE, which validates our implementation. We apply the extended SKIRT code on a 3D hydrodynamic simulation of a dusty AGN torus model and generate multi-wavelength images with CO rotational-line spectra against the underlying dust continuum. We find that the low-J CO emission traces the geometrically thick molecular torus, whereas the higher-J CO lines originate from the gas with high kinetic temperature located in the innermost regions of the torus. Comparing the calculations with and without dust radiative transfer, we find that higher-J CO lines are slightly attenuated by the surrounding cold dust when seen edge-on. This shows that atomic and molecular lines can experience attenuation, an effect that is particularly important for transitions at mid- and near-infrared wavelengths. Therefore, our self-consistent dust and non-LTE line radiative transfer calculations can help interpret the observational data from Herschel, ALMA, and JWST.
△ Less
Submitted 5 September, 2023;
originally announced September 2023.
-
3D simulations of AGB stellar winds -- II. Ray-tracer implementation and impact of radiation on the outflow morphology
Authors:
Mats Esseldeurs,
Lionel Siess,
Frederik De Ceuster,
Ward Homan,
Jolien Malfait,
Silke Maes,
Thomas Konings,
Thomas Ceulemans,
Leen Decin
Abstract:
Stars with an initial mass below ~ 8 Msun evolve through the asymptotic giant branch (AGB) phase, during which they develop strong stellar winds. Recent observations have revealed significant morphological complexities in their outflows, most likely caused by a companion. We study the impact of the radiation force on such companion-perturbed AGB outflows. We present the implementation of a ray tra…
▽ More
Stars with an initial mass below ~ 8 Msun evolve through the asymptotic giant branch (AGB) phase, during which they develop strong stellar winds. Recent observations have revealed significant morphological complexities in their outflows, most likely caused by a companion. We study the impact of the radiation force on such companion-perturbed AGB outflows. We present the implementation of a ray tracer for radiative transfer in smoothed particle hydrodynamics (SPH) and compared four different descriptions of radiative transfer: the free-wind, the geometrical, the Lucy, and the attenuation approximation. For both low and high mass-loss rates, the velocity profile of the outflow is modified when going from the free-wind to the geometrical approximation, also resulting in a different morphology. In the case of a low mass-loss rate, the effect of the Lucy and attenuation approximation is negligible due to the low densities but morphological differences appear in the high mass-loss rate regime. By comparing the radiative equilibrium temperature and radiation force to full 3D radiative transfer, we show that the Lucy approximation works best. Although, close to the companion, artificial heating occurs and it fails to simulate the shadow cast by the companion. The attenuation approximation produces a lower equilibrium temperature and weaker radiation force, but it produces the shadow cast by the companion. From the predictions of the 3D radiative transfer, we also conclude that a radially directed radiation force is a reasonable assumption. The radiation force thus plays a critical role in dust-driven AGB winds, impacting the velocity profile and morphological structures. For low mass-loss rates, the geometrical approximation suffices, while high mass-loss rates require a more rigorous method, where the Lucy approximation provides the most accurate results although not accounting for all effects.
△ Less
Submitted 19 April, 2023;
originally announced April 2023.
-
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.
-
The unusual 3D distribution of NaCl around the AGB star IK Tau
Authors:
A. Coenegrachts,
T. Danilovich,
F. De Ceuster,
L. Decin
Abstract:
NaCl is a diatomic molecule with a large dipole moment, which allows for its detection even at relatively small abundances. It has been detected towards several evolved stars, among which is the AGB star IK Tau, around which it is distributed in several clumps that lie off-center from the star. We aim to study the three-dimensional distribution of NaCl around the AGB star IK Tau, and to obtain the…
▽ More
NaCl is a diatomic molecule with a large dipole moment, which allows for its detection even at relatively small abundances. It has been detected towards several evolved stars, among which is the AGB star IK Tau, around which it is distributed in several clumps that lie off-center from the star. We aim to study the three-dimensional distribution of NaCl around the AGB star IK Tau, and to obtain the abundance of NaCl relative to H$_2$ for each of the clumps. First, a new value for the maximum expansion velocity is determined. The observed ALMA channel maps are then deprojected to create a three-dimensional model of the distribution of NaCl. This model is then used as input for the radiative transfer modelling code magritte, which is used to obtain the NaCl abundances of each of the clumps by comparing the observations with the results of the magritte simulations. Additionally, the rotational temperature of the clumps is determined using population diagrams. We derive an updated value for the maximum expansion velocity of IK Tau $\upsilon_\mathrm{exp}$ = 28.4 km/s. A spiral-like shape can be discerned in our three-dimensional distribution model of the NaCl. This spiral lies more or less in the plane of the sky. The distribution is also flatter in the line-of-sight direction than in the plane of the sky. We find clump abundances between $9 \times 10^{-9}$ and $5 \times 10^{-8}$ relative to H$_2$, where the relative abundance is typically lower for clumps closer to the star. For the first time, we used deprojection to understand the three-dimensional environment of an AGB star and calculated the fractional abundance of NaCl in clumps surrounding the star.
△ Less
Submitted 23 August, 2023; v1 submitted 13 February, 2023;
originally announced February 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.
-
Radiative Transfer as a Bayesian Linear Regression problem
Authors:
Frederik De Ceuster,
Thomas Ceulemans,
Jon Cockayne,
Leen Decin,
Jeremy Yates
Abstract:
Electromagnetic radiation plays a crucial role in various physical and chemical processes. Hence, almost all astrophysical simulations require some form of radiative transfer model. Despite many innovations in radiative transfer algorithms and their implementation, realistic radiative transfer models remain very computationally expensive, such that one often has to resort to approximate descriptio…
▽ More
Electromagnetic radiation plays a crucial role in various physical and chemical processes. Hence, almost all astrophysical simulations require some form of radiative transfer model. Despite many innovations in radiative transfer algorithms and their implementation, realistic radiative transfer models remain very computationally expensive, such that one often has to resort to approximate descriptions. The complexity of these models makes it difficult to assess the validity of any approximation and to quantify uncertainties on the model results. This impedes scientific rigour, in particular, when comparing models to observations, or when using their results as input for other models. We present a probabilistic numerical approach to address these issues by treating radiative transfer as a Bayesian linear regression problem. This allows us to model uncertainties on the input and output of the model with the variances of the associated probability distributions. Furthermore, this approach naturally allows us to create reduced-order radiative transfer models with a quantifiable accuracy. These are approximate solutions to exact radiative transfer models, in contrast to the exact solutions to approximate models that are often used. As a first demonstration, we derive a probabilistic version of the method of characteristics, a commonly-used technique to solve radiative transfer problems.
△ Less
Submitted 22 November, 2022;
originally announced November 2022.
-
3D hydrodynamical survey of the impact of a companion on the morphology and dynamics of AGB outflow
Authors:
Jolien Malfait,
Silke Maes,
Ward Homan,
Jan Bolte,
Lionel Siess,
Frederik De Ceuster,
Leen Decin
Abstract:
With the use of high-resolution ALMA observations, complex structures that resemble those observed in post-AGB stars and planetary nebulae are detected in the circumstellar envelopes of low-mass evolved stars. These deviations from spherical symmetry are believed to be caused primarily by the interaction with a companion star or planet. With the use of three-dimensional hydrodynamic simulations, w…
▽ More
With the use of high-resolution ALMA observations, complex structures that resemble those observed in post-AGB stars and planetary nebulae are detected in the circumstellar envelopes of low-mass evolved stars. These deviations from spherical symmetry are believed to be caused primarily by the interaction with a companion star or planet. With the use of three-dimensional hydrodynamic simulations, we study the impact of a binary companion on the wind morphology and dynamics of an AGB outflow. We classifiy the wind structures and morphology that form in these simulations with the use of a classification parameter, constructed with characteristic parameters of the binary configuration. Finally we conclude that the companion alters the wind expansion velocity through the slingshot mechanism, if it is massive enough.
△ Less
Submitted 28 July, 2022;
originally announced July 2022.
-
Route towards complete 3D hydro-chemical simulations of companion-perturbed AGB outflows
Authors:
Silke Maes,
Lionel Siess,
Ward Homan,
Jolien Malfait,
Frederik De Ceuster,
Thomas Ceulemans,
Dion Donne,
Mats Esseldeurs,
Leen Decin
Abstract:
Low- and intermediate mass stars experience a significant mass loss during the last phases of their evolution, which obscures them in a vast, dusty envelope. Although it has long been thought this envelope is generally spherically symmetric in shape, recent high-resolution observations find that most of these stars exhibit complex and asymmetrical morphologies, most likely resulting from binary in…
▽ More
Low- and intermediate mass stars experience a significant mass loss during the last phases of their evolution, which obscures them in a vast, dusty envelope. Although it has long been thought this envelope is generally spherically symmetric in shape, recent high-resolution observations find that most of these stars exhibit complex and asymmetrical morphologies, most likely resulting from binary interaction. In order to improve our understanding about these systems, theoretical studies are needed in the form of numerical simulations. Currently, a handful of simulations exist, albeit they mainly focus on the hydrodynamics of the outflow. Hence, we here present the pathway to more detailed and accurate modelling of companion-perturbed outflows with Phantom, by discussing the missing but crucial physical and chemical processes. With these state-of-the-art simulations we aim to make a direct comparison with observations to unveil the true identity on the embedded systems.
△ Less
Submitted 24 June, 2022;
originally announced June 2022.
-
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.
-
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.
-
SPH modelling of wind-companion interactions in eccentric AGB binary systems
Authors:
J. Malfait,
W. Homan,
S. Maes,
J. Bolte,
L. Siess,
F. De Ceuster,
L. Decin
Abstract:
The late evolutionary stages of low- and intermediate-mass stars are characterised by mass loss through a dust-driven stellar wind. Recent observations reveal complex structures within these winds, that are believed to be formed primarily via interaction with a companion. How these complexities arise, and which structures are formed in which type of systems, is still poorly understood. Particularl…
▽ More
The late evolutionary stages of low- and intermediate-mass stars are characterised by mass loss through a dust-driven stellar wind. Recent observations reveal complex structures within these winds, that are believed to be formed primarily via interaction with a companion. How these complexities arise, and which structures are formed in which type of systems, is still poorly understood. Particularly, there is a lack of studies investigating the structure formation in eccentric systems. We aim to improve our understanding of the wind morphology of eccentric AGB binary systems by investigating the mechanism responsible for the different small-scale structures and global morphologies that arise in a polytropic wind with different velocities. Using the smoothed particle hydrodynamics (SPH) code Phantom, we generate nine different high-resolution, 3D simulations of an AGB star with a solar-mass companion with various wind velocity and eccentricity combinations. The models assume a polytropic gas, with no additional cooling. We conclude that for models with a high wind velocity, the short interaction with the companion results in a regular spiral morphology, that is flattened. In the case of a lower wind velocity, the stronger interaction results in the formation of a high-energy region and bow-shock structure that can shape the wind into an irregular morphology if instabilities arise. High-eccentricity models show a complex, phase-dependent interaction leading to wind structures that are irregular in three dimensions. However, the significant interaction with the companion compresses matter into an equatorial density enhancement, irrespective of eccentricity.
△ Less
Submitted 2 July, 2021;
originally announced July 2021.
-
SPH modelling of companion-perturbed AGB outflows including a new morphology classification scheme
Authors:
S. Maes,
W. Homan,
J. Malfait,
L. Siess,
J. Bolte,
F. De Ceuster,
L. Decin
Abstract:
Asymptotic giant branch (AGB) stars are known to lose a significant amount of mass by a stellar wind, which controls the remainder of their stellar lifetime. High angular-resolution observations show that the winds of these cool stars typically exhibit mid- to small-scale density perturbations such as spirals and arcs, believed to be caused by the gravitational interaction with a (sub-)stellar com…
▽ More
Asymptotic giant branch (AGB) stars are known to lose a significant amount of mass by a stellar wind, which controls the remainder of their stellar lifetime. High angular-resolution observations show that the winds of these cool stars typically exhibit mid- to small-scale density perturbations such as spirals and arcs, believed to be caused by the gravitational interaction with a (sub-)stellar companion. We aim to explore the effects of the wind-companion interaction on the 3D density and velocity distribution of the wind, as a function of three key parameters: wind velocity, binary separation and companion mass. For the first time, we compare the impact on the outflow of a planetary companion to that of a stellar companion. We intend to devise a morphology classification scheme based on a singular parameter. With our grid of models we cover the prominent morphology changes in a companion-perturbed AGB outflow: slow winds with a close, massive binary companion show a more complex morphology. Additionally, we prove that massive planets are able to significantly impact the density structure of an AGB wind. We find that the interaction with a companion affects the terminal velocity of the wind, which can be explained by the gravitational slingshot mechanism. We distinguish between two types of wind focussing to the orbital plane resulting from distinct mechanisms: global flattening of the outflow as a result of the AGB star's orbital motion and the formation of an EDE as a consequence of the companion's gravitational pull. We investigate different morphology classification schemes and uncover that the ratio of the gravitational potential energy density of the companion to the kinetic energy density of the AGB outflow yields a robust classification parameter for the models presented in this paper.
△ Less
Submitted 1 July, 2021;
originally announced July 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.
-
Magritte, a modern software library for 3D radiative transfer: II. Adaptive ray-tracing, mesh construction and reduction
Authors:
Frederik De Ceuster,
Jan Bolte,
Ward Homan,
Silke Maes,
Jolien Malfait,
Leen Decin,
Jeremy Yates,
Peter Boyle,
James Hetherington
Abstract:
Radiative transfer is a notoriously difficult and computationally demanding problem. Yet, it is an indispensable ingredient in nearly all astrophysical and cosmological simulations. Choosing an appropriate discretization scheme is a crucial part of the simulation, since it not only determines the direct memory cost of the model but also largely determines the computational cost and the achievable…
▽ More
Radiative transfer is a notoriously difficult and computationally demanding problem. Yet, it is an indispensable ingredient in nearly all astrophysical and cosmological simulations. Choosing an appropriate discretization scheme is a crucial part of the simulation, since it not only determines the direct memory cost of the model but also largely determines the computational cost and the achievable accuracy. In this paper, we show how an appropriate choice of directional discretization scheme as well as spatial model mesh can help alleviate the computational cost, while largely retaining the accuracy. First, we discuss the adaptive ray-tracing scheme implemented in our 3D radiative transfer library Magritte, that adapts the rays to the spatial mesh and uses a hierarchical directional discretization based on HEALPix. Second, we demonstrate how the free and open-source software library Gmsh can be used to generate high quality meshes that can be easily tailored for Magritte. In particular, we show how the local element size distribution of the mesh can be used to optimise the sampling of both analytically and numerically defined models. Furthermore, we show that when using the output of hydrodynamics simulations as input for a radiative transfer simulation, the number of elements in the input model can often be reduced by an order of magnitude, without significant loss of accuracy in the radiation field. We demonstrate this for two models based on a hierarchical octree mesh resulting from adaptive mesh refinement (AMR), as well as two models based on smoothed-particle hydrodynamics (SPH) data.
△ Less
Submitted 30 November, 2020;
originally announced November 2020.
-
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.
-
(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.
-
Magritte, a modern software library for 3D radiative transfer: I. Non-LTE atomic and molecular line modelling
Authors:
Frederik De Ceuster,
Ward Homan,
Jeremy Yates,
Leen Decin,
Peter Boyle,
James Hetherington
Abstract:
Radiative transfer is a key component in almost all astrophysical and cosmological simulations. We present Magritte: a modern open-source software library for 3D radiative transfer. It uses a deterministic ray-tracer and formal solver, i.e. it computes the radiation field by tracing rays through the model and solving the radiative transfer equation in its second-order form along a fixed set of ray…
▽ More
Radiative transfer is a key component in almost all astrophysical and cosmological simulations. We present Magritte: a modern open-source software library for 3D radiative transfer. It uses a deterministic ray-tracer and formal solver, i.e. it computes the radiation field by tracing rays through the model and solving the radiative transfer equation in its second-order form along a fixed set of rays originating from each point. Magritte can handle structured and unstructured input meshes, as well as smoothed-particle hydrodynamics (SPH) particle data. In this first paper, we describe the numerical implementation, semi-analytic tests and cross-code benchmarks for the non-LTE line radiative transfer module of Magritte. This module uses the radiative transfer solver to self-consistently determine the populations of the quantised energy levels of atoms and molecules using an accelerated Lambda iteration (ALI) scheme. We compare Magritte with the established radiative transfer solvers Ratran (1D) and Lime (3D) on the van Zadelhoff benchmark and present a first application to a simple Keplerian disc model. Comparing with Lime, we conclude that Magritte produces more accurate and more precise results, especially at high optical depth, and that it is faster.
△ Less
Submitted 18 December, 2019;
originally announced December 2019.
-
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.
-
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.
