Showing 1–50 of 120 results for author: Richards, A M S

Maser Flares Driven by Isothermal Shock Waves

Authors: M. D. Gray, S. Etoka, B. Pimpanuwat, A. M. S. Richards

Abstract: We use 3D computer modelling to investigate the timescales and radiative output from maser flares generated by the impact of shock-waves on astronomical unit-scale clouds in interstellar and star-forming regions, and in circumstellar regions in some circumstances. Physical conditions are derived from simple models of isothermal hydrodynamic (single-fluid) and C-type (ionic and neutral fluid) shock… ▽ More We use 3D computer modelling to investigate the timescales and radiative output from maser flares generated by the impact of shock-waves on astronomical unit-scale clouds in interstellar and star-forming regions, and in circumstellar regions in some circumstances. Physical conditions are derived from simple models of isothermal hydrodynamic (single-fluid) and C-type (ionic and neutral fluid) shock-waves, and based on the ortho-H$_2$O 22-GHz transition. Maser saturation is comprehensively included, and we find that the most saturated maser inversions are found predominantly in the shocked material. We study the effect on the intensity, flux density and duration of flares of the following parameters: the pre-shock level of saturation, the observer's viewpoint, and the shock speed. Our models are able to reproduce observed flare rise times of a few times 10 days, specific intensities of up to 10$^5$ times the saturation intensity and flux densities of order $100(R/d)^2$Jy from a source of radius $R$ astronomical units at a distance of $d$ kiloparsec. We found that flares from C-type shocks are approximately 5 times more likely to be seen by a randomly placed observer than flares from hydrodynamically shocked clouds of similar dimensions. We computed intrinsic beaming patterns of the maser emission, finding substantial extension of the pattern parallel to the shock front in the hydrodynamic models. Beaming solid angles for hydrodynamic models can be as small as $1.3\times 10^{-5}$sr, but are an order of magnitude larger for C-type models. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: Accepted for publication in MNRAS - 27 pages, 21 figures, 2 Tables

Deep JWST/NIRCam imaging of Supernova 1987A

Authors: Mikako Matsuura, M. Boyer, Richard G. Arendt, J. Larsson, C. Fransson, A. Rest, A. P. Ravi, S. Park, P. Cigan, T. Temim, E. Dwek, M. J. Barlow, P. Bouchet, G. Clayton, R. Chevalier, J. Danziger, J. De Buizer, I. De Looze, G. De Marchi, O. Fox, C. Gall, R. D. Gehrz, H. L. Gomez, R. Indebetouw, T. Kangas , et al. (24 additional authors not shown)

Abstract: JWST/NIRCam obtained high angular-resolution (0.05-0.1''), deep near-infrared 1--5 micron imaging of Supernova (SN) 1987A taken 35 years after the explosion. In the NIRCam images, we identify: 1) faint H2 crescents, which are emissions located between the ejecta and the equatorial ring, 2) a bar, which is a substructure of the ejecta, and 3) the bright 3-5 micron continuum emission exterior to the… ▽ More JWST/NIRCam obtained high angular-resolution (0.05-0.1''), deep near-infrared 1--5 micron imaging of Supernova (SN) 1987A taken 35 years after the explosion. In the NIRCam images, we identify: 1) faint H2 crescents, which are emissions located between the ejecta and the equatorial ring, 2) a bar, which is a substructure of the ejecta, and 3) the bright 3-5 micron continuum emission exterior to the equatorial ring. The emission of the remnant in the NIRCam 1-2.3 micron images is mostly due to line emission, which is mostly emitted in the ejecta and in the hot spots within the equatorial ring. In contrast, the NIRCam 3-5 micron images are dominated by continuum emission. In the ejecta, the continuum is due to dust, obscuring the centre of the ejecta. In contrast, in the ring and exterior to the ring, synchrotron emission contributes a substantial fraction to the continuum. Dust emission contributes to the continuum at outer spots and diffuse emission exterior to the ring, but little within the ring. This shows that dust cooling and destruction time scales are shorter than the synchrotron cooling time scale, and the time scale of hydrogen recombination in the ring is even longer than the synchrotron cooling time scale. With the advent of high sensitivity and high angular resolution images provided by JWST/NIRCam, our observations of SN 1987A demonstrate that NIRCam opens up a window to study particle-acceleration and shock physics in unprecedented details, probed by near-infrared synchrotron emission, building a precise picture of how a SN evolves. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: Accepted for publication in MNRAS; 18 pages

Detection of Rydberg lines from the atmosphere of Betelgeuse

Authors: W. R. F. Dent, G. Harper, A. M. S. Richards, P. Kervella, L. D. Matthews

Abstract: Emission lines from Rydberg transitions are detected for the first time from a region close to the surface of Betelgeuse. The H30$α$ line is observed at 231.905 GHz, with a FWHM ~42 km/s and extended wings. A second line at 232.025 GHz (FWHM ~21 km/s), is modeled as a combination of Rydberg transitions of abundant low First Ionization Potential metals. Both H30$α$ and the Rydberg combined line X30… ▽ More Emission lines from Rydberg transitions are detected for the first time from a region close to the surface of Betelgeuse. The H30$α$ line is observed at 231.905 GHz, with a FWHM ~42 km/s and extended wings. A second line at 232.025 GHz (FWHM ~21 km/s), is modeled as a combination of Rydberg transitions of abundant low First Ionization Potential metals. Both H30$α$ and the Rydberg combined line X30$α$ are fitted by Voigt profiles, and collisional broadening with electrons may be partly responsible for the Lorentzian contribution, indicating electron densities of a few 10$^8$cm$^{-3}$. X30$α$ is located in a relatively smooth ring at a projected radius of 0.9x the optical photospheric radius R$_*$, whereas H30$α$ is more clumpy, reaching a peak at ~1.4R$_*$. We use a semi-empirical thermodynamic atmospheric model of Betelgeuse to compute the 232 GHz (1.29mm) continuum and line profiles making simple assumptions. Photoionized abundant metals dominate the electron density and the predicted surface of continuum optical depth unity at 232 GHz occurs at ~1.3R$_*$, in good agreement with observations. Assuming a Saha-Boltzmann distribution for the level populations of Mg, Si, and Fe, the model predicts that the X30$α$ emission arises in a region of radially-increasing temperature and turbulence. Inclusion of ionized C and non-LTE effects could modify the integrated fluxes and location of emission. These simulations confirm the identity of the Rydberg transition lines observed towards Betelgeuse, and reveal that such diagnostics can improve future atmospheric models. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: ApJLett (accepted). 4 figures

Resolved ALMA observations of water in the inner astronomical units of the HL Tau disk

Authors: Stefano Facchini, Leonardo Testi, Elizabeth Humphreys, Mathieu Vander Donckt, Andrea Isella, Ramon Wrzosek, Alain Baudry, Malcom D. Gray, Anita M. S. Richards, Wouter Vlemmmings

Abstract: The water molecule is a key ingredient in the formation of planetary systems, with the water snowline being a favourable location for the growth of massive planetary cores. Here we present Atacama Large Millimeter/ submillimeter Array data of the ringed protoplanetary disk orbiting the young star HL Tauri that show centrally peaked, bright emission arising from three distinct transitions of the ma… ▽ More The water molecule is a key ingredient in the formation of planetary systems, with the water snowline being a favourable location for the growth of massive planetary cores. Here we present Atacama Large Millimeter/ submillimeter Array data of the ringed protoplanetary disk orbiting the young star HL Tauri that show centrally peaked, bright emission arising from three distinct transitions of the main water isotopologue. The spatially and spectrally resolved water content probes gas in a thermal range down to the water sublimation temperature. Our analysis implies a stringent lower limit of 3.7 Earth oceans of water vapour available within the inner 17 astronomical units of the system. We show that our observations are limited to probing the water content in the atmosphere of the disk, due to the high dust column density and absorption, and indicate that the main water isotopologue is the best tracer to spatially resolve water vapour in protoplanetary disks. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: Published in Nature Astronomy (2024), arXiv version presents text and figures before proofing

The Hidden Clumps in VY CMa Uncovered by ALMA

Authors: Roberta M. Humphreys, A. M. S. Richards, Kris Davidson, A. P. Singh, L. Decin, L. M. Ziurys

Abstract: The red hypergiant VY CMa is famous for its very visible record of high mass loss events. Recent CO observations with ALMA revealed three previously unknown large scale outflows (Paper I). In this paper we use the CO maps to investigate the motions of a cluster of four clumps close to the star, not visible in the optical or infrared images. We present their proper motions measured from two epochs… ▽ More The red hypergiant VY CMa is famous for its very visible record of high mass loss events. Recent CO observations with ALMA revealed three previously unknown large scale outflows (Paper I). In this paper we use the CO maps to investigate the motions of a cluster of four clumps close to the star, not visible in the optical or infrared images. We present their proper motions measured from two epochs of ALMA images and determine the line of sight velocities of the gas in emission at the clumps. We estimate their masses and ages, or time since ejection, and conclude that all four were ejected during VY CMa's active period in the early 20th century. Together with two additional knots observed with HST, VY CMa experienced at least six massive outflows during a 30 year period with a total mass lost greater than 0.07 Msun. The position-velocity map of the $^{12}$CO emission reveals previously unnoticed attributes of the older outer ejecta. In a very narrow range of Doppler velocities, $^{12}$CO absorption and emission causes some of this outer material to be quite opaque. At those frequencies the inner structure is hidden and we see only emission from an extended outer region. This fact produces a conspicuous but illusory dark spot if one attempts to subtract the continuum in a normal way. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Comments: To appear in the Astronomical Journal

ATOMIUM: Molecular inventory of 17 oxygen-rich evolved stars observed with ALMA

