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.

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.

We discuss the science drivers for ALMA Band 2 which spans the frequency range from 67 to 90 GHz. The key science in this frequency range are the study of the deuterated molecules in cold, dense, quiescent gas and the study of redshifted emission from galaxies in CO and other species. However, Band 2 has a range of other applications which are also presented. The science enabled by a single receiver system which would combine ALMA Bands 2 and 3 covering the frequency range 67 to 116 GHz, as well as the possible doubling of the IF bandwidth of ALMA to 16 GHz, are also considered.

We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5\mu$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1.

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 role of massive stars is central to an understanding of galactic ecology. It is important to establish the details of how massive stars provide radiative, chemical, and mechanical feedback in galaxies. Central to these issues is an understanding of the evolution of massive stars, and the critical role of mass loss via strongly structured winds and stellar binarity. Ultimately, massive stellar clusters shape the structure and energetics of galaxies. We aim to conduct high-resolution, deep field mapping at 21cm of the core of the massive Cygnus OB2 association and to characterise the properties of the massive stars and colliding winds at this waveband. We used seven stations of the e-MERLIN radio facility, with its upgraded bandwidth and enhanced sensitivity to conduct a 21cm census of Cygnus OB2. Based on 42 hours of observations, seven overlapping pointings were employed over multiple epochs during 2014 resulting in 1 sigma sensitivities down to ~21microJy and a resolution of ~180mas. A total of 61 sources are detected at 21cm over a ~0.48deg x 0.48deg region centred on the heart of the Cyg OB2 association. Of these 61 sources, 33 are detected for the first time. We detect a number of previously identified sources including four massive stellar binary systems, two YSOs, and several known X-ray and radio sources. We also detect the LBV candidate (possible binary system) and blue hypergiant (BHG) star of Cyg OB2 #12. The 21cm observations secured in the COBRaS Legacy project provide data to constrain conditions in the outer wind regions of massive stars; determine the non-thermal properties of massive interacting binaries; examine evidence for transient sources, including those associated with young stellar objects; and provide unidentified sources that merit follow-up observations. The 21cm data are of lasting value and will serve in combination with other key surveys of Cyg OB2.

We present and analyse SO and SO$_2$, recently observed with high angular resolution and sensitivity in a spectral line survey with ALMA, for two oxygen-rich AGB stars: the low mass-loss rate R Dor and high mass-loss rate IK Tau. We analyse 8 lines of SO detected towards both stars, 78 lines of SO$_2$detected towards R Dor and 52 lines of SO$_2$ detected towards IK Tau. We detect several lines of $^{34}$SO, $^{33}$SO and $^{34}$SO$_2$ towards both stars, and tentatively S$^{18}$O towards R Dor, and hence derive isotopic ratios for these species. The spatially resolved observations show us that the two sulphur oxides are co-located towards R Dor and trace out the same wind structures in the circumstellar envelope (CSE). Much of the emission is well reproduced with a Gaussian abundance distribution spatially centred on the star. Emission from the higher energy levels of SO and SO$_2$ towards R Dor provide evidence in support of a rotating inner region of gas identified in earlier work. The new observations allow us to refine the abundance distribution of SO in IK Tau derived from prior observations with single antennas, and confirm the distribution is shell-like with the peak in the fractional abundance not centred on the star. The confirmation of different types of SO abundance distributions will help fine-tune chemical models and allows for an additional method to discriminate between low and high mass-loss rates for oxygen-rich AGB stars.

LUX-ZEPLIN (LZ) is a next generation dark matter direct detection experiment that will operate 4850 feet underground at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. Using a two-phase xenon detector with an active mass of 7 tonnes, LZ will search primarily for low-energy interactions with Weakly Interacting Massive Particles (WIMPs), which are hypothesized to make up the dark matter in our galactic halo. In this paper, the projected WIMP sensitivity of LZ is presented based on the latest background estimates and simulations of the detector. For a 1000 live day run using a 5.6 tonne fiducial mass, LZ is projected to exclude at 90% confidence level spin-independent WIMP-nucleon cross sections above $1.6 \times 10^{-48}$ cm$^{2}$ for a 40 $\mathrm{GeV}/c^{2}$ mass WIMP. Additionally, a $5\sigma$ discovery potential is projected reaching cross sections below the existing and projected exclusion limits of similar experiments that are currently operating. For spin-dependent WIMP-neutron(-proton) scattering, a sensitivity of $2.7 \times 10^{-43}$ cm$^{2}$ ($8.1 \times 10^{-42}$ cm$^{2}$) for a 40 $\mathrm{GeV}/c^{2}$ mass WIMP is expected. With underground installation well underway, LZ is on track for commissioning at SURF in 2020.

We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements.

We report the first detection of CO J = 3-2 around a truly metal-poor evolved star. RU Vulpeculae is modelled to have Teff ~ 3620 K, L ~ 3128 +/- 516 L_Sun, log(g) = 0.0 ~ 0.2 dex and [Fe/H] = -1.3 to -1.0 dex, and is modelled to have recently undergone a thermal pulse. Its infrared flux has approximately doubled over 35 years. ALMA observations show the 3-2 line is narrow (half-width ~1.8-3.5 km/s). The 2-1 line is much weaker: it is not confidently detected. Spectral-energy-distribution fitting indicates very little circumstellar absorption, despite its substantial mid-infrared emission. A VISIR mid-infrared spectrum shows features typical of previously observed metal-poor stars, dominated by a substantial infrared excess but with weak silicate and (possibly) Al2O3 emission. A lack of resolved emission, combined with weak 2!1 emission, indicates the dense circumstellar material is truncated at large radii. We suggest that rapid dust condensation is occurring, but with an aspherical geometry (e.g., a disc or clumps) that does not obscure the star. We compare with T UMi, a similar star which is currently losing its dust.

We present high angular resolution (~80 mas) ALMA continuum images of the SN 1987A system, together with CO $J$=2 $\!\rightarrow\!$ 1, $J$=6 $\!\rightarrow\!$ 5, and SiO $J$=5 $\!\rightarrow\!$ 4 to $J$=7 $\!\rightarrow\!$ 6 images, which clearly resolve the ejecta (dust continuum and molecules) and ring (synchrotron continuum) components. Dust in the ejecta is asymmetric and clumpy, and overall the dust fills the spatial void seen in H$\alpha$ images, filling that region with material from heavier elements. The dust clumps generally fill the space where CO $J$=6 $\!\rightarrow\!$ 5 is fainter, tentatively indicating that these dust clumps and CO are locationally and chemically linked. In these regions, carbonaceous dust grains might have formed after dissociation of CO. The dust grains would have cooled by radiation, and subsequent collisions of grains with gas would also cool the gas, suppressing the CO $J$=6 $\!\rightarrow\!$ 5 intensity. The data show a dust peak spatially coincident with the molecular hole seen in previous ALMA CO $J$=2 $\!\rightarrow\!$ 1 and SiO $J$=5 $\!\rightarrow\!$ 4 images. That dust peak, combined with CO and SiO line spectra, suggests that the dust and gas could be at higher temperatures than the surrounding material, though higher density cannot be totally excluded. One of the possibilities is that a compact source provides additional heat at that location. Fits to the far-infrared--millimeter spectral energy distribution give ejecta dust temperatures of 18--23K. We revise the ejecta dust mass to $\mathrm{M_{dust}} = 0.2-0.4$M$_\odot$ for carbon or silicate grains, or a maximum of $<0.7$M$_\odot$ for a mixture of grain species, using the predicted nucleosynthesis yields as an upper limit.

The Spectral Imaging of the Coronal Environment (SPICE) instrument is a high-resolution imaging spectrometer operating at extreme ultraviolet (EUV) wavelengths. In this paper, we present the concept, design, and pre-launch performance of this facility instrument on the ESA/NASA Solar Orbiter mission. The goal of this paper is to give prospective users a better understanding of the possible types of observations, the data acquisition, and the sources that contribute to the instrument's signal. The paper discusses the science objectives, with a focus on the SPICE-specific aspects, before presenting the instrument's design, including optical, mechanical, thermal, and electronics aspects. This is followed by a characterisation and calibration of the instrument's performance. The paper concludes with descriptions of the operations concept and data processing. The performance measurements of the various instrument parameters meet the requirements derived from the mission's science objectives. The SPICE instrument is ready to perform measurements that will provide vital contributions to the scientific success of the Solar Orbiter mission.

Red supergiant stars are surrounded by a gaseous and dusty circumstellar environment created by their mass loss which spreads heavy elements into the interstellar medium. The structure and the dynamics of this envelope are crucial to understand the processes driving the red supergiant mass loss and the shaping of the pre-supernova ejecta. We have observed the emission from the CO $J = 2-1$ line from the red supergiant star $\mu$~Cep with the NOEMA interferometer. In the line the synthesized beam was $0.92 \times 0.72$~arcsec ($590 \times 462$~au at 641~pc). The continuum map shows only the unresolved contribution of the free-free emission of the star chromosphere. The continuum-subtracted channel maps reveal a very inhomogeneous and clumpy circumstellar environment. In particular, we detected a bright CO clump, as bright as the central source in the line, at 1.80~arcsec south-west from the star, in the blue channel maps. After a deprojection of the radial velocity assuming two different constant wind velocities, the observations were modelled using the 3D radiative transfer code \textsc{lime} to derive the characteristics of the different structures. We determine that the gaseous clumps observed around $\mu$~Cep are responsible for a mass loss rate of $(4.9 \pm 1.0) \times 10^{-7}~{\rm M}_\odot\,{\rm yr}^{-1}$, in addition to a spatially unresolved wind component with an estimated mass-loss rate of $2.0 \times 10^{-6}~{\rm M}_\odot\,{\rm yr}^{-1}$. Therefore, the clumps have a significant role in $\mu$~Cep's mass loss ($\ge 25 \%$). We cannot exclude that the unresolved central outflow may be made of smaller unresolved clumps.

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.

In this paper we present new ALMA observations towards the proto-planet hosting transitional disc of Herbig Ae/Be star HD 100546. This includes resolved 1.3 mm continuum, $^{13}$CO and the first detection of C$^{18}$O in this disc, which displays azimuthal asymmetry in regions spatially coincident with structures previously identified in HST images related to spiral arms. The lower limit on the mass of the dust disc is calculated to be 9.6x10$^{-4}$M$_\odot$. A firm lower-limit on the total gas mass calculated from optically thin, mid-plane tracing C$^{18}$O (2-1) emission is 0.018M$_\odot$ assuming ISM abundances. These mass estimates provide an estimate of gas-to-dust ratio in the disc of 19, the ratio will increase if C$^{18}$O is relatively under-abundant in the disc compared to CO and H2. Through deprojection and azimuthal averaging of the image plane we detect 1.3 mm continuum emission out to 290+/-10 au,$^{13}$CO to 390+/-10 au and C$^{18}$O to 300+/-10au. We measure a radially increasing millimetre spectral index between wavelengths of 867$\mu$m and 1.3 mm, which shows that grain sizes increase towards the star, with solid particles growing to cm scales in the inner disc.

