Showing 1–50 of 132 results for author: Vlemmings, W H T

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)

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

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

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)

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

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)

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)

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

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

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.

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

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

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)

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)

Resolving the extended stellar atmospheres of Asymptotic Giant Branch stars at (sub-)millimetre wavelengths

Authors: W. H. T. Vlemmings, T. Khouri, H. Olofsson

Abstract: The initial conditions for the mass loss during the asymptotic giant branch (AGB) phase are set in their extended atmospheres, where, among others, convection and pulsation driven shocks determine the physical conditions. High resolution observations of AGB stars at (sub)millimetre wavelengths can now directly determine the morphology, activity, density, and temperature close to the stellar photos… ▽ More The initial conditions for the mass loss during the asymptotic giant branch (AGB) phase are set in their extended atmospheres, where, among others, convection and pulsation driven shocks determine the physical conditions. High resolution observations of AGB stars at (sub)millimetre wavelengths can now directly determine the morphology, activity, density, and temperature close to the stellar photosphere. We used the Atacama Large Millimeter/submillimeter Array (ALMA) high angular resolution observations to resolve the extended atmospheres of four of the nearest AGB stars: W Hya, Mira A, R Dor and R Leo. We interpreted the observations using a parameterised atmosphere model. We resolve all four AGB stars and determine the brightness temperature structure between $1$ and $2$ stellar radii. For W Hya and R Dor we confirm the existence of hotspots with brightness temperatures $>3000$ to $10000$~K. All four stars show deviations from spherical symmetry. We find variations on a timescale of days to weeks, and for R Leo we directly measure an outward motion of the millimetre wavelength surface with a velocity of at least $10.6\pm1.4$~km~s$^{-1}$. For all objects but W Hya we find that the temperature-radius and size-frequency relations require the existence of a (likely inhomogeneous) layer of enhanced opacity. The ALMA observations provide a unique probe of the structure of the extended AGB atmosphere. We find highly variable structures of hotspots and likely convective cells. In the future, these observations can be directly compared to multi-dimensional chromosphere and atmosphere models that determine the temperature, density, velocity, and ionisation structure between the stellar photosphere and the dust formation region. However, our results show that for the best interpretation, both very accurate flux calibration and near-simultaneous observations are essential. △ Less

Submitted 6 May, 2019; v1 submitted 12 April, 2019; originally announced April 2019.

Comments: 18 pages, 13 figures, Accepted to A&A, final version after language editing

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

Stringent limits on the magnetic field strength in the disc of TW Hya: ALMA observations of CN polarisation

Authors: W. H. T. Vlemmings, B. Lankhaar, P. Cazzoletti, C. Ceccobello, D. Dall'Olio, E. F. van Dishoeck, S. Facchini, E. M. L. Humphreys, M. V. Persson, L. Testi, J. P. Williams

Abstract: Despite their importance in the star formation process, measurements of magnetic field strength in proto-planetary discs remain rare. While linear polarisation of dust and molecular lines can give insight into the magnetic field structure, only observations of the circular polarisation produced by Zeeman splitting provide a direct measurement of magnetic field strenghts. One of the most promising… ▽ More Despite their importance in the star formation process, measurements of magnetic field strength in proto-planetary discs remain rare. While linear polarisation of dust and molecular lines can give insight into the magnetic field structure, only observations of the circular polarisation produced by Zeeman splitting provide a direct measurement of magnetic field strenghts. One of the most promising probes of magnetic field strengths is the paramagnetic radical CN. Here we present the first Atacama Large Millimeter/submillimeter Array (ALMA) observations of the Zeeman splitting of CN in the disc of TW~Hya. The observations indicate an excellent polarisation performance of ALMA, but fail to detect significant polarisation. An analysis of eight individual CN hyperfine components as well as a stacking analysis of the strongest (non-blended) hyperfine components yields the most stringent limits obtained so far on the magnetic field strength in a proto-planetary disc. We find that the vertical component of the magnetic field $|B_z|<0.8$~mG ($1σ$ limit). We also provide a $1σ$ toroidal field strength limit of $<30$~mG. These limits rule out some of the earlier accretion disc models, but remain consistent with the most recent detailed models with efficient advection. We detect marginal linear polarisation from the dust continuum, but the almost purely toroidal geometry of the polarisation vectors implies that his is due to radiatively aligned grains. △ Less

Submitted 2 April, 2019; originally announced April 2019.

Comments: 10 pages, 5 figures, accepted for publication in A&A letters

Detection of highly excited OH towards AGB stars: A new probe of shocked gas in the extended atmospheres

Authors: T. Khouri, L. Velilla-Prieto, E. De Beck, W. H. T. Vlemmings, H. Olofsson, B. Lankhaar, J. H. Black, A. Baudry

Abstract: We report the detection and investigate the properties of high-excitation lambda-doubling line emission of hydroxyl (OH) detected towards three asymptotic giant branch (AGB) stars (W Hya, R Dor, and IK Tau) using ALMA. The OH lines are produced very close to the central stars and seem optically thin and with no maser effect. We analyse the molecular excitation using population diagrams and find ro… ▽ More We report the detection and investigate the properties of high-excitation lambda-doubling line emission of hydroxyl (OH) detected towards three asymptotic giant branch (AGB) stars (W Hya, R Dor, and IK Tau) using ALMA. The OH lines are produced very close to the central stars and seem optically thin and with no maser effect. We analyse the molecular excitation using population diagrams and find rotational temperatures of $\sim 2500$ K and column densities of $\sim 10^{19}$ cm$^{-2}$ for both W Hya and R Dor. For W Hya, we observe emission from vibrationally excited H2O arising from the same region as the OH emission. Moreover, CO $v = 1; J = 3 - 2$ emission also shows a brightness peak in the same region. Considering optically thin emission and the rotational temperature derived for OH, we find a CO column density $\sim 15$ times higher than that of OH, within an area of ($92 \times 84$) mas$^2$ centred on the OH emission peak. These results should be considered tentative because of the simple methods employed. The observed OH line frequencies differ significantly from the predicted transition frequencies in the literature, and provide the possibility of using OH lines observed in AGB stars to improve the accuracy of the Hamiltonian used for the OH molecule. We predict stronger OH lambda-doubling lines at millimetre wavelengths than those we detected. These lines will be a good probe of shocked gas in the extended atmosphere and are possibly even suitable as probes of the magnetic field in the atmospheres of close-by AGB stars through the Zeeman effect. △ Less

