Showing 1–35 of 35 results for author: Homan, W

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)

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

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

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

Submitted 4 May, 2023; originally announced May 2023.

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

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

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.

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

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

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

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

Submitted 5 January, 2023; originally announced January 2023.

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

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

3D simulations of AGB stellar winds -- I. Steady winds and dust formation

Authors: L. Siess, W. Homan, S. Toupin, D. J. Price

Abstract: Aims. We present the implementation of the treatment of particle ejection and dust nucleation in the smoothed particle hydrodynamics (SPH) code phantom. These developments represent the first step toward a more complete modeling of dust-driven winds emanating from AGB stars. Methods. The AGB outflow is modeled by injecting the SPH particles from a spherical inner boundary. This boundary is a serie… ▽ More Aims. We present the implementation of the treatment of particle ejection and dust nucleation in the smoothed particle hydrodynamics (SPH) code phantom. These developments represent the first step toward a more complete modeling of dust-driven winds emanating from AGB stars. Methods. The AGB outflow is modeled by injecting the SPH particles from a spherical inner boundary. This boundary is a series of concentric shells, with the AGB star at its center, and the particles are positioned on these shells on the vertices of an isocahedron geodesic surface. The outermost shell is ejected with a predefined radial velocity, and subsequent lower shells replenish the ejected ones, all rotated randomly to improve the isotropy of the outflow. The physical properties of the particles on these shells are set by solving the 1D analytic steady wind equations. The formation of dust is calculated starting from a compact chemical network for carbon-rich material, which creates the building blocks of the solid-state particles. Subsequently, the theory of the moments is used to obtain dust growth rates, without requiring knowledge on the grain size distribution. Results. We tested our implementation against a series of 1D reference solutions. We demonstrate that our method is able to reproduce Parker-type wind solutions. For the trans-sonic solution, small oscillations are present in the vicinity of the sonic point, but these do not impact the trans-sonic passage or terminal wind velocity. Supersonic solutions always compare nicely with 1D analytic profiles. We also tested our implementation of dust using two formalisms: an analytic prescription for the opacity devised by Bowen and the full treatment of carbon-dust formation. Both simulations reproduce the 1D analytic solution displaying the expected additional acceleration when the gas temperature falls below the condensation temperature. △ Less

Submitted 29 August, 2022; originally announced August 2022.

Comments: 16 pages, 12 figures

Journal ref: A&A 667, A75 (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.

Comments: Proceedings IAU Symposium No. 366, 2021, The origin of outflows in evolved stars

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.

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

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)

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.

Comments: 23 pages, 22 figures