ALMA products are stored in the Science Archive in the form of FITS images. It is a common idea that the FITS image headers should collect in their keywords all the information that an archive User might want to search for in order to quickly select, compare, or discard datasets. With this perspective in mind, we first present a short description of the current status of the ALMA FITS archive and images. We realized that at the moment most of the parameters that could be useful for a general User are still missing in the archived data. We then provide a CASA task generating the image header keywords that we suggest to be relevant for the scientific exploitation of the ALMA archival data. The proposed tool could be also applied to several types of interferometer data and part of it is implemented in a web interface. An example of the scientific application of the keywords is also discussed.

We report the first detection of circumstellar CO in a globular cluster. Observations with ALMA have detected the CO J=3-2 and SiO v=1 J=8-7 transitions at 345 and 344 GHz, respectively, around V3 in 47 Tucanae (NGC 104; [Fe/H] = -0.72 dex), a star on the asymptotic giant branch. The CO line is detected at 7 sigma at a rest velocity v_LSR = -40.6 km/s and expansion velocity of 3.2 +/- ~0.4 km/s. The brighter, asymmetric SiO line may indicate a circumstellar maser. The stellar wind is slow compared to similar Galactic stars, but the dust opacity remains similar to Galactic comparisons. We suggest that the mass-loss rate is set by the levitation of material into the circumstellar environment by pulsations, but that the terminal wind-expansion velocity is determined by radiation pressure on the dust: a pulsation-enhanced dust-driven wind. We suggest the metal-poor nature of the star decreases the grain size, slowing the wind and increasing its density and opacity. Metallic alloys at high altitudes above the photosphere could also provide an opacity increase. The CO line is weaker than expected from Galactic AGB stars, but its strength confirms a model that includes CO dissociation by the strong interstellar radiation field present inside globular clusters.

We aim to determine the distributions of molecular SiS and CS in the circumstellar envelopes of oxygen-rich asymptotic giant branch stars and how these distributions differ between stars that lose mass at different rates. In this study we analyse ALMA observations of SiS and CS emission lines for three oxygen-rich galactic AGB stars: IK Tau, with a moderately high mass-loss rate of $5\times10^{-6}$M$_\odot$ yr$^{-1}$, and W Hya and R Dor with low mass loss rates of $\sim1\times10^{-7}$M$_\odot$ yr$^{-1}$. These molecules are usually more abundant in carbon stars but the high sensitivity of ALMA allows us to detect their faint emission in the low mass-loss rate AGB stars. The high spatial resolution of ALMA also allows us to precisely determine the spatial distribution of these molecules in the circumstellar envelopes. We run radiative transfer models to calculate the molecular abundances and abundance distributions for each star. We find a spread of peak SiS abundances with $\sim10^{-8}$ for R Dor, $\sim10^{-7}$ for W Hya, and $\sim3\times10^{-6}$ for IK Tau relative to H$_2$. We find lower peak CS abundances of $\sim7\times10^{-9}$ for R Dor, $\sim7\times10^{-8}$ for W Hya and $\sim4\times10^{-7}$ for IK Tau, with some stratifications in the abundance distributions. For IK Tau we also calculate abundances for the detected isotopologues: C$^{34}$S, $^{29}$SiS, $^{30}$SiS, Si$^{33}$S, Si$^{34}$S, $^{29}$Si$^{34}$S, and $^{30}$Si$^{34}$S. Overall the isotopic ratios we derive for IK Tau suggest a lower metallicity than solar.

We present a comparison of Doppler-shifted H-alpha line emission observed by the Global Jet Watch from freshly-launched jet ejecta at the nucleus of the Galactic microquasar SS433 with subsequent ALMA imaging at mm-wavelengths of the same jet ejecta. There is a remarkable similarity between the transversely-resolved synchrotron emission and the prediction of the jet trace from optical spectroscopy: this is an a priori prediction not an a posteriori fit, confirming the ballistic nature of the jet propagation. The mm-wavelength of the ALMA polarimetry is sufficiently short that the Faraday rotation is negligible and therefore that the observed E-vector directions are accurately orthogonal to the projected local magnetic field. Close to the nucleus the B-field vectors are perpendicular to the direction of propagation. Further out from the nucleus, the B-field vectors that are coincident with the jet instead become parallel to the ridge line; this occurs at a distance where the jet bolides are expected to expand into one another. X-ray variability has also been observed at this location; this has a natural explanation if shocks from the expanding and colliding bolides cause particle acceleration. In regions distinctly separate from the jet ridge line, the fractional polarisation approaches the theoretical maximum for synchrotron emission.

High-resolution interferometric observations of the circumstellar environments of AGB stars show a variety of morphologies. Guided by the unusual carbon monoxide line profile of the AGB star EP Aquarii, we have observed its circumstellar environment with ALMA band 6 in cycle 4. We describe the morphological complexity of the CO, SiO and SO2 molecular emission. The CO emission exhibits the characteristics of a bi-conical wind with a bright nearly face-on spiral feature around the systemic velocity. This is the first convincing detection of a spiral morphology in an O-rich wind. Based on the offsets of the centers of the two bi-conical wind hemispheres, we deduce the position angle of the inclination axis to be ~150 deg measured counterclockwise from north. Based on the velocity width of the spiral signature we estimate the inclination angle of the system to be between 4 deg and 18 deg. The central emission zone exhibits a morphology that resembles simulations modelling the spiral-inducing wind-Roche-lobe-overflow mechanism. Though the spiral may be a (companion-induced) density enhancement in the stellar outflow, the extremely narrow width of the spiral signature in velocity space suggests that it may be a hydrodynamical perturbation in a face-on differentially rotating disk. The SiO emission does not show the spiral, but exhibits a local emission void approximately 0.5'' west of the continuum brightness peak. We hypothesise that this may be a local environment caused by the presence of a (stellar) companion with a mass of at most 0.1 Msol, based on its non-detection in the continuum...

We use MasterCode to perform a frequentist analysis of the constraints on a phenomenological MSSM model with 11 parameters, the pMSSM11, including constraints from ~ 36/fb of LHC data at 13 TeV and PICO, XENON1T and PandaX-II searches for dark matter scattering, as well as previous accelerator and astrophysical measurements, presenting fits both with and without the $(g-2)_{\mu}$ constraint. The pMSSM11 is specified by the following parameters: 3 gaugino masses $M_{1,2,3}$, a common mass for the first-and second-generation squarks $m_{\tilde{q}}$ and a distinct third-generation squark mass $m_{\tilde{q}_3}$, a common mass for the first-and second-generation sleptons $m_{\tilde l}$ and a distinct third-generation slepton mass $m_{\tilde \tau}$, a common trilinear mixing parameter $A$, the Higgs mixing parameter $\mu$, the pseudoscalar Higgs mass $M_A$ and $\tan\beta$. In the fit including $(g-2)_{\mu}$, a Bino-like $\tilde\chi^0_1$ is preferred, whereas a Higgsino-like $\tilde \chi^0_1$ is favoured when the $(g-2)_{\mu}$ constraint is dropped. We identify the mechanisms that operate in different regions of the pMSSM11 parameter space to bring the relic density of the lightest neutralino, $\tilde\chi^0_1$, into the range indicated by cosmological data. In the fit including $(g-2)_{\mu}$, coannihilations with $\tilde \chi^0_2$ and the Wino-like $\tilde\chi^{\pm}_1$ or with nearly-degenerate first- and second-generation sleptons are favoured, whereas coannihilations with the $\tilde \chi^0_2$ and the Higgsino-like $\tilde\chi^{\pm}_1$ or with first- and second-generation squarks may be important when the $(g-2)_{\mu}$ constraint is dropped. Prospects remain for discovering strongly-interacting sparticles at the LHC as well as for discovering electroweakly-interacting sparticles at a future linear $e^+ e^-$ collider such as the ILC or CLIC.

Massive stars play an important role in both cluster and galactic evolution and the rate at which they lose mass is a key driver of both their own evolution and their interaction with the environment up to and including their SNe explosions. Young massive clusters provide an ideal opportunity to study a co-eval population of massive stars. We performed 3mm continuum observations with the Atacama Large Millimetre/submillimetre Array of the Galactic cluster Westerlund 1, to study the constituent massive stars and determine mass-loss rates for the diverse post-main sequence population. We detected emission from 50 stars in Westerlund 1, comprising all 21 Wolf-Rayets within the field of view, eight cool and 21 OB super-/hypergiants. Emission nebulae were associated with a number of the cool hypergiants while, unexpectedly, a number of hot stars also appear spatially resolved. We measured the mass-loss rates for a unique population of massive post-main sequence stars at every stage of evolution, confirming a significant increase as stars transition from OB supergiant to WR states. The range of spectral types exhibited provides a critical test of radiatively driven wind theory and the reality of the bi-stability jump. The extreme mass-loss rate inferred for the interacting binary Wd1-9 in comparison to other cluster members confirmed the key role binarity plays in massive stellar evolution. The presence of compact nebulae around a number of OB and WR stars is unexpected; by analogy to the cool super-/hypergiants we attribute this to confinement and sculpting of the stellar wind via interaction with the intra-cluster medium/wind. Given the morphology of core collapse SNe depend on the nature of the pre-explosion circumstellar environment, if this hypothesis is correct then the properties of the explosion depend not just on the progenitor, but also the environment in which it is located.