Submitted 1 March, 2019; v1 submitted 28 February, 2019; originally announced February 2019.

Comments: Published in Astronomy and Astrophysics Letters

Journal ref: A&A 623, L1 (2019)

EVN observations of 6.7 GHz methanol maser polarization in massive star-forming regions IV. Magnetic field strength limits and structure for 7 additional sources

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

Abstract: MHD simulations show that the magnetic field can drive molecular outflows during the formation of massive protostars. The best probe to observationally measure both the morphology and the strength of this magnetic field at scales of 10-100 au is maser polarization. We measure the direction of magnetic fields at mas resolution around a sample of massive star-forming regions to determine whether the… ▽ More MHD simulations show that the magnetic field can drive molecular outflows during the formation of massive protostars. The best probe to observationally measure both the morphology and the strength of this magnetic field at scales of 10-100 au is maser polarization. We measure the direction of magnetic fields at mas resolution around a sample of massive star-forming regions to determine whether there is a relation between the orientation of the magnetic field and of the outflows. In addition, by estimating the magnetic field strength via the Zeeman splitting measurements, the role of magnetic field in the dynamics of the massive star-forming region is investigated. We selected a flux-limited sample of 31 massive star-forming regions to perform a statistical analysis of the magnetic field properties with respect to the molecular outflows characteristics. We report the linearly and circularly polarized emission of 6.7 GHz CH3OH masers towards seven massive star-forming regions of the total sample with the EVN. The sources are: G23.44-0.18, G25.83-0.18, G25.71-0.04, G28.31-0.39, G28.83-0.25, G29.96-0.02, and G43.80-0.13. We identified a total of 219 CH3OH maser features, 47 and 2 of which showed linearly and circularly polarized emission, respectively. We measured well-ordered linear polarization vectors around all the massive YSOs and Zeeman splitting towards G25.71-0.04 and G28.83-0.25. Thanks to recent theoretical results, we were able to provide lower limits to the magnetic field strength from our Zeeman splitting measurements. We further confirm (based on ~80% of the total flux-limited sample) that the magnetic field on scales of 10-100 au is preferentially oriented along the outflow axes. The estimated magnetic field strength of |B_{||}|>61 mG and >21 mG towards G25.71-0.04 and G28.83-0.2, respectively, indicates that it dominates the dynamics of the gas in both regions. △ Less

Submitted 21 February, 2019; originally announced February 2019.

Comments: 19 pages, 8 figures, 12 tables, accepted by Astronomy & Astrophysics

Journal ref: A&A 623, A130 (2019)

HD101584: Circumstellar characteristics and evolutionary status

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

Abstract: We have performed a study of the characteristics of the circumstellar environment of the binary object HD101584, that provides information on a likely evolutionary scenario. We have obtained and analysed ALMA observations, complemented with observations using APEX, of a large number of molecular lines. An analysis of the spectral energy distribution has also been performed. Emissions from 12 molec… ▽ More We have performed a study of the characteristics of the circumstellar environment of the binary object HD101584, that provides information on a likely evolutionary scenario. We have obtained and analysed ALMA observations, complemented with observations using APEX, of a large number of molecular lines. An analysis of the spectral energy distribution has also been performed. Emissions from 12 molecular species (not counting isotopologues) have been observed, and most of them mapped with angular resolutions in the range 0.1" to 0.6". Four circumstellar components are identified: i) a central compact source of size 0.15", ii) an expanding equatorial density enhancement (a flattened density distribution in the plane of the orbit) of size 3", iii) a bipolar high-velocity outflow (150 km/s), and iv) an hourglass structure. The outflow is directed almost along the line of sight. There is evidence of a second bipolar outflow. The mass of the circumstellar gas is 0.5[D/1 kpc]^2 Msun, about half of it lies in the equatorial density enhancement. The dust mass is 0.01[D/1 kpc]^2 Msun, and a substantial fraction of this is in the form of large-sized, up to 1 mm, grains. The estimated kinetic age of the outflow is 770[D/1 kpc] yr. The kinetic energy and the scalar momentum of the accelerated gas are estimated to be 7x10^(45)[D/1 kpc]^2 erg and 10^(39)[D/1 kpc]^2 g cm/s, respectively. We provide good evidence that the binary system HD101584 is in a post-common-envelope-evolution phase, that ended before a stellar merger. Isotope ratios combined with stellar mass estimates suggest that the primary star's evolution was terminated already on the first red giant branch (RGB). Most of the energy required to drive the outflowing gas was probably released when material fell towards the companion. △ Less

Submitted 6 February, 2019; originally announced February 2019.

Comments: Accepted for publication in A&A

High-resolution observations of gas and dust around Mira using ALMA and SPHERE/ZIMPOL

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

