Showing 1–50 of 202 results for author: Vlemmings, W

Probing the Extended Atmospheres of AGB Stars: I. Synthetic imaging of 1D hydrodynamical models at radio and (sub-)millimeter wavelengths

Authors: Behzad Bojnordi Arbab, Wouter Vlemmings, Theo Khouri, Susanne Höfner

Abstract: We investigate the observable characteristics of the extended atmospheres of AGB stars across a wide range of radio and (sub-)mm wavelengths using state-of-the-art 1D dynamical atmosphere and wind models over one pulsation period. We also study the relationships between the observable features and model properties. We further study practical distance ranges for observable sources assuming the capa… ▽ More We investigate the observable characteristics of the extended atmospheres of AGB stars across a wide range of radio and (sub-)mm wavelengths using state-of-the-art 1D dynamical atmosphere and wind models over one pulsation period. We also study the relationships between the observable features and model properties. We further study practical distance ranges for observable sources assuming the capabilities of current and upcoming observatories. We present time-variable, frequency-dependent profiles of pulsating AGB stars' atmospheres, illustrating observable features in resolved and unresolved observations, including disc brightness temperature, photosphere radius, and resolved and unresolved spectral indices. Notably, temporal variations in disc brightness temperature closely mirror the temperature variability of the stellar atmosphere. We find that while the photospheric radius decreases due to gas dilution in the layers between consecutive shocks, the increase in the observed stellar radius reflects shock propagation through the atmosphere during the expansion phase, providing a direct measurement method for the shock velocity. Furthermore, our models indicate that enhanced gas temperatures after the passage of a strong shock might be observable in the high-frequency ALMA bands as a decrease in the brightness temperature with increasing frequency. We demonstrate that synthetic observations based on state-of-the-art dynamical atmosphere and wind models are necessary for proper interpretations of current (ALMA and VLA) and future (SKA and ngVLA) observations and that multi-wavelength observations of AGB stars are crucial for empirical studies of their extended atmospheres. △ Less

Submitted 9 October, 2024; originally announced October 2024.

Comments: Accepted for publication in ApJ. 14 pages, 7 figures, and 2 tables in the main text. Appendix including figures for additional models and extended frequency ranges

One month convection timescale on the surface of a giant evolved star

Authors: W. Vlemmings, T. Khouri, B. Bojnordi Arbab, E. De Beck, M. Maercker

Abstract: The transport of energy through convection is important during many stages of stellar evolution, and is best studied in our Sun or giant evolved stars. Features that are attributed to convection are found on the surface of massive red supergiant stars. Also for lower mass evolved stars, indications of convection are found, but convective timescales and sizes remain poorly constrained. Models indic… ▽ More The transport of energy through convection is important during many stages of stellar evolution, and is best studied in our Sun or giant evolved stars. Features that are attributed to convection are found on the surface of massive red supergiant stars. Also for lower mass evolved stars, indications of convection are found, but convective timescales and sizes remain poorly constrained. Models indicate that convective motions are crucial for the production of strong winds that return the products of stellar nucleosynthesis into the interstellar medium. Here we report a series of reconstructed interferometric images of the surface of the evolved giant star R Doradus. The images reveal a stellar disc with prominent small scale features that provide the structure and motions of convection on the stellar surface. We find that the dominant structure size of the features on the stellar disc is $0.72\pm0.05$ astronomical units (au). We measure the velocity of the surface motions to vary between $-18$ and $+20$ km s$^{-1}$, which means the convective timescale is approximately one month. This indicates a possible difference between the convection properties of low-mass and high-mass evolved stars. △ Less

Submitted 13 September, 2024; v1 submitted 10 September, 2024; originally announced September 2024.

Comments: This version of the article has been accepted for publication in Nature, but is not the Version of Record and does not reflect post-acceptance improvements, or any corrections. The Version of Record is available online at http://dx.doi.org/10.1038/s41586-024-07836-9

High-Speed Outflows and Dusty Disks during the AGB to PN Transition: The PANORAMA survey

Authors: Raghvendra Sahai, Javier Alcolea, Bruce Balick, Eric G. Blackman, Valentin Bujarrabal, Arancha Castro-Carrizo, Orsola De Marco, Joel Kastner, Hyosun Kim, Eric Lagadec, Chin-Fei Lee, Laurence Sabin, M. Santander-Garcia, Carmen Sánchez Contreras, Daniel Tafoya, Toshiya Ueta, Wouter Vlemmings, Albert Zijlstra

Abstract: As mass-losing asymptotic giant branch (AGB) stars evolve to planetary nebulae (PNe), the mass outflow geometries transform from nearly spherical to extreme aspherical. The physical mechanisms governing this transformation are widely believed to be linked to binarity and the associated production of disks and fast jets during transitional (post-AGB) evolutionary stages. We are carrying out a syste… ▽ More As mass-losing asymptotic giant branch (AGB) stars evolve to planetary nebulae (PNe), the mass outflow geometries transform from nearly spherical to extreme aspherical. The physical mechanisms governing this transformation are widely believed to be linked to binarity and the associated production of disks and fast jets during transitional (post-AGB) evolutionary stages. We are carrying out a systematic ALMA survey ($P$re-planet$A$ry $N$ebulae high-angular-res$O$lution su$R$vey with $A$L$MA$ or PANORAMA) of a representative sample of bipolar and multipolar post-AGB objects. We have obtained high angular-resolution (0".1-0".4) observations of the CO(3--2) and/or 6--5 emission in order to probe the spatio-kinematic structure of the collimated outflows and the central disk/torii. The results are remarkable, generally showing the presence of bipolar or multipolar high-velocity outflows, dense toroidal waists, and in one case, a geometrically-thin circular ring around the central bipolar nebula. A high degree of point-symmetry characterizes the morphology of the mass ejecta. In this contribution, we present these and other highlights from our survey. We aim to use 2D/3D radiative transfer modeling in order to derive accurate outflow momenta, masses and mass-loss rates for our sample, and build hydrodynamical models that can explain the observed spatio-kinematic structures. These results will then be used to distinguish between different classes of PN-shaping binary interaction models. △ Less

Submitted 9 September, 2024; originally announced September 2024.

An impressionist view of V Hydrae. When MATISSE paints Asymmetric Giant Blobs

Authors: L. Planquart, C. Paladini, A. Jorissen, A. Escorza, E. Pantin, J. Drevon, B. Aringer, F. Baron, A. Chiavassa, P. Cruzalèbes, W. Danchi, E. De Beck, M. A. T. Groenewegen, S. Höfner, J. Hron, T. Khouri, B. Lopez, F. Lykou, M. Montarges, N. Nardetto, K. Ohnaka, H. Olofsson, G. Rau, A. Rosales-Guzmán, J. Sanchez-Bermudez , et al. (7 additional authors not shown)

Abstract: Our purpose is to study the effect of binary companions located within the first 10 stellar radii from the primary AGB star. In this work, we target the mass-losing carbon star V Hydrae (V Hya), looking for signatures of its companion in the dust forming region of the atmosphere. The star was observed in the L- and N-bands with the VLTI/MATISSE instrument at low spectral resolution. We reconstruct… ▽ More Our purpose is to study the effect of binary companions located within the first 10 stellar radii from the primary AGB star. In this work, we target the mass-losing carbon star V Hydrae (V Hya), looking for signatures of its companion in the dust forming region of the atmosphere. The star was observed in the L- and N-bands with the VLTI/MATISSE instrument at low spectral resolution. We reconstructed images of V Hya's photosphere and surroundings using the two bands and compared our interferometric observables with VLTI/MIDI and VISIR archival data. To constrain the dust properties, we used DUSTY to model the spectral energy distribution. The star is dominated by dust emission in the L- and N- bands. The VISIR image confirms the presence of a large-scale dusty circumstellar envelope surrounding V Hya. The MATISSE reconstructed images show asymmetric and elongated structures in both infrared bands. In the L-band, we detected an elongated shape of approximately 15 mas, likely to be of photospheric origin. In the N-band, we found a 20 mas extension North-East from the star, and perpendicular to the L-band elongated axis. The position angle and the size of the N-band extension match the prediction of the companion position at MATISSE epoch. By comparing MATISSE N-band with MIDI data, we deduce that the elongation axis in the N-band has rotated since the previous interferometric measurements 13 years ago, supporting the idea that the particle enhancement is related to the dusty clump moving along with the companion. The MATISSE images unveil the presence of a dust enhancement at the companion position, opening new doors for further analysis on the binary interaction with an AGB component. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 687, A306 (2024)

Probing the dynamical and kinematical structures of detached shells around AGB stars

Authors: M. Maercker, E. De Beck, T. Khouri, W. H. T. Vlemmings, J. Gustafsson, H. Olofsson, D. Tafoya, F. Kerschbaum, M. Lindqvist

Abstract: Aims. We aim to resolve the spatial and kinematic sub-structures in five detached-shell sources to provide detailed constraints for hydrodynamic models that describe the formation and evolution of the shells. Methods. We use observations of the 12 CO (1-0) emission towards five carbon-AGB stars with ALMA. The data have angular resolutions of 0.3 arcsec to 1arcsec and a velocity resolution of 0.3 k… ▽ More Aims. We aim to resolve the spatial and kinematic sub-structures in five detached-shell sources to provide detailed constraints for hydrodynamic models that describe the formation and evolution of the shells. Methods. We use observations of the 12 CO (1-0) emission towards five carbon-AGB stars with ALMA. The data have angular resolutions of 0.3 arcsec to 1arcsec and a velocity resolution of 0.3 km/s . This enables us to quantify spatial and kinematic structures in the shells. Results. The observed emission is separated into two distinct components: a more coherent, bright outer shell and a more filamentary, fainter inner shell. The kinematic information shows that the inner sub-shells move at a higher velocity relative to the outer sub-shells. The observed sub-structures confirm the predictions from hydrodynamical models. However, the models do not predict a double-shell structure, and the CO emission likely only traces the inner and outer edges of the shell, implying a lack of CO in the middle layers of the detached shell. Previous estimates of the masses and temperatures are consistent with originating mainly from the brighter subshell, but the total shell masses are likely lower limits. Conclusions. The observed spatial and kinematical splittings of the shells appear consistent with results from hydrodynamical models, provided the CO emission does not trace the H2 density distribution in the shell but rather traces the edges of the shells. It is therefore not possible to constrain the total shell mass based on the CO observations alone. Complementary observations of, e.g., CI as a dissociation product of CO would be necessary to understand the distribution of CO compared to H2. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: 18 pages (incl. 5 pages Appendix), 13 Figures

Journal ref: A&A 687, A112 (2024)

Molecular line polarisation from the circumstellar envelopes of Asymptotic Giant Branch stars

Authors: W. H. T. Vlemmings, B. Lankhaar, L. Velilla-Prieto

Abstract: Polarisation observations of masers in the circumstellar envelopes (CSEs) around Asymptotic Giant Branch (AGB) stars have revealed strong magnetic fields. However, masers probe only specific lines-of-sight through the CSE. Non-masing molecular line polarisation observation can more directly reveal the large scale magnetic field morphology and hence probe the effect of the magnetic field on AGB mas… ▽ More Polarisation observations of masers in the circumstellar envelopes (CSEs) around Asymptotic Giant Branch (AGB) stars have revealed strong magnetic fields. However, masers probe only specific lines-of-sight through the CSE. Non-masing molecular line polarisation observation can more directly reveal the large scale magnetic field morphology and hence probe the effect of the magnetic field on AGB mass-loss and the shaping of the AGB wind. Observations and models of CSE molecular line polarisation can now be used to describe the magnetic field morphology and estimate its strength throughout the entire CSE. We use observations taken with ALMA of molecular line polarisation in the envelope of two AGB stars (CW~Leo and R~Leo). We model the observations using the multi-dimensional polarised radiative transfer tool PORTAL. We find linearly polarised emission, with maximum fractional polarisation on the order of a few percent, in several molecular lines towards both stars. We can explain the observed differences in polarisation structure between the different molecular lines by alignment of the molecules by a combination of the Goldreich-Kylafis effect and radiative alignment effects. We specifically show that the polarisation of CO traces the morphology of the magnetic field. Competition between the alignment mechanisms allows us to describe the behaviour of the magnetic field strength with radius throughout the circumstellar envelope of CW~Leo. The magnetic field strength derived using this method is inconsistent with the magnetic field strength derived using a structure function analysis of the CO polarisation and the strength previously derived using CN Zeeman observations. In contrast with CW~Leo, the magnetic field in the outer envelope of R~Leo appears to be advected outwards by the stellar wind. {abridged abstract} △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: 23 pages, 21 figures, accepted for publication in A&A

Journal ref: A&A 686, A274 (2024)

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 Vlemmings

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 6 August, 2024; v1 submitted 1 March, 2024; originally announced March 2024.

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

An empirical view of the extended atmosphere and inner envelope of the AGB star R Doradus I. Physical model based on CO lines

Authors: T. Khouri, H. Olofsson, W. H. T. Vlemmings, T. Schirmer, D. Tafoya, M. Maercker, E. De Beck, L. -Å. Nyman, M. Saberi

