Search

IRAS source 19312+1950 (hereafter I19312) is an infrared point source with maser emissions of SiO, H$_2$O, and OH molecules. Although initial observations suggested that I19312 might be an evolved star, its characteristics are not fully consistent with this classification. This study aims to further investigate the nature of I19312 by conducting long-term monitoring of its maser emissions and comparing the results with other known astrophysical objects. We conducted long-term monitoring of SiO, H$_2$O, and OH maser emissions using single-dish radio telescopes. The results were then compared with historical maser data and the characteristics of similar objects to infer the possible origin of I19312. The SiO maser emissions from I19312 were detected over a wide velocity range and exhibited significant time variability. The OH maser lines suggest characteristics of an evolved star, while the H$_2$O maser lines indicate molecular outflows. These features suggest that I19312 could be a candidate for a Water Fountain (WF) star, though there are inconsistencies, such as the large molecular gas mass, that challenge this hypothesis. The possibility of I19312 being a Red Nova Remnant (RNR) is also considered, but this remains speculative due to the lack of direct evidence. The evolutionary stage of I19312 remains unclear, but it shares multiple characteristics with both evolved stars with peculiar properties and RNRs. Further long-term monitoring and high-resolution interferometric observations are required to better constrain the nature of this object., Comment: 36 pages, 17 figures, Accepted for publication in ApJ

We analyse continuum and molecular emission, observed with ALMA, from the dust-enshrouded intermediate-mass AGB star OH 30.1 -0.7. We find a secondary peak in the continuum maps, "feature B", separated by 4.6" from the AGB star, which corresponds to a projected separation of $1.8 \times 10^4$ au, placing a lower limit on the physical separation. This feature is most likely composed of cold dust and is likely to be ejecta associated with the AGB star, though we cannot rule out that it is a background object. The molecular emission we detect includes lines of CO, SiS, CS, SO$_2$, NS, NaCl, and KCl. We find that the NS emission is off centre and arranged along an axis perpendicular to the direction of feature B, indicative of a UV-emitting binary companion (e.g. a G-type main sequence star or hotter), perhaps on an eccentric orbit, contributing to its formation. However, the NaCl and KCl emission constrain the nature of that companion to not be hotter than a late B-type main sequence star. We find relatively warm emission arising from the inner wind and detect several vibrationally excited lines of SiS (3 = 1), NaCl (up to 3 = 4) and KCl (up to 3 = 2), and emission from low energy levels in the mid to outer envelope, as traced by SO$_2$. The CO emission is abruptly truncated around 3.5" or 14,000 au from the continuum peak, suggesting that mass loss at a high rate may have commenced as little as 2800 years ago., Comment: Accepted for publication in MNRAS

Binary interactions have been proposed to explain a variety of circumstellar structures seen around evolved stars, including asymptotic giant branch (AGB) stars and planetary nebulae. Studies resolving the circumstellar envelopes of AGB stars have revealed spirals, discs and bipolar outflows, with shaping attributed to interactions with a companion. For the first time, we have used a combined chemical and dynamical analysis to reveal a highly eccentric and long-period orbit for W Aquilae, a binary system containing an AGB star and a main sequence companion. Our results are based on anisotropic SiN emission, the first detections of NS and SiC towards an S-type star, and density structures observed in the CO emission. These features are all interpreted as having formed during periastron interactions. Our astrochemistry-based method can yield stringent constraints on the orbital parameters of long-period binaries containing AGB stars, and will be applicable to other systems., Comment: This version of the article has been accepted for publication, after peer review (when applicable) and is subject to Springer Nature's AM terms of use, but is not the Version of Record and does not reflect post-acceptance improvements, or any corrections. The Version of Record is available online at: https://doi.org/10.1038/s41550-023-02154-y

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

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