Abstract: The outflows of oxygen-rich asymptotic giant branch (AGB) stars are thought to be driven by radiation pressure by photon scattering on grains with sizes of tenths of microns. The details of the formation of dust in the extended atmospheres of these stars and the mass-loss process is still not well understood. We obtained quasi-simultaneous observations of the AGB star Mira using ALMA and ZIMPOL to… ▽ More The outflows of oxygen-rich asymptotic giant branch (AGB) stars are thought to be driven by radiation pressure by photon scattering on grains with sizes of tenths of microns. The details of the formation of dust in the extended atmospheres of these stars and the mass-loss process is still not well understood. We obtained quasi-simultaneous observations of the AGB star Mira using ALMA and ZIMPOL to probe the distribution of gas and large dust grains, respectively. The polarized light images show dust around Mira~A, the companion (Mira~B) and in a trail that connects the two sources. ALMA reveals that dust around Mira~A is contained in a high-gas-density region with a significant fraction of the polarized light arising from its edge. We constrained the gas density, temperature, and velocity within a few stellar radii from the star by modelling the CO $v=1, J=3-2$ line. We find a mass $(\sim 3.8 \pm 1.3) \times 10^{-4}~M_\odot$ to be contained between the stellar millimetre photosphere, $R^{\rm 338~GHz}_\star$, and $4~R^{\rm 338~GHz}_\star$. Our best-fit models with lower masses also reproduce the $^{13}$CO $v=0, J=3-2$ line emission from this region. We find TiO$_2$ and AlO abundances corresponding to 4.5% and $< 0.1$% of the total titanium and aluminium expected for a solar-composition gas. The low abundance of AlO allows for efficient Al depletion into dust already very close to the star, as expected from thermal dust emission observations and theoretical calculations of Mira variables. We constrain the presence of aluminium oxide grains based on the scattered light observations and our gas-phase model. We find that aluminium oxide grains can account for a significant fraction of the total aluminium atoms in this region only if the grains have sizes $\lesssim 0.02~μ$m. This is an order of magnitude smaller than the maximum sizes predicted by dust-formation and wind-driving models. △ Less

Submitted 9 October, 2018; originally announced October 2018.

Comments: 16 pages, 10 figures

Journal ref: A&A 620, A75 (2018)

CO envelope of the symbiotic star R Aqr seen by ALMA

Authors: S. Ramstedt, S. Mohamed, T. Olander, W. H. T. Vlemmings, T. Khouri, S. Liljegren

Abstract: The symbiotic star R Aqr is part of a small sample of binary AGB stars observed with ALMA. The sample stars are: R Aqr, Mira, W Aql, and pi1 Gru. The sample covers a range in binary separation and wind properties, where R Aqr is the source with the smallest separation. The R Aqr binary pair consists of an M-type AGB star and a white dwarf at a separation of 45 mas, equivalent to about 10 AU at 218… ▽ More The symbiotic star R Aqr is part of a small sample of binary AGB stars observed with ALMA. The sample stars are: R Aqr, Mira, W Aql, and pi1 Gru. The sample covers a range in binary separation and wind properties, where R Aqr is the source with the smallest separation. The R Aqr binary pair consists of an M-type AGB star and a white dwarf at a separation of 45 mas, equivalent to about 10 AU at 218 pc. The aim of the ALMA study is to investigate the dependence of the wind shaping on the binary separation and to provide constraints for hydrodynamical binary interaction models. R Aqr is particularly interesting as the source with the smallest separation and a complex circumstellar environment that is strongly affected by the interaction between the two stars and by the high-energy radiation resulting from this interaction and from the hot white dwarf companion. The CO(J=3-2) line emission has been observed at 0.5" spatial resolution. The CO envelope around the binary pair is marginally resolved, showing what appears to be a rather complex distribution. The outer radius of the CO emitting region is estimated from the data and found to be about a factor of 10 larger than previously thought. This implies an average mass-loss rate during the past 100 yr of 2x10-7 Msun/yr, a factor of 45 less than previous estimates. The channel maps are presented and the molecular gas distribution is discussed and set into the context of what was previously known about the system. Additional molecular line emission detected within the bandwidth covered by the ALMA observations is also presented. Because of the limited extent of the emission, firm conclusions about the dynamical evolution of the system will have to wait for higher spatial resolution observations. However, the data presented here support the assumption that the mass-loss rate from the Mira star strongly varies and is focused on the orbital plane. △ Less

Submitted 13 June, 2018; originally announced June 2018.

Comments: 6 pages, 4 figures, accepted in A&A

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

Molecular line study of the S-type AGB star W Aquilae. ALMA observations of CS, SiS, SiO and HCN

Authors: Magdalena Brunner, Taissa Danilovich, Sofia Ramstedt, Ivan Marti-Vidal, Elvire De Beck, Wouter H. T. Vlemmings, Michael Lindqvist, Franz Kerschbaum

Abstract: Context. With the outstanding spatial resolution and sensitivity of the Atacama Large Millimeter/sub-millimeter Array (ALMA), molecular gas other than the abundant CO can be observed and resolved in circumstellar envelopes (CSEs) around evolved stars, such as the binary S-type Asymptotic Giant Branch (AGB) star W Aquilae. Aims. We aim to constrain the chemical composition of the CSE and determine… ▽ More Context. With the outstanding spatial resolution and sensitivity of the Atacama Large Millimeter/sub-millimeter Array (ALMA), molecular gas other than the abundant CO can be observed and resolved in circumstellar envelopes (CSEs) around evolved stars, such as the binary S-type Asymptotic Giant Branch (AGB) star W Aquilae. Aims. We aim to constrain the chemical composition of the CSE and determine the radial abundance distribution, the photospheric peak abundance, and isotopic ratios of a selection of chemically important molecular species in the innermost CSE of W Aql. The derived parameters are put into the context of the chemical evolution of AGB stars and are compared with theoretical models. Methods. We employ one-dimensional radiative transfer modeling - with the accelerated lambda iteration (ALI) radiative transfer code - of the radial abundance distribution of a total of five molecular species (CS, SiS, 30SiS, 29SiO and H13CN) and determine the best fitting model parameters based on high-resolution ALMA observations as well as archival single-dish observations. The additional advantage of the spatially resolved ALMA observations is that we can directly constrain the radial profile of the observed line transitions from the observations. Results. We derive abundances and e-folding radii for CS, SiS, 30SiS, 29SiO and H13CN and compare them to previous studies, which are based only on unresolved single-dish spectra. Our results are in line with previous results and are more accurate due to resolution of the emission regions. △ Less

