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 use the Met Office Unified Model to explore the potential of a tidally locked M dwarf planet, nominally Proxima Centauri b irradiated by a quiescent version of its host star, to sustain an atmospheric ozone layer. We assume a slab ocean surface layer, and an Earth-like atmosphere of nitrogen and oxygen with trace amounts of ozone and water vapour. We describe ozone chemistry using the Chapman mechanism and the hydrogen oxide (HO$_x$, describing the sum of OH and HO$_2$) catalytic cycle. We find that Proxima Centauri radiates with sufficient UV energy to initialize the Chapman mechanism. The result is a thin but stable ozone layer that peaks at 0.75 parts per million at 25 km. The quasi-stationary distribution of atmospheric ozone is determined by photolysis driven by incoming stellar radiation and by atmospheric transport. Ozone mole fractions are smallest in the lowest 15 km of the atmosphere at the sub-stellar point and largest in the nightside gyres. Above 15 km the ozone distribution is dominated by an equatorial jet stream that circumnavigates the planet. The nightside ozone distribution is dominated by two cyclonic Rossby gyres that result in localized ozone hotspots. On the dayside the atmospheric lifetime is determined by the HO$_x$ catalytic cycle and deposition to the surface, with nightside lifetimes due to chemistry much longer than timescales associated with atmospheric transport. Surface UV values peak at the substellar point with values of 0.01 W/m$^2$, shielded by the overlying atmospheric ozone layer but more importantly by water vapour clouds.

Radiative transfer is a key component in almost all astrophysical and cosmological simulations. We present Magritte: a modern open-source software library for 3D radiative transfer. It uses a deterministic ray-tracer and formal solver, i.e. it computes the radiation field by tracing rays through the model and solving the radiative transfer equation in its second-order form along a fixed set of rays originating from each point. Magritte can handle structured and unstructured input meshes, as well as smoothed-particle hydrodynamics (SPH) particle data. In this first paper, we describe the numerical implementation, semi-analytic tests and cross-code benchmarks for the non-LTE line radiative transfer module of Magritte. This module uses the radiative transfer solver to self-consistently determine the populations of the quantised energy levels of atoms and molecules using an accelerated Lambda iteration (ALI) scheme. We compare Magritte with the established radiative transfer solvers Ratran (1D) and Lime (3D) on the van Zadelhoff benchmark and present a first application to a simple Keplerian disc model. Comparing with Lime, we conclude that Magritte produces more accurate and more precise results, especially at high optical depth, and that it is faster.

Unravelling the composition and characteristics of gas and dust lost by asymptotic giant branch (AGB) stars is important as these stars play a vital role in the chemical life cycle of galaxies. The general hypothesis of their mass loss mechanism is a combination of stellar pulsations and radiative pressure on dust grains. However, current models simplify dust formation, which starts as a microscopic phase transition called nucleation. Various nucleation theories exist, yet all assume chemical equilibrium, growth restricted by monomers, and commonly use macroscopic properties for a microscopic process. Such simplifications for initial dust formation can have large repercussions on the type, amount, and formation time of dust. By abandoning equilibrium assumptions, discarding growth restrictions, and using quantum mechanical properties, we have constructed and investigated an improved nucleation theory in AGB wind conditions for four dust candidates, TiO$_2$, MgO, SiO and Al$_2$O$_3$. This paper reports the viability of these candidates as first dust precursors and reveals implications of simplified nucleation theories. Monomer restricted growth underpredicts large clusters at low temperatures and overpredicts formation times. Assuming the candidates are present, Al$_2$O$_3$ is the favoured precursor due to its rapid growth at the highest considered temperatures. However, when considering an initially atomic chemical mixture, only TiO$_2$-clusters form. Still, we believe Al$_2$O$_3$ to be the prime candidate due to substantial physical evidence in presolar grains, observations of dust around AGB stars at high temperatures, and its ability to form at high temperatures and expect the missing link to be insufficient quantitative data of Al-reactions.

In studies of the interstellar medium in galaxies, radiative transfer models of molecular emission are useful for relating molecular line observations back to the physical conditions of the gas they trace. However, doing this requires solving a highly degenerate inverse problem. In order to alleviate these degeneracies, the abundances derived from astrochemical models can be converted into column densities and fed into radiative transfer models. This enforces that the molecular gas composition used by the radiative transfer models be chemically realistic. However, because of the complexity and long running time of astrochemical models, it can be difficult to incorporate chemical models into the radiative transfer framework. In this paper, we introduce a statistical emulator of the UCLCHEM astrochemical model, built using neural networks. We then illustrate, through examples of parameter estimations, how such an emulator can be applied on real and synthetic observations.

Grain surface chemistry and its treatment in gas-grain chemical models is an area of large uncertainty. Whilst laboratory experiments are making progress, there is still much that is unknown about grain surface chemistry. Further, the results and parameters produced by experiment are often not easily translated to the rate equation approach most commonly used in astrochemical modelling. It is possible that statistical methods can reduce the uncertainty in grain surface chemical networks. In this work, a simple model of grain surface chemistry in a molecular cloud is developed and a Bayesian inference of the reactions rates is performed through MCMC sampling. Using observational data of the solid state abundances of major chemical species in molecular clouds, the posterior distributions for the rates of seven reactions producing CO, CO$_2$, CH$_3$OH and H$_2$O are calculated, in a form that is suitable for rate equation models. This represents a vital first step in the development of a method to infer reaction rates from observations of chemical abundances in astrophysical environments.

Future remote sensing of exoplanets will be enhanced by a thorough investigation of our solar system Ice Giants (Neptune-size planets). What can the configuration of the magnetic field tell us (remotely) about the interior, and what implications does that field have for the structure of the magnetosphere; energy input into the atmosphere, and surface geophysics (for example surface weathering of satellites that might harbour sub-surface oceans). How can monitoring of auroral emission help inform future remote observations of emission from exoplanets? Our Solar System provides the only laboratory in which we can perform in-situ experiments to understand exoplanet formation, dynamos, systems and magnetospheres.

We present 850 and 450 micron observations of the dense regions within the Auriga-California molecular cloud using SCUBA-2 as part of the JCMT Gould Belt Legacy Survey to identify candidate protostellar objects, measure the masses of their circumstellar material (disk and envelope), and compare the star formation to that in the Orion A molecular cloud. We identify 59 candidate protostars based on the presence of compact submillimeter emission, complementing these observations with existing Herschel/SPIRE maps. Of our candidate protostars, 24 are associated with young stellar objects (YSOs) in the Spitzer and Herschel/PACS catalogs of 166 and 60 YSOs, respectively (177 unique), confirming their protostellar nature. The remaining 35 candidate protostars are in regions, particularly around LkHalpha 101, where the background cloud emission is too bright to verify or rule out the presence of the compact 70 micron emission that is expected for a protostellar source. We keep these candidate protostars in our sample but note that they may indeed be prestellar in nature. Our observations are sensitive to the high end of the mass distribution in Auriga-Cal. We find that the disparity between the richness of infrared star forming objects in Orion A and the sparsity in Auriga-Cal extends to the submillimeter, suggesting that the relative star formation rates have not varied over the Class II lifetime and that Auriga-Cal will maintain a lower star formation efficiency.

We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta morphology is likely elongated. The difference of the CO and SiO line profiles is consistent with hydrodynamic simulations, which show that Rayleigh-Taylor instabilities cause mixing of gas, with heavier elements much more disturbed, making more elongated structure. We obtained isotopologue ratios of 28SiO/29SiO>13, 28SiO/30SiO>14, and 12CO/13CO>21, with the most likely limits of 28SiO/29SiO>128, 28SiO/30SiO>189. Low 29Si and 30Si abundances in SN 1987A are consistent with nucleosynthesis models that show inefficient formation of neutron-rich isotopes in a low metallicity environment, such as the Large Magellanic Cloud. The deduced large mass of HCO+ (~5x10^-6 Msun) and small SiS mass (<6x10^-5 Msun) might be explained by some mixing of elements immediately after the explosion. The mixing might have caused some hydrogen from the envelope to sink into carbon and oxygen-rich zones after the explosion, enabling the formation of a substantial mass of HCO+. Oxygen atoms may have penetrated into silicon and sulphur zones, suppressing formation of SiS. Our ALMA observations open up a new window to investigate chemistry, dynamics and explosive-nucleosynthesis in supernovae.

