Esseldeurs, M., Mathis, S., & Decin, L. (2024). Tidal Dissipation in Evolved Low and Intermediate Mass Stars. A&A, 690, A266.
@article{Esseldeurs2024,
author = {{Esseldeurs}, M. and {Mathis}, S. and {Decin}, L.},
title = {{Tidal Dissipation in Evolved Low and Intermediate Mass Stars}},
journal = {A&A},
keywords = {methods: numerical, planet-star interactions, binaries: close, stars: evolution, planetary systems, Astrophysics - Solar and Stellar Astrophysics, Astrophysics - Earth and Planetary Astrophysics},
year = {2024},
month = oct,
volume = {690},
eid = {A266},
pages = {A266},
doi = {10.1051/0004-6361/202449648},
archiveprefix = {arXiv},
eprint = {2407.10573},
primaryclass = {astro-ph.SR},
adsurl = {https://ui.adsabs.harvard.edu/abs/2024A&A...690A.266E},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
Context. As the observed occurrence for planets or stellar companions orbiting low- and intermediate-mass evolved stars is increasing, so is the importance of understanding and evaluating the strength of their interactions. This is important for the further evolution of both our own Earth-Sun system and most of the observed exoplanetary systems. One of the most fundamental mechanisms behind this interaction is the tidal dissipation in these stars, as it is one of the engines of the orbital and rotational evolution of star-planet and star-star systems.
Aims. This article builds upon previous works that studied the evolution of the tidal dissipation along the pre-main sequence and the main sequence of low- and intermediate-mass stars and found a strong link between the structural and rotational evolution of stars and tidal dissipation. This article provides, for the first time, a complete picture of tidal dissipation along the entire evolution of low- and intermediate-mass stars, including the advanced phases of evolution.
Methods. Using stellar evolutionary models, the internal structure of the star was computed from the pre-main sequence all the way up to the white dwarf phase for stars with initial masses between 1 and 4 Msun. Using this internal structure, the tidal dissipation was computed along the entire stellar evolution. Tidal dissipation was separated into two components: the dissipation of the equilibrium (non-wave-like) tide and the dissipation of the dynamical (wave-like) tide. For evolved stars, the dynamical tide is constituted by progressive internal gravity waves. The evolution of the tidal dissipation was investigated for both the equilibrium and dynamical tides, and the results were compared.
Results. The significance of both the equilibrium and dynamical tide dissipation becomes apparent within distinct domains of the parameter space. The dissipation of the equilibrium tide is dominant when the star is large or the companion is far from the star. Conversely, the dissipation of the dynamical tide is important when the star is small or the companion is close to the star. The size and location of these domains depend on the masses of both the star and the companion, as well as on the evolutionary phase.
Conclusions. Both the equilibrium and the dynamical tides are important in evolved stars, and therefore both need to be taken into account when studying the tidal dissipation in evolved stars and the evolution of the planetary and/or stellar companions orbiting them.
Esseldeurs, M., Mathis, S., & Decin, L. (2024). Towards a complete picture of the evolution of planetary systems around evolved stars. In A. Lemaitre & A.-S. Libert (Eds.), IAU Symposium (Vol. 382, pp. 98–100).
@inproceedings{Esseldeurs2023b,
author = {{Esseldeurs}, Mats and {Mathis}, St{\'e}phane and {Decin}, Leen},
title = {{Towards a complete picture of the evolution of planetary systems around evolved stars}},
keywords = {AGB stars, stellar winds and mass-loss, tides, star-planet interactions},
booktitle = {IAU Symposium},
year = {2024},
editor = {{Lemaitre}, Anne and {Libert}, Anne-Sophie},
series = {IAU Symposium},
volume = {382},
month = jan,
pages = {98-100},
doi = {10.1017/S1743921323004817},
adsurl = {https://ui.adsabs.harvard.edu/abs/2024IAUS..382...98E},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
Solar-like stars evolve through the Asymptotic Giant Branch (AGB) phase. This phase is characterized by increased radii, high luminosities, and significant mass loss. In order to understand the survival of companions during this phase, and explain the presence of planets orbiting white dwarfs, it is essential to examine the orbital evolution of these systems. Several physical mechanisms come into play for AGB stars, such as stellar mass-loss and tidal interactions between the star and its companion. On the one hand, evaluating mass-loss rates and accretion to the companion requires complex radiation-hydro-chemical simulations. On the other hand, the full history of tidal dissipation in low-mass stars during their late stages of evolution, which strongly depends on internal structure and boundary conditions, still requires dedicated studies. Finally a simultaneous treatment of winds and tides is required to predict a planet’s orbital evolution.
Esseldeurs, M., Siess, L., De Ceuster, F., Homan, W., Malfait, J., Maes, S., Konings, T., Ceulemans, T., & Decin, L. (2023). 3D simulations of AGB stellar winds. II. Ray-tracer implementation and impact of radiation on the outflow morphology. A&A, 674, A122.
@article{Esseldeurs2023a,
author = {{Esseldeurs}, M. and {Siess}, L. and {De Ceuster}, F. and {Homan}, W. and {Malfait}, J. and {Maes}, S. and {Konings}, T. and {Ceulemans}, T. and {Decin}, L.},
title = {{3D simulations of AGB stellar winds. II. Ray-tracer implementation and impact of radiation on the outflow morphology}},
journal = {A&A},
keywords = {stars: winds, outflows, methods: numerical, hydrodynamics, stars: AGB and post-AGB, radiative transfer, Astrophysics - Solar and Stellar Astrophysics, Astrophysics - Astrophysics of Galaxies},
year = {2023},
month = jun,
volume = {674},
eid = {A122},
pages = {A122},
doi = {10.1051/0004-6361/202346282},
archiveprefix = {arXiv},
eprint = {2304.09786},
primaryclass = {astro-ph.SR},
adsurl = {https://ui.adsabs.harvard.edu/abs/2023A&A...674A.122E},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
Context. Stars with an initial mass below ∼ 8 M⊙ will evolve through the Asymptotic Giant Branch (AGB) phase, during which they develop a strong stellar wind, due to radiation pressure on newly formed dust grains. Recent observations have revealed significant morphological complexities in AGB outflows, which are most probably caused by the interaction with a companion.
Aims. We aim for a more accurate description of AGB wind morphologies by accounting for both the radiation force in dust-driven winds and the impact of a companion on the AGB wind morphology.
Methods. We present the implementation of a ray-tracer for radiative transfer in the smoothed particle hydrodynamics (SPH) code Phantom. Our method allows for the creation of a 3D map of the optical depth around the AGB star. The effects of 4 different descriptions of radiative transfer, with different degrees of complexity, are compared: the free-wind approximation, the geometrical approximation, the Lucy approximation, and the attenuation approximation. Finally, we compare the Lucy and attenuation approximation to predictions with the 3D radiative transfer code Magritte.
Results. The effects of the different radiative transfer treatments are analysed considering both a low and high mass-loss rate regime and this both in the case of a single AGB star, as well as for an AGB binary system. 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 wind morphology for AGB binary systems. 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 the predictions from Magritte, we show that for most of the model the Lucy approximation works best, although close to the companion artificial heating occurs; and that it fails to simulate the shadow cast by the companion. The attenuation approximation leads to a stronger absorption of the radiation field yielding lower equilibrium temperature, weaker radiation force but produces the shadow cast by the companion. From the predictions of the 3D radiative transfer code Magritte, we also conclude that a radially directed radiation force is a reasonable assumption.
Conclusions. The radiation force 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, but for high mass-loss rates a more rigorous method is required. Among the studied approaches, the Lucy approximation provides the most accurate results, although it does not account for all effects.
Thesis
Esseldeurs, M. (2022). SPH Approach to Modelling Dust Attenuation of the Wind Acceleration in AGB Binaries. KU Leuven. Faculteit Wetenschappen.
@thesis{Esseldeurs2022,
language = {eng},
publisher = {KU Leuven. Faculteit Wetenschappen},
title = {SPH Approach to Modelling Dust Attenuation of the Wind Acceleration in AGB Binaries},
year = {2022},
author = {Esseldeurs, Mats},
address = {Leuven}
}
Circumstellar envelopes of asymptotic giant branch (AGB) stars have for a long time been modelled assuming a spherical symmetry. However, recent observations of the inner winds of these stars have shown that they exhibit a variety of complex structures, such as spirals, equatorial density enhancements, disks, bipolar outflows, etc. These structures are believed to originate from the presence of a companion, either stellar and/or planetary. Due to the inherently three-dimensional (3D) nature of this phenomenon, its investigation requires advanced 3D-simulations. It has already been shown that using hydrodynamic simulations, some of these wind morphologies can be obtained. However, the computational cost of truly selfconsistent calculations, including the crucial chemical and radiation processes, is currently still computationally prohibitive. Therefore, incremental modelling improvements using ever more refined approximations can significantly advance the quality and physical consistency of the simulations. One such problem revolves around the transfer of momentum from the stellar photons to the dust particles in the AGB wind. Until recently, this coupling has always either been ignored, or treated in the optically thin limit, which only requires knowledge on the local quantities. However, properly accounting for the dust opacity and its attenuation of the stellar radiation field can drastically affect the effective radiation pressure on the dust. In turn, this can significantly alter the dynamics and morphology of these stellar winds.
In the context of a large collaboration that revolves around upgrading the treatment of dusty winds within the smoothed-particle-hydrodynamics code PHANTOM, I have worked on levelling up the way in which dust acceleration is calculated. To calculate the attenuation of the stellar radiation field, knowledge on the distribution and optical properties of the matter in-between each particle and the star is required. To this end the ray-tracing algorithm of the radiative transfer code MAGRITTE was extracted, made compatible with the SPH philosophy, and coupled to PHANTOM. We investigated different options to speed up the ray-tracer with only a minimal loss in accuracy. We show that the best results are found when rays are traced outwards from the star in a uniform distribution, set by HEALPix. Because not all particles are struck by a ray, we also investigated different interpolation approaches. We find that interpolation scaling with the inverse square of the perpendicular distance to the four closest rays gives the most desirable result. Finally, we demonstrate the validity of our new approach in a simulation.
Co-author publications
Danilovich, T., Samaratunge, N., Mori, Y., Richards, A. M. S., Baudry, A., Etoka, S., Montargès, M., Kervella, P., McDonald, I., Gottlieb, C. A., Wallace, A., Price, D. J., Decin, L., Bolte, J., Ceulemans, T., De Ceuster, F., de Koter, A., Dionese, D., El Mellah, I., … Zijlstra, A. (2025). ATOMIUM: Dust and tracers of binarity in the continua. ArXiv e-Prints, arXiv:2504.00517.
@article{Danilovich2025,
author = {{Danilovich}, T. and {Samaratunge}, N. and {Mori}, Y. and {Richards}, A.~M.~S. and {Baudry}, A. and {Etoka}, S. and {Montarg{\`e}s}, M. and {Kervella}, P. and {McDonald}, I. and {Gottlieb}, C.~A. and {Wallace}, A. and {Price}, D.~J. and {Decin}, L. and {Bolte}, J. and {Ceulemans}, T. and {De Ceuster}, F. and {de Koter}, A. and {Dionese}, D. and {El Mellah}, I. and {Esseldeurs}, M. and {Gray}, M. and {Herpin}, F. and {Khouri}, T. and {Lagadec}, E. and {Landri}, C. and {Marinho}, L. and {Menten}, K.~M. and {Millar}, T.~J. and {M{\"u}ller}, H.~S.~P. and {Pimpanuwat}, B. and {Plane}, J.~M.~C. and {Sahai}, R. and {Siess}, L. and {Van de Sande}, M. and {Vermeulen}, O. and {Wong}, K.~T. and {Yates}, J. and {Zijlstra}, A.},
title = {{ATOMIUM: Dust and tracers of binarity in the continua}},
journal = {arXiv e-prints},
keywords = {Solar and Stellar Astrophysics, Astrophysics of Galaxies},
year = {2025},
month = apr,
eid = {arXiv:2504.00517},
pages = {arXiv:2504.00517},
doi = {10.48550/arXiv.2504.00517},
archiveprefix = {arXiv},
eprint = {2504.00517},
primaryclass = {astro-ph.SR},
adsurl = {https://ui.adsabs.harvard.edu/abs/2025arXiv250400517D},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
Low- and intermediate-mass stars on the asymptotic giant branch (AGB) account for a significant portion of the dust and chemical enrichment in their host galaxy. Here we present ALMA observations of the continuum emission at 1.24 mm around a sample of 17 stars from the ATOMIUM survey. From our analysis of the stellar contributions to the continuum flux, we find that the semi-regular variables all have smaller physical radii and fainter monochromatic luminosities than the Mira variables. Comparing these properties with pulsation periods, we find a positive trend between stellar radius and period only for the Mira variables with periods above 300 days and a positive trend between the period and the monochromatic luminosity only for the red supergiants and the most extreme AGB stars with periods above 500 days. We find that the continuum emission at 1.24 mm can be classified into four groups. "Featureless" continuum emission is confined to the (unresolved) regions close to the star for five stars in our sample, relatively uniform extended flux is seen for four stars, tentative bipolar features are seen for three stars, and the remaining five stars have unique or unusual morphological features in their continuum maps. These features can be explained by binary companions to 10 out of the 14 AGB stars in our sample. Based on our results we conclude that there are two modes of dust formation: well established pulsation-enhanced dust formation and our newly proposed companion-enhanced dust formation. If the companion is located close to the AGB star, in the wind acceleration region, then additional dust formed in the wake of the companion can increase the mass lost through the dust driven wind. This explains the different dust morphologies seen around our stars and partly accounts for a large scatter in literature mass-loss rates, especially among semiregular stars with small pulsation periods.
Willcox, R., Marchant, P., Vigna-Gómez, A., Sana, H., Bodensteiner, J., Deshmukh, K., Esseldeurs, M., Fabry, M., Hénault-Brunet, V., Janssens, S., Mahy, L., Patrick, L., Pauli, D., Renzo, M., Sander, A. A. C., Shenar, T., van Son, L. A. C., & Stoop, M. (2025). Binarity at LOw Metallicity (BLOeM): Bayesian inference of natal kicks from inert black hole binaries. ArXiv e-Prints, arXiv:2504.16669.
@article{Willcox2025,
author = {{Willcox}, R. and {Marchant}, P. and {Vigna-G{\'o}mez}, A. and {Sana}, H. and {Bodensteiner}, J. and {Deshmukh}, K. and {Esseldeurs}, M. and {Fabry}, M. and {H{\'e}nault-Brunet}, V. and {Janssens}, S. and {Mahy}, L. and {Patrick}, L. and {Pauli}, D. and {Renzo}, M. and {Sander}, A.~A.~C. and {Shenar}, T. and {van Son}, L.~A.~C. and {Stoop}, M.},
title = {{Binarity at LOw Metallicity (BLOeM): Bayesian inference of natal kicks from inert black hole binaries}},
journal = {arXiv e-prints},
keywords = {Solar and Stellar Astrophysics, Astrophysics of Galaxies, High Energy Astrophysical Phenomena},
year = {2025},
month = apr,
eid = {arXiv:2504.16669},
pages = {arXiv:2504.16669},
doi = {10.48550/arXiv.2504.16669},
archiveprefix = {arXiv},
eprint = {2504.16669},
primaryclass = {astro-ph.SR},
adsurl = {https://ui.adsabs.harvard.edu/abs/2025arXiv250416669W},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
Context. The emerging population of inert black hole binaries (BHBs) provides a unique opportunity to constrain black hole (BH) formation physics. These systems are composed of a stellar-mass BH in a wide orbit around a non-degenerate star with no observed Xray emission. Inert BHBs allow for narrow constraints to be inferred on the natal kick and mass loss during BH-forming core-collapse events.
Aims. In anticipation of the upcoming BLOeM survey, we aim to provide tight constraints on BH natal kicks by exploiting the full parameter space obtained from combined spectroscopic and astrometric data to characterize the orbits of inert BHBs. Multi-epoch spectroscopy from the BLOeM project will provide measurements of periods, eccentricities, and radial velocities for inert BHBs in the SMC, which complements Gaia astrometric observations of proper motions.
Methods. We present a Bayesian parameter estimation framework to infer natal kicks and mass loss during core-collapse from inert BHBs, accounting for all available observables, including the systemic velocity and its orientation relative to the orbital plane. The framework further allows for circumstances when some of the observables are unavailable, such as for the distant BLOeM sources which preclude resolved orbits.
Results. With our new framework, we are able to distinguish between BH formation channels, even in the absence of a resolved orbit. In cases when the pre-explosion orbit can be assumed to be circular, we precisely recover the parameters of the core-collapse, highlighting the importance of understanding the eccentricity landscape of pre-explosion binaries, both theoretically and observationally. Treating the near-circular, inert BHB, VFTS 243, as a representative of the anticipated BLOeM systems, we constrain the natal kick to less than 27 km/s and the mass loss to less than 2.9 Msun within a 90% credible interval.
Montargès, M., Malfait, J., Esseldeurs, M., de Koter, A., Baron, F., Kervella, P., Danilovich, T., Richards, A. M. S., Sahai, R., McDonald, I., Khouri, T., Shetye, S., Zijlstra, A., Van de Sande, M., El Mellah, I., Herpin, F., Siess, L., Etoka, S., Gobrecht, D., … Yates, aJ. (2025). An accreting dwarf star orbiting the S-type giant star pi1 Gru. ArXiv e-Prints, arXiv:2504.16845.
@article{Montarges2025,
author = {{Montarg{\`e}s}, M. and {Malfait}, J. and {Esseldeurs}, M. and {de Koter}, A. and {Baron}, F. and {Kervella}, P. and {Danilovich}, T. and {Richards}, A.~M.~S. and {Sahai}, R. and {McDonald}, I. and {Khouri}, T. and {Shetye}, S. and {Zijlstra}, A. and {Van de Sande}, M. and {El Mellah}, I. and {Herpin}, F. and {Siess}, L. and {Etoka}, S. and {Gobrecht}, D. and {Marinho}, L. and {Wallstr{\"o}m}, S.~H.~J. and {Wong}, K.~T. and {Yates}, aJ.},
title = {{An accreting dwarf star orbiting the S-type giant star pi1 Gru}},
journal = {arXiv e-prints},
keywords = {Solar and Stellar Astrophysics},
year = {2025},
month = apr,
eid = {arXiv:2504.16845},
pages = {arXiv:2504.16845},
doi = {10.48550/arXiv.2504.16845},
archiveprefix = {arXiv},
eprint = {2504.16845},
primaryclass = {astro-ph.SR},
adsurl = {https://ui.adsabs.harvard.edu/abs/2025arXiv250416845M},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
Aims. We aim to characterize the properties of the inner companion of the S-type AGB star pi1 Gru, and to identify plausible future evolution scenarios for this triple system.
Methods. We observed pi1 Gru with ALMA and VLT/SPHERE. In addition, we collected archival photometry data and used the Hipparcos-Gaia proper motion anomaly. We derive the best orbital parameters from Bayesian inference.
Results. The inner companion, pi1 Gru C was located at 37.4 +/- 2.0 mas from the primary in June-July 2019 (projected separation of 6.05 +/- 0.55 au at 161.7 +/- 11.7 pc). The best orbital solution gives a companion mass of 0.86 (+0.22/-0.20) Msun (using the derived mass of the primary), and a semi-major axis of 7.05 (+0.54/-0.57) au. This leads to an orbital period of 11.0 (+1.7/-1.5) yr. The best solution is an elliptical orbit with eccentricity e = 0.35 (+0.18/-0.17), but a circular orbit cannot be totally excluded. The close companion can either be a K1V (F9.5V/K7V) star or a white dwarf. The ultraviolet and millimeter continuum photometry are consistent with the presence of an accretion disk around the close companion. The ultraviolet emission could then either originate in hot spots in an overall cooler disk, or also from a hot disk in case the companion is a white dwarf.
Conclusions. Though the close companion and the AGB star are interacting, and an accretion disk is observed around the companion, the mass-accretion rate is too low to cause a Ia supernova but could produce novae every 900 yr. Short wavelength spatially resolved observations are needed to further constrain the nature of the C companion. Searches for close-in companions similar to this system will help to better understand the physics of mass- and angular-momentum transfer, and orbital evolution in the late evolutionary stages.
Malfait, J., Siess, L., Esseldeurs, M., De Ceuster, F., Wallström, S. H. J., de Koter, A., & Decin, L. (2024). Impact of H I cooling and study of accretion disks in asymptotic giant branch wind-companion smoothed particle hydrodynamic simulations. Astronomy & Astrophysics, 691, A84.
@article{Malfait2024a,
author = {{Malfait}, J. and {Siess}, L. and {Esseldeurs}, M. and {De Ceuster}, F. and {Wallstr{\"o}m}, S.~H.~J. and {de Koter}, A. and {Decin}, L.},
title = {{Impact of H I cooling and study of accretion disks in asymptotic giant branch wind-companion smoothed particle hydrodynamic simulations}},
journal = {Astronomy & Astrophysics},
keywords = {methods: numerical, stars: AGB and post-AGB, binaries: general, stars: winds, outflows, Astrophysics - Solar and Stellar Astrophysics},
year = {2024},
month = nov,
volume = {691},
eid = {A84},
pages = {A84},
doi = {10.1051/0004-6361/202450338},
archiveprefix = {arXiv},
eprint = {2408.13158},
primaryclass = {astro-ph.SR},
adsurl = {https://ui.adsabs.harvard.edu/abs/2024A&A...691A..84M},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
Context. High-resolution observations reveal that the outflows of evolved low- and intermediate-mass stars harbour complex morphological structures that are linked to the presence of one or multiple companions. Hydrodynamical simulations provide a way to study the impact of a companion on the shaping of the asymptotic giant branch (AGB) star out-flow.
Aims. Using smoothed particle hydrodynamic (SPH) simulations of an AGB star undergoing mass loss, which also has a binary companion, we study the impact of including H I atomic line cooling on the flow morphology. We also study how this affects the properties of the accretion disks that form around the companion.
Methods. We used the PHANTOM code to perform high-resolution 3D SPH simulations of the interaction of a solar-mass companion with the outflow of an AGB star, using different wind velocities and eccentricities. We compared the model properties, computed with and without the inclusion of H I cooling.
Results. The inclusion of H I cooling produces a sizeable decrease in the temperature, up to one order of magnitude, in the region closely surrounding the companion star. As a consequence, the morphological irregularities and relatively energetic (bipolar) outflows that were obtained without H I cooling no longer appear. In the case of an eccentric orbit and a low wind velocity, these morphologies are still highly asymmetric, but the same structures recur at every orbital period, making the morphology more regular. Flared accretion disks, with a (sub-)Keplerian velocity profile, are found to form around the companion in all our models with H I cooling, provided the accretion radius is small enough. The disks have radial sizes ranging from about 0.4 to 0.9 au and masses around 10‑7‑10‑8 M⊙. For the considered wind velocities, mass accretion onto the companion is up to a factor of 2 higher than predicted by the standard Bondi Hoyle Littleton rate, ranging between 4 to 21% of the AGB wind mass loss rate. The lower the wind velocity at the location of the companion, the larger and the more massive the disk and the higher the mass accretion efficiency. In eccentric systems, the disk size, disk mass, and mass accretion efficiency vary, depending on the orbital phase.
Conclusions. H I cooling is an essential ingredient to properly model the medium around the companion where density-enhanced wind structures form and it favours the formation of an accretion disk.
Malfait, J., Siess, L., Vermeulen, O., Esseldeurs, M., Wallström, S. H. J., Richards, A. M. S., De Ceuster, F., Maes, S., Bolte, J., & Decin, L. (2024). SPH modelling of AGB wind morphology in hierarchical triple systems and a comparison to observation of R Aql. Astronomy & Astrophysics, 691, A57.
@article{Malfait2024b,
author = {{Malfait}, J. and {Siess}, L. and {Vermeulen}, O. and {Esseldeurs}, M. and {Wallstr{\"o}m}, S.~H.~J. and {Richards}, A.~M.~S. and {De Ceuster}, F. and {Maes}, S. and {Bolte}, J. and {Decin}, L.},
title = {{SPH modelling of AGB wind morphology in hierarchical triple systems and a comparison to observation of R Aql}},
journal = {Astronomy & Astrophysics},
keywords = {hydrodynamics, methods: numerical, stars: AGB and post-AGB, binaries: close, stars: winds, outflows, ISM: jets and outflows, Astrophysics - Solar and Stellar Astrophysics},
year = {2024},
month = nov,
volume = {691},
eid = {A57},
pages = {A57},
doi = {10.1051/0004-6361/202450844},
archiveprefix = {arXiv},
eprint = {2408.16565},
primaryclass = {astro-ph.SR},
adsurl = {https://ui.adsabs.harvard.edu/abs/2024A&A...691A..57M},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
Context. Complex asymmetric 3D structures are observed in the outflows of evolved low- and intermediate-mass stars, and are believed to be shaped through the interaction of companions that remain hidden within the dense wind. One example is the AGB star R Aql, for which ALMA (Atacama Large Millimeter Array) observations have revealed complex wind structures that might originate from a higher-order stellar system.
Aims. We investigate how triple systems can shape the outflow of asymptotic giant branch (AGB) stars and characterise the different wind structures that form. For simplicity, we solely focus on co-planar systems in a hierarchical, stable orbit, consisting of an AGB star with one relatively close companion, and another one at a large orbital separation.
Methods. We modelled a grid of hierarchical triple systems including a wind-launching AGB star, with the smoothed-particle- hydrodynamic PHANTOM code. We varied the outer companion mass, the AGB wind velocity, and the orbital eccentricities to study the impact of these parameters on the wind morphology. To study the impact of adding a triple companion, we additionally modelled and analysed a small grid of binary sub-systems, for comparison. To investigate if R Aql could be shaped by a triple system, we postprocessed one of our triple models with a radiative transfer routine, and compared this to data of the ALMA ATOMIUM programme.
Results. The characteristic wind structures resulting from a hierarchical triple system are the following. A large two-edged spiral wake results behind the outer companion star. This structure lies on top of the spiral structure formed by the close binary, which is itself affected by the orbital motion around the system’s centre of mass, such that it resembles a snail-shell pattern. This dense inner spiral pattern interacts with, and strongly impacts, the spiral wake of the outer companion, resulting in a wave pattern on the outer edge of this spiral wake. The higher the mass of the outer companion, the larger the density enhancement and the more radially compressed the outer spiral. Lowering the wind velocity has a similar effect, and additionally results in an elongation of the global wind morphology. Introducing eccentricity in the inner and outer orbit of the hierarchical system results in complex phase-dependent wind-companion interactions, and consequently in asymmetries in the inner part of the wind and the global morphology, respectively. From the comparison of our models to the observations of R Aql, we conclude that this circumstellar environment might be shaped by a similar system to the ones modelled in this work, but an elaborate study of the observational data is needed to better determine the orbital parameters and characteristics of the central system.
Conclusions. The modelled outflow of an AGB star in a co-planar hierarchical systems is characterised by a large-scale spiral wake with a wavey outer edge, attached to the outer companion, on top of a compact inner spiral pattern that resembles a snail-shell pattern.
Wallström, S. H. J., Danilovich, T., Müller, H. S. P., Gottlieb, C. A., Maes, S., Van de Sande, M., Decin, L., Richards, A. M. S., Baudry, A., Bolte, J., Ceulemans, T., De Ceuster, F., de Koter, A., El Mellah, I., Esseldeurs, M., Etoka, S., Gobrecht, D., Gottlieb, E., Gray, M., … Zijlstra, A. (2024). ATOMIUM: Molecular inventory of 17 oxygen-rich evolved stars observed with ALMA. Astronomy & Astrophysics, 681, A50.
@article{Wallstrom2024,
author = {{Wallstr{\"o}m}, S.~H.~J. and {Danilovich}, T. and {M{\"u}ller}, H.~S.~P. and {Gottlieb}, C.~A. and {Maes}, S. and {Van de Sande}, M. and {Decin}, L. and {Richards}, A.~M.~S. and {Baudry}, A. and {Bolte}, J. and {Ceulemans}, T. and {De Ceuster}, F. and {de Koter}, A. and {El Mellah}, I. and {Esseldeurs}, M. and {Etoka}, S. and {Gobrecht}, D. and {Gottlieb}, E. and {Gray}, M. and {Herpin}, F. and {Jeste}, M. and {Kee}, D. and {Kervella}, P. and {Khouri}, T. and {Lagadec}, E. and {Malfait}, J. and {Marinho}, L. and {McDonald}, I. and {Menten}, K.~M. and {Millar}, T.~J. and {Montarg{\`e}s}, M. and {Nuth}, J.~A. and {Plane}, J.~M.~C. and {Sahai}, R. and {Waters}, L.~B.~F.~M. and {Wong}, K.~T. and {Yates}, J. and {Zijlstra}, A.},
title = {{ATOMIUM: Molecular inventory of 17 oxygen-rich evolved stars observed with ALMA}},
journal = {Astronomy & Astrophysics},
keywords = {stars: AGB and post-AGB, supergiants, circumstellar matter, line: identification, instrumentation: interferometers, astrochemistry, Astrophysics - Solar and Stellar Astrophysics, Astrophysics - Astrophysics of Galaxies},
year = {2024},
month = jan,
volume = {681},
eid = {A50},
pages = {A50},
doi = {10.1051/0004-6361/202347632},
archiveprefix = {arXiv},
eprint = {2312.03467},
primaryclass = {astro-ph.SR},
adsurl = {https://ui.adsabs.harvard.edu/abs/2024A&A...681A..50W},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
Context: 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.
Aims: 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.
Methods: 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.
Results: 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. Overall, we find strong correlations between the following molecular pairs: CS and SiS, CS and AlF, NaCl and KCl, AlO and SO, SO2 and SO, and SO2 and H2O; meaning both molecules tend to have more detected emission lines in the same sources. The measured isotopic ratios of Si and S are found to be consistent with previous measurements, except for an anomalously high 29Si/30Si ratio of 4 ± 1 in the RSG VX Sgr.
Conclusions: This paper presents the overall molecular inventory and an initial analysis of the large ATOMIUM dataset, laying the groundwork for future work deriving molecular abundances and abundance profiles using radiative transfer modeling which will provide more rigorous tests for chemical models.</pre>
</div>
</li>
Maes, S., Siess, L., Homan, W., Malfait, J., De Ceuster, F., Ceulemans, T., Donné, D., Esseldeurs, M., & Decin, L. (2022). Route towards complete 3D hydro-chemical simulations of companion-perturbed AGB outflows. The Origin of Outflows in Evolved Stars, 366, 227–233.
@inproceedings{Maes2022,
author = {{Maes}, Silke and {Siess}, Lionel and {Homan}, Ward and {Malfait}, Jolien and {De Ceuster}, Frederik and {Ceulemans}, Thomas and {Donn{\'e}}, Dion and {Esseldeurs}, Mats and {Decin}, Leen},
title = {{Route towards complete 3D hydro-chemical simulations of companion-perturbed AGB outflows}},
keywords = {Stars: AGB, Stars: winds, outflows, Methods: numerical, Astrophysics - Solar and Stellar Astrophysics, Astrophysics - Astrophysics of Galaxies},
booktitle = {The Origin of Outflows in Evolved Stars},
year = {2022},
volume = {366},
month = jan,
pages = {227-233},
doi = {10.1017/S1743921322000217},
archiveprefix = {arXiv},
eprint = {2206.12278},
primaryclass = {astro-ph.SR},
adsurl = {https://ui.adsabs.harvard.edu/abs/2022IAUS..366..227M},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
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, 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.