Here we aim to characterise submillimetre water masers at 437, 439, 471, and 474 GHz towards a sample of evolved stars. We used the Atacama Pathfinder Experiment (APEX) to observe submillimetre water transitions and the CO (4-3) line towards 11 evolved stars. The sample included semi-regular and Mira variables, plus a red supergiant star. We performed radiative transfer modelling for the water masers. We also used the CO observations to determine mass loss rates for the stars. From the sample of 11 evolved stars, 7 display one or more of the masers at 437, 439, 471, and 474 GHz. We therefore find that these masers are common in evolved star circumstellar envelopes. The fact that the maser lines are detected near the stellar velocity indicates that they are likely to originate from the inner circumstellar envelopes of our targets. We tentatively link the presence of masers to the degree of variability of the target star, that is, masers are more likely to be present in Mira variables than in semi-regular variables. We suggest that this indicates the importance of strong shocks in creating the necessary conditions for the masers. Typically, the 437 GHz line is the strongest maser line observed among those studied here. We cannot reproduce the above finding in our radiative transfer models. In general, we find that maser emission is very sensitive to dust temperature in the lines studied here. To produce strong maser emission, the dust temperature must be significantly lower than the gas kinetic temperature. In addition to running grids of models in order to determine the optimum physical conditions for strong masers in these lines, we performed smooth wind modelling for which we cannot reproduce the observed line shapes. This also suggests that the masers must originate predominantly from the inner envelopes.

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

The astrophysical agenda of the 21st century requires a very sharp view of celestial objects. High angular resolution studies are essential for fundamental studies of a broad variety of astrophysical phenomena ranging from relativistic physics of black holes, their gravitational and electromagnetic imprints, violent transient processes, including those producing detectable gravitational waves, birth and evolution of planetary systems. Over the past decades, radio astronomy made huge leap in achieving ground-breaking angular resolution measured in tens of microarcseconds (one tenth of nanoradian and better). Recently a global Event Horizon Telescope (EHT) collaboration obtained first direct images of the shadow of a super-massive black hole in the nucleus of the active galaxy M87. These observations were conducted at 230 GHz. The two first generation Space Very Long Baseline Interferometry (VLBI) missions, VSOP/HALCA led by the Japan Aerospace Exploration Agency (JAXA) and RadioAstron led by the Russia Roscosmos State Corporation and Russia Academy of Sciences, achieved the highest angular resolution at frequencies from 0.3 to 22 GHz in observations conducted in the period 1997 - 2019. The next step in advancing high angular resolution radio astronomy is in combining high frequency (millimeter and sub-millimeter wavelengths) and interferometric baselines exceeding the Earth diameter. The present THEZA White Paper describes a combination which would unify technology developments in giga-/tera-hertz instrumentation and space-borne radio astronomy. The current preprint version of the THEZA White Paper is slightly re-formatted and edited comparing to the official submitted version.

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.

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 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.

We present ALMA Atacama Compact Array (ACA) 870 micron continuum maps of 28 infrared-bright SDSS quasars with Herschel/SPIRE detections at redshifts 2-4, the largest such sample ever observed with ALMA. The ACA detections are centred on the SDSS coordinates to within 1 arcsec for about 80 per cent of the sample. Larger offsets indicate that the far-infrared (FIR) emission detected by Herschel might come from a companion source. The majority of the objects (about 70 per cent) have unique ACA counterparts within the SPIRE beam down to 3-4 arcsec resolution. Only 30 per cent of the sample shows clear evidence for multiple sources with secondary counterparts contributing to the total 870 micron flux within the SPIRE beam to at least 25 per cent. We discuss the limitations of the data based on simulated pairs of point-like sources at the resolution of the ACA and present an extensive comparison of our findings with recent works on the multiplicities of sub-millimetre galaxies. We conclude that, despite the coarse resolution of the ACA, our data support the idea that, for a large fraction of FIR-bright quasars, the sub-mm emission comes from single sources. Our results suggest that, on average, optically bright quasars with strong FIR emission are not triggered by early-stage mergers but are, instead, together with their associated star formation rates, the outcome of either late-stage mergers or secular processes.

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.

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 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.

We present the first images of the 28SiO v=1, J=2-1 maser emission around the closest known massive young stellar object Orion Source I observed at 86 GHz (3mm) with the VLBA. These images have high spatial (~0.3 mas) and spectral (~0.054 km/s) resolutions. We find that the 3mm masers lie in an X-shaped locus consisting of four arms, with blue-shifted emission in the south and east arms and red-shifted emission in the north and west arms. Comparisons with previous images of the 28SiO v=1,2, J=1-0 transitions at 7mm (observed in 2001-2002) show that the bulk of the J=2-1 transition emission follows the streamlines of the J=1-0 emission and exhibits an overall velocity gradient consistent with the gradient at 7mm. While there is spatial overlap between the 3mm and 7mm transitions, the 3mm emission, on average, lies at larger projected distances from Source I (~44 AU compared with ~35 AU for 7mm). The spatial overlap between the v=1, J=1-0 and J=2-1 transitions is suggestive of a range of temperatures and densities where physical conditions are favorable for both transitions of a same vibrational state. However, the observed spatial offset between the bulk of emission at 3mm and 7mm possibly indicates different ranges of temperatures and densities for optimal excitation of the masers. We discuss different maser pumping models that may explain the observed offset. We interpret the 3mm and 7mm masers as being part of a single wide-angle outflow arising from the surface of an edge-on disk rotating about a northeast-southwest axis, with a continuous velocity gradient indicative of differential rotation consistent with a Keplerian profile in a high-mass proto-binary.

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.

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.