Abstract: The mass loss experienced on the asymptotic giant branch (AGB) at the end of the lives of low- and intermediate-mass stars is widely accepted to rely on radiation pressure acting on dust grains formed in the extended AGB atmospheres. The interaction of convection, stellar pulsation, and heating and cooling processes cause the density, velocity and temperature distributions in the inner regions of… ▽ More The mass loss experienced on the asymptotic giant branch (AGB) at the end of the lives of low- and intermediate-mass stars is widely accepted to rely on radiation pressure acting on dust grains formed in the extended AGB atmospheres. The interaction of convection, stellar pulsation, and heating and cooling processes cause the density, velocity and temperature distributions in the inner regions of the envelope to be complex, making the dust-formation process difficult to calculate. Hence, characterising the extended atmospheres and inner outflow empirically is paramount to advance our understanding of the dust-formation and wind-driving processes. To this end, we observe the AGB star R Dor using ALMA and modelled the $^{12}$CO $v=0, J=2-1$, $v=1, J=2-1$ and $3-2$ and $^{13}$CO $v=0, J=3-2$ lines using the 3D radiative transfer code LIME up to a distance of $\sim 4$ times the radius of the star at sub-mm wavelengths. We find a complex velocity field with structure down to scales at least equal to the resolution of the observations. The observed maps are well reproduced assuming spherical symmetry for the gas temperature and density distributions. We find the radial profiles of these two quantities to be very steep close to the star and shallower for radii larger than $\sim 1.6$ times the stellar sub-mm radius. This change is consistent with the transition between extended atmosphere and outflow. We constrain the standard deviation of the stochastic velocity distribution in the large-scale outflow to be $\lesssim 0.4$ km/s. We observe two emission blobs in the CO $v=0, J=2-1$ line and find their gas densities and radial velocities to be substantially larger than those of the surrounding gas. Monitoring the evolution of these blobs will lead to a better understanding of the role of these structures in the mass-loss process of R Dor. △ Less

Submitted 21 February, 2024; originally announced February 2024.

MSC Class: 85-02

Maser polarization through anisotropic pumping

Authors: Boy Lankhaar, Gabriele Surcis, Wouter Vlemmings, Violette Impellizzeri

Abstract: (Abridged) Polarized emission from masers is an excellent tool to study magnetic fields in maser sources. The linear polarization of most masers is understood as an interplay of maser saturation and anisotropic pumping. However, for the latter mechanism, no quantitative modeling has been presented yet. We present a comprehensive model of maser polarization, including quantitative modeling of both… ▽ More (Abridged) Polarized emission from masers is an excellent tool to study magnetic fields in maser sources. The linear polarization of most masers is understood as an interplay of maser saturation and anisotropic pumping. However, for the latter mechanism, no quantitative modeling has been presented yet. We present a comprehensive model of maser polarization, including quantitative modeling of both anisotropic pumping and the effects of maser saturation on the polarization of masers. We extend regular maser excitation modeling with a dimension that describes the molecular population alignments, as well as including the linear polarization dimension to the radiative transfer. The results of the excitation analysis yield the anisotropic pumping and decay parameters, that are subsequently used in one-dimensional proper maser polarization radiative transfer modeling. We present the anisotropic pumping parameters for a variety of transitions from class I CH$_3$OH masers, H$_2$O masers and SiO masers. SiO masers are highly anisotropically pumped due to them occurring in the vicinity of a late-type star, that irradiates the maser region with a strong directional radiation field. Class I CH$_3$OH masers and H$_2$O masers occur in association with shocks, and they are modestly anisotropically pumped due to the anisotropy of the excitation region. Our modeling constitutes the first quantitative constraints on the anisotropic pumping of masers. We find that anisotropic pumping can explain the high polarization yields of SiO masers, as well as the modest polarization of unsaturated class I CH$_3$OH masers. We predict that the $183$ GHz H$_2$O maser is strongly anisotropically pumped. Finally, we outline a mechanism through which non-Zeeman circular polarization is produced, when the magnetic field changes direction along the propagation through an anisotropically pumped maser. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 22 pages, 10 figures, accepted to A&A

Charting Circumstellar Chemistry of Carbon-rich AGB Stars: I. ALMA 3 mm spectral surveys

Authors: R. Unnikrishnan, E. De Beck, L. A. Nyman, H. Olofsson, W. H. T. Vlemmings, D. Tafoya, M. Maercker, S. B. Charnley, M. A. Cordiner, I. de Gregorio, E. Humphreys, T. J. Millar, M. G. Rawlings

Abstract: AGB stars are major contributors to the chemical enrichment of the ISM through nucleosynthesis and extensive mass loss. Most of our current knowledge of AGB atmospheric and circumstellar chemistry, in particular in a C-rich environment, is based on observations of the carbon star IRC+10216. We aim to obtain a more generalised understanding of the chemistry in C-rich AGB CSEs by studying a sample o… ▽ More AGB stars are major contributors to the chemical enrichment of the ISM through nucleosynthesis and extensive mass loss. Most of our current knowledge of AGB atmospheric and circumstellar chemistry, in particular in a C-rich environment, is based on observations of the carbon star IRC+10216. We aim to obtain a more generalised understanding of the chemistry in C-rich AGB CSEs by studying a sample of three carbon stars, IRAS15194-5115, IRAS15082-4808, and IRAS07454-7112, and test the archetypal status often attributed to IRC+10216. We performed spatially resolved, unbiased spectral surveys in ALMA Band 3. We identify a total of 132 rotational transitions from 49 molecular species. There are two main morphologies of the brightness distributions: centrally-peaked (e.g. HCN) and shell-like (e.g. C$_2$H). We estimated the sizes of the molecular emitting regions using azimuthally-averaged radial profiles of the line brightness distributions, and derived abundance estimates. Of the shell distributions, the cyanopolyynes peak at slightly smaller radii than the hydrocarbons, and CN and HNC show the most extended emission. The emitting regions for each species are the smallest for IRAS07454-7112. We find that, within the uncertainties of the analysis, the three stars present similar abundances for most species, also compared to IRC+10216. We find that SiO is more abundant in our three stars compared to IRC+10216. Our estimated isotopic ratios match well the literature values for the sources. The observed circumstellar chemistry appears very similar across our sample and compared to that of IRC+10216, both in terms of the relative location of the emitting regions and molecular abundances. This implies that, to a first approximation, the chemical models tailored to IRC+10216 are able to reproduce the observed chemistry in C-rich envelopes across roughly an order of magnitude in wind density. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 21 pages. 13 figures and 8 tables in the main text. 5 appendices contain additional figures and tables. Appendix tables are available in electronic form at the CDS, along with the reprocessed ALMA cubes and spectra, at http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ or via anonymous ftp to cdsarc.u-strasbg.fr(130.79.128.5). Accepted for publication in Astronomy & Astrophysics

An SiO Toroid and Wide-angle Outflow associated with the Massive Protostar W75N(B)-VLA2

Authors: Jose F. Gomez, Jose M. Torrelles, Josep M. Girart, Gabriele Surcis, Jeong-Sook Kim, Jorge Canto, Guillem Anglada, Salvador Curiel, Wouter H. T. Vlemmings, Carlos Carrasco-Gonzalez, Adriana R. Rodriguez-Kamenetzky, Soon-Wook Kim, Ciriaco Goddi, Huib J. van Langevelde, Alvaro Sanchez-Monge

Abstract: We have carried out ALMA observations of the massive star-forming region W75N(B), which contains the massive protostars VLA1, VLA2, and VLA3. Particularly, VLA2 is an enigmatic protostar associated with a wind-driven H$_2$O maser shell, which has evolved from an almost isotropic outflow to a collimated one in just 20 years. The shell expansion seemed to be halted by an obstacle located to the nort… ▽ More We have carried out ALMA observations of the massive star-forming region W75N(B), which contains the massive protostars VLA1, VLA2, and VLA3. Particularly, VLA2 is an enigmatic protostar associated with a wind-driven H$_2$O maser shell, which has evolved from an almost isotropic outflow to a collimated one in just 20 years. The shell expansion seemed to be halted by an obstacle located to the northeast of VLA2. Here we present our findings from observing the 1.3 mm continuum and H$_2$CO and SiO emission lines. Within a region of $\sim 30"$ ($\sim 39,000$ au) diameter, we have detected 40 compact mm-continuum sources, three of them coinciding with VLA1, VLA2, and VLA3. While the H$_2$CO emission is mainly distributed in a fragmented structure around the three massive protostars, but without any of the main H$_2$CO clumps spatially coinciding with them, the SiO is highly concentrated on VLA2, indicating the presence of very strong shocks generated near this protostar. The SiO emission is clearly resolved into an elongated structure ($\sim 0.6"\times0.3"$; $\sim 780$ au$\times$390 au) perpendicular to the major axis of the wind-driven maser shell. The structure and kinematics of the SiO emission are consistent with a toroid and a wide-angle outflow surrounding a central mass of $\sim 10$ M$_{\odot}$, thus supporting previous theoretical predictions regarding the evolution of the outflow. Additionally, we have identified the expected location and estimated the gas density of the obstacle that is hindering the expansion of the maser shell. △ Less

Submitted 13 October, 2023; v1 submitted 12 October, 2023; originally announced October 2023.

Comments: To be published in The Astrophysical Journal Letters. Sixteen pages, seven figures. Updated metadata

Exploring the environment, magnetic fields, and feedback effects of massive high-redshift galaxies with [CII]

Authors: K. Kade, K. K. Knudsen, W. Vlemmings, F. Stanley, B. Gullberg, S. Konig

Abstract: Massive galaxies are expected to grow through different transformative evolutionary phases where high-redshift starburst galaxies and quasars are examples of such phases. The physical mechanisms driving these phases include companion galaxy interactions, active galactic nuclei feedback, and magnetic fields. Our aim is to characterize the physical properties and the environment of the submillimeter… ▽ More Massive galaxies are expected to grow through different transformative evolutionary phases where high-redshift starburst galaxies and quasars are examples of such phases. The physical mechanisms driving these phases include companion galaxy interactions, active galactic nuclei feedback, and magnetic fields. Our aim is to characterize the physical properties and the environment of the submillimeter galaxy AzTEC-3 at z = 5.3 and the lensed quasar BRI 0952-0115 at z = 4.4, to set a limit on the polarization properties, as well as placing both in the broader context of galaxy evolution. We used full polarization, sub-arcsecond-resolution, ALMA band-7 observations of both BRI 0952-0115 and AzTEC-3 and detect [CII] line emission towards both galaxies, along with companions in each field. We present an updated gravitational lensing model for BRI 0952-0115. We present infrared luminosities, star-formation rates, and [CII] line to infrared luminosity ratios for each source. The [CII] emission line profile for both BRI 0952-0115 and AzTEC-3 exhibit a broad, complex morphology, indicating the possible presence of outflows. We present evidence of a 'gas bridge' between AzTEC-3 and a companion source. Using a simple dynamical mass estimate for the sources, we suggest that both systems are undergoing minor or major mergers. No polarization is detected for the [CII], placing an upper limit below that of theoretical predictions. Our results show that high-velocity wings are detected, indicating possible signs of massive outflows; however, the presence of companion galaxies can affect the final interpretation. Furthermore, the results provide additional evidence in support of the hypothesis that massive galaxies form in overdense regions, growing through interactions. Finally, strong, ordered magnetic fields are unlikely to exist at the kiloparsec scale in the two studied sources. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Monitoring of the polarized $H_{2}O$ maser emission around the massive protostars W75N(B)-VLA1 and W75N(B)-VLA2

Authors: G. Surcis, W. H. T. Vlemmings, C. Goddi, J. M. Torrelles, J. F. Gómez, A. Rodríguez-Kamenetzky, C. Carrasco-González, S. Curiel, S. -W. Kim, J. -S. Kim, H. J. van Langevelde

Abstract: Several radio sources have been detected in the HMSFR W75N(B), among them the massive YSOs VLA1 and VLA2 are of great interest. These are thought to be in different evolutionary stages. In particular, VLA1 is at the early stage of the photoionization and it is driving a thermal radio jet, while VLA2 is a thermal, collimated ionized wind surrounded by a dusty disk or envelope. In both sources 22 GH… ▽ More Several radio sources have been detected in the HMSFR W75N(B), among them the massive YSOs VLA1 and VLA2 are of great interest. These are thought to be in different evolutionary stages. In particular, VLA1 is at the early stage of the photoionization and it is driving a thermal radio jet, while VLA2 is a thermal, collimated ionized wind surrounded by a dusty disk or envelope. In both sources 22 GHz water masers have been detected in the past. Those around VLA1 show a persistent distribution along the radio jet and those around VLA2 have instead traced the evolution from a non-collimated to a collimated outflow over a period of 20 years. By monitoring the polarized emission of the water masers around both VLA1 and VLA2 over a period of 6 years, we aim to determine whether the maser distributions show any variation over time and whether the magnetic field behaves accordingly. The EVN was used in full polarization and phase-reference mode to measure the absolute positions of the masers and to determine both the orientation and the strength of the magnetic field. We observed four epochs separated by two years from 2014 to 2020. We detected polarized emission from the water masers around both the YSOs in all the epochs. We find that the masers around VLA1 are tracing a nondissociative shock originating from the expansion of the thermal radio jet, while the masers around VLA2 are tracing an asymmetric expansion of the gas that is halted in the northeast where the gas likely encounters a very dense medium. We also found that the magnetic field inferred from the water masers in each epoch can be considered as a portion of a quasi-static magnetic field estimated in that location rather than in that time. This allowed us to study locally the morphology of the magnetic field around both YSOs in a larger area by considering the vectors estimated in all the epochs as a whole. △ Less