Authors: S. H. J. Wallstrom, T. Danilovich, H. S. P. Muller, C. A. Gottlieb, S. Maes, M. Van de Sande, L. Decin, A. M. S. Richards, A. Baudry, J. Bolte, T. Ceulemans, F. De Ceuster, A. de Koter, I. El Mellah, M. Esseldeurs, S. Etoka, D. Gobrecht, E. Gottlieb, M. Gray, F. Herpin, M. Jeste, D. Kee, P. Kervella, T. Khouri, E. Lagadec , et al. (13 additional authors not shown)

Abstract: The dusty winds of cool evolved stars are a major contributor of the newly synthesised material enriching the Galaxy and future generations of stars. However, the details of the physics and chemistry behind dust formation and wind launching have yet to be pinpointed. Recent spatially resolved observations show the importance of gaining a more comprehensive view of the circumstellar chemistry, but… ▽ More The dusty winds of cool evolved stars are a major contributor of the newly synthesised material enriching the Galaxy and future generations of stars. However, the details of the physics and chemistry behind dust formation and wind launching have yet to be pinpointed. Recent spatially resolved observations show the importance of gaining a more comprehensive view of the circumstellar chemistry, but a comparative study of the intricate interplay between chemistry and physics is still difficult because observational details such as frequencies and angular resolutions are rarely comparable. Aiming to overcome these deficiencies, ATOMIUM is an ALMA Large Programme to study the physics and chemistry of the circumstellar envelopes of a diverse set of oxygen-rich evolved stars under homogeneous observing conditions at three angular resolutions between ~0.02"-1.4". Here we summarize the molecular inventory of these sources, and the correlations between stellar parameters and molecular content. Seventeen oxygen-rich or S-type asymptotic giant branch (AGB) and red supergiant (RSG) stars have been observed in several tunings with ALMA Band 6, targeting a range of molecules to probe the circumstellar envelope and especially the chemistry of dust formation close to the star. We systematically assigned the molecular carriers of the spectral lines and measured their spectroscopic parameters and the angular extent of the emission of each line from integrated intensity maps. Across the ATOMIUM sample, we detect 291 transitions of 24 different molecules and their isotopologues. This includes several first detections in oxygen-rich AGB/RSG stars: PO v=1, SO2 v1=1 and v2=2, and several high energy H2O transitions. We also find several first detections in S-type AGB stars: vibrationally excited HCN v2=2,3 and SiS v=4,5,6, as well as first detections of the molecules SiC, AlCl, and AlF in W Aql... △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: 19 pages plus appendices, forthcoming publication in Astronomy & Astrophysics

ALMA Observations of Supernova Remnant N49 in the Large Magellanic Cloud. II. Non-LTE Analysis of Shock-heated Molecular Clouds

Authors: H. Sano, Y. Yamane, J. Th. van Loon, K. Furuya, Y. Fukui, R. Z. E. Alsaberi, A. Bamba, R. Enokiya, M. D. Filipović, R. Indebetouw, T. Inoue, A. Kawamura, M. Lakićević, C. J. Law, N. Mizuno, T. Murase, T. Onishi, S. Park, P. P. Plucinsky, J. Rho, A. M. S. Richards, G. Rowell, M. Sasaki, J. Seok, P. Sharda , et al. (6 additional authors not shown)

Abstract: We present the first compelling evidence of shock-heated molecular clouds associated with the supernova remnant (SNR) N49 in the Large Magellanic Cloud (LMC). Using $^{12}$CO($J$ = 2-1, 3-2) and $^{13}$CO($J$ = 2-1) line emission data taken with the Atacama Large Millimeter/Submillimeter Array, we derived the H$_2$ number density and kinetic temperature of eight $^{13}$CO-detected clouds using the… ▽ More We present the first compelling evidence of shock-heated molecular clouds associated with the supernova remnant (SNR) N49 in the Large Magellanic Cloud (LMC). Using $^{12}$CO($J$ = 2-1, 3-2) and $^{13}$CO($J$ = 2-1) line emission data taken with the Atacama Large Millimeter/Submillimeter Array, we derived the H$_2$ number density and kinetic temperature of eight $^{13}$CO-detected clouds using the large velocity gradient approximation at a resolution of 3.5$''$ (~0.8 pc at the LMC distance). The physical properties of the clouds are divided into two categories: three of them near the shock front show the highest temperatures of ~50 K with densities of ~500-700 cm$^{-3}$, while other clouds slightly distant from the SNR have moderate temperatures of ~20 K with densities of ~800-1300 cm$^{-3}$. The former clouds were heated by supernova shocks, but the latter were dominantly affected by the cosmic-ray heating. These findings are consistent with the efficient production of X-ray recombining plasma in N49 due to thermal conduction between the cold clouds and hot plasma. We also find that the gas pressure is roughly constant except for the three shock-engulfed clouds inside or on the SNR shell, suggesting that almost no clouds have evaporated within the short SNR age of ~4800 yr. This result is compatible with the shock-interaction model with dense and clumpy clouds inside a low-density wind bubble. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: 12 pages, 5 figures, 1 table, accepted for publication in The Astrophysical Journal (ApJ)

ALMA High-frequency Long Baseline Campaign in 2021: Highest Angular Resolution Submillimeter Wave Images for the Carbon-rich Star R Lep

Authors: Yoshiharu Asaki, Luke T. Maud, Harold Francke, Hiroshi Nagai, Dirk Petry, Edward B. Fomalont, Elizabeth Humphreys, Anita M. S. Richards, Ka Tat Wong, William Dent, Akihiko Hirota, Jose Miguel Fernandez, Satoko Takahashi, Antonio S. Hales

Abstract: The Atacama Large Millimeter/submillimeter Array (ALMA) was used in 2021 to image the carbon-rich evolved star R Lep in Bands 8-10 (397-908 GHz) with baselines up to 16 km. The goal was to validate the calibration, using band-to-band (B2B) phase referencing with a close phase calibrator J0504-1512, 1.2 deg from R Lep in this case, and the imaging procedures required to obtain the maximum angular r… ▽ More The Atacama Large Millimeter/submillimeter Array (ALMA) was used in 2021 to image the carbon-rich evolved star R Lep in Bands 8-10 (397-908 GHz) with baselines up to 16 km. The goal was to validate the calibration, using band-to-band (B2B) phase referencing with a close phase calibrator J0504-1512, 1.2 deg from R Lep in this case, and the imaging procedures required to obtain the maximum angular resolution achievable with ALMA. Images of the continuum emission and the hydrogen cyanide (HCN) maser line at 890.8 GHz, from the J=10-9 transition between the (1110) and (0400) vibrationally excited states, achieved angular resolutions of 13, 6, and 5 mas in Bands 8-10, respectively. Self-calibration (self-cal) was used to produce ideal images as to compare with the B2B phase referencing technique. The continuum emission was resolved in Bands 9 and 10, leaving too little flux for self-cal of the longest baselines, so these comparisons are made at coarser resolution. Comparisons showed that B2B phase referencing provided phase corrections sufficient to recover 92%, 83%, and 77% of the ideal image continuum flux densities. The HCN maser was sufficiently compact to obtain self-cal solutions in Band 10 for all baselines (up to 16 km). In Band 10, B2B phase referencing as compared to the ideal images recovered 61% and 70% of the flux density for the HCN maser and continuum, respectively. △ Less

Submitted 18 October, 2023; v1 submitted 14 October, 2023; originally announced October 2023.

Comments: 37 pages, 12 figures, 9 tables, accepted by ApJ (Aug 30, 2023)

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)

ALMA detection of CO rotational line emission in red supergiant stars of the massive young star cluster RSGC1 -- Determination of a new mass-loss rate prescription for red supergiants

Authors: Leen Decin, Anita M. S. Richards, Pablo Marchant, Hugues Sana

Abstract: [Abridged] Aim: We aim to derive a new mass-loss rate prescription for RSGs that is not afflicted with some uncertainties inherent in preceding studies. Methods: We have observed CO rotational line emission towards a sample of RSGs in the open cluster RSGC1 that all are of similar initial mass. The ALMA CO(2-1) line detections allow to retrieve the gas mass-loss rates (Mdot_CO). In contrast to mas… ▽ More [Abridged] Aim: We aim to derive a new mass-loss rate prescription for RSGs that is not afflicted with some uncertainties inherent in preceding studies. Methods: We have observed CO rotational line emission towards a sample of RSGs in the open cluster RSGC1 that all are of similar initial mass. The ALMA CO(2-1) line detections allow to retrieve the gas mass-loss rates (Mdot_CO). In contrast to mass-loss rates derived from the analysis of dust spectral features (Mdot_SED), the data allow a direct determination of the wind velocity and no uncertain dust-to-gas correction factor is needed. Results: Five RSGs in RSGC1 have been detected in CO(2-1). The retrieved Mdot_CO values are systematically lower than Mdot_SED. Although only five RSGs in RSGC1 have been detected, the data allow to propose a new mass-loss rate relation for M-type red supergiants that is dependent on luminosity and initial mass. The new mass-loss rate relation is based on the new Mdot_CO values for the RSGs in RSGC1 and on prior Mdot_SED values for RSGs in 4 clusters, including RSGC1. The new Mdot-prescription yields a good prediction for the mass-loss rate of some well-known Galactic RSGs that are observed in multiple CO rotational lines, including alpha Ori, mu Cep and VX Sgr. However, there are indications that a stronger, potentially eruptive, mass-loss process - different than captured by our new mass-loss rate prescription - is occurring during some fraction of the RSG lifetime. Implementing a lower mass-loss rate in evolution codes for massive stars has important consequences for the nature of their end-state. A reduction of the RSG mass-loss rate implies that quiescent RSG mass loss is not enough to strip a single star's hydrogen-rich envelope. Upon core-collapse such single stars would explode as RSG. △ Less

Submitted 13 November, 2023; v1 submitted 16 March, 2023; originally announced March 2023.

Comments: 26 pages, submitted to Astronomy & Astrophysics, minor revisions already implemented

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

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)

The Co-Ordinated Radio and Infrared Survey for High-Mass Star Formation. V. The CORNISH-South Survey and Catalogue