We present near-infrared interferometry of the carbon-rich asymptotic giant branch (AGB) star R Sculptoris. The visibility data indicate a broadly circular resolved stellar disk with a complex substructure. The observed AMBER squared visibility values show drops at the positions of CO and CN bands, indicating that these lines form in extended layers above the photosphere. The AMBER visibility values are best fit by a model without a wind. The PIONIER data are consistent with the same model. We obtain a Rosseland angular diameter of 8.9+-0.3 mas, corresponding to a Rosseland radius of 355+-55 Rsun, an effective temperature of 2640+-80 K, and a luminosity of log L/Lsun=3.74+-0.18. These parameters match evolutionary tracks of initial mass 1.5+-0.5 Msun and current mass 1.3+-0.7 Msun. The reconstructed PIONIER images exhibit a complex structure within the stellar disk including a dominant bright spot located at the western part of the stellar disk. The spot has an H-band peak intensity of 40% to 60% above the average intensity of the limb-darkening-corrected stellar disk. The contrast between the minimum and maximum intensity on the stellar disk is about 1:2.5. Our observations are broadly consistent with predictions by dynamic atmosphere and wind models, although models with wind appear to have a circumstellar envelope that is too extended compared to our observations. The detected complex structure within the stellar disk is most likely caused by giant convection cells, resulting in large-scale shock fronts, and their effects on clumpy molecule and dust formation seen against the photosphere at distances of 2-3 stellar radii.

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 present a line survey of the ultra-luminous infrared galaxy Arp 220, taken with the newly installed SEPIA Band 5 instrument on APEX. We illustrate the capacity of SEPIA to detect the 183.3 GHz H2O 31,3-22,0 line against the atmospheric H2O absorption feature. We confirm the previous detection of the HCN(2-1) line, and detect new transitions of standard dense gas tracers such as HNC(2-1), HCO+(2-1), CS(4-3), C34S(4-3), HC3N(20-19). We also detect HCN(2-1) v2=1 and the 193.5 GHz methanol (4-3) group for the first time. The absence of time variations in the megamaser water line compared to previous observations seems to rule out an AGN nuclear origin for the line. It could, on the contrary, favor a thermal origin instead, but also possibly be a sign that the megamaser emission is associated with star-forming cores washed-out in the beam. We finally discuss how the new transitions of HCN, HNC, HCO+ refine our knowledge of the ISM physical conditions in Arp 220.

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 present the results of models that were designed to study all possible water maser transitions in the frequency range 0-1.91THz, with particular emphasis on maser transitions that may be generated in evolved-star envelopes and observed with the ALMA and SOFIA telescopes. We used tens of thousands of radiative transfer models of both spin species of H2O, spanning a considerable parameter space in number density, kinetic temperature and dust temperature. Results, in the form of maser optical depths, have been summarized in a master table, Table 6. Maser transitions identified in these models were grouped according to loci of inverted regions in the density/kinetic temperature plane, a property clearly related to the dominant mode of pumping. A more detailed study of the effect of dust temperature on maser optical depth enabled us to divide the maser transitions into three groups: those with both collisional and radiative pumping schemes (22,96,209,321,325,395,941 and 1486\,GHz), a much larger set that are predominantly radiatively pumped, and another large group with a predominantly collisional pump. The effect of accelerative and decelerative velocity shifts of up to 5km/s was found to be generally modest, with the primary effect of reducing computed maser optical depths. More subtle asymmetric effects, dependent on line overlap, include maximum gains offset from zero shift by >1km/s, but these effects were predominantly found under conditions of weak amplification. These models will allow astronomers to use multi-transition water maser observations to constrain physical conditions down to the size of individual masing clouds (size of a few astronomical units).

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.

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.

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 present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum images of the asteroid 3 Juno obtained with an angular resolution of 0.042 arcseconds (60 km at 1.97 AU). The data were obtained over a single 4.4 hr interval, which covers 60% of the 7.2 hr rotation period, approximately centered on local transit. A sequence of ten consecutive images reveals continuous changes in the asteroid's profile and apparent shape, in good agreement with the sky projection of the three-dimensional model of the Database of Asteroid Models from Inversion Techniques. We measure a geometric mean diameter of 259pm4 km, in good agreement with past estimates from a variety of techniques and wavelengths. Due to the viewing angle and inclination of the rotational pole, the southern hemisphere dominates all of the images. The median peak brightness temperature is 215pm13 K, while the median over the whole surface is 197pm15 K. With the unprecedented resolution of ALMA, we find that the brightness temperature varies across the surface with higher values correlated to the subsolar point and afternoon areas, and lower values beyond the evening terminator. The dominance of the subsolar point is accentuated in the final four images, suggesting a reduction in the thermal inertia of the regolith at the corresponding longitudes, which are possibly correlated to the location of the putative large impact crater. These results demonstrate ALMA's potential to resolve thermal emission from the surface of main belt asteroids, and to measure accurately their position, geometric shape, rotational period, and soil characteristics.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations from the 2014 Long Baseline Campaign in dust continuum and spectral line emission from the HL Tau region. The continuum images at wavelengths of 2.9, 1.3, and 0.87 mm have unprecedented angular resolutions of 0.075 arcseconds (10 AU) to 0.025 arcseconds (3.5 AU), revealing an astonishing level of detail in the circumstellar disk surrounding the young solar analogue HL Tau, with a pattern of bright and dark rings observed at all wavelengths. By fitting ellipses to the most distinct rings, we measure precise values for the disk inclination (46.72pm0.05 degrees) and position angle (+138.02pm0.07 degrees). We obtain a high-fidelity image of the 1.0 mm spectral index ($\alpha$), which ranges from $\alpha\sim2.0$ in the optically-thick central peak and two brightest rings, increasing to 2.3-3.0 in the dark rings. The dark rings are not devoid of emission, we estimate a grain emissivity index of 0.8 for the innermost dark ring and lower for subsequent dark rings, consistent with some degree of grain growth and evolution. Additional clues that the rings arise from planet formation include an increase in their central offsets with radius and the presence of numerous orbital resonances. At a resolution of 35 AU, we resolve the molecular component of the disk in HCO+ (1-0) which exhibits a pattern over LSR velocities from 2-12 km/s consistent with Keplerian motion around a ~1.3 solar mass star, although complicated by absorption at low blue-shifted velocities. We also serendipitously detect and resolve the nearby protostars XZ Tau (A/B) and LkHa358 at 2.9 mm.

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.042 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using a very extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), corresponding to an un-magnified spatial scale of ~180 pc at z=3.042. The ALMA images clearly show two main gravitational arc components of an Einstein ring, with emission tracing a radius of ~1.5". We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO emission, at an angular resolution of ~170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies between the CO transitions and compared to the dust continuum. Our detection of H2O line emission, using only the shortest baselines, provides the most resolved detection to date of thermal H2O emission in an extragalactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing galaxy.

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.

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.

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.

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

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.