Submitted 28 February, 2023; originally announced March 2023.

Comments: 32 pages, 15 figures, 17 tables, accepted by Astronomy & Astrophysics

Journal ref: A&A 673, A10 (2023)

Polarisation of molecular lines in the circumstellar envelope of the post-Asymptotic Giant Branch star OH 17.7-2.0

Authors: W. H. T. Vlemmings, D. Tafoya

Abstract: (abridged) The role of magnetic field in the shaping of Planetary Nebulae (PNe), either directly or indirectly after being enhanced by binary interaction, has long been a topic of debate. Large scale magnetic fields around pre-PNe have been inferred from polarisation observations of masers. However, because masers probe very specific regions, it is still unclear if the maser results are representa… ▽ More (abridged) The role of magnetic field in the shaping of Planetary Nebulae (PNe), either directly or indirectly after being enhanced by binary interaction, has long been a topic of debate. Large scale magnetic fields around pre-PNe have been inferred from polarisation observations of masers. However, because masers probe very specific regions, it is still unclear if the maser results are representative of the intrinsic magnetic field in the circumstellar envelope (CSE). Molecular line polarisation can provide important information about the magnetic field. A comparison between the field morphology determined from maser observations and that observed in the more diffuse CO gas, can reveal if the two tracers probe the same magnetic field. We compare observations taken with ALMA of molecular line polarisation around the post-Asymptotic Giant Branch)/pre-PNe star OH~17.7-2.0 with previous observations of polarisation in the 1612~MHz OH maser region. We detect CO~$J=2-1$ molecular line polarisation at a level of $\sim4\%$ that displays an ordered linear polarisation structure. We find that, correcting for Faraday rotation of the OH~maser linear polarisation vectors, the OH and CO linearly polarised emission trace the same large scale magnetic field. A structure function analysis of the CO linear polarisation reveals a plane-of-the-sky magnetic field strength of $B_\perp\sim1$~mG in the CO region, consistent with previous OH Zeeman observations. The consistency of the ALMA CO molecular line polarisation with maser observations indicate that both can be used to determine the magnetic field in CSEs. The existence of a strong, ordered, magnetic-field around OH 17.7-2.0 indicates that magnetic fields are likely involved in the formation of this bipolar pre-PNe. △ Less

Submitted 2 February, 2023; originally announced February 2023.

Comments: 13 pages, 9 figures; accepted for publication in A&A

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

The messy death of a multiple star system and the resulting planetary nebula as observed by JWST

Authors: Orsola De Marco, Muhammad Akashi, Stavros Akras, Javier Alcolea, Isabel Aleman, Philippe Amram, Bruce Balick, Elvire De Beck, Eric G. Blackman, Henri M. J. Boffin, Panos Boumis, Jesse Bublitz, Beatrice Bucciarelli, Valentin Bujarrabal, Jan Cami, Nicholas Chornay, You-Hua Chu, Romano L. M. Corradi, Adam Frank, Guillermo Garcia-Segura, D. A. Garcia-Hernandez, Jorge Garcia-Rojas, Veronica Gomez-Llanos, Denise R. Goncalves, Martin A. Guerrero , et al. (44 additional authors not shown)

Abstract: Planetary nebulae (PNe), the ejected envelopes of red giant stars, provide us with a history of the last, mass-losing phases of 90 percent of stars initially more massive than the Sun. Here, we analyse James Webb Space Telescope (JWST) Early Release Observation (ERO) images of the PN NGC3132. A structured, extended H2 halo surrounding an ionised central bubble is imprinted with spiral structures,… ▽ More Planetary nebulae (PNe), the ejected envelopes of red giant stars, provide us with a history of the last, mass-losing phases of 90 percent of stars initially more massive than the Sun. Here, we analyse James Webb Space Telescope (JWST) Early Release Observation (ERO) images of the PN NGC3132. A structured, extended H2 halo surrounding an ionised central bubble is imprinted with spiral structures, likely shaped by a low-mass companion orbiting the central star at 40-60 AU. The images also reveal a mid-IR excess at the central star interpreted as a dusty disk, indicative of an interaction with another, closer companion. Including the previously known, A-type visual companion, the progenitor of the NGC3132 PN must have been at least a stellar quartet. The JWST images allow us to generate a model of the illumination, ionisation and hydrodynamics of the molecular halo, demonstrating the power of JWST to investigate complex stellar outflows. Further, new measurements of the A-type visual companion allow us to derive the value for the mass of the progenitor of a central star to date with excellent precision: 2.86+/-0.06 Mo. These results serve as path finders for future JWST observations of PNe providing unique insight into fundamental astrophysical processes including colliding winds, and binary star interactions, with implications for supernovae and gravitational wave systems. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: 32 pages, 5 figures for the main article. 12 pages 8 figures for the supplementary material

Journal ref: Nature Astronomy, 2022, Vol. 6, p. 1421

Distance estimates for AGB stars from parallax measurements

Authors: M. Andriantsaralaza, S. Ramstedt, W. H. T Vlemmings, E. De Beck

Abstract: Estimating the distances to asymptotic giant branch (AGB) stars using optical measurements of their parallaxes is not straightforward because of the large uncertainties introduced by their dusty envelopes, their large angular sizes, and their surface brightness variability. We compared the parallaxes from Gaia DR3 with parallaxes measured with maser observations with very long baseline interferome… ▽ More Estimating the distances to asymptotic giant branch (AGB) stars using optical measurements of their parallaxes is not straightforward because of the large uncertainties introduced by their dusty envelopes, their large angular sizes, and their surface brightness variability. We compared the parallaxes from Gaia DR3 with parallaxes measured with maser observations with very long baseline interferometry (VLBI) to determine a statistical correction factor for the DR3 parallaxes using a sub-sample of 33 maser-emitting oxygen-rich nearby AGB stars. We found that the nominal errors on the Gaia DR3 parallaxes are underestimated by a factor of 5.44 for the brightest sources ($G<8$ mag). Fainter sources ($8\leq G<12$) require a lower parallax error inflation factor of 2.74. We obtain a Gaia DR3 parallax zero-point offset of -0.077 mas for bright AGB stars. The offset becomes more negative for fainter AGB stars. We calculated the distances of about 200 AGB stars in the DEATHSTAR project using a Bayesian statistical approach on the corrected DR3 parallaxes and a prior based on the previously determined Galactic distribution of AGB stars. We found that the derived Gaia distances are associated with significant, asymmetrical errors for more than 40% of the sources in our sample. We performed radiative transfer modelling of the stellar and dust emission to determine the luminosity of the sources in the VLBI sub-sample based on the distances derived from maser parallaxes, and derived a new bolometric period-luminosity relation for Galactic oxygen-rich Mira variables. A new distance catalogue based on these results is provided for the sources in the DEATHSTAR sample. We found that a RUWE (re-normalised unit weight error) below 1.4 does not guarantee reliable distance estimates and we advise against the use of only the RUWE to measure the quality of Gaia DR3 astrometric data for individual AGB stars. △ Less

Submitted 7 November, 2022; v1 submitted 8 September, 2022; originally announced September 2022.

Comments: 15 pages, 12 figures, 3 appendices, accepted for publication in A&A, corrected typos, distance catalogue available on the CDS (and any updates that may occur)

Journal ref: A&A 667, A74 (2022)

First detection of AlF line emission towards M-type AGB stars

Authors: M. Saberi, T. Khouri, L. Velilla-Prieto, J. P. Fonfría, W. H. T. Vlemmings, S. Wedemeyer

Abstract: The nucleosynthesis production of fluorine (F) is still a matter of debate. Asymptotic giant branch (AGB) stars are one of the main candidates for F production. However, their contribution to the total F budget is not fully known due to the lack of observations. In this paper, we report the detection of AlF line emission, one of the two main carriers of F in the gas-phase in the outflow of evolved… ▽ More The nucleosynthesis production of fluorine (F) is still a matter of debate. Asymptotic giant branch (AGB) stars are one of the main candidates for F production. However, their contribution to the total F budget is not fully known due to the lack of observations. In this paper, we report the detection of AlF line emission, one of the two main carriers of F in the gas-phase in the outflow of evolved stars, towards five nearby oxygen-rich AGB stars, $o$ Ceti, R Leo, IK Tau, R Dor, and W Hya. From spatially resolved observations, we estimated the AlF emitting region with a radius $\sim11R_{\star}$ for $o$ Ceti and $\sim9R_{\star}$ for R Leo. From population diagram analysis, we report the AlF column densities of $\sim 5.8\times10^{15}$ cm$^{-2}$ and $\sim 3\times10^{15}$ cm$^{-2}$ for $o$ Ceti and R Leo, respectively, within these regions. For $o$ Ceti, we used the C$^{18}$O ($v=0$, $J=3-2$) observations to estimate the H$_2$ column density of the emitting region. We found a fractional abundance of $f_{\rm AlF/H_2}\sim(2.5\pm1.7)\times10^{-8}$. This gives a lower limit on the F budget in $o$ Ceti and is compatible with the solar F budget $f_{\rm F/H_2}=(5\pm2)\times10^{-8}$. For R Leo, a fractional abundance $f_{\rm AlF/H_2}=(1.2\pm0.5)\times10^{-8}$ is estimated. For other sources, we cannot precisely determine the emitting region based on the available data. Assuming an emitting region with a radius of $\sim 11R_{\star}$ and the rotational temperatures derived for $o$ Ceti and R Leo, we crudely approximated the AlF column density to be $\sim(1.2-1.5)\times10^{15}$ cm$^{-2}$ in W Hya, $\sim(2.5-3.0)\times10^{14}$ cm$^{-2}$ in R Dor, and $\sim(0.6-1.0)\times10^{16}$ cm$^{-2}$ in IK Tau. These result in fractional abundances within a range of $f_{\rm AlF/H_2}\sim(0.1-4)\times10^{-8}$ in W Hya, R Dor, and IK Tau. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: 11 pages, 9 figures

Journal ref: A&A 663, A54 (2022)

First images of the molecular gas around a born-again star revealed by ALMA

Authors: Daniel Tafoya, Jesús A. Toalá, Ramlal Unnikrishnan, Wouter H. T. Vlemmings, Martín A. Guerrero, Stefan Kimeswenger, Peter A. M. van Hoof, Luis A. Zapata, Sandra P. Treviño-Morales, Janis B. Rodríguez-González

Abstract: Born-again stars allow probing stellar evolution in human timescales and provide the most promising path for the formation of hydrogen-deficient post-asymptotic giant branch objects, but their cold and molecular components remain poorly explored. Here we present ALMA observations of V605 Aql that unveil for the first time the spatio-kinematic distribution of the molecular material associated to a… ▽ More Born-again stars allow probing stellar evolution in human timescales and provide the most promising path for the formation of hydrogen-deficient post-asymptotic giant branch objects, but their cold and molecular components remain poorly explored. Here we present ALMA observations of V605 Aql that unveil for the first time the spatio-kinematic distribution of the molecular material associated to a born-again star. Both the continuum and molecular line emission exhibit a clumpy ring-like structure with a total extent of $\approx$1$^{\prime\prime}$ in diameter. The bulk of the molecular emission is interpreted as being produced in a radially-expanding disk-like structure with an expansion velocity v$_{\rm exp}$$\sim$90 km s$^{-1}$ and an inclination $i$$\approx$60$^{\circ}$ with respect to the line-of-sight. The observations also reveal a compact high-velocity component, v$_{\rm exp}$$\sim$280 km s$^{-1}$, that is aligned perpendicularly to the expanding disk. This component is interpreted as a bipolar outflow with a kinematical age $τ$$\lesssim$20 yr, which could either be material that is currently being ejected from V605 Aql, or it is being dragged from the inner parts of the disk by a stellar wind. The dust mass of the disk is in the range $M_{\rm dust}$$\sim$0.2-8$\times$10$^{-3}$ M$_{\odot}$, depending on the dust absorption coefficient. The mass of the CO is $M_{\rm CO}$$\approx$1.1$\times10^{-5}$ $M_{\odot}$, which is more than three orders of magnitude larger than the mass of the other detected molecules. We estimate a $^{12}$C/$^{13}$C ratio of 5.6$\pm$0.6, which is consistent with the single stellar evolution scenario in which the star experienced a very late thermal pulse instead of a nova-like event as previously suggested. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 9 pages, 5 Figures; accepted to ApJL

Observational identification of a sample of likely recent Common-Envelope Events