Submitted 5 June, 2018; originally announced June 2018.

Journal ref: A&A 617, A23 (2018)

Magnetic field in a young circumbinary disk

Authors: Felipe O. Alves, Josep M. Girart, Marco Padovani, Daniele Galli, Gabriel A. P. Franco, Paola Caselli, Wouter H. T. Vlemmings, Qizhou Zhang, Helmut Wiesemeyer

Abstract: We use polarization observations of a circumbinary disk to investigate how the polarization properties change at distinct frequency bands. Our goal is to discern the main mechanism responsible for the polarization through comparison between our observations and model predictions. We used ALMA to perform full polarization observations at 97.5 GHz, 233 GHz and 343.5 GHz. The target is the Class I ob… ▽ More We use polarization observations of a circumbinary disk to investigate how the polarization properties change at distinct frequency bands. Our goal is to discern the main mechanism responsible for the polarization through comparison between our observations and model predictions. We used ALMA to perform full polarization observations at 97.5 GHz, 233 GHz and 343.5 GHz. The target is the Class I object BHB07-11, which is the youngest object in the Barnard 59 protocluster. Complementary VLA observations at 34.5 GHz revealed a binary system within the disk. We detect an extended and structured polarization pattern remarkably consistent among all three bands. The distribution of polarized intensity resembles a horseshoe shape with polarization angles following this morphology. From the spectral index between bands 3 and 7, we derive a dust opacity index $β\sim 1$ consistent with maximum grain sizes larger than expected to produce self-scattering polarization in each band. The polarization morphology do not match predictions from self-scattering. On the other hand, marginal correspondence is seen between our maps and predictions from radiation field assuming the brightest binary component as main radiation source. Molecular line data from BHB07-11 indicates disk rotation. We produced synthetic polarization maps from a rotating magnetized disk model assuming combined poloidal and toroidal magnetic field components. The magnetic field vectors (i. e., the polarization vectors rotated by 90$\degr$) are better represented by a model with poloidal magnetic field strength about 3 times the toroidal one. The similarity of our polarization patterns among the three bands provides a strong evidence against self-scattering and radiation fields. On the other hand, our data are reasonably well reproduced by a model of disk with toroidal magnetic field components slightly smaller than poloidal ones. △ Less

Submitted 6 August, 2018; v1 submitted 1 May, 2018; originally announced May 2018.

Comments: 8 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

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

Detection of CI line emission towards the oxygen-rich AGB star omi Cet

Authors: M. Saberi, W. H. T. Vlemmings, E. De Beck, R. Montez, S. Ramstedt

Abstract: We present the detection of neutral atomic carbon CI(3 P1 - 3 P0 ) line emission towards omi Cet. This is the first time that CI is detected in the envelope around an oxygen-rich M-type asymptotic giant branch (AGB) star. We also confirm the previously tentative CI detection around V Hya, a carbon-rich AGB star. As one of the main photodissociation products of parent species in the circumstellar e… ▽ More We present the detection of neutral atomic carbon CI(3 P1 - 3 P0 ) line emission towards omi Cet. This is the first time that CI is detected in the envelope around an oxygen-rich M-type asymptotic giant branch (AGB) star. We also confirm the previously tentative CI detection around V Hya, a carbon-rich AGB star. As one of the main photodissociation products of parent species in the circumstellar envelope (CSE) around evolved stars, CI can be used to trace sources of ultraviolet (UV) radiation in CSEs. The observed flux density towards omi Cet can be reproduced by a shell with a peak atomic fractional abundance of $2.4 \times 10^{-5}$ predicted based on a simple chemical model where CO is dissociated by the interstellar radiation field. However, the CI emission is shifted by $\sim$ 4 km/s from the stellar velocity. Based on this velocity shift, we suggest that the detected CI emission towards omi Cet potentially arises from a compact region near its hot binary companion. The velocity shift could, therefore, be the result of the orbital velocity of the binary companion around omi Cet. In this case, the CI column density is estimated to be $1.1 \times 10^{19}$ cm$^{-2}$. This would imply that strong UV radiation from the companion and/or accretion of matter between two stars is most likely the origin of the CI enhancement. However, this hypothesis can be confirmed by high-angular resolution observations. △ Less

Submitted 26 April, 2018; originally announced April 2018.

Journal ref: A&A 612, L11 (2018)

Rotation of the asymptotic giant branch star R Doradus

Authors: W. H. T. Vlemmings, T. Khouri, E. De Beck, H. Olofsson, G. Garcia-Segura, E. Villaver, A. Baudry, E. M. L. Humphreys, M. Maercker, S. Ramstedt