Authors: T. Irabor, M. G. Hoare, M. Burton, W. D. Cotton, P. Diamond, S. Dougherty, S. P. Ellingsen, R. Fender, G. A. Fuller, S. Garrington, P. F. Goldsmith, J. Green, A. G. Gunn, J. Jackson, S. Kurtz, S. L. Lumsden, J. Marti, I. McDonald, S. Molinari, T. J. Moore, M. Mutale, T. Muxlow, T. OBrien, R. D. Oudmaijer, R. Paladini , et al. (10 additional authors not shown)

Abstract: We present the first high spatial resolution radio continuum survey of the southern Galactic plane. The CORNISH project has mapped the region defined by $295^{\circ} < l < 350^{\circ}$; $|b| < 1^{\circ}$ at 5.5-GHz, with a resolution of 2.5$^{''}$ (FWHM). As with the CORNISH-North survey, this is designed to primarily provide matching radio data to the Spitzer GLIMPSE survey region. The CORNISH-So… ▽ More We present the first high spatial resolution radio continuum survey of the southern Galactic plane. The CORNISH project has mapped the region defined by $295^{\circ} < l < 350^{\circ}$; $|b| < 1^{\circ}$ at 5.5-GHz, with a resolution of 2.5$^{''}$ (FWHM). As with the CORNISH-North survey, this is designed to primarily provide matching radio data to the Spitzer GLIMPSE survey region. The CORNISH-South survey achieved a root mean square noise level of $\sim$ 0.11 mJy beam$^{-1}$, using the 6A configuration of the Australia Telescope Compact Array (ATCA). In this paper, we discuss the observations, data processing and measurements of the source properties. Above a 7$σ$ detection limit, 4701 sources were detected, and their ensemble properties show similar distributions with their northern counterparts. The catalogue is highly reliable and is complete to 90 per cent at a flux density level of 1.1 mJy. We developed a new way of measuring the integrated flux densities and angular sizes of non-Gaussian sources. The catalogue primarily provides positions, flux density measurements and angular sizes. All sources with IR counterparts at 8$μm$ have been visually classified, utilizing additional imaging data from optical, near-IR, mid-IR, far-IR and sub-millimetre galactic plane surveys. This has resulted in the detection of 524 H II regions of which 255 are ultra-compact H II regions, 287 planetary nebulae, 79 radio stars and 6 massive young stellar objects. The rest of the sources are likely to be extra-galactic. These data are particularly important in the characterization and population studies of compact ionized sources such as UCHII regions and PNe towards the Galactic mid-plane. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Characterising the high-mass star forming region IRAS 18144-1723 through methanol maser observations

Authors: Esraa Khafagy, khaled A. Edris, Osama M. Shalabiea, Anna Bartkiewicz, Anita M. S. Richards, Zainab Awad

Abstract: We introduce a study of the massive star forming region IRAS 18144--1723 using observations of the 6.7 GHz methanol maser line. Such regions are opaque at short wavelengths but can be observed through radio emission lines. In this study we traced the kinematics of the source on milliarcsecond scales using the Multi-Element-Radio-Interferometer-Network (MERLIN). We found 52 maser spots in the LSR v… ▽ More We introduce a study of the massive star forming region IRAS 18144--1723 using observations of the 6.7 GHz methanol maser line. Such regions are opaque at short wavelengths but can be observed through radio emission lines. In this study we traced the kinematics of the source on milliarcsecond scales using the Multi-Element-Radio-Interferometer-Network (MERLIN). We found 52 maser spots in the LSR velocity range 45--52 km s$^{-1}$, near the centre of the previously-detected CO range of 21.3--71.3 km s$^{-1}$, lying within $\sim$ 0$''$.5 of IRAS 18144--1723 `B', thought to be a young Class I protostar. Their distribution can be approximated as an ellipse, which, if it were rotating, would have its axis oriented south-east to north-west. The most probable morphology of the emitting regions is interaction between a disc and an outflow, possibly with a very large opening angle. The arcmin-scale CO outflow centred on source `B' is oriented East-West, and the methanol masers do show the highest dispersion of velocity gradients in approximately this direction, so the kinematics are complex and suggest that more than one source may be responsible. We also tested kinematic models for a Keplerian disc or a simple bipolar outflow, but neither are compatible with the kinematics of the maser clumps and the characteristics of their internal velocities. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Comments: 11 pages, 9 figures, 3 Tables, Accepted by MNRAS

Recovering Phosphine in Venus' Atmosphere from SOFIA Observations

Authors: Jane S. Greaves, Janusz J. Petkowski, Anita M. S. Richards, Clara Sousa-Silva, Sara Seager, David L. Clements

Abstract: Searches for phosphine in Venus' atmosphere have sparked a debate. Cordiner et al. 2022 analyse spectra from the Stratospheric Observatory For Infrared Astronomy (SOFIA) and infer <0.8 ppb of PH3. We noticed that some spectral artefacts arose from non-essential calibration-load signals. By-passing these signals allows simpler post-processing and a 5.7σ candidate detection, suggesting approx. 3 ppb… ▽ More Searches for phosphine in Venus' atmosphere have sparked a debate. Cordiner et al. 2022 analyse spectra from the Stratospheric Observatory For Infrared Astronomy (SOFIA) and infer <0.8 ppb of PH3. We noticed that some spectral artefacts arose from non-essential calibration-load signals. By-passing these signals allows simpler post-processing and a 5.7σ candidate detection, suggesting approx. 3 ppb of PH3 above the clouds. Compiling six phosphine results hints at an inverted abundance trend: decreasing above the clouds but rising again in the mesosphere from some unexplained source. However, no such extra source is needed if phosphine is undergoing destruction by sunlight (photolysis), to a similar degree as on Earth. Low phosphine values/limits are found where the viewed part of the super-rotating Venusian atmosphere had passed through sunlight, while high values are from views moving into sunlight. We suggest Venusian phosphine is indeed present, and so merits further work on models of its origins. △ Less

Submitted 2 March, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: v2 is an updated and revised version of the paper

Journal ref: Geophysical Research Letters, 50, e2023GL103539

Sciences with Thai National Radio Telescope

Authors: Phrudth Jaroenjittichai, Koichiro Sugiyama, Busaba H. Kramer, Boonrucksar Soonthornthum, Takuya Akahori, Kitiyanee Asanok, Willem Baan, Sherin Hassan Bran, Shari L. Breen, Se-Hyung Cho, Thanapol Chanapote, Richard Dodson, Simon P. Ellingsen, Sandra Etoka, Malcolm D. Gray, James A. Green, Kazuhiro Hada, Marcus Halson, Tomoya Hirota, Mareki Honma, Hiroshi Imai, Simon Johnston, Kee-Tae Kim, Michael Kramer, Di Li , et al. (22 additional authors not shown)

Abstract: This White Paper summarises potential key science topics to be achieved with Thai National Radio Telescope (TNRT). The commissioning phase has started in mid 2022. The key science topics consist of "Pulsars and Fast Radio Bursts (FRBs)", "Star Forming Regions (SFRs)", "Galaxy and Active Galactic Nuclei (AGNs)", "Evolved Stars", "Radio Emission of Chemically Peculiar (CP) Stars", and "Geodesy", cov… ▽ More This White Paper summarises potential key science topics to be achieved with Thai National Radio Telescope (TNRT). The commissioning phase has started in mid 2022. The key science topics consist of "Pulsars and Fast Radio Bursts (FRBs)", "Star Forming Regions (SFRs)", "Galaxy and Active Galactic Nuclei (AGNs)", "Evolved Stars", "Radio Emission of Chemically Peculiar (CP) Stars", and "Geodesy", covering a wide range of observing frequencies in L/C/X/Ku/K/Q/W-bands (1-115 GHz). As a single-dish instrument, TNRT is a perfect tool to explore time domain astronomy with its agile observing systems and flexible operation. Due to its ideal geographical location, TNRT will significantly enhance Very Long Baseline Interferometry (VLBI) arrays, such as East Asian VLBI Network (EAVN), Australia Long Baseline Array (LBA), European VLBI Network (EVN), in particular via providing a unique coverage of the sky resulting in a better complete "uv" coverage, improving synthesized-beam and imaging quality with reducing side-lobes. This document highlights key science topics achievable with TNRT in single-dish mode and in collaboration with VLBI arrays. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Comments: 75 pages, 16 figures, 15 tables. White Paper for Potential Key Sciences to be Achieved with Thai National Radio Telescope (TNRT)

Stringent limits on $^{28}$SiO maser emission from the recurrent nova T Coronae Borealis

Authors: A. Evans, B. Pimpanuwat, A. M. S. Richards, D. P. K. Banerjee, U. Munari, M. D. Gray, B. Hutawarakorn Kramer, A. Kraus

Abstract: There are indications that the third known eruption of the recurrent nova T CrB is imminent, and multi-wavelength observations prior to the eruption are important to characterise the system before it erupts. T CrB is known to display the SiO fundamental vibrational feature at 8$\,μ$m. When the anticipated eruption occurs, it is possible that the shock produced when the ejected material runs into t… ▽ More There are indications that the third known eruption of the recurrent nova T CrB is imminent, and multi-wavelength observations prior to the eruption are important to characterise the system before it erupts. T CrB is known to display the SiO fundamental vibrational feature at 8$\,μ$m. When the anticipated eruption occurs, it is possible that the shock produced when the ejected material runs into the wind of the red giant in the system may be traced using SiO maser emission. We have used the 100m Effelsberg Radio Telescope to search for $^{28}$SiO emission in the $\upsilon=1$, $\upsilon=2$, $J=1\rightarrow0$ transitions, at 43.122 GHz and 42.820~GHz respectively, while the system is in quiescence. We find no evidence for such emission. △ Less

Submitted 29 July, 2022; originally announced July 2022.

Comments: Published in MNRAS, 514, 4894 - 4897 (2022)

Journal ref: MNRAS, 514, 4894 - 4897 (2022)