We performed a kinematical analysis of the [CII] line emission of the BR 1202-0725 system at z~4,7 using ALMA observations. The most prominent sources of this system are a quasar and a submillimeter galaxy, separated by a projected distance of about 24 kpc and characterized by very high SFR, higher than 1000 Msun/yr. However, the ALMA observations reveal that these galaxies apparently have undisturbed rotating disks, which is at variance with the commonly accepted scenario in which strong star formation activity is induced by a major merger. We also detected faint components which, after spectral deblending, were spatially resolved from the main QSO and SMG emissions. The relative velocities and positions of these components are compatible with orbital motions within the gravitational potentials generated by the QSO host galaxy and the SMG, suggesting that they are smaller galaxies in interaction or gas clouds in accretion flows of tidal streams. We did not find any clear spectral evidence for outflows caused by AGN or stellar feedback. This suggests that the high star formation rates might be induced by interactions or minor mergers with these companions, which do not affect the large-scale kinematics of the disks, however. Our kinematical analysis also indicates that the QSO and the SMG have similar Mdyn, mostly in the form of molecular gas, and that the QSO host galaxy and the SMG are seen close to face-on with slightly different disk inclinations: the QSO host galaxy is seen almost face-on (i~15), while the SMG is seen at higher inclinations (i~25). Finally, the ratio between the black hole mass of the QSO, obtained from XShooter spectroscopy, and the Mdyn of the host galaxy is similar to value found in very massive local galaxies, suggesting that the evolution of black hole galaxy relations is probably better studied with dynamical than with stellar host galaxy masses.

In Orion BN/KL, proper motions of 7 mm vibrationally-excited SiO masers trace rotation of a nearly edge-on disk and a bipolar wide-angle outflow 10-100 AU from radio Source I, a binary young stellar object (YSO) of ~20 Msun. Here we map ground-state 7 mm SiO emission with the Very Large Array and track proper motions over 9 years. The innermost and strongest emission lies in two extended arcs bracketing Source I. The proper motions trace a northeast-southwest bipolar outflow 100-1000 AU from Source I with a median 3D motion of ~18 km/s. An overlying distribution of 1.3 cm H2O masers betrays similar flow characteristics. Gas dynamics and emission morphology traced by the masers suggest the presence of a magnetocentrifugal disk-wind. Reinforcing evidence lies in the colinearity of the flow, apparent rotation across the flow parallel to the disk rotation, and recollimation that narrows the flow opening angle ~120 AU downstream. The arcs of ground-state SiO emission may mark the transition point to a shocked super-Alfvenic outflow.

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.