Abstract: High resolution observations of the extended atmospheres of asymptotic giant branch (AGB) stars can now directly confront the theories that describe stellar mass loss. Using Atacama Large Millimeter/submillimeter Array (ALMA) high angular resolution ($30\times42$~mas) observations we have, for the first time, resolved stellar rotation of an AGB star, R~Dor. We measure an angular rotation velocity… ▽ More High resolution observations of the extended atmospheres of asymptotic giant branch (AGB) stars can now directly confront the theories that describe stellar mass loss. Using Atacama Large Millimeter/submillimeter Array (ALMA) high angular resolution ($30\times42$~mas) observations we have, for the first time, resolved stellar rotation of an AGB star, R~Dor. We measure an angular rotation velocity of $ω_R\sin{i}=(3.5\pm0.3)\times10^{-9}$~rad~s$^{-1}$ which indicates a rotational velocity of $|\upsilon_{\rm rot}\sin{i}|=1.0\pm0.1$~km~s$^{-1}$ at the stellar surface ($R_*=31.2$~mas at $214$~GHz). The rotation axis projected on the plane of the sky has a position angle $Φ=7\pm6^\circ$. We find that the rotation of R Dor is two orders of magnitude faster than expected for a solitary AGB star that will have lost most of its angular momentum. Its rotational velocity is consistent with angular momentum transfer from a close companion. As a companion has not been directly detected we thus suggest R~Dor has a low-mass, close-in, companion. The rotational velocity approaches the critical velocity, set by the local sound speed in the extended envelope, and is thus expected to affect the mass loss characteristics of R~Dor. △ Less

Submitted 12 April, 2018; originally announced April 2018.

Comments: 12 pages, 15 figures, Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 613, L4 (2018)

Vibrationally excited water emission at 658 GHz from evolved stars

Authors: A. Baudry, E. M. L. Humphreys, F. Herpin, K. Torstensson, W. H. T. Vlemmings, A. M. S. Richards, M. D. Gray, C. De Breuck, M. Olberg

Abstract: Several rotational transitions of water have been identified toward evolved stars in the ground vibrational state as well as in the first excited state of the bending mode. In the latter vibrational state of water, the 658 GHz J = 1_1,0-1_0,1 rotational transition is often strong and seems to be widespread in late-type stars. Our main goals are to better characterize the nature of the 658 GHz emis… ▽ More Several rotational transitions of water have been identified toward evolved stars in the ground vibrational state as well as in the first excited state of the bending mode. In the latter vibrational state of water, the 658 GHz J = 1_1,0-1_0,1 rotational transition is often strong and seems to be widespread in late-type stars. Our main goals are to better characterize the nature of the 658 GHz emission, compare the velocity extent of the 658 GHz emission with SiO maser emission to help locate the water layers and, more generally, investigate the physical conditions prevailing in the excited water layers of evolved stars. Another goal is to identify new 658 GHz emission sources and contribute in showing that this emission is widespread in evolved stars. Eleven evolved stars were extracted from our mini-catalog of existing and potential 658 GHz sources for observations with the APEX telescope equipped with the SEPIA receiver. The 13CO J=6-5 line was placed in the same receiver sideband for simultaneous observation with the 658 GHz line of water. We have compared the 658 GHz line properties with our H2O radiative transfer models in stars and we have compared the velocity ranges of the 658 GHz and SiO J=2-1, v=1 maser lines. All stars show 658 GHz emission with a peak flux density in the range 50-70 Jy to 2000-3000 Jy. We have shown that the 658 GHz line is masing and we found that the 658 GHz velocity extent tends to be correlated with that of the SiO maser suggesting that both emission lines are excited in circumstellar layers close to the central star. Broad and stable line profiles are observed at 658 GHz. This could indicate maser saturation although we have tentatively provided first information on time variability at 658 GHz. △ Less

Submitted 7 November, 2017; originally announced November 2017.

Comments: 9 pages, 5 figures. Accepted in A&A

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

The extended molecular envelope of the asymptotic giant branch star $π^{1}$ Gruis as seen by ALMA I. Large-scale kinematic structure and CO excitation properties

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

Abstract: The S-type asymptotic giant branch (AGB) star $π^{1}$ Gruis has a known companion at a separation of $\approx$400 AU. The envelope structure, including an equatorial torus and a fast bipolar outflow, is rarely seen in the AGB phase and is particularly unexpected in such a wide binary system. Therefore a second, closer companion has been suggested, but the evidence is not conclusive. The new ALMA… ▽ More The S-type asymptotic giant branch (AGB) star $π^{1}$ Gruis has a known companion at a separation of $\approx$400 AU. The envelope structure, including an equatorial torus and a fast bipolar outflow, is rarely seen in the AGB phase and is particularly unexpected in such a wide binary system. Therefore a second, closer companion has been suggested, but the evidence is not conclusive. The new ALMA $^{12}$CO and $^{13}$CO $J$=3-2 data, together with previously published $^{12}$CO $J$=2-1 data from the Submillimeter Array (SMA), and the $^{12}$CO $J$=5-4 and $J$=9-8 lines observed with Herschel/Heterodyne Instrument for the Far-Infrared (HIFI), is modeled with the 3D non-LTE radiative transfer code SHAPEMOL. The data analysis clearly confirms the torus-bipolar structure. The 3D model of the CSE that satisfactorily reproduces the data consists of three kinematic components: a radially expanding torus with velocity slowly increasing from 8 to 13 km s$^{-1}$ along the equator plane; a radially expanding component at the center with a constant velocity of 14 km s$^{-1}$; and a fast, bipolar outflow with velocity proportionally increasing from 14 km s$^{-1}$ at the base up to 100 km s$^{-1}$ at the tip, following a linear radial dependence. The results are used to estimate an average mass-loss rate during the creation of the torus of 7.7$\times$10$^{-7}$ M$_{\odot}$ yr$^{-1}$. The total mass and linear momentum of the fast outflow are estimated at 7.3$\times$10$^{-4}$ M$_{\odot}$ and 9.6$\times$10$^{37}$ g cm s$^{-1}$, respectively. The momentum of the outflow is in excess (by a factor of about 20) of what could be generated by radiation pressure alone, in agreement with recent findings for more evolved sources. The best-fit model also suggests a $^{12}$CO/$^{13}$CO abundance ratio of 50. Possible shaping scenarios for the gas envelope are discussed △ Less

Submitted 27 September, 2017; originally announced September 2017.