The role of massive ($\geq$ 8M$_{\odot}$) stars in defining the energy budget and chemical enrichment of the interstellar medium in their host galaxy is significant. In this first paper from the Tracing Evolution in Massive Protostellar Objects (TEMPO) project we introduce a colour-luminosity selected (L$_*$ $\sim$ 3$\times10^3$ to 1$\times10^5$ L$_{\odot}$) sample of 38 massive star forming regions observed with ALMA at 1.3mm and explore the fragmentation, clustering and flux density properties of the sample. The TEMPO sample fields are each found to contain multiple fragments (between 2-15 per field). The flux density budget is split evenly (53%-47%) between fields where emission is dominated by a single high flux density fragment and those in which the combined flux density of fainter objects dominates. The fragmentation scales observed in most fields are not comparable with the thermal Jeans length, $\lambda_J$, being larger in the majority of cases, suggestive of some non-thermal mechanism. A tentative evolutionary trend is seen between luminosity of the clump and the `spectral line richness' of the TEMPO fields; with 6.7GHz maser associated fields found to be lower luminosity and more line rich. This work also describes a method of line-free continuum channel selection within ALMA data and a generalised approach used to distinguishing sources which are potentially star-forming from those which are not, utilising interferometric visibility properties., Comment: 42 pages, 17 figures, 4 Tables including 3 Appendices (with additional tables and figures). Accepted for publication in MNRAS

Water and the hydroxyl radical are major constituents of the envelope of O-rich late-type stars. Transitions involving energy levels that are highly excited have been observed in both H$_2$O and OH. These and more recently discovered transitions can now be observed at a high sensitivity and angular resolution with the ALMA Array. Spectra and maps of H$_2$O and OH observed with an angular resolution of 20 to $\sim$200 mas were obtained at two epochs with the ALMA array. Observations with the Compact Array were also used to check for time variability of water transitions. Radiative transfer models of water were revisited to characterize masing conditions and up-to-date chemical models were used for comparison with our observations. Ten rotational transitions of H$_2$O with energies up to 9000 K were observed in various vibrational states. All but one are new detections in space, and from these we have derived accurate rest frequencies. Hyperfine split $\Lambda$-doubling transitions in v = 0, J = 27/2 and 29/2 levels of the $^2\Pi_{3/2}$ state and, $J = 33/2$ and 35/2 of the $^2\Pi_{1/2}$ state of OH with excitation energies up to 8900 K were also observed. Four of these transitions are new detections in space. Combining our measurements with earlier observations of OH, the v = 0 and v = 1 $\Lambda$-doubling frequencies have been improved. Our H$_2$O maps show compact emission and extensions up to twelve stellar radii or more. The 268.149 GHz emission line of water in the v$_2$ = 2 state is time variable, tends to be masing with dominant radiative pumping, and is widely excited. The widespread but weaker 262.898 GHz water line in v$_2$ = 1 also shows signs of maser emission. Emission and absorption of both H$_2$O and OH reveal an infall of matter and complex kinematics influenced by binarity. From our observed column densities, we derived OH/H$_2$O abundance ratios in a few stars., Comment: 41 pages with references and 25 figures in main text. 4 Tables in Appendix A. 10 figures in Appendix B. 2 figures in Appendix C. 14 figures in Appendix D. 5 figures in Appendix E. 2 figures in Appendix F

Aims. Through the ATOMIUM project, based on an ALMA large program, we aim to present a consistent view of a sample of 17 nearby cool evolved stars (Aymptotic Giant Branch and red supergiant stars). Methods. Here we present VLT/SPHERE-ZIMPOL polarimetric maps obtained in the visible of 14 out of the 17 ATOMIUM sources. They were obtained contemporaneously with the ALMA high spatial resolution data. To help interpret the polarized signal, we produced synthetic maps of light scattering by dust, through 3D radiative transfer simulations with the RADMC3D code. Results. The degree of linear polarization (DoLP) observed by ZIMPOL spreads across several optical filters. We infer that it primarily probes dust located just outside of the point spread function, and in or near the plane of the sky, with a total optical depth close to unity in the line of sight, representing only a fraction of the total circumstellar dust. The maximum DoLP ranges from 0.03-0.38 depending on the source, fractions that can be reproduced by our 3D pilot models for grains composed of common dust species. The spatial structure of the DoLP shows a diverse set of shapes. Only for three sources do we note a correlation between the ALMA CO and SiO lines, which trace the gas density, and the DoLP, which traces the dust. Conclusion. The clumpiness of the DoLP and the lack of a consistent correlation between the gas and the dust location show that, in the inner circumstellar environment (CSE), dust formation occurs at very specific sites. This has potential consequences for the derived mass-loss rates and dust-to-gas ratio in the inner region of the CSE. Except for $\pi^1$~Gru and perhaps GY Aql, we do not detect interactions between the circumstellar wind and the hypothesized companions that shape the wind at larger scales. This suggests that the orbits of any other companions are tilted out of the plane of the sky., Comment: Accepted for publication in Astronomy & Astrophysics. 22 pages, 15 figures, 5 tables

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

Context. To better understand the formation of high-mass stars, it is fundamental to investigate how matter accretes onto young massive stars, how it is ejected, and how all this differs from the low-mass case. The massive protocluster G31.41+0.31 is the ideal target to study all these processes because observations at millimeter and centimeter wavelengths have resolved the emission of the Main core into at least four massive dust continuum sources, named A, B, C, and D, within 1" or 0.018 pc, and have identified signatures of infall and several outflows associated with the core. Aims. We study the interplay between infall and outflow in G31.41+0.31 by investigating at a spatial resolution of a few 100 au their properties and their possible impact on the core. Methods. We carried out molecular line observations of typical high-density tracers, such as CH3CN or H2CO, and shock and outflow tracers, such as SiO, with ALMA at 1.4 mm that achieved an angular resolution of 0.09" (340 au). Results. The observations have revealed inverse P-Cygni profiles in CH3CN and H2CO toward the four sources embedded in the Main core, suggesting that all of them are undergoing collapse. The infall rates, estimated from the red-shifted absorption are on the order of 1E-2 Msun/yr. The individual infall rates imply that the accretion timescale of the Main core is an order of magnitude smaller than its rotation timescale. This confirms that rotating toroids such as the G31 Main core are non-equilibrium, transient collapsing structures that need to be constantly replenished with fresh material from a large-scale reservoir. For sources B, C, and D, the infall could be accelerating inside the sources, while for source A, the presence of a second emission component complicates the interpretation. The SiO observations have revealed the presence of at least six outflows in the G31.41+0.31 star-forming region, ..., Comment: 17 pages, 9 figures, 2 tables. Accepted by A&A

This overview paper presents ATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter-submillimeter Array (ALMA). The goal of ATOMIUM is to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical environments. 17 oxygen-rich AGB and RSG stars spanning a range in (circum)stellar parameters and evolutionary phases were observed in a homogeneous observing strategy allowing for an unambiguous comparison. Data were obtained between 213.83 and 269.71 GHz at high (0.025-0.050 arcsec), medium (0.13-0.24 arcsec), and low (about 1 arcsec) angular resolution. The sensitivity per 1.3 km/s channel was 1.5-5 mJy/beam. 13 molecules were designated as primary molecules in the survey: CO, SiO, AlO, AlOH, TiO, TiO2, HCN, SO, SO2, SiS, CS, H2O, and NaCl. The scientific motivation, survey design, sample properties, data reduction, and an overview of the data products are described; and we highlight one scientific result - the wind kinematics of the ATOMIUM sources. The ATOMIUM sources often have a slow wind acceleration, and a fraction of the gas reaches a velocity which can be up to a factor of two times larger than previously reported terminal velocities assuming isotropic expansion, and the wind kinematic profiles establish that the radial velocity described by the momentum equation for a spherical wind structure cannot capture the complexity of the velocity field. In 15 sources, some molecular transitions other than 12CO v=0 J=2-1 reach a higher outflow velocity, with a spatial emission zone that is often greater than 30 stellar radii, but much less than the extent of CO. Binary interaction with a (sub)stellar companion might (partly) explain the non-monotonic behaviour of the projected velocity field.

S-type asymptotic giant branch (AGB) stars are thought to be intermediates in the evolution of oxygen- to carbon-rich AGB stars. The chemical compositions of their circumstellar envelopes are also intermediate, but have not been studied in as much detail as their carbon- and oxygen-rich counterparts. We aim to determine the abundances of AlCl and AlF from rotational lines, which have been observed for the first time towards an S-type AGB star, W Aql. In combination with models based on PACS observations, we aim to update our chemical kinetics network based on these results. We analyse ALMA observations towards W Aql of AlCl in the ground and first two vibrationally excited states and AlF in the ground vibrational state. Using radiative transfer models, we determine the abundances and spatial abundance distributions of Al$^{35}$Cl, Al$^{37}$Cl, and AlF. We also model HCl and HF emission and compare these models to PACS spectra to constrain the abundances of these species. AlCl is found in clumps very close to the star, with emission confined within 0.1$^{\prime\prime}$ of the star. AlF emission is more extended, with faint emission extending 0.2$^{\prime\prime}$ to 0.6$^{\prime\prime}$ from the continuum peak. We find peak abundances, relative to H$_2$, of $1.7\times 10^{-7}$ for Al$^{35}$Cl, $7\times 10^{-8}$ for Al$^{37}$Cl and $1\times 10^{-7}$ for AlF. From the PACS spectra, we find abundances of $9.7\times 10^{-8}$ and $\leq 10^{-8}$, relative to H$_2$, for HCl and HF, respectively. The AlF abundance exceeds the solar F abundance, indicating that fluorine synthesised in the AGB star has already been dredged up to the surface of the star and ejected into the circumstellar envelope. From our analysis of chemical reactions in the wind, we conclude that AlF may participate in the dust formation process, but we cannot fully explain the rapid depletion of AlCl seen in the wind., Comment: Accepted for publication in Astronomy & Astrophysics

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

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

Context. ALMA observations at 1.4 mm and 0.2'' (750au) angular resolution of the Main core in the high-mass star forming region G31.41+0.31 have revealed a puzzling scenario: on the one hand, the continuum emission looks very homogeneous and the core appears to undergo solid-body rotation, suggesting a monolithic core stabilized by the magnetic field; on the other hand, rotation and infall speed up toward the core center, where two massive embedded free-free continuum sources have been detected, pointing to an unstable core having undergone fragmentation. Aims. To establish whether the Main core is indeed monolithic or its homogeneous appearance is due to a combination of large dust opacity and low angular resolution, we carried out millimeter observations at higher angular resolution and different wavelengths. Methods. We carried out ALMA observations at 1.4 mm and 3.5 mm that achieved angular resolutions of 0.1''(375 au) and 0.075'' (280 au), respectively. VLA observations at 7 mm and 1.3 cm at even higher angular resolution, 0.05'' (190 au) and 0.07'' (260 au), respectively, were also carried out to better study the nature of the free-free continuum sources detected in the core. Results. The millimeter continuum emission of the Main core has been clearly resolved into at least four sources, A, B, C, and D, within 1'', indicating that the core is not monolithic. The deconvolved radii of the dust emission of the sources, estimated at 3.5 mm, are 400-500au, their masses range from 15 to 26 Msun, and their number densities are several 1E9 cm-3. Sources A and B, located closer to the center of the core and separated by 750 au, are clearly associated with two free-free continuum sources, likely thermal radio jets, and are the brightest in the core. The spectral energy distribution of these two sources and their masses and sizes are similar and suggest a common origin., Comment: 18 pages,9 figures, 5 tables. Accepted by A&A

Binary interactions dominate the evolution of massive stars, but their role is less clear for low- and intermediate-mass stars. The evolution of a spherical wind from an asymptotic giant branch (AGB) star into a nonspherical planetary nebula (PN) could be due to binary interactions. We observed a sample of AGB stars with the Atacama Large Millimeter/submillimeter Array (ALMA) and found that their winds exhibit distinct nonspherical geometries with morphological similarities to planetary nebulae (PNe). We infer that the same physics shapes both AGB winds and PNe; additionally, the morphology and AGB mass-loss rate are correlated. These characteristics can be explained by binary interaction. We propose an evolutionary scenario for AGB morphologies that is consistent with observed phenomena in AGB stars and PNe., Comment: 19 pages main journal, 97 pages Supplementary Information

A 3D maser model has been used to perform an inverse problem on the light curves from three high-amplitude maser flares, selected on the basis of contemporaneous infra-red observations. Plots derived from the model recover the size of the maser cloud, and two parameters linked to saturation, from three observational properties of the light curve. Recovered sizes are consistent with independent interferometric measurements. Maser objects transition between weak and moderate saturation during a flare., Comment: 5 pages, 5 figures, accepted for publication in MNRAS

We simulate maser flares by varying either the pump rate or the background level of radiation in a 3D model of a maser cloud. We investigate the effect of different cloud shapes, saturation levels and viewpoints. Results are considered for clouds with both uniform and internally variable unsaturated inversion. Pumping and background variations are represented by several different driving functions, some of which are light curves drawn from observations. We summarise the pumping variability results in terms of three observable parameters, the maximum flux density achieved, a variability index and duty cycle. We demonstrate typical ranges of the flux density that may result from viewing an aspherical object from random viewpoints. The best object for a flare is a prolate cloud, viewed close to its long axis and driven from unsaturated conditions to at least modest saturation. Results for variation of the background level are qualitatively different from the variable pumping results in that they tend to produce short intervals of low flux density under conditions of moderate saturation and sufficient variability to be consistent with strong flaring. Variable background models typically have a significantly higher duty cycle than those with variable pumping., Comment: 19 pages, 15 figures, accepted for publication in MNRAS

We investigated, through numerical models, the flaring variability that may arise from the rotation of maser clouds of approximately spheroidal geometry, ranging from strongly oblate to strongly prolate examples. Inversion solutions were obtained for each of these examples over a range of saturation levels from unsaturated to highly saturated. Formal solutions were computed for rotating clouds with many randomly chosen rotation axes, and corresponding averaged maser light curves plotted with statistical information. The dependence of results on the level of saturation and on the degree of deformation from the spherical case were investigated in terms of a variability index and duty cycle. It may be possible to distinguish observationally between flares from oblate and prolate objects. Maser flares from rotation are limited to long timescales (at least a few years) and modest values of the variability index ($\lesssim 100$), and can be aperiodic or quasi-periodic. Rotation is therefore not a good model for H$_2$O variability on timescales of weeks to months, or of truly periodic flares., Comment: 11 pages, 12 figures, accepted for publication in MNRAS

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources., Comment: 23 pages, 13 fig (inc 3 Appendix figures)

As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core., Comment: 28 pages, 16 figures, and 7 tables. Accepted for publication in A&A

We set out the theory and discretization scheme for a new finite-element computer code, written specifically for the simulation of maser sources. The code was used to compute fractional inversions at each node of a 3-D domain for a range of optical thicknesses. Saturation behaviour of the nodes with regard to location and optical depth were broadly as expected. We have demonstrated via formal solutions of the radiative transfer equation that the apparent size of the model maser cloud decreases as expected with optical depth as viewed by a distant observer. Simulations of rotation of the cloud allowed the construction of light-curves for a number of observable quantities. Rotation of the model cloud may be a reasonable model for quasi-periodic variability, but cannot explain periodic flaring., Comment: Accepted for publication in MNRAS, 2018 February 23, 12 pages, 6 figures

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., Comment: 11 pages, accepted by A&A. The abstract here has been shortened

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., Comment: 12 pages, 9 figures, accepted for publication in Astronomy and Astrophysics

We present the analysis of the onset of the new 2010s' OH flaring event detected in the OH ground-state main line at 1665~MHz towards o Ceti and compare its characteristics with those of the 1990s' flaring event. This is based on a series of complementary single-dish and interferometric observations both in OH and H2O obtained with the Nancay Radio telescope (NRT), the Medicina and Effelsberg Telescopes, the European VLBI Network (EVN), and (e)Multi-Element Radio Linked Interferometer Network ((e)MERLIN). We compare the overall characteristics of o Ceti's flaring events with those which have been observed towards other thin-shell Miras, and explore the implication of these events with respect to the standard OH circumstellar-envelope model. The role of binarity in the specific characteristics of o Ceti's flaring events is also investigated. The flaring regions are found to be less than ~400$ +/- 40 mas (i.e., ~40 +/- 4$ AU) either side of o Ceti, with seemingly no preferential location with respect to the direction to the companion Mira B. Contrary to the usual expectation that the OH maser zone is located outside the H2O maser zone, the coincidence of the H2O and OH maser velocities suggests that both emissions arise at similar distances from the star. The OH flaring characteristics of Mira are similar to those observed in various Mira variables before, supporting the earlier results that the regions where the transient OH maser emission occurs are different from the standard OH maser zone., Comment: Preprint, 16 pages, 17 figures, 2 tables, accepted for publication in MNRAS (20 February, 2017)

We report on the measurement of the trigonometric parallaxes of 1612 MHz hydroxyl masers around two asymptotic giant branch stars, WX Psc and OH138.0+7.2, using the NRAO Very Long Baseline Array with in-beam phase referencing calibration. We obtained a 3-sigma upper limit of <=5.3 mas on the parallax of WX Psc, corresponding to a lower limit distance estimate of >~190 pc. The obtained parallax of OH138.0+7.2 is 0.52+/-0.09 mas (+/-18%), corresponding to a distance of 1.9(+0.4,-0.3) kpc, making this the first hydroxyl maser parallax below one milliarcsecond. We also introduce a new method of error analysis for detecting systematic errors in the astrometry. Finally, we compare our trigonometric distances to published phase-lag distances toward these stars and find a good agreement between the two methods., Comment: Preprint, accepted for publication in The Astronomical Journal (January 17, 2017)

A number of mechanisms to understand the periodic class II methanol masers associated with some high-mass star forming regions have been proposed in the past. Two recent proposals, ie. by Parfenov &Sobolev (2014) and Sanna et al. (2015) have been presented in order to explain the periodic masers in sources with light curves similar to the methanol masers in G9.62+0.20E. We evaluate to what extent the proposals and models presented by these authors can explain the light curve of the methanol masers in G9.62+0.20E. It is argued that neither of the proposed mechanisms can reproduce the light curves of the methanol masers in G9.62+0.20E., Comment: 6 pages, 3 figures. To appear in Astronomy and Astrophysics

Distances to AGB stars with optically thick circumstellar shells cannot be determined using optical parallaxes. However, for stars with OH 1612 MHz maser emission emanating from their circumstellar shells, distances can be determined by the phase-lag method. This method combines a linear diameter obtained from a phase-lag measurement with an angular diameter obtained from interferometry. The phase-lag of the variable emission from the back and front sides of the shells has been determined for 20 OH/IR stars in the galactic disk. These measurements are based on a monitoring program with the Nancay radio telescope ongoing for more than 6 years. The interferometric observations are continuing. We estimate that the uncertainties of the distance determination will be ~20%., Comment: 6 pages, 6 figures. Proceedings of the Conference "Why Galaxies Care About AGB Stars III", Vienna 2014, Eds. F. Kerschbaum, J. Hron, and R. Wing

Cool, evolved stars have copious, enriched winds. The structure of these winds and the way they are accelerated is not well known. We need to improve our understanding by studying the dynamics from the pulsating stellar surface to about 10 stellar radii, where radiation pressure on dust is fully effective. Some red supergiants have highly asymmetric nebulae, implicating additional forces. We retrieved ALMA Science Verification data providing images of sub-mm line and continuum emission from VY CMa. This enables us to locate water masers with milli-arcsec precision and resolve the dusty continuum. The 658-, 321- and 325-GHz masers lie in irregular, thick shells at increasing distances from the centre of expansion. For the first time this is confirmed as the stellar position, coinciding with a compact peak offset to the NW of the brightest continuum emission. The maser shells (and dust formation zone) overlap but avoid each other on tens-au scales. Their distribution is broadly consistent with excitation models but the conditions and kinematics appear to be complicated by wind collisions, clumping and asymmetries., Comment: Letter 4 pages, 5 figures plus appendix with 3 figures. Accepted by Astronomy and Astrophysics Letters

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants, Comment: Replaced figures 2 and 3: corrected position SRT. Corrected minor typo in 5.1

The present study aims at characterizing the massive star forming region G35.20N, which is found associated with at least one massive outflow and contains multiple dense cores, one of them recently found associated with a Keplerian rotating disk. We used ALMA to observe the G35.20N region in the continuum and line emission at 350 GHz. The observed frequency range covers tracers of dense gas (e.g. H13CO+, C17O), molecular outflows (e.g. SiO), and hot cores (e.g. CH3CN, CH3OH). The ALMA 870 um continuum emission map reveals an elongated dust structure (0.15 pc long and 0.013 pc wide) perpendicular to the large-scale molecular outflow detected in the region, and fragmented into a number of cores with masses 1-10 Msun and sizes 1600 AU. The cores appear regularly spaced with a separation of 0.023 pc. The emission of dense gas tracers such as H13CO+ or C17O is extended and coincident with the dust elongated structure. The three strongest dust cores show emission of complex organic molecules characteristic of hot cores, with temperatures around 200 K, and relative abundances 0.2-2x10^(-8) for CH3CN and 0.6-5x10^(-6) for CH3OH. The two cores with highest mass (cores A and B) show coherent velocity fields, with gradients almost aligned with the dust elongated structure. Those velocity gradients are consistent with Keplerian disks rotating about central masses of 4-18 Msun. Perpendicular to the velocity gradients we have identified a large-scale precessing jet/outflow associated with core B, and hints of an east-west jet/outflow associated with core A. The elongated dust structure in G35.20N is fragmented into a number of dense cores that may form massive stars. Based on the velocity field of the dense gas, the orientation of the magnetic field, and the regularly spaced fragmentation, we interpret this elongated structure as the densest part of a 1D filament fragmenting and forming massive stars., Comment: 24 pages, 26 figures, accepted for publication in Astronomy and Astrophysics (abstract modified to fit arXiv restrictions)

We report on ALMA observations of continuum and molecular line emission with 0.4" resolution towards the high-mass star forming region G35.20-0.74 N. Two dense cores are detected in typical hot-core tracers, such as CH3CN, which reveal velocity gradients. In one of these cores, the velocity field can be fitted with an almost edge-on Keplerian disk rotating about a central mass of 18 Msun. This finding is consistent with the results of a recent study of the CO first overtone bandhead emission at 2.3mum towards G35.20-0.74 N. The disk radius and mass are >2500 au and 3 Msun. To reconcile the observed bolometric luminosity (3x10^4 Lsun) with the estimated stellar mass of 18 Msun, we propose that the latter is the total mass of a binary system., Comment: 5 pages, 4 figures, accepted for publication in Astronomy and Astrophysics Letters

Convection, pulsation and magnetic fields have all been suggested as mechanisms for the transport of mass and energy from the optical photosphere of red supergiants, out to the region where the stellar wind is launched. We imaged the red supergiant Betelgeuse at 0.06-0.18 arcsec resolution, using e-MERLIN at 5.5--6.0 GHz, with a sensitivity of ~0.01 mJy/beam. Most of the radio emission comes from within an ellipse (0.235x0.218) arcsec^2 (~5x the optical radius), with a flux density of 1.62 mJy, giving an average brightness temperature ~1250 K. This radio photosphere contains two hotspots of 0.53 and 0.79 mJy/beam, separated by 90 milli-arcsec, with brightness temperatures 5400+/-600 K and 3800+/-500 K. Similar hotspots, at more than double the distance from the photosphere of those seen in any other regime, were detected by the less-sensitive `old' MERLIN in 1992, 1995 and 1996 and many exceed the photospheric temperature of 3600 K. Such brightness temperatures are high enough to emanate from pockets of chromospheric plasma. Other possibilities include local shock heating, the convective dredge-up of hot material or exceptionally cool, low density regions, transparent down to the hottest layer at ~40 milliarcsec radius. We also detect an arc 0.2--0.3 arcsec to the SW, brightness temperature ~150 K, in a similar direction to extensions seen on both smaller and larger scales in the infra-red and in CO at mm wavelengths. These preliminary results will be followed by further e-MERLIN, VLA and ALMA observations to help resolve the problem of mass elevation from 1 to 10 R* in red supergiants., Comment: 6pp (incl. appendix), 3 figs Accepted by MNRAS letters 2013 Mar 12