Using Science Verification data from the Atacama Large Millimeter/Submillimeter Array (ALMA), we have identified and imaged five rotational transitions (J=5-4 and J=6-5) of the three silicon monoxide isotopologues 28SiO v=0, 1, 2 and 29SiO v=0 and 28Si18O v=0 in the frequency range from 214 to 246 GHz towards the Orion BN/KL region. The emission of the ground-state 28SiO, 29SiO and 28Si18O shows an extended bipolar shape in the northeast-southwest direction at the position of Radio Source I, indicating that these isotopologues trace an outflow (~18 km/s, P.A. ~50deg, ~5000 AU in diameter) that is driven by this embedded high-mass young stellar object (YSO). Whereas on small scales (10-1000 AU) the outflow from Source I has a well-ordered spatial and velocity structure, as probed by Very Long Baseline Interferometry (VLBI) imaging of SiO masers, the large scales (500-5000 AU) probed by thermal SiO with ALMA reveal a complex structure and velocity field, most likely related to the effects of the environment of the BN/KL region on the outflow emanating from Source I. The emission of the vibrationally-excited species peaks at the position of Source I. This emission is compact and not resolved at an angular resolution of ~1.5" (~600 AU at a distance of 420 pc). 2-D Gaussian fitting to individual velocity channels locates emission peaks within radii of 100 AU, i.e. they trace the innermost part of the outflow. A narrow spectral profile and spatial distribution of the v=1 J=5-4 line similar to the masing v=1 J=1-0 transition, provide evidence for the most highly rotationally excited (frequency > 200 GHz) SiO maser emission associated with Source I known to date. The maser emission will enable studies of the Source I disk-outflow interface with future ALMA longest baselines.

Radio and infrared interferometry of SiO maser stars provide complementary information on the atmosphere and circumstellar environment at comparable spatial resolution. Here, we present the latest results on the atmospheric structure and the dust condensation region of AGB stars based on our recent infrared spectro-interferometric observations, which represent the environment of SiO masers. We discuss, as an example, new results from simultaneous VLTI and VLBA observations of the Mira variable AGB star R Cnc, including VLTI near- and mid-infrared interferometry, as well as VLBA observations of the SiO maser emission toward this source. We present preliminary results from a monitoring campaign of high-frequency SiO maser emission toward evolved stars obtained with the APEX telescope, which also serves as a precursor of ALMA images of the SiO emitting region. We speculate that large-scale long-period chaotic motion in the extended molecular atmosphere may be the physical reason for observed deviations from point symmetry of atmospheric molecular layers, and for the observed erratic variability of high-frequency SiO maser emission

We report the outcome of a 3-day workshop on the Hubble constant (H_0) that took place during February 6-8 2012 at the Kavli Institute for Particle Astrophysics and Cosmology, on the campus of Stanford University. The participants met to address the following questions. Are there compelling scientific reasons to obtain more precise and more accurate measurements of H_0 than currently available? If there are, how can we achieve this goal? The answers that emerged from the workshop are (1) better measurements of H_0 provide critical independent constraints on dark energy, spatial curvature of the Universe, neutrino physics, and validity of general relativity, (2) a measurement of H_0 to 1% in both precision and accuracy, supported by rigorous error budgets, is within reach for several methods, and (3) multiple paths to independent determinations of H_0 are needed in order to access and control systematics.

Using the Expanded Very Large Array, we have mapped the vicinity of the Orion BN/KL Hot Core with sub-arcsecond angular resolution in seven metastable inversion transitions of ammonia: (J,K)=(6,6) to (12,12). This emission comes from levels up to 1500 K above the ground state, enabling identification of source(s) responsible for heating the region. We used this multi-transition dataset to produce images of the rotational/kinetic temperature and the column density of ammonia for ortho and para species separately and on a position-by-position basis. We find rotational temperature and column density in the range 160-490 K and (1-4)x10^17 cm^-2, respectively. Our spatially-resolved images show that the highest (column) density and hottest gas is found in a northeast-southwest elongated ridge to the southeast of Source I. We have also measured the ortho-para ratio of ammonia, estimated to vary in the range 0.9-1.6. Enhancement of ortho with respect to para and the offset of hot ammonia emission peaks from known (proto)stellar sources provide evidence that the ammonia molecules have been released from dust grains into the gas-phase through the passage of shocks and not by stellar radiation. We propose that the combined effect of Source I's proper motion and its low-velocity outflow impinging on a pre-existing dense medium is responsible for the excitation of ammonia and the Orion Hot Core. Finally, we found for the first time evidence of a slow (5 km/s) and compact (1000 AU) outflow towards IRc7.