Journal ref: A&A, Volume 605, September 2017

The circumstellar envelope around the S-type AGB star W Aql Effects of an eccentric binary orbit

Authors: S. Ramstedt, S. Mohamed, W. H. T. Vlemmings, T. Danilovich, M. Brunner, E. De Beck, E. M. L. Humphreys, M. Lindqvist, M. Maercker, H. Olofsson, F. Kerschbaum, G. Quintana-Lacaci

Abstract: The CO(J=3-2) emission from the CSE of the binary S-type AGB star W Aql has been observed at subarcsecond resolution using ALMA. The aim of this paper is to investigate the wind properties of the AGB star and to analyse how the known companion has shaped the CSE. The average mass-loss rate during the creation of the detected CSE is estimated through modelling, using the ALMA brightness distributio… ▽ More The CO(J=3-2) emission from the CSE of the binary S-type AGB star W Aql has been observed at subarcsecond resolution using ALMA. The aim of this paper is to investigate the wind properties of the AGB star and to analyse how the known companion has shaped the CSE. The average mass-loss rate during the creation of the detected CSE is estimated through modelling, using the ALMA brightness distribution and previously published single-dish measurements as observational constraints. The ALMA observations are presented and compared to the results from a 3D smoothed particle hydrodynamics (SPH) binary interaction model with the same properties as the W Aql system and with two different orbital eccentricities. Three-dimensional radiative transfer modelling is performed and the response of the interferometer is modelled and discussed. The estimated average mass-loss rate of W~Aql agrees with previous results. The size of the emitting region is consistent with photodissociation models. The CO(J=3-2) emission is dominated by a smooth component overlayed with two weak arc patterns with different separations. The larger pattern is predicted by the binary interaction model with separations of 10" and therefore likely due to the known companion. It is consistent with a binary orbit with low eccentricity. The smaller separation pattern is asymmetric and coincides with the dust distribution, but the separation timescale (200 yrs) is not consistent with any known process of the system. The separation of the known companions of the system is large enough to not have a very strong effect on the circumstellar morphology. The density contrast across the envelope of a binary with an even larger separation will not be easily detectable, even with ALMA, unless the orbit is strongly asymmetric or the AGB star has a much larger mass-loss rate. △ Less

Submitted 21 September, 2017; originally announced September 2017.

Comments: 10 pages, 8 figures

Journal ref: A&A 605, A126 (2017)

Methanol masers reveal the magnetic field of the high-mass protostar IRAS 18089-1732

Authors: D. Dall'Olio, W. H. T. Vlemmings, G. Surcis, H. Beuther, B. Lankhaar, M. V. Persson, A. M. S. Richards, E. Varenius

Abstract: Context. The importance of the magnetic field in high-mass-star formation is not yet fully clear and there are still many open questions concerning its role in the accretion processes and generation of jets and outflows. In the past few years, masers have been successfully used to probe the magnetic field morphology and strength at scales of a few au around massive protostars, by measuring linear… ▽ More Context. The importance of the magnetic field in high-mass-star formation is not yet fully clear and there are still many open questions concerning its role in the accretion processes and generation of jets and outflows. In the past few years, masers have been successfully used to probe the magnetic field morphology and strength at scales of a few au around massive protostars, by measuring linear polarisation angles and Zeeman splitting. The massive protostar IRAS 18089-1732 is a well studied high-mass-star forming region, showing a hot core chemistry and a disc-outflow system. Previous SMA observations of polarised dust revealed an ordered magnetic field oriented around the disc of IRAS 18089-1732. Aims. We want to determine the magnetic field in the dense region probed by 6.7 GHz methanol maser observations and compare it with observations in dust continuum polarisation, to investigate how the magnetic field in the compact maser region relates to the large-scale field around massive protostars. Methods. We reduced MERLIN observations at 6.7 GHz of IRAS 18089-1732 and we analysed the polarised emission by methanol masers. Results. Our MERLIN observations show that the magnetic field in the 6.7 GHz methanol maser region is consistent with the magnetic field constrained by the SMA dust polarisation observations. A tentative detection of circularly polarised line emission is also presented. Conclusions. We found that the magnetic field in the maser region has the same orientation as in the disk. Thus the large-scale field component, even at the au scale of the masers, dominates over any small-scale field fluctuations. We obtained, from the circular polarisation tentative detection, a field strength along the line of sight of 5.5 mG which appeared to be consistent with the previous estimates. △ Less

Submitted 9 August, 2017; originally announced August 2017.

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

Journal ref: A&A 607, A111 (2017)

Rings and filaments: The remarkable detached CO shell of U Antliae

Authors: F. Kerschbaum, M. Maercker, M. Brunner, M. Lindqvist, H. Olofsson, M. Mecina, E. De Beck, M. A. T. Groenewegen, E. Lagadec, S. Mohamed, C. Paladini, S. Ramstedt, W. H. T. Vlemmings, M. Wittkowski

Abstract: Aims. Our goal is to characterize the intermediate age, detached shell carbon star U Antliae morphologically and physically in order to study the mass-loss evolution after a possible thermal pulse. Methods. High spatial resolution ALMA observations of unprecedented quality in thermal CO lines allow us to derive first critical spatial and temporal scales and constrain modeling efforts to estimate m… ▽ More Aims. Our goal is to characterize the intermediate age, detached shell carbon star U Antliae morphologically and physically in order to study the mass-loss evolution after a possible thermal pulse. Methods. High spatial resolution ALMA observations of unprecedented quality in thermal CO lines allow us to derive first critical spatial and temporal scales and constrain modeling efforts to estimate mass-loss rates for both the present day as well as the ejection period of the detached shell. Results. The detached shell is remarkably thin, overall spherically symmetric, and shows a barely resolved filamentary substructure possibly caused by instabilities in the interaction zone of winds with different outflow velocities. The expansion age of the detached shell is of the order of 2700 years and its overall width indicates a high expansion-velocity and high mass-loss period of only a few hundred years at an average mass-loss rate of $\approx$ 10$^{-5}$ $M_\odot$ yr$^{-1}$. The post-high-mass-loss-rate-epoch evolution of U Ant shows a significant decline to a substantially lower gas expansion velocity and a mass-loss rate amounting to 4$\times$10$^{-8}$ $M_\odot$ yr$^{-1}$, at present being consistent with evolutionary changes as predicted for the period between thermal pulses. △ Less