Cool, evolved stars undergo copious mass loss but the details of how the matter is returned to the ISM are still under debate. We investigated the structure and evolution of the wind at 5 to 50 stellar radii from Asymptotic Giant Branch and Red Supergiant stars. 22-GHz water masers around seven evolved stars were imaged using MERLIN, at sub-AU resolution. Each source was observed at between 2 and 7 epochs (several stellar periods). We compared our results with long-term Pushchino single dish monitoring. The 22-GHz emission is located in ~spherical, thick, unevenly filled shells. The outflow velocity doubles between the inner and outer shell limits. Water maser clumps could be matched at successive epochs separated by <2 years for AGB stars, or at least 5 years for RSG. This is much shorter than the decades taken for the wind to cross the maser shell, and comparison with spectral monitoring shows that some features fade and reappear. In 5 sources, most of the matched features brighten or dim in concert from one epoch to the next. One cloud in W Hya was caught in the act of passing in front of a background cloud leading to 50-fold, transient amplification. The masing clouds are 1-2 orders of magnitude denser than the wind average and contain a substantial fraction of the mass loss in this region, with a filling factor <1%. The RSG clouds are ~10x bigger than those round the AGB stars. Proper motions are dominated by expansion, with no systematic rotation. The maser clouds survive for decades (the shell crossing time) but the masers are not always beamed in our direction. Radiative effects cause changes in flux density throughout the maser shells on short timescales. Cloud size is proportional to parent star size; clouds have a similar radius to the star in the 22-GHz maser shell. Stellar properties such as convection cells must determine the clumping scale., Comment: Accepted by A&A 2012 July 10 Main text 29 pages, 62 figures Appendix 44 pages, 23 figures