Self-calibration and improving image fidelity for ALMA and other radio interferometers

Authors: A. M. S. Richards, E. Moravec, S. Etoka, E. B. Fomalont, A. F. Pérez-Sánchez, M. C. Toribio, R. A. Laing

Abstract: This manual is intended to help ALMA and other interferometer users improve images by recognising limitations and how to overcome them and deciding when and how to use self-calibration. The images provided by the ALMA Science Archive are calibrated using standard observing and data processing routines, including a quality assurance process to make sure that the observations meet the proposer's sci… ▽ More This manual is intended to help ALMA and other interferometer users improve images by recognising limitations and how to overcome them and deciding when and how to use self-calibration. The images provided by the ALMA Science Archive are calibrated using standard observing and data processing routines, including a quality assurance process to make sure that the observations meet the proposer's science requirements. This may not represent the full potential of the data, since any interferometry observation can be imaged with a range of resolutions and surface brightness sensitivity. The separation between phase calibration source and target usually limits the target dynamic range to a few hundred (or 50--100 for challenging conditions) but if the noise in the target field has not reached the thermal limit, improvements may be possible using self-calibration. This often requires judgements based on the target properties and is not yet automated for all situations. This manual provides background on the instrumental and atmospheric causes of visibility phase and amplitude errors, their effects on imaging and how to improve the signal to noise ratio and image fidelity by self-calibration. We introduce the conditions for self-calibration to be useful and how to estimate calibration parameter values for a range of observing modes (continuum, spectral line etc.). We also summarise more general error recognition and other techniques to tackle imaging problems. The examples are drawn from ALMA interferometric data processed using CASA, but the principles are generally applicable to most similar cm to sub-mm imaging. △ Less

Submitted 12 July, 2022; originally announced July 2022.

Comments: 76 pages, 55 figures, ALMA Memo series

Report number: 620

A MeerKAT, e-MERLIN, H.E.S.S. and Swift search for persistent and transient emission associated with three localised FRBs

Authors: James O. Chibueze, M. Caleb, L. Spitler, H. Ashkar, F. Schussler, B. W. Stappers, C. Venter, I. Heywood, A. M. S. Richards, D. R. A. Williams, M. Kramer, R. Beswick, M. C. Bezuidenhout, R. P. Breton, L. N. Driessen, F. Jankowski, E. F. Keane, M. Malenta, M. Mickaliger, V. Morello, H. Qiu, K. Rajwade, S. Sanidas, M. Surnis, T. W. Scragg , et al. (134 additional authors not shown)

Abstract: We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, an… ▽ More We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, and X-ray bands. For this FRB, we obtain a UV flux upper limit of 1.39x10^-16 erg/cm^-2/s/Amstrong, X-ray limit of ~ 6.6x10^-14 erg/cm^-2/s and a limit on the very-high-energy gamma-ray flux (Phi) (E > 120 GeV) < 1.7 x 10^-12 erg/cm^-2/s. We obtain a radio upper limit of ~15 microJy/beam for persistent emission at the locations of both FRBs 20190711A and 20171019A, but detect diffuse radio emission with a peak brightness of ~53 microJy/beam associated with FRB 20190714A at z = 0.2365. This represents the first detection of the radio continuum emission potentially associated with the host (galaxy) of FRB 20190714A, and is only the third known FRB to have such an association. Given the possible association of a faint persistent source, FRB 20190714A may potentially be a repeating FRB whose age lies between that of FRB 20121102A and FRB 20180916A. A parallel search for repeat bursts from these FRBs revealed no new detections down to a fluence of 0.08 Jy ms for a 1 ms duration burst. △ Less

Submitted 31 December, 2021; originally announced January 2022.

Comments: 15 pages, 9 figures

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)

Venusian phosphine: a 'Wow!' signal in chemistry?

Authors: William Bains, Janusz J. Petkowski, Sara Seager, Sukrit Ranjan, Clara Sousa-Silva, Paul B. Rimmer, Zhuchang Zhan, Jane S. Greaves, Anita M. S. Richards

Abstract: The potential detection of ppb levels phosphine (PH3) in the clouds of Venus through millimeter-wavelength astronomical observations is extremely surprising as PH3 is an unexpected component of an oxidized environment of Venus. A thorough analysis of potential sources suggests that no known process in the consensus model of Venus' atmosphere or geology could produce PH3 at anywhere near the observ… ▽ More The potential detection of ppb levels phosphine (PH3) in the clouds of Venus through millimeter-wavelength astronomical observations is extremely surprising as PH3 is an unexpected component of an oxidized environment of Venus. A thorough analysis of potential sources suggests that no known process in the consensus model of Venus' atmosphere or geology could produce PH3 at anywhere near the observed abundance. Therefore, if the presence of PH3 in Venus' atmosphere is confirmed, it is highly likely to be the result of a process not previously considered plausible for Venusian conditions. The source of atmospheric PH3 could be unknown geo- or photochemistry, which would imply that the consensus on Venus' chemistry is significantly incomplete. An even more extreme possibility is that strictly aerial microbial biosphere produces PH3. This paper summarizes the Venusian PH3 discovery and the scientific debate that arose since the original candidate detection one year ago. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Comments: A short overview of the Venusian PH3 discovery and the scientific debate that arose since the original candidate detection in September 2020. Additional discussion of possible non-canonical sources of PH3 on Venus is also included. arXiv admin note: text overlap with arXiv:2009.06499

Journal ref: GPSS

The Nearby Evolved Stars Survey II: Constructing a volume-limited sample and first results from the James Clerk Maxwell Telescope

Authors: P. Scicluna, F. Kemper, I. McDonald, S. Srinivasan, A. Trejo, S. H. J. Wallström, J. G. A. Wouterloot, J. Cami, J. Greaves, Jinhua He, D. T. Hoai, Hyosun Kim, O. C. Jones, H. Shinnaga, C. J. R. Clark, T. Dharmawardena, W. Holland, H. Imai, J. Th. van Loon, K. M. Menten, R. Wesson, H. Chawner, S. Feng, S. Goldman, F. C. Liu , et al. (67 additional authors not shown)

Abstract: The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a n… ▽ More The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a new metric for estimating the distances to evolved stars and compare its results to \emph{Gaia} EDR3. Replicating other studies, the most-evolved, highly enshrouded objects in the Galactic Plane dominate the dust returned by our sources, and we initially estimate a total DPR of $4.7\times 10^{-5}$ M$_\odot$ yr$^{-1}$ from our sample. Our sub-mm fluxes are systematically higher and spectral indices are typically shallower than dust models typically predict. The 450/850 $μ$m spectral indices are consistent with the blackbody Rayleigh--Jeans regime, suggesting a large fraction of evolved stars have unexpectedly large envelopes of cold dust. △ Less

Submitted 24 October, 2021; originally announced October 2021.

Comments: 20 pages, 13 figures, 5 tables, accepted for publications in MNRAS

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)

Low levels of sulphur dioxide contamination of Venusian phosphine spectra

Authors: Jane S. Greaves, Paul B. Rimmer, Anita M. S. Richards, Janusz J. Petkowski, William Bains, Sukrit Ranjan, Sara Seager, David L. Clements, Clara Sousa Silva, Helen J. Fraser

Abstract: New analysis is presented of the 1.1 mm wavelength absorption lines in Venus' atmosphere that suggested the presence of phosphine. We retrieve a sulphur dioxide observation from the JCMT archive that was simultaneous within a few days of the PH3 1-0 spectrum obtained in June 2017, and demonstrate via a radiative transfer calculation that contamination of PH3 by SO2 was ~10 per cent. We also presen… ▽ More New analysis is presented of the 1.1 mm wavelength absorption lines in Venus' atmosphere that suggested the presence of phosphine. We retrieve a sulphur dioxide observation from the JCMT archive that was simultaneous within a few days of the PH3 1-0 spectrum obtained in June 2017, and demonstrate via a radiative transfer calculation that contamination of PH3 by SO2 was ~10 per cent. We also present ALMA 2019 spectra of PH3 1-0 and an SO2 transition acquired simultaneously, and infer that SO2 line-contamination was ~2 percent (for the least-noisy half of the planetary disc). The contamination-subtracted ALMA and JCMT spectra (of 6-8 sigma confidence) are now consistent with similar absorption-depths at the two epochs. The two values span -1.9(+/-0.2) 10-4 of the continuum signal (which was re-estimated for ALMA), albeit for differing planetary areas. This suggests that the abundance attributed to phosphine in Venus' atmosphere was broadly similar in 2017 and 2019. △ Less

Submitted 21 June, 2022; v1 submitted 18 August, 2021; originally announced August 2021.

Comments: v3 is an proofread version of the paper that will appear in print in MNRAS

Journal ref: MNRAS 2022

The Additional Representative Images for Legacy (ARI-L) project for the ALMA Science Archive

Authors: M. Massardi, F. Stoehr, G. J. Bendo, M. Bonato, J. Brand, V. Galluzzi, F. Guglielmetti, E. Liuzzo, N. Marchili, A. M. S. Richards, K. L. J. Rygl, F. Bedosti, A. Giannetti, M. Stagni, C. Knapic, M. Sponza, G. A. Fuller, T. W. B. Muxlow