We describe a likelihood analysis using MasterCode of variants of the MSSM in which the soft supersymmetry-breaking parameters are assumed to have universal values at some scale $M_{in}$ below the supersymmetric grand unification scale $M_{GUT}$, as can occur in mirage mediation and other models. In addition to $M_{in}$, such `sub-GUT' models have the 4 parameters of the CMSSM, namely a common gaugino mass $m_{1/2}$, a common soft supersymmetry-breaking scalar mass $m_0$, a common trilinear mixing parameter $A$ and the ratio of MSSM Higgs vevs $\tan\beta$, assuming that the Higgs mixing parameter $\mu > 0$. We take into account constraints on strongly- and electroweakly-interacting sparticles from $\sim 36$/fb of LHC data at 13 TeV and the LUX and 2017 PICO, XENON1T and PandaX-II searches for dark matter scattering, in addition to the previous LHC and dark matter constraints as well as full sets of flavour and electroweak constraints. We find a preference for $M_{in} \sim 10^5$ to $10^9$ GeV, with $M_{in} \sim M_{GUT}$ disfavoured by $\Delta \chi^2 \sim 3$ due to the ${\rm BR}(B_{s, d} \to \mu^+\mu^-)$ constraint. The lower limits on strongly-interacting sparticles are largely determined by LHC searches, and similar to those in the CMSSM. We find a preference for the LSP to be a Bino or Higgsino with $\tilde{\chi^0_1} \sim 1$ TeV, with annihilation via heavy Higgs bosons $H/A$ and stop coannihilation, or chargino coannihilation, bringing the cold dark matter density into the cosmological range. We find that spin-independent dark matter scattering is likely to be within reach of the planned LUX-Zeplin and XENONnT experiments. We probe the impact of the $(g-2)_\mu$ constraint, finding similar results whether or not it is included.

The circumstellar environment of the AGB star L$_{\rm 2}$ Puppis was observed with ALMA in cycle 3, with a resolution of $15 \times 18 \rm\ mas$. The molecular emission shows a differentially rotating disk, inclined to a nearly edge-on position. In the first paper in this series (paper I) the molecular emission was analysed to accurately deduce the motion of the gas in the equatorial regions of the disk. In this work we model the optically thick $^{\rm 12}$CO $J=$3$-$2 and the optically thin $^{\rm 13}$CO $J=$3$-$2 rotational transition to constrain the physical conditions in the disk. To realise this effort we make use of the 3D NLTE radiative transfer code {\tt LIME}. The temperature structure and velocity structure show a high degree of complexity, both radially and vertically. The radial H$_{\rm 2}$ density profile in the disk plane is characterised by a power law with a slope of $-3.1$. We find a $^{\rm 12}$CO over $^{\rm 13}$CO abundance ratio of 10 inside the disk. Finally, estimations of the angular momentum in the disk surpass the expected available angular momentum of the star, strongly supporting the indirect detection of a compact binary companion reported in paper I. We estimate the mass of the companion to be around 1 Jupiter mass.

Using ALMA, we observed the stellar wind of two oxygen-rich Asymptotic Giant Branch (AGB) stars, IK Tau and R Dor, between 335 and 362 GHz. One aim was to detect metal oxides and metal hydroxides (AlO, AlOH, FeO, MgO, MgOH), some of which are thought to be direct precursors of dust nucleation and growth. We report on the potential first detection of FeO (v=0, Omega=4, J=11-10) in RDor (mass-loss rate, Mdot, ~1e-7 Msun/yr). The presence of FeO in IK Tau (Mdot~5e-6 Msun/yr) cannot be confirmed due to a blend with 29SiS, a molecule that is absent in R Dor. The detection of AlO in R Dor and of AlOH in IK Tau was reported earlier by Decin et al. (2017). All other metal oxides and hydroxides, as well as MgS, remain undetected. We derive a column density N(FeO) of 1.1+/-0.9e15 cm^{-2} in R Dor, or a fractional abundance [FeO/H]~1.5e-8 accounting for non-LTE effects. The derived fractional abundance [FeO/H] is a factor ~20 larger than conventional gas-phase chemical kinetic predictions. This discrepancy may be partly accounted for by the role of vibrationally excited OH in oxidizing Fe, or may be evidence for other currently unrecognised chemical pathways producing FeO. Assuming a constant fractional abundance w.r.t. H_2, the upper limits for the other metals are [MgO/H_2] <5.5e-10 (R Dor) and <7e-11 (IK Tau), [MgOH/H_2] <9e-9 (R Dor) and <1e-9 (IK Tau), [CaO/H_2] <2.5e-9 (R Dor) and <1e-10 (IK Tau), [CaOH/H_2] <6.5e-9 (R Dor) and <9e-10 (IK Tau), and [MgS/H_$] <4.5e-10 (R Dor) and <6e-11 (IK Tau). The retrieved upper limit abundances for these latter molecules are in accord with the chemical model predictions.

A spectral line and imaging survey of the low mass-loss rate AGB star R Dor (Mdot ~ 1e-7 Msun/yr) and the high mass-loss rate AGB star IK Tau (Mdot ~5e-6 Msun/yr) was made with ALMA between 335 and 362 GHz at a spatial resolution of ~150 mas, corresponding to the locus of the main dust formation region of both targets. Some 200 spectral features from 15 molecules (and their isotopologues) were observed, including rotational lines in both the ground and vibrationally excited states. Detected species include the gaseous precursors of dust grains such as SiO, AlO, AlOH, TiO, and TiO2. We present a spectral atlas for both stars and the parameters of all detected spectral features. A clear dichotomy for the sulphur chemistry is seen: while CS, SiS, SO, and SO2 are abundantly present in IK Tau, only SO and SO2 are detected in R Dor. Also other species such as NaCl, NS, AlO, and AlOH display a completely different behaviour. From some selected species, the minor isotopologues can be used to assess the isotopic ratios. The channel maps of many species prove that both large and small-scale inhomogeneities persist in the inner wind of both stars in the form of blobs, arcs, and/or a disk. The high sensitivity of ALMA allows us to spot the impact of these correlated density structures in the spectral line profiles. The spectral lines often display a half width at zero intensity much larger than expected from the terminal velocity, v_inf, previously derived for both objects (36 km/s versus v_inf ~17.7 km/s for IK Tau and 23 km/s versus v_inf ~5.5 km/s for R Dor). Both a more complex 3D morphology and a more forceful wind acceleration of the (underlying) isotropic wind can explain this trend. The formation of fractal grains in the region beyond ~400 mas can potentially account for the latter scenario. From the continuum map, we deduce a dust mass of ~3.7e-7 Msun for IK Tau and ~2e-8 Msun for R Dor.

We observed Betelgeuse using ALMA's extended configuration in band 7 (f~340 GHz, {\lambda}~0.88 mm), resulting in a very high angular resolution of 18 mas. Using a solid body rotation model of the 28SiO(v=2,J=8-7) line emission, we show that the supergiant is rotating with a projected equatorial velocity of v_eq sin i = 5.47 +/- 0.25 km/s at the equivalent continuum angular radius R_star = 29.50 +/- 0.14 mas. This corresponds to an angular rotation velocity of {\omega} sin i = (5.6 +/- 1.3) x 10^(-9) rad/s. The position angle of its north pole is PA = 48.0 +/- 3.5{\deg}. The rotation period of Betelgeuse is estimated to P/sin i = 36 +/- 8 years. The combination of our velocity measurement with previous observations in the ultraviolet shows that the chromosphere is co-rotating with the star up to a radius of ~10 au (45 mas or 1.5x the ALMA continuum radius). The coincidence of the position angle of the polar axis of Betelgeuse with that of the major ALMA continuum hot spot, a molecular plume, and a partial dust shell (from previous observations) suggests that focused mass loss is currently taking place in the polar region of the star. We propose that this hot spot corresponds to the location of a particularly strong "rogue" convection cell, which emits a focused molecular plume that subsequently condenses into dust at a few stellar radii. Rogue convection cells therefore appear to be an important factor shaping the anisotropic mass loss of red supergiants.

Variable OH/IR stars are Asymptotic Giant Branch (AGB) stars with an optically thick circumstellar envelope that emit strong OH 1612 MHz emission. They are commonly observed throughout the Galaxy but also in the LMC and SMC. Hence, the precise inference of the distances of these stars will ultimately result in better constraints on their mass range in different metallicity environments. Through a multi-year long-term monitoring program at the Nancay Radio telescope (NRT) and a complementary high-sensitivity mapping campaign at the eMERLIN and JVLA to measure precisely the angular diameter of the envelopes, we have been re-exploring distance determination through the phase-lag method for a sample of stars, in order to refine the poorly-constrained distances of some and infer the currently unknown distances of others. We present here an update of this project.

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.

Proxima Centauri, the star closest to our Sun, is known to host at least one terrestrial planet candidate in a temperate orbit. Here we report the ALMA detection of the star at 1.3 mm wavelength and the discovery of a belt of dust orbiting around it at distances ranging between 1 and 4 au, approximately. Given the low luminosity of the Proxima Centauri star, we estimate a characteristic temperature of about 40 K for this dust, which might constitute the dust component of a small-scale analog to our solar system Kuiper belt. The estimated total mass, including dust and bodies up to 50 km in size, is of the order of 0.01 Earth masses, which is similar to that of the solar Kuiper belt. Our data also show a hint of warmer dust closer to the star. We also find signs of two additional features that might be associated with the Proxima Centauri system, which, however, still require further observations to be confirmed: an outer extremely cold (about 10 K) belt around the star at about 30 au, whose orbital plane is tilted about 45 degrees with respect to the plane of the sky; and additionally, we marginally detect a compact 1.3 mm emission source at a projected distance of about 1.2 arcsec from the star, whose nature is still unknown.

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.

The formation of massive stars is still not well understood. Accumulating a large amount of mass infalling within a single entity in spite of radiation pressure is possible if, among several other conditions, enough thermal energy is released. Despite numerous water line observations, with the Herschel Space Observatory, in most of the sources observations were not able to trace the emission from the hot core around the newly forming protostellar object. We want to probe the physical conditions and water abundance in the inner layers of the host protostellar object NGC7538-IRS1 using a highly excited H2O line. Water maser models predict that several THz water masers should be detectable in these objects. We present SOFIA observations of the o-H2O 8(2,7)-7(3,4) line at 1296.41106 GHz and a 6(1,6)-5(2,3) 22 GHz e-MERLIN map of the region (first-ever 22 GHz images made after the e-MERLIN upgrade). In order to be able to constrain the nature of the emission - thermal or maser - we use near-simultaneous observations of the 22 GHz water maser performed with the Effelsberg radiotelescope and e-MERLIN. A thermal water model using the RATRAN radiative transfer code is presented based on HIFI pointed observations. Molecular water abundances are derived for the hot core. The H2O 8(2,7)- 7(3,4) line is detected toward NGC7538-IRS1 with one feature at the source velocity (-57.7 km/s) and another one at -48.4 km/s. We propose that the emission at the source velocity is consistent with thermal excitation and is excited in the innermost part of the IRS1a massive protostellar object's closest circumstellar environment. The other emission is very likely the first detection of a water THz maser line, pumped by shocks due to IRS1b outflow, in a star-forming region. Assuming thermal excitation of the THz line, the water abundance in NGC7538-IRS1's hot core is estimated to be 5.2x10^{-5} with respect to H2.

The mechanisms responsible for heating the extended atmospheres of early-M spectral-type supergiants are poorly understood. So too is the subsequent role these mechanisms play in driving the large mass-loss rates of these stars. Here we present ALMA long (i.e., $\sim$16 km) baseline 338 GHz (0.89 mm) continuum observations of the free-free emission in the extended atmosphere of the M2 spectral-type supergiant Betelgeuse. The spatial resolution of 14 mas exquisitely resolves the atmosphere, revealing it to have a mean temperature of 2760 K at $\sim$1.3 R$_{\star}$, which is below both the photospheric effective temperature ($T_{\textrm{eff}} = 3690$ K) and the temperatures at $\sim$$2 $R$_{\star}$. This is unambiguous proof for the existence of an inversion of the mean temperature in the atmosphere of a red supergiant. The emission is clearly not spherically symmetric with two notable deviations from a uniform disk detected in both the images and visibilities. The most prominent asymmetry is located in the north-east quadrant of the disk and is spatially resolved showing it to be highly elongated with an axis-ratio of 2.4 and occupying $\sim$$5\%$ of the disk projected area. Its temperature is approximately 1000 K above the measured mean temperature at 1.3 R$_{\star}$. The other main asymmetry is located on the disk limb almost due east of the disk center and occupies $\sim$$3\%$ of the disk projected area. Both emission asymmetries are clear evidence for localized heating taking place in the atmosphere of Betelgeuse. We suggest that the detected localized heating is related to magnetic activity generated by large-scale photospheric convection.

We provide an update for the astrometric solution for the Type II supernova progenitor Betelgeuse using the revised Hipparcos Intermediate Astrometric Data (HIAD) of van Leeuwen, combined with existing VLA and new e-MERLIN and ALMA positions. The 2007 Hipparcos refined abscissa measurements required the addition of so-called Cosmic Noise of 2.4 mas to find an acceptable 5-parameter stochastic solution. We find that a measure of radio Cosmic Noise should also be included for the radio positions because surface inhomogeneities exist at a level significant enough to introduce additional intensity centroid uncertainty. Combining the 2007 HIAD with the proper motions based solely on the radio positions leads to a parallax of $\pi = 5.27\pm 0.78$ mas ($190^{+33}_{-25}$ pc), smaller than the Hipparcos 2007 value of $6.56\pm 0.83$ mas ($152^{+22}_{-17}$ pc; van Leeuwen 2007). Furthermore, combining the VLA and new e-MERLIN and ALMA radio positions with the 2007 HIAD, and including radio Cosmic Noise of 2.4 mas, leads to a nominal parallax solution of $4.51 \pm 0.80$ mas ($222^{+48}_{-34}$ pc), which while only $0.7\sigma$ different from the 2008 solution of Harper et al. it is $2.6\sigma$ different from the solution of van Leeuwen. An accurate and precise parallax for Betelgeuse is always going to be difficult to obtain because it is small compared to the stellar angular diameter ($\theta=44$ mas). We outline an observing strategy, utilizing future mm and sub-mm high-spatial resolution interferometry that must be used if substantial improvements in precision and accuracy of the parallax and distance are to be achieved.

We perform a likelihood analysis of the constraints from accelerator experiments and astrophysical observations on supersymmetric (SUSY) models with SU(5) boundary conditions on soft SUSY-breaking parameters at the GUT scale. The parameter space of the models studied has 7 parameters: a universal gaugino mass $m_{1/2}$, distinct masses for the scalar partners of matter fermions in five- and ten-dimensional representations of SU(5), $m_5$ and $m_{10}$, and for the $\mathbf{5}$ and $\mathbf{\bar 5}$ Higgs representations $m_{H_u}$ and $m_{H_d}$, a universal trilinear soft SUSY-breaking parameter $A_0$, and the ratio of Higgs vevs $\tan \beta$. In addition to previous constraints from direct sparticle searches, low-energy and flavour observables, we incorporate constraints based on preliminary results from 13 TeV LHC searches for jets + MET events and long-lived particles, as well as the latest PandaX-II and LUX searches for direct Dark Matter detection. In addition to previously-identified mechanisms for bringing the supersymmetric relic density into the range allowed by cosmology, we identify a novel ${\tilde u_R}/{\tilde c_R} - \tilde{\chi}^0_1$ coannihilation mechanism that appears in the supersymmetric SU(5) GUT model and discuss the role of ${\tilde \nu_\tau}$ coannihilation. We find complementarity between the prospects for direct Dark Matter detection and SUSY searches at the LHC.

(abbreviated) We aim to constrain the dust formation histories in the winds of oxygen-rich AGB stars. We have obtained ALMA observations with a spatial resolution of 120x150 mas tracing the dust formation region of a low mass-loss rate and a high mass-loss rate AGB star, R Dor and IK Tau. Emission line profiles of AlO, AlOH and AlCl are detected and are used to derive a lower limit of atomic aluminium incorporated in molecules. We show that the gas-phase aluminium chemistry is completely different in both stars, with a remarkable difference in the AlO and AlOH abundance stratification. The amount of aluminium locked up in these 3 molecules is small, <=1.1e-7, for both stars, i.e. only <=2% of the total aluminium budget. This leaves ample of room for aluminium to be incorporated in grains. A fundamental result is that AlO and AlOH, being the direct precursors of alumina grains, are detected well beyond the onset of the dust condensation proving that the aluminium oxide condensation cycle is not fully efficient. The ALMA observations allow us to quantitatively assess the current generation of theoretical dynamical-chemical models for AGB winds. We discuss how the current proposed scenario of aluminium dust condensation for low mass-loss rate AGB stars at a distance of ~1.5 Rstar, in particular for the stars R Dor and W Hya, poses a challenge if one wishes to explain both the dust spectral features in the spectral energy distribution (SED), in interferometric data, and in polarized light signal. In particular, the estimated grain temperature of Al2O3 is too high for the grains to retain their amorphous structure. We propose that large gas-phase (Al2O3)n-clusters (n>34) can be the potential agents of the broad 11 micron feature in the SED and in the interferometric data and we explain how these large clusters can be formed.

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.

In this Technical Design Report (TDR) we describe the LZ detector to be built at the Sanford Underground Research Facility (SURF). The LZ dark matter experiment is designed to achieve sensitivity to a WIMP-nucleon spin-independent cross section of three times ten to the negative forty-eighth square centimeters.

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

We perform a likelihood analysis of the minimal Anomaly-Mediated Supersymmetry Breaking (mAMSB) model using constraints from cosmology and accelerator experiments. We find that a wino-like or a Higgsino-like neutralino LSP, $m_{\tilde \chi^0_{1}}$, may provide the cold dark matter (DM) with similar likelihood. The upper limit on the DM density from Planck and other experiments enforces $m_{\tilde \chi^0_{1}} \lesssim 3~TeV$ after the inclusion of Sommerfeld enhancement in its annihilations. If most of the cold DM density is provided by the $\tilde \chi_0^1$, the measured value of the Higgs mass favours a limited range of $\tan \beta \sim 5$ (or for $\mu > 0$, $\tan \beta \sim 45$) but the scalar mass $m_0$ is poorly constrained. In the wino-LSP case, $m_{3/2}$ is constrained to about $900~TeV$ and ${m_{\tilde \chi^0_{1}}}$ to $2.9\pm0.1~TeV$, whereas in the Higgsino-LSP case $m_{3/2}$ has just a lower limit $\gtrsim 650TeV$ ($\gtrsim 480TeV$) and $m_{\tilde \chi^0_{1}}$ is constrained to $1.12 ~(1.13) \pm0.02~TeV$ in the $\mu>0$ ($\mu<0$) scenario. In neither case can the anomalous magnetic moment of the muon, ${(g-2)_\mu}$, be improved significantly relative to its Standard Model (SM) value, nor do flavour measurements constrain the model significantly, and there are poor prospects for discovering supersymmetric particles at the LHC, {though there} are some prospects for direct DM detection. On the other hand, if the ${m_{\tilde \chi^0_{1}}}$ contributes only a fraction of the cold DM density, {future LHC $E_T$-based searches for gluinos, squarks and heavier chargino and neutralino states as well as disappearing track searches in the wino-like LSP region will be relevant}, and interference effects enable ${\rm BR}(B_{s, d} \to \mu^+\mu^-)$ to agree with the data better than in the SM in the case of wino-like DM with $\mu > 0$.

Six billion years from now, while evolving on the asymptotic giant branch (AGB), the Sun will metamorphose from a red giant into a beautiful planetary nebula. This spectacular evolution will impact the Solar System planets, but observational confirmations of the predictions of evolution models are still elusive as no planet orbiting an AGB star has yet been discovered. The nearby AGB red giant L2 Puppis (d = 64 pc) is surrounded by an almost edge-on circumstellar dust disk. We report new observations with ALMA at very high angular resolution (18 x 15 mas) in band 7 (f ~ 350 GHz) that allow us to resolve the velocity profile of the molecular disk. We establish that the gas velocity profile is Keplerian within the central cavity of the dust disk, allowing us to derive the mass of the central star L2 Pup A, mA = 0.659 +/- 0.011 +/- 0.041 Msun (+/- 6.6%). From evolutionary models, we determine that L2 Pup A had a near-solar main sequence mass, and is therefore a close analog of the future Sun in 5 to 6 Gyr. The continuum map reveals the presence of a secondary source (B) at a radius of 2 AU contributing fB/ fA = 1.3 +/- 0.1% of the flux of the AGB star. L2 Pup B is also detected in CO emission lines at a radial velocity of vB = 12.2 +/- 1.0 km/s. The close coincidence of the center of rotation of the gaseous disk with the position of the continuum emission from the AGB star allows us to constrain the mass of the companion to mB = 12 +/- 16 MJup. L2 Pup B is most likely a planet or low mass brown dwarf with an orbital period around 5 years. Its continuum brightness and molecular emission suggest that it may be surrounded by an extended molecular atmosphere or an accretion disk. L2 Pup therefore emerges as a promising vantage point on the distant future of our Solar System.

(abreviated) We aim to study the inner-wind structure (R<250 Rstar) of the well-known red supergiant VY CMa. We analyse high spatial resolution (~0".24x0".13) ALMA Science Verification (SV) data in band 7 in which four thermal emission lines of gaseous sodium chloride (NaCl) are present at high signal-to-noise ratio. For the first time, the NaCl emission in the inner wind region of VY CMa is spatially resolved. The ALMA observations reveal the contribution of up to four different spatial regions. The NaCl emission pattern is different compared to the dust continuum and TiO2 emission already analysed from the ALMA SV data. The emission can be reconciled with an axisymmetric geometry, where the lower density polar/rotation axis has a position angle of ~50 degrees measured from north to east. However, this picture can not capture the full morphological diversity, and discrete mass ejection events need to be invoked to explain localized higher-density regions. The velocity traced by the gaseous NaCl line profiles is significantly lower than the average wind terminal velocity, and much slower than some of the fastest mass ejections, signalling a wide range of characteristic speeds for the mass loss. Gaseous NaCl is detected far beyond the main dust condensation region. Realising the refractory nature of this metal halide, this hints at a chemical process preventing all NaCl from condensing onto dust grains. We show that in the case of the ratio of the surface binding temperature to the grain temperature being ~50, only some 10% of NaCl remains in gaseous form, while for lower values of this ratio thermal desorption efficiently evaporates NaCl. Photodesorption by stellar photons seems not to be a viable explanation for the detection of gaseous NaCl at 220 Rstar from the central star, and instead, we propose shock-induced sputtering driven by localized mass ejection events as alternative.

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

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.

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

Distances to AGB stars with optically thick circumstellar shells cannot be determined using optical parallaxes. However, for stars with OH 1612 MHz maser emission emanating from their circumstellar shells, distances can be determined by the phase-lag method. This method combines a linear diameter obtained from a phase-lag measurement with an angular diameter obtained from interferometry. The phase-lag of the variable emission from the back and front sides of the shells has been determined for 20 OH/IR stars in the galactic disk. These measurements are based on a monitoring program with the Nancay radio telescope ongoing for more than 6 years. The interferometric observations are continuing. We estimate that the uncertainties of the distance determination will be ~20%.

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.

The SKA will be a transformational instrument in the study of our local Universe. In particular, by virtue of its high sensitivity (both to point sources and diffuse low surface brightness emission), angular resolution and the frequency ranges covered, the SKA will undertake a very wide range of astrophysical research in the field of nearby galaxies. By surveying vast numbers of nearby galaxies of all types with $\mu$Jy sensitivity and sub-arcsecond angular resolutions at radio wavelengths, the SKA will provide the cornerstone of our understanding of star-formation and accretion activity in the local Universe. In this chapter we outline the key continuum and molecular line science areas where the SKA, both during phase-1 and when it becomes the full SKA, will have a significant scientific impact.

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.

(abbreviated) We aim to study the inner wind of the well-known AGB star CW Leo. Different diagnostics probing different geometrical scales have pointed toward a non-homogeneous mass-loss process: dust clumps are observed at milli-arcsec scale, a bipolar structure is seen at arcsecond-scale and multi-concentric shells are detected beyond 1". We present the first ALMA Cycle 0 band 9 data around 650 GHz. The full-resolution data have a spatial resolution of 0".42x0".24, allowing us to study the morpho-kinematical structure within ~6". Results: We have detected 25 molecular lines. The emission of all but one line is spatially resolved. The dust and molecular lines are centered around the continuum peak position. The dust emission has an asymmetric distribution with a central peak flux density of ~2 Jy. The molecular emission lines trace different regions in the wind acceleration region and suggest that the wind velocity increases rapidly from about 5 R* almost reaching the terminal velocity at ~11 R*. The channel maps for the brighter lines show a complex structure; specifically for the 13CO J=6-5 line different arcs are detected within the first few arcseconds. The curved structure present in the PV map of the 13CO J=6-5 line can be explained by a spiral structure in the inner wind, probably induced by a binary companion. From modeling the ALMA data, we deduce that the potential orbital axis for the binary system lies at a position angle of ~10-20 deg to the North-East and that the spiral structure is seen almost edge-on. We infer an orbital period of 55 yr and a binary separation of 25 au (or ~8.2 R*). We tentatively estimate that the companion is an unevolved low-mass main-sequence star. The ALMA data hence provide us for the first time with the crucial kinematical link between the dust clumps seen at milli-arcsecond scale and the almost concentric arcs seen at arcsecond scale.

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 outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science.

Multi-wavelength centimeter continuum observations of non-dusty, non-pulsating K spectral-type red giants directly sample their chromospheres and wind acceleration zones. Such stars are feeble emitters at these wavelengths however, and previous observations have provided only a small number of modest S/N measurements slowly accumulated over three decades. We present multi-wavelength Karl G. Jansky Very Large Array thermal continuum observations of the wind acceleration zones of two dust-free red giants, Arcturus (Alpha Boo: K2 III) and Aldebaran (Alpha Tau: K5 III). Importantly, most of our observations of each star were carried out over just a few days, so that we obtained a snapshot of the different stellar atmospheric layers sampled at different wavelengths, independent of any long-term variability. We report the first detections at several wavelengths for each star including a detection at 10 cm (3.0 GHz: S band) for both stars and a 20 cm (1.5 GHz: L band) detection for Alpha Boo. This is the first time single luminosity class III red giants have been detected at these continuum wavelengths. Our long-wavelength data sample the outer layers of Alpha Boo's atmosphere where its wind velocity is approaching its terminal value and the ionization balance is becoming frozen-in. For Alpha Tau, however, our long-wavelength data are still sampling its inner atmosphere, where the wind is still accelerating probably due to its lower mass-loss rate. We compare our data with published semi-empirical models based on ultraviolet data, and the marked deviations highlight the need for new atmospheric models to be developed. Spectral indices are used to discuss the possible properties of the stellar atmospheres, and we find evidence for a rapidly cooling wind in the case of Alpha Boo. Finally, we develop a simple analytical wind model for Alpha Boo based on our new long-wavelength flux measurements.

We present a comprehensive view of the W51B HII region complex and the W51C supernova remnant (SNR) using new radio observations from the VLA, VLBA, MERLIN, JCMT, and CSO along with archival data from Spitzer, ROSAT, ASCA, and Chandra. Our VLA data include the first 400 cm (74 MHz) continuum image of W51 at high resolution (88 arcsec). The 400 cm image shows non-thermal emission surrounding the G49.2-0.3 HII region, and a compact source of non-thermal emission (W51B_NT) coincident with the previously-identified OH (1720 MHz) maser spots, non-thermal 21 and 90 cm emission, and a hard X-ray source. W51B_NT falls within the region of high likelihood for the position of TeV gamma-ray emission. Using the VLBA three OH (1720 MHz) maser spots are detected in the vicinity of W51B_NT with sizes of 60 to 300 AU and Zeeman effect magnetic field strengths of 1.5 to 2.2 mG. The multiwavelength data demonstrate that the northern end of the W51B HII region complex has been partly enveloped by the advancing W51C SNR and this interaction explains the presence of W51B_NT and the OH masers. This interaction also appears in the thermal molecular gas which partially encircles W51B_NT and exhibits narrow pre-shock (DeltaV 5 km/s) and broad post-shock (DeltaV 20 km/s) velocity components. RADEX radiative transfer modeling of these two components yield physical conditions consistent with the passage of a non-dissociative C-type shock. Confirmation of the W51B/W51C interaction provides additional evidence in favor of this region being one of the best candidates for hadronic particle acceleration known thus far.

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

Silicon monoxide maser emission has been detected in the circumstellar envelopes of many evolved stars in various vibrationally-excited rotational transitions. It is considered a good tracer of the wind dynamics close to the photosphere of the star. We have investigated the polarization morphology in the circumstellar envelope of an AGB star, R Cas. We mapped the linear and circular polarization of SiO masers in the v=1, J=1-0 transition. The linear polarization is typically a few tens of percent while the circular polarization is a few percent. The fractional polarization tends to be higher for emission of lower total intensity. We found that, in some isolated features the fractional linear polarization appears to exceed 100%. We found the Faraday rotation is not negligible but is ~15 deg., which could produce small scale structure in polarized emission whilst total intensity is smoother and partly resolved out. The polarization angles vary considerably from feature to feature but there is a tendency to favour the directions parallel or perpendicular to the radial direction with respect to the star. In some features, the polarization angle abruptly flips 90 deg. We found that our data are in the regime where the model of Goldreich et al (1973) can be applied and the polarization angle flip is caused when the magnetic field is at close to 55 deg. to the line of sight. The polarization angle configuration is consistent with a radial magnetic field although other configurations are not excluded.

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

This document defines the high level metadata necessary to describe the physical parameter space of observed or simulated astronomical data sets, such as 2D-images, data cubes, X-ray event lists, IFU data, etc.. The Characterisation data model is an abstraction which can be used to derive a structured description of any relevant data and thus to facilitate its discovery and scientific interpretation. The model aims at facilitating the manipulation of heterogeneous data in any VO framework or portal. A VO Characterisation instance can include descriptions of the data axes, the range of coordinates covered by the data, and details of the data sampling and resolution on each axis. These descriptions should be in terms of physical variables, independent of instrumental signatures as far as possible. Implementations of this model has been described in the IVOA Note available at: http://www.ivoa.net/Documents/latest/ImplementationCharacterisation.html Utypes derived from this version of the UML model are listed and commented in the following IVOA Note: http://www.ivoa.net/Documents/latest/UtypeListCharacterisationDM.html An XML schema has been build up from the UML model and is available at: http://www.ivoa.net/xml/Characterisation/Characterisation-v1.11.xsd

The star-forming region W75N hosts bright OH masers that are observed to be variable. We present observations taken in 2008 of the ground-state OH maser transitions with the Very Long Baseline Array (VLBA) and the Multi-Element Radio-Linked Interferometer Network (MERLIN) and with the Nancay Radio Telescope in 2011. Several of the masers in W75N were observed to be flaring, with the brightest 1720-MHz maser in excess of 400 Jy. The 1720-MHz masers appear to be associated with the continuum source VLA 1, unlike the bright flaring 1665- and 1667-MHz masers, which are associated with VLA 2. The 1720-MHz masers are located in an outflow traced by water masers and are indicative of very dense molecular material near the H II region. The magnetic field strengths are larger in the 1720-MHz maser region than in most regions hosting only main-line OH masers. The density falls off along the outflow, and the order of appearance of different transitions of OH masers is consistent with theoretical models. The 1665- and 1667-MHz VLBA data are compared against previous epochs over a time baseline of over 7 years. The median maser motion is 3.5 km/s, with a scatter that is comparable to thermal turbulence. The general pattern of maser proper motions observed in the 1665- and 1667-MHz transitions is consistent with previous observations.

We present multi-epoch, total intensity, high-resolution images of 43GHz, v=1, J=1-0 SiO maser emission toward the Mira variable R Cas. In total we have 23 epochs of data for R Cas at approximate monthly intervals over an optical pulsation phase range from 0.158 to 1.78. These maps show a ring-like distribution of the maser features in a shell, which is assumed to be centred on the star at a radius of 1.6 to 2.3 times the stellar radii. It is clear from these images that the maser emission is significantly extended around the star. At some epochs a faint outer arc can be seen at 2.2 stellar radii. The intensity of the emission waxes and wanes during the stellar phase. Some maser features are seen infalling as well as outflowing. We have made initial comparisons of our data with models by Gray et. al. (2009).

The nature of maser emission means that the apparent angular size of an individual maser spot is determined by the amplification process as well as by the instrinsic size of the emitting cloud. Highly sensitive MERLIN radio interferometry images spatially and spectrally resolve water maser clouds around evolved stars. We measured the properties of clouds around the red supergiant S Per and the AGB stars IK Tau, RT Vir, U Her and U Ori, to test maser beaming theory. Spherical clouds are expected to produce an inverse relationship between maser intensity and apparent size, which would not be seen from cylindrical or slab-like regions. We analysed the maser properties, in order to estimate the saturation state, and investigated the variation of observed spot size with intensity and across the spectral line profiles. Circumstellar masers emanate from discrete clouds from about one to 20 AU in diameter depending on the star. Most of the maser features have negative excitation temperatures close to zero and modest optical depths, showing that they are mainly unsaturated. Around S Per and (at most epochs) RT Vir and IK Tau, the maser component size shrinks with increasing intensity. In contrast, the masers around U Ori and U Her tend to increase in size, with a larger scatter. The water masers from S Per, RT Vir and IK Tau are mainly beamed into spots with an observed angular size much smaller than the emitting clouds and smallest of all at the line peaks. This suggests that the masers are amplification-bounded, emanating from approximately spherical clouds. Many of the masers around U Her and U Ori have apparent sizes which are more similar to the emitting clouds and have less or no dependence on intensity, suggesting that these masers are matter-bounded. This is consistent with an origin in flattened clouds and these two stars have shown other behaviour indicating the presence of shocks.

The aim of this study is to investigate the physical properties of molecular envelopes of planetary nebulae in their earliest stages of evolution. Using the 100m telescope at Effelsberg, we have undertaken a high sensitivity discrete source survey for the first excited state of OH maser emission (J=5/2, 2PI3/2 at 6GHz) in the direction of planetary and proto-planetary nebulae exhibiting 18cm OH emission (main and/or satellite lines), and we further validate our detections using the Nan\c{c}ay radio telescope at 1.6-1.7GHz and MERLIN interferometer at 1.6-1.7 and 6GHz. Two sources have been detected at 6035MHz (5cm), both of them are young (or very young) planetary nebulae. The first one is a confirmation of the detection of a weak 6035MHz line in Vy 2-2. The second one is a new detection, in K 3-35, which was already known to be an exceptional late type star because it exhibits 1720MHz OH emission. The detection of 6035MHz OH maser emission is confirmed by subsequent observations made with the MERLIN interferometer. These lines are very rarely found in evolved stars. The 1612MHz masers surround but are offset from the 1720 and 6035MHz masers which in turn lie close to a compact 22GHz continuum source embedded in the optical nebula.

We report results of an extensive observational campaign of the 6 cm formaldehyde maser in the young massive stellar object IRAS18566+0408 (G37.55+0.20) conducted from 2002 to 2009. Using Arecibo, VLA, and GBT, we discovered quasi-periodic formaldehyde flares (P ~ 237 days). Based on Arecibo observations, we also discovered correlated variability between formaldehyde (H2CO) and methanol (CH3OH) masers. The H2CO and CH3OH masers are not spatially coincident, as demonstrated by different line velocities and high angular resolution MERLIN observations. The flares could be caused by variations in the infrared radiation field, possibly modulated by periodic accretion onto a young binary system.

We report sub-arcsecond scale radio continuum observations of a field of 35 by 22 arcmin centred in NGC3351 obtained with the Multi-Element Radio Linked Interferometer Network (MERLIN). We found 23 radio sources in this field, 6 of which are projected within the D$_{25}$ isophote of the galaxy, and 3 are located inside the central 100 arcsec in radius. Two of these three are significantly extended, while the third one is relatively compact. This one is the only source with a previously detected counterpart at other wavelengths and could constitute the radio counterpart of a young supernova remnant. The other two are probably related to jets from a background AGN. We are not able to detect individual supernovae or SNRs in the central region ($r<600$ pc) of the galaxy. This could imply that the ionising populations of the circumnuclear star-forming regions (CNSFRs) are too young (less than a few Myr) to host supernovae. Also the diffusion length of the relativistic electrons in the ISM associated with the SN from previous events of star formation seems to be larger than our maximum resolution of 50 pc in radius. Detecting the thermal bremsstrahlung emission from the circumnuclear HII regions probably requires deeper observations.

We investigate which structures the 6.7 GHz methanol masers trace in the environment of high-mass protostar candidates by observing a homogenous sample of methanol masers selected from Torun surveys. We also probed their origins by looking for associated H II regions and IR emission. We selected 30 methanol sources with improved position accuracies achieved using MERLIN and another 3 from the literature. We imaged 31 of these using the European VLBI Network's expanded array of telescopes with 5-cm (6-GHz) receivers. We used the VLA to search for 8.4 GHz radio continuum counterparts and inspected Spitzer GLIMPSE data at 3.6-8 um from the archive. High angular resolution images allowed us to analyze the morphology and kinematics of the methanol masers in great detail and verify their association with radio continuum and mid-infrared emission. A new class of "ring-like" methanol masers in star--forming regions appeared to be suprisingly common, 29 % of the sample. The new morphology strongly suggests that methanol masers originate in the disc or torus around a proto- or a young massive star. However, the maser kinematics indicate the strong influence of outflow or infall. This suggests that they form at the interface between the disc/torus and a flow. This is also strongly supported by Spitzer results because the majority of the masers coincide with 4.5 um emission to within less than 1 arcsec. Only four masers are associated with the central parts of UC H II regions. This implies that 6.7 GHz methanol maser emission occurs before H II region observable at cm wavelengths is formed.

We have imaged the disc of the young star HL Tau using the VLA at 1.3 cm, with 0.08" resolution (as small as the orbit of Jupiter). The disc is around half the stellar mass, assuming a canonical gas-mass conversion from the measured mass in large dust grains. A simulation shows that such discs are gravitationally unstable, and can fragment at radii of a few tens of AU to form planets. The VLA image shows a compact feature in the disc at 65 AU radius (confirming the `nebulosity' of Welch et al. 2004), which is interpreted as a localised surface density enhancement representing a candidate proto-planet in its earliest accretion phase. If correct, this is the first image of a low-mass companion object seen together with the parent disc material out of which it is forming. The object has an inferred gas plus dust mass of approximately 14 M(Jupiter), similar to the mass of a proto-planet formed in the simulation. The disc instability may have been enhanced by a stellar flyby: the proper motion of the nearby star XZ Tau shows it could have recently passed the HL Tau disc as close as ~600 AU.