We present MERLIN observations of OH, water and methanol masers towards the young high mass stellar object IRAS 20126+4104. Emission from the 1665-MHz OH, 22-GHz H_2O and 6.7-GHz CH_3OH masers is detected and all originates very close to the central source. The OH and methanol masers appear to trace part of the circumstellar disk around the central source. The positions and velocities of the OH and methanol masers are consistent with Keplerian rotation around a central mass of ~ 5 Msun. The water masers are offset from the OH and methanol masers and have significantly changed since they were last observed, but still appear to be associated to the outflow from the source. All the OH masers components are circularly polarised, in some cases reaching 100 percent while some OH components also have linear polarisation. We identify one Zeeman pair of OH masers and the splitting of this pair indicates a magnetic field of strength ~ 11 mG within ~ 0.5" (850 AU) of the central source. The OH and methanol maser emission suggest that the disk material is dense, n > 10^6 cm^-3, and warm, T > 125 K and the high abundance of methanol required by the maser emission is consistent with the evaporation of the mantles on dust grains in the disk as a result of heating or shocking of the disk material, Comment: 9 pages, 7 figures and 6 tables

Context. Maser emission is a strong tool for studying high mass star forming regions and their evolutionary stages. OH masers in particular can trace the circumstellar material around protostars and determine their magnetic field strengths at milliarcsecond resolution. Aims. Imaging OH maser mission towards high mass protostellar objects to determine their evolutionary stages and to locate the detected maser emission in the process of high mass star formation. Methods. In 2007, we surveyed OH maser towards 217 high mass protostellar objects to study its presence. In this paper, we present a follow up MERLIN observations of a ground state OH maser emission towards one of these objects, IRAS 19092+0841. Results. Emission from the two OH main spectral lines, 1665 and 1667 MHz, were detected close to the central object. The positions and velocities of the OH maser features have been determined. The masers are distributed over a region of ~ 500 corresponding to 22400 AU (or ~ 0.1 pc) at a distance of 4.48 kpc. The polarization properties of the OH maser features were determined as well. We identify three Zeeman pairs from which we inferred a magnetic field strength of ~ 4:4mG pointing towards the observer. Conclusions. The relatively small velocity spread and the relatively wide spacial distribution of the OH maser features support the suggestion that this object could be in an early evolutionary state before the presence of disk and/or jet/outfows., Comment: 6 pages, 2 figures and 3 tables