Journal ref: A&A 652, A51 (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.

Comments: 19 pages, 11 figures

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

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.

Comments: 12 pages, published in MNRAS

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

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

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

Submitted 12 October, 2020; originally announced October 2020.

Comments: 19 pages, 20 figures

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

(Sub)stellar companions shape the winds of evolved stars

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

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

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

Comments: 19 pages main journal, 97 pages Supplementary Information

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

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

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

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

Submitted 19 August, 2020; originally announced August 2020.

Comments: 24 pages, 25 Figures

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

Wind morphology around cool evolved stars in binaries: the case of slowly accelerating oxygen-rich outflows

Authors: I. El Mellah, J. Bolte, L. Decin, W. Homan, R. Keppens

Abstract: The late stellar evolutionary phases of low and intermediate-mass stars are strongly constrained by their mass-loss rates. The wind surrounding cool evolved stars frequently shows non-spherical features, thought to be due to an unseen companion orbiting the donor star. We study the morphology of the circumbinary envelope, in particular around oxygen-rich asymptotic giant branch (AGB) stars. We run… ▽ More The late stellar evolutionary phases of low and intermediate-mass stars are strongly constrained by their mass-loss rates. The wind surrounding cool evolved stars frequently shows non-spherical features, thought to be due to an unseen companion orbiting the donor star. We study the morphology of the circumbinary envelope, in particular around oxygen-rich asymptotic giant branch (AGB) stars. We run a grid of 70 3D hydrodynamics simulations of a progressively accelerating wind propagating in the Roche potential formed by a mass-loosing evolved star in orbit with a main sequence companion. We resolve the flow structure both in the immediate vicinity of the secondary, where bow shocks, outflows and wind-captured disks form, and up to 40 orbital separations, where spiral arms, arcs and equatorial density enhancements develop. When the companion is deeply engulfed in the wind, the lower terminal wind speeds and more progressive wind acceleration around oxygen-rich AGB stars make them more prone than carbon-rich AGB stars to display more disturbed outflows, a disk-like structure around the companion and a wind concentrated in the orbital plane. In these configurations, a large fraction of the wind is captured by the companion which leads to a significant shrinking of the orbit over the mass-loss timescale, if the donor star is at least a few times more massive than its companion. Provided the companion has a mass of at least a tenth of the mass of the donor star, it can compress the wind in the orbital plane up to large distances. Our grid of models covers a wide scope of configurations function of the dust chemical content, the terminal wind speed relative to the orbital speed, the extension of the dust condensation region around the cool evolved star and the mass ratio. It provides a frame of reference to interpret high-resolution maps of the outflows surrounding cool evolved stars. △ Less

Submitted 13 January, 2020; originally announced January 2020.

Journal ref: A&A 637, A91 (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.

Comments: 16 pages, accepted in MNRAS

NOEMA maps the CO $J = 2-1$ environment of the red supergiant $μ$ Cep

Authors: M. Montargès, W. Homan, D. Keller, N. Clementel, S. Shetye, L. Decin, G. M. Harper, P. Royer, J. M. Winters, T. Le Bertre, A. M. S. Richards

Abstract: Red supergiant stars are surrounded by a gaseous and dusty circumstellar environment created by their mass loss which spreads heavy elements into the interstellar medium. The structure and the dynamics of this envelope are crucial to understand the processes driving the red supergiant mass loss and the shaping of the pre-supernova ejecta. We have observed the emission from the CO $J = 2-1$ line fr… ▽ More Red supergiant stars are surrounded by a gaseous and dusty circumstellar environment created by their mass loss which spreads heavy elements into the interstellar medium. The structure and the dynamics of this envelope are crucial to understand the processes driving the red supergiant mass loss and the shaping of the pre-supernova ejecta. We have observed the emission from the CO $J = 2-1$ line from the red supergiant star $μ$~Cep with the NOEMA interferometer. In the line the synthesized beam was $0.92 \times 0.72$~arcsec ($590 \times 462$~au at 641~pc). The continuum map shows only the unresolved contribution of the free-free emission of the star chromosphere. The continuum-subtracted channel maps reveal a very inhomogeneous and clumpy circumstellar environment. In particular, we detected a bright CO clump, as bright as the central source in the line, at 1.80~arcsec south-west from the star, in the blue channel maps. After a deprojection of the radial velocity assuming two different constant wind velocities, the observations were modelled using the 3D radiative transfer code \textsc{lime} to derive the characteristics of the different structures. We determine that the gaseous clumps observed around $μ$~Cep are responsible for a mass loss rate of $(4.9 \pm 1.0) \times 10^{-7}~{\rm M}_\odot\,{\rm yr}^{-1}$, in addition to a spatially unresolved wind component with an estimated mass-loss rate of $2.0 \times 10^{-6}~{\rm M}_\odot\,{\rm yr}^{-1}$. Therefore, the clumps have a significant role in $μ$~Cep's mass loss ($\ge 25 \%$). We cannot exclude that the unresolved central outflow may be made of smaller unresolved clumps. △ Less

Submitted 29 April, 2019; v1 submitted 17 March, 2019; originally announced March 2019.

Comments: 15 pages, 4 tables, 9 figures. 2nd version : one co-author removed and acknowledgement updated (consistent with erratum https://doi.org/10.1093/mnras/stz1006)

Journal ref: MNRAS 485, 2417 (2019)

Reduction of the maximum mass-loss rate of OH/IR stars due to unnoticed binary interaction

Authors: L. Decin, W. Homan, T. Danilovich, A. de Koter, D. Engels, L. B. F. M. Waters, S. Muller, C. Gielen, D. A. García-Hernández, R. J. Stancliffe, M. Vande Sande, G. Molenberghs, F. Kerschbaum, A. A. Zijlstra, I. El Mellah

Abstract: In 1981, the idea of a superwind that ends the life of cool giant stars was proposed. Extreme OH/IR-stars develop superwinds with the highest mass-loss rates known so far, up to a few 10^(-4) Msun/yr, informing our understanding of the maximum mass-loss rate achieved during the Asymptotic Giant Branch (AGB) phase. A condundrum arises whereby the observationally determined duration of the superwind… ▽ More In 1981, the idea of a superwind that ends the life of cool giant stars was proposed. Extreme OH/IR-stars develop superwinds with the highest mass-loss rates known so far, up to a few 10^(-4) Msun/yr, informing our understanding of the maximum mass-loss rate achieved during the Asymptotic Giant Branch (AGB) phase. A condundrum arises whereby the observationally determined duration of the superwind phase is too short for these stars to become white dwarfs. Here, we report on the detection of spiral structures around two cornerstone extreme OH/IR-stars, OH26.5+0.6 and OH30.1-0.7, identifying them as wide binary systems. Hydrodynamical simulations show that the companion's gravitational attraction creates an equatorial density enhancement mimicking a short extreme superwind phase, thereby solving the decades-old conundrum. This discovery restricts the maximum mass-loss rate of AGB stars around the single-scattering radiation-pressure limit of a few 10^(-5) Msun/yr. This brings about crucial implications for nucleosynthetic yields, planet survival, and the wind-driving mechanism. △ Less

Submitted 25 February, 2019; originally announced February 2019.

Comments: Publication date: 25 February 2019 at 16:00h GMT

Journal ref: Nature Astronomy, 2019

An unusual face-on spiral in the wind of the M-type AGB star EP Aquarii

Authors: Ward Homan, Anita Richards, Leen Decin, Alex de Koter, Pierre Kervella

Abstract: High-resolution interferometric observations of the circumstellar environments of AGB stars show a variety of morphologies. Guided by the unusual carbon monoxide line profile of the AGB star EP Aquarii, we have observed its circumstellar environment with ALMA band 6 in cycle 4. We describe the morphological complexity of the CO, SiO and SO2 molecular emission. The CO emission exhibits the characte… ▽ More High-resolution interferometric observations of the circumstellar environments of AGB stars show a variety of morphologies. Guided by the unusual carbon monoxide line profile of the AGB star EP Aquarii, we have observed its circumstellar environment with ALMA band 6 in cycle 4. We describe the morphological complexity of the CO, SiO and SO2 molecular emission. The CO emission exhibits the characteristics of a bi-conical wind with a bright nearly face-on spiral feature around the systemic velocity. This is the first convincing detection of a spiral morphology in an O-rich wind. Based on the offsets of the centers of the two bi-conical wind hemispheres, we deduce the position angle of the inclination axis to be ~150 deg measured counterclockwise from north. Based on the velocity width of the spiral signature we estimate the inclination angle of the system to be between 4 deg and 18 deg. The central emission zone exhibits a morphology that resembles simulations modelling the spiral-inducing wind-Roche-lobe-overflow mechanism. Though the spiral may be a (companion-induced) density enhancement in the stellar outflow, the extremely narrow width of the spiral signature in velocity space suggests that it may be a hydrodynamical perturbation in a face-on differentially rotating disk. The SiO emission does not show the spiral, but exhibits a local emission void approximately 0.5'' west of the continuum brightness peak. We hypothesise that this may be a local environment caused by the presence of a (stellar) companion with a mass of at most 0.1 Msol, based on its non-detection in the continuum... △ Less

Submitted 3 May, 2018; v1 submitted 16 April, 2018; originally announced April 2018.

Comments: 20 pages, 16 Figures

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

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)

Simplified models of circumstellar morphologies for interpreting high-resolution data. Analytical approach to the equatorial density enhancement

Authors: Ward Homan, Jels Boulangier, Leen Decin, Alex de Koter

Abstract: Equatorial density enhancements (EDEs) are a very common astronomical phenomenon. Studies of the circumstellar environments (CSE) of young stellar objects and of evolved stars have shown that these objects often possess these features. These are believed to originate from different mechanisms, ranging from binary interactions to the gravitational collapse of interstellar material. Quantifying the… ▽ More Equatorial density enhancements (EDEs) are a very common astronomical phenomenon. Studies of the circumstellar environments (CSE) of young stellar objects and of evolved stars have shown that these objects often possess these features. These are believed to originate from different mechanisms, ranging from binary interactions to the gravitational collapse of interstellar material. Quantifying the effect of the presence of this type of EDE on the observables is essential for a correct interpretation of high-resolution data. We seek to investigate the manifestation in the observables of a circumstellar EDE, to assess which properties can be constrained, and to provide an intuitive bedrock on which to compare and interpret upcoming high-resolution data (e.g. \emph{ALMA} data) using 3D models. We develop a simplified analytical parametrised description of a 3D EDE, with a possible substructure such as warps, gaps, and spiral instabilities. In addition, different velocity fields (Keplerian, radial, super-Keplerian, sub-Keplerian and rigid rotation) are considered. The effect of a bipolar outflow is also investigated. The geometrical models are fed into the 3D radiative transfer code {\tt LIME}, that produces 3D intensity maps throughout velocity space. We investigate the spectral signature of the $J$=3$-$2 up to $J$=7$-$6 rotational transitions of CO in the models, as well as the spatial aspect of this emission by means of channel maps, wide-slit position-velocity (PV) diagrams, stereograms, and spectral lines. Additionally, we discuss methods of constraining the geometry of the EDE, the inclination, the mass-contrast between the EDE and the bipolar outflow, and the global velocity field. Finally, we simulated \emph{ALMA} observations... △ Less

Submitted 14 March, 2018; originally announced March 2018.

Journal ref: Astronomy & Astrophysics, 12/2016, Volume 596, id.A91, 29 pp

ALMA observations of the nearby AGB star L$_{\rm 2}$ Puppis. II. Gas disk properties derived from $^{\rm 12}$CO and $^{\rm13}$CO $J=$3$-$2 emission

Authors: Ward Homan, Anita Richards, Leen Decin, Pierre Kervella, Alex de Koter, Iain McDonald, Keiichi Ohnaka

Abstract: The circumstellar environment of the AGB star L$_{\rm 2}$ Puppis was observed with ALMA in cycle 3, with a resolution of $15 \times 18 \rm\ mas$. The molecular emission shows a differentially rotating disk, inclined to a nearly edge-on position. In the first paper in this series (paper I) the molecular emission was analysed to accurately deduce the motion of the gas in the equatorial regions of th… ▽ More The circumstellar environment of the AGB star L$_{\rm 2}$ Puppis was observed with ALMA in cycle 3, with a resolution of $15 \times 18 \rm\ mas$. The molecular emission shows a differentially rotating disk, inclined to a nearly edge-on position. In the first paper in this series (paper I) the molecular emission was analysed to accurately deduce the motion of the gas in the equatorial regions of the disk. In this work we model the optically thick $^{\rm 12}$CO $J=$3$-$2 and the optically thin $^{\rm 13}$CO $J=$3$-$2 rotational transition to constrain the physical conditions in the disk. To realise this effort we make use of the 3D NLTE radiative transfer code {\tt LIME}. The temperature structure and velocity structure show a high degree of complexity, both radially and vertically. The radial H$_{\rm 2}$ density profile in the disk plane is characterised by a power law with a slope of $-3.1$. We find a $^{\rm 12}$CO over $^{\rm 13}$CO abundance ratio of 10 inside the disk. Finally, estimations of the angular momentum in the disk surpass the expected available angular momentum of the star, strongly supporting the indirect detection of a compact binary companion reported in paper I. We estimate the mass of the companion to be around 1 Jupiter mass. △ Less

Submitted 14 March, 2018; originally announced March 2018.

Comments: 21 pages; 19 figures

Journal ref: Astronomy & Astrophysics, 05/2017, Volume 601, id.A5, 19 pp

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)

ALMA spectral line and imaging survey of a low and a high mass-loss rate AGB star between 335 and 362 GHz

Authors: L. Decin, A. M. S. Richards, T. Danilovich, W. Homan, J. A. Nuth

Abstract: A spectral line and imaging survey of the low mass-loss rate AGB star R Dor (Mdot ~ 1e-7 Msun/yr) and the high mass-loss rate AGB star IK Tau (Mdot ~5e-6 Msun/yr) was made with ALMA between 335 and 362 GHz at a spatial resolution of ~150 mas, corresponding to the locus of the main dust formation region of both targets. Some 200 spectral features from 15 molecules (and their isotopologues) were obs… ▽ More A spectral line and imaging survey of the low mass-loss rate AGB star R Dor (Mdot ~ 1e-7 Msun/yr) and the high mass-loss rate AGB star IK Tau (Mdot ~5e-6 Msun/yr) was made with ALMA between 335 and 362 GHz at a spatial resolution of ~150 mas, corresponding to the locus of the main dust formation region of both targets. Some 200 spectral features from 15 molecules (and their isotopologues) were observed, including rotational lines in both the ground and vibrationally excited states. Detected species include the gaseous precursors of dust grains such as SiO, AlO, AlOH, TiO, and TiO2. We present a spectral atlas for both stars and the parameters of all detected spectral features. A clear dichotomy for the sulphur chemistry is seen: while CS, SiS, SO, and SO2 are abundantly present in IK Tau, only SO and SO2 are detected in R Dor. Also other species such as NaCl, NS, AlO, and AlOH display a completely different behaviour. From some selected species, the minor isotopologues can be used to assess the isotopic ratios. The channel maps of many species prove that both large and small-scale inhomogeneities persist in the inner wind of both stars in the form of blobs, arcs, and/or a disk. The high sensitivity of ALMA allows us to spot the impact of these correlated density structures in the spectral line profiles. The spectral lines often display a half width at zero intensity much larger than expected from the terminal velocity, v_inf, previously derived for both objects (36 km/s versus v_inf ~17.7 km/s for IK Tau and 23 km/s versus v_inf ~5.5 km/s for R Dor). Both a more complex 3D morphology and a more forceful wind acceleration of the (underlying) isotropic wind can explain this trend. The formation of fractal grains in the region beyond ~400 mas can potentially account for the latter scenario. From the continuum map, we deduce a dust mass of ~3.7e-7 Msun for IK Tau and ~2e-8 Msun for R Dor. △ Less

Submitted 28 January, 2018; originally announced January 2018.

Comments: 39 pages, Astronomy & Astrophysics, in press

Journal ref: A&A 615, A28 (2018)

The close circumstellar environment of Betelgeuse - V. Rotation velocity and molecular envelope properties from ALMA

Authors: Pierre Kervella, Leen Decin, Anita M. S. Richards, Graham M. Harper, Iain McDonald, Eamon O'Gorman, Miguel Montargès, Ward Homan, Keiichi Ohnaka

Abstract: We observed Betelgeuse using ALMA's extended configuration in band 7 (f~340 GHz, λ~0.88 mm), resulting in a very high angular resolution of 18 mas. Using a solid body rotation model of the 28SiO(v=2,J=8-7) line emission, we show that the supergiant is rotating with a projected equatorial velocity of v_eq sin i = 5.47 +/- 0.25 km/s at the equivalent continuum angular radius R_star = 29.50 +/- 0.14… ▽ More We observed Betelgeuse using ALMA's extended configuration in band 7 (f~340 GHz, λ~0.88 mm), resulting in a very high angular resolution of 18 mas. Using a solid body rotation model of the 28SiO(v=2,J=8-7) line emission, we show that the supergiant is rotating with a projected equatorial velocity of v_eq sin i = 5.47 +/- 0.25 km/s at the equivalent continuum angular radius R_star = 29.50 +/- 0.14 mas. This corresponds to an angular rotation velocity of ω sin i = (5.6 +/- 1.3) x 10^(-9) rad/s. The position angle of its north pole is PA = 48.0 +/- 3.5°. The rotation period of Betelgeuse is estimated to P/sin i = 36 +/- 8 years. The combination of our velocity measurement with previous observations in the ultraviolet shows that the chromosphere is co-rotating with the star up to a radius of ~10 au (45 mas or 1.5x the ALMA continuum radius). The coincidence of the position angle of the polar axis of Betelgeuse with that of the major ALMA continuum hot spot, a molecular plume, and a partial dust shell (from previous observations) suggests that focused mass loss is currently taking place in the polar region of the star. We propose that this hot spot corresponds to the location of a particularly strong "rogue" convection cell, which emits a focused molecular plume that subsequently condenses into dust at a few stellar radii. Rogue convection cells therefore appear to be an important factor shaping the anisotropic mass loss of red supergiants. △ Less

Submitted 28 November, 2017; v1 submitted 21 November, 2017; originally announced November 2017.

Comments: 18 pages, 19 figures, accepted for publication in A&A, update of bibliography and acknowledgements

Journal ref: A&A 609, A67 (2018)

Radiation-pressure-driven sub-Keplerian rotation of the disc around the AGB star L2 Pup

Authors: Thomas J. Haworth, Richard A. Booth, Ward Homan, Leen Decin, Cathie J. Clarke, Subhanjoy Mohanty

Abstract: We study the sub-Keplerian rotation and dust content of the circumstellar material around the asymptotic giant branch (AGB) star L$_2$ Puppis. We find that the thermal pressure gradient alone cannot explain the observed rotation profile. We find that there is a family of possible dust populations for which radiation pressure can drive the observed sub-Keplerian rotation. This set of solutions is f… ▽ More We study the sub-Keplerian rotation and dust content of the circumstellar material around the asymptotic giant branch (AGB) star L$_2$ Puppis. We find that the thermal pressure gradient alone cannot explain the observed rotation profile. We find that there is a family of possible dust populations for which radiation pressure can drive the observed sub-Keplerian rotation. This set of solutions is further constrained by the spectral energy distribution (SED) of the system, and we find that a dust-to-gas mass ratio of $\sim10^{-3}$ and a maximum grain size that decreases radially outwards can satisfy both the rotation curve and SED. These dust populations are dynamically tightly coupled to the gas azimuthally. However grains larger than $\sim0.5μ$m are driven outward radially by radiation pressure at velocities $\sim5$kms$^{-1}$, which implies a dust replenishment rate of $\sim3\times10^{-9}$M$_\odot$yr$^{-1}$. This replenishment rate is consistent with observational estimates to within uncertainties. Coupling between the radial motion of the dust and gas is weak and hence the gas does not share in this rapid outward motion. Overall we conclude that radiation pressure is a capable and necessary mechanism to explain the observed rotation profile of L$_2$ Pup, and offers other additional constraints on the dust properties. △ Less

Submitted 8 September, 2017; originally announced September 2017.

Comments: Resubmitted to MNRAS after minor revisions

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)

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)