We present MERLIN observations of Galactic 21-cm HI absorption at an angular resolution of c. 0.1-0.2 arcsec and a velocity resolution of 0.5 km/s, in the direction of three moderately low latitude (-8< b <-12 deg) extragalactic radio sources, 3C111, 3C123 and 3C161, all of which are heavily reddened. HI absorption is observed against resolved background emission sources up to c. 2 arcsec in extent and we distinguish details of the opacity distribution within 1-1.5 arcsec regions towards 3C~123 and 3C~161. This study is the second MERLIN investigation of small scale structure in interstellar HI (earlier work probed Galactic HI in the directions of the compact sources 3C138 and 3C147). The 0.1-arcsec scale is intermediate between HI absorption studies made with other fixed element interferometers with resolution of 1 to 10 arcsec and VLBI studies with resolutions of 10-20 milli-arcsec. At a scale of 1 arcsec (about 500 AU), prominent changes in Galactic HI opacity in excess of 1-1.5 are determined in the direction of 3C161 with a signal-to-noise ratio of at least 10 sigma. Possible fluctuations in the HI opacity at the level of about 1 are detected at the 2.5-3 sigma level in the direction of 3C123.

In this paper we investigate the radio-MIR correlation at very low flux densities using extremely deep 1.4 GHz sub-arcsecond angular resolution MERLIN+VLA observations of a 8'.5 by 8'.5 field centred upon the Hubble Deep Field North, in conjunction with Spitzer 24micron data. From these results the MIR-radio correlation is extended to the very faint (~microJy) radio source population. Tentatively we detect a small deviation from the correlation at the faintest IR flux densities. We suggest that this small observed change in the gradient of the correlation is the result of a suppression of the MIR emission in faint star-forming galaxies. This deviation potentially has significant implications for using either the MIR or non-thermal radio emission as a star-formation tracer of very low luminosity galaxies.