We present an overview on our project to study the extended atmospheres and dust formation zones of Mira stars using coordinated observations with the Very Large Telescope Interferometer (VLTI), the Very Long Baseline Array (VLBA), and the Atacama Pathfinder Experiment (APEX). The data are interpreted using an approach of combining recent dynamic model atmospheres with a radiative transfer model of the dust shell, and combining the resulting model structure with a maser propagation model.

An astrophysical MASER (Microwave Amplification by Stimulated Emission of Radiation) is a source of stimulated spectral line emission. Maser emission is observed from the circumstellar envelopes of evolved stars, molecular clouds/star-forming regions, active galactic nuclei, supernova remnants, comets, and the Saturnian moons. It arises from molecules such as water (H2O), hydroxyl radicals (OH), methanol (CH3OH), formaldehyde (CH2O), silicon monoxide (SiO), ammonia (NH3), silicon sulphide (SiS), hydrogen cyanide (HCN), and from atomic hydrogen recombination lines. Masers are compact, of high brightness temperature, and often display narrow spectral line widths, polarized emission and variability. Free electron-cyclotron astrophysical masers additionally exist.

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.

We present new 7mm continuum observations of Orion BN/KL with the VLA. We resolve the emission from the protostar radio Source I and BN at several epochs. Source I is highly elongated NW-SE, and remarkably stable in flux density, position angle, and overall morphology over nearly a decade. This favors the extended emission component arising from an ionized disk rather than a jet. We have measured the proper motions of Source I and BN for the first time at 43 GHz. We confirm that both sources are moving at high speed (12 and 26 km/s, respectively) approximately in opposite directions, as previously inferred from measurements at lower frequencies. We discuss dynamical scenarios that can explain the large motions of both BN and Source I and the presence of disks around both. Our new measurements support the hypothesis that a close (~50 AU) dynamical interaction occurred around 500 years ago between Source I and BN as proposed by Gomez et al. From the dynamics of encounter we argue that Source I today is likely to be a binary with a total mass on the order of 20 Msun, and that it probably existed as a softer binary before the close encounter. This enables preservation of the original accretion disk, though truncated to its present radius of ~50 AU. N-body numerical simulations show that the dynamical interaction between a binary of 20 Msun total mass (I) and a single star of 10 Msun mass (BN) may lead to the ejection of both and binary hardening. The gravitational energy released in the process would be large enough to power the wide-angle flow traced by H2 and CO emission in the BN/KL nebula. Assuming the proposed dynamical history is correct, the smaller mass for Source I recently estimated from SiO maser dynamics (>7 Msun) by Matthews et al., suggests that non-gravitational forces (e.g. magnetic) must play an important role in the circumstellar gas dynamics.

The Megamaser Cosmology Project (MCP) aims to determine H0 by measuring angular-diameter distances to galaxies in the Hubble flow using observations of water vapor megamasers in the circumnuclear accretion disks of active galaxies. The technique is based only on geometry and determines H0 in one step, independent of standard candles and the extragalactic distance ladder. In Paper I we presented a VLBI map of the maser emission from the Seyfert 2 galaxy UGC 3789. The map reveals an edge-on, sub-parsec disk in Keplerian rotation, analogous to the megamaser disk in NGC 4258. Here we present 3.2 years of monthly GBT observations of the megamaser disk in UGC 3789. We use these observations to measure the centripetal accelerations of both the systemic and high-velocity maser components. The measured accelerations suggest that maser emission lines near the systemic velocity originate on the front side of the accretion disk, primarily from segments of two narrow rings. Adopting a two-ring model for the systemic features, we determine the angular-diameter distance to UGC 3789 to be 49.9 +/- 7.0 Mpc. This is the most accurate geometric distance yet obtained to a galaxy in the Hubble flow. Based on this distance, we determine H0 = 69 +/- 11 km/s/Mpc. We also measure the mass of the central black hole to be 1.09 x 10^7 solar masses +/- 14%. With additional observations the uncertainty in the distance to this galaxy can be reduced to under 10%. Observations of megamaser disks in other galaxies will further reduce the uncertainty in H0 as measured by the MCP.

We present multi-epoch VLBA imaging of the 28SiO v=1 & v=2, J=1-0 maser emission toward the massive YSO Orion Source I. Both SiO transitions were observed simultaneously with an angular resolution of ~0.5 mas (~0.2 AU for d=414 pc). Here we explore the global properties and kinematics of the emission through two 19-epoch animated movies spanning 21 months (2001 March 19 to 2002 December 10). These movies provide the most detailed view to date of the dynamics and temporal evolution of molecular material within ~20-100 AU of a massive (~>8M_sun) YSO. The bulk of the SiO masers surrounding Source I lie in an X-shaped locus; emission in the South/East arms is predominantly blueshifted and in the North and West is predominantly redshifted. In addition, bridges of intermediate-velocity emission connect the red and blue sides of the emission distribution. We have measured proper motions of >1000 maser features and find a combination of radially outward migrations along the four arms and motions tangent to the bridges. We interpret the SiO masers as arising from a wide-angle bipolar wind emanating from a rotating, edge-on disk. The detection of maser features along extended, curved filaments suggests that magnetic fields may play a role in launching and/or shaping the wind. Our observations appear to support a picture in which stars with M ~>8 M_sun form via disk-mediated accretion. However, we cannot rule out that the Source I disk may have been formed or altered following a close encounter. (Abridged).