We present 450 and 850 micron submillimetre continuum observations of the IC5146 star-forming region taken as part of the JCMT Gould Belt Survey. We investigate the location of bright submillimetre (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual submillimetre clumps are identified using FellWalker and their physical properties are examined. These clumps are found to be relatively massive, ranging from 0.5to 116 MSun with a mean mass of 8 MSun and a median mass of 3.7 MSun. A stability analysis for the clumps suggest that the majority are (thermally) Jeans stable, with M/M_J < 1. We further compare the locations of known protostars with the observed submillimetre emission, finding that younger protostars, i.e., Class 0 and I sources, are strongly correlated with submillimetre peaks and that the clumps with protostars are among the most Jeans unstable. Finally, we contrast the evolutionary conditions in the two major star-forming regions within IC5146: the young cluster associated with the Cocoon Nebula and the more distributed star formation associated with the Northern Streamer filaments. The Cocoon Nebula appears to have converted a higher fraction of its mass into dense clumps and protostars, the clumps are more likely to be Jeans unstable, and a larger fraction of these remaining clumps contain embedded protostars. The Northern Streamer, however, has a larger number of clumps in total and a larger fraction of the known protostars are still embedded within these clumps.

We use a simple organism lifecycle model to explore the viability of an atmospheric habitable zone (AHZ), with temperatures that could support Earth-centric life, which sits above an environment that does not support life. To illustrate our model we use a cool Y dwarf atmosphere, such as $\mathrm{WISE~J}085510.83-0714442.5$ whose $4.5-5.2$ micron spectrum shows absorption features consistent with water vapour and clouds. We allow organisms to adapt to their atmospheric environment (described by temperature, convection, and gravity) by adopting different growth strategies that maximize their chance of survival and proliferation. We assume a constant upward vertical velocity through the AHZ. We found that the organism growth strategy is most sensitive to the magnitude of the atmospheric convection. Stronger convection supports the evolution of more massive organisms. For a purely radiative environment we find that evolved organisms have a mass that is an order of magnitude smaller than terrestrial microbes, thereby defining a dynamical constraint on the dimensions of life that an AHZ can support. Based on a previously defined statistical approach we infer that there are of order $10^9$ cool Y brown dwarfs in the Milky Way, and likely a few tens of these objects are within ten parsecs from Earth. Our work also has implications for exploring life in the atmospheres of temperate gas giants. Consideration of the habitable volumes in planetary atmospheres significantly increases the volume of habitable space in the galaxy.

We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that, if modelled as thermally-supported Bonnor-Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival $^{13}$CO velocity dispersion measurements and find that our cores are typically pressure-confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields.

Using the PACS and SPIRE spectrometers on-board the Herschel Space Observatory, we obtained spectra of two red supergiants (RSGs) in the Large Magellanic Cloud (LMC). Multiple rotational CO emission lines (J=6-5 to 15-14) and 15 H2O lines were detected from IRAS 05280-6910, and one CO line was detected from WOH G64. This is the first time CO rotational lines have been detected from evolved stars in the LMC. Their CO line intensities are as strong as those of the Galactic RSG, VY CMa. Modelling the CO lines and the spectral energy distribution results in an estimated mass-loss rate for IRAS 05280-6910 of 3x10^-4 Msun per yr. The model assumes a gas-to-dust ratio and a CO-to-H2 abundance ratio is estimated from the Galactic values scaled by the LMC metallicity ([Fe/H]~-0.3), i.e., that the CO-to-dust ratio is constant for Galactic and LMC metallicities within the uncertainties of the model. The key factor determining the CO line intensities and the mass-loss rate found to be the stellar luminosity.

We present e-MERLIN 21 cm (L-band) observations of single luminous OB stars in the Cygnus OB2 association, from the COBRaS Legacy programme. The radio observations potentially offer the most straightforward, least model-dependent, determinations of mass-loss rates, and can be used to help resolve current discrepancies in mass-loss rates via clumped and structured hot star winds. We report here that the 21 cm flux densities of O3 to O6 supergiant and giant stars are less than ~ 70 microJy. These fluxes may be translated to `smooth' wind mass-loss upper limits of ~ 4.4 - 4.8 x 10^(-6) M_sol/yr for O3 supergiants and < 2.9 x 10^(-6) M_sol/yr for B0 to B1 supergiants. The first ever resolved 21 cm detections of the hypergiant (and LBV candidate) Cyg OB2 #12 are discussed; for multiple observations separated by 14 days, we detect a ~ 69% increase in its flux density. Our constraints on the upper limits for the mass-loss rates of evolved OB stars in Cyg OB2 support the model that the inner wind region close to the stellar surface (where Halpha forms) is more clumped than the very extended geometric region sampled by our radio observations.

We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($\delta \leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $\mu$m continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only $\sim$72\% of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 $\mu$m emission structures (column densities $>3.7\times10^{21}\mathrm{\:cm}^{-2}$). The remaining 28\% of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of young stellar objects such that sources which are more evolved are associated with higher velocities.

The dust emissivity spectral index, $\beta$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $\beta$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $\beta$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, 500 $\mu$m, and 850 $\mu$m bands. Most of the derived $\beta$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $\beta$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $\beta$ variations within individual clumps and find low $\beta$ regions correlate with local temperature peaks, hinting at the possible origins of low $\beta$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape.

We present SCUBA-2 450{\mu}m and 850{\mu}m observations of the W40 complex in the Serpens-Aquila region as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) of nearby star-forming regions. We investigate radiative heating by constructing temperature maps from the ratio of SCUBA-2 fluxes using a fixed dust opacity spectral index, {\beta} = 1.8, and a beam convolution kernel to achieve a common 14.8" resolution. We identify 82 clumps ranging between 10 and 36K with a mean temperature of 20{\pm}3K. Clump temperature is strongly correlated with proximity to the external OB association and there is no evidence that the embedded protostars significantly heat the dust. We identify 31 clumps that have cores with densities greater than 105cm{^{-3}}. Thirteen of these cores contain embedded Class 0/I protostars. Many cores are associated with bright-rimmed clouds seen in Herschel 70 {\mu}m images. From JCMT HARP observations of the 12CO 3-2 line, we find contamination of the 850{\mu}m band of up to 20 per cent. We investigate the free-free contribution to SCUBA-2 bands from large-scale and ultracompact H ii regions using archival VLA data and find the contribution is limited to individual stars, accounting for 9 per cent of flux per beam at 450 {\mu}m or 12 per cent at 850 {\mu}m in these cases. We conclude that radiative heating has potentially influenced the formation of stars in the Dust Arc sub-region, favouring Jeans stable clouds in the warm east and fragmentation in the cool west.

We present a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 micron peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 x 10^23 cm^-2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 x 10^23 cm^-2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023 / 2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars.

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

The interaction of ionizing and far-ultraviolet radiation with the interstellar medium is of great importance. It results in the formation of regions in which the gas is ionized, beyond which are photodissociation regions (PDRs) in which the gas transitions to its atomic and molecular form. Several numerical codes have been implemented to study these two main phases of the interstellar medium either dynamically or chemically. In this paper we present TORUS-3DPDR, a new self-consistent code for treating the chemistry of three-dimensional photoionization and photodissociation regions. It is an integrated code coupling the two codes TORUS, a hydrodynamics and Monte Carlo radiation transport code, and 3D-PDR, a photodissociation regions code. The new code uses a Monte Carlo radiative transfer scheme to account for the propagation of the ionizing radiation including the diffusive component as well as a ray-tracing scheme based on the HEALPix package in order to account for the escape probability and column density calculations. Here, we present the numerical techniques we followed and we show the capabilities of the new code in modelling three-dimensional objects including single or multiple sources. We discuss the effects introduced by the diffusive component of the UV field in determining the thermal balance of PDRs as well as the effects introduced by a multiple sources treatment of the radiation field. We find that diffuse radiation can positively contribute to the formation of CO. With this new code, three-dimensional synthetic observations for the major cooling lines are possible, for making feasible a detailed comparison between hydrodynamical simulations and observations.

ALMA observations show a non-detection of carbon monoxide around the four most luminous asymptotic giant branch (AGB) stars in the globular cluster 47 Tucanae. Stellar evolution models and star counts show that the mass-loss rates from these stars should be ~1.2-3.5 x 10^-7 solar masses per year. We would naively expect such stars to be detectable at this distance (4.5 kpc). By modelling the ultraviolet radiation field from post-AGB stars and white dwarfs in 47 Tuc, we conclude CO should be dissociated abnormally close to the stars. We estimate that the CO envelopes will be truncated at a few hundred stellar radii from their host stars and that the line intensities are about two orders of magnitude below our current detection limits. The truncation of CO envelopes should be important for AGB stars in dense clusters. Observing the CO (3-2) and higher transitions and targeting stars far from the centres of clusters should result in the detections needed to measure the outflow velocities from these stars.

Aim: The late stages of stellar evolution are mainly governed by the mass of the stars. Low- and intermediate-mass stars lose copious amounts of mass during the asymptotic giant branch (AGB) which obscure the central star making it difficult to study the stellar spectra and determine the stellar mass. In this study, we present observational data that can be used to determine lower limits to the stellar mass. Method: Spectra of nine heavily reddened AGB stars taken by the Herschel Space Observatory display numerous molecular emission lines. The strongest emission lines are due to H2O. We search for the presence of isotopologues of H2O in these objects. Result: We detected the 16O and 17O isotopologues of water in these stars, but lines due to H2^{18}O are absent. The lack of 18O is predicted by a scenario where the star has undergone hot-bottom burning which preferentially destroys 18O relative to 16O and 17O. From stellar evolution calculations, this process is thought to occur when the stellar mass is above 5 Msun for solar metallicity. Hence, observations of different isotopologues of H2O can be used to help determine the lower limit to the initial stellar mass. Conclusion: From our observations, we deduce that these extreme OH/IR stars are intermediate-mass stars with masses of >= 5 Msun. Their high mass-loss rates of ~ 1.0e-4 Msun/yr may affect the enrichment of the interstellar medium and the overall chemical evolution of our Galaxy.

Aims. We observe the dense gas tracer CS in two nearby starburst galaxies to determine how the conditions of the dense gas varies across the circumnuclear regions in starburst galaxies. Methods. Using the IRAM-30m telescope, we mapped the distribution of the CS(2-1) and CS(3-2) lines in the circumnuclear regions of the nearby starburst galaxies NGC 3079 and NGC 6946. We also detected the formaldehyde (H2CO) and methanol (CH3OH) in both galaxies. We marginally detect the isotopologue C34S. Results. We calculate column densities under LTE conditions for CS and CH3OH. Using the detections accumulated here to guide our inputs, we link a time and depth dependent chemical model with a molecular line radiative transfer model; we reproduce the observations, showing how conditions where CS is present are likely to vary away from the galactic centres. Conclusions. Using the rotational diagram method for CH3OH, we obtain a lower limit temperature of 14 K. In addition to this, by comparing the chemical and radiative transfer models to observations, we determine the properties of the dense gas as traced by CS (and CH3OH). We also estimate the quantity of the dense gas. We find that, provided that there are a between 10^5 and 10^6 dense cores in our beam, for both target galaxies, emission of CS from warm (T = 100 - 400 K), dense (n(H2) = 10^5-6 cm-3) cores, possibly with a high cosmic ray ionisation rate (zeta = 100 zeta0) best describes conditions for our central pointing. In NGC 6946, conditions are generally cooler and/or less dense further from the centre, whereas in NGC 3079, conditions are more uniform. The inclusion of shocks allows for more efficient CS formation, leading to an order of magnitude less dense gas being required to replicate observations in some cases.

We present 850$\mu$m and 450$\mu$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks with single stellar hosts of spectral type K have higher masses than those of spectral type M, supporting the hypothesis that higher mass stars have more massive disks. Variations in disk masses calculated at the two wavelengths suggests there may be differences in dust opacity and/or dust temperature between disks with hosts of spectral types K to those with spectral type M.

In this paper we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N$_{2}$H$^{+}$ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialised. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N$_{2}$H$^{+}$ is a good tracer of the bound material of prestellar cores, although we find some evidence for N$_{2}$H$^{+}$ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C$^{18}$O and that traced by N$_{2}$H$^{+}$, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor-Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from southwest to northeast, indicating sequential star formation across the region.

We present SCUBA-2 450micron and 850micron observations of the Serpens MWC 297 region, part of the JCMT Gould Belt Survey of nearby star-forming regions. Simulations suggest that radiative feedback influences the star-formation process and we investigate observational evidence for this by constructing temperature maps. Maps are derived from the ratio of SCUBA-2 fluxes and a two component model of the JCMT beam for a fixed dust opacity spectral index of beta = 1.8. Within 40 of the B1.5Ve Herbig star MWC 297, the submillimetre fluxes are contaminated by free-free emission with a spectral index of 1.03+-0.02, consistent with an ultra-compact HII region and polar winds/jets. Contamination accounts for 73+-5 per cent and 82+-4 per cent of peak flux at 450micron and 850micron respectively. The residual thermal disk of the star is almost undetectable at these wavelengths. Young Stellar Objects are confirmed where SCUBA-2 850micron clumps identified by the fellwalker algorithm coincide with Spitzer Gould Belt Survey detections. We identify 23 objects and use Tbol to classify nine YSOs with masses 0.09 to 5.1 Msun. We find two Class 0, one Class 0/I, three Class I and three Class II sources. The mean temperature is 15+-2K for the nine YSOs and 32+-4K for the 14 starless clumps. We observe a starless clump with an abnormally high mean temperature of 46+-2K and conclude that it is radiatively heated by the star MWC 297. Jeans stability provides evidence that radiative heating by the star MWC 297 may be suppressing clump collapse.

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 micron data and improved imaging quality at 100 and 160 micron compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 micron [O I] line flux, eliminating the possibility that line contaminations distort the previously estimated dust mass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 micron flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5+-0.1 Msun of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3 Msun of amorphous carbon and 0.5 Msun of silicates, totalling 0.8 Msun of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies.

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.

The evolution of low- and intermediate-mass stars on the asymptotic giant branch (AGB) is mainly controlled by the rate at which these stars lose mass in a stellar wind. Understanding the driving mechanism and strength of the stellar winds of AGB stars and the processes enriching their surfaces with products of nucleosynthesis are paramount to constraining AGB evolution and predicting the chemical evolution of galaxies. In a previous paper we have constrained the structure of the outflowing envelope of W Hya using spectral lines of the $^{12}$CO molecule. Here we broaden this study by modelling an extensive set of H$_{2}$O and $^{28}$SiO lines observed by the three instruments on board Herschel using a state-of-the-art molecular excitation and radiative transfer code. The oxygen isotopic ratios and the $^{28}$SiO abundance profile can be connected to the initial stellar mass and to crucial aspects of dust formation at the base of the stellar wind, respectively. The modelling of H$_{2}$O and $^{28}$SiO confirms the properties of the envelope model of W Hya derived from $^{12}$CO lines. We find an H$_2$O ortho-to-para ratio of 2.5\,$^{+2.5}_{-1.0}$, consistent with what is expected for an AGB wind. The O$^{16}$/O$^{17}$ ratio indicates that W Hya has an initial mass of about 1.5 M$_\odot$. Although the ortho- and para-H$_{2}$O lines observed by HIFI appear to trace gas of slightly different physical properties, a turbulence velocity of $0.7\pm0.1$ km s$^{-1}$ fits the HIFI lines of both spin isomers and those of $^{28}$SiO well. The ortho- and para-H$_2^{16}$O and $^{28}$SiO abundances relative to H$_{2}$ are $(6^{+3}_{-2}) \times 10^{-4}$, $(3^{+2}_{-1}) \times 10^{-4}$, and $(3.3\pm 0.8)\times 10^{-5}$, respectively. Assuming a solar silicon-to-carbon ratio, the $^{28}$SiO line emission model is consistent with about one-third of the silicon atoms being locked up in dust particles.

The discovery of a hyper metal-poor star with total metallicity of $\le 10^{-5}$ Z$_\odot$, has motivated new investigations of how such objects can form from primordial gas polluted by a single supernova. In this paper we present a shock-cloud model which simulates a supernova remnant interacting with a cloud in a metal-free environment at redshift $z=10$. Pre-supernova conditions are considered, which include a multiphase neutral medium and H II region. A small dense clump ($n=100$ cm$^{-3}$), located 40 pc from a 40 M$_\odot$ metal-free star, embedded in a $n=10$ cm$^{-3}$ ambient cloud. The evolution of the supernova remnant (explosion energy $10^{52}$ erg) and its subsequent interaction with the dense clump is examined. This is the first study to include a comprehensive treatment of the non-equilibrium chemistry and associated radiative cooling that is occurring at all stages of the shock-cloud model. We have included a primordial chemistry network that covers the temperature range $10-10^9$ K, and is coupled to thermal models of atomic & molecular cooling. We find $\times10^{3}$ density enhancement of the clump (i.e maximum density $\sim 78000$ cm$^{-3}$) within this metal-free model. This is consistent with Galactic shock-cloud models considering solar metallicity gas with equilibrium cooling functions. Despite this strong compression, the cloud does not become gravitationally unstable. We find that the small cloud modelled here is destroyed for shock velocities $\gtrsim 50\,$km s$^{-1}$, and not significantly affected by shocks with velocity $\lesssim 30\,$km s$^{-1}$. Rather specific conditions are required to make such a cloud collapse.

We present a model of the photodissociation regions of NGC 4038, which is part of the Antennae galaxies. We have considered one-dimensional slabs of uniform density all having a maximum $A_V=10\,{\rm mag}$, interacting with plane-parallel radiation. The density range in our simulations spans four orders of magnitude ($100\le n\le 10^6\,{\rm cm}^{-3}$) and the UV field strength spans more than three orders of magnitude ($10\le\chi\le10^{4.5}$ multiples of the ${\rm Draine}$ field), from which we generated a grid of about 1400 simulations. We compare our results with Herschel SPIRE-FTS, CSO and ISO-LWS observations of eight CO transition lines ($J=1-0$ to $8-7$) and the [CI] 609\mu m and [OI] 146\mu m fine structure lines. We find that the molecular and atomic emission lines trace different gas components of NGC 4038, thus single emission models are insufficient to reproduce the observed values. In general, low-$J$ CO transition lines correspond to either low density regions interacting with low UV field strengths, or high density regions interacting with high UV field strengths. Higher $J$ CO transition lines are less dependent on the UV field strength and are fitted by gas with density $n\sim10^{4.5}-10^{5.2}\,{\rm cm}^{-3}$. We find that the observed fine structure line ratio of [CI] 609\mu m/[OI] 146\mu m is reproduced by clouds subject to weaker UV fields compared to the CO lines. We make estimates of the $X_{\rm CO}$ factor which relates the CO emission with the column density of molecular hydrogen, and find that it is less than the canonical Milky Way value.

Asymptotic giant branch (AGB) stars lose their envelopes by means of a stellar wind whose driving mechanism is not understood well. Characterizing the composition and thermal and dynamical structure of the outflow provides constraints that are essential for understanding AGB evolution, including the rate of mass loss and isotopic ratios. We modeled the CO emission from the wind of the low mass-loss rate oxygen-rich AGB star W Hya using data obtained by the HIFI, PACS, and SPIRE instruments onboard the Herschel Space Observatory and ground-based telescopes. $^{12}$CO and $^{13}$CO lines are used to constrain the intrinsic $^{12}$C/$^{13}$C ratio from resolved HIFI lines. The acceleration of the outflow up to about 5.5 km/s is quite slow and can be represented by a beta-type velocity law with index 5. Beyond this point, acceleration up the terminal velocity of 7 km/s is faster. Using the J=10-9, 9-8, and 6-5 transitions, we find an intrinsic $^{12}$C/$^{13}$C ratio of $18\pm10$ for W Hya, where the error bar is mostly due to uncertainties in the $^{12}$CO abundance and the stellar flux around 4.6 $\mu$m. To match the low-excitation CO lines, these molecules need to be photo-dissociated at about 500 stellar radii. The radial dust emission intensity profile measured by PACS images at 70 $\mu$m shows substantially stronger emission than our model predicts beyond 20 arcsec. The initial slow acceleration of the wind implies inefficient wind driving in the lower part of the envelope. The final injection of momentum in the wind might be the result of an increase in the opacity thanks to the late condensation of dust species. The derived intrinsic isotopologue ratio for W Hya is consistent with values set by the first dredge-up and suggestive of an initial mass of 2 M$_\odot$ or more. However, the uncertainty in the main-sequence mass derived based on this isotopologic ratio is large.

Previously, we have presented the first study to investigate the response of the Jovian thermosphere to transient variations in solar wind dynamic pressure, using a coupled, azimuthally symmetric global circulation model coupled with a simple magnetosphere model. This work (Yates et al., 2013, submitted) described the response of thermospheric flows, momentum sources, and the magnetosphere-ionosphere coupling currents to transient compressions and expansions in the magnetosphere. The present study describes the response of thermospheric heating, cooling and the auroral emissions to the aforementioned transient events. We find that transient compressions and expansions, on time scales <= 3 hours, cause at least a factor of two increase in Joule heating per unit volume. Ion drag significantly changes the kinetic energy of the thermospheric neutrals depending on whether the magnetosphere is compressed or expanded. These processes lead to local temperature variations >= 25 K and a ~2000 TW increase in the total power dissipated in the thermosphere. In terms of auroral processes, transient compressions increase main oval UV emission by a factor of ~4.5 whilst transient expansions increase this main emission by a more modest 37%. Both types of transient event cause shifts in the position of the main oval, of up to 1 deg latitude.

The importance of the Jovian thermosphere with regard to magnetosphere-ionosphere coupling is often neglected in magnetospheric physics. We present the first study to investigate the response of the Jovian thermosphere to transient variations in solar wind dynamic pressure, using an azimuthally symmetric global circulation model coupled to a simple magnetosphere and fixed auroral conductivity model. In our simulations, the Jovian magnetosphere encounters a solar wind shock or rarefaction region and is subsequently compressed or expanded. We present the ensuing response of the coupling currents, thermospheric flows, heating and cooling terms, and the aurora to these transient events. Transient compressions cause the reversal, with respect to steady state, of magnetosphere-ionosphere coupling currents and momentum transfer between the thermosphere and magnetosphere. They also cause at least a factor of two increase in the Joule heating rate. Ion drag significantly changes the kinetic energy of the thermospheric neutrals depending on whether the magnetosphere is compressed or expanded. Local temperature variations appear between ~ -45 and 175 K for the compression scenario and ~ -20 and 50 K for the expansion case. Extended regions of equatorward flow develop in the wake of compression events - we discuss the implications of this behaviour for global energy transport. Both compressions and expansions lead to a ~ 2000 TW increase in the total power dissipated or deposited in the thermosphere. In terms of auroral processes, transient compressions increase main oval UV emission by a factor of 4.5 whilst transient expansions increase this main emission by a more modest 37%. Both types of transient event cause shifts in the position of the main oval, of up to 1 deg latitude.

We present an analysis of the far-infrared and submillimetre molecular emission line spectrum of the luminous M-supergiant VY CMa, observed with the SPIRE and PACS spectrometers aboard the Herschel Space Observatory. Over 260 emission lines were detected in the 190-650-micron SPIRE FTS spectra, with one-third of the observed lines being attributable to H2O. Other detected species include CO, 13CO, H2^18O, SiO, HCN, SO, SO2, CS, H2S, and NH3. Our model fits to the observed 12CO and 13CO line intensities yield a 12C/13C ratio of 5.6+-1.8, consistent with measurements of this ratio for other M supergiants, but significantly lower than previously estimated for VY CMa from observations of lower-J lines. The spectral line energy distribution for twenty SiO rotational lines shows two temperature components: a hot component at 1000 K, which we attribute to the stellar atmosphere and inner wind, plus a cooler ~200 K component, which we attribute to an origin in the outer circumstellar envelope. We fit the line fluxes of 12CO, 13CO, H2O and SiO, using the SMMOL non-LTE line transfer code, with a mass-loss rate of 1.85x10^-4 Msun yr^-1 between 9 R* and 350 R*. To fit the observed line fluxes of 12CO, 13CO, H2O and SiO with SMMOL non-LTE line radiative transfer code, along with a mass-loss rate of 1.85x10^-4 Msun yr^-1. To fit the high rotational lines of CO and H2O, the model required a rather flat temperature distribution inside the dust condensation radius, attributed to the high H2O opacity. Beyond the dust condensation radius the gas temperature is fitted best by an r^-0.5 radial dependence, consistent with the coolant lines becoming optically thin. Our H2O emission line fits are consistent with an ortho:para ratio of 3 in the outflow.

We report spectroscopic and imaging observations of rotational transitions of cold CO and SiO in the ejecta of SN1987A, the first such emission detected in a supernova remnant. In addition to line luminosities for the CO J=1-0, 2-1, 6-5, and 7-6 transitions, we present upper limits for all other transitions up to J=13-12, collectively measured from the Atacama Large Millimeter Array (ALMA), the Atacama Pathfinder EXperiment (APEX), and the Herschel Spectral and Photometric Imaging REceiver (SPIRE). Simple models show the lines are emitted from at least 0.01 solar masses of CO at a temperature > 14 K, confined within at most 35% of a spherical volume expanding at ~ 2000 km/s. Moreover, we locate the emission within 1'' of the central debris. These observations, along with a partial observation of SiO, confirm the presence of cold molecular gas within supernova remnants and provide insight into the physical conditions and chemical processes in the ejecta. Furthermore, we demonstrate the powerful new window into supernova ejecta offered by submillimeter observations.

Aim : In order to study the history of mass loss in extreme OH/IR stars, we observed a number of these objects using CO as a tracer of the density and temperature structure of their circumstellar envelopes. Method : Combining CO observations from the Herschel Space Observatory with those from the ground, we trace mass loss rates as a function of radius in five extreme OH/IR stars. Using radiative transfer modelling, we modelled the dusty envelope as well as the CO emission. The high-rotational transitions of CO indicate that they originate in a dense superwind region close to the star while the lower transitions tend to come from a more tenuous outer wind which is a result of the mass loss since the early AGB phase. Result : The models of the circumstellar envelopes around these stars suggest that they have entered a superwind phase in the past 200 - 500 years. The low 18O/17O (~ 0.1 compared to the solar abundance ratio of ~ 5) and 12C/13C (3-30 cf. the solar value of 89) ratios derived from our study support the idea that these objects have undergone hot-bottom burning and hence that they are massive M >= 5 solar-mass AGB stars.

Photodissociation regions (PDRs) define the transition zone between an ionized and a dark molecular region. They consist of neutral gas which interacts with far-ultraviolet radiation and are characterized by strong infrared line emission. Various numerical codes treating one-dimensional PDRs have been developed in the past, simulating the complexity of chemical reactions occurring and providing a better understanding of the structure of a PDR. In this paper we present the three-dimensional code, 3D-PDR, which can treat PDRs of arbitrary density distribution. The code solves the chemistry and the thermal balance self-consistently within a given three-dimensional cloud. It calculates the total heating and cooling functions at any point in a given PDR by adopting an escape probability method. It uses a HEALPix-based ray-tracing scheme to evaluate the attenuation of the far-ultraviolet radiation in the PDR and the propagation of the far-infrared/submm line emission out of the PDR. We present benchmarking results and apply 3D-PDR to i) a uniform-density spherical cloud interacting with a plane-parallel external radiation field, ii) a uniform-density spherical cloud interacting with a two-component external radiation field, and iii) a cometary globule interacting with a plane-parallel external radiation field. We find that the code is able to reproduce the benchmarking results of various other one-dimensional numerical codes treating PDRs. We also find that the accurate treatment of the radiation field in the fully three-dimensional treatment of PDRs can in some cases leads to different results when compared to a standard one-dimensional treatment.

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.

CW Leo has been observed six times between October 2009 and June 2012 with the SPIRE instrument on board the Herschel satellite. Variability has been detected in the flux emitted by the central star with a period of 639 \pm 4 days, in good agreement with determinations in the literature. Variability is also detected in the bow shock around CW Leo that had previously been detected in the ultraviolet and Herschel PACS/SPIRE data. Although difficult to prove directly, our working hypothesis is that this variability is directly related to that of the central star. In this case, fitting a sine curve with the period fixed to 639 days results in a time-lag in the variability between bow shock and the central star of 402 \pm 37 days. The orientation of the bow shock relative to the plane of the sky is unknown (but see below). For an inclination angle of zero degrees, the observed time-lag translates into a distance to CW Leo of 130 \pm 13 pc, and for non-zero inclination angles the distance is smaller. Fitting the shape of the bow shock with an analytical model (Wilkin 1996), the effect of the inclination angle on the distance may be estimated. Making the additional assumption that the relative peculiar velocity between the interstellar medium (ISM) and CW Leo is determined entirely by the star space velocity with respect to the local standard of rest (i.e. a stationary ISM), the inclination angle is found to be (-33.3 \pm 0.8) degrees based on the observed proper motion and radial velocity. Using the Wilkin model, our current best estimate of the distance to CW Leo is 123 \pm 14 pc. For a distance of 123 pc, we derive a mean luminosity of 7790 \pm 150 Lsol (internal error).

The Gould Belt Legacy Survey will map star-forming regions within 500 pc, using HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common-User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope (JCMT). This paper describes HARP observations of the J = 3-2 transitions of 13CO and C18O towards Orion A. The 1500-resolution observations cover 5 pc of the Orion filament, including OMC1 (inc. BN-KL and Orion Bar), OMC 2/3 and OMC 4, and allow a comparative study of the molecular gas properties throughout the star-forming cloud. The filament shows a velocity gradient of ~1 km/s /pc between OMC 1, 2 and 3, and high velocity emission is detected in both isotopologues. The Orion Nebula and Bar have the largest masses and line widths, and dominate the mass and energetics of the high velocity material. Compact, spatially resolved emission from CH3CN, 13CH3OH, SO, HCOOCH3, C2H5OH, CH3CHO and CH3OCHO is detected towards the Orion Hot Core. The cloud is warm, with a median excitation temperature of ~24 K; the Orion Bar has the highest excitation temperature gas, at >80 K. The C18O excitation temperature correlates well with the dust temperature (to within 40%). The C18O emission is optically thin, and the 13CO emission is marginally optically thick; despite its high mass, OMC 1 shows the lowest opacities. A virial analysis indicates that Orion A is too massive for thermal or turbulent support, but is consistent with a model of a filamentary cloud that is threaded by helical magnetic fields. The variation of physical conditions across the cloud is reflected in the physical characteristics of the dust cores....continued

We investigate the origin and nature of the profiles of water and ammonia observ ed toward the L1157-B1 clump as part of the HIFI CHESS survey (Ceccarelli et al. 2010) using a new code coupling a gas-grain chemical model with a parametric shock model. Fir st results from the unbiased survey (Lefloch et al. 2010, Codella et al. 2010) reveal different molecular components at different excitation conditions coexisting in the B1 bow shock structure, with NH$_3$, H$_2$CO and CH$_3$OH emitting only at relatively low outflow velocities whereas H$_2$O shows bright emission at high velocities. Our model suggests that these differences are purely chemical and can be explained by the presence of a C-type shock whose maximum temperature must be close to 4000 K along the B1 clump.

The coupling of Jupiter's magnetosphere and ionosphere plays a vital role in creating its auroral emissions. The strength of these emissions is dependent on the difference in speed of the rotational flows within Jupiter's high-latitude thermosphere and the planet's magnetodisc. Using an azimuthally symmetric global circulation model, we have simulated how upstream solar wind conditions affect the energy and direction of atmospheric flows. In order to simulate the effect of a varying dynamic pressure in the upstream solar wind, we calculated three magnetic field profiles representing compressed, averaged and expanded `middle' magnetospheres. These profiles were then used to solve for the angular velocity of plasma in the magnetosphere. This angular velocity determines the strength of currents flowing between the ionosphere and magnetosphere. We examine the influence of variability in this current system upon the global winds and energy inputs within the Jovian thermosphere. We find that the power dissipated by Joule heating and ion drag increases by ~190% and ~185% from our compressed to expanded model respectively. We investigated the effect of exterior boundary conditions on our models and found that by reducing the radial current at the outer edge of the magnetodisc, we also limit the thermosphere's ability to transmit angular momentum to this region.

In 2001, the discovery of circumstellar water vapour around the ageing carbon star IRC+10216 was announced. This detection challenged the current understanding of chemistry in old stars, since water vapour was predicted to be absent in carbon-rich stars. Several explanations for the occurrence of water vapour were postulated, including the vaporization of icy bodies (comets or dwarf planets) in orbit around the star, grain surface reactions, and photochemistry in the outer circumstellar envelope. However, the only water line detected so far from one carbon-rich evolved star can not discriminate, by itself, between the different mechanisms proposed. Here we report on the detection by the Herschel satellite of dozens of water vapour lines in the far-infrared and sub-millimetre spectrum of IRC+10216, including some high-excitation lines with energies corresponding to ~1000 K. The emission of these high-excitation water lines can only be explained if water vapour is present in the warm inner region of the envelope. A plausible explanation for the formation of warm water vapour appears to be the penetration of ultraviolet (UV) photons deep into a clumpy circumstellar envelope. This mechanism triggers also the formation of other molecules such as ammonia, whose observed abundances are much higher than hitherto predicted.

We present a theoretical study of CS line profiles in archetypal hot cores. We provide estimates of line fluxes from the CS(1-0) to the CS(15-14) transitions and present the temporal variation of these fluxes. We find that \textit{i)} the CS(1-0) transition is a better tracer of the Envelope of the hot core whereas the higher-J CS lines trace the ultra-compact core; \textit{ii)} the peak temperature of the CS transitions is a good indicator of the temperature inside the hot core; \textit{iii)} in the Envelope, the older the hot core the stronger the self-absorption of CS; \textit{iv)} the fractional abundance of CS is highest in the innermost parts of the ultra-compact core, confirming the CS molecule as one of the best tracers of very dense gas.

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 Gould Belt Legacy Survey (GBS) on the JCMT has observed a region of 260 square arcminutes in 12CO J=3--2 emission, and a 190 square arcminute subset of this in 13CO and C18O towards the Serpens molecular cloud. We examine the global velocity structure of the non-outflowing gas, and calculate excitation temperatures and opacities. The large scale mass and energetics of the region are evaluated, with special consideration for high velocity gas. We find the cloud to have a mass of 203 solar masses, and to be gravitationally bound, and that the kinetic energy of the outflowing gas is approximately seventy percent of the turbulent kinetic energy of the cloud. We identify compact outflows towards some of the submillimetre Class 0/I sources in the region