Using extremely deep (rms 3.3 microJy/bm) 1.4GHz sub-arcsecond resolution MERLIN + VLA radio observations of a 8'.5 by 8'.5 field centred upon the Hubble Deep Field North, in conjunction with Spitzer 24 micron data we present an investigation of the radio-MIR correlation at very low flux densities. By stacking individual sources within these data we are able to extend the MIR-radio correlation to the extremely faint (~microJy and even sub-microJy) radio source population. Tentatively we demonstrate a small deviation from the correlation for the faintest MIR sources. We suggest that this small observed change in the gradient of the correlation is the result of a suppression of the MIR emission in faint star-forming galaxies. This deviation potentially has significant implications for using either the MIR or non-thermal radio emission as a star-formation tracer at low luminosities.

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.

A 10-arcmin field around the HDF(N) contains 92 radio sources >40 uJy, resolved by MERLIN+VLA at 0".2-2".0 resolution. 55 have Chandra X-ray counterparts including 18 with a hard X-ray photon index and high luminosity characteristic of a type-II (obscured) AGN. >70% of the radio sources have been classified as starbursts or AGN using radio morphologies, spectral indices and comparisons with optical appearance and MIR emission. Starbursts outnumber radio AGN 3:1. This study extends the VO methods previously used to identify X-ray-selected obscured type-II AGN to investigate whether very luminous radio and X-ray emission originates from different phenomena in the same galaxy. The high-redshift starbursts have typical sizes of 5--10 kpc and star formation rates of ~1000 Msun/yr. There is no correlation between radio and X-ray luminosities nor spectral indices at z>~1.3. ~70% of both the radio-selected AGN and the starburst samples were detected by Chandra. The X-ray luminosity indicates the presence of an AGN in at least half of the 45 cross-matched radio starbursts, of which 11 are type-II AGN including 7 at z>1.5. This distribution overlaps closely with the X-ray detected radio sources which were also detected by SCUBA. Stacked 1.4-GHz emission at the positions of radio-faint X-ray sources is correlated with X-ray hardness. Most extended radio starbursts at z>1.3 host X-ray selected obscured AGN. Radio emission from most of these ultra-luminous objects is dominated by star formation but it contributes less than 1/3 of their X-ray luminosity. Our results support the inferences from SCUBA and IR data, that at z>1.5, star formation is an order of magnitude more extended and more copious, it is closely linked to AGN activity and it is triggered differently, compared with star formation at lower redshifts.