ALMA observations of the nearby AGB star L2 Puppis - I. Mass of the central star and detection of a candidate planet

Authors: P. Kervella, W. Homan, A. M. S. Richards, L. Decin, I. McDonald, M. Montargès, K. Ohnaka

Abstract: Six billion years from now, while evolving on the asymptotic giant branch (AGB), the Sun will metamorphose from a red giant into a beautiful planetary nebula. This spectacular evolution will impact the Solar System planets, but observational confirmations of the predictions of evolution models are still elusive as no planet orbiting an AGB star has yet been discovered. The nearby AGB red giant L2… ▽ More Six billion years from now, while evolving on the asymptotic giant branch (AGB), the Sun will metamorphose from a red giant into a beautiful planetary nebula. This spectacular evolution will impact the Solar System planets, but observational confirmations of the predictions of evolution models are still elusive as no planet orbiting an AGB star has yet been discovered. The nearby AGB red giant L2 Puppis (d = 64 pc) is surrounded by an almost edge-on circumstellar dust disk. We report new observations with ALMA at very high angular resolution (18 x 15 mas) in band 7 (f ~ 350 GHz) that allow us to resolve the velocity profile of the molecular disk. We establish that the gas velocity profile is Keplerian within the central cavity of the dust disk, allowing us to derive the mass of the central star L2 Pup A, mA = 0.659 +/- 0.011 +/- 0.041 Msun (+/- 6.6%). From evolutionary models, we determine that L2 Pup A had a near-solar main sequence mass, and is therefore a close analog of the future Sun in 5 to 6 Gyr. The continuum map reveals the presence of a secondary source (B) at a radius of 2 AU contributing fB/ fA = 1.3 +/- 0.1% of the flux of the AGB star. L2 Pup B is also detected in CO emission lines at a radial velocity of vB = 12.2 +/- 1.0 km/s. The close coincidence of the center of rotation of the gaseous disk with the position of the continuum emission from the AGB star allows us to constrain the mass of the companion to mB = 12 +/- 16 MJup. L2 Pup B is most likely a planet or low mass brown dwarf with an orbital period around 5 years. Its continuum brightness and molecular emission suggest that it may be surrounded by an extended molecular atmosphere or an accretion disk. L2 Pup therefore emerges as a promising vantage point on the distant future of our Solar System. △ Less