We have used the Very Large Array (VLA) at 7 mm wavelength to image five rotational transitions (J=1-0) from three SiO isotopologues towards Orion BN/KL: 28SiO v=0,1,2; and 29SiO and 30SiO v=0. For the first time, we have mapped the 29SiO and 30SiO J=1-0 emission, established the maser nature of the emission, and confirmed association with the deeply embedded high-mass young stellar object commonly denoted radio Source I. The 28SiO v=0 maser emission shows a bipolar structure that extends over ~700 AU along a northeast-southwest axis, and we propose that it traces a bipolar outflow driven by Source I. The high-brightness isotopic SiO maser emission imaged with a ~0.2 arcsec resolution has a more compact distribution, generally similar to that of the 28SiO v=1,2 emission, and it probably traces bulk gas flows in a region of diameter <100 AU centered on Source I. On small scales of <10 AU, however, compact 29SiO/30SiO v=0 and 28SiO v=1,2 emission features may be offset from one another in position and line-of-sight velocity. From a radiative transfer analysis based on a large velocity gradient (LVG) pumping model, we derive similar temperatures and densities for the optimum excitation of both 29SiO/30SiO v=0 and 28SiO v=1,2 masers, significantly higher than required for 28SiO v=0 maser excitation. In order to account for the small-scale differences among the isotopologues (v=0) and the main species (v=1,2), follow-up radiative transfer modeling that incorporates non-local line overlap among transitions of all SiO isotopic species may be required.

Joint analysis of Cosmic Microwave Background, Baryon Acoustic Oscillation, and supernova data has enabled precision estimation of cosmological parameters. New programs will push to 1% uncertainty in the dark energy equation of state and tightened constraint on curvature, requiring close attention to systematics. Direct 1% measurement of the Hubble constant (H0) would provide a new constraint. It can be obtained without overlapping systematics directly from recessional velocities and geometric distance estimates for galaxies via the mapping of water maser emission that traces the accretion disks of nuclear black holes. We identify redshifts 0.02<z<0.06 as best for small samples, e.g., 10 widely distributed galaxies, each with 3% distance uncertainty. Knowledge of peculiar radial motion is also required. Mapping requires very long baseline interferometry (VLBI) with the finest angular resolution, sensitivity to individual lines of a few mJy-km/s, and baselines that can detect a complex of ~10 mJy lines (peak) in < 1 min. For 2010-2020, large ground apertures (50-100m diameter) augmenting the VLBA are critical, such as EVLA, GBT, Effelsberg, and the Large Millimeter Telescope, for which we propose a 22 GHz receiver and VLBI instrumentation. A space-VLBI aperture may be required, thus motivating US participation in the Japanese VSOP-2 mission (launch c.2013). This will provide 3-4x longer baselines and ~5x improvement in distance uncertainty. There are now 5 good targets at z>0.02, out of ~100 known masers. A single-dish discovery survey of >10,000 nuclei (>2500 hours on the GBT) would build a sample of tens of potential distance anchors. Beyond 2020, a high-frequency SKA could provide larger maser samples, enabling estimation of H0 from individually less accurate distances, and possibly without the need for peculiar motion corrections.

We describe a combined dynamic atmosphere and maser propagation model of SiO maser emission in Mira variables. This model rectifies many of the defects of an earlier model of this type, particularly in relation to the infra-red (IR) radiation field generated by dust and various wavelength-dependent, optically thick layers. Modelled masers form in rings with radii consistent with those found in VLBI observations and with earlier models. This agreement requires the adoption of a radio photosphere of radius approximately twice that of the stellar photosphere, in agreement with observations. A radio photosphere of this size renders invisible certain maser sites with high amplification at low radii, and conceals high-velocity shocks, which are absent in radio continuum observations. The SiO masers are brightest at an optical phase of 0.1 to 0.25, which is consistent with observed phase-lags. Dust can have both mild and profound effects on the maser emission. Maser rings, a shock and the optically thick layer in the SiO pumping band at 8.13\micron appear to be closely associated in three out of four phase samples.

Using the Expanded Very Large Array, we have conducted a search for 22.2 GHz H2O megamaser emission in the strongly lensed submm galaxy, SMM J16359+6612 at z=2.517. This object is lensed into three components, and after a correction for magnification is applied to its submm-wavelength flux density, it is typical of the bulk of the high-redshift, submm galaxy population responsible for the 850 um extragalactic background (S(850um)~1 mJy). We do not detect any H2O megamaser emission, but the lensing allows us to place an interesting constraint on the luminosity of any megamasers present, L(H2O) < 5305 solar luminosities for an assumed linewidth of 80 km/s. Because the far-infrared luminosity in submm galaxies is mainly powered by star formation, and very luminous H2O megamasers are more commonly associated with quasar activity, it could be that blind searches for H2O megamasers will not be an effective means of determining redshifts for less luminous members of the submm galaxy population.