The interstellar medium is enriched primarily by matter ejected from evolved low and intermediate mass stars. The outflows from these stars create a circumstellar envelope in which a rich gas-phase and dust-nucleation chemistry takes place. We observed the nearest carbon-rich evolved star, IRC+10216, using the PACS (55-210 {\mu}m) and SPIRE (194-672 {\mu}m) spectrometers on board Herschel. We find several tens of lines from SiS and SiO, including lines from the v=1 vibrational level. For SiS these transitions range up to J=124-123, corresponding to energies around 6700K, while the highest detectable transition is J=90-89 for SiO, which corresponds to an energy around 8400K. Both species trace the dust formation zone of IRC+10216, and the broad energy ranges involved in their detected transitions permit us to derive the physical properties of the gas and the particular zone in which each species has been formed. This allows us to check the accuracy of chemical thermodynamical equilibrium models and the suggested depletion of SiS and SiO due to accretion onto dust grains.

We present far-infrared and submillimetre spectra of three carbon-rich evolved objects, AFGL 2688, AFGL 618 and NGC 7027. The spectra were obtained with the SPIRE Fourier transform spectrometer on board the Herschel Space Observatory, and cover wavelengths from 195-670 um, a region of the electromagnetic spectrum hitherto difficult to study in detail. The far infrared spectra of these objects are rich and complex, and we measure over 150 lines in each object. Lines due to 18 different species are detected. We determine physical conditions from observations of the rotational lines of several molecules, and present initial large velocity gradient models for AFGL 618. We detect water in AFGL 2688 for the first time, and confirm its presence in AFGL 618 in both ortho and para forms. In addition, we report the detection of the J=1-0 line of CH+ in NGC 7027.