Submitted 18 November, 2016; originally announced November 2016.

Comments: 18 pages, 12 figures, Astronomy & Astrophysics, in press

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)

Simplified models of stellar wind anatomy for interpreting high-resolution data: Analytical approach to embedded spiral geometries

Authors: Ward Homan, Leen Decin, Alex de Koter, Allard Jan van Marle, Robin Lombaert, Wouter Vlemmings

Abstract: Recent high-resolution observations have shown stellar winds to harbour complexities which strongly deviate from spherical symmetry, generally assumed as standard wind model. One such morphology is the archimedean spiral, generally believed to be formed by binary interactions, which has been directly observed in multiple sources. We seek to investigate the manifestation in the observables of spira… ▽ More Recent high-resolution observations have shown stellar winds to harbour complexities which strongly deviate from spherical symmetry, generally assumed as standard wind model. One such morphology is the archimedean spiral, generally believed to be formed by binary interactions, which has been directly observed in multiple sources. We seek to investigate the manifestation in the observables of spiral structures embedded in the spherical outflows of cool stars. We aim to provide an intuitive bedrock with which upcoming ALMA data can be compared and interpreted. By means of an extended parameter study, we model rotational CO emission from the stellar outflow of asymptotic giant branch stars. To this end, we develop a simplified analytical parametrised description of a 3D spiral structure. This model is embedded into a spherical wind, and fed into the 3D radiative transfer code LIME, which produces 3D intensity maps throughout velocity space. Subsequently, we investigate the spectral signature of rotational transitions of CO of the models, as well as the spatial aspect of this emission by means of wide-slit PV diagrams. Additionally, the potential for misinterpretation of the 3D data in a 1D context is quantified. Finally, we simulate ALMA observations to explore the impact of interefrometric noise and artifacts on the emission signatures. The spectral signatures of the CO rotational transition v=0 J=3-2 are very efficient at concealing the dual nature of the outflow. Only a select few parameter combinations allow for the spectral lines to disclose the presence of the spiral structure. The inability to disentangle the spiral from the spherical signal can result in an incorrect interpretation in a 1D context. Consequently, erroneous mass loss rates would be calculated... △ Less