We have surveyed molecular line emission from Orion BN/KL from 42.3 to 43.6 GHz with the Green Bank Telescope. Sixty-seven lines were identified and ascribed to 13 different molecular species. The spectrum at 7 mm is dominated by SiO, SO2, CH3OCH3, and C2H5CN. Five transitions have been detected from the SiO isotopologues 28SiO, 29SiO, and 30SiO. We report here for the first time the spectra of the 29SiO and 30SiO v=0 J=1-0 emission in Orion BN/KL, and we show that they have double-peaked profiles with velocity extents similar to the main isotopologue. The main motivation for the survey was the search of high-velocity (100-1000 km/s) outflows in the BN/KL region as traced by SiO Doppler components. Some of the unidentified lines in principle could be high-velocity SiO features, but without imaging data their location cannot be established. Wings of emission are present in the v=0 28SiO, 29SiO and 30SiO profiles, and we suggest that the v=0 emission from the three isotopologues might trace a moderately high-velocity (~30-50 km/s) component of the flows around the high-mass protostar Source I in the Orion BN/KL region. We also confirm the 7 mm detection of a complex oxygen-bearing species, acetone (CH3COCH3), which has been recently observed towards the hot core at 3 mm, and we have found further indications of the presence of long cyanopolyynes (HC5N and HC7N) in the quiescent cold gas of the extended ridge.

We report measurements of centripetal accelerations of maser spectral components of NGC 4258 for 51 epochs spanning 1994 to 2004. This is the second paper of a series, in which the goal is determination of a new geometric maser distance to NGC 4258 accurate to possibly ~3%. We measure accelerations using a formal analysis method that involves simultaneous decomposition of maser spectra for all epochs into multiple, Gaussian components. Components are coupled between epochs by linear drifts (accelerations) from their centroid velocities at a reference epoch. For high-velocity emission, accelerations lie in the range -0.7 to +0.7 km/s/yr indicating an origin within 13 degrees of the disk midline (the perpendicular to the line-of-sight to the black hole). Comparison of high-velocity emission projected positions in VLBI images, with those derived from acceleration data, provides evidence that masers trace real gas dynamics. High-velocity emission accelerations do not support a model of trailing shocks associated with spiral arms in the disk. However, we find strengthened evidence for spatial periodicity in high-velocity emission, of wavelength 0.75 mas. This supports suggestions of spiral structure due to density waves in the nuclear accretion disk of an active galaxy. Accelerations of low-velocity (systemic) emission lie in the range 7.7 to 8.9 km/s/yr, consistent with emission originating from a concavity where the thin, warped disk is tangent to the line-of-sight. A trend in accelerations of low-velocity emission as a function of Doppler velocity may be associated with disk geometry and orientation, or with the presence of spiral structure.

Despite theoretical predictions of the existence of many submillimeter masers, and some pioneering observational discoveries over the past few decades, these lines have remained relatively unstudied due to (i) challenges associated with observing at shorter wavelength; and, (ii) lack of possibility of high (< 14'' at 345 GHz) angular resolution observations. With the advent of the SMA, the first submillimeter imaging array capable of sub-arcsecond resolution, APEX, and the promise of ALMA, opportunities are opening for performing new science with millimeter/submillmeter masers. In this talk, I will review recent work in the field - including extragalactic water millimeter masers, hydrogen recombination masers, submillimeter masers in star-forming regions, and in the envelopes of evolved stars - and discuss prospects for the future.

A comprehensive picture of high-mass star formation has remained elusive, in part because examples of high-mass YSOs tend to be relatively distant, deeply embedded, and confused with other emission sources. These factors have impeded dynamical investigations within tens of AU of high-mass YSOs--scales that are critical for probing the interfaces where outflows from accretion disks are launched and collimated. Using observations of SiO masers obtained with the VLA and the VLBA, the KaLYPSO project is overcoming these limitations by mapping the structure and dynamical/temporal evolution of the material 10-1000 AU from the nearest high-mass YSO: Radio Source I in the Orion BN/KL region. Our data include ~40 epochs of VLBA observations over a several-year period, allowing us to track the proper motions of individual SiO maser spots and to monitor changes in the physical conditions of the emitting material with time. Ultimately these data will provide 3-D maps of the outflow structure over approximately 30% of the outflow crossing time. Here we summarize recent results from the KaLYPSO project, including evidence that high-mass star formation is occurring via disk-mediated accretion.