With a luminosity > 10^5 Lsun and a mass-loss rate of about 2.10-4 Msun/yr, the red supergiant VY CMa truly is a spectacular object. Because of its extreme evolutionary state, it could explode as supernova any time. Studying its circumstellar material, into which the supernova blast will run, provides interesting constraints on supernova explosions and on the rich chemistry taking place in such complex circumstellar envelopes. We have obtained spectroscopy of VYCMa over the full wavelength range offered by the PACS and SPIRE instruments of Herschel, i.e. 55 to 672 micron. The observations show the spectral fingerprints of more than 900 spectral lines, of which more than half belong to water. In total, we have identified 13 different molecules and some of their isotopologues. A first analysis shows that water is abundantly present, with an ortho-to-para ratio as low as 1.3:1, and that chemical non-equilibrium processes determine the abundance fractions in the inner envelope.

We present chemical and radiative transfer models for the many far-IR ortho- and para-H2O lines that were observed from the Orion-KL region in high resolution Fabry-Perot (FP) mode by the Long Wavelength Spectrometer (LWS) on board the Infrared Space Observatory (ISO). The chemistry of the region was first studied by simulating the conditions in the different known components of Orion-KL: chemical models for a hot core, a plateau and a ridge were coupled with an accelerated Lambda$-iteration (ALI) radiative transfer model to predict H2O line fluxes and profiles. Our models include the first 45 energy levels of ortho- and para-H2O. We find that lines arising from energy levels below 560 K were best reproduced by a gas of density 3x10^5 cm^-3 at a temperature of 70-90 K, expanding at a velocity of 30 km s^-1 and with a H2O/H2 abundance ratio of the order of 2 - 3 x 10^-5, similar to the abundance derived by Cernicharo et al. (2006). However, the model that best reproduces the fluxes and profiles of H2O lines arising from energy levels above 560 K has a significantly higher H2O/H2 abundance, 1 - 5 x 10^-4, originating from gas of similar density, in the Plateau region, that has been heated to 300 K, relaxing to 90-100 K. We conclude that the observed water lines do not originate from high temperature shocks.