Submitted 5 May, 2015; v1 submitted 20 April, 2015; originally announced April 2015.

Journal ref: A&A 579, A118 (2015)

Gravity-mode period spacings as seismic diagnostic for a sample of gamma Doradus stars from Kepler space photometry and high-resolution ground-based spectroscopy

Authors: T. Van Reeth, A. Tkachenko, C. Aerts, P. I. Papics, S. A. Triana, K. Zwintz, P. Degroote, J. Debosscher, S. Bloemen, V. S. Schmid, K. De Smedt, Y. Fremat, A. S. Fuentes, W. Homan, M. Hrudkova, R. Karjalainen, R. Lombaert, P. Nemeth, R. Oestensen, G. Van De Steene, J. Vos, G. Raskin, H. Van Winckel

Abstract: Gamma Doradus stars (hereafter gamma Dor stars) are gravity-mode pulsators of spectral type A or F. Such modes probe the deep stellar interior, offering a detailed fingerprint of their structure. Four-year high-precision space-based Kepler photometry of gamma Dor stars has become available, allowing us to study these stars with unprecedented detail. We selected, analysed, and characterized a sampl… ▽ More Gamma Doradus stars (hereafter gamma Dor stars) are gravity-mode pulsators of spectral type A or F. Such modes probe the deep stellar interior, offering a detailed fingerprint of their structure. Four-year high-precision space-based Kepler photometry of gamma Dor stars has become available, allowing us to study these stars with unprecedented detail. We selected, analysed, and characterized a sample of 67 gamma Dor stars for which we have Kepler observations available. For all the targets in the sample we assembled high-resolution spectroscopy to confirm their F-type nature. We found fourteen binaries, among which four single-lined binaries, five double-lined binaries, two triple systems and three binaries with no detected radial velocity variations. We estimated the orbital parameters whenever possible. For the single stars and the single-lined binaries, fundamental parameter values were determined from spectroscopy. We searched for period spacing patterns in the photometric data and identified this diagnostic for 50 of the stars in the sample, 46 of which are single stars or single-lined binaries. We found a strong correlation between the spectroscopic vsini and the period spacing values, confirming the influence of rotation on gamma Dor-type pulsations as predicted by theory. We also found relations between the dominant g-mode frequency, the longest pulsation period detected in series of prograde modes, vsini, and log Teff. △ Less

Submitted 8 April, 2015; originally announced April 2015.

Comments: 61 pages, 61 figures, 6 tables, accepted for publication in ApJS