Authors: Theo Khouri, Wouter Vlemmings, Daniel Tafoya, Andrés F. Pérez-Sánchez, Carmen Sánchez Contreras, José F. Gómez, Hiroshi Imai, Raghvendra Sahai

Abstract: One of the most poorly understood stellar evolutionary paths is that of binary systems undergoing common-envelope evolution, when the envelope of a giant star engulfs the orbit of a companion. Although this interaction leads to a great variety of astrophysical systems, direct empirical studies are difficult because few objects experiencing common-envelope evolution are known. We present ALMA obser… ▽ More One of the most poorly understood stellar evolutionary paths is that of binary systems undergoing common-envelope evolution, when the envelope of a giant star engulfs the orbit of a companion. Although this interaction leads to a great variety of astrophysical systems, direct empirical studies are difficult because few objects experiencing common-envelope evolution are known. We present ALMA observations towards sources known as water fountains that reveal they had low initial masses ($<4~{\rm M}_\odot$) and ejected a significant fraction of it over less than a few hundred years. The only mechanism able to explain such rapid mass ejection is common-envelope evolution. Our calculations show that the water-fountain sample accounts for a large fraction of the systems in our Galaxy which have just experienced the common-envelope phase. Since water-fountain sources show characteristic fast bipolar outflows, outflows and jets likely play an important role right before, during or immediately after the common-envelope phase. △ Less

Submitted 17 December, 2021; originally announced December 2021.

Comments: Published on Nature Astronomy on Dec 16, 2021: https://rdcu.be/cDlX8

EVN observations of 6.7 GHz methanol maser polarization in massive star-forming regions. V. Completion of the flux-limited sample

Authors: G. Surcis, W. H. T. Vlemmings, H. J. van Langevelde, B. Hutawarakorn Kramer, A. Bartkiewicz