We have conducted a Galactic plane survey of methanol masers at 6668 MHz using a 7-beam receiver on the Parkes telescope. Here we present results from the first part, which provides sensitive unbiased coverage of a large region around the Galactic Centre. Details are given for 183 methanol maser sites in the longitude range 345$^{\circ}$ through the Galactic Centre to 6$^{\circ}$. Within 6$^{\circ}$ of the Centre, we found 88 maser sites, of which more than half (48) are new discoveries. The masers are confined to a narrow Galactic latitude range, indicative of many sources at the Galactic Centre distance and beyond, and confined to a thin disk population; there is no high latitude population that might be ascribed to the Galactic Bulge. Within 2$^{\circ}$ of the Galactic Centre the maser velocities all lie between -60 and +77 \kms, a range much smaller than the 540 \kms range observed in CO. Elsewhere, the maser with highest positive velocity (+107 \kms) occurs, surprisingly, near longitude 355$^{\circ}$ and is probably attributable to the Galactic bar. The maser with the most negative velocity (-127 \kms) is near longitude 346$^{\circ}$, within the longitude-velocity locus of the near side of the `3-kpc arm'. It has the most extreme velocity of a clear population of masers associated with the near and far sides of the 3-kpc arm. Closer to the Galactic Centre the maser space density is generally low, except within 0.25 kpc of the Centre itself, the `Galactic Centre Zone', where it is 50 times higher, which is hinted at by the longitude distribution, and confirmed by the unusual velocities.