Submitted 9 August, 2017; originally announced August 2017.

Journal ref: A&A 605, A116 (2017)

Molecular outflow launched beyond the disk edge

Authors: Felipe O. Alves, Josep M. Girart, Paola Caselli, Gabriel A. P. Franco, Bo Zhao, Wouter H. T. Vlemmings, Marc G. Evans, Luca Ricci

Abstract: One of the long-standing problems of star formation is the excess of angular momentum of the parent molecular cloud. In the classical picture, a fraction of angular momentum of the circumstellar material is removed by the magneto-centrifugally driven disk wind that is launched from a wide region throughout the disk. In this work, we investigate the kinematics in the envelope-disk transition zone o… ▽ More One of the long-standing problems of star formation is the excess of angular momentum of the parent molecular cloud. In the classical picture, a fraction of angular momentum of the circumstellar material is removed by the magneto-centrifugally driven disk wind that is launched from a wide region throughout the disk. In this work, we investigate the kinematics in the envelope-disk transition zone of the Class I object BHB07-11, in the B59 core. For this purpose, we used the Atacama Large Millimeter/submillimeter Array in extended configuration to observe the thermal dust continuum emission ($λ_0 \sim$ 1.3 mm) and molecular lines (CO, C$^{18}$O and H$_2$CO), which are suitable tracers of disk, envelope, and outflow dynamics at a spatial resolution of $\sim 30$ AU. We report a bipolar outflow that was launched at symmetric positions with respect to the disk ($\sim$80~AU in radius), but was concentrated at a distance of 90--130~AU from the disk center. The two outflow lobes had a conical shape and the gas inside was accelerating. The large offset of the launching position coincided with the landing site of the infall material from the extended spiral structure (seen in dust) onto the disk. This indicates that bipolar outflows are efficiently launched within a narrow region outside the disk edge. We also identify a sharp transition in the gas kinematics across the tip of the spiral structure, which pinpoints the location of the so-called centrifugal barrier. △ Less

Submitted 5 July, 2017; originally announced July 2017.

Comments: 5 pages, 5 figures, Accepted for publication in A&A Letters

Journal ref: A&A 603, L3 (2017)

First detection of methanol towards a post-AGB object, HD101584

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

Abstract: The circumstellar environments of objects on the asymptotic giant branch and beyond are rich in molecular species. Nevertheless, methanol has never been detected in such an object, and is therefore often taken as a clear signpost for a young stellar object. However, we report the first detection of CH3OH in a post-AGB object, HD101584, using ALMA. Its emission, together with emissions from CO, SiO… ▽ More The circumstellar environments of objects on the asymptotic giant branch and beyond are rich in molecular species. Nevertheless, methanol has never been detected in such an object, and is therefore often taken as a clear signpost for a young stellar object. However, we report the first detection of CH3OH in a post-AGB object, HD101584, using ALMA. Its emission, together with emissions from CO, SiO, SO, CS, and H2CO, comes from two extreme velocity spots on either side of the object where a high-velocity outflow appears to interact with the surrounding medium. We have derived molecular abundances, and propose that the detected molecular species are the effect of a post-shock chemistry where circumstellar grains play a role. We further provide evidence that HD101584 was a low-mass, M-type AGB star. △ Less

Submitted 26 June, 2017; originally announced June 2017.

Simultaneous 183 GHz H2O Maser and SiO Observations Towards Evolved Stars Using APEX SEPIA Band 5

Authors: E. M. L. Humphreys, K. Immer, M. D. Gray, E. De Beck, W. H. T. Vlemmings, A. Baudry, A. M. S. Richards, M. Wittkowski, K. Torstensson, C. De Breuck, P. Moller, S. Etoka, M. Olberg

Abstract: We investigate the use of 183 GHz H2O masers for characterization of the physical conditions and mass loss process in the circumstellar envelopes of evolved stars. We used APEX SEPIA Band 5 to observe the 183 GHz H2O line towards 2 Red Supergiant and 3 Asymptotic Giant Branch stars. Simultaneously, we observed lines in 28SiO v0, 1, 2 and 3, and for 29SiO v0 and 1. We detected the 183 GHz H2O line… ▽ More We investigate the use of 183 GHz H2O masers for characterization of the physical conditions and mass loss process in the circumstellar envelopes of evolved stars. We used APEX SEPIA Band 5 to observe the 183 GHz H2O line towards 2 Red Supergiant and 3 Asymptotic Giant Branch stars. Simultaneously, we observed lines in 28SiO v0, 1, 2 and 3, and for 29SiO v0 and 1. We detected the 183 GHz H2O line towards all the stars with peak flux densities greater than 100 Jy, including a new detection from VY CMa. Towards all 5 targets, the water line had indications of being due to maser emission and had higher peak flux densities than for the SiO lines. The SiO lines appear to originate from both thermal and maser processes. Comparison with simulations and models indicate that 183 GHz maser emission is likely to extend to greater radii in the circumstellar envelopes than SiO maser emission and to similar or greater radii than water masers at 22, 321 and 325 GHz. We speculate that a prominent blue-shifted feature in the W Hya 183 GHz spectrum is amplifying the stellar continuum, and is located at a similar distance from the star as mainline OH maser emission. From a comparison of the individual polarizations, we find that the SiO maser linear polarization fraction of several features exceeds the maximum fraction allowed under standard maser assumptions and requires strong anisotropic pumping of the maser transition and strongly saturated maser emission. The low polarization fraction of the H2O maser however, fits with the expectation for a non-saturated maser. 183 GHz H2O masers can provide strong probes of the mass loss process of evolved stars. Higher angular resolution observations of this line using ALMA Band 5 will enable detailed investigation of the emission location in circumstellar envelopes and can also provide information on magnetic field strength and structure. △ Less