We characterize the molecular environment of candidate massive young stellar objects (MYSOs) signposted by methanol masers. Single pixel observations of 10 transitions of HCO^+, CO and CS isotopomers were carried out, using the IRAM 30m telescope. We studied a sample of 28 targets for which the 6.7GHz maser emission positions are known with a sub-arcsecond accuracy. The systemic velocity inferred from the optically thin lines agrees within 3km/s with the central velocity of the maser emission for most of the sources. About 64% of the sources show line wings in one or more transitions of CO, HCO^+ and CS species, indicating the presence of molecular outflows. Comparison of the widths of line wings and methanol maser emission suggests that the 6.7GHz maser line traces the environment of MYSO of various kinematic regimes. Therefore conditions conducive for the methanol maser can exist in the inner parts of molecular clouds or circumstellar discs as well as in the outer parts associated with molecular outflows. Calculations of the physical conditions based on the CO and HCO^+ lines and the CS line intensity ratios refine the input parameters for maser models.

We present recent results from very deep MERLIN/VLA 1.4GHz observations of an 8.5-by-8.5 arcminute field centred on the the Hubble Deep Field North, in conjunction with BVIz catalogues of GOODS/ACS sources detected in this field. In an extension of work by Muxlow et al. (also presented at this conference) on the statistical properties of over 8000 sources, the results presented here take advantage of the deep (rms: 3.3 microJy/bm), sub-arcsecond radio imaging of this field to characterise the properties at 1.4GHz of a number of individual radio counterparts to GOODS/ACS sources.