As part of a JCMT Legacy Survey of star formation in the Gould Belt, we present early science results for Taurus. CO J=3-2 maps have been secured along the north-west ridge and bowl, collectively known as L 1495, along with deep 13CO and C18O J=3-2 maps in two sub-regions. With these data we search for molecular outflows, and use the distribution of flows, HH objects and shocked H2 line emission features, together with the population of young stars, protostellar cores and starless condensations to map star formation across this extensive region. In total 21 outflows are identified. It is clear that the bowl is more evolved than the ridge, harbouring a greater population of T Tauri stars and a more diffuse, more turbulent ambient medium. By comparison, the ridge contains a much younger, less widely distributed population of protostars which, in turn, is associated with a greater number of molecular outflows. We estimate the ratio of the numbers of prestellar to protostellar cores in L 1495 to be ~ 1.3-2.3, and of gravitationally unbound starless cores to (gravitationally bound) prestellar cores to be ~ 1. If we take previous estimates of the protostellar lifetime of ~ 5 x 10^5 yrs, this indicates a prestellar lifetime of 9(+/-3) x 10^5 yrs. From the number of outflows we also crudely estimate the star formation efficiency in L 1495, finding it to be compatible with a canonical value of 10-15 %. We note that molecular outflow-driving sources have redder near-IR colours than their HH jet-driving counterparts. We also find that the smaller, denser cores are associated with the more massive outflows, as one might expect if mass build-up in the flow increases with the collapse and contraction of the protostellar envelope.

Far-infrared [C ii] 158 micron, [O i] 145 micron and [O i] 63 micron fine structure emission line fluxes were measured from archival Infrared Space Observatory Long Wavelength Spectrometer spectra of 46 extragalactic sources, with 28 sources providing detections in all three lines. For 12 of the sources, the contribution to the [C ii] 158 micron line flux from H ii regions could be estimated from their detected [N ii] 122 micron line fluxes. The measured [C ii]/[O i] and [O i] 63/145 micron line flux ratios were compared with those from a grid of PDR models previously computed using the UCL PDR code. Persistent offsets between the observed and modelled line ratios could be partly attributed to the effects of [O i] 63 micron self-absorption. Using the SMMOL code, we calculated model [O i] line profiles and found that the strength of the [O i] 63 micron line was reduced by 20-80%, depending on the PDR parameters. We conclude that high PDR densities and radiation field strengths, coupled with the effects of [O i] 63 micron self-absorption, are likely to provide the best match to the observed line flux ratios.

The Gould Belt Legacy Survey will survey nearby star-forming regions (within 500 pc), using HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common- User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope (JCMT). This paper describes the initial data obtained using HARP to observe 12CO, 13CO and C18O J = 3 - 2 towards two regions in Orion B, NGC 2024 and NGC 2071. We describe the physical characteristics of the two clouds, calculating temperatures and opacities utilizing all three isotopologues. We find good agreement between temperatures calculated from CO and from dust emission in the dense, energetic regions. We determine the mass and energetics of the clouds, and of the high-velocity material seen in 12CO emission, and compare the relative energetics of the high- and low-velocity material in the two clouds. We present a CLUMPFIND analysis of the 13CO condensations. The slope of the condensation mass functions, at the high-mass ends, is similar to the slope of the initial mass function.