We present the first view of the magnetic field structure in the OH shell of the extreme OH/IR star OH 26.5+0.6. MERLIN interferometric observations of this object were obtained in December 1993 in full polarisation, at 1612, 1665 and 1667 MHz. The maser spots show a spheroidal distribution both at 1612 and 1667 MHz, while at 1665 MHz emission from the blue-shifted maser peak is concentrated on the stellar position, and the red-shifted peak emission exhibits a filamentary structure oriented on a SE-NW axis. The linear polarisation in both main lines is rather faint, ranging from 9 to 20% at 1665 MHz and from 0 to 30% at 1667 MHz. At 1612 MHz most maser spots exhibit a similar range of linear polarisation although those in the outermost parts of the envelope reach values as high as 66%. This is particularly apparent in the southern part of the shell. The detailed distribution of the polarisation vectors could only be obtained at 1612 MHz. The polarisation vectors show a highly structured distribution indicative of a poloidal magnetic field inclined by 40-60$^\circ$ to the line of sight. The velocity distribution of the maser spots with respect to the radial distance is well explained by an isotropic outflow at constant velocity in the case of a prolate shaped spheroid envelope, also tilted about 45-65$^\circ$ to the line of sight., Comment: 20 pages, 16 figures, accepted for publication in MNRAS