Abstract: Although the role of magnetic fields in launching molecular outflows in massive YSOs has been convincingly demonstrated by theoretical arguments, observationally, the alignment of the magnetic field lines with the molecular outflows is still under debate. We aim to complete the measurements of the direction of the magnetic fields at mas resolution around a sample of massive star-forming regions (M… ▽ More Although the role of magnetic fields in launching molecular outflows in massive YSOs has been convincingly demonstrated by theoretical arguments, observationally, the alignment of the magnetic field lines with the molecular outflows is still under debate. We aim to complete the measurements of the direction of the magnetic fields at mas resolution around a sample of massive star-forming regions (MSFRs) to determine whether the magnetic field and outflows are aligned. In 2012, we started a large VLBI campaign with the EVN to measure the magnetic field orientation and strength toward a sample of 31 MSFRs (the flux-limited sample) by analyzing the polarized emission of 6.7GHz CH3OH masers. In the previous papers of the series, we have presented 80% of the sample. Here, we report the linearly and circularly polarized emission of 6.7GHz CH3OH masers toward the last five MSFRs of the flux-limited sample. The sources are G30.70-0.07, G30.76-0.05, G31.28+0.06, G32.03+0.06, and G69.52-0.97. We detected a total of 209 masers, 15% of which show linearly polarized emission (0.07%-16.7%), and 2% of which show circularly polarized emission (0.2%-4.2%). Zeeman splitting was measured toward G30.70-0.07, G32.03+0.06, and G69.52-0.97. The statistical analysis of the entire flux-limited sample shows that the observations are consistent with a bimodal distribution in the difference between the 3D magnetic field direction and the outflow axis, with half the magnetic field directions being perpendicular and the other half being parallel to the outflow. In addition, we determined that typical values of the linear and circular polarization fractions for 6.7 GHz CH3OH masers are Pl=1.0%-2.5% and Pv=0.5%-0.75%, respectively. We found that a typical Zeeman splitting is in the range between 0.5 m/s and 2.0 m/s. This would correspond to 9 mG<$|B_{||}|$<40 mG if F=3->4 is the most favored hyperfine transition. △ Less

Submitted 15 November, 2021; originally announced November 2021.

Comments: 22 pages, 12 figures, 9 tables, accepted by Astronomy & Astrophysics. arXiv admin note: text overlap with arXiv:1902.08210

Journal ref: A&A 658, A78 (2022)

Milliarcsecond localisation of the repeating FRB 20201124A

Authors: K. Nimmo, D. M. Hewitt, J. W. T. Hessels, F. Kirsten, B. Marcote, U. Bach, R. Blaauw, M. Burgay, A. Corongiu, R. Feiler, M. P. Gawroński, M. Giroletti, R. Karuppusamy, A. Keimpema, M. A. Kharinov, M. Lindqvist, G. Maccaferri, A. Melnikov, A. Mikhailov, O. S. Ould-Boukattine, Z. Paragi, M. Pilia, A. Possenti, M. P. Snelders, G. Surcis , et al. (6 additional authors not shown)

Abstract: Very long baseline interferometric (VLBI) localisations of repeating fast radio bursts (FRBs) have demonstrated a diversity of local environments: from nearby star-forming regions to globular clusters. Here we report the VLBI localisation of FRB 20201124A using an ad-hoc array of dishes that also participate in the European VLBI Network (EVN). In our campaign, we detected 18 total bursts from FRB… ▽ More Very long baseline interferometric (VLBI) localisations of repeating fast radio bursts (FRBs) have demonstrated a diversity of local environments: from nearby star-forming regions to globular clusters. Here we report the VLBI localisation of FRB 20201124A using an ad-hoc array of dishes that also participate in the European VLBI Network (EVN). In our campaign, we detected 18 total bursts from FRB 20201124A at two separate epochs. By combining the visibilities from both observing epochs, we were able to localise FRB 20201124A with a 1-$σ$ error of 4.5 milliarcseconds (mas). We use the relatively large burst sample to investigate astrometric accuracy, and find that for $\gtrsim20$ baselines ($\gtrsim7$ dishes) that we can robustly reach milliarcsecond precision even using single-burst data sets. Sub-arcsecond precision is still possible for single bursts, even when only $\sim$ six baselines (four dishes) are available. We explore two methods for determining the individual burst positions: the peaks of the dirty maps and a Gaussian fit to the cross fringe pattern on the dirty maps. We found the latter to be more reliable due to the lower mean and standard deviation in the offsets from the FRB position. Our VLBI work places FRB 20201124A 705$\pm$26 mas (1-$σ$ errors) from the optical centre of the host galaxy, and consistent with originating from within the recently-discovered extended radio structure associated with star-formation in the host galaxy. Future high-resolution optical observations, e.g. with Hubble Space Telescope, can determine the proximity of our FRB 20201124A VLBI position to nearby knots of star formation. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: submitted, comments welcome

Discovery of Molecular Line Polarization in the Disk of TW Hya

Authors: Richard Teague, Chat L. H. Hull, Stéphane Guilloteau, Edwin A. Bergin, Anne Dutrey, Thomas Henning, Rolf Kuiper, Dmitry Semenov, Ian W. Stephens, Wouter H. T. Vlemmings

Abstract: We report observations of polarized line and continuum emission from the disk of TW~Hya using the Atacama Large Millimeter/submillimeter Array. We target three emission lines, $^{12}$CO (3-2), $^{13}$CO (3-2) and CS (7-6), to search for linear polarization due to the Goldreich-Kylafis effect, while simultaneously tracing the continuum polarization morphology at 332\,GHz (900\,\micron{}), achieving… ▽ More We report observations of polarized line and continuum emission from the disk of TW~Hya using the Atacama Large Millimeter/submillimeter Array. We target three emission lines, $^{12}$CO (3-2), $^{13}$CO (3-2) and CS (7-6), to search for linear polarization due to the Goldreich-Kylafis effect, while simultaneously tracing the continuum polarization morphology at 332\,GHz (900\,\micron{}), achieving a spatial resolution of 0.5\arcsec{} (30~au). We detect linear polarization in the dust continuum emission; the polarization position angles show an azimuthal morphology, and the median polarization fraction is $\sim$\,0.2\%, comparable to previous, lower frequency observations. Adopting a `shift-and-stack' technique to boost the sensitivity of the data, combined with a linear combination of the $Q$ and $U$ components to account for their azimuthal dependence, we detect weak linear polarization of $^{12}$CO and $^{13}$CO line emission at a $\sim 10σ$ and $\sim 5σ$ significance, respectively. The polarization was detected in the line wings, reaching a peak polarization fraction of $\sim 5\%$ and $\sim 3\%$ for the two molecules between disk radii of 0.5" and 1". The sign of the polarization was found to flip from the blue-shifted side of the emission to the red-shifted side, suggesting a complex, asymmetric polarization morphology. Polarization is not robustly detected for the CS emission; however, a tentative signal, comparable in morphology to that found for the $^{12}$CO and $^{13}$CO emission, is found at a $\lesssim 3σ$ significance. We are able to reconstruct a polarization morphology, consistent with the azimuthally averaged profiles, under the assumption that this is also azimuthally symmetric, which can be compared with future higher-sensitivity observations. △ Less

Submitted 22 September, 2021; v1 submitted 19 September, 2021; originally announced September 2021.

Comments: Accepted by ApJ

Maser emission from the CO envelope of the asymptotic giant branch star W Hydrae

Authors: W. H. T. Vlemmings, T. Khouri, D. Tafoya

Abstract: Observation of CO emission around asymptotic giant branch (AGB) stars is the primary method to determine gas mass-loss rates. While radiative transfer models have shown that molecular levels of CO can become mildly inverted, causing maser emission, CO maser emission has yet to be confirmed observationally. High-resolution observations of the CO emission around AGB stars now have the brightness tem… ▽ More Observation of CO emission around asymptotic giant branch (AGB) stars is the primary method to determine gas mass-loss rates. While radiative transfer models have shown that molecular levels of CO can become mildly inverted, causing maser emission, CO maser emission has yet to be confirmed observationally. High-resolution observations of the CO emission around AGB stars now have the brightness temperature sensitivity to detect possible weak CO maser emission. We used high angular resolution observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA) to observe the small-scale structure of CO $J=3-2$ emission around the oxygen-rich AGB star W Hya. We find CO maser emission amplifying the stellar continuum with an optical depth $τ\approx-0.55$. The maser predominantly amplifies the limb of the star because CO $J=3-2$ absorption from the extended stellar atmosphere is strongest towards the centre of the star. The CO maser velocity corresponds to a previously observed variable component of high-frequency H$_2$O masers and with the OH maser that was identified as the amplified stellar image. This implies that the maser originates beyond the acceleration region and constrains the velocity profile since we find the population inversion primarily in the inner circumstellar envelope. We find that inversion can be explained by the radiation field at 4.6 $μ$m and that the existence of CO maser emission is consistent with the estimated mass-loss rates for W Hya. However, the pumping mechanism requires a complex interplay between absorption and emission lines in the extended atmosphere. Excess from dust in the circumstellar envelope of W Hya is not sufficient to contribute significantly to the required radiation field at 4.6 $μ$m. The interplay between molecular lines that cause the pumping can be constrained by future multi-level CO observations. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: Accepted for publication in A\&A, 14 pages, 14 figures

Journal ref: A&A 654, A18 (2021)

DEATHSTAR: Nearby AGB stars with the Atacama Compact Array II. CO envelope sizes and asymmetries: The S-type stars

Authors: M. Andriantsaralaza, S. Ramstedt, W. H. T. Vlemmings, T. Danilovich, E. De Beck, M. A. T. Groenewegen, S. Höfner, F. Kerschbaum, T. Khouri, M. Lindqvist, M. Maercker, H. Olofsson, G. Quintana-Lacaci, M. Saberi, R. Sahai, A. Zijlstra

Abstract: We aim to constrain the sizes of the CO circumstellar envelopes (CSEs) of 16 S-type stars, along with an additional 7 and 4 CSEs of C-type and M-type AGB stars, respectively. We map the emission from the CO J=2-1 and 3-2 lines observed with the Atacama Compact Array (ACA) and its total power (TP) antennas, and fit with a Gaussian distribution in the uv- and image planes for ACA-only and TP observa… ▽ More We aim to constrain the sizes of the CO circumstellar envelopes (CSEs) of 16 S-type stars, along with an additional 7 and 4 CSEs of C-type and M-type AGB stars, respectively. We map the emission from the CO J=2-1 and 3-2 lines observed with the Atacama Compact Array (ACA) and its total power (TP) antennas, and fit with a Gaussian distribution in the uv- and image planes for ACA-only and TP observations, respectively. The major axis of the fitted Gaussian for the CO(2-1) line data gives a first estimate of the size of the CO-line-emitting CSE. We investigate possible signs of deviation from spherical symmetry by analysing the line profiles, the results from visibility fitting, and by investigating the deconvolved images. The sizes of the CO-line-emitting CSEs of low-mass-loss-rate (low-MLR) S-stars fall between the sizes of the CSEs of C-stars, which are larger, and those of M-stars, which are smaller, as expected because of the differences in their respective CO abundances. The sizes of the low-MLR S-type stars show no dependence on circumstellar density, while a steeper density dependence is observed at high MLR. Furthermore, our results show that the CO CSEs of most of the S-stars in our sample are consistent with a spherically symmetric and smooth outflow. The CO envelope sizes obtained in this paper will be used to constrain detailed radiative transfer modelling to directly determine more accurate MLR estimates for the stars in our sample. For several of our sources that present signs of deviation from spherical symmetry, further high-resolution observations would be necessary to investigate the nature of, and the physical processes behind, these asymmetrical structures. This will provide further insight into the mass-loss process and its related chemistry in S-type AGB stars. △ Less

Submitted 24 August, 2021; originally announced August 2021.

Comments: 8 pages, 2 figures, 4 appendices, accepted in A&A

Journal ref: A&A 653, A53 (2021)

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

Burst timescales and luminosities link young pulsars and fast radio bursts

Authors: K. Nimmo, J. W. T. Hessels, F. Kirsten, A. Keimpema, J. M. Cordes, M. P. Snelders, D. M. Hewitt, R. Karuppusamy, A. M. Archibald, V. Bezukovs, M. Bhardwaj, R. Blaauw, S. T. Buttaccio, T. Cassanelli, J. E. Conway, A. Corongiu, R. Feiler, E. Fonseca, O. Forssen, M. Gawronski, M. Giroletti, M. A. Kharinov, C. Leung, M. Lindqvist, G. Maccaferri , et al. (23 additional authors not shown)

Abstract: Fast radio bursts (FRBs) are extragalactic radio flashes of unknown physical origin. Their high luminosities and short durations require extreme energy densities, like those found in the vicinity of neutron stars and black holes. Studying the burst intensities and polarimetric properties on a wide range of timescales, from milliseconds down to nanoseconds, is key to understanding the emission mech… ▽ More Fast radio bursts (FRBs) are extragalactic radio flashes of unknown physical origin. Their high luminosities and short durations require extreme energy densities, like those found in the vicinity of neutron stars and black holes. Studying the burst intensities and polarimetric properties on a wide range of timescales, from milliseconds down to nanoseconds, is key to understanding the emission mechanism. However, high-time-resolution studies of FRBs are limited by their unpredictable activity levels, available instrumentation and temporal broadening in the intervening ionised medium. Here we show that the repeating FRB 20200120E can produce isolated shots of emission as short as about 60 nanoseconds in duration, with brightness temperatures as high as $3\times 10^{41}$ K (excluding relativistic effects), comparable to `nano-shots' from the Crab pulsar. Comparing both the range of timescales and luminosities, we find that FRB 20200120E observationally bridges the gap between known Galactic young pulsars and magnetars, and the much more distant extragalactic FRBs. This suggests a common magnetically powered emission mechanism spanning many orders of magnitude in timescale and luminosity. In this work, we probe a relatively unexplored region of the short-duration transient phase space; we highlight that there likely exists a population of ultra-fast radio transients at nanosecond to microsecond timescales, which current FRB searches are insensitive to. △ Less

Submitted 29 September, 2021; v1 submitted 24 May, 2021; originally announced May 2021.

Comments: Submitted. Comments welcome

A repeating fast radio burst source in a globular cluster

Authors: F. Kirsten, B. Marcote, K. Nimmo, J. W. T. Hessels, M. Bhardwaj, S. P. Tendulkar, A. Keimpema, J. Yang, M. P. Snelders, P. Scholz, A. B. Pearlman, C. J. Law, W. M. Peters, M. Giroletti, Z. Paragi, C. Bassa, D. M. Hewitt, U. Bach, V. Bezrukovs, M. Burgay, S. T. Buttaccio, J. E. Conway, A. Corongiu, R. Feiler, O. Forssén , et al. (41 additional authors not shown)

Abstract: Fast radio bursts (FRBs) are exceptionally luminous flashes of unknown physical origin, reaching us from other galaxies (Petroff et al. 2019). Most FRBs have only ever been seen once, while others flash repeatedly, though sporadically (Spitler et al. 2016, CHIME/FRB Collaboration et al. 2021). Many models invoke magnetically powered neutron stars (magnetars) as the engines producing FRB emission (… ▽ More Fast radio bursts (FRBs) are exceptionally luminous flashes of unknown physical origin, reaching us from other galaxies (Petroff et al. 2019). Most FRBs have only ever been seen once, while others flash repeatedly, though sporadically (Spitler et al. 2016, CHIME/FRB Collaboration et al. 2021). Many models invoke magnetically powered neutron stars (magnetars) as the engines producing FRB emission (Margalit & Metzger 2018, CHIME/FRB Collaboration et al. 2020). Recently, CHIME/FRB announced the discovery (Bhardwaj et al. 2021) of the repeating FRB 20200120E, coming from the direction of the nearby grand design spiral galaxy M81. Four potential counterparts at other observing wavelengths were identified (Bhardwaj et al. 2021) but no definitive association with these sources, or M81, could be made. Here we report an extremely precise localisation of FRB 20200120E, which allows us to associate it with a globular cluster (GC) in the M81 galactic system and to place it ~2pc offset from the optical center of light of the GC. This confirms (Bhardwaj et al. 2021) that FRB 20200120E is 40 times closer than any other known extragalactic FRB. Because such GCs host old stellar populations, this association strongly challenges FRB models that invoke young magnetars formed in a core-collapse supernova as powering FRB emission. We propose, instead, that FRB 20200120E is a highly magnetised neutron star formed via either accretion-induced collapse of a white dwarf or via merger of compact stars in a binary system (Margalit et al. 2019). Alternative scenarios involving compact binary systems, efficiently formed inside globular clusters, could also be responsible for the observed bursts. △ Less

Submitted 29 September, 2021; v1 submitted 24 May, 2021; originally announced May 2021.

Comments: Submitted. Comments welcome

Tracing the large-scale magnetic field morphology in protoplanetary disks using molecular line polarization

Authors: Boy Lankhaar, Wouter Vlemmings, Per Bjerkeli

Abstract: (abridged) Magnetic fields are fundamental to the accretion dynamics of protoplanetary disks and they likely affect planet formation. Typical methods to study the magnetic field morphology observe the polarization of dust or spectral lines. However, it has recently become clear that dust-polarization in ALMA's spectral regime does not always faithfully trace the magnetic field structure of protopl… ▽ More (abridged) Magnetic fields are fundamental to the accretion dynamics of protoplanetary disks and they likely affect planet formation. Typical methods to study the magnetic field morphology observe the polarization of dust or spectral lines. However, it has recently become clear that dust-polarization in ALMA's spectral regime does not always faithfully trace the magnetic field structure of protoplanetary disks, which leaves spectral line polarization as a promising method for mapping the magnetic field morphologies of such sources. We aim to model the emergent polarization of different molecular lines in the ALMA wavelength regime that are excited in protoplanetary disks. We explore a variety of disk models and molecules to identify those properties that are conducive to the emergence of polarization in spectral lines and may therefore be viably used for magnetic field measurements in protoplanetary disks. We used PORTAL in conjunction with LIME. Together, they allowed us to treat the polarized line radiative transfer of complex three-dimensional physical and magnetic field structures. We present simulations of the emergence of spectral line polarization of different molecules and molecular transitions in the ALMA wavelength regime. We find that molecules that thermalize at high densities, such as HCN, are also the most susceptible to polarization. We find that such molecules are expected to be significantly polarized in protoplanetary disks, while molecules that thermalize at low densities, such as CO, are only significantly polarized in the outer disk regions. We present the simulated polarization maps at a range of inclinations and magnetic field morphologies, and we comment on the observational feasibility of ALMA linear polarization observations of protoplanetary disks. △ Less

Submitted 9 November, 2021; v1 submitted 13 May, 2021; originally announced May 2021.

Comments: Accepted to A&A

Journal ref: A&A 657, A106 (2022)

Heavy-element Rydberg transition line emission from the post-giant-evolution star HD101584

Authors: H. Olofsson, J. H. Black, T. Khouri, W. H. T. Vlemmings, E. M. L. Humphreys, M. Lindqvist, M. Maercker, L. Nyman, S. Ramstedt, D. Tafoya

Abstract: We report the detection of two lines at millimetre wavelengths towards the immediate surroundings of the post-giant and most likely post-common-envelope star HD101584 using high-angular-resolution ALMA observations. The circumstellar environment of this object is rich in different molecular species, but we find no viable identifications in terms of molecular lines. We aim to determine whether or… ▽ More We report the detection of two lines at millimetre wavelengths towards the immediate surroundings of the post-giant and most likely post-common-envelope star HD101584 using high-angular-resolution ALMA observations. The circumstellar environment of this object is rich in different molecular species, but we find no viable identifications in terms of molecular lines. We aim to determine whether or not these lines can be attributed to the Rydberg transitions -- X30alpha and X26alpha -- of neutral atoms of elements heavier than carbon. A simple model in strict local thermodynamic equilibrium for a warm-gas environment of the moderate-temperature star (T_eff about 8500 K) was constructed to corroborate our findings. A geometrically thin, disc-like geometry seen face-on was chosen and a distance of 1 kpc. The observed flux densities of the lines and the continuum at 232 and 354 GHz can be reproduced using 10^(-3) M_sun of gas at a temperature of about 2800 K and a hydrogen density of about 10^(12) cm(-3), assuming solar abundances for the elements. The gas lies within a distance of about 5 au from the star (assuming a distance of 1 kpc). The ionisation fraction is low, about 3x10^(-5). The origin of such a region is not clear, but it may be related to a common-envelope-evolution phase. With these conditions, the line emissions are dominated by Rydberg transitions within the stable isotopes of Mg. A turbulent velocity field in the range 5.5 - 7.5 km s^(-1) is required to fit the Gaussian line shapes. An upper limit to the average magnetic field in the line-emitting region of 1G is set using the Zeeman effect in these lines. We speculate that Rydberg transitions of heavy elements may be an interesting probe for the close-in environments of other moderate-temperature objects like AGB stars, red supergiants, yellow hypergiants, and binaries of various types. △ Less

Submitted 3 May, 2021; originally announced May 2021.

Comments: 19 pages, 6 figures

Multiwavelength Observations of the RV Tauri Variable System U Monocerotis: Long-Term Variability Phenomena That Can Be Explained by Binary Interactions with a Circumbinary Disk

Authors: Laura D. Vega, Keivan G. Stassun, Rodolfo Montez Jr., Tomasz Kamiński, Laurence Sabin, Eric M. Schlegel, Wouter H. T. Vlemmings, Joel H. Kastner, Sofia Ramstedt, Patricia T. Boyd

Abstract: We present X-ray through submillimeter observations of the classical RV Tauri (RVb-type) variable U Mon, a post-AGB binary with a circumbinary disk (CBD). Our SMA observations indicate a CBD diameter of $\lesssim$550 au. Our XMM-Newton observations make U Mon the first RV Tauri variable detected in X-rays. The X-ray emission is characteristic of a hot plasma ($\sim$10 MK), with L… ▽ More We present X-ray through submillimeter observations of the classical RV Tauri (RVb-type) variable U Mon, a post-AGB binary with a circumbinary disk (CBD). Our SMA observations indicate a CBD diameter of $\lesssim$550 au. Our XMM-Newton observations make U Mon the first RV Tauri variable detected in X-rays. The X-ray emission is characteristic of a hot plasma ($\sim$10 MK), with L$_{X}=5\times10^{30}{\rm erg}~{\rm s}^{-1}$, and we consider its possible origin from U Mon, its companion, and/or binary system interactions. Combining DASCH and AAVSO data, we extend the time-series photometric baseline back to the late 1880s and find evidence that U Mon has secular changes that appear to recur on a timescale of $\sim$60 yr, possibly caused by a feature in the CBD. From literature radial velocities we find that the binary companion is a $\sim$2 M$_{\odot}$ A-type main-sequence star. The orientation of the binary's orbit lies along our line of sight ($ω= 95^\circ$), such that apastron corresponds to photometric RVb minima, consistent with the post-AGB star becoming obscured by the near side of the CBD. In addition, we find the size of the inner-CBD hole ($\sim$4.5-9 au) to be comparable to the binary separation, implying that one or both stars may interact with the CBD at apastron. The obscuration of the post-AGB star implicates the companion as the likely source of the enhanced H$α$ observed at RVb minima and of the X-ray emission that may arise from accreted material. △ Less

Submitted 12 March, 2021; originally announced March 2021.

Comments: 19 pages, 9 figures, 4 tables, in aastex62 format, published in ApJ, NASA news feature available upon request. Light curve in Figure 2 available as "Data behind the Figure"

Journal ref: The Astrophysical Journal, 909:138 (14pp), 2021 March 12

A study of radial self-similar non-relativistic MHD outflow models: parameter space exploration and application to the water fountain W43A

Authors: C. Ceccobello, M. H. M. Heemskerk, Y. Cavecchi, W. H. T. Vlemmings, D. Tafoya

Abstract: Outflows, spanning a wide range of dynamical properties and spatial extensions, have now been associated with a variety of accreting astrophysical objects, from supermassive black holes at the core of active galaxies to young stellar objects. The role of such outflows is key to the evolution of the system that generates them, for they extract a fraction of the orbiting material and angular momentu… ▽ More Outflows, spanning a wide range of dynamical properties and spatial extensions, have now been associated with a variety of accreting astrophysical objects, from supermassive black holes at the core of active galaxies to young stellar objects. The role of such outflows is key to the evolution of the system that generates them, for they extract a fraction of the orbiting material and angular momentum from the region close to the central object and release them in the surroundings. The details of the launching mechanism and their impact on the environment are fundamental to understand the evolution of individual sources and the similarities between different types of outflow-launching systems. We solve semi-analytically the non-relativistic, ideal, magnetohydrodynamics (MHD) equations describing outflows launched from a rotating disk threaded with magnetic fields using our new numerical scheme. We present here a parameter study of a large sample of new solutions. We study the different combinations of forces that lead to a successfully launched jet and discuss their global properties. We show how these solutions can be applied to the outflow of the water fountain W43A for which we have observational constraints on magnetic field, density and velocity of the flow at the location of two symmetrical water maser emitting regions. △ Less

Submitted 5 November, 2020; originally announced November 2020.

Comments: 22 pages, 20 figures, the version includes changes based on referee comments

Gravity and rotation drag the magnetic field in high-mass star formation

Authors: H. Beuther, J. D. Soler, H. Linz, Th. Henning, C. Gieser, R. Kuiper, W. Vlemmings, P. Hennebelle, S. Feng, R. Smith, A. Ahmadi

Abstract: The formation of hot stars out of the cold interstellar medium lies at the heart of astrophysical research. Understanding the importance of magnetic fields during star formation remains a major challenge. With the advent of the Atacama Large Millimeter Array, the potential to study magnetic fields by polarization observations has tremendously progressed. However, the major question remains how muc… ▽ More The formation of hot stars out of the cold interstellar medium lies at the heart of astrophysical research. Understanding the importance of magnetic fields during star formation remains a major challenge. With the advent of the Atacama Large Millimeter Array, the potential to study magnetic fields by polarization observations has tremendously progressed. However, the major question remains how much magnetic fields shape the star formation process or whether gravity is largely dominating. Here, we show that for the high-mass star-forming region G327.3 the magnetic field morphology appears to be dominantly shaped by the gravitational contraction of the central massive gas core where the star formation proceeds. We find that in the outer parts of the region, the magnetic field is directed toward the gravitational center of the region. Filamentary structures feeding the central core exhibit U-shaped magnetic field morphologies directed toward the gravitational center as well, again showing the gravitational drag toward the center. The inner part then shows rotational signatures, potentially associated with an embedded disk, and there the magnetic field morphology appears to be rotationally dominated. Hence, our results demonstrate that for this region gravity and rotation are dominating the dynamics and shaping the magnetic field morphology. △ Less

Submitted 12 October, 2020; originally announced October 2020.

Comments: 10 pages, 4 figures, accepted for the Astrophysical Journal, also available at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

Polarization properties of methanol masers

Authors: Daria Dall'Olio, Wouter Vlemmings, Boy Lankhaar, Gabriele Surcis

Abstract: (Abridged) Astronomical masers have been effective tools to study magnetic fields for many years. In particular, methanol can be used to probe different parts of protostars such as accretion discs and outflows, since it produces one of the strongest and the most commonly observed masers in massive star-forming regions. We investigate the polarization properties of selected methanol maser transitio… ▽ More (Abridged) Astronomical masers have been effective tools to study magnetic fields for many years. In particular, methanol can be used to probe different parts of protostars such as accretion discs and outflows, since it produces one of the strongest and the most commonly observed masers in massive star-forming regions. We investigate the polarization properties of selected methanol maser transitions in light of newly calculated methanol Landé g-factors and considering hyperfine components. We compare our results with previous observations and we evaluate the effect of preferred hyperfine pumping and non-Zeeman effects. We run simulations using the radiative transfer code CHAMP. We find a dependence of linear and circular polarization fractions on the hyperfine transitions. Preferred hyperfine pumping can explain some high levels of linear and circular polarization and some of the peculiar features seen in the S-shape of observed V-profiles. Methanol masers are not significantly affected by non-Zeeman effects. Our models show that for methanol maser emission, both the linear and circular polarization percentages depend on which hyperfine transition is masing and the degree to which it is being pumped. Since non-Zeeman effects become more relevant at high values of brightness temperatures, it is important to obtain good estimates of these quantities and on maser beaming angles. Better constraints on the brightness temperature will help in understand about the extent to which non-Zeeman effects contribute to the observed polarization percentages. In order to detect separate hyperfine components, an intrinsic thermal line width significantly smaller than the hyperfine separation is required. △ Less

Submitted 9 October, 2020; originally announced October 2020.

Comments: Accepted for publication in Astronomy & Astrophysics

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

VLTI/PIONIER reveals the close environment of the evolved system HD101584

Authors: J. Kluska, H. Olofsson, H. Van Winckel, T. Khouri, M. Wittkowski, W. J. de Wit, E. M. L. Humphreys, M. Lindqvist, M. Maercker, S. Ramstedt, D. Tafoya, W. H. T. Vlemmings

Abstract: Context: The observed orbital characteristics of post-asymptotic giant branch (post-AGB) and post-red giant branch (post-RGB) binaries are not understood. We suspect that the missing ingredients to explain them probably lie in the continuous interaction of the central binary with its circumstellar environment. Aims: We aim at studying the circumbinary material in these complex systems by investiga… ▽ More Context: The observed orbital characteristics of post-asymptotic giant branch (post-AGB) and post-red giant branch (post-RGB) binaries are not understood. We suspect that the missing ingredients to explain them probably lie in the continuous interaction of the central binary with its circumstellar environment. Aims: We aim at studying the circumbinary material in these complex systems by investigating the connection between the innermost and large-scale structures. Methods: We perform high-angular resolution observations in the near-infrared continuum of HD101584, which has a complex structure as seen at millimeter wavelengths with a disk-like morphology and a bipolar outflow due to an episode of strong binary interaction. To account for the complexity of the target we first perform an image reconstruction and use this result to fit a geometrical model to extract the morphological and thermal features of the environment. Results: The image reveals an unexpected double-ring structure. We interpret the inner ring to be produced by emission from dust located in the plane of the disk and the outer ring to be produced by emission from dust that is located 1.6[D/1kpc] au above the disk plane. The inner ring diameter (3.94[D/1kpc] au), and temperature (T=1540$\pm$10K) are compatible with the dust sublimation front of the disk. The origin of the out-of-plane ring (with a diameter of 7.39[D/1kpc] au and a temperature of 1014$\pm10$K) could be due to episodic ejection or a dust condensation front in the outflow. Conclusion: The observed outer ring is possibly linked with the blue-shifted side of the large scale outflow seen by ALMA and is tracing its launching location to the central star. Such observations give morphological constraints on the ejection mechanism. Additional observations are needed to constrain the origin of the out-of-plane structure. △ Less

Submitted 26 August, 2020; originally announced August 2020.

Comments: Accepted to A&A. 14 pages, 13 figures

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

DEATHSTAR: Nearby AGB stars with the Atacama Compact Array I. CO envelope sizes and asymmetries: A new hope for accurate mass-loss-rate estimates

Authors: S. Ramstedt, W. H. T. Vlemmings, L. Doan, T. Danilovich, M. Lindqvist, M. Saberi, H. Olofsson, E. De Beck, M. A. T. Groenewegen, S. Höfner, J. H. Kastner, F. Kerschbaum, T. Khouri, M. Maercker, R. Montez, G. Quintana-Lacaci, R. Sahai, D. Tafoya, A. Zijlstra

Abstract: This is the first publication of the DEATHSTAR project. The goal of the project is to reduce the uncertainties of observational estimates of mass-loss rates from Asymptotic Giant Branch (AGB) stars. Line emission from 12CO J=2-1 and 3-2 were mapped using the ACA. In this initial analysis, the emission distribution was fit to a Gaussian distribution in the uv-plane. Detailed radiative transfer anal… ▽ More This is the first publication of the DEATHSTAR project. The goal of the project is to reduce the uncertainties of observational estimates of mass-loss rates from Asymptotic Giant Branch (AGB) stars. Line emission from 12CO J=2-1 and 3-2 were mapped using the ACA. In this initial analysis, the emission distribution was fit to a Gaussian distribution in the uv-plane. Detailed radiative transfer analysis will be presented in the future. The axes of the best-fit Gaussian at the line center of the 12CO J=2-1 emission gives a first indication of the size of the emitting region. Furthermore, the fitting results, such as the major and minor axis, center position, and the goodness of fit across both lines, constrain the symmetry of the emission distribution. We find that the CO envelope sizes are, in general, larger for C-type than for M-type AGB stars, which is expected if the CO/H2 ratio is larger in C-type stars. Furthermore, a relation between the 12CO J=2-1 size and circumstellar density is shown that, while in broad agreement with photodissociation calculations, reveals large scatter and systematic differences between the stellar types. The majority of the sources have CO envelopes that are consistent with a spherically symmetric, smooth outflow. For about a third of the sources, indications of strong asymmetries are found. This is consistent with previous interferometric investigations of northern sources. Smaller scale asymmetries are found in a larger fraction of sources. These results for CO envelope radii and shapes can be used to constrain detailed radiative transfer modeling of the same stars so as to determine mass-loss rates that are independent of photodissociation models. For a large fraction of the sources, observations at higher spatial resolution will be necessary to further investigate the complex circumstellar dynamics revealed by our ACA observations. △ Less

Submitted 18 August, 2020; originally announced August 2020.

Comments: 10 pages, 5 figures, 5 appendices, accepted in A&A

ALMA resolves the remarkable molecular jet and rotating wind in the extremely radio-quiet galaxy NGC 1377

Authors: S. Aalto, N. Falstad, S. Muller, K. Wada, J. S. Gallagher, S. König, K. Sakamoto, W. Vlemmings, C. Ceccobello, K. Dasyra, F. Combes, S. García-Burillo, Y. Oya, S. Martín, P. van der Werf, A. S. Evans, J. Kotilainen

Abstract: Submillimetre and millimetre observations are important in probing the properties of the molecular gas and dust around obscured active galactic nuclei (AGNs) and their feedback. With very high-resolution (0."02x0."03 (2x3 pc)) ALMA 345 GHz observations of CO 3-2, HCO$^+$ 4-3, HCN 4-3 $ν_2$=1$f$, and continuum we have studied the molecular outflow and nucleus of the extremely radio-quiet lenticular… ▽ More Submillimetre and millimetre observations are important in probing the properties of the molecular gas and dust around obscured active galactic nuclei (AGNs) and their feedback. With very high-resolution (0."02x0."03 (2x3 pc)) ALMA 345 GHz observations of CO 3-2, HCO$^+$ 4-3, HCN 4-3 $ν_2$=1$f$, and continuum we have studied the molecular outflow and nucleus of the extremely radio-quiet lenticular galaxy NGC1377. The outflow is resolved, revealing a 150 pc long, clumpy, high-velocity, collimated molecular jet. The molecular emission is emerging from the spine of the jet with an average diameter of 3-7 pc. A narrow-angle, rotating molecular wind surrounds the jet and is enveloped by a larger-scale, slower CO-emitting structure. The jet and narrow wind are turbulent ($σ>$40 kms$^{-1}$) and have steep radial gas excitation gradients. The jet shows velocity reversals that we propose are caused by precession, or episodic directional changes. We suggest that an important process powering the outflow is magneto-centrifugal driving. In contrast, the large-scale CO-envelope may be a slow wind, or cocoon that stems from jet-wind interactions. An asymmetric, nuclear r$\sim$2 pc and hot (>180 K) dust structure with a high molecular column density, N(H$_2$)$\sim1.8 \times 10^{24}$ cm$^{-2}$, is detected in continuum and vibrationally excited HCN. Its luminosity is likely powered by a buried AGN. The mass of the supermassive black hole (SMBH) is estimated to $\sim9\times10^6$ M$_\odot$ and the SMBH of NGC1377 appears to be at the end of an intense phase of accretion. The nuclear growth may be fuelled by low-angular momentum gas inflowing from gas ejected in the molecular jet and wind. Such a feedback-loop of cyclic accretion and outflows would be an effective process in growing the nuclear SMBH. This result invites new questions as to SMBH growth processes in obscured, dusty galaxies. △ Less

Submitted 26 July, 2020; v1 submitted 15 July, 2020; originally announced July 2020.

Comments: 16 pages, 12 figures, accepted for publication in Astronomy and Astrophysics. Updated affiliations, added reference

VLBI20-30: a scientific roadmap for the next decade -- The future of the European VLBI Network

Authors: Tiziana Venturi, Zsolt Paragi, Michael Lindqvist, Anna Bartkiewicz, Rob Beswick, Tamara Bogdanović, Walter Brisken, Patrick Charlot, Francisco Colomer, John Conway, Sándor Frey, José Carlos Guirado, Leonid Gurvits, Huib van Langevelde, Andrei Lobanov, John McKean, Raffaella Morganti, Tom Muxlow, Miguel Pérez-Torres, Kazi Rygl, Robert Schulz, Arpad Szomoru, Pablo de Vicente, Tao An, Guillem Anglada , et al. (55 additional authors not shown)

Abstract: This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, dem… ▽ More This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end. △ Less

Submitted 5 July, 2020; originally announced July 2020.

Comments: Editors: Tiziana Venturi, Zsolt Paragi, Michael Lindqvist. "EVN Vision Document (2020)", 7 chapters, 3 appendices, 191 pages, 64 figures. Chapter coordinators and all the contributors are listed in the document

Mapping circumstellar magnetic fields of late-type evolved stars with the Goldreich-Kylafis effect: CARMA observations at $λ1.3$ mm of R Crt and R Leo

Authors: K. -Y. Huang, A. J. Kemball, W. H. T. Vlemmings, S. -P. Lai, L. Yang, I. Agudo

Abstract: Mapping magnetic fields is the key to resolving what remains an unclear physical picture of circumstellar magnetic fields in late-type evolved stars. Observations of linearly polarized emission from thermal molecular line transitions due to the Goldreich-Kylafis (G-K) effect provides valuable insight into the magnetic field geometry in these sources that is complementary to other key studies. In t… ▽ More Mapping magnetic fields is the key to resolving what remains an unclear physical picture of circumstellar magnetic fields in late-type evolved stars. Observations of linearly polarized emission from thermal molecular line transitions due to the Goldreich-Kylafis (G-K) effect provides valuable insight into the magnetic field geometry in these sources that is complementary to other key studies. In this paper, we present the detection of spectral-line polarization from both the thermal $J=2-1$ CO line and the $v=1, J=5-4$ SiO maser line toward two thermal-pulsating (TP-) AGB stars, R Crt and R Leo. The observed fractional linear polarization in the CO emission is measured as $m_l\sim 3.1\%$ and $m_l\sim9.7\%$ for R Crt and R Leo respectively. A circumstellar envelope (CSE) model profile and the associated parameters are estimated and used as input to a more detailed modeling of the predicted linear polarization expected from the G-K effect. The observed thermal line polarization level is consistent with the predicted results from the G-K model for R Crt; additional effects need to be considered for R Leo. △ Less

Submitted 1 July, 2020; originally announced July 2020.

Collisional polarization of molecular ions: a signpost of ambipolar diffusion

Authors: Boy Lankhaar, Wouter Vlemmings

Abstract: Magnetic fields play a role in the dynamics of many astrophysical processes, but they are hard to detect. In a partially ionized plasma, a magnetic field works directly on the ionized medium but not on the neutral medium, which gives rise to a velocity drift between them: ambipolar diffusion. This process is suggested to be important in the process of star formation, but has never been directly ob… ▽ More Magnetic fields play a role in the dynamics of many astrophysical processes, but they are hard to detect. In a partially ionized plasma, a magnetic field works directly on the ionized medium but not on the neutral medium, which gives rise to a velocity drift between them: ambipolar diffusion. This process is suggested to be important in the process of star formation, but has never been directly observed. We introduce a method that could be used to detect ambipolar diffusion and the magnetic field that gives rise to it, where we exploit the velocity drift between the charged and neutral medium. By using a representative classical model of the collision dynamics, we show that molecular ions partially align themselves when a velocity drift is present between the molecular ion and its main collision partner H2. We demonstrate that ambipolar diffusion potently aligns molecular ions in regions denser than their critical density, which subsequently leads to partially polarized emission from these species. We include a model for HCO+ and show that collisional polarization could be detectable for the ambipolar drifts predicted by numerical simulations of the inner protostellar disk regions. The polarization vectors are aligned perpendicular to the magnetic field direction projected on the plane of the sky. △ Less

Submitted 22 June, 2020; originally announced June 2020.

Comments: 5 pages, 2 figures. Accepted and published in A&A

Journal ref: A&A 638, L7 (2020)

Characterizing the radio continuum nature of sources in the massive star-forming region W75N (B)

Authors: A. Rodríguez-Kamenetzky, C. Carrasco-González, J. M. Torrelles, W. H. T. Vlemmings, L. F. Rodríguez, G. Surcis, J. F. Gómez, J. Cantó, C. Goddi, J. S. Kim, S. -W. Kim, N. Añez-López, S. Curiel, H. J. van Langevelde

Abstract: The massive star-forming region W75N~(B) is thought to host a cluster of massive protostars (VLA~1, VLA~2, and VLA~3) undergoing different evolutionary stages. In this work, we present radio continuum data with the highest sensitivity and angular resolution obtained to date in this region, using the VLA-A and covering a wide range of frequencies (4-48~GHz), which allowed us to study the morphology… ▽ More The massive star-forming region W75N~(B) is thought to host a cluster of massive protostars (VLA~1, VLA~2, and VLA~3) undergoing different evolutionary stages. In this work, we present radio continuum data with the highest sensitivity and angular resolution obtained to date in this region, using the VLA-A and covering a wide range of frequencies (4-48~GHz), which allowed us to study the morphology and the nature of the emission of the different radio continuum sources. We also performed complementary studies with multi-epoch VLA data and ALMA archive data at 1.3 mm wavelength. We find that VLA~1 is driving a thermal radio jet at scales of $\approx$0.1 arcsec ($\approx$130 au), but also shows signs of an incipient hyper-compact HII region at scales of $\lesssim$ 1 arcsec ($\lesssim$ 1300~au). VLA~3 is also driving a thermal radio jet at scales of a few tenths of arcsec (few hundred of au). We conclude that this jet is shock-exciting the radio continuum sources Bc and VLA~4 (obscured HH objects), which show proper motions moving outward from VLA~3 at velocities of $\approx$112--118~km/s. We have also detected three new weak radio continuum sources, two of them associated with millimeter continuum cores observed with ALMA, suggesting that these two sources are also embedded YSOs in this massive star-forming region. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: 14 pages, 7 figures. Accepted for publication in MNRAS

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

ALMA reveals the coherence of the magnetic field geometry in OH 231.8+4.2

Authors: L. Sabin, R. Sahai, W. H. T. Vlemmings, Q. Zhang, A. A. Zijlstra, T. Gledhill, M. Huarte-Espinosa, A. F. Pérez Sánchez, E. Lagadec, S. G. Navarro

Abstract: In a continuing effort to investigate the role of magnetic fields in evolved low and intermediate mass stars (principally regarding the shaping of their envelopes), we present new ALMA high resolution polarization data obtained for the nebula OH 231.8+4.2. We found that the polarized emission likely arises from aligned grains in the presence of magnetic fields rather than radiative alignment and s… ▽ More In a continuing effort to investigate the role of magnetic fields in evolved low and intermediate mass stars (principally regarding the shaping of their envelopes), we present new ALMA high resolution polarization data obtained for the nebula OH 231.8+4.2. We found that the polarized emission likely arises from aligned grains in the presence of magnetic fields rather than radiative alignment and self scattering. The ALMA data show well organized electric field orientations in most of the nebula and the inferred magnetic field vectors (rotated by 90 degrees) trace an hourglass morphology centred on the central system of the nebula. One region in the southern part of OH 231.8+4.2 shows a less organized distribution probably due to the shocked environment. These findings, in conjunction with earlier investigations (maser studies and dust emission analysis at other scales and wavelengths) suggest an overall magnetic hourglass located inside a toroidal field. We propose the idea that the magnetic field structure is closely related to the architecture of a magnetic tower and that the outflows were therefore magnetically launched. While the current dynamical effect of the fields might be weak in the equatorial plane principally due to the evolution of the envelope, it would still be affecting the outflows. In that regard, the measurement of the magnetic field at the stellar surface, which is still missing, combined with a full MHD treatment are required to better understand and constrain the events occurring in OH 231.8+4.2. △ Less

Submitted 20 May, 2020; originally announced May 2020.

Comments: 9 pages, 7 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

CO and HCN isotopologue ratios in the outflows of AGB stars

Authors: M. Saberi, H. Olofsson, W. H. T. Vlemmings, E. De Beck, T. Khouri, S. Ramstedt

Abstract: Isotopologue line intensity ratios of circumstellar molecules have been widely used to trace the photospheric elemental isotopic ratios of evolved stars. However, depending on the molecular species and the physical conditions of the environment, the circumstellar isotopologue ratio may deviate considerably from the stellar atmospheric value. In this paper, we aim to examine how the CO and HCN abun… ▽ More Isotopologue line intensity ratios of circumstellar molecules have been widely used to trace the photospheric elemental isotopic ratios of evolved stars. However, depending on the molecular species and the physical conditions of the environment, the circumstellar isotopologue ratio may deviate considerably from the stellar atmospheric value. In this paper, we aim to examine how the CO and HCN abundance ratios vary radially due to chemical reactions in the outflows of AGB stars and the effect of excitation and optical depth on the resulting line intensity ratios. We find that the circumstellar 12CO/13CO can deviate from its atmospheric value by up to 25-94% and 6-60% for C- and O-type CSEs, respectively. We show that variations of the intensity of the ISRF and the gas kinetic temperature can significantly influence the CO isotopologue ratio in the outer CSEs. On the contrary, the H12CN/H13CN ratio is stable for all tested mass-loss rates. The RT modeling shows that the integrated line intensity ratio of CO of different rotational transitions varies significantly for stars with intermediate mass-loss rates due to combined chemical and excitation effects. In contrast, the excitation conditions for the both HCN isotopologues are the same. We demonstrate the importance of using the isotopologue abundance profiles from chemical models as inputs to RT models in the interpretation of isotopologue observations. Previous studies of CO isotopologue ratios are based on multi-transition data for individual sources and it is difficult to estimate the errors in the reported values due to assumptions that are not entirely correct according to this study. If anything, previous studies may have overestimated the circumstellar 12CO/13CO abundance ratio. The use of the HCN as a tracer of C isotope ratios is affected by fewer complicating problems, provided one accounts corrections for high optical depths. △ Less

Submitted 12 May, 2020; originally announced May 2020.

Comments: 14 pages, 11 figures

Journal ref: A&A 638, A99 (2020)

Inner dusty envelope of the AGB stars W~Hydrae, SW\,Virginis, and R~Crateris using SPHERE/ZIMPOL

Authors: T. Khouri, W. H. T. Vlemmings, C. Paladini, C. Ginski, E. Lagadec, M. Maercker, P. Kervella, E. De Beck, L. Decin, A. de Koter, L. B. F. M. Waters

Abstract: Stars with initial masses between $\sim0.8$ and 8~$M_\odot$ present copious mass loss during the asymptotic giant branch (AGB) at the end of their lives. The accepted mass-loss mechanism requires radiation pressure acting on dust grains that form in the extended AGB stellar atmospheres. The details of this process are not yet well understood, however. Using the extreme-adaptive-optics imager and p… ▽ More Stars with initial masses between $\sim0.8$ and 8~$M_\odot$ present copious mass loss during the asymptotic giant branch (AGB) at the end of their lives. The accepted mass-loss mechanism requires radiation pressure acting on dust grains that form in the extended AGB stellar atmospheres. The details of this process are not yet well understood, however. Using the extreme-adaptive-optics imager and polarimeter SPHERE/ZIMPOL, we observed light polarised by grains around W\,Hya, SW\,Vir, and R\,Crt, which have mass-loss rates between 10$^{-7}$ and 10$^{-6}~M_\odot~{\rm yr^{-1}}$. We find the distribution of dust to be asymmetric around the three targets. A biconical morphology is seen for R Crt, with a position angle that is very similar to those inferred from interferometric observations of maser emission and of mid-infrared continuum emission. The cause of the biconical outflow cannot be directly inferred from the ZIMPOL data. The dust grains polarise light more efficiently at 0.65~$μ$m for R\,Crt and SW\,Vir and at 0.82~$μ$m for W\,Hya. This indicates that at the time of the observations, the grains around SW\,Vir and R\,Crt had sizes $< 0.1~μ$m, while those around W\,Hya were larger, with sizes $\gtrsim 0.1~μ$m. The asymmetric distribution of dust around R\,Crt makes the interpretation more uncertain for this star, however. We find that polarised light is produced already from within the visible photosphere of W~Hya, which we reproduce using models with an inner dust shell that is optically thick to scattering. The radial profile of the polarised light observed around W\,Hya reveal a steep dust density profile. We find the wind-acceleration region of W\,Hya to extend to at least $\sim 7~R_\star$, in agreement with theoretical predictions of acceleration up to $\sim 12~R_\star$. △ Less

Submitted 13 March, 2020; originally announced March 2020.

Comments: 16 pages, 9 figures, 3 tables

PORTAL--three-dimensional polarized (sub)millimeter line radiative transfer

Authors: Boy Lankhaar, Wouter Vlemmings

Abstract: Context. Magnetic fields are important to the dynamics of many astrophysical processes and can typically be studied through polarization observations. Polarimetric interferometry capabilities of modern (sub)millimeter telescope facilities have made it possible to obtain detailed velocity resolved maps of molecular line polarization. To properly analyze these for the information they carry regardin… ▽ More Context. Magnetic fields are important to the dynamics of many astrophysical processes and can typically be studied through polarization observations. Polarimetric interferometry capabilities of modern (sub)millimeter telescope facilities have made it possible to obtain detailed velocity resolved maps of molecular line polarization. To properly analyze these for the information they carry regarding the magnetic field, the development of adaptive three-dimensional polarized line radiative transfer models is necessary. Aims. We aim to develop an easy-to-use program to simulate the polarization maps of molecular and atomic (sub)millimeter lines in magnetized astrophysical regions, such as protostellar disks, circumstellar envelopes, or molecular clouds. Methods. By considering the local anisotropy of the radiation field as the only alignment mechanism, we can model the alignment of molecular or atomic species inside a regular line radiative transfer simulation by only making use of the converged output of this simulation. Calculations of the aligned molecular or atomic states can subsequently be used to ray trace the polarized maps of the three-dimensional simulation. Results. We present a three-dimensional radiative transfer code, POlarized Radiative Transfer Adapted to Lines (PORTAL), that can simulate the emergence of polarization in line emission through a magnetic field of arbitrary morphology. Our model can be used in stand-alone mode, assuming LTE excitation, but it is best used when processing the output of regular three-dimensional (nonpolarized) line radiative transfer modeling codes. We present the spectral polarization map of test cases of a collapsing sphere and protoplanetary disk for multiple three-dimensional magnetic field morphologies. △ Less

Submitted 9 March, 2020; originally announced March 2020.

Comments: 13 pages, 8 figures, Accepted for publication in A&A

Journal ref: A&A 636, A14 (2020)

The extended molecular envelope of the asymptotic giant branch star $π^{1}$ Gruis as seen by ALMA II. The spiral-outflow observed at high-angular resolution

Authors: L. Doan, S. Ramstedt, W. H. T. Vlemmings, S. Mohamed, S. Höfner, E. De Beck, F. Kerschbaum, M. Lindqvist, M. Maercker, C. Paladini, M. Wittkowski

Abstract: The AGB star $π^{1}$ Gruis has a known companion (at a separation of ~400 AU) which cannot explain the strong deviations from the spherical symmetry of the CSE. Recently, hydrodynamic simulations of mass transfer in closer binary systems have successfully reproduced the spiral-shaped CSEs found around a handful of sources. There is growing evidence for an even closer, undetected companion complica… ▽ More The AGB star $π^{1}$ Gruis has a known companion (at a separation of ~400 AU) which cannot explain the strong deviations from the spherical symmetry of the CSE. Recently, hydrodynamic simulations of mass transfer in closer binary systems have successfully reproduced the spiral-shaped CSEs found around a handful of sources. There is growing evidence for an even closer, undetected companion complicating the case of $π^{1}$ Gruis further. The improved spatial resolution allows for the investigation of the complex circumstellar morphology and the search for imprints on the CSE of the third component. We have observed the 12CO J=3-2 line emission from $π^{1}$ Gruis using both the compact and extended array of Atacama Large Millimeter/submillimeter Array (ALMA). The interferometric data has furthermore been combined with data from the ALMA total power (TP) array. The imaged brightness distribution has been used to constrain a non-local, non-LTE 3D radiative transfer model of the CSE. The high-angular resolution ALMA data have revealed the first example of a source on the AGB where both a faster bipolar outflow and a spiral pattern along the orbital plane can be seen in the gas envelope. The spiral can be traced in the low- to intermediate velocity, from 13 to 25 km s$^{-1}$, equatorial torus. The largest spiral-arm separation is $\approx$5".5 and consistent with a companion with an orbital period of $\approx$330 yrs and a separation of less than 70 AU. The kinematics of the bipolar outflow is consistent with it being created during a mass-loss eruption where the mass-loss rate from the system increased by at least a factor of 5 during 10-15 yrs. The spiral pattern is the result of an undetected companion. The bipolar outflow is the result of a rather recent mass-loss eruption event. △ Less

Submitted 27 November, 2019; v1 submitted 25 November, 2019; originally announced November 2019.

Comments: 12 pages, 11 figures

VLA cm-wave survey of young stellar objects in the Oph A cluster: constraining extreme UV- and X-ray-driven disk photo-evaporation -- A pathfinder for Square Kilometre Array studies

Authors: A. Coutens, H. B. Liu, I. Jiménez-Serra, T. L. Bourke, J. Forbrich, M. Hoare, L. Loinard, L. Testi, M. Audard, P. Caselli, A. Chacón-Tanarro, C. Codella, J. Di Francesco, F. Fontani, M. Hogerheijde, A. Johansen, D. Johnstone, S. Maddison, O. Panić, L. M. Pérez, L. Podio, A. Punanova, J. M. C. Rawlings, D. Semenov, M. Tazzari , et al. (6 additional authors not shown)

Abstract: Observations of young stellar objects (YSOs) in centimeter bands can probe the continuum emission from growing dust grains, ionized winds, and magnetospheric activity, which are intimately connected to the evolution of protoplanetary disks and the formation of planets. We have carried out sensitive continuum observations toward the Ophiuchus A star-forming region using the Karl G. Jansky Very Larg… ▽ More Observations of young stellar objects (YSOs) in centimeter bands can probe the continuum emission from growing dust grains, ionized winds, and magnetospheric activity, which are intimately connected to the evolution of protoplanetary disks and the formation of planets. We have carried out sensitive continuum observations toward the Ophiuchus A star-forming region using the Karl G. Jansky Very Large Array (VLA) at 10 GHz over a field-of-view of 6$'$ with a spatial resolution of $θ_{maj}$ $\times$ $θ_{min}$ $\sim$ 0.4$''$ $\times$ 0.2$''$. We achieved a 5 $μ$Jy beam$^{-1}$ root-mean-square noise level at the center of our mosaic field of view. Among the eighteen sources we detected, sixteen are YSOs (three Class 0, five Class I, six Class II, and two Class III) and two are extragalactic candidates. We find that thermal dust emission generally contributes less that 30% of the emission at 10 GHz. The radio emission is dominated by other types of emission such as gyro-synchrotron radiation from active magnetospheres, free-free emission from thermal jets, free-free emission from the outflowing photo-evaporated disk material, and/or synchrotron emission from accelerated cosmic-rays in jet or protostellar surface shocks. These different types of emission could not be clearly disentangled. Our non-detections towards Class II/III disks suggest that extreme UV-driven photoevaporation is insufficient to explain the disk dispersal, assuming that the contribution of UV photoevaporating stellar winds to radio flux does not evolve with time. The sensitivity of our data cannot exclude photoevaporation due to X-ray photons as an efficient mechanism for disk dispersal. Deeper surveys with the Square Kilometre Array will be able to provide strong constraints on disk photoevaporation. △ Less

Submitted 8 September, 2019; originally announced September 2019.

Comments: Accepted in A&A

Journal ref: A&A 631, A58 (2019)

Spatio-kinematical model of the collimated molecular outflow in the water-fountain nebula IRAS 16342-3814

Authors: D. Tafoya, G. Orosz, W. H. T. Vlemmings, R. Sahai, A. F. Pérez-Sánchez

Abstract: In this work we aimed to describe the three-dimensional morphology and kinematics of the molecular gas of the water-fountain nebula IRAS 16342-3814. In order to do this, we retrieved data from the ALMA archive to analyse it using a simple spatio-kinematical model. We used the software SHAPE to construct a three-dimensional spatio-kinematical model of the molecular gas in IRAS 16342-3814. By reprod… ▽ More In this work we aimed to describe the three-dimensional morphology and kinematics of the molecular gas of the water-fountain nebula IRAS 16342-3814. In order to do this, we retrieved data from the ALMA archive to analyse it using a simple spatio-kinematical model. We used the software SHAPE to construct a three-dimensional spatio-kinematical model of the molecular gas in IRAS 16342-3814. By reproducing the intensity distribution and position-velocity diagram of the CO emission from the ALMA observations we derived the morphology and velocity field of the gas. We used CO(1-0) data to support the physical interpretation of the model. A spatio-kinematical model that includes a high-velocity collimated outflow embedded within material expanding at relatively lower velocity reproduces the images and position-velocity diagrams from the observations. The high-velocity collimated outflow exhibits deceleration across its length, while the velocity of the surrounding component increases with distance. The morphology of the emitting region; the velocity field and the mass of the gas as function of velocity are in excellent agreement with the properties predicted for a molecular outflow driven by a jet. The timescale of the molecular outflow is estimated to be ~70-100 years. An oscillating pattern was found associated to the high-velocity collimated outflow. The oscillation period of the pattern is T~60-90 years and its opening angle is ~2 degrees. The CO (3-2) emission in IRAS 16342-3814 is interpreted in terms of a jet-driven molecular outflow expanding along an elongated region. The position-velocity diagram and the mass spectrum reveal a feature due to entrained material that is associated to the driving jet. It is likely that the jet in those objects has already disappeared since it is expected to last only for a couple of hundred years. △ Less

Submitted 14 June, 2019; originally announced June 2019.

Comments: 11 pages, 8 figues. Accepted for publication in A&A

Journal ref: A&A 629, A8 (2019)

Characterizing maser polarization: effects of saturation, anisotropic pumping and hyperfine structure

Authors: Boy Lankhaar, Wouter Vlemmings

Abstract: The polarization of masers contains information on the magnetic field strength and direction of the regions they occur in. Many maser polarization observations have been performed over the last 30 years. However, versatile maser polarization models that can aide in the interpretation of these observations are not available. We aim to develop a program suite that can compute the polarization by a m… ▽ More The polarization of masers contains information on the magnetic field strength and direction of the regions they occur in. Many maser polarization observations have been performed over the last 30 years. However, versatile maser polarization models that can aide in the interpretation of these observations are not available. We aim to develop a program suite that can compute the polarization by a magnetic field of any non-paramagnetic maser specie at arbitrarily high maser saturation. Furthermore, we aim to investigate the polarization of masers by non-Zeeman polarizing effects. We aim to present a general interpretive structure for maser polarization observations. We expand existing maser polarization theories of non-paramagnetic molecules and incorporate these in a numerical modeling program suite. We present a modeling program that CHAracterizes Maser Polarization (CHAMP) that can examine the polarization of masers of arbitrarily high maser saturation and high angular momentum. Also, hyperfine multiplicity of the maser-transition can be incorporated. The user is able to investigate non-Zeeman polarizing mechanisms such as anisotropic pumping and polarized incident seed radiation. We present an analysis of the polarization of v = 1 SiO masers and the 22 GHz water maser. We comment on the underlying polarization mechanisms, and also investigate non-Zeeman effects. We identify the regimes where different polarizing mechanisms will be dominant and present the polarization characteristics of the SiO and water masers. From the results of our calculations, we identify markers to recognize alternative polarization mechanisms. △ Less

Submitted 13 May, 2019; originally announced May 2019.

Comments: 67 pages, 27 figures. Accepted to be published in A&A

ALMA reveals the magnetic field evolution in the high-mass star forming complex G9.62+0.19

Authors: Daria Dall'Olio, W. H. T. Vlemmings, M. V. Persson, F. O. Alves, H. Beuther, J. M. Girart, G. Surcis, J. M. Torrelles, H. J. Van Langevelde

Abstract: Context. The role of magnetic fields during the formation of high-mass stars is not yet fully understood, and the processes related to the early fragmentation and collapse are largely unexplored today. The high-mass star forming region G9.62+0.19 is a well known source, presenting several cores at different evolutionary stages. Aims. We determine the magnetic field morphology and strength in the h… ▽ More Context. The role of magnetic fields during the formation of high-mass stars is not yet fully understood, and the processes related to the early fragmentation and collapse are largely unexplored today. The high-mass star forming region G9.62+0.19 is a well known source, presenting several cores at different evolutionary stages. Aims. We determine the magnetic field morphology and strength in the high-mass star forming region G9.62+0.19, to investigate its relation to the evolutionary sequence of the cores. Methods. We use Band 7 ALMA observations in full polarisation mode and we analyse the polarised dust emission. We estimate the magnetic field strength via the Davis-Chandrasekhar-Fermi and the Structure Function methods. Results. We resolve several protostellar cores embedded in a bright and dusty filamentary structure. The polarised emission is clearly detected in six regions. Moreover the magnetic field is oriented along the filament and appears perpendicular to the direction of the outflows. We suggest an evolutionary sequence of the magnetic field, and the less evolved hot core exhibits a magnetic field stronger than the more evolved one. We detect linear polarisation from thermal line emission and we tentatively compared linear polarisation vectors from our observations with previous linearly polarised OH masers observations. We also compute the spectral index, the column density and the mass for some of the cores. Conclusions. The high magnetic field strength and the smooth polarised emission indicate that the magnetic field could play an important role for the fragmentation and the collapse process in the star forming region G9.62+019 and that the evolution of the cores can be magnetically regulated. On average, the magnetic field derived by the linear polarised emission from dust, thermal lines and masers is pointing in the same direction and has consistent strength. △ Less

Submitted 13 May, 2019; v1 submitted 1 May, 2019; originally announced May 2019.

Comments: accepted by A&A, version after language editing

Journal ref: A&A 626, A36 (2019)