A new 7-beam 6-7 GHz receiver has been built to survey the Galaxy and the Magellanic Clouds for newly forming high-mass stars that are pinpointed by strong methanol maser emission at 6668 MHz. The receiver was jointly constructed by Jodrell Bank Observatory (JBO) and the Australia Telescope National Facility (ATNF) and allows simultaneous coverage at 6668 and 6035 MHz. It was successfully commissioned at Parkes in January 2006 and is now being used to conduct the Parkes-Jodrell multibeam maser survey of the Milky Way. This will be the first systematic survey of the entire Galactic plane for masers of not only 6668-MHz methanol, but also 6035-MHz excited-state hydroxyl. The survey is two orders of magnitude faster than most previous systematic surveys and has an rms noise level of ~0.17 Jy.This paper describes the observational strategy, techniques and reduction procedures of the Galactic and Magellanic Cloud surveys, together with deeper, pointed, follow-up observations and complementary observations with other instruments. It also includes an estimate of the survey detection efficiency. The 111 days of observations with the Parkes telescope have so far yielded >800 methanol sources, of which ~350 are new discoveries. The whole project will provide the first comprehensive Galaxy-wide catalogue of 6668-MHz and 6035-MHz masers.

The evolution of star-forming core analogues undergoing inside-out collapse is studied with a multi-point chemodynamical model which self-consistently computes the abundance distribution of chemical species in the core. For several collapse periods the output chemistry of infall tracer species such as HCO+, CS, and N2H+, is then coupled to an accelerated Lambda-iteration radiative transfer code, which predicts the emerging molecular line profiles using two different input gas/dust temperature distributions. We investigate the sensitivity of the predicted spectral line profiles and line asymmetry ratios to the core temperature distribution, the time-dependent model chemistry, as well as to ad hoc abundance distributions. The line asymmetry is found to be strongly dependent on the adopted chemical abundance distribution. In general, models with a warm central region show higher values of blue asymmetry in optically thick HCO+ and CS lines than models with a starless core temperature profile. We find that in the formal context of Shu-type inside-out infall, and in the absence of rotation or outflows, the relative blue asymmetry of certain HCO+ and CS transitions is a function of time and, subject to the foregoing caveats, can act as a collapse chronometer. The sensitivity of simulated HCO+ line profiles to linear radial variations, subsonic or supersonic, of the internal turbulence field is investigated in the separate case of static cores.

As part of the first high resolution far-IR spectral survey of the Orion KL region (Lerate et al. 2006), we observed 20 CO emission lines with Jup=16 to Jup=39 (upper levels from approx 752 K to 4294 K above the ground state). Observations were taken using the Long Wavelength Spectrometer (LWS) on board the Infrared Space Observatory (ISO), in its high resolution Fabry-Perot (FP) mode (approx 33 km s$^{-1}$). We present here an analysis of the final calibrated CO data, performed with a more sophisticated modelling technique than hitherto, including a detailed analysis of the chemistry, and discuss similarities and differences with previous results. The inclusion of chemical modelling implies that atomic and molecular abundances are time-predicted by the chemistry. This provides one of the main differences with previous studies in which chemical abundances needed to be assumed as initial condition. The chemistry of the region is studied by simulating the conditions of the different known components of the KL region: chemical models for a hot core, a plateau and a ridge are coupled with an accelerated Lambda-iteration (ALI)radiative transfer model to predict line fluxes and profiles. We conclude that the CO transitions with 18<Jup<25 mainly arise from a hot core of diameter 0.02 pc and a density of 10$^{7}$ cm$^{-3}$ rather from the plateau as previous studies had indicated.

We present the results of the first complete survey of the Large and Small Magellanic Clouds for 6668-MHz methanol and 6035-MHz excited-state hydroxyl masers. In addition to the survey, higher-sensitivity targeted searches towards known star-formation regions were conducted. The observations yielded the discovery of a fourth 6668-MHz methanol maser in the Large Magellanic Cloud (LMC), found towards the star-forming region N160a, and a second 6035-MHz excited-state hydroxyl maser, found towards N157a. We have also re-observed the three previously known 6668-MHz methanol masers and the single 6035-MHz hydroxyl maser. We failed to detect emission from either transition in the Small Magellanic Cloud. All observations were initially made using the Methanol Multibeam (MMB) survey receiver on the 64-m Parkes telescope as part of the MMB project and accurate positions have been measured with the Australia Telescope Compact Array (ATCA). We compare the maser populations in the Magellanic Clouds with those of our Galaxy and discuss their implications for the relative rates of massive star-formation, heavy metal abundance, and the abundance of complex molecules. The LMC maser populations are demonstrated to be smaller than their Milky Way counterparts. Methanol masers are under-abundant by a factor of ~45, whilst hydroxyl and water masers are a factor of ~10 less abundant than our Galaxy.

We present the scientific motivation and observing plan for an upcoming detection survey for debris disks using the James Clerk Maxwell Telescope. The SCUBA-2 Unbiased Nearby Stars (SUNS) Survey will observe 500 nearby main sequence and sub-giant stars (100 of each of the A, F, G, K and M spectral classes) to the 850 micron extragalactic confusion limit to search for evidence of submillimeter excess, an indication of circumstellar material. The survey distance boundaries are 8.6, 16.5, 22, 25 and 45 pc for M, K, G, F and A stars, respectively, and all targets lie between the declinations of -40 deg to 80 deg. In this survey, no star will be rejected based on its inherent properties: binarity, presence of planetary companions, spectral type or age. This will be the first unbiased survey for debris disks since IRAS. We expect to detect ~125 debris disks, including ~50 cold disks not detectable in current shorter wavelength surveys. A substantial amount of complementary data will be required to constrain the temperatures and masses of discovered disks. High resolution studies will likely be required to resolve many of the disks. Therefore, these systems will be the focus of future observational studies using a variety of observatories to characterize their physical properties. For non-detected systems, this survey will set constraints (upper limits) on the amount of circumstellar dust, of typically 200 times the Kuiper Belt mass, but as low as 10 times the Kuiper Belt mass for the nearest stars in the sample (approximately 2 pc).

This paper describes a James Clerk Maxwell Telescope (JCMT) legacy survey that has been awarded roughly 500 hrs of observing time to be carried out from 2007 to 2009. In this survey we will map with SCUBA-2 (Submillimetre Common User Bolometer Array 2) almost all of the well-known low-mass and intermediate-mass star-forming regions within 0.5 kpc that are accessible from the JCMT. Most of these locations are associated with the Gould Belt. From these observations we will produce a flux-limited snapshot of star formation near the Sun, providing a legacy of images, as well as point-source and extended-source catalogues, over almost 700 square degrees of sky. The resulting images will yield the first catalogue of prestellar and protostellar sources selected by submillimetre continuum emission, and should increase the number of known sources by more than an order of magnitude. We will also obtain CO maps with the array receiver HARP (Heterodyne Array Receiver Programme), in three CO isotopologues, of a large typical sample of prestellar and protostellar sources. We will then map the brightest hundred sources with the SCUBA-2 polarimeter (POL-2), producing the first statistically significant set of polarization maps in the submillimetre. The images and source catalogues will be a powerful reference set for astronomers, providing a detailed legacy archive for future telescopes, including ALMA, Herschel and JWST.

Stars form in the densest, coldest, most quiescent regions of molecular clouds. Molecules provide the only probes which can reveal the dynamics, physics, chemistry and evolution of these regions, but our understanding of the molecular inventory of sources and how this is related to their physical state and evolution is rudimentary and incomplete. The Spectral Legacy Survey (SLS) is one of seven surveys recently approved by the JCMT Board. Starting in 2007, the SLS will produce a spectral imaging survey of the content and distribution of all the molecules detected in the 345 GHz atmospheric window (between 332 GHz and 373 GHz) towards a sample of 5 sources. Our intended targets are: a low mass core (NGC1333 IRAS4), 3 high mass cores spanning a range of star forming environments and evolutionary states (W49, AFGL2591, and IRAS20126), and a PDR (the Orion Bar). The SLS will use the unique spectral imaging capabilities of HARP-B/ACSIS to study the molecular inventory and the physical structure of these objects, which span different evolutionary stages and physical environments, to probe their evolution during the star formation process. As its name suggests, the SLS will provide a lasting data legacy from the JCMT that is intended to benefit the entire astronomical community. As such, the entire data set (including calibrated spectral datacubes, maps of molecular emission, line identifications, and calculations of the gas temperature and column density) will be publicly available.

We present the detection and modeling of more than 70 far-IR pure rotational lines of water vapor, including the 18O and 17O isotopologues, towards Orion KL. Observations were performed with the Long Wavelength Spectrometer Fabry-Perot (LWS/FP; R~6800-9700) on board the Infrared Space Observatory (ISO) between ~43 and ~197 um. The water line profiles evolve from P-Cygni type profiles (even for the H2O18 lines) to pure emission at wavelengths above ~100 um. We find that most of the water emission/absorption arises from an extended flow of gas expanding at 25+-5 kms^-1. Non-local radiative transfer models show that much of the water excitation and line profile formation is driven by the dust continuum emission. The derived beam averaged water abundance is 2-3x10^-5. The inferred gas temperature Tk=80-100 K suggests that: (i) water could have been formed in the "plateau" by gas phase neutral-neutral reactions with activation barriers if the gas was previously heated (e.g. by shocks) to >500 K and/or (ii) H2O formation in the outflow is dominated by in-situ evaporation of grain water-ice mantles and/or (iii) H2O was formed in the innermost and warmer regions (e.g. the hot core) and was swept up in ~1000 yr, the dynamical timescale of the outflow.

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.