Abstract: The Additional Representative Images for Legacy (ARI-L) project is a European Development project for ALMA Upgrade approved by the Joint ALMA Observatory (JAO) and the European Southern Observatory (ESO), started in June 2019. It aims to increase the legacy value of the ALMA Science Archive (ASA) by bringing the reduction level of ALMA data from Cycles 2-4 close to that of data from more recent Cy… ▽ More The Additional Representative Images for Legacy (ARI-L) project is a European Development project for ALMA Upgrade approved by the Joint ALMA Observatory (JAO) and the European Southern Observatory (ESO), started in June 2019. It aims to increase the legacy value of the ALMA Science Archive (ASA) by bringing the reduction level of ALMA data from Cycles 2-4 close to that of data from more recent Cycles processed for imaging with the ALMA Pipeline. As of mid-2021 more than 150000 images have been returned to the ASA for public use. At its completion in 2022, the project will have provided enhanced products for at least 70% of the observational data from Cycles 2-4 processable with the ALMA Pipeline. In this paper we present the project rationale, its implementation, and the new opportunities offered to ASA users by the ARI-L products. The ARI-L cubes and images complement the much limited number of archival image products generated during the data quality assurance stages (QA2), which cover only a small fraction of the available data for those Cycles. ARI-L imaging products are highly relevant for many science cases and significantly enhance the possibilities for exploiting archival data. Indeed, ARI-L products facilitate archive access and data usage for science purposes even for non-expert data miners, provide a homogeneous view of all data for better dataset comparisons and download selections, make the archive more accessible to visualization and analysis tools, and enable the generation of preview images and plots similar to those possible for subsequent Cycles. △ Less

Submitted 23 July, 2021; originally announced July 2021.

Comments: 15 pages. Accepted for publication in PASP

SiO masers and 7 mm Continuum in Mira and R Aqr

Authors: W. D. Cotton, E. Humphreys, M. Wittkowski, A. Baudry, A. M. S. Richards, W. Vlemmings, T. Khouri, S. Etoka

Abstract: Interactions between AGB stars and a secondary in a close orbit are one possible explanation of why some AGB stars develop into aspherical planetary nebulae. This study uses millimeter observations of the continuum and SiO maser emission in several symbiotic Miras looking for evidence of an interaction between the two stars. New JVLA observations at ~45 mas resolution are analyzed, imaging continu… ▽ More Interactions between AGB stars and a secondary in a close orbit are one possible explanation of why some AGB stars develop into aspherical planetary nebulae. This study uses millimeter observations of the continuum and SiO maser emission in several symbiotic Miras looking for evidence of an interaction between the two stars. New JVLA observations at ~45 mas resolution are analyzed, imaging continuum and SiO masers. Two of the three targets were detected and accurately registered continuum and line images were derived. No clear evidence of an interaction was found between components B and A in Mira. R Aqr has a well known jet arising from the secondary star. The jet may be disturbing the circumstellar envelop of the AGB star or possibly just nearly aligned with it. △ Less

Submitted 20 July, 2021; originally announced July 2021.

Comments: 7 pages, 7 figures

Methanol and water maser observations separate disc and outflow sources in IRAS 19410+2336

Authors: M. S. Darwish, K. A. Edris, A. M. S. Richards, S. Etoka, M. S. Saad, M. M. Beheary, G. A Fuller

Abstract: We investigate the kinematics of high mass protostellar objects within the high mass star forming region IRAS 19410+2336. We performed high angular resolution observations of 6.7-GHz methanol and 22 GHz water masers using the MERLIN (Multi-Element Radio Linked Interferometer Network) and e-MERLIN interferometers. The 6.7-GHz methanol maser emission line was detected within the $\sim$ 16--27 km s… ▽ More We investigate the kinematics of high mass protostellar objects within the high mass star forming region IRAS 19410+2336. We performed high angular resolution observations of 6.7-GHz methanol and 22 GHz water masers using the MERLIN (Multi-Element Radio Linked Interferometer Network) and e-MERLIN interferometers. The 6.7-GHz methanol maser emission line was detected within the $\sim$ 16--27 km s$^{-1}$ velocity range with a peak flux density $\sim$50 Jy. The maser spots are spread over $\sim$1.3 arcsec on the sky, corresponding to $\sim$2800 au at a distance of 2.16 kpc. These are the first astrometric measurements at 6.7 GHz in IRAS 19410+2336. The 22-GHz water maser line was imaged in 2005 and 2019 (the latter with good astrometry). Its velocities range from 13 to $\sim$29 km s$^{-1}$. The peak flux density was found to be 18.7 Jy and 13.487 Jy in 2005, and 2019, respectively. The distribution of the water maser components is up to 165 mas, $\sim$350 au at 2.16 kpc. We find that the Eastern methanol masers most probably trace outflows from the region of millimetre source mm1. The water masers to the West lie in a disc (flared or interacting with outflow/infall) around another more evolved millimetre source (13-s). The maser distribution suggests that the disc lies at an angle of 60$^{\circ}$ or more to the plane of the sky and the observed line of sight velocities then suggest an enclosed mass between 44 M$_{\odot}$ and as little as 11 M$_{\odot}$ if the disc is edge-on. The Western methanol masers may be infalling. △ Less

Submitted 24 May, 2021; originally announced May 2021.

Journal ref: 2020MNRAS.493.4442D

Recovery of Spectra of Phosphine in Venus' Clouds

Authors: Jane S. Greaves, Anita M. S. Richards, William Bains, Paul B. Rimmer, David L. Clements, Sara Seager, Janusz J. Petkowski, Clara Sousa-Silva, Sukrit Ranjan, Helen J. Fraser

Abstract: We recover PH3 in the atmosphere of Venus in data taken with ALMA, using three different calibration methods. The whole-planet signal is recovered with 5.4σ confidence using Venus bandpass self-calibration, and two simpler approaches are shown to yield example 4.5-4.8σ detections of the equatorial belt. Non-recovery by Villanueva et al. is attributable to (a) including areas of the planet with hig… ▽ More We recover PH3 in the atmosphere of Venus in data taken with ALMA, using three different calibration methods. The whole-planet signal is recovered with 5.4σ confidence using Venus bandpass self-calibration, and two simpler approaches are shown to yield example 4.5-4.8σ detections of the equatorial belt. Non-recovery by Villanueva et al. is attributable to (a) including areas of the planet with high spectral-artefacts and (b) retaining all antenna baselines which raises the noise by a factor ~2.5. We release a data-processing script that enables our whole-planet result to be reproduced. The JCMT detection of PH3 remains robust, with the alternative SO2 attribution proposed by Villanueva et al. appearing inconsistent both in line-velocity and with millimetre-wavelength SO2 monitoring. SO2 contamination of the ALMA PH3-line is minimal. Net abundances for PH3, in the gas column above ~55 km, are up to ~20 ppb planet-wide with JCMT, and ~7 ppb with ALMA (but with signal-loss possible on scales approaching planetary size). Derived abundances will differ if PH3 occupies restricted altitudes - molecules in the clouds will contribute significantly less absorption at line-centre than equivalent numbers of mesospheric molecules - but in the latter zone, PH3 lifetime is expected to be short. Given we recover phosphine, we suggest possible solutions (requiring substantial further testing): a small collisional broadening coefficient could give narrow lines from lower altitude, or a high eddy diffusion coefficient could allow molecules to survive longer at higher altitudes. Alternatively, PH3 could be actively produced by an unknown mechanism in the mesosphere, but this would need to be in addition to cloud-level PH3 detected retrospectively by Pioneer-Venus. △ Less

Submitted 19 April, 2021; originally announced April 2021.

Comments: Revised and updated re-analysis of the phosphine spectra in Venus' clouds; Submitted to Nature Astronomy "Matters Arising" on Apr/19/2021

Journal ref: Nature Astronomy 5.7 (2021): 636-639

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)

OH maser toward IRAS 06056+2131: polarization parameters and evolution status

Authors: Mohamed Said darwish, Anita M. S. Richards, Sandra Etoka, Khaled Edris, Somaya Saad, Mohamed Beheary, Gary Fullar

Abstract: We present high angular resolution observations of OH maser emission towards the high-mass star forming region IRAS 06056+2131. The observations were carried out using the UK radio interferometer array, Multi-Element Radio Linked Interferometer Network (MERLIN) in the OH main lines at 1665- and 1667-MHz, in addition to the OH satellite line at 1720-MHz. The results of this study revealed emission… ▽ More We present high angular resolution observations of OH maser emission towards the high-mass star forming region IRAS 06056+2131. The observations were carried out using the UK radio interferometer array, Multi-Element Radio Linked Interferometer Network (MERLIN) in the OH main lines at 1665- and 1667-MHz, in addition to the OH satellite line at 1720-MHz. The results of this study revealed emission in the 1665 MHz line with an estimated total intensity of $\sim 4$ Jy. We did not detect any emission from the 1667-MHz and 1720-MHz lines. The full polarization mode of MERLIN enables us to investigate the magnetic field in the OH maser region. Our results show that IRAS 06056+2131 is a highly circularly polarized source. In this transition, a Zeeman pair is identified from which a magnetic strength of $\sim -1.5$ mG is inferred. The orientation of the linear polarization vectors suggests that the magnetic field lines at the location of the OH maser emission \textbf{might be} in agreement with the orientation of the outflow thought to be associated with this source. The star forming evolutionary status of the embedded proto-stellar object is discussed. △ Less

Submitted 17 January, 2021; originally announced January 2021.

Comments: 10 pages, 5 figures

Journal ref: 2020MNRAS.499.1441D

The Mass-Loss History of the Red Hypergiant VY CMa

Authors: Roberta M. Humphreys, Kris Davidson, A. M. S. Richards, L. M. Ziurys, Terry J. Jones, Kazunori Ishibashi

Abstract: Imaging and spectroscopy of the knots, clumps, and extended arcs in the complex ejecta of VY CMa confirm a record of high mass loss events over the past few hundred years. HST/STIS spectroscopy of numerous small knots close to the star allow us to measure their radial velocities from the strong K I emission and determine their separate motions, spatial orientations, and time since ejecta. Their ag… ▽ More Imaging and spectroscopy of the knots, clumps, and extended arcs in the complex ejecta of VY CMa confirm a record of high mass loss events over the past few hundred years. HST/STIS spectroscopy of numerous small knots close to the star allow us to measure their radial velocities from the strong K I emission and determine their separate motions, spatial orientations, and time since ejecta. Their ages concentrate around 70, 120, 200 and 250 years ago. A K I emission knot only 50 mas from the star ejected as recently as 1985 -- 1995 may coincide with an H2O maser. Comparison with VY CMa's historic light curve from 1800 to the present, shows several knots with ejection times that correspond with extended periods of variability and deep minima. The similarity of this correspondence in VY CMa with the remarkable recent dimming of Betelgeuse and an outflow of gas is apparent. The evidence for similar outflows from the surface of a more typical red supergiant suggests that discrete ejections are more common and surface or convective activity is a major source of mass loss for red supergiants. △ Less