We report a three year, 18 epoch, VLBI monitoring study of H2O masers in the sub-parsec, warped, accretion disk within the NGC4258 AGN. Our immediate goals are to trace the geometry of the underlying disk, track rotation via measurement of proper motion, and ascertain the radii of masers for which centripetal acceleration may be measured separately. The monitoring includes ~ 4 times as many epochs, ~ 3 times denser sampling, and tighter control over sources of systematic error than earlier VLBI investigations. Coverage of a ~ 2400 km/s bandwidth has also enabled mapping of molecular material ~ 30% closer to the black hole than accomplished previously, which will strengthen geometric and dynamical disk models. Through repeated observation we have also measured for the first time a 5 microarsecond (1 sigma) thickness of the maser medium. Assuming this corresponds to the thickness of the accretion disk, hydrostatic equilibrium requires a disk plane temperature of ~ 600 K. Our long-term goal is a geometric distance to NGC4258 that is accurate to ~ 3%, a ~ 2 times improvement over the current best estimate. A geometric estimate of distance can be compared to distances obtained from analysis of Cepheid light curves, with the intent to recalibrate the extragalactic distance scale with reduced systematic uncertainties. This is the first paper in a series. We present here VLBI observations, data reduction, and temporal and spatial characteristics of the maser emission. Later papers will report estimation of orbital acceleration and proper motion, modeling of disk 3-D geometry and dynamics, and estimation of a "maser distance." Estimation of a "Cepheid distance" is presented in a parallel paper series.

We report the first detection of an extragalactic millimeter wavelength H2O maser at 183 GHz towards NGC 3079 using the Submillimeter Array (SMA), and a tentative submillimeter wave detection of the 439 GHz maser towards the same source using the James Clerk Maxwell Telescope (JCMT). These H2O transitions are known to exhibit maser emission in star-forming regions and evolved stars. NGC 3079 is a well-studied nuclear H2O maser source at 22 GHz with a time-variable peak flux density in the range 3 -- 12 Jy. The 183 GHz H2O maser emission, with peak flux density $\sim$0.5 Jy (7$\sigma$ detection), also originates from the nuclear region of NGC 3079 and is spatially coincident with the dust continuum peak at 193 GHz (53 mJy integrated). Peak emission at both 183 and 439 GHz occurs in the same range of velocity as that covered by the 22 GHz spectrum. We estimate the gas to dust ratio of the nucleus of NGC 3079 to be $\approx$150, comparable to the Galactic value of 160. Discovery of maser emission in an active galactic nucleus beyond the long-known 22 GHz transition opens the possibility of future position-resolved radiative transfer modeling of accretion disks and outflows $<1$ pc from massive black holes.

We report on our ongoing, high-resolution study of H2O masers in the innermost parsec of NGC 4258. Over thirty epochs of VLBA and VLA data, taken over six years, are being used to monitor the velocities, accelerations, positions and proper motions of water masers rotating in a warped, Keplerian disk about a supermassive central object. Our extensive monitoring results in an improved accuracy in the distance determination to this galaxy (i) via a reduction in experimental random errors due to a longer time-baseline than in previous work; and (ii) via a better modeling of sources of systematic error, such as disk eccentricity. These data can therefore yield an extremely accurate geometric distance to NGC 4258.

A stellar hydrodynamic pulsation model has been combined with a SiO maser model in an attempt to calculate the temporal variability of SiO maser emission in the circumstellar envelope (CE) of a model AGB star. This study investigates whether the variations in local physical conditions brought about by shocks are the predominant contributing factor to SiO maser variability because, in this work, the radiative part of the pump is constant. We find that some aspects of the variability are not consistent with a pump provided by shock-enhanced collisions alone. In these simulations, gas parcels of relatively enhanced SiO abundance are distributed in a model CE by a Monte Carlo method, at a single epoch of the stellar cycle. From this epoch on, Lagrangian motions of individual parcels are calculated according to the velocity fields encountered in the model CE during the stellar pulsation cycle. The potentially masing gas parcels therefore experience different densities and temperatures, and have varying line-of-sight velocity gradients throughout the stellar cycle, which may or may not be suitable to produce maser emission. At each epoch (separated by 16.6 days), emission lines from the parcels are combined to produce synthetic spectra and VLBI-type images. We report here the results for v=1, J=1-0 (43-GHz) and J=2-1 (86-GHz) masers.

We report SIMBA 1.2 mm dust continuum observations of the environments of eight methanol maser sources, all discovered during spatially fully-sampled, untargeted surveys of the galactic plane. We summarise our search for possible associations of the masers with IR sources (IRAS and MSX) and find that it is not always possible to make definite associations. A preliminary characterisation of the IR sources found in the maser neighbourhood is given according their position in the [60-25] -- [25-12] colour-colour diagram.

Building on the detection of J=7-6 SiO maser emission in both the v=1 and v=2 vibrational states towards the symbiotic Mira, R Aquarii, we have used the James Clerk Maxwell Telescope to study the changes in the SiO maser features from R Aqr over a stellar pulsational period. The observations, complemented by contemporaneous data taken at 86 GHz, represent a test of the popular thermal-instability clump models of SiO masers. The `clump' model of SiO maser emission considers the SiO masers to be discrete emitting regions which differ from their surroundings in the values of one or more physical variables (SiO abundance, for example). We find that our observational data are consistent with a clump model in which the appearance of maser emission in the J=7-6 transitions coincides with an outward-moving shock impinging on the inner edge of the maser zone.