We have observed the bipolar post-AGB candidate OH 231.8+4.2, using the mid-infrared interferometer MIDI and the infrared camera with the adaptive optics system NACO on the Very Large Telescope. An unresolved core (<200 mas in FWHM) is found at the center of the OH 231.8+4.2 in the 3.8 micron image. This compact source is resolved with the interferometer. We used two 8-meter telescopes with four different baselines, which cover projected baseline lengths from 62 to 47 meters, and projected position angles from 112 to 131 degrees that are almost perpendicular to the bipolar outflow. Fringes from 8 to 9 micron and from 12 to 13.5 micron were clearly detected, whilst the strong silicate self-absorption allows only marginal detection of visibilities between 9 and 12 micron. The fringes from the four baselines consistently show the presence of a compact circumstellar object with an inner radius of 30-40 mas, which is equivalent to 40-50 AU at 1.3 kpc. This clearly shows that the mid-infrared compact source is not the central star (3 AU) but circumstellar material. The measured size of the circumstellar material is consistent with the size of such disks calculated by hydrodynamic models, implying the circumstellar material may have a disk configuration., Comment: Accepted by ApJL High-resolution jpg image is available from http://optik2.mtk.nao.ac.jp/~mikako/photo/oh231_fig1.jpg

We present the preliminary results of the first alignment to milliarsecond accuracy of the 6 GHz maser emission in OH and methanol maser lines toward W3(OH). The identifications of Zeeman pairs allowed us to infer the actual velocity of the OH material and the magnetic field strength at the location where 6.7 GHz methanol maser emission arises., Comment: 2 pages. 2 figures. Proceedings of the 7th European VLBI Network Symposium held in Toledo, Spain on October 12-15, 2004. Editors: R. Bachiller, F. Colomer, J.-F. Desmurs, P. de Vicente (Observatorio Astronomico Nacional), p. 199-200. Needs evn2004.cls