Submitted 15 December, 2020; originally announced December 2020.

Comments: Accepted for publication in the Astronomical Journal

Planet formation in intermediate-separation binary systems

Authors: O. Panić, T. J. Haworth, M. G. Petr-Gotzens, J. Miley, M. van den Ancker, M. Vioque, L. Siess, R. Parker, C. J. Clarke, I. Kamp, G. Kennedy, R. D. Oudmaijer, I. Pascucci, A. M. S. Richards, T. Ratzka, C. Qi

Abstract: We report the first characterisation of the individual discs in the intermediate separation binary systems KK Oph and HD 144668 at millimetre wavelengths. In both systems the circum-primary and the circum-secondary discs are detected in the millimetre continuum emission, but not in $^{13}$CO nor C$^{18}$O lines. Even though the disc structure is only marginally resolved, we find indications of lar… ▽ More We report the first characterisation of the individual discs in the intermediate separation binary systems KK Oph and HD 144668 at millimetre wavelengths. In both systems the circum-primary and the circum-secondary discs are detected in the millimetre continuum emission, but not in $^{13}$CO nor C$^{18}$O lines. Even though the disc structure is only marginally resolved, we find indications of large-scale asymmetries in the outer regions of the primary discs, most likely due to perturbation by the companion. The derived dust masses are firmly above debris disc level for all stars. The primaries have about three times more dust in their discs than the secondaries. In the case of HD 144668 the opacity spectral index of the primary and secondary differ by the large margin of 0.69 which may be a consequence of the secondary disc being more compact. Upper limits on the gas masses imply less than 0.1 M$_{\textrm{jup}}$ in any of these discs, meaning that giant planets can no longer form in them. Considering that there have been no massive gas discs identified to date in intermediate separation binaries (i.e., binaries at a few hundred au separation), this opens space for speculation whether their binarity causes the removal of gas, with tidal interaction truncating the discs and hence shortening the accretion timescale. More systematic studies in this respect are sorely needed. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: 12 pages. Accepted for publication in MNRAS

Re-analysis of Phosphine in Venus' Clouds

Authors: Jane S. Greaves, Anita M. S. Richards, William Bains, Paul B. Rimmer, David L. Clements, Sara Seager, Janusz J. Petkowski, Clara Sousa-Silva, Sukrit Ranjan, Helen J. Fraser