We use Virtual Observatory methods to investigate the association between radio and X-ray emission at high redshifts. Fifty-five of the 92 HDF(N) sources resolved by combining MERLIN+VLA data were detected by Chandra, of which 18 are hard enough and bright enough to be obscured AGN. The high-z population of microJy radio sources is dominated by starbursts an order of magnitude more active and more extended than any found at z<1 and at least a quarter of these simultaneously host highly X-ray-luminous obscured AGN.

We present recent and ongoing results from extremely deep 18 day MERLIN + VLA 1.4GHz observations (rms: 3.3microJy/bm) of an 8.5-by-8.5 arcminute field centred upon the Hubble Deep Field North. This area of sky has been the subject of some of the deepest observations ever made over a wide range of frequencies, from X-rays to the radio. The results presented here use our deep, sub-arcsecond radio imaging of this field to characterise the radio structures of the several hundred GOODS Spitzer MIR sources in this field. These MIR sources primarily trace the luminous starburst sources. A significant proportion of the MIR sources are detected and resolved by our radio observations, allowing these observations to trace the IR/Radio correlation for galaxies over ~7 orders of magnitude, extending it to ever lower luminosities.

Deep combination radio observations at 1.4GHz with the VLA and MERLIN have imaged a region 10 arcminutes square surrounding the Hubble Deep Field North (HDF-N). Initial studies of the weak radio source population have shown that the proportion of starburst systems increases with decreasing radio flux density with more than 70% of radio sources being starburst in nature at flux densities less than S1.4GHz)~70microJy. The recently published GOODS ACS field overlaps this area, and here we present the results of a follow-up statistical study of the very weak radio sources (S(1.4GHz)<40microJy) in an 8.5 arcminute square field centred on the HDF-N which contains the region of overlap. Radio emission at the level of a few microJy are statistically detected associated with ACS galaxies brighter than a z-band magnitude of 25. These very faint radio sources are extended starburst systems with average radii in the range 0.6 to 0.8 arcseconds and for those with measured redshifts, radio luminosities typically several times that of the nearby well-studies starburst galaxy M82.

Star-formation and the Starburst phenomenon are presented with respect to a number of nearby star-forming galaxies where our understanding of the process can be calibrated. Methods of estimating star-formation rates are discussed together with the role played in the investigation of the process by multi-wavelength studies of a few selected starburst galaxies (especially the well studied galaxy M82). Our understanding of nearby systems allows us to study the star-formation history of the Universe by observing high-redshift starburst galaxies. These begin to dominate the radio source populations at centimetric wavelengths at flux densities below a few 10s of Jy. New very sensitive, high resolution telescopes in the sub-mm and radio will revolutionize our understanding of these distant star-forming systems, some of which may contain embedded AGN.

We present a summary of the major contributions to the Special Session on Data Management held at the IAU General Assembly in Prague in 2006. While recent years have seen enormous improvements in access to astronomical data, and the Virtual Observatory aims to provide astronomers with seamless access to on-line resources, more attention needs to be paid to ensuring the quality and completeness of those resources. For example, data produced by telescopes are not always made available to the astronomical community, and new instruments are sometimes designed and built with insufficient planning for data management, while older but valuable legacy data often remain undigitised. Data and results published in journals do not always appear in the data centres, and astronomers in developing countries sometimes have inadequate access to on-line resources. To address these issues, an 'Astronomers Data Manifesto' has been formulated with the aim of initiating a discussion that will lead to the development of a 'code of best practice' in astronomical data management.

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

The silicate carbon star V778 Cyg is a source of 22 GHz water maser emission which was recently resolved by MERLIN. Observations revealed an elongated S-like structure along which the velocities of the maser features show a linear dependence on the impact parameter. This is consistent with a doubly-warped m=2 disc observed edge-on. Water masers and silicate dust emission (detected by IRAS and ISO) have a common origin in O-rich material and are likely to be co-located in the disc. We propose a detailed self-consistent model of a masing gas-dust disc around a companion to the carbon star in a binary system, which allows us to estimate the companion mass of 1.7 +- 0.1 M_sun, the disc radius of 40 +-3 AU and the distance between companions of about 80 AU. Using a dust-gas coupling model for water masing, we calculate the maser power self-consistently, accounting for both the gas and the dust energy balances. Comparing the simulation results with the observational data, we deduce the main physical parameters of the masing disc, such as the gas and dust temperatures and their densities. We also present an analysis of the stability of the disc.

Various scenarios have been proposed to explain the presence of silicate features associated with carbon stars, such as V778 Cyg. We have attempted to constrain these theories by means of mapping water maser mission from V778 Cyg. The 22 GHz water maser emission from this star has been mapped using MERLIN with an astrometric accuracy of 25 mas. The spatially- and kinematically-resolved maser complex is displaced by about 190 mas from the position of the C-star as measured 10 years earlier using Tycho. Our simulations and analysis of available data show that this position difference is unlikely to be due to proper motion if V778 Cyg is at the assumed distance of 1.4 kpc. The maser components seem to form a distorted S-shaped structure extended over 18 mas with a clear velocity gradient. We propose a model which explains the observed water maser structure as an O-rich warped disk around a companion of the C-star in V 778 binary system, which is seen almost edge-on. Analysis of observational data, especially those obtained with MERLIN, suggests that V778 Cyg (and, by implication, other silicate carbon stars) are binary systems composed of a C-rich star and a companion which stores circumstellar O-rich material.

