Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

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

Abstract

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

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Observed SiN emission towards W Aql.
Fig. 2: Plots of the orbit of the W Aql system as seen in the plane of the sky.
Fig. 3: Observed SiO, SiS, CS and HCN emission towards W Aql.
Fig. 4: Observed emission of CN-bearing molecules towards W Aql.
Fig. 5: Observations and simulations of CO emission towards W Aql.

Similar content being viewed by others

Data availability

The observational data used here are openly available through the data archives for ALMA (https://almascience.nrao.edu/aq/), ESO for the APEX and SPHERE data (http://archive.eso.org), and HST (https://hla.stsci.edu). Custom ALMA data products that were produced for this study are available from T.D. or A.M.S.R. (a.m.s.richards@manchester.ac.uk) upon reasonable request. The Phantom input and output files from our main model can be downloaded from https://doi.org/10.26180/24240001.v1.

Code availability

Phantom is open source under the GPLv3 license and can be downloaded via https://github.com/danieljprice/phantom. MCFOST is open source under the GPLv3 license and can be downloaded via https://mcfost.readthedocs.io/en/latest/overview.html. ALI, the one-dimensional radiative transfer code, is available from the corresponding author upon reasonable request.

References

  1. Höfner, S. & Olofsson, H. Mass loss of stars on the asymptotic giant branch. Mechanisms, models and measurements. Astron. Astrophys. Rev. 26, 1 (2018).

    ADS  Google Scholar 

  2. Kobayashi, C., Karakas, A. I. & Lugaro, M. The origin of elements from carbon to uranium. Astrophys. J. 900, 179 (2020).

    ADS  CAS  Google Scholar 

  3. Decin, L. et al. Reduction of the maximum mass-loss rate of OH/IR stars due to unnoticed binary interaction. Nat. Astron. 3, 408–415 (2019).

  4. Van de Sande, M. & Millar, T. J. The impact of stellar companion UV photons on the chemistry of the circumstellar environments of AGB stars. Mon. Not. R. Astron. Soc. 510, 1204–1222 (2022).

    ADS  Google Scholar 

  5. De Marco, O. et al. The messy death of a multiple star system and the resulting planetary nebula as observed by JWST. Nat. Astron. 6, 1421–1432 (2022).

    ADS  Google Scholar 

  6. Decin, L. et al. (Sub)stellar companions shape the winds of evolved stars. Science 369, 1497–1500 (2020).

    ADS  CAS  PubMed  Google Scholar 

  7. Jorissen, A. in Asymptotic Giant Branch Stars (eds. H. J. Habing & H. Olofsson) 461–518 (Springer, 2004).

  8. Karovska, M., Hack, W., Raymond, J. & Guinan, E. First Hubble Space Telescope observations of Mira AB wind-accreting binary systems. Astrophys. J. Lett. 482, 175–178 (1997).

    ADS  Google Scholar 

  9. Ramstedt, S. et al. The wonderful complexity of the Mira AB system. Astron. Astrophys. 570, 14 (2014).

    Google Scholar 

  10. Kervella, P. et al. ALMA observations of the nearby AGB star L2 Puppis. I. Mass of the central star and detection of a candidate planet. Astron. Astrophys. 596, 92 (2016).

    Google Scholar 

  11. Homan, W. et al. ATOMIUM: a high-resolution view on the highly asymmetric wind of the AGB star π1Gruis. I. First detection of a new companion and its effect on the inner wind. Astron. Astrophys. 644, 61 (2020).

    Google Scholar 

  12. Mauron, N. & Huggins, P. J. Imaging the circumstellar envelopes of AGB stars. Astron. Astrophys. 452, 257–268 (2006).

    ADS  CAS  Google Scholar 

  13. Kim, H. et al. The large-scale nebular pattern of a superwind binary in an eccentric orbit. Nat. Astron. 1, 0060 (2017).

    ADS  Google Scholar 

  14. Maercker, M. et al. Unexpectedly large mass loss during the thermal pulse cycle of the red giant star R Sculptoris. Nature 490, 232–234 (2012).

    ADS  CAS  PubMed  Google Scholar 

  15. Sahai, R. et al. The rapidly evolving asymptotic giant branch star, V Hya: ALMA finds a multiring circus with high-velocity outflows. Astrophys. J. 929, 59 (2022).

    ADS  Google Scholar 

  16. Ramos-Larios, G. et al. Rings and arcs around evolved stars—I. Fingerprints of the last gasps in the formation process of planetary nebulae. Mon. Not. R. Astron. Soc. 462, 610–635 (2016).

    ADS  CAS  Google Scholar 

  17. Herbig, G. H. Physical companions to long-period variables. Veroeffentlichungen der Remeis-Sternwarte zu Bamberg 27, 164 (1965).

    ADS  Google Scholar 

  18. Danilovich, T., Olofsson, G., Black, J. H., Justtanont, K. & Olofsson, H. Classifying the secondary component of the binary star W Aquilae. Astron. Astrophys. 574, 23 (2015).

    ADS  Google Scholar 

  19. Ramstedt, S., Maercker, M., Olofsson, G., Olofsson, H. & Schöier, F. L. Imaging the circumstellar dust around AGB stars with PolCor. Astron. Astrophys. 531, 148 (2011).

    ADS  Google Scholar 

  20. Mayer, A. et al. Large-scale environments of binary AGB stars probed by Herschel. I. Morphology statistics and case studies of R Aquarii and W Aquilae. Astron. Astrophys. 549, 69 (2013).

    Google Scholar 

  21. Danilovich, T. et al. Detailed modelling of the circumstellar molecular line emission of the S-type AGB star W Aquilae. Astron. Astrophys. 569, 76 (2014).

    Google Scholar 

  22. Ramstedt, S. et al. The circumstellar envelope around the S-type AGB star W Aql. Effects of an eccentric binary orbit. Astron. Astrophys. 605, 126 (2017).

    Google Scholar 

  23. Brunner, M. et al. Molecular line study of the S-type AGB star W Aquilae. ALMA observations of CS, SiS, SiO and HCN. Astron. Astrophys. 617, 23 (2018).

    Google Scholar 

  24. Danilovich, T. et al. ATOMIUM: halide molecules around the S-type AGB star W Aquilae. Astron. Astrophys. 655, 80 (2021).

    Google Scholar 

  25. Richichi, A., Percheron, I. & Khristoforova, M. CHARM2: an updated catalog of high angular resolution measurements. Astron. Astrophys. 431, 773–777 (2005).

    ADS  Google Scholar 

  26. Turner, B. E. Detection of SiN in IRC+10216. Astrophys. J. Lett. 388, 35 (1992).

    ADS  Google Scholar 

  27. Huggins, P. J., Glassgold, A. E. & Morris, M. CN and C2H in IRC+10216. Astrophys. J. 279, 284–290 (1984).

    ADS  CAS  Google Scholar 

  28. Tokovinin, A. & Kiyaeva, O. Eccentricity distribution of wide binaries. Mon. Not. R. Astron. Soc. 456, 2070–2079 (2016).

    ADS  Google Scholar 

  29. Boffin, H. M. J., Cerf, N. & Paulus, G. Statistical analysis of a sample of spectroscopic binaries containing late-type giants. Astron. Astrophys. 271, 125–138 (1993).

    ADS  Google Scholar 

  30. Moe, M. & Di Stefano, R. Mind your Ps and Qs: the interrelation between period (P) and mass-ratio (Q) distributions of binary stars. Astrophys. J. Suppl. Ser. 230, 15 (2017).

    ADS  Google Scholar 

  31. Maercker, M. et al. A detailed view of the gas shell around R Sculptoris with ALMA. Astron. Astrophys. 586, 5 (2016).

    Google Scholar 

  32. Lykou, F. et al. The curious case of II Lup: a complex morphology revealed with SAM/NACO and ALMA. Mon. Not. R. Astron. Soc. 480, 1006–1021 (2018).

    ADS  Google Scholar 

  33. Doan, L. et al. The extended molecular envelope of the asymptotic giant branch star π1 Gruis as seen by ALMA. I. Large-scale kinematic structure and CO excitation properties. Astron. Astrophys. 605, 28 (2017).

    Google Scholar 

  34. Sahai, R., Scibelli, S. & Morris, M. R. High-speed bullet ejections during the AGB-to-planetary nebula transition: HST observations of the carbon star, V Hydrae. Astrophys. J. 827, 92 (2016).

    ADS  Google Scholar 

  35. Oomen, G.-M. et al. Orbital properties of binary post-AGB stars. Astron. Astrophys. 620, 85 (2018).

    Google Scholar 

  36. Jorissen, A., Van Eck, S., Mayor, M. & Udry, S. Insights into the formation of barium and Tc-poor S stars from an extended sample of orbital elements. Astron. Astrophys. 332, 877–903 (1998).

    ADS  CAS  Google Scholar 

  37. Ramstedt, S., Schöier, F. L. & Olofsson, H. Circumstellar molecular line emission from S-type AGB stars: mass-loss rates and SiO abundances. Astron. Astrophys. 499, 515–527 (2009).

    ADS  CAS  Google Scholar 

  38. Knapp, G. R. & Morris, M. Mass loss from evolved stars. III—Mass loss rates for fifty stars from CO J = 1–0 observations. Astrophys. J. 292, 640–669 (1985).

    ADS  CAS  Google Scholar 

  39. Jura, M. Mass loss from S stars. Astrophys. J. Suppl. Ser. 66, 33 (1988).

    ADS  CAS  Google Scholar 

  40. Groenewegen, M. A. T. & De Jong, T. CO observations and mass loss of MS- and S-stars. Astron. Astrophys. 337, 797–807 (1998).

    ADS  Google Scholar 

  41. Gaia Collaboration The Gaia mission. Astron. Astrophys. 595, 1 (2016).

    Google Scholar 

  42. Andriantsaralaza, M., Ramstedt, S., Vlemmings, W. H. T. & De Beck, E. Distance estimates for AGB stars from parallax measurements. Astron. Astrophys. 667, 74 (2022).

    ADS  Google Scholar 

  43. Gaia Collaboration Gaia Early Data Release 3. Summary of the contents and survey properties. Astron. Astrophys. 649, 1 (2021).

    Google Scholar 

  44. Habets, G. M. H. J. & Heintze, J. R. W. Empirical bolometric corrections for the main-sequence. Astron. Astrophys. Suppl. Ser. 46, 193–237 (1981).

    ADS  Google Scholar 

  45. De Nutte, R. et al. Nucleosynthesis in AGB stars traced by oxygen isotopic ratios. I. Determining the stellar initial mass by means of the 17O/18O ratio. Astron. Astrophys. 600, 71 (2017).

    Google Scholar 

  46. De Beck, E. & Olofsson, H. The surprisingly carbon-rich environment of the S-type star W Aql. Astron. Astrophys. 642, 20 (2020).

    Google Scholar 

  47. Vassiliadis, E. & Wood, P. R. Evolution of low- and intermediate-mass stars to the end of the asymptotic giant branch with mass loss. Astrophys. J. 413, 641–657 (1993).

    ADS  Google Scholar 

  48. Montargès, M. et al. The VLT/SPHERE view of the ATOMIUM cool evolved star sample. I. Overview: sample characterization through polarization analysis. Astron. Astrophys. 671, 96 (2023).

    Google Scholar 

  49. Maire, A.-L. et al. High-precision astrometric studies in direct imaging with SPHERE. The Messenger 183, 7–12 (2021).

    ADS  Google Scholar 

  50. Gottlieb, C. A. et al. ATOMIUM: ALMA tracing the origins of molecules in dust forming oxygen rich M-type stars. Motivation, sample, calibration, and initial results. Astron. Astrophys. 660, 94 (2022).

    Google Scholar 

  51. The CASA Team et al. CASA, the Common Astronomy Software Applications for radio astronomy. Publ. Astron. Soc. Pacific 134, 114501 (2022).

  52. Thompson, A. R., Moran, J. M. & Swenson, G. W. Jr Interferometry and Synthesis in Radio Astronomy 3rd edn (Springer, 2017).

  53. Olofsson, H., Eriksson, K., Gustafsson, B. & Carlstrom, U. A study of circumstellar envelopes around bright carbon stars. I—Structure, kinematics, and mass-loss rate. Astrophys. J. Suppl. Ser. 87, 267–304 (1993).

    ADS  CAS  Google Scholar 

  54. Homan, W. et al. 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. Astron. Astrophys. 651, 82 (2021).

    Google Scholar 

  55. Vlemmings, W. H. T., Khouri, T. & Olofsson, H. Resolving the extended stellar atmospheres of asymptotic giant branch stars at (sub)millimetre wavelengths. Astron. Astrophys. 626, 81 (2019).

    ADS  Google Scholar 

  56. Baudry, A. et al. ATOMIUM: probing the inner wind of evolved O-rich stars with new, highly excited H2O and OH lines. Astron. Astrophys. 674, 125 (2023).

    Google Scholar 

  57. Cernicharo, J., Gottlieb, C. A., Guelin, M., Thaddeus, P. & Vrtilek, J. M. Astronomical and laboratory detection of the SiC radical. Astrophys. J. Lett. 341, 25 (1989).

    ADS  Google Scholar 

  58. Massalkhi, S. et al. Abundance of SiC2 in carbon star envelopes. Astron. Astrophys. 611, 29 (2018).

    Google Scholar 

  59. Cernicharo, J., Marcelino, N., Agúndez, M. & Guélin, M. Molecular shells in IRC+10216: tracing the mass loss history. Astron. Astrophys. 575, 91 (2015).

    ADS  Google Scholar 

  60. Decin, L. et al. ALMA data suggest the presence of spiral structure in the inner wind of CW Leonis. Astron. Astrophys. 574, 5 (2015).

    Google Scholar 

  61. Siebert, M. A., Van de Sande, M., Millar, T. J. & Remijan, A. J. Investigating anomalous photochemistry in the inner wind of IRC+10216 through interferometric observations of HC3N. Astrophys. J. 941, 90 (2022).

    ADS  Google Scholar 

  62. Patel, N., Gottlieb, C. & Young, K. Probing the dust formation zone in IRC+10216 with the SMA. In Proc. The Life Cycle of Dust in the Universe: Observations (eds Kemper et al.) 98 (2013).

  63. Velilla-Prieto, L. et al. Circumstellar chemistry of Si-C bearing molecules in the C-rich AGB star IRC+10216. IAU Symp. 343, 535–537 (2019).

    ADS  Google Scholar 

  64. Lee, S. K., Ozeki, H. & Saito, S. Microwave spectrum of the NS radical in the 2Π ground electronic state. Astrophys. J. Suppl. Ser. 98, 351 (1995).

    ADS  CAS  Google Scholar 

  65. Velilla Prieto, L. et al. The millimeter IRAM-30 m line survey toward IK Tauri. Astron. Astrophys. 597, 25 (2017).

    Google Scholar 

  66. Decin, L., Richards, A. M. S., Danilovich, T., Homan, W. & Nuth, J. A. ALMA spectral line and imaging survey of a low and a high mass-loss rate AGB star between 335 and 362 GHz. Astron. Astrophys. 615, 28 (2018).

    ADS  Google Scholar 

  67. Decin, L. et al. Discovery of multiple dust shells beyond 1 arcmin in the circumstellar envelope of IRC+10216 using Herschel/PACS. Astron. Astrophys. 534, 1 (2011).

    Google Scholar 

  68. Guélin, M. et al. IRC +10 216 in 3D: morphology of a TP-AGB star envelope. Astron. Astrophys. 610, 4 (2018).

    Google Scholar 

  69. Schöier, F. L. et al. The abundance of HCN in circumstellar envelopes of AGB stars of different chemical type. Astron. Astrophys. 550, 78 (2013).

    Google Scholar 

  70. Danilovich, T. et al. Sulphur-bearing molecules in AGB stars. II. Abundances and distributions of CS and SiS. Astron. Astrophys. 617, 132 (2018).

    Google Scholar 

  71. Massalkhi, S., Agúndez, M. & Cernicharo, J. Study of CS, SiO, and SiS abundances in carbon star envelopes: assessing their role as gas-phase precursors of dust. Astron. Astrophys. 628, 62 (2019).

    ADS  Google Scholar 

  72. Morris, M. & Jura, M. Molecular self-shielding in the outflows from late-type stars. Astrophys. J. 264, 546–553 (1983).

    ADS  CAS  Google Scholar 

  73. Bachiller, R. et al. A survey of CN in circumstellar envelopes. Astron. Astrophys. 319, 235–243 (1997).

    ADS  CAS  Google Scholar 

  74. Rybicki, G. B. & Hummer, D. G. An accelerated lambda iteration method for multilevel radiative transfer. I—Non-overlapping lines with background continuum. Astron. Astrophys. 245, 171–181 (1991).

    ADS  CAS  Google Scholar 

  75. Pickett, H. M. The fitting and prediction of vibration–rotation spectra with spin interactions. J. Mol. Spectrosc. 148, 371–377 (1991).

    ADS  CAS  Google Scholar 

  76. Lique, F. et al. Rotational excitation of CN(X 2Σ+) by He: theory and comparison with experiments. J. Chem. Phys. 132, 024303 (2010).

    ADS  PubMed  Google Scholar 

  77. Montez Jr, R., Ramstedt, S., Kastner, J. H., Vlemmings, W. & Sanchez, E. A catalog of GALEX ultraviolet emission from asymptotic giant branch stars. Astrophys. J. 841, 33 (2017).

    ADS  Google Scholar 

  78. Schrijver, C. J. Magnetic structure in cool stars. XI. Relations between radiative fluxes mesuring stellar activity, and evidence for two components in stellar chromospheres. Astron. Astrophys. 172, 111–123 (1987).

    ADS  CAS  Google Scholar 

  79. Gálvez, M. C. et al. Multiwavelength optical observations of chromospherically active binary systems. IV. The X-ray/EUV selected binary BK Psc (2RE J0039+103). Astron. Astrophys. 389, 524–536 (2002).

    ADS  Google Scholar 

  80. Gray, R. O., Corbally, C. J. & Burgasser, A. J. Stellar Spectral Classification Princeton Series in Astrophysics (Princeton Univ. Press, 2009); http://books.google.com.au/books?id=S_Sh1i226wwC

  81. Asplund, M., Amarsi, A. M. & Grevesse, N. The chemical make-up of the Sun: a 2020 vision. Astron. Astrophys. 653, 141 (2021).

    ADS  Google Scholar 

  82. Agúndez, M., Martínez, J. I., de Andres, P. L., Cernicharo, J. & Martín-Gago, J. A. Chemical equilibrium in AGB atmospheres: successes, failures, and prospects for small molecules, clusters, and condensates. Astron. Astrophys. 637, 59 (2020).

    Google Scholar 

  83. Agúndez, M. et al. Growth of carbon chains in IRC+10216 mapped with ALMA. Astron. Astrophys. 601, 4 (2017).

    Google Scholar 

  84. Cordiner, M. A. & Millar, T. J. Density-enhanced gas and dust shells in a new chemical model for IRC+10216. Astrophys. J. 697, 68–78 (2009).

    ADS  CAS  Google Scholar 

  85. Kerkines, I. S. K. & Mavridis, A. On the electron affinity of SiN and spectroscopic constants of SiN. J. Chem. Phys. 123, 124301–124301 (2005).

    ADS  PubMed  Google Scholar 

  86. Thomson, R. & Dalby, F. W. Experimental determination of the dipole moments of the X(2Σ+) and B(2Σ+) states of the CN molecule. Can. J. Phys. 46, 2815 (1968).

    ADS  CAS  Google Scholar 

  87. Monaghan, J. J. Smoothed particle hydrodynamics. Rep. Prog. Phys. 68, 1703–1759 (2005).

    ADS  MathSciNet  Google Scholar 

  88. Price, D. J. Smoothed particle hydrodynamics and magnetohydrodynamics. J. Comput. Phys. 231, 759–794 (2012).

    ADS  MathSciNet  Google Scholar 

  89. Price, D. J. et al. Phantom: a smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics. Publ. Astron. Soc. Aust. 35, 031 (2018).

    ADS  Google Scholar 

  90. Siess, L., Homan, W., Toupin, S. & Price, D. J. 3D simulations of AGB stellar winds. I. Steady winds and dust formation. Astron. Astrophys. 667, 75 (2022).

    ADS  Google Scholar 

  91. Maes, S. et al. SPH modelling of companion-perturbed AGB outflows including a new morphology classification scheme. Astron. Astrophys. 653, 25 (2021).

    Google Scholar 

  92. Malfait, J. et al. SPH modelling of wind-companion interactions in eccentric AGB binary systems. Astron. Astrophys. 652, 51 (2021).

    Google Scholar 

  93. Mastrodemos, N. & Morris, M. Bipolar preplanetary nebulae: hydrodynamics of dusty winds in binary systems. I. Formation of accretion disks. Astrophys. J. 497, 303–329 (1998).

    ADS  CAS  Google Scholar 

  94. Mastrodemos, N. & Morris, M. Bipolar pre-planetary nebulae: hydrodynamics of dusty winds in binary systems. II. Morphology of the circumstellar envelopes. Astrophys. J. 523, 357–380 (1999).

    ADS  Google Scholar 

  95. Kim, H., Liu, S.-Y. & Taam, R. E. Templates of binary-induced spiral-shell patterns around mass-losing post-main-sequence stars. Astrophys. J. Suppl. Ser. 243, 35 (2019).

    ADS  CAS  Google Scholar 

  96. Pinte, C., Ménard, F., Duchêne, G. & Bastien, P. Monte Carlo radiative transfer in protoplanetary disks. Astron. Astrophys. 459, 797–804 (2006).

    ADS  Google Scholar 

  97. Pinte, C. et al. Benchmark problems for continuum radiative transfer. High optical depths, anisotropic scattering, and polarisation. Astron. Astrophys. 498, 967–980 (2009).

    ADS  Google Scholar 

  98. Weingartner, J. C. & Draine, B. T. Dust grain-size distributions and extinction in the Milky Way, Large Magellanic Cloud, and Small Magellanic Cloud. Astrophys. J. 548, 296–309 (2001).

    ADS  Google Scholar 

  99. Pinte, C. et al. Direct mapping of the temperature and velocity gradients in discs. Imaging the vertical CO snow line around IM Lupi. Astron. Astrophys. 609, 47 (2018).

    Google Scholar 

  100. Winnewisser, G., Belov, S. P., Klaus, T. & Schieder, R. Sub-Doppler measurements on the rotational transitions of carbon monoxide. J. Mol. Spectrosc. 184, 468–472 (1997).

    ADS  CAS  Google Scholar 

  101. Bogey, M., Demuynck, C. & Destombes, J. L. Laboratory measurement of the submillimeter wave spectrum of SiC and isotopomers. Astron. Astrophys. 232, 19 (1990).

    ADS  Google Scholar 

  102. Yamada, K. M. T., Moravec, A. & Winnewisser, G. Sub-millimeter wave spectra of cyanoacetylene and revised ground state constants. Z. Naturforsch. A 50, 1179–1181 (1995).

    ADS  CAS  Google Scholar 

  103. Müller, H. S. P. et al. Rotational spectroscopy of isotopologues of silicon monoxide, SiO, and spectroscopic parameters from a combined fit of rotational and rovibrational data. J. Phys. Chem. A 117, 13843–13854 (2013).

    PubMed  Google Scholar 

  104. Müller, H. S. P. et al. Rotational spectroscopy of the isotopic species of silicon monosulfide, SiS. Phys. Chem. Chem. Phys. 9, 1579 (2007).

    PubMed  Google Scholar 

  105. Gottlieb, C. A., Myers, P. C. & Thaddeus, P. Precise millimeter-wave laboratory frequencies for CS and C34S. Astrophys. J. 588, 655–661 (2003).

    ADS  CAS  Google Scholar 

  106. Ahrens, V. et al. Sub-Doppler saturation spectroscopy of HCN up to 1 THz and detection of J = 3 → 2 (4 → 3) emission from TMC1. Z. Naturforsch. A 57, 669–681 (2002).

    ADS  CAS  Google Scholar 

  107. Fuchs, U. et al. High resolution spectroscopy of HCN isotopomers: H13CN, HC15N, and H13C15N in the ground and first excited bending vibrational state. Z. Naturforsch. A 59, 861–872 (2004).

    ADS  CAS  Google Scholar 

  108. Bogey, M., Demuynck, C. & Destombes, J. L. The millimetre wave spectrum of the 13C14N radical in its ground state. Can. J. Phys. 62, 1248–1253 (1984).

    ADS  CAS  Google Scholar 

  109. Müller, H. S. P., Thorwirth, S., Roth, D. A. & Winnewisser, G. The Cologne Database for Molecular Spectroscopy, CDMS. Astron. Astrophys. 370, 49–52 (2001).

    ADS  Google Scholar 

  110. Müller, H. S. P., Schlöder, F., Stutzki, J. & Winnewisser, G. The Cologne Database for Molecular Spectroscopy, CDMS: a useful tool for astronomers and spectroscopists. J Mol. Struct. 742, 215–227 (2005).

    ADS  Google Scholar 

Download references

Acknowledgements

We thank S.-H. Cho of the Korean VLBI Network for KVN observations of W Aql to confirm consistency with our ALMA results. T.D. is supported in part by the Australian Research Council through a Discovery Early Career Researcher Award (DE230100183). T.D., F.D.C. and S.H.J.W. acknowledge support from the Research Foundation Flanders (FWO) through grants 12N9920N, 1253223N and 1285221N, respectively. J.M. and S.M. acknowledge support from the Research Foundation Flanders (FWO) grant G099720N. M.V.d.S. acknowledges support from European Union’s Horizon 2020 Research and Innovation programme under the Marie Skłodowska-Curie grant agreement number 882991. M.M. acknowledges funding form the Programme Paris Region fellowship supported by the Région Ile-de-France. P.K. acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation programme (synergy grant project UniverScale, grant agreement 951549). T.J.M. is grateful to the STFC for support through grant ST/P000312/1 and thanks the Leverhulme Trust for the award of an Emeritus Fellowship. J.M.C.P. was supported by STFC grant number ST/T000287/1. L.D., J.M.C.P., S.H.J.W., S.M. and D.G. acknowledge support from ERC consolidator grant 646758 AEROSOL. E.D.B. acknowledges support from the Swedish National Space Agency. D.G. was funded by the project grant ‘The Origin and Fate of Dust in Our Universe’ from the Knut and Alice Wallenberg Foundation. K.T.W. acknowledges support from the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation programme (grant agreement number 883867, project EXWINGS). F.H., A.B. and L.M. acknowledge funding from the French National Research Agency (ANR) project PEPPER (ANR-20-CE31-0002). H.S.P.M. acknowledges support by the Deutsche Forschungsgemeinschaft through the collaborative research grant SFB 956 (project ID 184018867). R.S.’s contribution to the research described here was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with NASA, and funded in part by NASA via ADAP awards, and multiple HST GO awards from the Space Telescope Science Institute. A.Z. is funded by STFC/UKRI through grant ST/T000414/1. This research was supported in part by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. This project has received funding from the Framework Program for Research and Innovation ‘Horizon 2020’ under the convention Marie Skłodowska-Curie No 945298. Computational resources and services used in this work were provided by the VSC (Flemish Supercomputer Center), funded by the Research Foundation Flanders (FWO) and the Flemish Government, department EWI. This research was undertaken with the assistance of resources and services from the National Computational Infrastructure (NCI), which is supported by the Australian Government. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2018.1.00659.L. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programme 0103.D-0772(A). We acknowledge excellent support from the UK ALMA Regional Centre (UK ARC), which is hosted by the Jodrell Bank Centre for Astrophysics (JBCA) at the University of Manchester. The UK ARC Node is supported by STFC Grant ST/P000827/1.

Author information

Authors and Affiliations

Authors

Contributions

T.D. conceived of and led this publication, analysed and interpreted data, performed the radiative transfer models, wrote the paper, and created most of the figures. J.M. performed and interpreted the hydrodynamics model and made Fig. 5c, and Supplementary Figs. 9a,b and 10. M.V.d.S. led the chemical interpretation and made Supplementary Figs. 7 and 8. M.M. and P.K. contributed the analysis of the resolved imaging. A.M.S.R. performed the ALMA data reduction. F.D.C. and A.C. contributed to the three-dimensional interpretation of the data. T.J.M. and J.M.C.P. contributed to the chemical interpretation. C.A.G. contributed to the line identifications and interpretation. C.P. assisted in the MCFOST modelling. D.J.P. assisted in the Phantom modelling and interpretation. E.D.B. contributed the fully reduced APEX data. The ALMA proposal was led by L.D. and C.A.G., with contributions from M.M., T.D., A.d.K., K.M.M., R.S., A.M.S.R., J.M.C.P., H.S.P.M., E.D.B., P.K., A.B., K.T.W., M.V.d.S., E.L., D.G., J.Y. and D.J.P. All authors contributed to discussions of the results, and commented on the paper and analysis.

Corresponding author

Correspondence to T. Danilovich.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Astronomy thanks the anonymous reviewers for their contribution to the peer review of this work.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data

Extended Data Fig. 1 A series of sketches illustrating the formation of SiN during the periastron passage of the W Aql system.

The orbit (black line) is shown face on in the frame of the AGB star and the F9 star is assumed to be moving clockwise. Relative to our observations, the observer is located to the left, represented by the radio dish. (a) The F9 star (yellow) approaches the AGB star (red) and enters the dense inner wind region (\({n}_{{{{{\rm{H}}}}}_{2}} \sim 1{0}^{8}\) to 1010 cm−3). (b) The rapid periastron passage is completed and SiN (and, similarly, SiC and NS) has formed in the wake of the F9 star (cyan region), with formation initiated by the F9 UV flux (see ‘Chemical modelling’ in the Methods). (c) As the F9 star continues on its orbit, the arc of SiN expands away from the AGB star, along with the stellar wind in which it is embedded. The present-day configuration of SiN can be seen in Fig. 1, where the PV diagram is a good approximation of the final arc shape that would be seen around the AGB star were the orbit viewed face-on.

Extended Data Fig. 2 Observed SiC emission towards W Aql.

(a) Zeroth moment map of SiC towards W Aql with contours at levels of 3 and 5σ. Transition details are given in Table 1. North is up and east is to the left. The position of the AGB star is indicated by the red star at (0,0) and the location of the F9 companion is indicated by the yellow star to the south-west. North is up and east is left. The white ellipse in the bottom left corner indicates the size of the synthesized beam. (b) Position-velocity diagram of SiC towards W Aql, taken with the same wide slit as used for SiN (Fig. 1), with a position angle of north 33 east. Dashed black contours are at levels of 3 and 5σ, a dotted white parabola is fit to the data (see ‘SiN and SiC’ in the ‘Data reduction’ section in the Methods), and a dash-dotted pink ellipse is plotted to emphasize the shape of the emission in the PV diagram. The position and LSR velocity of the AGB star is indicated by the red star and the horizontal yellow dotted line indicates the present offset of the F9 star.

Extended Data Fig. 3 Plots of NS towards W Aql.

(a) Zeroth moment map of NS towards W Aql with contours at levels of 3 and 5σ. North is up and east is to the left. The position of the AGB star is indicated by the red star at (0,0) and the location of the F9 companion is indicated by the yellow star to the south-west. The white ellipse in the bottom left corner indicates the size of the synthesized beam. (b) Position-velocity diagram of NS taken with the same wide slit as used for SiN (Fig. 1). The position and LSR velocity of the AGB star is indicated by the red star and the horizontal dotted yellow line indicates the present offset of the F9 star. (c) Spectra of the NS, SiN and SiC lines given in Table 1. All lines were extracted for circular apertures with radii 0.25, centred on the continuum peak. The flux of the SiC spectrum is multiplied by 5 to allow for a more direct comparison to SiN and NS.

Extended Data Fig. 4 Central channels of molecular emission towards W Aql.

A channel of SiS (left) and CS (right) are plotted, showing the asymmetric distribution of these molecules caused by the flux from the F9 star. The positions of the AGB and F9 stars are indicated by the red and yellow stars, respectively. The channel velocities are given in the top right corners and the beam is shown in the bottom left corner. Contours are plotted for levels of 3, 5, and 10σ. North is up and east is left.

Extended Data Fig. 5 A comparison of spectral lines at different locations in the AGB wind.

Plots are shown of CS, HCN, SiO and H13CN emission extracted from circular apertures with 100 mas radii centred on the F9 star (blue), on the AGB star (orange) and at the same separation as the F9 star but on the opposite side of the AGB (Opp. F9, brown, dashed). (See Table 1 for line frequencies.) The AGB and Opp. F9 line profiles are scaled by the factor given in the legend to facilitate comparison with the F9 line profiles. The vertical grey line indicates υLSR = − 23km s−1.

Extended Data Fig. 6 Zeroth moment maps of HC3N towards W Aql.

The transition for each map is given in the top right, with further details given in Table 1. Contours are plotted at levels of 3 and 5σ. North is up and east is to the left. The position of the AGB star is indicated by the red star at (0,0), also corresponding to the continuum peak, and the location of the F9 companion is indicated by the yellow star to the south-west. The white ellipse in the bottom left corner indicates the size of the synthesized beam.

Extended Data Fig. 7 CO emission towards W Aql.

Channel maps of CO (J = 2 → 1) are shown towards W Aql, obtained by combining observations from three configurations of ALMA. The AGB star is located at (0,0) and is marked by a red cross. The LSR velocity of each channel is given in the top right hand corner and the three channels closest to the W Aql υLSR = − 23km s−1 are highlighted with red borders and summed for Fig. 5. The synthetic beam is given by the white ellipse in the bottom left corner of each channel. North is up and east is left.

Extended Data Fig. 8 Variations in the radial emission of CO with angle.

The plots show the radial emission distribution against angle for the summed central three channels of CO (Fig. 5) with a full revolution shown in the centre (0 to 2π) and half a revolution is shown on either side ( − π to 0 and 2π to 3π) to show how the structures extend onwards. The location of the F9 star is indicated by the yellow star and a yellow dotted line which passes through both stars is plotted in the central winding to guide the eye. The black, red and white curves correspond to the same features highlighted in Fig. 5. The top plot shows the full observed extent of the CO emission (out to 15) and the bottom plot focuses on the regions out to 5 from the AGB star. These plots are reproduced without the additional curves in Supplementary Figure 4.

Extended Data Fig. 9 Observations and models of blue and red CO emission towards W Aql.

The plots show that blue (a and c) and red (b and d) channels equidistant from the stellar LSR velocity (υLSR = − 23km s−1) in velocity space do not exhibit identical CO emission patterns. The ALMA observations (a and b) show an elongated emission region on the blue side (a) and an approximately round emission region on the red side (b). The same pattern is mimicked in the red (c) and blue (d) channels of the hydrodynamic model processed with MCFOST. The red and black crosses correspond to the locations of the AGB and F9 stars. Note that the modelled and observed positions do not exactly correspond. Details are given in ‘Hydrodynamic simulations’ in the Methods.

Extended Data Table 1 Possible orbital solutions for the W Aql system

Supplementary information

Supplementary Information

Supplementary Discussion (‘Radiation pressure on dust’ and ‘Anisotropic mass loss’), Table 1 and Figs. 1–13.

Rights and permissions

Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Danilovich, T., Malfait, J., Van de Sande, M. et al. Chemical tracers of a highly eccentric AGB–main-sequence star binary. Nat Astron 8, 308–327 (2024). https://doi.org/10.1038/s41550-023-02154-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41550-023-02154-y

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing