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

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.

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

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.

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_{\nu}} \propto {\nu}^{1.42}$ between 0.7 mm and 1.4 cm, which we associate with chromosphere-dominated emission, while the flux density ${S_{\nu}} \propto {\nu}^{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.

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.

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.

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~$\mu$m for R\,Crt and SW\,Vir and at 0.82~$\mu$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~\mu$m, while those around W\,Hya were larger, with sizes $\gtrsim 0.1~\mu$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$.

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

The AGB star $\pi^{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 $\pi^{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 $\pi^{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.

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

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.

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

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.

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.

Ultraviolet (UV) photodissociation of CO controls the abundances and distribution of CO and its photodissociation products. This significantly influences the gas-phase chemistry in the circumstellar envelope (CSE) around evolved stars. A better understanding of CO photodissociation in outflows also provides a more precise estimate of mass-loss rates. We aim to update the CO photodissociation rate in an expanding spherical envelope assuming that the interstellar radiation field (ISRF) photons penetrate through the envelope. This will allow us to precisely estimate the CO abundance distributions in CSEs around evolved stars. We used the most recent CO spectroscopic data to precisely calculate the depth dependency of the photodissociation rate of each CO dissociating line. We calculated the CO self- and mutual-shielding functions in an expanding envelope. We investigated the dependence of the CO profile on the five fundamental parameters mass-loss rate, the expansion velocity, the CO initial abundance, the CO excitation temperature, and the strength of the ISRF. Our derived CO envelope size is smaller than the commonly used radius derived by Mamon et al. 1988. The difference between results varies from 1% to 39% and depends on the H2 and CO densities of the envelope. We list two fitting parameters for a large grid of models to estimate the CO abundance distribution. We demonstrate that the CO envelope size can differ between outflows with the same effective content of CO, but different CO abundance, mass-loss rate, and the expansion velocity as a consequence of differing amounts of shielding by H2 and CO. We find that the abundance of CO close to the star and the outflow density both can have a significant effect on the size of the molecular envelope. We also demonstrate that modest variations in the ISRF can cause measurable differences in the envelope extent.

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\sigma$ limit). We also provide a $1\sigma$ 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.

Models of the chemical evolution of the interstellar medium, galaxies, and the Universe rely on our understanding of the amounts and chemical composition of the material returned by stars and supernovae. Stellar yields are obtained from stellar-evolution models, which currently lack predictive prescriptions of stellar mass loss, although it significantly affects stellar lifetimes, nucleosynthesis, and chemical ejecta. Galaxy properties are derived from observations of the integrated light of bright member stars. Stars in the late stages of their evolution are among the infrared-brightest objects in galaxies. An unrealistic treatment of the mass-loss process introduces significant uncertainties in galaxy properties derived from their integrated light. We describe current efforts and future needs and opportunities to characterize AGB outflows: driving mechanisms, outflow rates, underlying fundamental physical and chemical processes such as dust grain formation, and dependency of these on metallicity.

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

There is ample evidence for strong magnetic fields in the envelopes of (Post-)Asymptotic Giant Branch (AGB) stars as well as supergiant stars. The origin and role of these fields are still unclear. This paper updates the current status of magnetic field observations around AGB, post-AGB stars and describes their possible role during these stages of evolution. The discovery of magnetically aligned dust around a supergiant star is also highlighted. In our search for the origin of the magnetic fields, recent observations show the signatures of possible magnetic activity and rotation, indicating that the magnetic fields might be intrinsic to the AGB stars.

Evolved stars dominate galactic spectra, enrich the galactic medium, expand to change their planetary systems, eject winds of a complex nature, produce spectacular nebulae and illuminate them, and transfer material between binary companions. While doing this, they fill the HR diagram with diagnostic loops that write the story of late stellar evolution. Evolved stars sometimes release unfathomable amounts of energy in neutrinos, light, kinetic flow, and gravitational waves. During these late-life times, stars evolve complexly, with expansion, convection, mixing, pulsation, mass loss. Some processes have virtually no spatial symmetries, and are poorly addressed with low-resolution measurements and analysis. Even a "simple" question as how to model mass loss resists solution. However, new methods offer increasingly diagnostic tools. Astrometry reveals populations and groupings. Pulsations/oscillations support study of stellar interiors. Optical/radio interferometry enable 2-3d imagery of atmospheres and shells. Bright stars with rich molecular spectra and velocity fields are a ripe opportunity for imaging with high spatial and spectral resolution, giving insight into the physics and modeling of later stellar evolution.

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.

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.

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.

With the publication of Gaia DR2, 1.3 billion stars now have public parallax and proper motion measurements. In this contribution, we compare the results for sources that have both optical and radio measurements, focusing on circumstellar masers. For these large, variable and bright AGB stars, the VLBI astrometry results can be more robust. Moreover, there are a number of applications where VLBI astrometry provides unique data for studying stellar populations and Galactic structure. The BeSSel project not only provides parallax and proper motions at much larger distances than Gaia can reach, but it also uniquely samples the spiral arms of the Galaxy. The evolved stars in the BAaDE sample can potentially constrain the dynamics and stellar content of the inner bulge and bar of the Milky Way, not reachable in the optical.

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~\mu$m. This is an order of magnitude smaller than the maximum sizes predicted by dust-formation and wind-driving models.

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 $\beta \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.

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.

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.

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.

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 $\omega_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 $\Phi=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.

Magnetic fields play an important role during star formation. Direct magnetic field strength observations have proven specifically challenging in the extremely dynamic protostellar phase. Because of their occurrence in the densest parts of star forming regions, masers, through polarization observations, are the main source of magnetic field strength and morphology measurements around protostars. Of all maser species, methanol is one of the strongest and most abundant tracers of gas around high-mass protostellar disks and in outflows. However, as experimental determination of the magnetic characteristics of methanol has remained largely unsuccessful, a robust magnetic field strength analysis of these regions could hitherto not be performed. Here we report a quantitative theoretical model of the magnetic properties of methanol, including the complicated hyperfine structure that results from its internal rotation. We show that the large range in values of the Land\'{e} g-factors of the hyperfine components of each maser line lead to conclusions which differ substantially from the current interpretation based on a single effective g-factor. These conclusions are more consistent with other observations and confirm the presence of dynamically important magnetic fields around protostars. Additionally, our calculations show that (non-linear) Zeeman effects must be taken into account to further enhance the accuracy of cosmological electron-to-proton mass ratio determinations using methanol.

Context: The importance of magnetic fields at the onset of star formation related to the early fragmentation and collapse processes is largely unexplored today. Aims: We want to understand the magnetic field properties at the earliest evolutionary stages of high-mass star formation. Methods: The Atacama Large Millimeter Array is used at 1.3mm wavelength in full polarization mode to study the polarized emission and by that the magnetic field morphologies and strengths of the high-mass starless region IRDC18310-4. Results: The polarized emission is clearly detected in four sub-cores of the region. In general it shows a smooth distribution, also along elongated cores. Estimating the magnetic field strength via the Davis-Chandrasekhar-Fermi method and following a structure function analysis, we find comparably large magnetic field strengths between ~0.6 and 3.7mG. Comparing the data to spectral line observations, the turbulent-to-magnetic energy ratio is low, indicating that turbulence does not significantly contribute to the stability of the gas clump. A mass-to-flux ratio around the critical value 1.0 - depending on column density - indicates that the region starts to collapse which is consistent with the previous spectral line analysis of the region. Conclusions: While this high-mass region is collapsing and thus at the verge of star formation, the high magnetic field values and the smooth spatial structure indicate that the magnetic field is important for the fragmentation and collapse process. This single case study can only be the starting point for larger sample studies of magnetic fields at the onset of star formation.

The majority of searches for radio emission from exoplanets have to date focused on short period planets, i.e., the so-called hot Jupiter type planets. However, these planets are likely to be tidally locked to their host stars and may not generate sufficiently strong magnetic fields to emit electron cyclotron maser emission at the low frequencies used in observations (typically >150 MHz). In comparison, the large mass-loss rates of evolved stars could enable exoplanets at larger orbital distances to emit detectable radio emission. Here, we first show that the large ionized mass-loss rates of certain evolved stars relative to the solar value could make them detectable with the Low Frequency Array (LOFAR) at 150 MHz ($\lambda$ = 2 m), provided they have surface magnetic field strengths >50 G. We then report radio observations of three long period (>1 au) planets that orbit the evolved stars $\beta$ Gem, $\iota$ Dra, and $\beta$ UMi using LOFAR at 150 MHz. We do not detect radio emission from any system but place tight 3$\sigma$ upper limits of 0.98, 0.87, and 0.57 mJy on the flux density at 150 MHz for $\beta$ Gem, $\iota$ Dra, and $\beta$ UMi, respectively. Despite our non-detections these stringent upper limits highlight the potential of LOFAR as a tool to search for exoplanetary radio emission at meter wavelengths.

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.

Our current understanding of the chemistry and mass-loss processes in solar-like stars at the end of their evolution depends critically on the description of convection, pulsations and shocks in the extended stellar atmosphere. Three-dimensional hydrodynamical stellar atmosphere models provide observational predictions, but so far the resolution to constrain the complex temperature and velocity structures seen in the models has been lacking. Here we present submillimeter continuum and line observations that resolve the atmosphere of the asymptotic giant branch star W Hya. We show that hot gas with chromospheric characteristics exists around the star. Its filling factor is shown to be small. The existence of such gas requires shocks with a cooling time larger than commonly assumed. A shocked hot layer will be an important ingredient in the models of stellar convection, pulsation and chemistry that underlie our current understanding of the late stages of stellar evolution.

The S-type asymptotic giant branch (AGB) star $\pi^{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

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.

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.

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.

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 ($\lambda_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.

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.

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.

We have performed a comprehensive study of the UV emission detected from AGB stars by the Galaxy Evolution Explorer (GALEX). Of the 468 AGB stars in our sample, 316 were observed by GALEX. In the NUV bandpass ($\lambda_{\rm eff} \sim 2310$z A), 179 AGB stars were detected and 137 were not detected. Only 38 AGB stars were detected in the FUV bandpass ($\lambda_{\rm eff} \sim1528$ A). We find that NUV emission is correlated with optical to near infrared emission leading to higher detection fractions among the brightest, hence closest, AGB stars. Comparing the AGB time-variable visible phased light curves to corresponding GALEX NUV phased light curves we find evidence that for some AGB stars the NUV emission varies in phase with the visible light curves. We also find evidence that the NUV emission, and possibly, the FUV emission are anti-correlated with the circumstellar envelope density. These results suggest that the origin of the GALEX-detected UV emission is an inherent characteristic of the AGB stars that can most likely be traced to a combination of photospheric and chromospheric emission. In most cases, UV detections of AGB stars are not likely to be indicative of the presence of binary companions.

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.

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.

Extragalactic observations of water emission can provide valuable insights into the excitation of the interstellar medium. In addition, extragalactic megamasers are powerful probes of kinematics close to active nuclei. Therefore, it is paramount to determine the true origin of the water emission, whether it is excited by processes close to an AGN or in star-forming regions. We use ALMA Band 5 science verification observations to analyse the emission of the 183 GHz water line in Arp 220 on sub-arcsecond scales, in conjunction with new ALMA Band 7 data at 325 GHz. Specifically, the nature of the process leading to the excitation of emission at these water lines is studied in this context. Supplementary 22 GHz VLA observations are used to better constrain the parameter space in the excitation modelling of the water lines. We detect 183 GHz H2O and 325 GHz water emission towards the two compact nuclei at the center of Arp 220, being brighter in Arp 220 West. The emission at these two frequencies is compared to previous single-dish data and does not show evidence of variability. The 183 and 325 GHz lines show similar spectra and kinematics, but the 22 GHz profile is significantly different in both nuclei due to a blend with an NH3 absorption line. Our findings suggest that the most likely scenario to cause the observed water emission in Arp 220 is a large number of independent masers originating from numerous star-forming regions.

(Abridged) The formation of large-scale (hundreds to few thousands of AU) bipolar structures in the circumstellar envelopes (CSEs) of post-Asymptotic Giant Branch (post-AGB) stars is poorly understood. The shape of these structures, traced by emission from fast molecular outflows, suggests that the dynamics at the innermost regions of these CSEs does not depend only on the energy of the radiation field of the central star. Deep into the Water Fountains is an observational project based on the results of programs carried out with three telescope facilities: The Karl G. Jansky Very Large Array (JVLA), The Australia Telescope Compact Array (ATCA), and the Very Large Telescope (SINFONI-VLT). Here we report the results of the observations towards the WF nebula IRAS 18043$-$2116: Detection of radio continuum emission in the frequency range 1.5GHz - 8.0GHz; H$_{2}$O maser spectral features and radio continuum emission detected at 22GHz, and H$_{2}$ ro-vibrational emission lines detected at the near infrared. The high-velocity H$_{2}$O maser spectral features, and the shock-excited H$_{2}$ emission detected could be produced in molecular layers which are swept up as a consequence of the propagation of a jet-driven wind. Using the derived H$_{2}$ column density, we estimated a molecular mass-loss rate of the order of $10^{-9}$M$_{\odot}$yr$^{-1}$. On the other hand, if the radio continuum flux detected is generated as a consequence of the propagation of a thermal radio jet, the mass-loss rate associated to the outflowing ionized material is of the order of 10$^{-5}$M$_{\odot}$yr$^{-1}$. The presence of a rotating disk could be a plausible explanation for the mass-loss rates estimated.

Class I methanol masers are found to be good tracers of the interaction between outflows from massive young stellar objects with their surrounding media. Although polarization observations of Class II methanol masers have been able to provide information about magnetic fields close to the central (proto)stars, polarization observations of Class I methanol masers are rare, especially at 44 and 95 GHz. We present the results of linear polarization observations of 39 Class I methanol maser sources at 44 and 95 GHz. These two lines are observed simultaneously with one of the 21 m Korean VLBI Network telescopes in single-dish mode. Approximately 60% of the observed sources have fractional polarizations of a few percent in at least one transition. This is the first reported detection of linear polarization of the 44 GHz methanol maser. The two maser transitions show similar polarization properties, indicating that they trace similar magnetic environments, although the fraction of the linear polarization is slightly higher at 95 GHz. We discuss the association between the directions of polarization angles and outflows. We also discuss some targets having different polarization properties at both lines, including DR21(OH) and G82.58+0.20, which show the 90{\deg} polarization angle flip at 44 GHz.

The detection and study of molecular gas in born-again stars would be of great importance to understand their composition and chemical evolution. In addition, the molecular emission would be an invaluable tool to explore the physical conditions, kinematics and formation of asymmetric structures in the circumstellar envelopes of these evolved stars. However, until now, all attempts to detect molecular emission from the cool material around born-again stars have failed. We carried out observations using the APEX and IRAM 30m telescopes to search for molecular emission toward four well studied born-again stars, V4334 Sgr, V605 Aql, A30 and A78, that are thought to represent an evolutionary sequence. We detected for the first time emission from HCN and H$^{13}$CN molecules toward V4334 Sgr, and CO emission in V605 Aql. No molecular emission was detected above the noise level toward A30 and A78. A first estimate of the H$^{12}$CN/H$^{13}$CN abundance ratio in the circumstellar environment of V4334 Sgr is $\approx$3, which is similar to the value of the $^{12}$C/$^{13}$C ratio measured from other observations. We derived a rotational temperature of $T_{\rm rot}$=13$\pm1$ K, and a total column density of $N_{{\rm HCN}}$=1.6$\pm0.1\times$10$^{16}$ cm$^{-2}$ for V4334 Sgr. This result sets a lower limit on the amount of hydrogen that was ejected into the wind during the born-again event of this source. For V605 Aql, we obtained a lower limit for the integrated line intensities $I_{^{12}\rm C}$/$I_{^{13}\rm C}$>4.

Abridged. The 12CO/13CO ratio in the circumstellar envelope (CSE) of asymptotic giant branch (AGB) stars has been extensively used as the tracer of the photospheric 12C/13C ratio. However, spatially-resolved ALMA observations of R Scl, a carbon rich AGB star, have shown that the 12CO/13CO ratio is not consistent over the entire CSE. Hence, it can not necessarily be used as a tracer of the 12C/13C ratio. The most likely hypothesis to explain the observed discrepancy between the 12CO/13CO and 12C/13C ratios is CO isotopologue selective photodissociation by UV radiation. Unlike the CO isotopologue ratio, the HCN isotopologue ratio is not affected by UV radiation. Therefore, HCN isotopologue ratios can be used as the tracer of the atomic C ratio in UV irradiated regions. We have performed a detailed non-LTE excitation analysis of circumstellar H12CN and H13CN line emission around R Scl, observed with ALMA and APEX, using a radiative transfer code, ALI. The spatial extent of the molecular distribution for both isotopologues is constrained based on the spatially resolved H13CN(4-3) ALMA observations. We find fractional abundances of H12CN/H2 = (5.0 +\- 2.0) x 10^{-5} and H13CN/H2 = (1.9 +\- 0.4) x 10^{-6} in the inner wind (r < (2.0 +\- 0.25) x 10^{15} cm) of R Scl. The derived circumstellar isotopologue ratio of H12CN/H13CN = 26.3 +\- 11.9 is consistent with the photospheric ratio of 12C/13C ~ 19 \pm 6. We show that the circumstellar H12CN/H13CN ratio traces the photospheric 12C/13C ratio. These results support the previously proposed explanation that CO isotopologue selective-shielding is the main factor responsible for the observed discrepancy between 12C/13C and 12CO/13CO ratios in the inner CSE of R Scl. This indicates that UV radiation impacts on the CO isotopologue ratio.

We have mapped 12CO J=3-2 and other molecular lines from the "water-fountain" bipolar pre-planetary nebula (PPN) IRAS 16342-3814 with ~0."35 resolution using ALMA. We find (i) two very high-speed knotty, jet-like molecular outflows, (ii) a central high-density (> few x 10^6 cm^{-3}), expanding torus of diameter 1300 AU, and (iii) the circumstellar envelope of the progenitor AGB, generated by a sudden, very large increase in the mass-loss rate to >3.5 x 10^{-4} Msun/yr in the past ~455 yr. Strong continuum emission at 0.89 mm from a central source (690 mJy), if due to thermally-emitting dust, implies a substantial mass (0.017 Msun) of very large (~mm-sized) grains. The measured expansion ages of the above structural components imply that the torus (age~160 yr) and the younger high-velocity outflow (age~110 yr) were formed soon after the sharp increase in the AGB mass-loss rate. Assuming a binary model for the jets in IRAS 16342, the high momentum rate for the dominant jet-outflow in IRAS 16342 implies a high minimum accretion rate, ruling out standard Bondi-Hoyle-Lyttleton wind accretion and wind Roche lobe overflow (RLOF) models with white-dwarf or main-sequence companions. Most likely, enhanced RLOF from the primary or accretion modes operating within common envelope evolution are needed.

Yellow hypergiants are rare and represent a fast evolutionary stage of massive evolved stars. That evolutionary phase is characterised by a very intense mass loss, the understanding of which is still very limited. Here we report ALMA Compact Array observations of a 50$"$-mosaic toward the Fried Egg nebula, around one of the few Galactic yellow hypergiants IRAS 17163-3907. The emission from the $^{12}$CO J=2-1 line, H30$\alpha$ recombination line, and continuum is imaged at a resolution of $\sim$8$"$, revealing the morphology of the molecular environment around the star. The continuum emission is unresolved and peaks at the position of the star. The radio recombination line H30$\alpha$ shows unresolved emission at the star, with an approximately gaussian spectrum centered on a velocity of 21$\pm$3~km/s with a width of 57$\pm$6~km/s. In contrast, the CO 2-1 emission is complex and decomposes into several components beyond the contamination from interstellar gas in the line of sight. The CO spectrum toward the star is a broad plateau, centered at the systemic velocity of +18 km/s and with an expansion velocity of 100$\pm$10 km/s. Assuming isotropic and constant mass-loss, we estimate a mass-loss rate of 8$\pm$1.5 $\times10^{-5}$~M$_\odot$ yr$^{-1}$. At a radius of 25$"$ from the star, we detect CO emission associated with the dust ring previously imaged by {\it Herschel}. The kinematics of this ring, however, is not consistent with an expanding shell, but show a velocity gradient of $v_{sys} \pm$20 km/s. In addition, we find a puzzling bright feature radially connecting the star to the CO ring, at a velocity of +40 km/s relative to the star. This spur feature may trace a unidirectional ejection event from the star. Our ACA observations reveal the complex morphology around IRAS 17163 and illustrate the breakthroughs that ALMA will bring to the field of massive stellar evolution.

We report the detection of thermal continuum radio emission from the K0 III coronal giant Pollux ($\beta$ Gem) with the Karl G. Jansky Very Large Array (VLA). The star was detected at 21 and 9 GHz with flux density values of $150\pm21$ and $43\pm8\,\mu$Jy, respectively. We also place a $3\sigma_{\mathrm{rms}}$ upper limit of $23\,\mu$Jy for the flux density at 3 GHz. We find the stellar disk-averaged brightness temperatures to be approximately 9500, 15000, and $<71000\,$K, at 21, 9, and 3 GHz, respectively, which are consistent with the values of the quiet Sun. The emission is most likely dominated by optically thick thermal emission from an upper chromosphere at 21 and 9 GHz. We discuss other possible additional sources of emission at all frequencies and show that there may also be a small contribution from gyroresonance emission above active regions, coronal free-free emission and free-free emission from an optically thin stellar wind, particularly at the lower frequencies. We constrain the maximum mass-loss rate from Pollux to be less than $3.7\times 10^{-11}\,M_{\odot}$ yr$^{-1}$ (assuming a wind terminal velocity of 215 km s$^{-1}$), which is about an order of magnitude smaller than previous constraints for coronal giants and is in agreement with existing predictions for the mass-loss rate of Pollux. These are the first detections of thermal radio emission from a single (i.e., non-binary) coronal giant and demonstrate that low activity coronal giants like Pollux have atmospheres at radio frequencies akin to the quiet Sun.

Aluminium monoxide, AlO, is likely efficiently depleted from the gas around oxygen-rich evolved stars to form alumina clusters and dust seeds. Its presence in the extended atmospheres of evolved stars has been derived from optical spectroscopy. More recently, AlO gas was also detected at long wavelengths around the supergiant VY CMa and the oxygen-rich asymptotic giant branch (AGB) star o Cet (Mira A). In search of AlO, we mined data obtained with APEX, the IRAM 30m telescope, Herschel/HIFI, SMA, and ALMA, which were primarily aimed at studying other molecular species. We report here on observations of AlO towards a sample of eight oxygen-rich AGB stars in different rotational transitions, up to seven for some stars. We present definite detections of one rotational transition of AlO for o Cet and R Aqr, and tentative detections of one transition for R Dor and o Cet, and two for IK Tau and W Hya. The presented spectra of WX Psc, R Cas, and TX Cam show no signature of AlO. For o Cet, R Aqr, and IK Tau, we find that the AlO(N=9-8) emission likely traces the inner parts of the wind, out to only a few tens of AU, where the gas has not yet reached its terminal velocity. The conclusive detections of AlO emission in the case of o Cet and R Aqr confirm the presence of AlO gas in outflows of AGB stars. The tentative detections further support this. Since most of the observations presented in this study were obtained with stronger emission from other species than AlO in mind, observations with higher sensitivity in combination with high angular resolution will improve our understanding of the presence and behaviour of AlO. From the current data sets we cannot firmly conclude whether there is a direct correlation between the wind properties and the detection rate of AlO emission. We hope that this study can serve as a stimulus to perform sample studies in search of AlO in oxygen-rich outflows.

Context: Low- and intermediate-mass stars lose most of their stellar mass at the end of their lives on the asymptotic giant branch (AGB). Determining gas and dust mass-loss rates (MLRs) is important in quantifying the contribution of evolved stars to the enrichment of the interstellar medium. Aims: Attempt to, for the first time, spectrally resolve CO thermal line emission in a small sample of AGB stars in the Large Magellanic Cloud. Methods: ALMA was used to observe 2 OH/IR stars and 4 carbon stars in the LMC in the CO J= 2-1 line. Results: We present the first measurement of expansion velocities in extragalactic carbon stars. All four C-stars are detected and wind expansion velocities and stellar velocities are directly measured. Mass-loss rates are derived from modelling the spectral energy distribution and Spitzer/IRS spectrum with the DUSTY code. Gas-to-dust ratios are derived that make the predicted velocities agree with the observed ones. The expansion velocities and MLRs are compared to a Galactic sample of well-studied relatively low MLRs stars supplemented with "extreme" C-stars that have properties more similar to the LMC targets. Gas MLRs derived from a simple formula are significantly smaller than derived from the dust modelling, indicating an order of magnitude underestimate of the estimated CO abundance, time-variable mass loss, or that the CO intensities in LMC stars are lower than predicted by the formula derived for Galactic objects. This could be related to a stronger interstellar radiation field in the LMC. Conclusions: Although the LMC sample is small and the comparison to Galactic stars is non-trivial because of uncertainties in their distances it appears that for C stars the wind expansion velocities in the LMC are lower than in the solar neighbourhood, while the MLRs appear similar. This is in agreement with dynamical dust-driven wind models.

We report the serendipitous detection with ALMA of the vibrationally-excited pure-rotational CO transition $v=1, J=3-2$ towards five asymptotic giant branch (AGB) stars, $o$ Cet, R Aqr, R Scl, W Aql, and $\pi^1$ Gru. The observed lines are formed in the poorly-understood region located between the stellar surface and the region where the wind starts, the so-called warm molecular layer. We successfully reproduce the observed lines profiles using a simple model. We constrain the extents, densities, and kinematics of the region where the lines are produced. R Aqr and R Scl show inverse P-Cygni line profiles which indicate infall of material onto the stars. The line profiles of $o$ Cet and R Scl show variability. The serendipitous detection towards these five sources shows that vibrationally-excited rotational lines can be observed towards a large number of nearby AGB stars using ALMA. This opens a new possibility for the study of the innermost regions of AGB circumstellar envelopes.

We report full polarimetric VLBA observations of water masers towards the Turner-Welch Object in the W3(OH) high-mass star forming complex. This object drives a synchrotron jet, which is quite exceptional for a high-mass protostar, and is associated with a strongly polarized water maser source, W3(H$_2$O), making it an optimal target to investigate the role of magnetic fields on the innermost scales of protostellar disk-jet systems. The linearly polarized emission from water masers provides clues on the orientation of the local magnetic field, while the measurement of the Zeeman splitting from circular polarization provides its strength. The water masers trace a bipolar, biconical outflow at the center of the synchrotron jet. Although on scales of a few thousand AU the magnetic field inferred from the masers is on average orientated along the flow axis, on smaller scales (10s to 100s of AU), we have revealed a misalignment between the magnetic field and the velocity vectors, which arises from the compression of the field component along the shock front. Our measurements support a scenario where the magnetic field would evolve from having a dominant component parallel to the outflow velocity in the pre-shock gas, with field strengths of the order of a few tens of mG (at densities of $10^7$ cm$^{-3}$), to being mainly dominated by the perpendicular component of order of a few hundred of mG in the post-shock gas where the water masers are excited (at densities of $10^9$ cm$^{-3}$). The general implication is that in the undisturbed (i.e. not-shocked) circumstellar gas, the flow velocities would follow closely the magnetic field lines, while in the gas shocked by the prostostellar jet the magnetic field would be re-configured to be parallel to the shock front.

Current mm/submm interferometers, like the Atacama Large mm/submm Array (ALMA), use receivers that register the sky signal in a linear polarization basis. In the case of observations performed in full-polarization mode (where the cross-correlations are computed among all the polarization channels) it is possible to reconstruct the full-polarization brightness distribution of the observed sources, as long as a proper calibration of delay offsets and leakage among polarization channels can be performed. Observations of calibrators, preferably with some linear polarization, with a good parallactic angle coverage are usually needed for such a calibration. In principle, dual-polarization observations only allow us to recover the Stokes $I$ intensity distribution of the sources, regardless of the parallactic angle coverage of the observations. In this paper, we present a novel technique of dual differential polarimetry that makes it possible to obtain information related to the full-polarization brightness distribution of the observed sources from dual-polarization observations. This technique is inspired in the Earth-rotation polarization synthesis and can be applied even to sources with complex structures.

Aim: The aim of this work is to search Seyfert 2 galaxy NGC 4945, a well-known 22 GHz water megamaser galaxy, for water (mega)maser emission at 183 GHz. Method: We used APEX SEPIA Band 5 to perform the observations. Results: We detected 183 GHz water maser emission towards NGC 4945 with a peak flux density of ~3 Jy near the galactic systemic velocity. The emission spans a velocity range of several hundred km/s. We estimate an isotropic luminosity of > 1000 Lsun, classifying the emission as a megamaser. A comparison of the 183 GHz spectrum with that observed at 22 GHz suggests that 183 GHz emission also arises from the active galactic nucleus (AGN) central engine. If the 183 GHz emission originates from the circumnuclear disk, then we estimate that a redshifted feature at 1084 km/s in the spectrum should arise from a distance of 0.022 pc from the supermassive black hole (1.6 x 10(5) Schwarzschild radii), i.e. closer than the water maser emission previously detected at 22 GHz. This is only the second time 183 GHz maser emission has been detected towards an AGN central engine (the other galaxy being NGC 3079). It is also the strongest extragalactic millimetre/submillimetre water maser detected to date. Conclusions: Strong millimetre 183 GHz water maser emission has now been shown to occur in an external galaxy. For NGC 4945, we believe that the maser emission arises, or is dominated by, emission from the AGN central engine. Emission at higher velocity, i.e. for a Keplerian disk closer to the black hole, has been detected at 183 GHz compared with that for the 22 GHz megamaser. This indicates that millimetre/submillimetre water masers can indeed be useful probes for tracing out more of AGN central engine structures and dynamics than previously probed. Future observations using ALMA Band 5 should unequivocally determine the origin of the emission in this and other galaxies.

We use high-angular-resolution images obtained with SPHERE/ZIMPOL to study the photosphere, the warm molecular layer, and the inner wind of the close-by oxygen-rich AGB star R Doradus. We present observations in filters V, cntH$\alpha$, and cnt820 and investigate the surface brightness distribution of the star and of the polarised light produced in the inner envelope. Thanks to second-epoch observations in cntH$\alpha$, we are able to see variability on the stellar photosphere. We find that in the first epoch the surface brightness of R Dor is asymmetric in V and cntH$\alpha$, the filters where molecular opacity is stronger, while in cnt820 the surface brightness is closer to being axisymmetric. The second-epoch observations in cntH$\alpha$ show that the morphology of R Dor changes completely in a timespan of 48 days to a more axisymmetric and compact configuration. The polarised intensity is asymmetric in all epochs and varies by between a factor of 2.3 and 3.7 with azimuth for the different images. We fit the radial profile of the polarised intensity using a spherically symmetric model and a parametric description of the dust density profile, $\rho(r)=\rho_\circ r^{-n}$. On average, we find exponents of $- 4.5 \pm 0.5$ that correspond to a much steeper density profile than that of a wind expanding at constant velocity. The dust densities we derive imply an upper limit for the dust-to-gas ratio of $\sim 2\times10^{-4}$ at 5.0 $R_\star$. Given the uncertainties in observations and models, this value is consistent with the minimum values required by wind-driving models for the onset of a wind, of $\sim 3.3\times10^{-4}$. However, if the steep density profile we find extends to larger distances from the star, the dust-to-gas ratio will quickly become too small for the wind of R Dor to be driven by the grains that produce the scattered light.

The observed relation between the X-ray and radio properties of low-luminosity accreting black holes has enabled the identification of multiple candidate black hole X-ray binaries (BHXBs) in globular clusters. Here we report an identification of the radio source VLA J213002.08+120904 (aka M15 S2), recently reported in Kirsten et al. 2014, as a BHXB candidate. They showed that the parallax of this flat-spectrum variable radio source indicates a 2.2$^{+0.5}_{-0.3}$ kpc distance, which identifies it as lying in the foreground of the globular cluster M15. We determine the radio characteristics of this source, and place a deep limit on the X-ray luminosity of $\sim4\times10^{29}$ erg s$^{-1}$. Furthermore, we astrometrically identify a faint red stellar counterpart in archival Hubble images, with colors consistent with a foreground star; at 2.2 kpc its inferred mass is 0.1-0.2 $M_{\odot}$. We rule out that this object is a pulsar, neutron star X-ray binary, cataclysmic variable, or planetary nebula, concluding that VLA J213002.08+120904 is the first accreting black hole X-ray binary candidate discovered in quiescence outside a globular cluster. Given the relatively small area over which parallax studies of radio sources have been performed, this discovery suggests a much larger population of quiescent BHXBs in our Galaxy, $2.6\times10^4-1.7\times10^8$ BHXBs at $3\sigma$ confidence, than has been previously estimated ($\sim10^2-10^4$) through population synthesis.

Active galactic nuclei (AGN) host some of the most energetic phenomena in the Universe. AGN are thought to be powered by accretion of matter onto a rotating disk that surrounds a supermassive black hole. Jet streams can be boosted in energy near the event horizon of the black hole and then flow outward along the rotation axis of the disk. The mechanism that forms such a jet and guides it over scales from a few light-days up to millions of light-years remains uncertain, but magnetic fields are thought to play a critical role. Using the Atacama large mm/submm array (ALMA), we have detected a polarization signal (Faraday rotation) related to the strong magnetic field at the jet base of a distant AGN, PKS1830-211. The amount of Faraday rotation (rotation measure) is proportional to the magnetic field strength along the line of sight times the density of electrons. Although it is impossible to precisely infer the magnetic fields in the region of Faraday rotation, the high rotation measures derived suggest magnetic fields of at least tens of Gauss (and possibly considerably higher) on scales of the order of light days (0.01 pc) from the black hole.

The unprecedented sensitivity of the Atacama Large millimeter/submillimeter array (ALMA) is providing many new discoveries. Several of these are serendipitous to the original goal of the observations. We report the discovery of previously unknown continuum sources, or a single fast moving new source, in our ALMA observations. Here we aim to determine the nature of the detections. The detections, at $>5.8\sigma$ in the image plane and $>14\sigma$ in the $(u,v)-$plane, were made in two epochs of ALMA observations of a $25$ arc second region around the asymptotic giant branch star W Aql in the continuum around 345 GHz. At a third epoch, covering $50x50$ arcseconds, the source(s) were not seen. We have investigated if the detections could be spurious, if they could constitute a population of variable background sources, or if the observations revealed a fast moving single object. Based on our analysis, we conclude that a single object (with a flux of $\sim3.0$ mJy) exhibiting a large proper motion ($\sim87$ arcsec/yr) is the most likely explanation. Until the nature of the source becomes clear, we have named it Gna. Unless there are yet unknown, but significant, issues with ALMA observations, we have detected a previously unknown objects in our solar system. Based on proper motion analysis we find that, if it is gravitationally bound, Gna is currently located at $12-25$ AU distance and has a size of $\sim220-880$ km. Alternatively it is a much larger, planet-sized, object, gravitationally unbound, and located within $\sim4000$ AU, or beyond (out to $\sim0.3$~pc) if it is strongly variable. Our observations highlight the power of ALMA in detecting possible solar system objects, but also show how multiple epoch observations are crucial to identify what are otherwise probably assumed to be extra-galactic sources.

In recent years there has been a paradigm shift from centralised to geographically distributed resources. Individual entities are no longer able to host or afford the necessary expertise in-house, and, as a consequence, society increasingly relies on widespread collaborations. Although such collaborations are now the norm for scientific projects, more technical structures providing support to a distributed scientific community without direct financial or other material benefits are scarce. The network of European ALMA Regional Centre (ARC) nodes is an example of such an internationally distributed user support network. It is an organised effort to provide the European ALMA user community with uniform expert support to enable optimal usage and scientific output of the ALMA facility. The network model for the European ARC nodes is described in terms of its organisation, communication strategies and user support.

We recently announced the detection of an unknown submillimeter source in our ALMA observations of alpha Cen AB. The source was detected in two epochs, a strong detection at 445~GHz and one at lower significance at 343.5~GHz. After valuable feedback of the community, it turns out that the detection at 343.5~GHz could not be reproduced with a different reduction software nor with fitting within the $(u,v)$-plane. The detection at 445~GHz has been further confirmed with modeling of the $(u,v)$-data and was shown to be robust at $>12\sigma$, confirming our detection of this unknown source. However, based on only one epoch, further analysis and preferably new data are needed, before publication of an article in which the nature of the new source can be discussed. The analysis has indicated the importance of both $(u,v)$-plane fitting and alternative data reduction when dealing with low signal to noise source detections.

Thermal pulses are fundamental to the chemical evolution of AGB stars and their circumstellar envelopes. A further consequence of thermal pulses is the formation of detached shells of gas and dust around the star. We aim to determine the physical properties of the detached gas shell around R Sculptoris, in particular the shell mass and temperature, and to constrain the evolution of the mass-loss rate during and after a thermal pulse. We analyse CO(1-0), CO(2-1), and CO(3-2) emission, observed by. The spatial resolution of the ALMA data allows us to separate the detached shell emission from the extended emission inside the shell. We perform radiative transfer modelling of both components to determine the shell properties and the post-pulse mass-loss properties. The ALMA data show a gas shell with a radius of 19.5" expanding at 14.3km/s. The different scales probed by the ALMA Cycle 0 array show that the shell must be entirely filled with gas, contrary to the idea of a detached shell. The comparison to single-dish spectra and radiative transfer modelling confirms this. We derive a shell mass of 4.5e-3 Msun with a temperature of 50K. Typical timescales for thermal pulses imply a pulse mass-loss rate of 2.3e-5 Msun/yr. For the post-pulse mass-loss rate, we find evidence for a gradual decline of the mass-loss rate, with an average value of 1.6e-5 Msun/yr. The total amount of mass lost since the last thermal pulse is 0.03 Msun, a factor four higher compared to classical models, with a sharp decline in mass-loss rate immediately after the pulse. We find that the mass-loss rate after a thermal pulse has to decline more slowly than generally expected from models of thermal pulses. This may cause the star to lose significantly more mass during a thermal pulse cycle, which affects the chemical evolution of the star and the interstellar medium.

We want to study the velocity and magnetic field morphology in the vicinity (<1000 AU) of a massive young stellar object (YSO), at very high spatial resolution (10-100 AU). We performed milli-arcsecond polarimetric observations of the strong CH3OH maser emission observed in the vicinity of an O-type YSO, in G023.01-00.41. We have combined this information with the velocity field of the CH3OH masing gas previously measured at the same angular resolution. We analyse the velocity and magnetic fields in the reference system defined by the direction of the molecular outflow and the equatorial plane of the hot molecular core at its base, as recently observed on sub-arcsecond scales. We provide a first detailed picture of the gas dynamics and magnetic field configuration within a radius of 2000 AU from a massive YSO. We have been able to reproduce the magnetic field lines for the outer regions (>600 AU) of the molecular envelope, where the magnetic field orientation shows a smooth change with the maser cloudlets position (0.2 degree/AU). Overall, the velocity field vectors well accommodate with the local, magnetic field direction, but still show an average misalignment of 30 degrees. We interpret this finding as the contribution of a turbulent velocity field of about 3.5 km/s, responsible for braking up the alignment between the velocity and magnetic field vectors. We do resolve different gas flows which develop both along the outflow axis and across the disk plane, with an average speed of 7 km/s. In the direction of the outflow axis, we establish a collimation of the gas flow, at a distance of about 1000 AU from the disk plane. In the disk region, gas appears to stream outward along the disk plane for radii greater than 500-600 AU, and inward for shorter radii.

The current paradigm of star formation through accretion disks, and magnetohydrodynamically driven gas ejections, predicts the development of collimated outflows, rather than expansion without any preferential direction. We present radio continuum observations of the massive protostar W75N(B)-VLA 2, showing that it is a thermal, collimated ionized wind and that it has evolved in 18 years from a compact source into an elongated one. This is consistent with the evolution of the associated expanding water-vapor maser shell, which changed from a nearly circular morphology, tracing an almost isotropic outflow, to an elliptical one outlining collimated motions. We model this behavior in terms of an episodic, short-lived, originally isotropic, ionized wind whose morphology evolves as it moves within a toroidal density stratification.

We use the Very Large Array (VLA) in the A configuration with the Pie Town (PT) Very Long Baseline Array (VLBA) antenna to spatially resolve the extended atmosphere of Betelgeuse over multiple epochs at 0.7, 1.3, 2.0, 3.5, and 6.1 cm. The extended atmosphere deviates from circular symmetry at all wavelengths while at some epochs we find possible evidence for small pockets of gas significantly cooler than the mean global temperature. We find no evidence for the recently reported e-MERLIN radio hotspots in any of our multi-epoch VLA/PT data, despite having sufficient spatial resolution and sensitivity at short wavelengths, and conclude that these radio hotspots are most likely interferometric artefacts. The mean gas temperature of the extended atmosphere has a typical value of 3000 K at 2 $R_{\star}$ and decreases to 1800 K at 6 $R_{\star}$, in broad agreement with the findings of the single epoch study from Lim et al. (1998). The overall temperature profile of the extended atmosphere between $2 R_{\star} \lesssim r \lesssim 6 R_{\star}$ can be described by a power law of the form $T_{\mathrm{gas}}(r) \propto r^{-0.6}$, with temporal variability of a few 100 K evident at some epochs. Finally, we present over 12 years of V band photometry, part of which overlaps our multi-epoch radio data. We find a correlation between the fractional flux density variability at V band with most radio wavelengths. This correlation is likely due to shock waves induced by stellar pulsations, which heat the inner atmosphere and ionize the more extended atmosphere through radiative means. Stellar pulsations may play an important role in exciting Betelgeuse's extended atmosphere.

Titanium dioxide, TiO$_2$, is a refractory species that could play a crucial role in the dust-condensation sequence around oxygen-rich evolved stars. To date, gas phase TiO$_2$ has been detected only in the complex environment of the red supergiant VY CMa. We aim to constrain the distribution and excitation of TiO$_2$ around VY CMa in order to clarify its role in dust formation. We analyse spectra and channel maps for TiO$_2$ extracted from ALMA science verification data. We detect 15 transitions of TiO$_2$, and spatially resolve the emission for the first time. The maps demonstrate a highly clumpy, anisotropic outflow in which the TiO$_2$ emission likely traces gas exposed to the stellar radiation field. A roughly east-west oriented, accelerating bipolar-like structure is found, of which the blue component runs into and breaks up around a solid continuum component. A distinct tail to the south-west is seen for some transitions, consistent with features seen in the optical and near-infrared. We find that a significant fraction of TiO$_2$ remains in the gas phase outside the dust-formation zone and suggest that this species might play only a minor role in the dust-condensation process around extreme oxygen-rich evolved stars like VY CMa.

Recent high-resolution observations have shown stellar winds to harbour complexities which strongly deviate from spherical symmetry, generally assumed as standard wind model. One such morphology is the archimedean spiral, generally believed to be formed by binary interactions, which has been directly observed in multiple sources. We seek to investigate the manifestation in the observables of spiral structures embedded in the spherical outflows of cool stars. We aim to provide an intuitive bedrock with which upcoming ALMA data can be compared and interpreted. By means of an extended parameter study, we model rotational CO emission from the stellar outflow of asymptotic giant branch stars. To this end, we develop a simplified analytical parametrised description of a 3D spiral structure. This model is embedded into a spherical wind, and fed into the 3D radiative transfer code LIME, which produces 3D intensity maps throughout velocity space. Subsequently, we investigate the spectral signature of rotational transitions of CO of the models, as well as the spatial aspect of this emission by means of wide-slit PV diagrams. Additionally, the potential for misinterpretation of the 3D data in a 1D context is quantified. Finally, we simulate ALMA observations to explore the impact of interefrometric noise and artifacts on the emission signatures. The spectral signatures of the CO rotational transition v=0 J=3-2 are very efficient at concealing the dual nature of the outflow. Only a select few parameter combinations allow for the spectral lines to disclose the presence of the spiral structure. The inability to disentangle the spiral from the spherical signal can result in an incorrect interpretation in a 1D context. Consequently, erroneous mass loss rates would be calculated...

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

Theoretical simulations and observations at different angular resolutions have shown that magnetic fields have a central role in massive star formation. Like in low-mass star formation, the magnetic field in massive young stellar objects can either be oriented along the outflow axis or randomly. Measuring the magnetic field at milliarcsecond resolution (10-100 au) around a substantial number of massive young stellar objects permits determining with a high statistical significance whether the direction of the magnetic field is correlated with the orientation of the outflow axis or not. In late 2012, we started a large VLBI campaign with the European VLBI Network to measure the linearly and circularly polarized emission of 6.7 GHz methanol masers around a sample of massive star-forming regions. This paper focuses on the first seven observed sources, G24.78+0.08, G25.65+1.05, G29.86-0.04, G35.03+0.35, G37.43+1.51, G174.20-0.08, and G213.70-12.6. For all these sources, molecular outflows have been detected in the past. We detected a total of 176 methanol masing cloudlets toward the seven massive star-forming regions, 19% of which show linearly polarized emission. The methanol masers around the massive young stellar object MM1 in G174.20-0.08 show neither linearly nor circularly polarized emission. The linear polarization vectors are well ordered in all the other massive young stellar objects. We measured significant Zeeman splitting toward both A1 and A2 in G24.78+0.08, and toward G29.86-0.04 and G213.70-12.6. By considering all the 19 massive young stellar objects reported in the literature for which both the orientation of the magnetic field at milliarcsecond resolution and the orientation of outflow axes are known, we find evidence that the magnetic field (on scales 10-100 au) is preferentially oriented along the outflow axes.

We present the size, shape and flux densities at millimeter continuum wavelengths, based on ALMA science verification observations in Band 3 (~94.6 GHz) and Band 6 (~228.7 GHz), from the binary Mira A (o Ceti) and Mira B. The Mira AB system has been observed with ALMA at a spatial resolution of down to ~25 mas. The extended atmosphere of Mira A and the wind around Mira B sources are resolved and we derive the size of Mira A and of the ionized region around Mira B. The spectral indices within Band 3 (between 89-100 GHz) and between Band 3 and Band 6 are also derived. The spectral index of Mira A is found to change from 1.71+-0.05 within Band 3 to 1.54+-0.04 between Band 3 and 6. The spectral index of Mira B is 1.3+-0.2 in Band 3, in good agreement with measurements at longer wavelengths. However it rises to 1.72+-0.11 between the bands. For the first time the extended atmosphere of a star is resolved at these frequencies and for Mira A the diameter is ~3.8x3.2 AU in Band 3 (with brightness temperature Tb~5300 K) and ~4.0x3.6 AU in Band 6 (Tb~2500 K). Additionally, a bright hotspot of ~0.4 AU and with Tb~10000 K is found on the stellar disc of Mira A. The size of the ionized region around the accretion disk of Mira B is found to be ~2.4 AU. The emission around Mira B is consistent with that from a partially ionized wind of gravitationally bound material from Mira A close to the accretion disk of Mira B. The Mira A atmosphere does not fully match predictions, with brightness temperatures in Band 3 significantly higher than expected, potentially due to shock heating. The hotspot is likely due to magnetic activity and could be related to the previously observed X-ray flare of Mira A.

We discuss the unique opportunities for maser astrometry with the inclusion of the Square Kilometre Array (SKA) in Very Long Baseline Interferometry (VLBI) networks. The first phase of the SKA will enable observations of hydroxyl and methanol masers, positioning the latter to an accuracy of 5 microarcseconds, and the second phase may allow water maser observations. These observations will provide trigonometric distances with errors as small as 1%. The unrivalled sensitivity of the SKA will enable large-scale surveys and, through joint operations, will turn any VLBI network into a fast astrometry device. Both evolved stars and high mass star formation regions will be accessible throughout the (Southern) Milky Way, completing our understanding of the content, dynamics and history of our Galaxy. Maser velocities and proper motions will be measurable in the Local Group of galaxies and beyond, providing new insights into their kinematics and evolution.

We present new proper motion and parallax measurements obtained with the European VLBI Network (EVN) at 5$\,$GHz for the three isolated pulsars B1929+10, B2020+28, and B2021+51. For B1929+10 we combined our data with earlier VLBI measurements and confirm the robustness of the astrometric parameters of this pulsar. For pulsars B2020+28 and B2021+51 our observations indicate that both stars are almost a factor of two closer to the solar system than previously thought, placing them at a distance of $1.39_{-0.06}^{+0.05}$ and $1.25_{-0.17}^{+0.14}\,$kpc. Using our new astrometry, we simulated the orbits of all three pulsars in the Galactic potential with the aim to confirm or reject previously proposed birth locations. Our observations ultimately rule out a claimed binary origin of B1929+10 and the runaway star $\zeta$ Ophiuchi in Upper Scorpius. A putative common binary origin of B2020+28 and B2021+51 in the Cygnus Superbubble is also very unlikely.

We study the circumstellar evolution of the binary HD101584, consisting of a post-AGB star and a low-mass companion, which is most likely a post-common-envelope-evolution system. We used ALMA observations of the 12CO, 13CO, and C18O J=2-1 lines and the 1.3mm continuum to determine the morphology, kinematics, masses, and energetics of the circumstellar environment. The circumstellar medium has a bipolar hour-glass structure, seen almost pole-on, formed by an energetic jet, about 150 km/s. We conjecture that the circumstellar morphology is related to an event that took place about 500 year ago, possibly a capture event where the companion spiraled in towards the AGB star. However, the kinetic energy of the accelerated gas exceeds the released orbital energy, and, taking into account the expected energy transfer efficiency of the process, the observed phenomenon does not match current common-envelope scenarios. This suggests that another process must augment, or even dominate, the ejection process. A significant amount of material resides in an unresolved region, presumably in the equatorial plane of the binary system.

Magnetic fields have only recently been included in theoretical simulations of high-mass star formation. The simulations show that magnetic fields can play a crucial role not only in the formation and dynamics of molecular outflows, but also in the evolution of circumstellar disks. Therefore, new measurements of magnetic fields at milliarcsecond resolution close to massive young stellar objects (YSOs) are fundamental for providing new input for numerical simulations and for understanding the formation process of massive stars. The polarized emission of 6.7 GHz CH3OH masers allows us to investigate the magnetic field close to the massive YSO where the outflows and disks are formed. Recently, we have detected with the EVN CH3OH maser polarized emission towards 10 massive YSOs. From a first statistical analysis we have found evidence that magnetic fields are primarily oriented along the molecular outflows. To improve our statistics we are carrying on a large observational EVN campaign for a total of 19 sources, the preliminary results of the first seven sources are presented in this contribution. Furthermore, we also describe our efforts to estimate the Lande' g-factors of the CH3OH maser transition to determine the magnetic field strength from our Zeeman-splitting measurements.

The measurement of Zeeman splitting in spectral lines---both in emission and absorption---can provide direct estimates of the magnetic field strength and direction in atomic and molecular clouds, both in our own Milky Way and in external galaxies. This method will probe the magnetic field in the warm and cold neutral components of the interstellar medium, providing a complement to the extensive SKA Faraday studies planning to probe the field in the ionized components.

Stacking as a tool for studying objects that are not individually detected is becoming popular even for radio interferometric data, and will be widely used in the SKA era. Stacking is typically done using imaged data rather than directly using the visibilities (the uv-data). We have investigated and developed a novel algorithm to do stacking using the uv-data. We have performed exten- sive simulations comparing to image-stacking, and summarize the results of these simulations. Furthermore, we disuss the implications in light of the vast data volume produced by the SKA. Having access to the uv-stacked data provides a great advantage, as it allows the possibility to properly analyse the result with respect to calibration artifacts as well as source properties such as size. For SKA the main challenge lies in archiving the uv-data. For purposes of robust stacking analysis, it would be strongly desirable to either keep the calibrated uv-data at least in an aver- age form, or implement a stacking queue where stacking positions could be provided prior to the observations and the uv-stacking is done almost in real time.

Adding VLBI capability to the SKA arrays will greatly broaden the science of the SKA, and is feasible within the current specifications. SKA-VLBI can be initially implemented by providing phased-array outputs for SKA1-MID and SKA1-SUR and using these extremely sensitive stations with other radio telescopes, and in SKA2 by realising a distributed configuration providing baselines up to thousands of km, merging it with existing VLBI networks. The motivation for and the possible realization of SKA-VLBI is described in this paper.

The precise mechanisms that provide the non-radiative energy for heating the chromosphere and the corona of the Sun and those of other stars constitute an active field of research. By studying stellar chromospheres one aims at identifying the relevant physical processes. Defining the permittable extent of the parameter space can also serve as a template for the Sun-as-a-star. Earlier observations with Herschel and APEX have revealed the temperature minimum of alpha Cen, but these were unable to spatially resolve the binary into individual components. With the data reported here, we aim at remedying this shortcoming. Furthermore, these earlier data were limited to the wavelength region between 100 and 870mu. In the present context, we intend to extend the spectral mapping to longer wavelengths, where the contrast between stellar photospheric and chromospheric emission becomes increasingly evident. ALMA is particularly suited to point sources, such as unresolved stars. ALMA provides the means to achieve our objectives with both its high sensitivity of the collecting area for the detection of weak signals and the high spatial resolving power of its adaptable interferometer for imaging close multiple stars. This is the first detection of main-sequence stars at a wavelength of 3mm. Furthermore, the individual components of the binary alpha CenAB are clearly detected and spatially well resolved at all ALMA wavelengths. The high S/N of these data permit accurate determination of their relative flux ratios. The previously obtained flux ratio of 0.44, which was based on measurements in the optical and at 70mu, is consistent with the present ALMA results, albeit with a large error bar. Given the distinct difference in their cyclic activity, the similarity of their submm SEDs appears surprising.

The processes leading to dust formation and the subsequent role it plays in driving mass loss in cool evolved stars is an area of intense study. Here we present high resolution ALMA Science Verification data of the continuum emission around the highly evolved oxygen-rich red supergiant VY CMa. These data enable us to study the dust in its inner circumstellar environment at a spatial resolution of 129 mas at 321 GHz and 59 mas at 658 GHz, thus allowing us to trace dust on spatial scales down to 11 R$_{\star}$ (71 AU). Two prominent dust components are detected and resolved. The brightest dust component, C, is located 334 mas (61 R$_{\star}$) South East of the star and has a dust mass of at least $2.5\times 10^{-4}$ M$_{\odot}$. It has a dust emissivity spectral index of $\beta =-0.1$ at its peak, implying that it is optically thick at these frequencies with a cool core of $T_{d}\lesssim 100$ K. Interestingly, not a single molecule in the ALMA data has emission close to the peak of this massive dust clump. The other main dust component, VY, is located at the position of the star and contains a total dust mass of $4.0 \times 10^{-5} $M$_{\odot}$. It also contains a weaker dust feature extending over $60$ R$_{\star}$ to the North with the total component having a typical dust emissivity spectral index of $\beta =0.7$. We find that at least $17\%$ of the dust mass around VY CMa is located in clumps ejected within a more quiescent roughly spherical stellar wind, with a quiescent dust mass loss rate of $5 \times 10^{-6}$ M$_{\odot} $yr$^{-1}$. The anisotropic morphology of the dust indicates a continuous, directed mass loss over a few decades, suggesting that this mass loss cannot be driven by large convection cells alone.

We present a new algorithm for stacking radio interferometric data in the uv-domain. The performance of uv-stacking is compared to the stacking of fully imaged data using simulated Atacama Large Millimeter/sub-millimeter Array (ALMA) and the Karl G. Jansky Very Large Array (VLA) deep extragalactic surveys. We find that image- and uv-stacking produce similar results, however, uv-stacking is typically the more robust method. An advantage of the uv-stacking algorithm is the availability of uv-data post stacking, which makes it possible to identify and remove problematic baselines. For deep VLA surveys uv-stacking yields a signal-to-noise ratio that is up to 20 per cent higher than image-stacking. Furthermore, we have investigated stacking of resolved sources with a simulated VLA data set where 1.5" (10-12 kpc at z ~ 1-4) sources are stacked. We find that uv-stacking, where a model is fitted directly to the visibilities, significantly improves the accuracy and robustness of the size estimates. While scientific motivation for this work is studying faint, high-z galaxies, the algorithm analysed here would also be applicable in other fields of astronomy. Stacking of radio interferometric data is also expected to play a big role for future surveys with telescopes such as LOFAR and Square Kilometre Array (SKA).

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.

We have mapped the CO(3-2) line emission around the Mira AB system at 0.5 resolution using the Atacama Large Millimeter/submillimeter Array (ALMA). The CO map shows amazing complexity. The circumstellar gas has been shaped by different dynamical actors during the evolution of the system and several morphological components can be identified. The companion is marginally resolved in continuum emission and is currently at 0.487$\pm$0.006 separation. In the main line component, centered on the stellar velocity, spiral arcs around Mira A are found. The spiral appears to be relatively flat and oriented in the orbital plane. An accretion wake behind the companion is clearly visible and the projected arc separation is of order 5''. In the blue wing of the line emission, offset from the main line, several large ($\sim$5-10''), opposing arcs are found. We tentatively suggest that this structure is created by the wind of Mira B blowing a bubble in the expanding envelope of Mira A.

[abridged] The class of water fountain nebulae is thought to represent the stage of the earliest onset of collimated bipolar outflows during the post-Asymptotic Giant Branch phase. They thus play a crucial role in the study of the formation of bipolar Planetary Nebulae (PNe). To date, 14 water fountain nebulae have been identified. The identification of more sources in this unique stage of stellar evolution will enable us to study the origin of bipolar PNe morphologies in more detail. We present the results of seven sources observed with the JVLA that were identified as water fountain candidates in an Effelsberg 100m telescope survey of 74 AGB and early post-AGB stars. We find that our sample of water fountain candidates displays strong variability in their 22 GHz H2O maser spectra. The JVLA observations show an extended bipolar H2O maser outflow for one source, the OH/IR star IRAS 18455+0448. This source was previously classified as a 'dying' OH/IR star based on the exponential decrease of its 1612 MHz OH maser and the lack of H2O masers. We therefore also re-observed the 1612, 1665, and 1667 MHz OH masers. We confirm that the 1612 MHz masers have not reappeared and find that the 1665/1667 MHz masers have decreased in strength by several orders of magnitude during the last decade. The OH/IR star IRAS 18455+0448 is confirmed to be a new addition to the class of water fountain nebulae. Its kinematic age is approximately 70 yr, but could be lower, depending on the distance and inclination. Previous observations indicate, with significant uncertainty, that IRAS 18455+0448 has a surprisingly low mass compared to available estimates for other water fountain nebulae. The available historical OH maser observations make IRAS 18455+0448 unique for the study of water fountain nebulae and the launch of post-AGB bipolar outflows...

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

We report the first detections of OH$^+$ emission in planetary nebulae (PNe). As part of an imaging and spectroscopy survey of 11 PNe in the far-IR using the PACS and SPIRE instruments aboard the Herschel Space Observatory, we performed a line survey in these PNe over the entire spectral range between 51 and 672$\mu$m to look for new detections. OH$^+$ rotational emission lines at 152.99, 290.20, 308.48, and 329.77$\mu$m were detected in the spectra of three planetary nebulae: NGC 6445, NGC 6720, and NGC 6781. Excitation temperatures and column densities derived from these lines are in the range of 27 to 47 K and 2$\times$10$^{10}$ to 4 $\times$10$^{11}$ cm$^{-2}$, respectively. In PNe, the OH+ rotational line emission appears to be produced in the photodissociation region (PDR) in these objects. The emission of OH+ is observed only in PNe with hot central stars (T$_{eff}$ > 100000 K), suggesting that high-energy photons may play a role in the OH+ formation and its line excitation in these objects, as it seems to be the case for ultraluminous galaxies.

W75N(B) is a massive star-forming region that contains three radio continuum sources (VLA 1, VLA 2, and VLA 3), which are thought to be three massive young stellar objects at three different evolutionary stages. VLA 1 is the most evolved and VLA 2 the least evolved source. The 22 GHz H2O masers associated with VLA 1 and VLA 2 have been mapped at several epochs over eight years. While the H2O masers in VLA 1 show a persistent linear distribution along a radio jet, those in VLA 2 are distributed around an expanding shell. Furthermore, H2O maser polarimetric measurements revealed magnetic fields aligned with the two structures. Using new polarimetric observations of H2O masers, we aim to confirm the elliptical expansion of the shell-like structure around VLA 2 and, at the same time, to determine if the magnetic fields around the two sources have changed. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman-splitting of the 22 GHz H2O masers towards the massive star-forming region W75N(B). The H2O maser distribution around VLA 1 is unchanged from that previously observed. We made an elliptical fit of the H2O masers around VLA 2. We find that the shell-like structure is still expanding along the direction parallel to the thermal radio jet of VLA 1. While the magnetic field around VLA 1 has not changed in the past 7 years, the magnetic field around VLA 2 has changed its orientation according to the new direction of the major-axis of the shell-like structure and it is now aligned with the magnetic field in VLA 1.

The globular cluster (GC) M15 (NGC 7078) is host to at least eight pulsars and two low mass X-ray binaries (LMXBs) one of which is also visible in the radio regime. Here we present the results of a multi-epoch global very long baseline interferometry (VLBI) campaign aiming at i) measuring the proper motion of the known compact radio sources, ii) finding and classifying thus far undetected compact radio sources in the GC, and iii) detecting a signature of the putative intermediate mass black hole (IMBH) proposed to reside at the core of M15. We measure the sky motion in right ascension ($\mu_{\alpha}$) and declination ($\mu_{\delta}$) of the pulsars M15A and M15C and of the LMXB AC211 to be $(\mu_{\alpha},\,\mu_{\delta})_{\text{M15A}}=(-0.54\pm0.14,\,-4.33\pm0.25)\,$mas$\,$yr$^{-1}$, $(\mu_{\alpha},\,\mu_{\delta})_{\text{M15C}}=(-0.75\pm0.09,\,-3.52\pm0.13)\,$mas$\,$yr$^{-1}$, and $(\mu_{\alpha},\,\mu_{\delta})_{\text{AC211}}=(-0.46\pm0.08,\,-4.31\pm0.20)\,$mas$\,$yr$^{-1}$, respectively. Based on these measurements we estimate the global proper motion of M15 to be $(\mu_{\alpha},\,\mu_{\delta})=(-0.58\pm0.18,\,-4.05\pm0.34)\,$mas$\,$yr$^{-1}$. We detect two previously known but unclassified compact sources within our field of view. Our observations indicate that one them is of extragalactic origin while the other one is a foreground source, quite likely an LMXB. The double neutron star system M15C became fainter during the observations, disappeared for one year and is now observable again---an effect possibly caused by geodetic precession. The LMXB AC211 shows a double lobed structure in one of the observations indicative of an outburst during this campaign. With the inclusion of the last two of a total of seven observations we confirm the upper mass limit for a putative IMBH to be M$_{\bullet}<500$ M$_{\odot}$.

This is the first of a series of investigations into far-IR characteristics of 11 planetary nebulae (PNs) under the Herschel Space Observatory Open Time 1 program, Herschel Planetary Nebula Survey (HerPlaNS). Using the HerPlaNS data set, we look into the PN energetics and variations of the physical conditions within the target nebulae. In the present work, we provide an overview of the survey, data acquisition and processing, and resulting data products. We perform (1) PACS/SPIRE broadband imaging to determine the spatial distribution of the cold dust component in the target PNs and (2) PACS/SPIRE spectral-energy-distribution (SED) and line spectroscopy to determine the spatial distribution of the gas component in the target PNs. For the case of NGC 6781, the broadband maps confirm the nearly pole-on barrel structure of the amorphous carbon-richdust shell and the surrounding halo having temperatures of 26-40 K. The PACS/SPIRE multi-position spectra show spatial variations of far-IR lines that reflect the physical stratification of the nebula. We demonstrate that spatially-resolved far-IR line diagnostics yield the (T_e, n_e) profiles, from which distributions of ionized, atomic, and molecular gases can be determined. Direct comparison of the dust and gas column mass maps constrained by the HerPlaNS data allows to construct an empirical gas-to-dust mass ratio map, which shows a range of ratios with the median of 195+-110. The present analysis yields estimates of the total mass of the shell to be 0.86 M_sun, consisting of 0.54 M_sun of ionized gas, 0.12 M_sun of atomic gas, 0.2 M_sun of molecular gas, and 4 x 10^-3 M_sun of dust grains. These estimates also suggest that the central star of about 1.5 M_sun initial mass is terminating its PN evolution onto the white dwarf cooling track.

The analysis of astronomical interferometric data is often performed on the images obtained after deconvolution of the interferometer's point spread function (PSF). This strategy can be understood (especially for cases of sparse arrays) as fitting models to models, since the deconvolved images are already non-unique model representations of the actual data (i.e., the visibilities). Indeed, the interferometric images may be affected by visibility gridding, weighting schemes (e.g., natural vs. uniform), and the particulars of the (non-linear) deconvolution algorithms. Fitting models to the direct interferometric observables (i.e., the visibilities) is preferable in the cases of simple (analytical) sky intensity distributions. In this paper, we present UVMULTIFIT, a versatile library for fitting visibility data, implemented in a Python-based framework. Our software is currently based on the CASA package, but can be easily adapted to other analysis packages, provided they have a Python API. We have tested the software with synthetic data, as well as with real observations. In some cases (e.g., sources with sizes smaller than the diffraction limit of the interferometer), the results from the fit to the visibilities (e.g., spectra of close by sources) are far superior to the output obtained from the mere analysis of the deconvolved images. UVMULTIFIT is a powerful improvement of existing tasks to extract the maximum amount of information from visibility data, especially in cases close to the sensitivity/resolution limits of interferometric observations.

IRAS20126+4104 is a well studied B0.5 protostar that is surrounded by a ~1000 au Keplerian disk and is where a large-scale outflow originates. Both 6.7-GHz CH3OH masers and 22-GHz H2O masers have been detected toward this young stellar object. The CH3OH masers trace the Keplerian disk, while the H2O masers are associated with the surface of the conical jet. Recently, observations of dust polarized emission (350 um) at an angular resolution of 9 arcseconds (~15000 au) have revealed an S-shaped morphology of the magnetic field around IRAS20126+4104. The observations of polarized maser emissions at milliarcsecond resolution (~20 au) can make a crucial contribution to understanding the orientation of the magnetic field close to IRAS20126+4104. This will allow us to determine whether the magnetic field morphology changes from arcsecond resolution to milliarcsecond resolution. The European VLBI Network was used to measure the linear polarization and the Zeeman splitting of the 6.7-GHz CH3OH masers toward IRAS20126+4104. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman splitting of the 22-GHz H2O masers toward the same region. We detected 26 CH3OH masers and 5 H2O masers at high angular resolution. Linear polarization emission was observed toward three CH3OH masers and toward one H2O maser. Significant Zeeman splitting was measured in one CH3OH maser (\Delta V_{Z}=-9.2 +/- 1.4 m/s). No significant (5 sigma) magnetic field strength was measured using the H2O masers. We found that in IRAS20126+4104 the rotational energy is less than the magnetic energy.

We report the first detection of submillimeter water maser emission toward water-fountain nebulae, which are post-AGB stars that exhibit high-velocity water masers. Using APEX we found emission in the ortho-H2O (10_29-9_36) transition at 321.226 GHz toward three sources: IRAS 15445-5449, IRAS 18043-2116 and IRAS 18286-0959. Similarly to the 22 GHz masers, the submillimeter water masers are expanding with a velocity larger than that of the OH masers, suggesting that these masers also originate in fast bipolar outflows. In IRAS 18043-2116 and IRAS 18286-0959, which figure among the sources with the fastest water masers, the velocity range of the 321 GHz masers coincides with that of the 22 GHz masers, indicating that they likely coexist. Towards IRAS 15445-5449 the submillimeter masers appear in a different velocity range, indicating that they are tracing different regions. The intensity of the submillimeter masers is comparable to that of the 22 GHz masers, implying that the kinetic temperature of the region where the masers originate should be Tk > 1000 K. We propose that the passage of two shocks through the same gas can create the conditions necessary to explain the presence of strong high-velocity 321 GHz masers coexisting with the 22 GHz masers in the same region.

Galactic chemical evolution (GCE) is important for understanding the composition of the present-day interstellar medium (ISM) and of our solar system. In this paper, we aim to track the GCE by using the 29Si/30Si ratios in evolved stars and tentatively relate this to presolar grain composition. We used the APEX telescope to detect thermal SiO isotopologue emission toward four oxygen-rich M-type stars. Together with the data retrieved from the Herschel science archive and from the literature, we were able to obtain the 29Si/30Si ratios for a total of 15 evolved stars inferred from their optically thin 29SiO and 30SiO emission. These stars cover a range of masses and ages, and because they do not significantly alter 29Si/30Si during their lifetimes, they provide excellent probes of the ISM metallicity (or 29Si/30Si ratio) as a function of time. The 29Si/30Si ratios inferred from the thermal SiO emission tend to be lower toward low-mass oxygen-rich stars (e.g., down to about unity for W Hya), and close to an interstellar or solar value of 1.5 for the higher-mass carbon star IRC+10216 and two red supergiants. There is a tentative correlation between the 29Si/30Si ratios and the mass-loss rates of evolved stars, where we take the mass-loss rate as a proxy for the initial stellar mass or current stellar age. This is consistent with the different abundance ratios found in presolar grains. We found that older objects (up to possibly 10 Gyr old) in our sample trace a previous, lower 29Si/30Si value of about 1. Material with this isotopic ratio is present in two subclasses of presolar grains, providing independent evidence of the lower ratio. Therefore, the 29Si/30Si ratio derived from the SiO emission of evolved stars is a useful diagnostic tool for the study of the GCE and presolar grains.

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

(Abridged) The evolution of low- and intermediate-initial-mass stars beyond the asymptotic giant branch (AGB) remains poorly understood. High-velocity outflows launched shortly after the AGB phase are thought to be the primary shaping mechanism of bipolar and multipolar planetary nebulae (PNe). However, little is known about the launching and driving mechanism for these jets, whose momentum and energy often far exceed the energy that can be provided by radiation pressure alone. Here, we report direct evidence of a magnetically collimated jet shaping the bipolar morphology of the circumstellar envelope of a post-AGB star. We present radio continuum observations of the post-AGB star IRAS 15445-5449 which has water masers tracing a fast bipolar outflow. Our observations confirm the earlier observed steep negative spectral index of the spectral energy distribution (SED) above 3 GHz, and resolve, for the first time, the emission to originate from a radio jet, proving the existence of such jets around a post-AGB star. The SED is consistent with a synchrotron jet embedded in a sheath of thermal electrons. We find a close correspondence between the extent and direction of the synchrotron jet and the bipolar shape of the object observed at other wavelenghts, suggesting that the jet is responsible for the source morphology. The jet is collimated by a magnetic field of the order of mG at almost 7000 AU from the central star. We recover observations from the ATCA archive that indicate that the emission measure of the thermal component has increased by a factor of three between 1998 and 2005 after which it has remained constant. The short timescale evolution of the radio emission suggests a short lifetime for the jet.

The Boomerang Nebula is the coldest known object in the Universe, and an extreme member of the class of Pre-Planetary Nebulae, objects which represent a short-lived transitional phase between the AGB and Planetary Nebula evolutionary stages. Previous single-dish CO (J=1-0) observations (with a 45 arcsec beam) showed that the high-speed outflow in this object has cooled to a temperature significantly below the temperature of the cosmic background radiation. Here we report the first observations of the Boomerang with ALMA in the CO J=2-1 and J=1-0 lines to resolve the structure of this ultra-cold nebula. We find a central hourglass-shaped nebula surrounded by a patchy, but roughly round, cold high-velocity outflow. We compare the ALMA data with visible-light images obtained with HST and confirm that the limb-brightened bipolar lobes seen in these data represent hollow cavities with dense walls of molecular gas and dust producing both the molecular-emission-line and scattered-light structures seen at millimeter and visible wavelengths. The large diffuse biconical shape of the nebula seen in the visible is likely due to preferential illumination of the cold high-velocity outflow. We find a compact source of millimeter-wave continuum in the nebular waist -- these data, together with sensitive upper limits on the radio continuum using observations with ATCA, indicate the presence of a substantial mass of very large (mm-sized) grains in the waist of the nebula. Another unanticipated result is the detection of CO emission regions beyond the ultracold region which indicate the re-warming of the cold gas, most likely due to photoelectric grain heating.

[abridged] The 12CO/13CO ratio is often used as a measure of the 12C/13C ratio in the circumstellar environment, carrying important information about the stellar nucleosynthesis. External processes can change the 12CO and 13CO abundances, and spatially resolved studies of the 12CO/13CO ratio are needed to quantify the effect of these processes on the globally determined values. Additionally, such studies provide important information on the conditions in the circumstellar environment. The detached-shell source R Scl, displaying CO emission from recent mass loss, in a binary-induced spiral structure as well as in a clumpy shell produced during a thermal pulse, provides a unique laboratory for studying the differences in CO isotope abundances throughout its recent evolution. We observed both the 12CO(J=3-2) and the 13CO(J=3-2) line using ALMA. We find significant variations in the 12CO/13CO intensity ratios and consequently in the abundance ratios. The average CO isotope abundance ratio is at least a factor three lower in the shell (~19) than that in the present-day (<300 years) mass loss (>60). Additionally, variations in the ratio of more than an order of magnitude are found in the shell itself. We attribute these variations to the competition between selective dissociation and isotope fractionation in the shell, of which large parts cannot be warmer than ~35 K. However, we also find that the 12CO/13CO ratio in the present-day mass loss is significantly higher than the 12C/13C ratio determined in the stellar photosphere from molecular tracers (~19). The origin of this discrepancy is still unclear, but we speculate that it is due to an embedded source of UV-radiation that is primarily photo-dissociating 13CO. This radiation source could be the hitherto hidden companion. Alternatively, the UV-radiation could originate from an active chromosphere of R Scl itself....

Magnetic fields have only recently been included in theoretical simulations of high-mass star formation. The simulations show that magnetic fields play an important role in the formation and dynamics of molecular outflows. Masers, in particular 6.7-GHz CH3OH masers, are the best probes of the magnetic field morphologies around massive young stellar objects on the smallest scales of 10-100 AU. This paper focuses on 4 massive young stellar objects, IRAS06058+2138-NIRS1, IRAS22272+6358A, S255-IR, and S231, which complement our previous 2012 sample (the first EVN group). From all these sources, molecular outflows have been detected in the past. Seven of the European VLBI Network antennas were used to measure the linear polarization and Zeeman-splitting of the 6.7-GHz CH3OH masers in the star-forming regions in this second EVN group. We detected a total of 128 CH3OH masing cloudlets. Fractional linear polarization (0.8%-11.3%) was detected towards 18% of the CH3OH masers in our sample. The linear polarization vectors are well ordered in all the massive young stellar objects. We measured significant Zeeman-splitting in IRAS06058+2138-NIRS1 (DVz=3.8+/-0.6 m/s) and S255-IR (DVz=3.2+/-0.7 m/s). By considering the 20 massive young stellar objects towards which the morphology of magnetic fields was determined by observing 6.7-GHz CH3OH masers in both hemispheres, we find no evident correlation between the linear distributions of CH3OH masers and the outflows or the linear polarization vectors. On the other hand, we present first statistical evidence that the magnetic field (on scales 10-100 AU) is primarily oriented along the large-scale outflow direction. Moreover, we empirically find that the linear polarization fraction of unsaturated CH3OH masers is P_l<4.5%.

CRL 618 is an object that exhibits characteristics of both AGB and post-AGB star. It also displays a spectacular array of bipolar lobes with a dense equatorial region, which makes it an excellent object to study the development of asymmetries in evolved stars. In the recent decades, an elliptical compact HII region located in the center of the nebula has been seen to be increasing in size and flux. This seems to be due to the ionization of the circumstellar envelope by the central star, and it would be indicating the beginning of the planetary nebula phase for CRL 618. We analyzed interferometric radio continuum data at ~5 and 22 GHz from observations carried out at seven epochs with the VLA. We traced the increase of the flux of the ionized region over a period of ~26 years. We measured the dimensions of the HII region directly from the brightness distribution images to determine the increase of its size over time. For one of the epochs we analyzed observations at six frequencies from which we estimated the electron density distribution. We carried out model calculations of the spectral energy distribution at two different epochs to corroborate our observational results. We found that the radio continuum flux and the size of the ionized region have been increasing monotonically in the last three decades. The size of the major axis of the HII region shows a dependance with frequency, which has been interpreted as a result of a gradient of the electron density in this direction. The growth of the HII region is due to the expansion of an ionized wind whose mass-loss rate increased continuously for a period of ~100 years until a few decades ago, when the mass-loss rate experienced a sudden decline. Our results indicate that the beginning of the ionization of the circumstellar envelope began around 1971, which marks the start of the planetary nebula phase of CRL 618.

We aim to detect the magnetic field and infer its properties around four AGB stars using H$_2$O maser observations. The sample we observed consists of the following sources: the semi-regular variable RT Vir and the Mira variables AP Lyn, IK Tau, and IRC+60370. We observed the 6$_{1,6}-5_{2,3}$ H$_2$O maser rotational transition, in full-polarization mode, to determine its linear and circular polarization. Based on the Zeeman effect, one can infer the properties of the magnetic field from the maser polarization analysis. We detected a total of 238 maser features, in three of the four observed sources. No masers were found toward AP Lyn. The observed masers are all located between 2.4 and 53.0 AU from the stars. Linear and circular polarization was found in 18 and 11 maser features, respectively. We more than doubled the number of AGB stars in which magnetic field has been detected from H$_2$O maser polarization, as our results confirm the presence of fields around IK Tau, RT Vir and IRC+60370. The strength of the field along the line of sight is found to be between 47 and 331 mG in the H$_2$O maser region. Extrapolating this result to the surface of the stars, assuming a toroidal field ($\propto$ r$^{-1}$), we find magnetic fields of 0.3-6.9 G on the stellar surfaces. If, instead of a toroidal field, we assume a poloidal field ($\propto$ r$^{-2}$), then the extrapolated magnetic field strength on the stellar surfaces are in the range between 2.2 and $\sim$115 G. Finally, if a dipole field ($\propto$ r$^{-3}$) is assumed, the field strength on the surface of the star is found to be between 15.8 and $\sim$1945 G. The magnetic energy of our sources is higher than the thermal and kinetic energy in the H$_2$O maser region of this class of objects. This leads us to conclude that, indeed, magnetic fields probably play an important role in shaping the outflows of evolved stars. (abridged)

(Abridged) The goal of this work is to investigate if SiO, H2O and HCN maser emission within the ALMA frequency range can be detected with observable levels of fractional linear polarization in the regime where the Zeeman frequency (gOmega) is greater than the stimulated emission rate (R) and the decay rate of the molecular state (Gamma). We have used a radiative transfer code to calculate the fractional linear polarization as a function of the emerging brightness temperature for a number of rotational transition of SiO, H2O and HCN which have been observed to display maser emission at submillimetre wavelengths. We assume typical magnetic field strengths measured towards galactic star forming regions and circumstellar envelopes of late-type stars from previous VLBI observations. Since the Land\'e g-factors have not been reported for the different rotational transitions we have modeled, we performed our calculations assuming conservative values of the Zeeman frequency for the different molecular species. Setting a lower limit for the Zeeman frequency which still satisfies the criteria gOmega>R and gOmega>Gamma, we find fractional polarization levels of up to 13%, 14% and 19% for the higher J transitions analysed for SiO, H2O and HCN, respectively, without considering anisotropic pumping or any other non-Zeeman effect. Such upper limits were calculated assuming a magnetic field oriented perpendicular to the direction of propagation of the maser radiation. According to our results SiO, H2O and HCN maser emission within the ALMA frequency range can be detected with suitable linear polarization to trace the magnetic field structure towards star forming regions and late-type stars, even if the detected polarization has been enhanced by non-Zeeman effects.

The asymptotic giant branch star R Sculptoris is surrounded by a detached shell of dust and gas. The shell originates from a thermal pulse during which the star undergoes a brief period of increased mass loss. It has hitherto been impossible to constrain observationally the timescales and mass-loss properties during and after a thermal pulse - parameters that determine the lifetime on the asymptotic giant branch and the amount of elements returned by the star. Here we report observations of CO emission from the circumstellar envelope and shell around R Sculptoris with an angular resolution of 1.3 arcsec. What was hitherto thought to be only a thin, spherical shell with a clumpy structure, is revealed to contain a spiral structure. Spiral structures associated with circumstellar envelopes have been seen previously, from which it was concluded that the systems must be binaries. Using the data, combined with hydrodynamic simulations, we conclude that R Sculptoris is a binary system that underwent a thermal pulse approximately 1800 years ago, lasting approximately 200 years. About 0.003 Msun of mass was ejected at a velocity of 14.3 km s-1 and at a rate approximately 30 times higher than the prepulse mass-loss rate. This shows that approximately 3 times more mass is returned to the interstellar medium during and immediately after a pulse than previously thought.

Magnetic fields have been measured around Asymptotic Giant Branch (AGB) stars of all chemical types using maser polarization observations. If present, a large-scale magnetic field would lead to X-ray emission, which should be observable using current X-ray observatories. The aim of this work is to search the archival data for AGB stars that are intrinsic X-ray emitters. We have searched the ROSAT, CXO, and XMM archives for serendipitous X-ray observations of a sample of ~500 AGB stars. We specifically searched for the AGB stars detected with GALEX. The data is calibrated, analyzed and the X-ray luminosities and temperatures are estimated as functions of the circumstellar absorption. We identify 13 AGB stars as having either serendipitous or targeted observations in the X-ray data archives, however for a majority of the sources the detailed analysis show that the detections are questionable. Two new sources are detected by ROSAT: T Dra and R UMa. The spectral analysis suggests that the emission associated with these sources could be due to coronal activity or interaction across a binary system. Further observations of the detected sources are necessary to clearly determine the origin of the X-ray emission. Moreover, additional objects should be subject to targeted X-ray observations in order to achieve better constraints for the magnetic fields around AGB stars.

We present a new publicly available tool (DustPol) aimed to model the polarised thermal dust emission. The module DustPol, which is publicly available, is part of the ARTIST (Adaptable Radiative Transfer Innovations for Submillimetre Telescopes) package, which also offers tools for modelling the polarisation of line emission together with a model library and a Python-based user interface. DustPol can easily manage analytical as well as pre-gridded models to generate synthetic maps of the Stokes I, Q, and U parameters. These maps are stored in FITS format which is straightforwardly read by the data reduction software used, e.g., by the Atacama Large Millimeter Array (ALMA). This turns DustPol into a powerful engine for the prediction of the expected polarisation features of a source observed with ALMA or the Planck satellite as well as for the interpretation of existing submillimetre observations obtained with other telescopes. DustPol allows the parameterisation of the maximum degree of polarisation and we find that, in a prestellar core, if there is depolarisation, this effect should happen at densities of 10^6 cm-3 or larger. We compare a model generated by DustPol with the observational polarisation data of the low-mass Class 0 object NGC 1333 IRAS 4A, finding that the total and the polarised emission are consistent.

We present spectro-polarimetric observations of several molecular lines obtained with the Submillimeter Array (SMA) toward the carbon rich AGB star IRC+10216. We have detected and mapped the linear polarization of the CO 3-2, SiS 19-18 and CS 7-6 lines. The polarization arises at a distance of ~450 AU from the star and is blueshifted with respect the Stokes I. The SiS 19-18 polarization pattern appears to be consistent with a locally radial magnetic field configuration. However, the CO 3-2 and CS 7-6 line polarization suggests an overall complex magnetic field morphology within the envelope. This work demonstrates the feasibility of using spectro-polarimetric observations to carry out tomographic imaging of the magnetic field in circumstellar envelopes.

Intermediate mass black holes (IMBHs) with expected masses M_BH ~ 10^4 M_sun are thought to bridge the gap between stellar mass black holes (M_BH ~ 3 - 100 M_sun) and supermassive black holes found at the centre of galaxies (M_BH > 10^6 M_sun). Until today, no IMBH has been confirmed observationally. The most promising objects to host an IMBH as their central mass are globular clusters. Here, we present high sensitivity multi-epoch 1.6 GHz very long baseline interferometry observations of the globular cluster M15 that has been suggested to host an IMBH. Assuming the IMBH to be accreting matter from its surrounding we expect to detect it as a point source moving with the global motion of the cluster. However, we do not detect any such object within a radius of 6000 AU of the cluster centre in any of the five observations spread over more than one year. This rules out any variability of the putative IMBH on the time scale of one to two months. To get the most stringent upper limit for the flux density of the putative IMBH we concatenate the data of all five epochs. In this data we measure a 3{\sigma} upper flux limit of 10 {\mu}Jy for a central source. We employ the fundamental plane of black hole activity to estimate the mass of the central IMBH candidate. Based on previous X-ray observations of M15 our measurements indicate a 3{\sigma} upper mass limit of ~500 M_sun.

Shock-induced H2O masers are important magnetic field tracers at very high density gas. Water masers are found in both high- and low-mass star-forming regions, acting as a powerful tool to compare magnetic field morphologies in both mass regimes. In this paper, we show one of the first magnetic field determinations in the low-mass protostellar core IRAS 16293-2422 at volume densities as high as 10^(8-10) cm^-3. Our goal is to discern if the collapsing regime of this source is controlled by magnetic fields or other factors like turbulence. We used the Very Large Array (VLA) to carry out spectro-polarimetric observations in the 22 GHz Zeeman emission of H2O masers. From the Stokes V line profile, we can estimate the magnetic field strength in the dense regions around the protostar. A blend of at least three maser features can be inferred from our relatively high spatial resolution data set (~ 0.1"), which is reproduced in a clear non-Gaussian line profile. The emission is very stable in polarization fraction and position angle across the channels. The maser spots are aligned with some components of the complex outflow configuration of IRAS 16293-2422, and they are excited in zones of compressed gas produced by shocks. The post-shock particle density is in the range of 1-3 x 10^9 cm^-3, consistent with typical water masers pumping densities. Zeeman emission is produced by a very strong line-of-sight magnetic field (B ~ 113 mG). The magnetic field pressure derived from our data is comparable to the ram pressure of the outflow dynamics. This indicates that the magnetic field is energetically important in the dynamical evolution of IRAS 16293-2422.

The role of magnetic fields in the formation of high-mass stars is still under debate, and recent measurements of their orientation and strength by using polarized maser emissions are contributing new insights. Masers polarization, in particular of the 6.7-GHz methanol masers, are one of the best probes of the magnetic field morphologies around massive protostars. Determining the magnetic field morphology around an increasing number of massive protostars at milliarcsecond resolution by observing 6.7-GHz methanol masers is crucial to better understand the role of magnetic fields in massive star formation.The First EVN Group consists of 4 massive star-forming complexes: W51, W48, IRAS18556+0138, and W3(OH). These contain well-studied \hii ~regions from some of which molecular bipolar outflows were also detected (W51-e2, G35.20-0.74N). Nine of the European VLBI Network antennas were used to measure the linear polarization and Zeeman-splitting of the 6.7-GHz methanol masers in the star-forming regions of the First EVN Group. We detected a total of 154 CH3OH masers, one third of these towards W3(OH). Fractional linear polarization (1.2-11.5%) was detected towards 55 masers. The linear polarization vectors are well-ordered in all the massive star-forming regions. We measured significant Zeeman-splitting in 3 massive star-forming regions (W51, W48, and W3(OH)) revealing a range of separations -3.5 m/s<\Delta V_{z}<3.8 m/s with the smallest |\Delta V_{z}|=0.4m/s. We were also able to compare our magnetic field results with those obtained from submillimeter wavelength dust observation in W51 and show that the magnetic field at low and high resolutions are in perfect agreement.

Molecular line polarization is a unique source of information about the magnetic fields and anisotropies in the circumstellar envelopes of evolved stars. Here we present the first detection of thermal CO(J=2-1) and SiO(J=5-4, v=0) polarization, in the envelope of the asymptotic giant branch star IK Tau. The observed polarization direction does not match predictions for circumstellar envelope polarization induced only by an anisotropic radiation field. Assuming that the polarization is purely due to the Goldreich-Kylafis effect, the linear polarization direction is defined by the magnetic field as even the small Zeeman splitting of CO and SiO dominates the molecular collisional and spontaneous emission rates. The polarization was mapped using the Submillimeter Array (SMA) and is predominantly north-south. There is close agreement between the CO and SiO observations, even though the CO polarization arises in the circumstellar envelope at ~800 AU and the SiO polarization at <250 AU. If the polarization indeed traces the magnetic field, we can thus conclude that it maintains a large-scale structure throughout the circumstellar envelope. We propose that the magnetic field, oriented either east-west or north-south is responsible for the east-west elongation of the CO distribution and asymmetries in the dust envelope. In the future, the Atacama Large Millimeter/submillimeter Array will be able to map the magnetic field using CO polarization for a large number of evolved stars.

Most of PNe are not spherical. The loss of spherical symmetry occurs somewhere between the AGB and PN phase. The cause of this change of morphology is not yet well known, but magnetic fields are one of the possible agents. Its origin remains to be determined, and potentially requires the presence of a massive companion to the AGB star. Therefore, further detections of the magnetic field around evolved stars (in particular those thought to be part of a binary system) are crucial to improve our understanding of the origin and role of magnetism on evolved stars. One such binaries is the pre-PN OH231.8, around which a magnetic field was detected in the OH maser region of the outer circumstellar envelope. We aim to detect and infer the properties of the magnetic field of this source in the water maser region. We observed the 6_{1,6}-5_{2,3} water maser rotational transition to determine its linear and circular polarization. These emissions are located within the inner regions of OH231.8 (at few tens of AU). We detected 30 water maser features, which occur in two distinct regions that are moving apart with a velocity on the sky of 2.3 mas/year. Taking into account the inclination angle of the source, this corresponds to an average separation velocity of 21 km/s. Based on the velocity gradient of the maser emission, the masers appear to be dragged along the direction of the nebula jet. Linear polarization is present in 3 of the features, and circular polarization was detected in the 2 brightest ones. We found that the strength of the magnetic field is |B_{||}|~45 mG which, when assuming a toroidal magnetic field, implies B~2.5 G on the stellar surface. The morphology of the field is not yet determined, but the high scatter found on the directions of the linear polarization vectors could indicate that the masers occur near the tangent points of a toroidal field.

(ABRIDGED) We present multi-epoch Very Long Baseline Array (VLBA) observations of V773 Tau A, the 51-day binary subsystem in the multiple young stellar system V773 Tau. Combined with previous interferometric and radial velocity measurements, these new data enable us to improve the characterization of the physical orbit of the A subsystem. In particular, we infer updated dynamical masses for the primary and the secondary components of 1.55 pm 0.11 Msun, and 1.293 pm 0.068 Msun, respectively, and an updated orbital parallax distance to the system of 135.7 pm 3.2 pc, all consistent with previous estimates. Using the improved orbit, we can calculate the absolute coordinates of the barycenter of V773 Tau A at each epoch of our VLBA observations, and fit for its trigonometric parallax and proper motion. This provides a direct measurement of the distance to the system almost entirely independent of the orbit modeling. The best fit yields a distance of 129.9 pm 3.2 pc, in good agreement (i.e. within 1 sigma) with the distance estimate based on the orbital fit. Taking the mean value of the orbital and trigonometric parallaxes, we conclude that V773 Tau is located at d =132.8 pm 2.3 pc. The accuracy of this determination is nearly one order of magnitude better than that of previous estimates. In projection, V773 Tau and two other young stars (Hubble 4 and HDE 283572) recently observed with the VLBA are located toward the dark cloud Lynds 1495, in the central region of Taurus. These three stars appear to have similar trigonometric parallaxes, radial velocities, and proper motions, and we argue that the weighted mean and dispersion of their distances (d = 131.4 pc and sigma_d = 2.4 pc) provide a good estimate of the distance to and depth of Lynds 1495 and its associated stellar population.

SiO maser emission occurs in the extended atmosphere of evolved stars and can be studied at high angular resolution. Very long baseline interferometry (VLBI) observations of Mira variables indicate that SiO masers are significantly linearly polarized with linear polarization fraction up to 100%. However, no information is available at high angular resolution for SiO masers in higher mass loss OH/IR stars. We extend the VLBI SiO maser studies to OH/IR stars. The observations enable us to understand the SiO pumping mechanisms in higher mass loss evolved objects and compare those with Mira variables. Additionally, polarimetric observations of SiO masers help us to understand the magnetic field strength and morphology and to distinguish between conflicting polarization theories. The 43 GHz SiO maser observations of the OH/IR star OH 44.8-2.3 were performed with the VLBA in full polarization spectral line mode. Auxiliary EVLA observations were performed to allow for the absolute calibration of the polarization angle. The Zeeman splitting was measured by cross correlating the right and left circular polarization spectra as well as the S-curve fitting. Additionally, we analyzed the 1612 MHz OH maser observations of OH 44.8-2.3 from the VLA archive. The SiO masers of OH 44.8-2.2 form a ring located at ~5.4 AU around the star. The masers appear to be highly linearly polarized with fractional linear polarization up to 100%. The linear polarization vectors are consistent with a dipole field morphology in this star. We report a tentative detection of circular polarization of ~0.7% for the brightest maser feature. The magnetic field measured for this feature corresponds to 1.5\pm0.3 G. Additionally, the distribution of the 1612 MHz OH maser emission could indicate an elongated morphology.

We present the morphology and linear polarization of the 22-GHz H2O masers in the high-velocity outflow of two post-AGB sources, d46 (IRAS 15445-5449) and b292 (IRAS 18043-2116). The observations were performed using The Australia Telescope Compact Array. Different levels of saturated maser emission have been detected for both sources. We also present the mid-infrared image of d46 overlaid with the distribution of the maser features that we have observed in the red-shifted lobe of the bipolar structure. The relative position of the observed masers and a previous radio continuum observation suggests that the continnum is produced along the blue-shifted lobe of the jet. It is likely due to synchrontron radiation, implying the presence of a strong magnetic field in the jet. The fractional polarization levels measured for the maser features of d46 indicate that the polarization vectors are tracing the poloidal component of the magnetic field in the emitting region. For the H2O masers of b292 we have measured low levels of fractional linear polarization. The linear polarization in the H2O maser region of this source likely indicates a dominant toroidal or poloidal magnetic field component. Since circular polarization was not detected it is not possible to determine the magnetic field strength. However, we present a 3-sigma evaluation of the upper limit intensity of the magnetic field in the maser emitting regions of both observed sources.

NGC7538 is a complex massive star-forming region. The region is composed of several radio continuum sources, one of which is IRS1, a high-mass protostar, from which a 0.3 pc molecular bipolar outflow was detected. Several maser species have been detected around IRS1. The CH3OH masers have been suggested to trace a Keplerian-disk, while the H2O masers are almost aligned to the outflow. More recent results suggested that the region hosts a torus and potentially a disk, but with a different inclination than the Keplerian-disk that is supposed to be traced by the CH3OH masers. Tracing the magnetic field close to protostars is fundamental for determining the orientation of the disk/torus. Recent studies showed that during the protostellar phase of high-mass star formation the magnetic field is oriented along the outflows and around or on the surfaces of the disk/torus. The observations of polarized maser emissions at milliarcsecond resolution can make a crucial contribution to understanding the orientation of the magnetic field and, consequently, the orientation of the disk/torus in NGC7538-IRS1. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman-splitting of the 22GHz H2O masers toward NGC7538-IRS1. The European VLBI Network and the MERLIN telescopes were used to measure the linear polarization and the Zeeman-splitting of the 6.7GHz CH3OH masers toward the same region. We detected 17 H2O masers and 49 CH3OH masers at high angular resolution. We detected linear polarization emission toward two H2O masers and toward twenty CH3OH masers. The CH3OH masers, most of which only show a core structure, seem to trace rotating and potentially infalling gas in the inner part of a torus. Significant Zeeman-splitting was measured in three CH3OH masers. [...] We also propose a new description of the structure of the NGC7538-IRS1 maser region.

The aim of this study is to investigate the distribution and excitation of the methanol gas around Cep A and to describe the physical conditions of the region. In addition the large scale abundance distribution is determined in order to understand the morphology and kinematics of star forming regions in which methanol masers occur. The spatial distribution of the methanol is studied by mapping line emission, as well as the column density and excitation temperature, which are estimated using rotation diagrams. For a limited number of positions the parameters are checked with non-LTE models. Furthermore, the distribution of the methanol abundance is derived in comparison with archival dust continuum maps. Methanol is detected over a 0.3x0.15 pc area centred on the Cep A HW2 source, showing an outflow signature. Most of the gas can be characterized by a moderately warm rotation temperature (30-60K). At the central position two velocity components are detected with different excitation characteristics, the first related to the large-scale outflow. The second component, uniquely detected at the central location, is probably associated with the maser emission on much smaller scales of 2". Detailed analysis reveals that the highest densities and temperatures occur for these inner components. In the inner region the dust and gas are shown to have different physical parameters. Abundances of methanol in the range 10E-9 - 10E-7 are inferred, with the abundance peaking at the maser position. The geometry of the large-scale methanol is in accordance with previous determinations of the Cep A geometry, in particular those from methanol masers. The dynamical and chemical time-scales are consistent with the methanol originating from a single driving source associated with the HW2 object and the masers in its equatorial region.

In recent studies of methanol masers, a substantial fraction of the objects show maser components aligned in large-scale elliptical configurations. These can be readily interpreted as rings centred on a high mass star in formation, seen in projection. Remarkably, most of these rings do not show signs of rotation, but rather the radial motions dominate. This must mean that their dynamics are governed by other than gravitational forces. In particular, we have studied the methanol masers around Cep A in detail, where it can be argued that the methanol masers show signs of infall. In this paper we discuss the dynamics of the Cep A methanol maser and sources from the Torun blind survey to argue that at least in a fraction of sources methanol masers could be associated with the shock interface between the large scale accretion, regulated by the magnetic field, and a 1000-AU scale circumstellar disk. We discuss the validity of such a model for the overall population of methanol maser sources.

We present Submillimeter Array observations of high frequency SiO masers around the supergiant VX Sgr and the semi-regular variable star W Hya. The J=5-4, v=1 28SiO and v=0 29SiO masers of VX Sgr are shown to be highly linearly polarized with a polarization from ~5-60%. Assuming the continuum emission peaks at the stellar position, the masers are found within ~60 mas of the star, corresponding to ~100 AU at a distance of 1.57 kpc. The linear polarization vectors are consistent with a large scale magnetic field, with position and inclination angles similar to that of the dipole magnetic field inferred in the H2O and OH maser regions at much larger distances from the star. We thus show for the first time that the magnetic field structure in a circumstellar envelope can remain stable from a few stellar radii out to ~1400 AU. This provides further evidence supporting the existence of large scale and dynamically important magnetic fields around evolved stars. Due to a lack of parallactic angle coverage, the linear polarization of masers around W Hya could not be determined. For both stars we observed the 28SiO and 29SiO isotopologues and find that they have a markedly different distribution and that they appear to avoid each other. Additionally, emission from the SO 5_5-4_4 line was imaged for both sources. Around W Hya we find a clear offset between the red- and blue-shifted SO emission. This indicates that W Hya is likely host to a slow bipolar outflow or a rotating disk-like structure.

A debated topic in star formation theory is the role of magnetic fields during the protostellar phase of high-mass stars. It is still unclear how magnetic fields influence the formation and dynamics of massive disks and outflows. Most current information on magnetic fields close to high-mass protostars comes from polarized maser emissions, which allows us to investigate the magnetic field on small scales by using very long-baseline interferometry. The massive star-forming region W75N contains three radio continuum sources (VLA1, VLA2, and VLA3), at three different evolutionary stages, and associated masers, while a large-scale molecular bipolar outflow is also present. Very recently, polarization observations of the 6.7 GHz methanol masers at milliarsecond resolution have been able to probe the strength and structure of the magnetic field over more than 2000 AU around VLA1. The magnetic field is parallel to the outflow, suggesting that VLA1 is its powering source. The observations of water masers at 22 GHz can give more information about the gas dynamics and the magnetic fields around VLA1 and VLA2. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman-splitting of the 22 GHz water masers in the star-forming region W75N. We detected 124 water masers, 36 around VLA1 and 88 around VLA2 of W75N, which indicate two different physical environments around the two sources, where VLA1 is in a more evolved state. The linear polarization of the masers confirms the tightly ordered magnetic field around VLA1, which is aligned with the large-scale molecular outflow, and also reveals an ordered magnetic field around VLA2, which is not parallel to the outflow. [abridged]

The 6.7 GHz methanol maser is exclusively associated with high-mass star formation. However, it remains unclear what structures harbour the methanol masers. Cepheus A is one of the closest regions of massive star formation, making it an excellent candidate for detailed studies. We determine the dynamics of maser spots in the high-mass star-forming region Cepheus A in order to infer where and when the maser emission occurs. Very long baseline interferometry (VLBI) observations of the 6.7 and 12.2 GHz methanol masers allows for mapping their spatial and velocity distribution. Phase-referencing is used to determine the astrometric positions of the maser emission, and multi-epoch observations can reveal 3D motions. The 6.7 GHz methanol masers are found in a filamentary structure over ~1350 AU, straddling the waist of the radio jet HW2. The positions agree well with previous observations of both the 6.7 and 12.2 GHz methanol masers. The velocity field of the maser spots does not show any sign of rotation, but is instead consistent with an infall signature. The 12.2 GHz methanol masers are closely associated with the 6.7 GHz methanol masers, and the parallax that we derive confirms previous measurements. We show that the methanol maser emission very likely arises in a shock interface in the equatorial region of Cepheus A HW2 and presents a model in which the maser emission occurs between the infalling gas and the accretion disk/process.

To characterize the magnetic field structure of the outflow and core region within a prototypical high-mass star-forming region, we analyzed polarized CO(3-2) - for the first time observed with the Submillimeter Array -- as well as 880mum submm continuum emission from the high-mass outflow/disk system IRAS18089-1732. Both emission features with polarization degrees at a few percent level indicate that the magnetic field structure is largely aligned with the outflow/jet orientation from the small core scales to the larger outflow scales. Although quantitative estimates are crude, the analysis indicates that turbulent energy dominates over magnetic energy. The data also suggest a magnetic field strength increase from the lower-density envelope to the higher-density core.

Observational evidence for strong magnetic fields throughout the envelopes of evolved stars is increasing. Many of the instruments coming on line in the near-future will be able to make further contributions to this field. Specifically, maser polarization observations and dust/line polarization in the sub-mm regime has the potential to finally provide a definite picture of the magnetic field strength and configuration from the Asymptotic Giant Branch (AGB) all the way to the Planetary Nebula phase. While current observations are limited in sample size, strong magnetic fields appear ubiquitous at all stages of (post-)AGB evolution. Recent observations also strongly support a field structure that is maintained from close to the star to several thousands of AU distance. While its origin is still unclear, the magnetic field is thus a strong candidate for shaping the stellar outflows on the path to the planetary nebula phase and might even play a role in determining the stellar mass-loss.

We present five epochs of VLBI water maser observations around the massive protostar Cepheus A HW2 with 0.4 mas (0.3 AU) resolution. The main goal of these observations was to follow the evolution of the remarkable water maser linear/arcuate structures found in earlier VLBI observations. Comparing the data of our new epochs of observation with those observed five years before, we find that at "large" scales of > 1" (700 AU) the main regions of maser emission persist, implying that both the surrounding medium and the exciting sources of the masers have been relatively stable during that time span. However, at smaller scales of < 0.1" (70 AU) we see large changes in the maser structures, particularly in the expanding arcuate structures R4 and R5. R4 traces a nearly elliptical patchy ring of ~ 70 mas size (50 AU) with expanding motions of ~ 5 mas/yr (15 km/s). This structure is probably driven by the wind of a still unidentified YSO located at the centre of the ring (~ 0.18" south of HW2). On the other hand, the R5 expanding bubble structure (driven by the wind of a previously identified YSO located ~ 0.6" south of HW2) is currently dissipating in the circumstellar medium and losing its previous degree of symmetry, indicating a very short-lived event. In addition, our results reveal, at scales of ~ 1" (700 AU), the simultaneous presence of a relatively slow (~ 10-70 km/s) wide-angle outflow (opening angle of ~ 102 deg, traced by the masers, and the fast (~ 500~km/s) highly collimated radio jet associated with HW2 (opening angle of ~ 18 deg, previously observed with the VLA. This simultaneous presence of a wide-angle outflow and a highly collimated jet associated with a massive protostar is similar to what is found in some low-mass YSOs. The implications of these results in the study of the formation of high-mass stars are discussed.

We present polarization observations of the 6.7-GHz methanol masers around the massive protostar Cepheus A HW2 and its associated disc. The data were taken with the Multi-Element Radio Linked Interferometer Network. The maser polarization is used to determine the full three-dimensional magnetic field structure around Cepheus A HW2. The observations suggest that the masers probe the large scale magnetic field and not isolated pockets of a compressed field. We find that the magnetic field is predominantly aligned along the protostellar outflow and perpendicular to the molecular and dust disc. From the three-dimensional magnetic field orientation and measurements of the magnetic field strength along the line of sight, we are able to determine that the high density material, in which the masers occurs, is threaded by a large scale magnetic field of ~23 mG. This indicates that the protostellar environment at ~1000 AU from Cepheus A HW2 is slightly supercritical (lambda approximately 1.7) and the relation between density and magnetic field is consistent with collapse along the magnetic field lines. Thus, the observations indicate that the magnetic field likely regulates accretion onto the disc. The magnetic field dominates the turbulent energies by approximately a factor of three and is sufficiently strong to be the crucial component stabilizing the massive accretion disc and sustaining the high accretion rates needed during massive star-formation.

The majority of the observed planetary nebulae exhibit elliptical or bipolar structures. Theoretical modeling has indicated that magnetically collimated jets may be responsible for the formation of the non-spherical planetary nebulae. The aim of this project is to measure the Zeeman splitting caused by the magnetic field in the OH and H2O maser regions occurring in the circumstellar envelope and bipolar outflow of the evolved star W43A. We report a measured magnetic field of approximately 100 micro-gauss in the OH maser region of the circumstellar envelope around W43A. The GBT observations reveal a magnetic field strength B|| of ~30 mG changing sign across the H2O masers at the tip of the red-shifted lobe of the bipolar outflow. We also find that the OH maser shell shows no sign of non-spherical expansion and that it probably has an expansion velocity that is typical for the shells of regular OH/IR stars. The GBT observations confirm that the magnetic field collimates the H2O maser jet, while the OH maser observations show that a strong large scale magnetic field is present in the envelope surrounding the W43A central star. The magnetic field in the OH maser envelope is consistent with the one extrapolated from the H2O measurements, confirming that magnetic fields play an important role in the entire circumstellar environment of W43A.

The role of magnetic fields during the protostellar phase of high-mass star-formation is a debated topic. In particular, it is still unclear how magnetic fields influence the formation and dynamic of disks and outflows. Most current information on magnetic fields close to high-mass protostars comes from H2O and OH maser observations. Recently, the first 6.7 GHz methanol maser polarization observations were made, and they reveal strong and ordered magnetic fields. The morphology of the magnetic field during high-mass star-formation needs to be investigated on small scales, which can only be done using very long baseline interferometry observations. The massive star-forming regionW75N contains three radio sources and associated masers, while a large-scale molecular bipolar outflow is also present. Polarization observations of the 6.7 GHz methanol masers at high angular resolution probe the strength and structure of the magnetic field and determine its relation to the outflow. Eight of the European VLBI network antennas were used to measure the linear polarization and Zeeman-splitting of the 6.7 GHz methanol masers in the star-forming region W75N. We detected 10 methanol maser features, 4 of which were undetected in previous work. All arise near the source VLA1 of W75N. The linear polarization of the masers reveals a tightly ordered magnetic field over more than 2000 AU around VLA1 that is exactly aligned with the large-scale molecular outflow. This is consistent with the twisted magnetic field model proposed for explaining dust polarization observations. The Zeeman-splitting measured on 3 of the maser features indicates a dynamically important magnetic field in the maser region of the order of 50mG. We suggest VLA1 is the powering sources of the bipolar outflow.

(Abridged) Recently, the magnetic field induced Zeeman splitting was measured for the strongest known 6.7 GHz methanol maser, which arises in the massive star forming region G09.62+0.20. This maser is one of a handful of periodically flaring methanol masers. The 100-m Effelsberg telescope was used to monitor the 6.7 GHz methanol masers of G09.62+0.20. With the exception of a two week period during the peak of the maser flare, we measure a constant magnetic field of B_||~11+-2 mG in the two strongest maser components of G09.62+0.20 that are separated by over 200 AU. In the two week period that coincides exactly with the peak of the maser flare of the strongest maser feature, we measure a sharp decrease and possible reversal of the Zeeman splitting. The exact cause of both maser and polarization variability is still unclear, but it could be related to either background amplification of polarized emission or the presence of a massive protostar with a close-by companion. Alternatively, the polarization variability could be caused by non-Zeeman effects related to the radiative transfer of polarized maser emission.

Astrometry can bring powerful constraints to bear on a variety of scientific questions about neutron stars, including their origins, astrophysics, evolution, and environments. Using phase-referenced observations at the VLBA, in conjunction with pulsar gating and in-beam calibration, we have measured the parallaxes and proper motions for 14 pulsars. The smallest measured parallax in our sample is 0.13+-0.02 mas for PSR B1541+09, which has a most probable distance of 7.2+1.3-1.1 kpc. We detail our methods, including initial VLA surveys to select candidates and find in-beam calibrators, VLBA phase-referencing, pulsar gating, calibration, and data reduction. The use of the bootstrap method to estimate astrometric uncertainties in the presence of unmodeled systematic errors is also described. Based on our new model-independent estimates for distance and transverse velocity, we investigate the kinematics and birth sites of the pulsars and revisit models of the Galactic electron density distribution. We find that young pulsars are moving away from the Galactic plane, as expected, and that age estimates from kinematics and pulsar spindown are generally in agreement, with certain notable exceptions. Given its present trajectory, the pulsar B2045-16 was plausibly born in the open cluster NGC 6604. For several high-latitude pulsars, the NE2001 electron density model underestimates the parallax distances by a factor of two, while in others the estimates agree with or are larger than the parallax distances, suggesting that the interstellar medium is irregular on relevant length scales. The VLBA astrometric results for the recycled pulsar J1713+0747 are consistent with two independent estimates from pulse timing, enabling a consistency check between the different reference frames.

Observations of H$_2$O masers towards the post-AGB star and water fountain source OH 009.1--0.4 were made as part of HOPS (The H$_2$O southern galactic Plane Survey), with the Mopra radiotelescope. Together with followup observations using the Australia Telescope Compact Array (ATCA), we have identified H$_2$O maser emission over a velocity spread of nearly 400km/s (--109 to +289km/s). This velocity spread appears to be the largest of any known maser source in our Galaxy. High resolution observations with the ATCA indicate the maser emission is confined to a region $0\farcs3 \times 0\farcs3$ and shows weak evidence for a separation of the red- and blueshifted maser spots. We are unable to determine if the water fountain is projected along the line of sight, or is inclined, but either way OH 009.1--0.4 is an interesting source, worthy of followup observations.

Context: Magnetic fields are thought to be one of the possible mechanisms responsible for shaping the generally spherical outflow of evolved stars into often aspherical planetary nebulae. However, direct measurements of magnetic fields during the transition to the planetary nebula phase are rare. Aims: The aim of this project is to expand the number of magnetic field measurements of stars in the (proto-)planetary nebula phase and find if the magnetic field strength is sufficient to affect the stellar outflow. Methods: We used Very Long Baseline Array observations to measure the circular polarization due to the Zeeman splitting of 22 GHz water masers in the envelope of the proto-planetary nebula candidate star IRAS 19296+2227 and the planetary nebula K3-35. Results: A strong magnetic field of B||=-135+-28 is detected in the water maser region of the proto-planetary nebula candidate IRAS 19296+2227. The water masers of K3-35 are too weak to detect circular polarization although we do present the measurements of weak linear polarization in those masers. Conclusions: The field measured in the masers of IRAS 19296+2227 is dynamically important and, if it is representative of the large scale field, is an important factor in driving the stellar mass loss and shaping the stellar outflow.

Binary pulsar systems are superb probes of stellar and binary evolution and the physics of extreme environments. In a survey with the Arecibo telescope, we have found PSR J1903+0327, a radio pulsar with a rotational period of 2.15 ms in a highly eccentric (e = 0.44) 95-day orbit around a solar mass companion. Infrared observations identify a possible main-sequence companion star. Conventional binary stellar evolution models predict neither large orbital eccentricities nor main-sequence companions around millisecond pulsars. Alternative formation scenarios involve recycling a neutron star in a globular cluster then ejecting it into the Galactic disk or membership in a hierarchical triple system. A relativistic analysis of timing observations of the pulsar finds its mass to be 1.74+/-0.04 Msun, an unusually high value.

Context: The role of magnetic fields during high-mass star formation is a matter of fierce debate, yet only a few direct probes of magnetic field strengths are available. Aims: The magnetic field is detected in a number of massive star-forming regions through polarization observations of 6.7 GHz methanol masers. Although these masers are the most abundant of the maser species occurring during high-mass star formation, most magnetic field measurements in the high-density gas currently come from OH and H2 O maser observations. Methods: The 100-m Effelsberg telescope was used to measure the Zeeman splitting of 6.7 GHz methanol masers for the first time. The observations were performed on a sample of 24 bright northern maser sources. Results: Significant Zeeman splitting is detected in 17 of the sources with an average magnitude of 0.56 m/s . Using the current best estimate of the 6.7 GHz methanol maser Zeeman splitting coefficient and a geometrical correction, this corresponds to an absolute magnetic field strength of 23 mG in the methanol maser region. Conclusions: The magnetic field is dynamically important in the dense maser regions. No clear relation is found with the available OH maser magnetic field measurements. The general sense of direction of the magnetic field is consistent with other Galactic magnetic field measurements, although a few of the masers display a change of direction between different maser features. Due to the abundance of methanol masers, measuring their Zeeman splitting provides the opportunity to construct a comprehensive sample of magnetic fields in high-mass star-forming regions.

MERLIN observations of 6.668-GHz methanol and both 6.031- and 6.035-GHz hydroxyl (OH) emission from the massive star-formation region ON1 are presented. These are the first methanol observations made in full polarization using 5 antennas of MERLIN, giving high resolution and sensitivity to extended emission. Maser features are found to lie at the southern edge of the ultra-compact HII region, following the known distribution of ground-state OH masers. The masers cover a region ~1 arcsec in extent, lying perpendicular to the H13CO+ bipolar outflow. Excited-state OH emission demonstrates consistent polarization angles across the strongest linearly polarized features which are parallel to the overall distribution. The linear polarizations vary between 10.0 and 18.5 per cent, with an average polarization angle of -60 deg +/- 28 deg. The strongest 6.668-GHz methanol features provide an upper limit to linear polarization of ~1 per cent. Zeeman splitting of OH shows magnetic fields between -1.1 to -5.8 mG, and a tentative methanol magnetic field strength of -18 mG is measured.

Aims: Accurate distances to evolved stars with high mass loss rates are needed for studies of many of their fundamental properties. However, as these stars are heavily obscured and variable, optical and infrared astrometry is unable to provide enough accuracy. Methods: Astrometry using masers in the circumstellar envelopes can be used to overcome this problem. We have observed the OH masers of a number of Asymptotic Giant Branch (AGB) stars for approximately 1 year with the Very Long Baseline Array (VLBA). We have used the technique of phase referencing with in-beam calibrators to test the improvements this technique can provide to Very Long Baseline Interferometry (VLBI) OH maser astrometric observations. Results: We have significantly improved the parallax and proper motion measurements of the Mira variable stars U Her, S CrB and RR Aql. Conclusions: It is shown that both in-beam phase-referencing and a decrease in solar activity during the observations significantly improves the accuracy of the astrometric observations. The improved distances to S CrB (418 +21 -18 pc) and RR Aql (633 +214 -128 pc) are fully consistent with published P-L relations, but the distance to U Her (266 +32 -28 pc) is significantly smaller. We conclude that for sources that are bright and have a nearby in-beam calibrator, VLBI OH maser astrometry can be used to determine distances to OH masing stars of up to ~2 kpc.

Through polarization observations masers are unique probes of the magnetic field in a variety of different astronomical objects, with the different maser species tracing different physical conditions. In recent years maser polarization observations have provided insights in the magnetic field strength and morphology in, among others, the envelopes around evolved stars, Planetary Nebulae (PNe), massive star forming regions and supernova remnants. More recently, maser observations have even been used to determine the magnetic field in megamaser galaxies. This review will present an overview of maser polarization observations and magnetic field determinations of the last several years and discuss the implications of the magnetic field measurements for several important fields of study, such as aspherical PNe creation and massive star formation.

We present observations of the 22 GHz water megamasers in the circumnuclear disk of NGC 3079 obtained with the Green Bank Telescope. The data are analyzed for circular polarization due to the Zeeman-induced splitting of the water maser lines. No circular polarization is detected and we derive a 1 sigma upper limit of 11 mG for the toroidal magnetic field at ~0.64 pc from the central black hole. This is the tightest upper limit for the magnetic field around a black hole to date. We use the magnetic field limit to derive an estimate of the mass accretion onto the central black hole. In addition to the polarimetric results, we present an observation of rapid variability in the maser lines, which we explain as weak interstellar scintillation. From the scintillation parameters, we estimate an intrinsic size of the mostly saturated maser features of ~12 microarcseconds. This is consistent with models assuming a thick, clumpy accretion disk.

Water maser polarization observations in the precessing jet of W43A have revealed that it is magnetically collimated. Here we present a detailed description of the physical properties of the water maser environment in the jet. We discuss the maser saturation level and beaming angle as well as the intrinsic temperatures and densities. Additionally, we show that the polarization angle of the strongest red-shifted maser feature undergoes a fast rotation of 90 degrees across the maser. Along with the variation of linear polarization fraction, this strongly supports the current theoretical description of maser linear polarization.

We present the first 6.7 GHz methanol maser linear polarization map of the extended filamentary maser structure around the compact HII region W3(OH). The methanol masers show linear polarization up to 8 per cent and the polarization angles indicate a magnetic field direction along the North-South maser structure. The polarization angles are consistent with those measured for the OH masers, taking into account external Faraday rotation toward W3(OH), and confirm that the OH and methanol masers are found in similar physical conditions. Additionally we discuss the Zeeman splitting of the 6.7 GHz methanol transition and present an upper limit of ~22 mG for the magnetic field strength in the maser region. The upper limit is fully consistent with the field strengths derived from OH maser Zeeman splitting.

We present the first direct measurements of the magnetic field strength and direction in a collimated jet from an evolved star on its way to become a planetary nebula. Very Long Baseline Array (VLBA) observations of the linear and circular polarization of the water masers in the collimated jet of W43A reveal a strong toroidal magnetic field, indicating that the jet is magnetically collimated. The magnetic field strength in the jet extrapolated back to the stellar surface yields a surface field of several Gauss, consistent with the measurements of maser polarization in a large sample of evolved stars. The origin of the magnetic field is yet unknown, although the jet precession might point to the existence of a heavy planet or stellar companion. This is the first direct observational evidence for magnetic collimation in the jets, that likely play an important role in shaping planetary nebulae.

The recently started Arecibo L-band Feed Array (ALFA) pulsar survey aims to find ~1000 new pulsars. Due to its high time and frequency resolution the survey is especially sensitive to millisecond pulsars, which have the potential to test gravitational theories, detect gravitational waves and probe the neutron-star equation of state. Here we report the results of our preliminary analysis: in the first months we have discovered 21 new pulsars. One of these, PSR J1906+0746, is a young 144-ms pulsar in a highly relativistic 3.98-hr low-eccentricity orbit. The 2.61 +- 0.02 solar-mass system is expected to coalesce in ~300 Myr and contributes significantly to the computed cosmic inspiral rate of compact binary systems.

Planetary nebulae often have asymmetric shapes, which could arise due to collimated jets from evolved stars before evolution to the planetary nebula phase. The source of jet collimation in these stars is unknown. Magnetic fields are thought to collimate outflows that are observed in many other astrophysical sources, such as active galactic nuclei and proto-stars, although hitherto there are no direct observations of both the magnetic field direction and strength in any collimated jet. Theoretical models have shown that magnetic fields could also be the dominant source of collimation of jet in evolved stars. Here we report measurements of the polarization of water vapour masers that trace the precessing jet emanating from the asymptotic giant branch star W43A at 2.6 kpc from the Sun, which is undergoing rapid evolution into a planetary nebula. The masers occur in two clusters at opposing tips of the jets, ~1,000 AU from the star. We find direct evidence that the magnetic field is collimating the jet.

We report the discovery of PSR J1906+0746, a young 144-ms pulsar in a highly relativistic 3.98-hr orbit with an eccentricity of 0.085 and expected gravitational wave coalescence time of 300 Myr. The new pulsar was found during precursor survey observations with the Arecibo 1.4-GHz feed array system and retrospectively detected in the Parkes Multibeam plane pulsar survey data. From radio follow-up observations with Arecibo, Jodrell Bank, Green Bank, and Parkes, we have measured the spin-down and binary parameters of the pulsar and its basic spectral and polarization properties. We also present evidence for pulse profile evolution, which is likely due to geodetic precession, a relativistic effect caused by the misalignment of the pulsar spin and total angular momentum vectors. Our measurements show that PSR J1906+0746 is a young object with a characteristic age of 112 kyr. From the measured rate of orbital periastron advance 7.57+/-0.03 deg/yr, we infer a total system mass of 2.61+/-0.02 Msun. While these parameters suggest that the PSR J1906+0746 binary system might be a younger version of the double pulsar system, intensive searches for radio pulses from the companion have so far been unsuccessful. It is therefore not known whether the companion is another neutron star or a massive white dwarf. Regardless of the nature of the companion, a simple calculation suggests that the Galactic birth rate of binaries similar to PSR J1906+0746 is 60/Myr. This implies that PSR J1906+0746 will make a significant contribution to the computed cosmic inspiral rate of compact binary systems.

We present linear and circular polarization observations of the water masers in 4 distinct regions spread over 1x2 arcseconds around the HW2 high-mass young stellar object in the Cepheus A star-forming region. We find magnetic fields between 100-500 mG in the central maser region, which has been argued to trace a circumstellar disk. The masers further from HW2 have field strengths between 30-100 mG. In all cases the magnetic field pressure is found to be similar to the dynamic pressure, indicating that the magnetic field is capable of controlling the outflow dynamics around HW2. In addition to several water maser complexes observed before, we also detect a new maser filament, ~1 arcsec (~690 AU) East of HW2, which we interpret as a shocked region between the HW2 outflow and the surrounding medium. We detect a linear polarization gradient along the filament as well as a reversal of the magnetic field direction. This is thought to mark the transition between the magnetic field associated with the outflow and that found in the surrounding molecular cloud. In addition to the magnetic field we determine several other physical properties of the maser region, including density and temperatures as well as the maser beaming angles.

We report results from the initial stage of a long-term pulsar survey of the Galactic plane using the Arecibo L-band Feed Array (ALFA), a seven-beam receiver operating at 1.4 GHz with 0.3 GHz bandwidth. The search targets Galactic latitudes |b| < 5 deg in the longitude ranges 32 deg < l < 77 deg and 168 deg < l < 77 deg. Data discussed here were collected over a 100 MHz passband centered on 1.42 GHz using a spectrometer that recorded 256 channels every 64 microsec. In a preliminary, standard period-DM analysis, we have detected 29 previously known pulsars and discovered 11 new ones. One of these, with a period of 69 ms and a low characteristic age of 82 kyr, is a plausible candidate for association with the unidentified EGRET source 3EG J1928+1733. Another is a non-recycled pulsar in a relativistic binary with orbital period of 3.98 hr. We also search the data for isolated dispersed pulses, a technique that yielded discovery of an extremely sporadic radio emitter with a spin period of 1.2 s. Simulations we have carried out indicate that about 1000 new pulsars will be found in the ALFA survey. In addition to providing a large sample for use in population analyses and for probing the magnetoionic interstellar medium, the survey maximizes the chances of finding rapidly spinning millisecond pulsars and pulsars in compact binary systems. Our search algorithms will exploit the multiple data streams from ALFA to discriminate between radio frequency interference and celestial signals, including pulsars and possibly new classes of transient radio sources.

Through polarization observations water masers are excellent probes of magnetic fields in the maser region. Magnetic field strengths, such as those in the water masers regions of the envelopes of late-type stars and star-forming regions, are typically determined using a direct relation between the field strength and the observed circular polarization. Here it is shown that velocity and magnetic field gradients along the maser have a significant effect on the field strengths obtained from circular polarization observations. Due to velocity gradients the actual magnetic field strength could be up to 100% higher than the field strength derived from the observations. Additionally, when a magnetic field gradient is present, the resulting circular polarization derived is caused predominantly by the average magnetic field in the unsaturated maser core. Measurements of the fractional linear polarization are not affected by velocity or magnetic field strength gradients, though changes in the magnetic field angle along the maser do quench the linear polarization intensity when the maser saturates.

Determination of pulsar parallaxes and proper motions addresses fundamental astrophysical questions. We have recently finished a VLBI astrometry project to determine the proper motions and parallaxes of 27 pulsars, thereby doubling the total number of pulsar parallaxes. Here we summarise our astrometric technique and present the discovery of a pulsar moving in excess of 1000 km/s. As an example of the application of high precision pulsar astrometry we also infer the identification of 2 pulsars originating from a disrupted binary in the Cygnus Superbubble.

The highest velocity neutron stars establish stringent constraints on natal kicks, asymmetries in supernova core collapse, and the evolution of close binary systems. Here we present the first results of a long-term pulsar astrometry program using the VLBA. We measure a proper motion and parallax for the pulsar B1508+55, leading to model-independent estimates of its distance (2.37+0.23-0.20 kpc) and transverse velocity (1083+103-90 km/s), the highest velocity directly measured for a neutron star. We trace the pulsar back from its present Galactic latitude of 52.3 degrees to a birth site in the Galactic plane near the Cyg OB associations, and find that it will inevitably escape the Galaxy. Binary disruption alone is insufficient to impart the required birth velocity, and a natal kick is indicated. A composite scenario including a large kick along with binary disruption can plausibly account for the high velocity.

We present the analysis of the circular polarization, due to Zeeman splitting, of the water masers around a sample of late-type stars to determine the magnetic fields in their circumstellar envelopes. The magnetic field strengths in the water maser regions around the Mira variable stars U Ori and U Her are shown to be several Gauss while those of the supergiants S Per, NML Cyg and VY CMa are several hundred mG. We also show that large scale magnetic fields permeate the CSE of an evolved star; the polarization of the water masers around VX Sgr reveals a dipole field structure. We shortly discuss the coupling of the magnetic field with the stellar outflow, as such fields could possibly be the cause of distinctly aspherical mass-loss and the resulting aspherical planetary nebulae.

Through polarization observations, circumstellar masers are excellent probes of the magnetic field in the envelopes of late-type stars. Whereas observations of the polarization of the SiO masers close to the star and on the OH masers much further out were fairly commonplace, observations of the magnetic field strength in the intermediate density and temperature region where the 22 GHz water masers occur have only recently become possible. Here we present the analysis of the circular polarization, due to Zeeman splitting, of the water masers around the Mira variable stars U Her and U Ori and the supergiant VX Sgr. We present an upper limit of the field around U Her that is lower but consistent with previous measurements, reflecting possible changes in the circumstellar envelope. The field strengths around U Ori and VX Sgr are shown to be of the order of several Gauss. Moreover, we show for the first time that large scale magnetic fields permeate the circumstellar envelopes of an evolved star; the polarization of the water masers around VX Sgr reveals a dipole field structure. We discuss the coupling of the magnetic field with the stellar outflow, as such fields could possibly be the cause of distinctly aspherical mass-loss.

The changes in the spectral line profile of the 22GHz H2O maser are calculated as a function of emerging maser flux. We address not only the narrowing and re-broadening of the maser lines, but also the deviations of Gaussian symmetry as a result of the various hyperfine components of the 22 GHz maser line. Non-LTE models of the H2O maser transition, including hyperfine structure and cross-relaxation, are compared with high spectral resolution Very Long Baseline Interferometry (VLBI) observations of the H2O masers in the circumstellar envelopes (CSEs) of a sample of late-type stars. This yields estimates on the thermal width in the maser region as well as the emerging maser flux and thus the level of saturation. We also discuss the effect of a velocity gradient along the maser path on the line widths and shapes of the line profile and the effect of the geometry of the maser region. We find that the typical velocity shift along the maser path is of the order of 1.0 km/s. The effect of this shift on the shape of the maser spectrum is difficult to distinguish from the effect of the hyperfine components.

Many observations of late-type M stars show large near-spherical circumstellar envelopes, yet planetary nebulae and supernova remnants are frequently axisymmetric. We present VLBI and MERLIN observations of masers around the red supergiant S Per which show varying degrees of axisymmetry and a dynamically significant magnetic field. There is no evidence for rotation here or in most similar objects. We examine possible origins of the magnetic field.

High precision astrometric data have enabled us to determine the trajectories through the Galactic potential for a growing number of pulsars. This has resulted in the discovery of a pulsar pair (B2020+28 and B2021+51) that has a common origin in the Cygnus Superbubble or in one of its related OB associations at an epoch which is comparable with the spin-down ages of the pulsars. Analysis of the Galactic orbits indicates that the progenitors of the pulsars had similar masses and were in a binary system, which was disrupted after the second supernova explosion. The implied pulsar birth velocities are consistent with the high velocities of neutron stars in general. The initial spin period of the pulsar that was formed in the second supernova explosion was approximately 200 ms. A further increase in astrometric accuracy will allow us to more tightly constrain the birth velocities and the kick velocities that were imparted by the two respective supernova explosions. Two additional pulsars in the sample of 24 with parallax measurements may also have originated in the Cygnus Superbubble.

We present the first pulsar parallaxes measured with phase-referenced pulsar VLBI observations at 5 GHz. Due to the steep spectra of pulsars, previous astrometric measurements have been at lower frequencies. However, the strongest pulsars can be observed at 5 GHz, offering the benefit of lower combined ionospheric and tropospheric phase errors, which usually limit VLBI astrometric accuracy. The pulsars B0329+54, B0355+54 and B1929+10 were observed for 7 epochs spread evenly over 2 years. For B0329+54, large systematic errors lead to only an upper limit on the parallax (pi < 1.5 mas). A new proper motion and parallax were measured for B0355+54 (pi = 0.91 +- 0.16 mas), implying a distance of 1.04+0.21-0.16 kpc and a transverse velocity of 61+12-9 km/s. The parallax and proper motion for B1929+10 were significantly improved (pi = 2.77 +- 0.07 mas), yielding a distance of 361+10-8 pc and a transverse velocity of 177+4-5 km/s. We demonstrate that the astrometric errors are correlated with the angular separation between the phase reference calibrator and the target source, with significantly lower errors at 5 GHz compared to 1.6 GHz. Finally, based on our new distance determinations for B1929+10 and B0355+54, we derive or constrain the luminosities of each pulsar at high energies. We show that, for thermal emission models, the emitting area for X-rays from PSR B1929+10 is roughly consistent with the canonical size for a heated polar cap, and that the conversion of spin-down power to gamma-ray luminosity in B0355+54 must be low. The new proper motion for B1929+10 also implies that its progenitor is unlikely to have been the binary companion of the runaway O-star zeta-Ophiuchi.

The OH main line masers in circumstellar shells were observed with astrometric VLBI in order to accurately determine the position of the star, as well as its motion, including the parallax. Results are presented on 4 Mira stars monitored for 2 - 8 years with the VLBA. The data show that in some stars the VLBI detection is dominated by blue-shifted emission that is associated with the stellar image amplified by the maser shell in front of the star. In other cases the maser is not directly tied to the stellar position, but still persistent enough to measure proper motion and parallax.

We present circular polarization measurements of circumstellar water masers around 2 Mira variables and 4 supergaint stars. Both an LTE and non-LTE method are used to analyze the circular polarization and total intensity spectra and to obtain the magnetic field strength. The non-LTE method is found to be able to reproduce the observations best. Using this, we find fields from a few hundred milliGauss up to a few Gauss, indicating a solar-type r^(-2) dependence of the magnetic field on the distance, which leads to stellar surface magnetic fields of up to several hundred Gauss. No linear polarization is detected to less than 1.5%.

The main-line OH masers around 4 AGB stars have been observed with the NRAO Very Long Baseline Array (VLBA) at 8 epochs over a period of 2.5 years. Using a phase referencing technique, the position of the most compact maser spot of each star was monitored with respect to two extragalactic reference sources. For U Her and W Hya, we observe the most blue-shifted maser spot, while for R Cas and S CrB we only detect a compact red-shifted maser spot. We managed to determine an accurate proper motion and parallax for U Her, R Cas and S CrB, while additional motion of the compact blue-shifted maser of W Hya is shown to possibly be related to the stellar pulsation. The motion and radio position are compared with the stellar trajectory and absolute optical position determined by the Hipparcos satellite. For U Her and W Hya, the most blue-shifted maser is consistent with the amplified stellar image. The new distances are compared with several published P--L relations, and in this respect the VLBI distances seem an improvement upon the Hipparcos distances.

The circumstellar masers around evolved stars offer an interesting possibility to measure stellar parameters through VLBI astrometry. In this paper the application of this technique is discussed, including the accuracy and the uncertainties of the method. The different maser species (OH, H_2O, SiO) have slightly different characteristics and applications. This paper does not concern astrometry of maser spots to study the kinematics of the envelope, but concentrates on attempting to measure the motion of the underlying star.

The 22 GHz H_2O masers in the circumstellar envelope of the Mira variable star U Her have been observed with MERLIN using a phase referencing technique to determine accurate astrometric positions. The positions were compared with the optical positions obtained with the Hipparcos satellite to an accuracy of 18 mas. The absolute radio position of the brightest H_2O maser spot is found to match the optical position, indicating that this spot is the stellar image amplified by the maser screen in front of it. The occurrence of an amplified image in the 22 GHz maser can be used to accurately determine the positions of the H_2O with respect to the star as well as with respect to the SiO and OH masers. Our observations seem to indicate that the star is not in the centre of the distribution of maser spots, which has been interpreted as a ring.

We present circular polarization measurements of circumstellar H_2O masers. The circular polarization detected in the (6_{16}-5_{23}) rotational transition of the H_{2}O maser can be attributed to Zeeman splitting in the intermediate temperature and density regime. The magnetic fields are derived using a general, LTE Zeeman analysis as well as a full radiative transfer method (non-LTE), which includes a treatment of all hyperfine components simultaneously as well as the effects of saturation and unequal populations of the magnetic substates. The differences and relevances of these interpretations are discussed extensively. The field strengths are compared with previous detections of the magnetic field on the SiO and OH masers. We show that the magnetic pressure dominates the thermal pressure by a factor of 20 or more.

We present the first results of our 3D Monte Carlo maser radiative transfer code, used to model the 1612 MHz OH maser shell and the amplification of emission from the stellar radio-photosphere.

We present the first circular polarization measurements of circumstellar water masers around a sample of late-type stars. These observations are used to obtain the magnetic field strength in the water maser region with both an LTE and non-LTE analysis. We find fields from a few hundred milliGauss up to a few Gauss, indicating a solar-type r^(-2) dependence of the magnetic field on the distance. No linear polarization is detected to less than 1%.

We present the first circular polarization measurements of circumstellar H2O masers. Previously the magnetic field in circumstellar envelopes has been estimated using polarization observations of SiO and OH masers. SiO masers are probes of the high temperature and density regime close to the central star. OH masers are found at much lower densities and temperatures, generally much further out in the circumstellar envelope. The detection of the circular polarization of the (6_16 - 5_23) rotational transition of the H2O maser could be attributed to Zeeman splitting due to the magnetic field in the intermediate temperature and density regime. The fields inferred here agree well with predicted values for a combination of the r^{-2} dependence of a solar-type magnetic field, and the coupling of the field to the high density masing regions. We also discuss the unexpected narrowing of the circular polarization spectrum.

The OH 1667 MHz maser in the circumstellar shell around the Mira variable U Herculis has been observed with the NRAO Very Long Baseline Array (VLBA) at 6 epochs, spread over 4 years. Using phase referencing techniques the position of the most blue-shifted maser spot was monitored with respect to two extra-galactic radio sources. The absolute radio positions of the maser can be compared with the stellar optical position measured by the Hipparcos satellite to 15 mas accuracy. This confirms the model in which one of the maser spots corresponds to the stellar continuum, amplified by the maser. The stellar proper motion and the annual parallax (5.3 +/- 2.1 mas) were measured.