Abstract: We first respond to two points raised by Villanueva et al. We show the JCMT discovery spectrum of PH3 can not be re-attributed to SO2, as the line width is larger than observed for SO2 features, and the required abundance would be an extreme outlier. The JCMT spectrum is also consistent with our simple model, constant PH3-abundance with altitude, with no discrepancy in line profile (within data li… ▽ More We first respond to two points raised by Villanueva et al. We show the JCMT discovery spectrum of PH3 can not be re-attributed to SO2, as the line width is larger than observed for SO2 features, and the required abundance would be an extreme outlier. The JCMT spectrum is also consistent with our simple model, constant PH3-abundance with altitude, with no discrepancy in line profile (within data limits); reconciliation with a full photochemical model is the subject of future work. Section 2 presents initial results from re-processed ALMA data. Villanueva et al. noted an issue with bandpass calibration. They have worked on a partially re-processed subset of the ALMA data, so we note where their conclusions, and those of Greaves et al., are now superseded. To summarise: we recover PH3 in Venus' atmosphere with ALMA (~5σ confidence). Localised abundance appears to peak at ~5-10 parts-per-billion (ppb), with suggestions of spatial variation. Advanced data-products suggest a planet-averaged PH3 abundance ~1-4 ppb, lower than from the earlier ALMA processing (which indicated 7+ ppb). The ALMA data are reconcilable with the JCMT detection (~20 ppb) if there is order-of-magnitude temporal variation; more advanced processing of the JCMT data is underway to check methods. Independent PH3 measurements suggest possible altitude dependence (under ~5 ppb at 60+ km, up to ~100 ppb at 50+ km; see Section 2: Conclusions.). Given that both ALMA and JCMT were working at the limit of observatory capabilities, new spectra should be obtained. The ALMA data in-hand are no longer limited by calibration, but spectral ripples still exist, probably due to size and brightness of Venus in relation to the primary beam. Further, spatial ripples are present, potentially reducing significance of real narrow spectral features. △ Less

Submitted 10 December, 2020; v1 submitted 16 November, 2020; originally announced November 2020.

Comments: Response to the Matters Arising: "No phosphine in the atmosphere of Venus" by G. Villanueva et al., preprint at arXiv:2010.14305; v2 contains updated values and additional technical details

Journal ref: Nature Astronomy 5.7 (2021): 636-639

Rotational spectra of vibrationally excited AlO and TiO in oxygen rich stars

Authors: T. Danilovich, C. A. Gottlieb, L. Decin, A. M. S. Richards, K. L. K. Lee, T. Kaminski, N. A. Patel, K. H. Young, K. M. Menten

Abstract: Rotational transitions in vibrationally excited AlO and TiO -- two possible precursors of dust -- were observed in the 300 GHz range (1 mm wavelength) towards the oxygen rich AGB stars R Dor and IK Tau with ALMA, and vibrationally excited AlO was observed towards the red supergiant VY CMa with the SMA. The $J=11 \to 10$ transition of TiO in the $v=1~{\rm{and}}~2$ levels, and the $N = 9 \to 8$ tran… ▽ More Rotational transitions in vibrationally excited AlO and TiO -- two possible precursors of dust -- were observed in the 300 GHz range (1 mm wavelength) towards the oxygen rich AGB stars R Dor and IK Tau with ALMA, and vibrationally excited AlO was observed towards the red supergiant VY CMa with the SMA. The $J=11 \to 10$ transition of TiO in the $v=1~{\rm{and}}~2$ levels, and the $N = 9 \to 8$ transition in the $v=2$ level of AlO were identified towards R Dor; the $J=11 \to 10$ line of TiO was identified in the $v=1$ level towards IK Tau; and two transitions in the $v=1~{\rm{and}}~2$ levels of AlO were identified towards VY CMa. The newly-derived high vibrational temperature of TiO and AlO in R Dor of $1800 \pm 200$ K, and prior measurements of the angular extent confirm that the majority of the emission is from a region within $\lesssim2R_{\star}$ of the central star. A full radiative transfer analysis of AlO in R Dor yielded a fractional abundance of $\sim$3% of the solar abundance of Al. From a similar analysis of TiO a fractional abundance of $\sim78$% of the solar abundance of Ti was found. The observations provide indirect evidence that TiO is present in a rotating disk close to the star. Further observations in the ground and excited vibrational levels are needed to determine whether AlO, TiO, and TiO$_2$ are seeds of the Al$_2$O$_3$ dust in R Dor, and perhaps in the gravitationally bound dust shells in other AGB stars with low mass loss rates. △ Less

Submitted 13 October, 2020; originally announced October 2020.

Comments: Accepted for publication in ApJ

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

Phosphine Gas in the Cloud Decks of Venus

Authors: Jane S. Greaves, Anita M. S. Richards, William Bains, Paul B. Rimmer, Hideo Sagawa, David L. Clements, Sara Seager, Janusz J. Petkowski, Clara Sousa-Silva, Sukrit Ranjan, Emily Drabek-Maunder, Helen J. Fraser, Annabel Cartwright, Ingo Mueller-Wodarg, Zhuchang Zhan, Per Friberg, Iain Coulson, E'lisa Lee, Jim Hoge

Abstract: Measurements of trace-gases in planetary atmospheres help us explore chemical conditions different to those on Earth. Our nearest neighbor, Venus, has cloud decks that are temperate but hyper-acidic. We report the apparent presence of phosphine (PH3) gas in Venusian atmosphere, where any phosphorus should be in oxidized forms. Single-line millimeter-waveband spectral detections (quality up to ~15… ▽ More Measurements of trace-gases in planetary atmospheres help us explore chemical conditions different to those on Earth. Our nearest neighbor, Venus, has cloud decks that are temperate but hyper-acidic. We report the apparent presence of phosphine (PH3) gas in Venusian atmosphere, where any phosphorus should be in oxidized forms. Single-line millimeter-waveband spectral detections (quality up to ~15 sigma) from the JCMT and ALMA telescopes have no other plausible identification. Atmospheric PH3 at ~20 parts-per-billion abundance is inferred. The presence of phosphine is unexplained after exhaustive study of steady-state chemistry and photochemical pathways, with no currently-known abiotic production routes in Venusian atmosphere, clouds, surface and subsurface, or from lightning, volcanic or meteoritic delivery. Phosphine could originate from unknown photochemistry or geochemistry, or, by analogy with biological production of phosphine on Earth, from the presence of life. Other PH3 spectral features should be sought, while in-situ cloud and surface sampling could examine sources of this gas. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Journal ref: Greaves, J. S., Richards, A. M., Bains, W., Rimmer, P. B., Sagawa, H., Clements, D. L., ... & Hoge, J. (2020). Phosphine gas in the cloud decks of Venus. Nature Astronomy, 1-10

Phosphine on Venus Cannot be Explained by Conventional Processes

Authors: William Bains, Janusz J. Petkowski, Sara Seager, Sukrit Ranjan, Clara Sousa-Silva, Paul B. Rimmer, Zhuchang Zhan, Jane S. Greaves, Anita M. S. Richards

Abstract: The recent candidate detection of ~1 ppb of phosphine in the middle atmosphere of Venus is so unexpected that it requires an exhaustive search for explanations of its origin. Phosphorus-containing species have not been modelled for Venus' atmosphere before and our work represents the first attempt to model phosphorus species in the Venusian atmosphere. We thoroughly explore the potential pathways… ▽ More The recent candidate detection of ~1 ppb of phosphine in the middle atmosphere of Venus is so unexpected that it requires an exhaustive search for explanations of its origin. Phosphorus-containing species have not been modelled for Venus' atmosphere before and our work represents the first attempt to model phosphorus species in the Venusian atmosphere. We thoroughly explore the potential pathways of formation of phosphine in a Venusian environment, including in the planet's atmosphere, cloud and haze layers, surface, and subsurface. We investigate gas reactions, geochemical reactions, photochemistry, and other non-equilibrium processes. None of these potential phosphine production pathways are sufficient to explain the presence of ppb phosphine levels on Venus. If PH3's presence in Venus' atmosphere is confirmed, it therefore is highly likely to be the result of a process not previously considered plausible for Venusian conditions. The process could be unknown geochemistry, photochemistry, or even aerial microbial life, given that on Earth phosphine is exclusively associated with anthropogenic and biological sources. The detection of phosphine adds to the complexity of chemical processes in the Venusian environment and motivates in situ follow up sampling missions to Venus. Our analysis provides a template for investigation of phosphine as a biosignature on other worlds. △ Less

Submitted 12 July, 2021; v1 submitted 14 September, 2020; originally announced September 2020.

Comments: v2 is in press in Astrobiology, Special Collection: Venus; v2 also expands on the potential of phosphides from the deep mantle volcanism as a source of PH3 (as suggested by Truong and Lunine 2021: https://www.pnas.org/content/118/29/e2021689118) and shows the volcanic source of PH3 to be unlikely

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)

Spatially Resolved Ultraviolet Spectroscopy of the Great Dimming of Betelgeuse

Authors: Andrea K. Dupree, Klaus G. Strassmeier, Lynn D. Matthews, Han Uitenbroek, Thomas Calderwood, Thomas Granzer, Edward F Guinan, Reimar Leike, Miguel Montargès, Anita M. S. Richards, Richard Wasatonic, Michael Weber

Abstract: The bright supergiant, Betelgeuse (Alpha Orionis, HD 39801) experienced a visual dimming during 2019 December and the first quarter of 2020 reaching an historic minimum 2020 February 7$-$13. During 2019 September-November, prior to the optical dimming event, the photosphere was expanding. At the same time, spatially resolved ultraviolet spectra using the Hubble Space Telescope/Space Telescope Imag… ▽ More The bright supergiant, Betelgeuse (Alpha Orionis, HD 39801) experienced a visual dimming during 2019 December and the first quarter of 2020 reaching an historic minimum 2020 February 7$-$13. During 2019 September-November, prior to the optical dimming event, the photosphere was expanding. At the same time, spatially resolved ultraviolet spectra using the Hubble Space Telescope/Space Telescope Imaging Spectrograph revealed a substantial increase in the ultraviolet spectrum and Mg II line emission from the chromosphere over the southern hemisphere of the star. Moreover, the temperature and electron density inferred from the spectrum and C II diagnostics also increased in this hemisphere. These changes happened prior to the Great Dimming Event. Variations in the Mg II k-line profiles suggest material moved outwards in response to the passage of a pulse or acoustic shock from 2019 September through 2019 November. It appears that this extraordinary outflow of material from the star, likely initiated by convective photospheric elements, was enhanced by the coincidence with the outward motions in this phase of the $\sim$400 day pulsation cycle. These ultraviolet observations appear to provide the connecting link between the known large convective cells in the photosphere and the mass ejection event that cooled to form the dust cloud in the southern hemisphere imaged in 2019 December, and led to the exceptional optical dimming of Betelgeuse in 2020 February. △ Less

Submitted 11 August, 2020; originally announced August 2020.

Comments: 11 pages, 8 figures, Astrophysical Journal, accepted

ALMA and VLA reveal the lukewarm chromospheres of the nearby red supergiants Antares and Betelgeuse

Authors: E. O'Gorman, G. M. Harper, K. Ohnaka, A. Feeney-Johansson, K. Wilkeneit-Braun, A. Brown, E. F. Guinan, J. Lim, A. M. S. Richards, N. Ryde, W. H. T. Vlemmings

Abstract: We first present spatially resolved ALMA and VLA continuum observations of the early-M red supergiant Antares to search for the presence of a chromosphere at radio wavelengths. We resolve the free-free emission of the Antares atmosphere at 11 unique wavelengths between 0.7 mm (ALMA band 8) and 10 cm (VLA S band). The projected angular diameter is found to continually increase with increasing wavel… ▽ More We first present spatially resolved ALMA and VLA continuum observations of the early-M red supergiant Antares to search for the presence of a chromosphere at radio wavelengths. We resolve the free-free emission of the Antares atmosphere at 11 unique wavelengths between 0.7 mm (ALMA band 8) and 10 cm (VLA S band). The projected angular diameter is found to continually increase with increasing wavelength, from a low of 50.7 mas at 0.7 mm up to a diameter of 431 mas at 10 cm, which corresponds to 1.35 and 11.6 times the photospheric angular diameter, respectively. All four ALMA measurements show that the shape of the atmosphere is elongated, with a flattening of 15% at a similar position angle. The disk-averaged gas temperature of the atmosphere initially rises from a value of 2700 K at 1.35 $R_{\star}$ (i.e., 0.35 $R_{\star}$ above the photosphere) to a peak value of 3800 K at $\sim$2.5 $R_{\star}$, after which it then more gradually decreases to 1650 K at 11.6 $R_{\star}$. The rise in gas temperature between 1.35 $R_{\star}$ and $\sim$2.5 $R_{\star}$ is evidence for a chromospheric temperature rise above the photosphere of a red supergiant. We detect a clear change in the spectral index across the sampled wavelength range, with the flux density ${S_ν} \propto ν^{1.42}$ between 0.7 mm and 1.4 cm, which we associate with chromosphere-dominated emission, while the flux density ${S_ν} \propto ν^{0.8}$ between 4.3 cm and 10 cm, which we associate with wind-dominated emission. We then perform nonlocal thermal equilibrium modeling of the far-ultraviolet radiation field of another early-M red supergiant, Betelgeuse, and find that an additional hot (i.e., $>7000$ K) chromospheric photoionization component with a much smaller filling factor must also exist throughout the chromospheres of these stars. △ Less

Submitted 14 June, 2020; originally announced June 2020.

Comments: 12 pages, 6 figures, accepted for publication in A&A

The e-MERLIN Galaxy Evolution Survey (e-MERGE): Overview and Survey Description

Authors: T. W. B. Muxlow, A. P. Thomson, J. F. Radcliffe, N. H. Wrigley, R. J. Beswick, Ian Smail, I. M. McHardy, S. T. Garrington, R. J. Ivison, M. J. Jarvis, I. Prandoni, M. Bondi, D. Guidetti, M. K. Argo, David Bacon, P. N. Best, A. D. Biggs, S. C. Chapman, K. Coppin, H. Chen, T. K. Garratt, M. A. Garrett, E. Ibar, Jean-Paul Kneib, Kirsten K. Knudsen , et al. (15 additional authors not shown)

Abstract: We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively clo… ▽ More We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5μ$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1. △ Less

Submitted 5 May, 2020; originally announced May 2020.

Comments: 21 pages, 12 colour figures, 3 tables -- accepted for publication in MNRAS

ALMA High-frequency Long Baseline Campaign in 2017: Band-to-band Phase Referencing in Submillimeter Waves

Authors: Yoshiharu Asaki, Luke T. Maud, Edward B. Fomalont, Neil M. Phillips, Akihiko Hirota, Tsuyoshi Sawada, Loreto Barcos-Muñoz, Anita M. S. Richards, William R. F. Dent, Satoko Takahashi, Stuartt Corder, John M. Carpenter, Eric Villard, Elizabeth M. Humphreys

Abstract: In 2017, an Atacama Large Millimeter/submillimeter Array (ALMA) high-frequency long baseline campaign was organized to test image capabilities with baselines up to 16 km at submillimeter (submm) wavelengths. We investigated image qualities using ALMA receiver Bands 7, 8, 9, and 10 (285-875 GHz) by adopting band-to-band (B2B) phase referencing in which a phase calibrator is tracked at a lower frequ… ▽ More In 2017, an Atacama Large Millimeter/submillimeter Array (ALMA) high-frequency long baseline campaign was organized to test image capabilities with baselines up to 16 km at submillimeter (submm) wavelengths. We investigated image qualities using ALMA receiver Bands 7, 8, 9, and 10 (285-875 GHz) by adopting band-to-band (B2B) phase referencing in which a phase calibrator is tracked at a lower frequency. For B2B phase referencing, it is expected that a closer phase calibrator to a target can be used, comparing to standard in-band phase referencing. In the first step, it is ensured that an instrumental phase offset difference between low- and high-frequency Bands can be removed using a differential gain calibration in which a phase calibrator is certainly detected while frequency switching. In the next step, comparative experiments are arranged to investigate the image quality between B2B and in-band phase referencing with phase calibrators at various separation angles. In the final step, we conducted long baseline imaging tests for a quasar at 289 GHz in Band 7 and 405 GHz in Band 8 and complex structure sources of HL Tau and VY CMa at ~670 GHz in Band 9. The B2B phase referencing was successfully applied, allowing us to achieve an angular resolution of 14x11 and 10x8 mas for HL Tau and VY CMa, respectively. There is a high probability of finding a low-frequency calibrator within 5.4 deg in B2B phase referencing, bright enough to use an 8 s scan length combined with a 7.5 GHz bandwidth. △ Less

Submitted 23 March, 2020; v1 submitted 16 March, 2020; originally announced March 2020.

Comments: 61 pages, 17 figures, 8 tables

Journal ref: 2020ApJS..247...23A

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

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

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

Submitted 9 March, 2020; originally announced March 2020.

Comments: Accepted for publication in MNRAS

COBRaS: The e-MERLIN 21 cm Legacy survey of Cygnus OB2

Authors: J. C. Morford, D. M. Fenech, R. K. Prinja, R. Blomme, J. A. Yates, J. J. Drake, S. P. S. Eyres, A. M. S. Richards, I. R. Stevens, N. J. Wright, J. S. Clark, S. Dougherty, J. M. Pittard, H. Smith, J. S. Vink

Abstract: The role of massive stars is central to an understanding of galactic ecology. It is important to establish the details of how massive stars provide radiative, chemical, and mechanical feedback in galaxies. Central to these issues is an understanding of the evolution of massive stars, and the critical role of mass loss via strongly structured winds and stellar binarity. Ultimately, massive stellar… ▽ More The role of massive stars is central to an understanding of galactic ecology. It is important to establish the details of how massive stars provide radiative, chemical, and mechanical feedback in galaxies. Central to these issues is an understanding of the evolution of massive stars, and the critical role of mass loss via strongly structured winds and stellar binarity. Ultimately, massive stellar clusters shape the structure and energetics of galaxies. We aim to conduct high-resolution, deep field mapping at 21cm of the core of the massive Cygnus OB2 association and to characterise the properties of the massive stars and colliding winds at this waveband. We used seven stations of the e-MERLIN radio facility, with its upgraded bandwidth and enhanced sensitivity to conduct a 21cm census of Cygnus OB2. Based on 42 hours of observations, seven overlapping pointings were employed over multiple epochs during 2014 resulting in 1 sigma sensitivities down to ~21microJy and a resolution of ~180mas. A total of 61 sources are detected at 21cm over a ~0.48deg x 0.48deg region centred on the heart of the Cyg OB2 association. Of these 61 sources, 33 are detected for the first time. We detect a number of previously identified sources including four massive stellar binary systems, two YSOs, and several known X-ray and radio sources. We also detect the LBV candidate (possible binary system) and blue hypergiant (BHG) star of Cyg OB2 #12. The 21cm observations secured in the COBRaS Legacy project provide data to constrain conditions in the outer wind regions of massive stars; determine the non-thermal properties of massive interacting binaries; examine evidence for transient sources, including those associated with young stellar objects; and provide unidentified sources that merit follow-up observations. The 21cm data are of lasting value and will serve in combination with other key surveys of Cyg OB2. △ Less

Submitted 20 March, 2020; v1 submitted 17 January, 2020; originally announced January 2020.

Comments: 41 pages, 12 figures, accepted in A&A

Journal ref: A&A 637, A64 (2020)

Circumstellar CO J = 3-2 detected around the evolving metal-poor ([Fe/H] ~ -1.15 dex) AGB star RU Vulpeculae

Authors: Iain McDonald, Stefan Uttenthaler, Albert A. Zijlstra, Anita M. S. Richards, Eric Lagadec

Abstract: We report the first detection of CO J = 3-2 around a truly metal-poor evolved star. RU Vulpeculae is modelled to have Teff ~ 3620 K, L ~ 3128 +/- 516 L_Sun, log(g) = 0.0 ~ 0.2 dex and [Fe/H] = -1.3 to -1.0 dex, and is modelled to have recently undergone a thermal pulse. Its infrared flux has approximately doubled over 35 years. ALMA observations show the 3-2 line is narrow (half-width ~1.8-3.5 km/… ▽ More We report the first detection of CO J = 3-2 around a truly metal-poor evolved star. RU Vulpeculae is modelled to have Teff ~ 3620 K, L ~ 3128 +/- 516 L_Sun, log(g) = 0.0 ~ 0.2 dex and [Fe/H] = -1.3 to -1.0 dex, and is modelled to have recently undergone a thermal pulse. Its infrared flux has approximately doubled over 35 years. ALMA observations show the 3-2 line is narrow (half-width ~1.8-3.5 km/s). The 2-1 line is much weaker: it is not confidently detected. Spectral-energy-distribution fitting indicates very little circumstellar absorption, despite its substantial mid-infrared emission. A VISIR mid-infrared spectrum shows features typical of previously observed metal-poor stars, dominated by a substantial infrared excess but with weak silicate and (possibly) Al2O3 emission. A lack of resolved emission, combined with weak 2!1 emission, indicates the dense circumstellar material is truncated at large radii. We suggest that rapid dust condensation is occurring, but with an aspherical geometry (e.g., a disc or clumps) that does not obscure the star. We compare with T UMi, a similar star which is currently losing its dust. △ Less

Submitted 23 October, 2019; originally announced October 2019.

Comments: Accepted MNRAS, 16 pages, 7 figures

The Unexpected Spectrum of the Innermost Ejecta of the Red Hypergiant VY CMa

Authors: Roberta M. Humphreys, L. M. Ziurys, J. J. Bernal, Michael S. Gordon, L. Andrew Helton, Kazunori Ishibashi, Terry J. Jones, A. M. S. Richards, Wouter Vlemmings

Abstract: HST/STIS spectra of the small clumps and filaments closest to the central star in VY CMa reveal that the very strong K I emission and TiO and VO molecular emission, long thought to form in a dusty circumstellar shell, actually originate in a few small clumps 100's of AU from the star. The K I lines are 10 to 20 times stronger in these nearest ejecta than on the star. The observations also confirm… ▽ More HST/STIS spectra of the small clumps and filaments closest to the central star in VY CMa reveal that the very strong K I emission and TiO and VO molecular emission, long thought to form in a dusty circumstellar shell, actually originate in a few small clumps 100's of AU from the star. The K I lines are 10 to 20 times stronger in these nearest ejecta than on the star. The observations also confirm VO as a circumstellar molecule. In this letter we discuss the spectra of the features, their motions and ages, and the identification of the molecular emission. The strength of the atomic and molecular features in the small clumps present an astrophysical problem for the excitation process. We show that the clumps must have a nearly clear line of sight to the star's radiation. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: To appear in the Astrophysical Journal

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)

Asymmetric mid-plane gas in ALMA images of HD~100546

Authors: J. Miley, O. Panić, T. J. Haworth, I. Pascucci, M. Wyatt, C. Clarke, A. M. S. Richards, T. Ratzka

Abstract: In this paper we present new ALMA observations towards the proto-planet hosting transitional disc of Herbig Ae/Be star HD 100546. This includes resolved 1.3 mm continuum, $^{13}$CO and the first detection of C$^{18}$O in this disc, which displays azimuthal asymmetry in regions spatially coincident with structures previously identified in HST images related to spiral arms. The lower limit on the ma… ▽ More In this paper we present new ALMA observations towards the proto-planet hosting transitional disc of Herbig Ae/Be star HD 100546. This includes resolved 1.3 mm continuum, $^{13}$CO and the first detection of C$^{18}$O in this disc, which displays azimuthal asymmetry in regions spatially coincident with structures previously identified in HST images related to spiral arms. The lower limit on the mass of the dust disc is calculated to be 9.6x10$^{-4}$M$_\odot$. A firm lower-limit on the total gas mass calculated from optically thin, mid-plane tracing C$^{18}$O (2-1) emission is 0.018M$_\odot$ assuming ISM abundances. These mass estimates provide an estimate of gas-to-dust ratio in the disc of 19, the ratio will increase if C$^{18}$O is relatively under-abundant in the disc compared to CO and H2. Through deprojection and azimuthal averaging of the image plane we detect 1.3 mm continuum emission out to 290+/-10 au,$^{13}$CO to 390+/-10 au and C$^{18}$O to 300+/-10au. We measure a radially increasing millimetre spectral index between wavelengths of 867$μ$m and 1.3 mm, which shows that grain sizes increase towards the star, with solid particles growing to cm scales in the inner disc. △ Less

Submitted 20 February, 2019; originally announced February 2019.

Circumstellar CO in metal-poor stellar winds: the highly irradiated globular cluster star 47 Tucanae V3

Authors: Iain McDonald, Martha L. Boyer, Martin A. T. Groenewegen, Eric Lagadec, Anita M. S. Richards, Gregory C. Sloan, Albert A. Zijlstra

Abstract: We report the first detection of circumstellar CO in a globular cluster. Observations with ALMA have detected the CO J=3-2 and SiO v=1 J=8-7 transitions at 345 and 344 GHz, respectively, around V3 in 47 Tucanae (NGC 104; [Fe/H] = -0.72 dex), a star on the asymptotic giant branch. The CO line is detected at 7 sigma at a rest velocity v_LSR = -40.6 km/s and expansion velocity of 3.2 +/- ~0.4 km/s. T… ▽ More We report the first detection of circumstellar CO in a globular cluster. Observations with ALMA have detected the CO J=3-2 and SiO v=1 J=8-7 transitions at 345 and 344 GHz, respectively, around V3 in 47 Tucanae (NGC 104; [Fe/H] = -0.72 dex), a star on the asymptotic giant branch. The CO line is detected at 7 sigma at a rest velocity v_LSR = -40.6 km/s and expansion velocity of 3.2 +/- ~0.4 km/s. The brighter, asymmetric SiO line may indicate a circumstellar maser. The stellar wind is slow compared to similar Galactic stars, but the dust opacity remains similar to Galactic comparisons. We suggest that the mass-loss rate is set by the levitation of material into the circumstellar environment by pulsations, but that the terminal wind-expansion velocity is determined by radiation pressure on the dust: a pulsation-enhanced dust-driven wind. We suggest the metal-poor nature of the star decreases the grain size, slowing the wind and increasing its density and opacity. Metallic alloys at high altitudes above the photosphere could also provide an opacity increase. The CO line is weaker than expected from Galactic AGB stars, but its strength confirms a model that includes CO dissociation by the strong interstellar radiation field present inside globular clusters. △ Less

Submitted 16 January, 2019; originally announced January 2019.

Comments: 5 pages, accepted MNRAS Letters

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

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

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

Submitted 31 December, 2018; originally announced January 2019.

Comments: 16 pages