High resolution MERLIN 5 GHz observations (0.04") of 7 Seyfert galaxies, selected as the ones previously showing evidence of collimated ejection, have been compared with high resolution archive HST data. The radio maps reveal rich structures in all the galaxies. NGC 2639 and TXFS 2226-184 have multiple knot parsec-scale extended structures, Mrk 1034, Mrk 1210, NGC 4922C and NGC 5506 reveal one-sided jets, while IC 1481 exhibits a jet-like extension. The close correlation between the radio-emitting relativistic plasma and the ionized gas in the inner regions of these galaxies allows us to study in detail the physics close to the center of low luminosity AGN.

We present MERLIN observations of OH maser and radio continuum emission from the Ultra Luminous IR Galaxy Markarian 231. The 1665- and 1667-MHz transitions have a combined velocity extent of 720 km/s and show a similar position-velocity structure including a gradient of 1.7 km/s/pc from NW to SE along the 420-pc major axis, steeper in the inner few tens of pc. The maser distribution is modelled as a torus rotating about an axis inclined at ~45deg. We estimate the enclosed mass density to be 320(90) Msun in a flattened distribution, including a central unresolved mass of </=8E+06 Msun. All the maser emission is projected against a region with a radio continuum brightness temperature >/=1E+05 K, giving a maser gain of </=2.2. The 1667:1665-MHz line ratio is close to the LTE ratio of 1.8 consistent with radiatively pumped, unsaturated masers. The size of individual masing regions is in the range 0.25-4 pc with a covering factor close to unity. There are no very bright compact masers, in contrast to galaxies such as the Seyfert 2 Markarian 273 where the masing torus is viewed nearer edge-on. The comparatively modest maser amplification seen from Markarian 231 is consistent with its classification as a Seyfert 1. Most of the radio continuum emission on 50-500 pc scales is probably of starburst origin but the compact peak is 0.4 per cent polarized by a magnetic field running north-south, similar to the jet direction on these scales. There is no close correlation between maser and continuum intensity. Comparisons with other data show that the jet changes direction close the nucleus and suggest that the sub-kpc disc hosting the masers and starburst activity is severely warped.

We present the results of MERLIN polarization mapping of OH masers at 1665 and 1667 MHz towards the Cepheus A star-forming region. The maser emission is spread over a region of 6 arcsec by 10 arcsec, twice the extent previously detected. In contrast to the 22 GHz water masers, the OH masers associated with H II regions show neither clear velocity gradients nor regular structures. We identified ten Zeeman pairs which imply a magnetic field strength along the line-of-sight from -17.3 to +12.7 mG. The magnetic field is organised on the arcsecond scale, pointing towards us in the west and away from us in the east side. The linearly polarized components, detected for the first time, show regularities in the polarization position angles depending on their position. The electric vectors of OH masers observed towards the outer parts of H II regions are consistent with the interstellar magnetic field orientation, while those seen towards the centres of H II regions are parallel to the radio-jets. A Zeeman quartet inside a southern H II region has now been monitored for 25 years; we confirm that the magnetic field decays monotonically over that period.

Eighteen days of MERLIN data and 42 hours of A-array VLA data at 1.4 GHz have been combined to image a 10-arcmin field centred on the Hubble Deep and Flanking Fields (HDF and HFF). A complete sample of 92 radio sources with 1.4-GHz flux densities above 40 microJy has been imaged using MERLIN+VLA. The images are amongst the most sensitive yet made at 1.4 GHz, with rms noise levels of 3.3 microJy/beam in the 0.2-arcsec images. Virtually all the sources are resolved, with angular sizes in the range 0.2 to 3 arcsec. No additional sources were detected down to 23 microJy in the central 3 arcmin, indicating that sources fainter than 40 microJy are heavily resolved with MERLIN and must have typical angular sizes greater than 0.5 arcsec. Compact radio sources were used to align the optical data to the ICRF, to <50 mas in the HDF. We find a statistical association of very faint (2 microJy and above) radio sources with optically bright HDF galaxies down to about 23 mag. Of the 92 radio sources above 40 microJy, about 85 percent are identified with galaxies brighter than about I = 25 mag; the remaining 15 percent are associated with optically faint systems. We identify several very red, optically faint systems including the the strongest sub-mm source in the HDF, HDF850.1. 72 percent of the radio sources are starburst or AGN-type systems; the remainder are unclassified. The proportion of starburst systems increases with decreasing flux density; below 100 microJy 70 percent of the sources are starburst-type systems in the redshift range 0.3 -- 1.3. Chandra detections are associated with 55 of the 92 radio sources but their X-ray flux densities do not appear to be correlated with the radio flux densities or morphologies.

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.

We present the first maps of 4765-MHz OH masers in two star-forming regions Cepheus A and W75N, made with Multi-Element Radio Linked Interferometer Network. In Cep A the emission has an arc-like structure of size 40 mas with a clear velocity gradient, and lies at the edge of H II region 3b Over a period of 8 weeks the maser diminished in intensity by a factor of 7 but its structure remained stable. This structure coincides with a newly mapped 1720-MHz maser in Cep A within the positional error, and matches it in velocity. The 1720-MHz line exhibits Zeeman splitting that corresponds to a magnetic field strength of -17.3 mG. In W75N the excited 4765-MHz OH maser has a linear structure of size 45 mas with a well defined velocity gradient, and lies at the edge of H II region VLA 1. This structure coincides in position and velocity with the 1720-MHz masers. We conclude that in both sources the 4765-MHz emission coexists with 1720-MHz emission in the same volume of gas. In such a case the physical conditions in these regions are tightly constrained by the maser-pumping models.

As part of a multi-wavelength campaign to observe the 2000--2002 outburst of the prototypical symbiotic star Z Andromedae, we observed this object six times each with the Multi-Element Radio Linked Interferometer Network (MERLIN) and Very Large Array (VLA). The radio flux varied significantly during the course of the optical outburst at all three observation frequencies (1.4, 5, and 15 GHz). A jet-like extension was present in the 2001 September MERLIN image and appeared to be aligned perpendicularly to the plane of the binary orbit. Assuming that the ejection took place at the beginning of the optical outburst, the 0.06 arcsec separation between the peak of the extended emission and the central core implies that the ejected material was moving with a velocity of \~400 km/s. This extended emission faded on a timescale of ~months and was not detected at any other epoch. We consider the implications of jets being a component of a ``prototypical'' symbiotic system and compare properties of the observed jet of Z And with those of the jets in X-ray binary systems.

The 6 cm formaldehyde (H2CO) maser sources in the compact HII regions NGC 7538-IRS1 and G29.96-0.02 have been imaged at high resolution (beams < 50 mas). Using the VLBA and MERLIN, we find the angular sizes of the NGC 7538 masers to be ~10 mas (30 AU) corresponding to brightness temperatures ~10^8 K. The angular sizes of the G29.96-0.02 masers are ~20 mas (130 AU) corresponding to brightness temperatures ~10^7 K. Using the VLA, we detect 2 cm formaldehyde absorption from the maser regions. We detect no emission in the 2 cm line, indicating the lack of a 2 cm maser and placing limits on the 6 cm excitation process. We find that both NGC 7538 maser components show an increase in intensity on 5-10 year timescales while the G29.96-0.02 masers show no variability over 2 years. A search for polarization provides 3-sigma upper limits of 1% circularly polarized and 10% linearly polarized emission in NGC 7538 and of 15% circularly polarized emission in G29.96-0.02. A pronounced velocity gradient of 28 km/s/arcsecond (1900 km/s/pc) is detected in the NGC 7538 maser gas.

We present MERLIN maps of the 22-GHz H2O masers around four low-mass late-type stars (IK Tau U Ori, RT Vir and U Her), made with an angular resolution of ~ 15 milliarcsec and a velocity resolution of 0.1 km s-1. The H2O masers are found in thick expanding shells with inner radii ~ 6 to 16 au and outer radii four times larger. The expansion velocity increases radially through the H2O maser regions, with logarithmic velocity gradients of 0.5--0.9. IK Tau and RT Vir have well-filled H2O maser shells with a spatial offset between the near and far sides of the shell, which suggests that the masers are distributed in oblate spheroids inclined to the line of sight. U Ori and U Her have elongated poorly-filled shells with indications that the masers at the inner edge have been compressed by shocks; these stars also show OH maser flares. MERLIN resolves individual maser clouds, which have diameters of 2 -- 4 au and filling factors of only ~ 0.01 with respect to the whole H2O maser shells. The CSE velocity structure gives additional evidence the maser clouds are density bounded. Masing clouds can be identified over a similar timescale to their sound crossing time (~2 yr) but not longer. The sizes and observed lifetimes of these clouds are an order of magnitude smaller than those around red supergiants, similar to the ratio of low-mass:high-mass stellar masses and sizes. This suggests that cloud size is determined by stellar properties, not local physical phenomena in the wind.

MERLIN and VLBA observations of the 1720 MHz maser emission from the OH molecule in the supernova remnant IC443 are presented. Based on MERLIN data with a resolution of 160 mas, the deconvolved sizes of the maser sources are in the range 90 to 180 mas (135 to 270 AU). The 12 mas resolution VLBA images show compact cores with sizes in the range 15 to 55 mas. The maser brightness temperatures are (2-34)x10^6 K for the MERLIN sources and (5-19)x10^8 K for the VLBA cores, in agreement with theory. Unlike the Zeeman Stokes V profiles observed in other OH (1720 MHz) SNR masers, single-handed circular polarization line profiles are observed in IC443 on all angular scales from 1000 to 10 mas resolution. For one line component, the observed line width is 0.24+/-0.07 km/s, compared to an estimated Doppler width of 0.49 km/s. This discrepancy in line widths can be accounted for if the maser emission arises from an elongated ellipsoidal region of masing gas.

Despite extensive observational efforts, the brightest sub--mm source in the Hubble Deep Field, HDF850.1, has failed to yield a convincing optical/infrared identification almost 4 years after its discovery. This failure is all the more notable given the availability of supporting multi-frequency data of unparalleled depth, and sub-arcsec positional accuracy for the sub-mm/mm source. Consequently, HDF850.1 has become a test case of the possibility that the most violently star-forming objects in the universe are too red and/or distant to be seen in the deepest optical images. Here we report the discovery of the host galaxy of HDF850.1. This object has been revealed by careful analysis of a new, deep K-prime image of the HDF obtained with the Subaru 8.2-m telescope. Its reality is confirmed by a similar analysis of the HST NICMOS F160W image of the same region. This object is extremely faint (K=23.5), clumpy (on sub-arcsec scales) and very red (I-K > 5.2; H-K = 1.4 +- 0.35). The likelihood that it is the correct identification is strongly reinforced by a reanalysis of the combined MERLIN+VLA 1.4-GHz map of the field which yields a new radio detection of HDF850.1 only 0.1 arcsec from the new near-ir counterpart, and with sufficient positional accuracy to exclude all previously considered alternative optical candidates. We have calculated new confidence limits on the estimated redshift of HDF850.1 and find z = 4.1 +- 0.5. We also calculate that the flux density of HDF850.1 has been boosted by a factor of ~3 through lensing by the intervening elliptical 3-586.0, consistent with predictions that a small but significant fraction of blank-field sub-mm sources are lensed by foreground galaxies. We discuss the wider implications of these results for the sub-mm population and cosmic star-formation history.

We have made a deep search for radio emission from all the northern hemisphere supersoft X-ray sources using the VLA and MERLIN telescopes, at 5 and 8.4 GHz. Three previously undetected sources: T Pyx, V1974 Cygni and RX J0019.8+2156 were imaged in quiescence using the VLA in order to search for any persistent emission. No radio emission was detected in any of the VLA fields down to a typical 1 sigma RMS noise of 20 uJy/beam, however, 17 new point sources were detected in the fields with 5 GHz fluxes between 100 and 1500 uJy giving an average 100 uJy-source density of around 200 per square degree, comparable to what was found in the MERLIN HDF survey. The persistent source AG Draconis was observed by MERLIN to provide a confirmation of previous VLA observations and to investigate the source at a higher resolution. The core is resolved at the milliarcsec scale into two components which have a combined flux of around 1 mJy. It is possible that we are detecting nebulosity which is becoming resolved out by the higher MERLIN resolution. We have investigated possible causes of radio emission from a wind environment, both directly from the secondary star, and also as a consequence of the high X-ray luminosity from the white dwarf. There is an order of magnitude discrepancy between observed and modelled values which can be explained by the uncertainty in fundamental quantities within these systems.

MERLIN images of Mrk 348 at 22 GHz show water maser emission at 0.02 - 0.11 Jy, within approximately 0.8 pc of the nucleus. This is the first direct confirmation that molecular material exists close to the Seyfert 2 nucleus. Mrk 348 was observed in 2000 May one month after Falcke et al. (2000) first identified the maser in single-dish spectra. The peak maser flux density has increased about threefold. The masing region is < 0.6 pc in radius. The flux density of radio continuum emission from the core has been rising for about 2 years. The maser-core separation is barely resolved but at the 3 sigma significance level they are not coincident along the line of sight. The masers lie in the direction of the northern radio lobes and probably emanate from material shocked by a jet with velocity close to c. The correlation between the radio continuum increase and maser flare is explained as arising from high level nuclear activity through a common excitation mechanism although direct maser amplification of the core by masers tracing a Keplerian disc is not completely ruled out.

We present adaptive optics imaging of the CLASS gravitational lens system B1359+154 obtained with the Canada-France-Hawaii Telescope (CFHT) in the infrared K band. The observations show at least three brightness peaks within the ring of lensed images, which we identify as emission from multiple lensing galaxies. The results confirm the suspected compound nature of the lens, as deduced from preliminary mass modeling. The detection of several additional nearby galaxies suggests that B1359+154 is lensed by the compact core of a small galaxy group. We attempted to produce an updated lens model based on the CFHT observations and new 5-GHz radio data obtained with the MERLIN array, but there are too few constraints to construct a realistic model at this time. The uncertainties inherent with modeling compound lenses make B1359+154 a challenging target for Hubble constant determination through the measurement of differential time delays. However, time delays will offer additional constraints to help pin down the mass model. This lens system therefore presents a unique opportunity to directly measure the mass distribution of a galaxy group at intermediate redshift.

We have studied the core morphology of the ultraluminous infrared galaxy Mrk 273 by combining a high-resolution adaptive optics near-infrared image with an optical image from the Hubble Space Telescope and interferometric radio continuum data, all at spatial resolutions of 150 mas or better. The near-infrared image reveals that the nucleus has two main components, both of which have radio counterparts. The strongest component (N) shows very similar extended structure in the radio and near-infrared. It has a flat radio spectrum and is resolved into a double-lobed structure (Ne; Nw), with a separation of 90\pm5 mas (70 parsec). A similar structure is detected in the near-infrared. We identify this component as the location of the active nucleus. The second component (SW), strong in the near-infrared but relatively weak in the radio, is located $\sim1$ arcsecond to the southwest. We interpret this as an obscured starburst region associated with the merger. The radio continuum images show a third, strong, component (SE) which has previously been interpreted as a second nucleus. However, it shows no associated optical or near-infrared emission, suggesting that it is in fact a background source.