We investigated the main-line spectral evolution with shell thickness of oxygen rich AGB stars. The study is based on a sample of 30 sources distributed along the IRAS colour-colour diagram. The sources were chosen to trace the Miras with thick shells and the whole range of OH/IR stars. The Miras exhibit a 1665 MHz emission strength comparable to that at 1667 MHz. Even though the Miras of the study have quite thick shells, their spectral characteristics in both main lines attest to a strong heterogeneity in their OH shell with, in particular, the presence of significant turbulence and acceleration. The expansion velocity has been found to be about the same at 1665 and 1667 MHz, taking into account a possible velocity turbulence of 1-2km/s at the location of the main-line maser emission. An increase in the intensity ratio 1667/1665 with shell thickness has been found. A plausible explanation for such a phenomenon is that competitive gain in favour of the 1667 MHz line increases when the shell is getting thicker. There is an evolution in the spectral profile shape with the appearance of a substantial inter-peak signal when the shell is getting thicker. Also, inter-peak components are found and can be as strong as the external standard peaks when the shell is very thick. This trend for an increase of the signal in between the two main peaks is thought to be the result of an increase of the saturation with shell thickness. All sources but two - a Mira and an OH/IR star from the lower part of the colour-colour diagram - are weakly polarized. The strong polarization observed for those two particular objects is thought to be the result of perturbations in their shells., Comment: 19 pages, 12 figures, accepted for publication in A&A