Submitted 7 April, 2017; originally announced April 2017.

Comments: 11 pages, accepted by A&A. The abstract here has been shortened

The Coldest Place in the Universe: Probing the Ultra-Cold Outflow and Dusty Disk in the Boomerang Nebula

Authors: R. Sahai, W. H. T. Vlemmings, L-A. Nyman

Abstract: Our Cycle 0 ALMA observations confirmed that the Boomerang Nebula is the coldest known object in the Universe, with a massive high-speed outflow that has cooled significantly below the cosmic background temperature. Our new CO 1-0 data reveal heretofore unseen distant regions of this ultra-cold outflow, out to $\gtrsim120,000$ AU. We find that in the ultra-cold outflow, the mass-loss rate (dM/dt)… ▽ More Our Cycle 0 ALMA observations confirmed that the Boomerang Nebula is the coldest known object in the Universe, with a massive high-speed outflow that has cooled significantly below the cosmic background temperature. Our new CO 1-0 data reveal heretofore unseen distant regions of this ultra-cold outflow, out to $\gtrsim120,000$ AU. We find that in the ultra-cold outflow, the mass-loss rate (dM/dt) increases with radius, similar to its expansion velocity ($V$) - taking $V\propto r$, we find $dM/dt \propto r^{0.9-2.2}$. The mass in the ultra-cold outflow is $\gtrsim3.3$ Msun, and the Boomerang's main-sequence progenitor mass is $\gtrsim4$ Msun. Our high angular resolution ($\sim$0".3) CO J=3-2 map shows the inner bipolar nebula's precise, highly-collimated shape, and a dense central waist of size (FWHM) $\sim$1740 AU$\times275$ AU. The molecular gas and the dust as seen in scattered light via optical HST imaging show a detailed correspondence. The waist shows a compact core in thermal dust emission at 0.87-3.3 mm, which harbors $(4-7)\times10^{-4}$ Msun~of very large ($\sim$mm-to-cm sized), cold ($\sim20-30$ K) grains. The central waist (assuming its outer regions to be expanding) and fast bipolar outflow have expansion ages of $\lesssim1925$ yr and $\le1050$ yr: the "jet-lag" (i.e., torus age minus the fast-outflow age) in the Boomerang supports models in which the primary star interacts directly with a binary companion. We argue that this interaction resulted in a common-envelope configuration while the Boomerang's primary was an RGB or early-AGB star, with the companion finally merging into the primary's core, and ejecting the primary's envelope that now forms the ultra-cold outflow. △ Less

Submitted 25 May, 2017; v1 submitted 20 March, 2017; originally announced March 2017.

Comments: accepted ApJ, 12 Apr, 2017

Magnetically aligned dust and SiO maser polarization in the envelope of the red supergiant VY CMa

Authors: W. H. T. Vlemmings, T. Khouri, I. Marti-Vidal, D. Tafoya, A. Baudry, S. Etoka, E. M. L. Humphreys, T. J. Jones, A. Kemball, E. O'Gorman, A. F. Perez-Sanchez, A. M. S. Richards

Abstract: We use Atacama Large Millimeter/submillimeter Array Band 5 science verification observations of the red supergiant VY CMa to study the polarization of SiO thermal/masers lines and dust continuum at ~1.7 mm wavelength. We analyse both linear and circular polarization and derive the magnetic field strength and structure, assuming the polarization of the lines originates from the Zeeman effect, and t… ▽ More We use Atacama Large Millimeter/submillimeter Array Band 5 science verification observations of the red supergiant VY CMa to study the polarization of SiO thermal/masers lines and dust continuum at ~1.7 mm wavelength. We analyse both linear and circular polarization and derive the magnetic field strength and structure, assuming the polarization of the lines originates from the Zeeman effect, and that of the dust originates from aligned dust grains. We also discuss other effects that could give rise to the observed polarization. We detect, for the first time, significant polarization (~3%) of the circumstellar dust emission at millimeter wavelengths. The polarization is uniform with an electric vector position angle of $\sim8^\circ$. Varying levels of linear polarization are detected for the J=4-3 28SiO v=0, 1, 2, and 29SiO v=0, 1 lines, with the strongest polarization fraction of ~30% found for the 29SiO v=1 maser. The linear polarization vectors rotate with velocity, consistent with earlier observations. We also find significant (up to ~1%) circular polarization in several lines, consistent with previous measurements. We conclude that the detection is robust against calibration and regular instrumental errors, although we cannot yet fully rule out non-standard instrumental effects. Emission from magnetically aligned grains is the most likely origin of the observed continuum polarization. This implies that the dust is embedded in a magnetic field >13 mG. The maser line polarization traces the magnetic field structure. The magnetic field in the gas and dust is consistent with an approximately toroidal field configuration, but only higher angular resolution observations will be able to reveal more detailed field structure. If the circular polarization is due to Zeeman splitting, it indicates a magnetic field strength of ~1-3 Gauss, consistent with previous maser observations. △ Less

Submitted 11 April, 2017; v1 submitted 6 March, 2017; originally announced March 2017.

Comments: 12 pages, 9 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 603, A92 (2017)