We have carried out a high spectral resolution line survey towards the Orion Kleinmann-Low (KL) cluster from 44-188 um. The observations were taken with the Long Wavelength Spectrometer (LWS) in Fabry-Perot mode, on board the Infrared Space Observatory (ISO). A total of 152 lines are clearly detected and a further 34 features are present as possible detections. The spectrum is dominated by the molecular species H2O, OH and CO, along with [OI] and [CII] lines from PDR or shocked gas and [OIII], [NIII] lines from the foreground M42 HII region. Several isotopic species, as well as NH3, are also detected. HDO and H3O+ are tentatively detected for the first time in the far-infrared range towards Orion-KL. A basic analysis of the line observations is carried out, by comparing with previous measurements and published models and deriving rotational temperatures and column densities in the case of the molecular species. The complexity of the region requires more sophisticated models for the interpretation of all the line observations.

In order to investigate the origin and the structure of the low velocity, chemically rich clumps observed along the lobes of low- and intermediate-mass outflows, we construct a detailed model of the S1 clump along the CB3 outflow. We use a time-dependent chemical model coupled with a radiative transfer model to reproduce the observed line profile for a direct comparison with previous observations of this clump. We find that the simultaneous fitting of multiple species and transitions is a powerful tool in constraining the physical parameters of the gas. Different scenarios for the clump formation have been investigated. The models that better reproduce all the observed lines are those where the clump is formed, at least partially, before the advent of the outflow; with the advent of the outflow the clump undergoes a short period of non-dissociative shock and the consequent release of the icy mantle together with the high temperature chemistry leads to the observed chemical enrichment. Our results also suggest the presence of substructure within the clump: a more extended component traced by CS, SO and the lower energy transitions (3_k-2_k and 2_k-1_k) of CH_3OH, and a more compact component traced by H_2CO, SO_2 and the higher energy transitions (5_k-4_k) of CH_3OH.

As part of the first far-IR line survey towards Orion KL, we present the detection of seven new rotationally excited OH Lambda-doublets (at 48, 65, 71, 79, 98 and 115 um). Observations were performed with the Long Wavelength Spectrometer (LWS) Fabry-Perots on board the Infrared Space Observatory (ISO). In total, more than 20 resolved OH rotational lines, with upper energy levels up to 620 K, have been detected at an angular and velocity resolutions of 80$'' and 33 km s^-1 respectively. OH line profiles show a complex behavior evolving from pure absorption, P-Cygni type to pure emission. We also present a large scale 6' declination raster in the OH ^2\Pi_3/2 J=5/2^+-3/2^- and ^2\Pi_3/2 J=7/2^-5/2^+ lines (at 119.441 and 84.597 um) revealing the decrease of excitation outside the core of the cloud. From the observed profiles, mean intrinsic line widths and velocity offsets between emission and absorption line peaks we conclude that most of the excited OH arises from Orion outflow(s), i.e. the ``plateau'' component. We determine an averaged OH abundance relative to H_2 of X(OH)=(0.5-1.0)x10^-6, a kinetic temperature of 100 K and a density of n(H_2)=5x10^5 cm^-3. Even with these conditions, the OH excitation is heavily coupled with the strong dust continuum emission from the inner hot core regions and from the expanding flow itself.

The effect of a newly born star cluster inside a giant molecular cloud (GMC) is to produce a hot bubble and a thin, dense shell of interstellar gas and dust swept up by the H II expansion, strong stellar winds, and repeated supernova explosions. Lying at the inner side of the shell is the photodissociation region (PDR), the origin of much of the far-infrared/sub-millimeter/millimeter (FIR/sub-mm/mm) radiation from the interstellar medium (ISM). We present a model for the expanding shell at different stages of its expansion which predict mm/sub-mm and far-IR emission line intensities from a series of key molecular and atomic constituents in the shell. The kinematic properties of the swept-up shell predicted by our model are in very good agreement with the measurements of the supershell detected in the nearby starburst galaxy M 82. We compare the modeling results with the ratio-ratio plots of the FIR/sub-mm/mm line emission in the central 1.0 kpc region to investigate the mechanism of star forming activity in M 82. Our model has yielded appropriate gas densities, temperatures, and structure scales compared to those measured in M 82, and the total H2 content is compatible with the observations. This implies that the neutral ISM of the central star-forming region is a product of fragments of the evolving shells.

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.

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.

In cold molecular clouds submillimetre emission lines are excited by the ambient radiation field. The pumping is dominated by the cosmic microwave background (CMB). It is usual in molecular line radiative transfer modelling to simply assume that this is the only incident radiation field. In this paper, a molecular line transport code and a dust radiative transfer code are used to explore the effects of the inclusion of a full interstellar radiation field (ISRF) on a simple test molecular cloud. It is found that in many galactic situations, the shape and strength of the line profiles that result are robust to variations in the ISRF and thus that in most cases, it is safe to adopt the CMB radiation field for the molecular line transport calculations. However, we show that in two examples, the inclusion of a plausible radiation field can have a significant effect on the the line profiles. Firstly, in the vicinity of an embedded massive star, there will be an enhanced far infared component to the radiation field. Secondly, for molecular clouds at large redshift, the CMB temperature increases and this of course also alters the radiation field. In both of these cases, the line profiles are weakened significantly compared to a cloud exposed to a standard radiation field. Therefore this effect should be accounted for when investigating prestellar cores in massive star forming regions and when searching for molecular clouds at high redshift.

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.

Comparison is made between a number of independent computer programs for radiative transfer in molecular rotational lines. The test models are spherically symmetric circumstellar envelopes with a given density and temperature profile. The first two test models have a simple power law density distribution, constant temperature and a fictive 2-level molecule, while the other two test models consist of an inside-out collapsing envelope observed in rotational transitions of HCO+. For the 2-level molecule test problems all codes agree well to within 0.2%, comparable to the accuracy of the individual codes, for low optical depth and up to 2% for high optical depths (tau=4800). The problem of the collapsing cloud in HCO+ has a larger spread in results, ranging up to 12% for the J=4 population. The spread is largest at the radius where the transition from collisional to radiative excitation occurs. The resulting line profiles for the HCO+ J=4-3 transition agree to within 10%, i.e., within the calibration accuracy of most current telescopes. The comparison project and the results described in this paper provide a benchmark for future code development, and give an indication of the typical accuracy of present day calculations of molecular line transfer.

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

We use ground-based imaging polarimetry to detect and image the dusty circumstellar envelopes of a sample of proto-planetary nebulae (PPNe) at near-infrared wavelengths. We detect extended (up to 9 arcsec diameter) circumstellar envelopes around 15 out of 16 sources with a range of morphologies including bipolars and shells. The distribution of scattered light in combination with its polarization (up to 40 per cent) provides unambiguous evidence for axisymmetry in 14 objects showing this to be a common trait of PPNe. We suggest that the range of observed envelope morphologies results from the development of an axisymmetric dust distribution during the superwind phase at the end of the AGB. We identify shells seen in polarized light with scattering from these superwind dust distributions which allows us to provide constraints on the duration of the superwind phase. In one object (IRAS 19475+3119) the circumstellar envelope has a two-armed spiral structure which we suggest results from the interaction of the mass losing star with a binary companion.

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

We report high spatial and spectral resolution measurements of masers towards R Aqr and H1-36, both of which are examples of the sub-class of symbiotic stars that contain a long-period Mira-type variable. Our observations have resulted in the first detection of 321-GHz H2O maser action towards a symbiotic Mira - R Aqr. Comparison with simultaneous 22-GHz H2O maser data suggests the masers do not have the same properties as those in the circumstellar envelopes of field Miras. R Aqr's 22-/321-GHz peak flux density and luminosity ratios are low, as is the line-width ratio. Continuum and spectral-line maps indicate that the 22-GHz maser and free-free emission are aligned. Three mechanisms can reproduce the data with varying degrees of success. All three lead naturally to normal levels of maser emission in SiO and 321-GHz H2O and anomalously weak OH and 22-GHz H2O masers. In the most convincing model, UV radiation and a fast wind from the companion remove the Mira's envelope of dusty, molecular gas, leaving a relatively small cavity of dense, neutral material within a large, ionized nebula. Excitation temperatures suggest that 321-GHz masers are normally excited close to the Mira whilst 22-GHz masers are more remote; in R Aqr, therefore, the 22-GHz masers do not form under optimum conditions. Instead, we see weak and narrow lines that form closer to the Mira, consistent with our high-resolution maps.

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.