We present a detailed analysis of the disk and circumstellar environment of the forming O-type star AFGL 4176 mm1, placing results from the Atacama Large Millimeter/submillimeter Array (ALMA) into context with multiwavelength data. With ALMA, we detect seventeen 1.2 mm continuum sources within 5$''$ (21,000 au) of AFGL 4176 mm1. We find that mm1 has a spectral index of 3.4$\pm$0.2 across the ALMA band, with $>$87% of its 1.2 mm continuum emission from dust. The source mm2, projected 4200 au from mm1, may be a companion or a blueshifted knot in a jet. We also explore the morphological differences between the molecular lines detected with ALMA, finding 203 lines from 25 molecules, which we categorize into several morphological types. Our results show that AFGL 4176 mm1 provides an example of a forming O-star with a large and chemically complex disk, which is mainly traced by nitrogen-bearing molecules. Lines that show strong emission on the blueshifted side of the disk are predominantly oxygen-bearing, which we suggest are tracing a disk accretion shock. The molecules C$^{34}$S, H$_2$CS and CH$_{3}$CCN trace a slow wide-angle wind or dense structures in the outflow cavity walls. With the Australia Telescope Compact Array (ATCA), we detect a compact continuum source ($<$2000 $\times$ 760 au) at 1.2 cm, associated with mm1, of which $>$96% is from ionized gas. The ATCA NH$_3$(1,1) and (2,2) emission traces a large-scale (r$\sim$0.5 pc) rotating toroid with the disk source mm1 in the blueshifted part of this structure offset to the northwest.

We present high-resolution (30 mas or 130 au at 4.2 kpc) Atacama Large Millimeter/submillimeter Array observations at 1.2 mm of the disc around the forming O-type star AFGL 4176 mm1. The disc (AFGL 4176 mm1-main) has a radius of ~1000 au and contains significant structure, most notably a spiral arm on its redshifted side. We fitted the observed spiral with logarithmic and Archimedean spiral models. We find that both models can describe its structure, but the Archimedean spiral with a varying pitch angle fits its morphology marginally better. As well as signatures of rotation across the disc, we observe gas arcs in CH$_3$CN that connect to other millimetre continuum sources in the field, supporting the picture of interactions within a small cluster around AFGL 4176 mm1-main. Using local thermodynamic equilibrium modelling of the CH$_3$CN K-ladder, we determine the temperature and velocity field across the disc, and thus produce a map of the Toomre stability parameter. Our results indicate that the outer disc is gravitationally unstable and has already fragmented or is likely to fragment in the future, possibly producing further companions. These observations provide evidence that disc fragmentation is one possible pathway towards explaining the high fraction of multiple systems around high-mass stars.

Very Long Baseline Interferometry (VLBI) offers unrivalled resolution in studies of celestial radio sources. The subjects of interest of the IAU Symposium No. 356, the Active Galactic Nuclei (AGN) of all types, constitute the major observing sample of modern VLBI networks. At present, the largest in the world in terms of the number of telescopes and geographical coverage is the European VLBI Network (EVN), which operates under the open sky policy via peer-reviewed observing proposals. Recent EVN observations cover a broad range of science themes from high-sensitivity monitoring of structural changes in inner AGN areas to observations of tidal eruptions in AGN cores and investigation of redshift-dependent properties of parsec-scale radio structures of AGN. All the topics above should be considered as potentially rewarding scientific activities of the prospective African VLBI Network (AVN), a natural scientific ally of EVN. This contribution briefly describes the status and near-term strategy for the AVN development as a southern extension of the EVN-AVN alliance and as an eventual bridge to the Square Kilometre Array (SKA) with its mid-frequency core in South Africa.

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

It is important to determine if massive stars form via disc accretion, like their low-mass counterparts. Theory and observation indicate that protostellar jets are a natural consequence of accretion discs and are likely to be crucial for removing angular momentum during the collapse. However, massive protostars are typically rarer, more distant and more dust enshrouded, making observational studies of their jets more challenging. A fundamental question is whether the degree of ionisation in jets is similar across the mass spectrum. Here we determine an ionisation fraction of $\sim5-12\%$ in the jet from the massive protostar G35.20-0.74N, based on spatially coincident infrared and radio emission. This is similar to the values found in jets from lower-mass young stars, implying a unified mechanism of shock ionisation applies in jets across most of the protostellar mass spectrum, up to at least $\sim10$ solar masses.

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

Massive young stellar objects (MYSOs) have recently been shown to drive jets whose particles can interact with either the magnetic fields of the jet or ambient medium to emit non-thermal radiation. We report a search for non-thermal radio emission from a sample of 15 MYSOs to establish the prevalence of the emission in the objects. We used their spectra across the L-, C- and Q-bands along with spectral index maps to characterise their emission. We find that about 50% of the sources show evidence for non-thermal emission with 40% showing clear non-thermal lobes, especially sources of higher bolometric luminosity. The common or IRAS names of the sources that manifest non-thermal lobes are; V645Cyg, IRAS 22134+5834, NGC 7538 IRS 9, IRAS 23262+640, AFGL 402d and AFGL 490. All the central cores of the sources are thermal with corresponding mass-loss rates that lie in the range 3X10^{-7} to 7X10^{-6} solar masses per year. Given the presence of non-thermal lobes in some of the sources and the evidence of non-thermal emission from some spectral index maps, it seems that magnetic fields play a significant role in the jets of massive protostars. Also noted is that some of the sources show evidence of binarity and variability.

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

We present sensitive ($\sigma<10{\rm \mu Jy \, beam^{-1}}$), radio continuum observations using the Australian telescope compact array (ATCA) at frequencies of $6$ and $9$ GHz towards 4 MYSOs. From a previous, less sensitive work, these objects are known to harbour ionized jets associated with radio lobes, which result from shock processes. In comparison with that work, further emission components are detected towards each MYSO. These include extended, direct, thermal emission from the ionized jet's stream, new radio lobes indicative of shocks close ($<10^5\, {\rm au}$) to the MYSO, 3 radio Herbig-Haro objects separated by up to $3.8 {\rm pc}$ from the jet's launching site and an IR-dark source coincident with CH$_3$OH maser emission. No significant, integrated flux variability is detected towards any jets or shocked lobes, and only one proper motion is observed ($1806\pm596\,{\rm km \, s^{-1}}$ parallel to the jet axis of G310.1420+00.7583A). Evidence for precession is detected in all four MYSOs with precession-periods and angles ranging from $66-15480 \, {\rm yr}$ and $6-36 ^\circ$ respectively. Should precession be the result of the influence from a binary companion, we infer orbital radii of $30-1800\, {\rm au} $.

Radio continuum surveys of the Galactic plane can find and characterize HII regions, supernova remnants (SNRs), planetary nebulae (PNe), and extragalactic sources. A number of surveys at high angular resolution (<25") at different wavelengths exist to study the interstellar medium (ISM), but no comparable high-resolution and high-sensitivity survey exists at long radio wavelengths around 21cm. We observed a large fraction of the Galactic plane in the first quadrant of the Milky Way (l=14.0-67.4deg and |b| < 1.25deg) with the Karl G. Jansky Very Large Array (VLA) in the C-configuration covering six continuum spectral windows. These data provide a detailed view on the compact as well as extended radio emission of our Galaxy and thousands of extragalactic background sources. We used the BLOBCAT software and extracted 10916 sources. After removing spurious source detections caused by the sidelobes of the synthesised beam, we classified 10387 sources as reliable detections. We smoothed the images to a common resolution of 25" and extracted the peak flux density of each source in each spectral window (SPW) to determine the spectral indices $\alpha$ (assuming $I(\nu)\propto\nu^\alpha$). By cross-matching with catalogs of HII regions, SNRs, PNe, and pulsars, we found radio counterparts for 840 HII regions, 52 SNRs, 164 PNe, and 38 pulsars. We found 79 continuum sources that are associated with X-ray sources. We identified 699 ultra-steep spectral sources ($\alpha < -1.3$) that could be high-redshift galaxies. Around 9000 of the sources we extracted are not classified specifically, but based on their spatial and spectral distribution, a large fraction of them is likely to be extragalactic background sources. More than 7750 sources do not have counterparts in the SIMBAD database, and more than 3760 sources do not have counterparts in the NED database.

We present a new radio-selected sample of PNe from the CORNISH survey. This is a radio continuum survey of the inner Galactic plane covering Galactic longitude, $10^\circ <l< 65^\circ$ and latitude, $|b| < 1^\circ$ with a resolution of 1.5$^{"}$ and sensitivity better than 0.4 mJy/beam. The radio regime, being unbiased by dust extinction, allows for a more complete sample selection, especially towards the Galactic mid-plane. Visual inspection of the CORNISH data, in combination with data from multi-wavelength surveys of the Galactic plane, allowed the CORNISH team to identify 169 candidate PNe. Here, we explore the use of multi-wavelength diagnostic plots and analysis to verify and classify the candidate PNe. We present the multi-wavelength properties of this new PNe sample. We find 90 new PNe, of which 12 are newly discovered and 78 are newly classified as PN. A further 47 previously suspected PNe are confirmed as such from the analysis presented here and 24 known PNe are detected. Eight sources are classified as possible PNe or other source types. Our sample includes a young sub-sample, with physical diameters < 0.12 pc, brightness temperatures (> 1000 K) and located closer than 7 kpc. Within this sample is a water-maser PN with a spectral index of $-0.55\pm 0.08$, which indicates non-thermal radio emission. Such a radio-selected sample, unaffected by extinction, will be particularly useful to compare with population synthesis models and should contribute to the understanding of the formation and evolution of PNe.

What role does Galactic structure play in star formation? We have used the Herschel Hi-GAL compact-clump catalogue to examine trends in evolutionary stage over large spatial scales in the inner Galaxy. We examine the relationship between the fraction of clumps with embedded star formation (the star-forming fraction, or SFF) and other measures of star-formation activity. Based on a positive correlation between SFF and evolutionary indicators such as the luminosity-to-mass ratio, we assert that the SFF principally traces the average evolutionary state of a sample and must depend on the local fraction of rapidly-evolving, high-mass young stellar objects. The spiral-arm tangent point longitudes show small excesses in the SFF, though these can be accounted for by a small number of the most massive clusters, just 7.6% of the total number of clumps in the catalogue. This suggests that while the arms tend to be home to the Galaxy's massive clusters, the remaining 92.4% of Hi-GAL clumps in our catalogue do not show an enhancement of star formation within arms. Globally, the SFF is highest at the Galactic midplane and inner longitudes. We find no significant trend in evolutionary stage as a function of position across spiral arms at the tangent-point longitudes. This indicates that the angular offset observed between gas and stars, if coordinated by a density wave, is not evident at the clump phase; alternatively, the onset of star formation is not triggered by the spiral density wave.

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation.

We present a composite model and radiative transfer simulations of the massive star forming core W33A MM1. The model was tailored to reproduce the complex features observed with ALMA at $\approx 0.2$ arcsec resolution in CH$_3$CN and dust emission. The MM1 core is fragmented into six compact sources coexisting within $\sim 1000$ au. In our models, three of these compact sources are better represented as disc-envelope systems around a central (proto)star, two as envelopes with a central object, and one as a pure envelope. The model of the most prominent object (Main) contains the most massive (proto)star ($M_\star\approx7~M_\odot$) and disc+envelope ($M_\mathrm{gas}\approx0.4~M_\odot$), and is the most luminous ($L_\mathrm{Main} \sim 10^4~L_\odot$). The model discs are small (a few hundred au) for all sources. The composite model shows that the elongated spiral-like feature converging to the MM1 core can be convincingly interpreted as a filamentary accretion flow that feeds the rising stellar system. The kinematics of this filament is reproduced by a parabolic trajectory with focus at the center of mass of the region. Radial collapse and fragmentation within this filament, as well as smaller filamentary flows between pairs of sources are proposed to exist. Our modelling supports an interpretation where what was once considered as a single massive star with a $\sim 10^3$ au disc and envelope, is instead a forming stellar association which appears to be virialized and to form several low-mass stars per high-mass object.

A catalogue of 239 ultra-compact HII regions (UCHIIs) found in the CORNISH survey at 5 GHz and 1.5" resolution in the region $10^{\circ} < l < 65^{\circ}, ~|b| < 1^{\circ}$ is presented. This is the largest complete and well-selected sample of UCHIIs to date and provides the opportunity to explore the global and individual properties of this key state in massive star formation at multiple wavelengths. The nature of the candidates was validated, based on observational properties and calculated spectral indices, and the analysis is presented in this work. The physical sizes, luminosities and other physical properties were computed by utilising literature distances or calculating the distances whenever a value was not available. The near- and mid-infrared extended source fluxes were measured and the extinctions towards the UCHIIs were computed. The new results were combined with available data at longer wavelengths and the spectral energy distributions (SEDs) were reconstructed for 177 UCHIIs. The bolometric luminosities obtained from SED fitting are presented. By comparing the radio flux densities to previous observational epochs, we find about 5% of the sources appear to be time variable. This first high-resolution area survey of the Galactic plane shows that the total number of UCHIIs in the Galaxy is ~ 750 - a factor of 3-4 fewer than found in previous large area radio surveys. It will form the basis for future tests of models of massive star formation.

We present medium-resolution (R~7000) near-infrared echelle spectroscopic data for 36 MYSOs drawn from the Red MSX Source (RMS) survey. This is the largest sample observed at this resolution at these wavelengths of MYSOs to date. The spectra are characterized mostly by emission from hydrogen recombination lines and accretion diagnostic lines. One MYSO shows photospheric HI absorption, a comparison with spectral standards indicates that the star is an A type star with a low surface gravity, implying that the MYSOs are probably swollen, as also suggested by evolutionary calculations. An investigation of the Br$\gamma$ line profiles shows that most are in pure emission, while 13$\pm$5 % display P Cygni profiles, indicative of outflow, while less than 8$\pm$4 % have inverse P Cygni profiles, indicative of infall. These values are comparable with investigations into the optically bright Herbig Be stars, but not with those of Herbig Ae and T Tauri stars, consistent with the notion that the more massive stars undergo accretion in a different fashion than lower mass objects which are undergoing magnetospheric accretion. Accretion luminosities and rates as derived from the Br$\gamma$ line luminosities agree with results for lower mass sources, providing tentative evidence for massive star formation theories based on scaling of low-mass scenarios. We present Br$\gamma$/Br12 line profile ratios exploiting the fact that optical depth effects can be traced as a function of Doppler shift across the lines. These show that the winds of MYSOs in this sample are nearly equally split between constant, accelerating, and decelerating velocity structures. There are no trends between the types of features we see and bolometric luminosities or near-infrared colours.

Hi-GAL is a large-scale survey of the Galactic plane, performed with Herschel in five infrared continuum bands between 70 and 500 $\mu$m. We present a band-merged catalogue of spatially matched sources and their properties derived from fits to the spectral energy distributions (SEDs) and heliocentric distances, based on the photometric catalogs presented in Molinari et al. (2016a), covering the portion of Galactic plane $-71.0^{\circ}< \ell < 67.0^{\circ}$. The band-merged catalogue contains 100922 sources with a regular SED, 24584 of which show a 70 $\mu$m counterpart and are thus considered proto-stellar, while the remainder are considered starless. Thanks to this huge number of sources, we are able to carry out a preliminary analysis of early stages of star formation, identifying the conditions that characterise different evolutionary phases on a statistically significant basis. We calculate surface densities to investigate the gravitational stability of clumps and their potential to form massive stars. We also explore evolutionary status metrics such as the dust temperature, luminosity and bolometric temperature, finding that these are higher in proto-stellar sources compared to pre-stellar ones. The surface density of sources follows an increasing trend as they evolve from pre-stellar to proto-stellar, but then it is found to decrease again in the majority of the most evolved clumps. Finally, we study the physical parameters of sources with respect to Galactic longitude and the association with spiral arms, finding only minor or no differences between the average evolutionary status of sources in the fourth and first Galactic quadrants, or between "on-arm" and "inter-arm" positions.

We present a comprehensive statistical analysis of star-forming objects located in the vicinities of 1 360 bubble structures throughout the Galactic Plane and their local environments. The compilation of ~70 000 star-forming sources, found in the proximity of the ionized (Hii) regions and detected in both Hi-GAL and GLIMPSE surveys, provided a broad overview of the different evolutionary stages of star-formation in bubbles, from prestellar objects to more evolved young stellar objects (YSOs). Surface density maps of star-forming objects clearly reveal an evolutionary trend where more evolved star-forming objects are found spatially located near the center, while younger star-forming objects are found at the edge of the bubbles. We derived dynamic ages for a subsample of 182 Hii regions for which kinematic distances and radio continuum flux measurements were available. We detect ~80% more star-forming sources per unit area in the direction of bubbles than in the surrounding fields. We estimate ~10% clump formation efficiency (CFE) of Hi-GAL clumps in bubbles, twice the CFE in fields not affected by feedback. We find higher CFE of protostellar clumps in younger bubbles, whose density of the bubble shells is higher. We argue that the formation rate from prestellar to protostellar phase is probably higher during the early stages of the bubble expansion. Evaluation of the fragmentation time inside the shell of bubbles advocates the preexistence of clumps in the medium before the bubble, as supported by numerical simulations. Approximately 23% of the Hi-GAL clumps are found located in the direction of a bubble, with 15% for prestellar clumps and 41% for protostellar clumps. We argue that the high fraction of protostellar clumps may be due to the acceleration of the star-formation process cause by the feedback of the (Hii) bubbles.

The H2O Southern Galactic Plane Survey (HOPS) has mapped 100 square degrees of the Galactic plane for water masers and thermal molecular line emission using the 22-m Mopra telescope. We describe the automated spectral-line fitting pipelines used to determine the properties of emission detected in HOPS datacubes, and use these to derive the physical and kinematic properties of gas in the survey. A combination of the angular resolution, sensitivity, velocity resolution and high critical density of lines targeted make the HOPS data cubes ideally suited to finding precursor clouds to the most massive and dense stellar clusters in the Galaxy. We compile a list of the most massive HOPS ammonia regions and investigate whether any may be young massive cluster progenitor gas clouds. HOPS is also ideally suited to trace the flows of dense gas in the Galactic Centre. We find the kinematic structure of gas within the inner 500pc of the Galaxy is consistent with recent predictions for the dynamical evolution of gas flows in the centre of the Milky Way. We confirm a recent finding that the dense gas in the inner 100pc has an oscillatory kinematic structure with characteristic length scale of ~20pc, and also identify similar oscillatory kinematic structure in the gas at radii larger than 100pc. Finally, we make all of the above fits and the remaining HOPS data cubes across the 100 square degrees of the survey available to the community.

We present the first data release of the James Clerk Maxwell Telescope (JCMT) Plane Survey (JPS), the JPS Public Release 1 (JPSPR1). JPS is an 850-um continuum survey of six fields in the northern inner Galactic Plane in a longitude range of l=7-63, made with the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2). This first data release consists of emission maps of the six JPS regions with an average pixel-to-pixel noise of 7.19 mJy beam^-1, when smoothed over the beam, and a compact-source catalogue containing 7,813 sources. The 95 per cent completeness limits of the catalogue are estimated at 0.04 Jy beam^-1 and 0.3 Jy for the peak and integrated flux densities, respectively. The emission contained in the compact-source catalogue is 42 +- 5 per cent of the total and, apart from the large-scale (greater than 8') emission, there is excellent correspondence with features in the 500-um Herschel maps. We find that, with two-dimensional matching, 98 +- 2 per cent of sources within the fields centred at l=20, 30, 40 and 50 are associated with molecular clouds, with 91 +- 3 per cent of the l=30 and 40 sources associated with dense molecular clumps. Matching the JPS catalogue to Herschel 70-um sources, we find that 38 +- 1 per cent of sources show evidence of ongoing star formation. The images and catalogue will be a valuable resource for studies of star formation in the Galaxy and the role of environment and spiral arms in the star formation process.

We have modelled the evolution of cometary HII regions produced by zero-age main-sequence stars of O and B spectral types, which are driving strong winds and are born off-centre from spherically symmetric cores with power-law ($\alpha = 2$) density slopes. A model parameter grid was produced that spans stellar mass, age and core density. Exploring this parameter space we investigated limb-brightening, a feature commonly seen in cometary HII regions. We found that stars with mass $M_\star \geq 12\, \mathrm{M}_\odot$ produce this feature. Our models have a cavity bounded by a contact discontinuity separating hot shocked wind and ionised ambient gas that is similar in size to the surrounding HII region. Due to early pressure confinement we did not see shocks outside of the contact discontinuity for stars with $M_\star \leq 40\, \mathrm{M}_\odot$, but the cavities were found to continue to grow. The cavity size in each model plateaus as the HII region stagnates. The spectral energy distributions of our models are similar to those from identical stars evolving in uniform density fields. The turn-over frequency is slightly lower in our power-law models due to a higher proportion of low density gas covered by the HII regions.

We targeted the massive star forming region W33A using the Atacama Large Sub/Millimeter Array (ALMA) in band 6 (230 GHz) and 7 (345 GHz) to search for a sub-1000au disc around the central O-type massive young stellar object (MYSO) W33A MM1-Main. Our data achieves a resolution of ~0.2" (~500au) and resolves the central core, MM1, into multiple components and reveals complex and filamentary structures. There is strong molecular line emission covering the entire MM1 region. The kinematic signatures are inconsistent with only Keplerian rotation although we propose that the shift in the emission line centroids within ~1000au of MM1-Main could hint at an underlying compact disc with Keplerian rotation. We cannot however rule out the possibility of an unresolved binary or multiple system. A putative smaller disc could be fed by the large scale spiral `feeding filament' we detect in both gas and dust emission. We also discuss the nature of the now-resolved continuum sources.

Until recently, there have been few studies of the protostellar evolution of intermediate-mass (IM) stars, which may bridge the low- and high-mass regimes. This paper aims to investigate whether the properties of an IM protostar within the IRAS 05373+2349 embedded cluster are similar to that of low- and/or high-mass protostars. We carried out Very Large Array as well as Combined Array for Research in Millimeter Astronomy continuum and 12CO(J=1-0) observations, which uncover seven radio continuum sources (VLA 1-7). The spectral index of VLA 2, associated with the IM protostar is consistent with an ionised stellar wind or jet. The source VLA 3 is coincident with previously observed H2 emission line objects aligned in the north-south direction (P.A. -20 to -12 deg), which may be either an ionised jet emanating from VLA 2 or (shock-)ionised cavity walls in the large-scale outflow from VLA 2. The position angle between VLA 2 and 3 is slightly misaligned with the large-scale outflow we map at ~5-arcsec resolution in 12CO (P.A. ~30 deg), which in the case of a jet suggests precession. The emission from the mm core associated with VLA 2 is also detected; we estimate its mass to be 12-23 Msun, depending on the contribution from ionised gas. Furthermore, the large-scale outflow has properties intermediate between outflows from low- and high-mass young stars. Therefore, we conclude that the IM protostar within IRAS 05373+2349 is phenomenologically as well as quantitatively intermediate between the low- and high-mass domains.

We present large-scale trends in the distribution of star-forming objects revealed by the Hi-GAL survey. As a simple metric probing the prevalence of star formation in Hi-GAL sources, we define the fraction of the total number of Hi-GAL sources with a 70-micron counterpart as the "star-forming fraction" or SFF. The mean SFF in the inner galactic disc (3.1 kpc < R_GC < 8.6 kpc) is 25%. Despite an apparent pile-up of source numbers at radii associated with spiral arms, the SFF shows no significant deviations at these radii, indicating that the arms do not affect the star-forming productivity of dense clumps either via physical triggering processes or through the statistical effects of larger source samples associated with the arms. Within this range of Galactocentric radii, we find that the SFF declines with R_GC at a rate of -0.026 +/- 0.002 per kiloparsec, despite the dense gas mass fraction having been observed to be constant in the inner Galaxy. This suggests that the SFF may be weakly dependent on one or more large-scale physical properties of the Galaxy, such as metallicity, radiation field, pressure or shear, such that the dense sub-structures of molecular clouds acquire some internal properties inherited from their environment.

Radio continuum observations using the Australia telescope compact array at 5.5, 9.0, 17.0 and 22.8 GHz have detected free-free emission associated with 45 of 49 massive young stellar objects and HII regions. Of these, 26 sources are classified as ionized jets (12 of which are candidates), 2 as ambiguous jets or disc winds, 1 as a disc-wind, 14 as HII regions and 2 were unable to be categorised. Classification as ionized jets is based upon morphology, radio flux and spectral index, in conjunction with previous observational results at other wavelengths. Radio-luminosity and momentum are found to scale with bolometric luminosity in the same way as low-mass jets, indicating a common mechanism for jet production across all masses. In 13 of the jets, we see associated non-thermal/optically-thin lobes resulting from shocks either internal to the jet and/or at working surfaces. Ten jets display non-thermal (synchrotron emission) spectra in their lobes, with an average spectral index of -0.55 consistent with Fermi acceleration in shocks. This shows that magnetic fields are present, in agreement with models of jet formation incorporating magnetic fields. Since the production of collimated radio jets is associated with accretion processes, the results presented in this paper support the picture of disc-mediated accretion for the formation of massive stars with an upper-limit on the jet phase lasting approximately $6.5 \times 10^4 yr$. Typical mass loss rates in the jet are found to be $1.4 \times 10^{-5} M_\odot yr^{-1}$ with associated momentum rates of the order $(1-2) \times 10^{-2} M_\odot km s^{-1} yr^{-1}$.

(Abridged) We present the first public release of high-quality data products (DR1) from Hi-GAL, the {\em Herschel} infrared Galactic Plane Survey. Hi-GAL is the keystone of a suite of continuum Galactic Plane surveys from the near-IR to the radio, and covers five wavebands at 70, 160, 250, 350 and 500 micron, encompassing the peak of the spectral energy distribution of cold dust for 8 < T < 50K. This first Hi-GAL data release covers the inner Milky Way in the longitude range 68{\deg} > l > -70{\deg} in a |b|<1{\deg} latitude strip. Photometric maps have been produced with the ROMAGAL pipeline, that optimally capitalizes on the excellent sensitivity and stability of the bolometer arrays of the {\em Herschel} PACS and SPIRE photometric cameras, to deliver images of exquisite quality and dynamical range, absolutely calibrated with {\em Planck} and {\em IRAS}, and recovering extended emission at all wavelengths and all spatial scales. The compact source catalogues have been generated with the CuTEx algorithm, specifically developed to optimize source detection and extraction in the extreme conditions of intense and spatially varying background that are found in the Galactic Plane in the thermal infrared. Hi-GAL DR1 images will be accessible via a dedicated web-based image cutout service. The DR1 Compact Source Catalogues are delivered as single-band photometric lists containing, in addition to source position, peak and integrated flux and source sizes, a variety of parameters useful to assess the quality and reliability of the extracted sources, caveats and hints to help this assessment are provided. Flux completeness limits in all bands are determined from extensive synthetic source experiments and depend on the specific line of sight along the Galactic Plane. Hi-GAL DR1 catalogues contain 123210, 308509, 280685, 160972 and 85460 compact sources in the five bands, respectively.

We present Atacama Large Millimeter/submillimeter Array (ALMA) line and continuum observations at 1.2mm with ~0.3" resolution that uncover a Keplerian-like disk around the forming O-type star AFGL 4176. The continuum emission from the disk at 1.21 mm (source mm1) has a deconvolved size of 870+/-110 AU x 330+/-300 AU and arises from a structure ~8 M_sun in mass, calculated assuming a dust temperature of 190 K. The first-moment maps, pixel-to-pixel line modeling, assuming local thermodynamic equilibrium (LTE), and position-velocity diagrams of the CH3CN J=13-12 K-line emission all show a velocity gradient along the major axis of the source, coupled with an increase in velocity at small radii, consistent with Keplerian-like rotation. The LTE line modeling shows that where CH3CN J=13-12 is excited, the temperatures in the disk range from ~70 to at least 300 K and that the H2 column density peaks at 2.8x10^24 cm^-2. In addition, we present Atacama Pathfinder Experiment (APEX) 12CO observations which show a large-scale outflow from AFGL 4176 perpendicular to the major axis of mm1, supporting the disk interpretation. Finally, we present a radiative transfer model of a Keplerian disk surrounding an O7 star, with a disk mass and radius of 12 M_sun and 2000 AU, that reproduces the line and continuum data, further supporting our conclusion that our observations have uncovered a Keplerian disk around an O-type star.

We have observed 99 mid-infrared-bright, massive young stellar objects and compact HII regions drawn from the Red MSX source (RMS) survey in the J=3$-$2 transition of $^{12}$CO and $^{13}$CO, using the James Clerk Maxwell Telescope. 89 targets are within 6 kpc of the Sun, covering a representative range of luminosities and core masses. These constitute a relatively unbiased sample of bipolar molecular outflows associated with massive star formation. Of these, 59, 17 and 13 sources (66, 19 and 15 percent) are found to have outflows, show some evidence of outflow, and have no evidence of outflow, respectively. The time-dependent parameters of the high-velocity molecular flows are calculated using a spatially variable dynamic timescale. The canonical correlations between the outflow parameters and source luminosity are recovered and shown to scale with those of low-mass sources. For coeval star formation we find the scaling is consistent with all the protostars in an embedded cluster providing the outflow force, with massive stars up to $\sim$30 M$_{\odot}$ generating outflows. Taken at face value, the results support the model of a scaled-up version of the accretion-related outflow-generation mechanism associated with discs and jets in low-mass objects with time-averaged accretion rates of $\sim$10$^{-3}$ M$_{\odot}$ yr$^{-1}$ onto the cores. However, we also suggest an alternative model, in which the molecular outflow dynamics are dominated by the entrained mass and are unrelated to the details of the acceleration mechanism. We find no evidence that outflows contribute significantly to the turbulent kinetic energy of the surrounding dense cores.

We present early results from the JCMT Plane Survey (JPS), which has surveyed the northern inner Galactic plane between longitudes l=7 and l=63 degrees in the 850-{\mu}m continuum with SCUBA-2, as part of the James Clerk Maxwell Telescope Legacy Survey programme. Data from the l=30 degree survey region, which contains the massive star-forming regions W43 and G29.96, are analysed after approximately 40% of the observations had been completed. The pixel-to-pixel noise is found to be 19 mJy/beam, after a smooth over the beam area, and the projected equivalent noise levels in the final survey are expected to be around 10 mJy/beam. An initial extraction of compact sources was performed using the FellWalker method resulting in the detection of 1029 sources above a 5-{\sigma} surface-brightness threshold. The completeness limits in these data are estimated to be around 0.2 Jy/beam (peak flux density) and 0.8 Jy (integrated flux density) and are therefore probably already dominated by source confusion in this relatively crowded section of the survey. The flux densities of extracted compact sources are consistent with those of matching detections in the shallower ATLASGAL survey. We analyse the virial and evolutionary state of the detected clumps in the W43 star-forming complex and find that they appear younger than the Galactic-plane average.

We used Herschel Hi-GAL survey data to determine whether massive young stellar objects (MYSOs) are resolved at 70$\mu$m and to study their envelope density distribution. Our analysis of three relatively isolated sources in the l=30{\deg} and l=59{\deg} Galactic fields show that the objects are partially resolved at 70$\mu$m. The Herschel Hi-GAL survey data have a high scan velocity which makes unresolved and partially resolved sources appear elongated in the 70$\mu$m images. We analysed the two scan directions separately and examine the intensity profile perpendicular to the scan direction. Spherically symmetric radiative transfer models with a power law density distribution were used to study the circumstellar matter distribution. Single dish sub-mm data were also included to study how different spatial information affects the fitted density distribution. The density distribution which best fits both the 70$\mu$m intensity profile and SED has an average index of ~0.5. This index is shallower than expected and is probably due to the dust emission from bipolar outflow cavity walls not accounted for in the spherical models. We conclude that 2D axisymmetric models and Herschel images at low scan speeds are needed to better constrain the matter distribution around MYSOs.

Very Long Baseline Interferometry (VLBI) at radio wavelengths can provide astrometry accurate to 10 micro-arcseconds or better (i.e. better than the target GAIA accuracy) without being limited by dust obscuration. This means that unlike GAIA, VLBI can be applied to star-forming regions independently of their internal and line-of-sight extinction. Low-mass young stellar objects (particularly T Tauri stars) are often non-thermal compact radio emitters, ideal for astrometric VLBI radio continuum experiments. Existing observations for nearby regions (e.g. Taurus, Ophiuchus, or Orion) demonstrate that VLBI astrometry of such active T Tauri stars enables the reconstruction of both the regions' 3D structure (through parallax measurements) and their internal kinematics (through proper motions, combined with radial velocities). The extraordinary sensitivity of the SKA telescope will enable similar "tomographic mappings" to be extended to regions located several kpc from Earth, in particular to nearby spiral arm segments. This will have important implications for Galactic science, galactic dynamics and spiral structure theories.

The SKA will be a state of the art radiotelescope optimized for both large area surveys as well as for deep pointed observations. In this paper we analyze the impact that the SKA will have on Galactic studies, starting from the immense legacy value of the all-sky survey proposed by the continuum SWG but also presenting some areas of Galactic Science that particularly benefit from SKA observations both surveys and pointed. The planned all-sky survey will be characterized by unique spatial resolution, sensitivity and survey speed, providing us with a wide-field atlas of the Galactic continuum emission. Synergies with existing, current and planned radio Galactic Plane surveys will be discussed. SKA will give the opportunity to create a sensitive catalog of discrete Galactic radio sources, most of them representing the interaction of stars at various stages of their evolution with the environment: complete census of all stage of HII regions evolution; complete census of late stages of stellar evolution such as PNe and SNRs; detection of stellar winds, thermal jets, Symbiotic systems, Chemically Peculiar and dMe stars, active binary systems in both flaring and quiescent states. Coherent emission events like Cyclotron Maser in the magnetospheres of different classes of stars can be detected. Pointed, deep observations will allow new insights into the physics of the coronae and plasma processes in active stellar systems and single stars, enabling the detection of flaring activity in larger stellar population for a better comprehension of the mechanism of energy release in the atmospheres of stars with different masses and age.

The vast collecting area of the Square Kilometre Array (SKA), harnessed by sensitive receivers, flexible digital electronics and increased computational capacity, could permit the most sensitive and exhaustive search for technologically-produced radio emission from advanced extraterrestrial intelligence (SETI) ever performed. For example, SKA1-MID will be capable of detecting a source roughly analogous to terrestrial high-power radars (e.g. air route surveillance or ballistic missile warning radars, EIRP (EIRP = equivalent isotropic radiated power, ~10^17 erg sec^-1) at 10 pc in less than 15 minutes, and with a modest four beam SETI observing system could, in one minute, search every star in the primary beam out to ~100 pc for radio emission comparable to that emitted by the Arecibo Planetary Radar (EIRP ~2 x 10^20 erg sec^-1). The flexibility of the signal detection systems used for SETI searches with the SKA will allow new algorithms to be employed that will provide sensitivity to a much wider variety of signal types than previously searched for. Here we discuss the astrobiological and astrophysical motivations for radio SETI and describe how the technical capabilities of the SKA will explore the radio SETI parameter space. We detail several conceivable SETI experimental programs on all components of SKA1, including commensal, primary-user, targeted and survey programs and project the enhancements to them possible with SKA2. We also discuss target selection criteria for these programs, and in the case of commensal observing, how the varied use cases of other primary observers can be used to full advantage for SETI.

By matching infrared-selected, massive young stellar objects (MYSOs) and compact HII regions in the RMS survey to massive clumps found in the submillimetre ATLASGAL survey, we have identified ~1000 embedded young massive stars between 280\degr < $\ell$ < 350\degr and 10degr < $\ell$ < 60\degr with |b|<1.5degr. Combined with an existing sample of radio-selected methanol masers and compact HII regions, the result is a catalogue of ~1700 massive stars embedded within ~1300 clumps located across the inner Galaxy, containing three observationally distinct subsamples, methanol-maser, MYSO and HII-region associations, covering the most important tracers of massive star formation, thought to represent key stages of evolution. We find that massive star formation is strongly correlated with the regions of highest column density in spherical, centrally condensed clumps. We find no significant differences between the three samples in clump structure or the relative location of the embedded stars, which suggests that the structure of a clump is set before the onset of star formation, and changes little as the embedded object evolves towards the main sequence. There is a strong linear correlation between clump mass and bolometric luminosity, with the most massive stars forming in the most massive clumps. We find that the MYSO and HII-region subsamples are likely to cover a similar range of evolutionary stages and that the majority are near the end of their main accretion phase. We find few infrared-bright MYSOs associated with the most massive clumps, probably due to very short pre-main sequence lifetimes in the most luminous sources.

The VISTA Variables in the V\'ia L\'actea (VVV) survey is one of six public ESO surveys, and is now in its 4th year of observing. Although far from being complete, the VVV survey has already delivered many results, some directly connected to the intended science goals (detection of variables stars, microlensing events, new star clusters), others concerning more exotic objects, e.g. novae. Now, at the end of the fourth observing period, and comprising roughly 50% of the proposed observations, the actual status of the survey, as well some of the results based on the VVV data, are presented.

Abridged: We have used the well-selected sample of ~1750 embedded, young, massive stars identified by the RMS survey to investigate the Galactic distribution of recent massive star formation. We describe the various methods used to assign distances extracted from the literature, and solve the distance ambiguities towards ~200 sources located within the Solar circle using archival HI data. These distances are used to calculate bolometric luminosities and estimate the survey completeness (~2x10^4 lsun). In total, we calculate the distance and luminosity of ~1650 sources, one third of which are above the survey's completeness threshold. Examination of the sample's longitude, latitude, radial velocities and mid-infrared images has identified ~120 small groups of sources, many of which are associated with well known star formation complexes, such as W43, W49 and W51. We compare the positional distribution of the sample with the expected locations of the spiral arms, assuming a model of the Galaxy consisting of four gaseous arms. The distribution of young massive stars in the Milky Way is spatially correlated with the spiral arms, with strong peaks in the source position and luminosity distributions at the arms' Galactocentric radii. The overall source and luminosity surface densities are both well correlated with the surface density of the molecular gas, which suggests that the massive star formation rate per unit molecular mass is approximately constant across the Galaxy. We estimate the total luminosity of the embedded massive star population to be ~0.76x10^8 lsun, 30% of which is associated with the ten most active star forming complexes. We measure the scale height as a function of Galactocentric distance and find that it increases only modestly from ~20-30 pc between 4 and 8 kpc, but much more rapidly at larger distances.

We present medium resolution ($R\approx5300$) $K$-band integral field spectroscopy of six MYSOs. The targets are selected from the RMS survey, and we used the NIFS on the Gemini North telescope. The data show various spectral line features including Br$\gamma$, CO, H$_2$, and \mbox{He\,{\sc i}}. The Br$\gamma$ line is detected in emission in all objects with $v_\mathrm{FWHM}\sim100$ -- 200 kms$^{-1}$. V645 Cyg shows a high-velocity P-Cygni profile between -800 kms$^{-1}$ and -300 kms$^{-1}$. We performed three-dimensional spectroastrometry to diagnose the circumstellar environment in the vicinity of the central stars using the Br$\gamma$ line. We measured the centroids of the velocity components with sub-mas precision. The centroids allow us to discriminate the blueshifted and redshifted components in a roughly east--west direction in both IRAS 18151--1208 and S106 in Br$\gamma$. This lies almost perpendicular to observed larger scale outflows. We conclude, given the widths of the lines and the orientation of the spectroastrometric signature, that our results trace a disc wind in both IRAS 18151--1208 and S106. The CO $\nu=2-0$ absorption lines at low $J$ transitions are detected in IRAS 18151--1208 and AFGL 2136. We analysed the velocity structure of the neutral gas discs. In IRAS 18151--1208, the absorption centroids of the blueshifted and redshifted components are separated in a direction of north-east to south-west, nearly perpendicular to that of the larger scale $H_2$ jet. The position-velocity relations of these objects can be reproduced with central masses of 30 $M_{\sun}$ for IRAS 18151--1208 and 20 $M_{\sun}$ for AFGL 2136. We also detect CO $\nu=2-0$ bandhead emission in IRAS 18151--1208, S106 and V645 Cyg. The results can be fitted reasonably with a Keplerian rotation model, with masses of 15, 20 and 20 $M_{\sun}$ respectively.

We present the Red MSX Source (RMS) Survey, the largest statistically selected catalog of young massive protostars and HII regions to date. We outline the construction of the catalog using mid and near infrared color selection, as well as the detailed follow up work at other wavelengths, and at higher spatial resolution in the infrared. We show that within the adopted selection bounds we are more than 90% complete for the massive protostellar population, with a positional accuracy of the exciting source of better than 2 arcseconds. We briefly summarize some of the results that can be obtained from studying the properties of the objects in the catalog as a whole, and find evidence that the most massive stars form: (i) preferentially nearer the Galactic centre than the anti-centre; (ii) in the most heavily reddened environments, suggestive of high accretion rates; and (iii) from the most massive cloud cores.

We present a complete sample of molecular clumps containing compact and ultra-compact (UC) HII regions between \ell=10\degr and 60\degr\ and $|b|<1\degr, identified by combining the the ATLASGAL submm and CORNISH radio continuum surveys with visual examination of archival infrared data. Our sample is complete to optically thin, compact and UCHII regions driven by a zero age main sequence star of spectral type B0 or earlier embedded within a 1,000 Msun clump. In total we identify 213 compact and UCHII regions, associated with 170 clumps. Unambiguous kinematic distances are derived for these clumps and used to estimate their masses and physical sizes, as well as the Lyman continuum fluxes and sizes of their embedded HII regions. We find a clear lower envelope for the surface density of molecular clumps hosting massive star formation of 0.05 g cm^{-2}, which is consistent with a similar sample of clumps associated with 6.7 GHz masers. The mass of the most massive embedded stars is closely correlated with the mass of their natal clump. Young B stars appear to be significantly more luminous in the ultraviolet than predicted by current stellar atmosphere models. The properties of clumps associated with compact and UCHII regions are very similar to those associated with 6.7 GHz methanol masers and we speculate that there is little evolution in the structure of the molecular clumps between these two phases. Finally, we identify a significant peak in the surface density of compact and UCHII regions associated with the W49A star-forming complex, noting that this complex is truly one of the most massive and intense regions of star formation in the Galaxy.

Near-infrared H- and K-band spectra are presented for 247 objects, selected from the Red MSX Source (RMS) survey as potential young stellar objects (YSOs). 195 (~80%) of the targets are YSOs, of which 131 are massive YSOs (L_BOL > 5x10^3 L_solar), M > 8M_solar. This is the largest spectroscopic study of massive YSOs to date, providing a valuable resource for the study of massive star formation. In this paper we present our exploratory analysis of the data. The YSOs observed have a wide range of embeddedness (2.7 < A_V < 114), demonstrating that this study covers minimally obscured objects right through to very red, dusty sources. Almost all YSOs show some evidence for emission lines, though there is a wide variety of observed properties. The most commonly detected lines are Brgamma, H_2, fluorescent FeII, CO bandhead, [FeII] and HeI 2-1 2^1S-2^1P, in order of frequency of occurrence. In total, ~40% of the YSOs display either fluorescent FeII 1.6878um or CO bandhead emission (or both), indicative of a circumstellar disc; however, no correlation of the strength of these lines with bolometric luminosity was found. We also find that ~60% of the sources exhibit [FeII] or H_2 emission, indicating the presence of an outflow. Three quarters of all sources have Brgamma in emission. A good correlation with bolometric luminosity was observed for both the Brgamma and H_2 emission line strengths, covering 1 L_solar< L_BOL < 3.5x10^5 L_solar. This suggests that the emission mechanism for these lines is the same for low-, intermediate-, and high-mass YSOs, i.e. high-mass YSOs appear to resemble scaled-up versions of low-mass YSOs.

The CORNISH project is the highest resolution radio continuum survey of the Galactic plane to date. It is the 5 GHz radio continuum part of a series of multi-wavelength surveys that focus on the northern GLIMPSE region (10 deg < l < 65 deg), observed by the Spitzer satellite in the mid-infrared. Observations with the Very Large Array in B and BnA configurations have yielded a 1.5" resolution Stokes I map with a root-mean-squared noise level better than 0.4 mJy/beam. Here we describe the data-processing methods and data characteristics, and present a new, uniform catalogue of compact radio-emission. This includes an implementation of automatic deconvolution that provides much more reliable imaging than standard CLEANing. A rigorous investigation of the noise characteristics and reliability of source detection has been carried out. We show that the survey is optimised to detect emission on size scales up to 14" and for unresolved sources the catalogue is more than 90 percent complete at a flux density of 3.9 mJy. We have detected 3,062 sources above a 7-sigma detection limit and present their ensemble properties. The catalogue is highly reliable away from regions containing poorly-sampled extended emission, which comprise less than two percent of the survey area. Imaging problems have been mitigated by down-weighting the shortest spacings and potential artefacts flagged via a rigorous manual inspection with reference to the Spitzer infrared data. We present images of the most common source types found: regions, planetary nebulae and radio-galaxies. The CORNISH data and catalogue are available online at http://cornish.leeds.ac.uk

Massive stars play an important role in many areas of astrophysics, but numerous details regarding their formation remain unclear. In this paper we present and analyse high resolution (R ~ 30,000) near-infrared 2.3 micron spectra of 20 massive young stellar objects from the RMS database, in the largest such study of CO first overtone bandhead emission to date. We fit the emission under the assumption it originates from a circumstellar disc in Keplerian rotation. We explore three approaches to modelling the physical conditions within the disc - a disc heated mainly via irradiation from the central star, a disc heated mainly via viscosity, and a disc in which the temperature and density are described analytically. We find that the models described by heating mechanisms are inappropriate because they do not provide good fits to the CO emission spectra. We therefore restrict our analysis to the analytic model, and obtain good fits to all objects that possess sufficiently strong CO emission, suggesting circumstellar discs are the source of this emission. On average, the temperature and density structure of the discs correspond to geometrically thin discs, spread across a wide range of inclinations. Essentially all the discs are located within the dust sublimation radius, providing strong evidence that the CO emission originates close to the central protostar, on astronomical unit scales. In addition, we show that the objects in our sample appear no different to the general population of MYSOs in the RMS database, based on their near- and mid-infrared colours. The combination of observations of a large sample of MYSOs with CO bandhead emission and our detailed modelling provide compelling evidence of the presence of small scale gaseous discs around such objects, supporting the scenario in which massive stars form via disc accretion.

We describe the motivation, design and implementation of the CORNISH survey, an arcsecond resolution radio continuum survey of the inner Galactic plane at 5 GHz using the Karl G. Jansky Very Large Array (VLA). It is a blind survey co-ordinated with the northern Spitzer GLIMPSE I region covering 10 deg < l < 65 deg and |b| < 1 deg at similar resolution. We discuss in detail the strategy that we employed to control the shape of the synthesised beam across this survey that covers a wide range of fairly low declinations. Two snapshots separated by 4 hours in hour angle kept the beam elongation to less that 1.5 over 75% of the survey area and less than 2 over 98% of the survey. The prime scientific motivation is to provide an unbiased survey for ultra-compact HII regions to study this key phase in massive star formation. A sensitivity around 2 mJy will allow the automatic distinction between radio loud and quiet mid-IR sources found in the Spitzer surveys. This survey has many legacy applications beyond star formation including evolved stars, active stars and binaries, and extragalactic sources. The CORNISH survey for compact ionized sources complements other Galactic plane surveys that target diffuse and non-thermal sources as well as atomic and molecular phases to build up a complete picture of the ISM in the Galaxy.

The H2O Southern Galactic Plane Survey (HOPS) has mapped a 100 degree strip of the Galactic plane (-70deg > l > 30deg, |b| < 0.5deg) using the 22-m Mopra antenna at 12-mm wavelengths. Observations were conducted in on-the-fly mode using the Mopra spectrometer (MOPS), targeting water masers, thermal molecular emission and radio-recombination lines. Foremost among the thermal lines are the 23 GHz transitions of NH3 J,K = (1,1) and (2,2), which trace the densest parts of molecular clouds (n > 10^4 cm^{-3}). In this paper we present the NH3 (1,1) and (2,2) data, which have a resolution of 2 arcmin and cover a velocity range of +/-200 km/s. The median sensitivity of the NH3 data-cubes is sigma_Tmb = 0.20 +/1 0.06 K. For the (1,1) transition this sensitivity equates to a 3.2 kpc distance limit for detecting a 20 K, 400 Msun cloud at the 5-sigma level. Similar clouds of mass 5,000 Msun would be detected as far as the Galactic centre, while 30,000 Msun clouds would be seen across the Galaxy. We have developed an automatic emission finding procedure based on the ATNF DUCHAMP software and have used it to create a new catalogue of 669 dense molecular clouds. The catalogue is 100 percent complete at the 5-sigma detection limit (Tmb = 1.0 K). A preliminary analysis of the ensemble cloud properties suggest that the near kinematic distances are favoured. The cloud positions are consistent with current models of the Galaxy containing a long bar. Combined with other Galactic plane surveys this new molecular-line dataset constitutes a key tool for examining Galactic structure and evolution. Data-cubes, spectra and catalogues are available to the community via the HOPS website.

A survey of the Milky Way disk and the Magellanic System at the wavelengths of the 21-cm atomic hydrogen (HI) line and three 18-cm lines of the OH molecule will be carried out with the Australian Square Kilometre Array Pathfinder telescope. The survey will study the distribution of HI emission and absorption with unprecedented angular and velocity resolution, as well as molecular line thermal emission, absorption, and maser lines. The area to be covered includes the Galactic plane (|b|< 10deg) at all declinations south of delta = +40deg, spanning longitudes 167deg through 360deg to 79deg at b=0deg, plus the entire area of the Magellanic Stream and Clouds, a total of 13,020 square degrees. The brightness temperature sensitivity will be very good, typically sigma_T ~ 1 K at resolution 30arcsec and 1 km/s. The survey has a wide spectrum of scientific goals, from studies of galaxy evolution to star formation, with particular contributions to understanding stellar wind kinematics, the thermal phases of the interstellar medium, the interaction between gas in the disk and halo, and the dynamical and thermal states of gas at various positions along the Magellanic Stream.

We present a study of the kinematical properties of a small sample of nearby near-infrared bright massive and intermediate mass young stellar objects using emission lines sensitive to discs and winds. We show for the first time that the broad ($\sim500$kms$^{-1}$) symmetric line wings on the HI Brackett series lines are due to Stark broadening or electron scattering, rather than pure Doppler broadening due to high speed motion. The results are consistent with the presence of a very dense circumstellar environment. In addition, many of these lines show evidence for weak line self-absorption, suggestive of a wind or disc-wind origin for that part of the absorbing material. The weakness of the self-absorption suggests a large opening angle for such an outflow. We also study the fluorescent 1.688$\mu$m FeII line, which is sensitive to dense material. We fitted a Keplerian disc model to this line, and find reasonable fits in all bar one case, in agreement with previous finding for classical Be stars that fluorescent iron transitions are reasonable disc tracers. Overall the picture is one in which these stars still have accretion discs, with a very dense inner circumstellar environment which may be tracing either the inner regions of a disc, or of a stellar wind, and in which ionised outflow is also present. The similarity with lower mass stars is striking, suggesting that at least in this mass range they form in a similar fashion.

We test two different methods of using near-infrared extinction to estimate distances to dark clouds in the first quadrant of the Galaxy using large near infrared (2MASS and UKIDSS) surveys. VLBI parallax measurements of masers around massive young stars provide the most direct and bias-free measurement of the distance to these dark clouds. We compare the extinction distance estimates to these maser parallax distances. We also compare these distances to kinematic distances, including recent re-calibrations of the Galactic rotation curve. The extinction distance methods agree with the maser parallax distances (within the errors) between 66% and 100% of the time (depending on method and input survey) and between 85% and 100% of the time outside of the crowded Galactic center. Although the sample size is small, extinction distance methods reproduce maser parallax distances better than kinematic distances; furthermore, extinction distance methods do not suffer from the kinematic distance ambiguity. This validation gives us confidence that these extinction methods may be extended to additional dark clouds where maser parallaxes are not available.

Massive stars form whilst they are still embedded in dense envelopes. As a result, the roles of rotation, mass loss and accretion in massive star formation are not well understood. This study evaluates the source of the Q-band, lambda=19.5 microns, emission of massive young stellar objects (MYSOs). This allows us to determine the relative importance of rotation and outflow activity in shaping the circumstellar environments of MYSOs on 1000 AU scales. We obtained diffraction limited mid-infrared images of a sample of 20 MYSOs using the VLT/VISIR and Subaru/COMICS instruments. For these 8 m class telescopes and the sample selected, the diffraction limit, ~0.6", corresponds to approximately 1000 AU. We compare the images and the spectral energy distributions (SEDs) observed to a 2D, axis-symmetric dust radiative transfer model that reproduces VLTI/MIDI observations of the MYSO W33A. We vary the inclination, mass infall rate, and outflow opening angle to simultaneously recreate the behaviour of the sample of MYSOs in the spatial and spectral domains. The mid-IR emission of 70 percent of the MYSOs is spatially resolved. In the majority of cases, the spatial extent of their emission and their SEDs can be reproduced by the W33A model featuring an in-falling, rotating dusty envelope with outflow cavities. There is independent evidence that most of the sources which are not fit by the model are associated with ultracompact HII regions and are thus more evolved. We find that, in general, the diverse 20 micron morphology of MYSOs can be attributed to warm dust in the walls of outflow cavities seen at different inclinations. This implies that the warm dust in the outflow cavity walls dominates the Q-band emission of MYSOs. In turn, this emphasises that outflows are an ubiquitous feature of massive star formation.

The ESO Public Survey VISTA Variables in the Via Lactea (VVV) started in 2010. VVV targets 562 sq. deg in the Galactic bulge and an adjacent plane region and is expected to run for ~5 years. In this paper we describe the progress of the survey observations in the first observing season, the observing strategy and quality of the data obtained. The observations are carried out on the 4-m VISTA telescope in the ZYJHKs filters. In addition to the multi-band imaging the variability monitoring campaign in the Ks filter has started. Data reduction is carried out using the pipeline at the Cambridge Astronomical Survey Unit. The photometric and astrometric calibration is performed via the numerous 2MASS sources observed in each pointing. The first data release contains the aperture photometry and astrometric catalogues for 348 individual pointings in the ZYJHKs filters taken in the 2010 observing season. The typical image quality is ~0.9-1.0". The stringent photometric and image quality requirements of the survey are satisfied in 100% of the JHKs images in the disk area and 90% of the JHKs images in the bulge area. The completeness in the Z and Y images is 84% in the disk, and 40% in the bulge. The first season catalogues contain 1.28x10^8 stellar sources in the bulge and 1.68x10^8 in the disk area detected in at least one of the photometric bands. The combined, multi-band catalogues contain more than 1.63x10^8 stellar sources. About 10% of these are double detections due to overlapping adjacent pointings. These overlapping multiple detections are used to characterise the quality of the data. The images in the JHKs bands extend typically ~4 mag deeper than 2MASS. The magnitude limit and photometric quality depend strongly on crowding in the inner Galactic regions. The astrometry for Ks=15-18 mag has rms ~35-175 mas.

We present high-resolution HI data obtained using the Australia Telescope Compact Array to resolve the near/far distance ambiguities towards a sample of compact HII regions from the Red MSX Source (RMS) survey. The high resolution data are complemented with lower resolution archival HI data extracted from the Southern and VLA Galactic Plane surveys. We resolve the distance ambiguity for nearly all of the 105 sources where the continuum was strong enough to allow analysis of the HI absorption line structure. This represents another step in the determination of distances to the total RMS sample, which with over 1,000 massive young stellar objects and compact HII regions, is the largest and most complete sample of its kind. The full sample will allow the distribution of massive star formation in the Galaxy to be examined.

A new Monte Carlo method has been developed in order to derive ages of young embedded clusters within massive star forming regions where there is strong differential reddening. After foreground and infrared excess source candidates are removed, each cluster candidate star is individually dereddened. Simulated clusters are constructed using isochrones, an IMF, realistic photometric errors, simulated background field populations and extinction distributions. These synthetic clusters are then dereddened in the same way as the real data, obtained from a deep near infrared survey, and used to derive the ages of 3 embedded clusters. Results were found to be consistent with those determined using spectrophotometric methods. This new method provides way to determine the ages of embedded clusters when only photometric data are available and there is strong differential reddening.

The Red MSX Source (RMS) survey has identified a sample of ~1200 massive young stellar objects (MYSOs), compact and ultra compact HII regions from a sample of ~2000 MSX and 2MASS colour selected sources. We have used the 100 m Green Bank telescope to search for 22-24 GHz water maser and ammonia (1,1), (2,2) and (3,3) emission towards ~600 RMS sources located within the northern Galactic plane. We have identified 308 H2O masers which corresponds to an overall detection rate of ~50%. Abridged: We detect ammonia emission towards 479 of these massive young stars, which corresponds to ~80%. Ammonia is an excellent probe of high density gas allowing us to measure key parameters such as gas temperatures, opacities, and column densities, as well as providing an insight into the gas kinematics. The average kinetic temperature, FWHM line width and total NH3 column density for the sample are approximately 22 K, 2 km/s and 2x10^{15} cm^{-2}, respectively. We find that the NH3 (1,1) line width and kinetic temperature are correlated with luminosity and finding no underlying dependence of these parameters on the evolutionary phase of the embedded sources, we conclude that the observed trends in the derived parameters are more likely to be due to the energy output of the central source and/or the line width-clump mass relationship. The velocities of the peak H2O masers and the NH3 emission are in excellent agreement with each other, which would strongly suggest an association between the dense gas and the maser emission. Moreover, we find the bolometric luminosity of the embedded source and the isotropic luminosity of the H2O maser are also correlated. We conclude from the correlations of the cloud and water maser velocities and the bolometric and maser luminosity that there is a strong dynamical relationship between the embedded young massive star and the H2O maser.

We present first results of the H_2O Southern Galactic Plane Survey (HOPS), using the Mopra radiotelescope with a broad band backend and a beam size of about 2'. We have observed 100 square degrees of the southern Galactic plane at 12mm (19.5 to 27.5GHz), including spectral line emission from water masers, multiple metastable transitions of ammonia, cyanoacetylene, methanol and radio recombination lines. In this paper, we report on the characteristics of the survey and water maser emission. We find 540 water masers, of which 334 are new detections. The strongest maser is 3933Jy and the weakest is 0.7Jy, with 62 masers over 100Jy. In 14 maser sites, the spread in velocity of the water maser emission exceeds 100km/s. In one region, the water maser velocities are separated by 351.3km/s. The rms noise levels are typically between 1-2Jy, with 95% of the survey under 2Jy. We estimate completeness limits of 98% at around 8.4Jy and 50% at around 5.5Jy. We estimate that there are between 800 and 1500 water masers in the Galaxy that are detectable in a survey with similar completeness limits to HOPS. We report possible masers in NH_3 (11,9) and (8,6) emission towards G19.61-0.23 and in the NH_3 (3,3) line towards G23.33-0.30.

There is currently no accepted theoretical framework for the formation of the most massive stars, and the manner in which protostars continue to accrete and grow in mass beyond \sim10Msun is still a controversial topic. In this study we use several prescriptions of stellar accretion and a description of the Galactic gas distribution to simulate the luminosities and spatial distribution of massive protostellar population of the Galaxy. We then compare the observables of each simulation to the results of the Red MSX Source (RMS) survey, a recently compiled database of massive young stellar objects. We find that the observations are best matched by accretion rates which increase as the protostar grows in mass, such as those predicted by the turbulent core and competitive accretion (i.e. Bondi-Hoyle) models. These 'accelerating accretion' models provide very good qualitative and quantitative fits to the data, though we are unable to distinguish between these two models on our simulations alone. We rule out models with accretion rates which are constant with time, and those which are initially very high and which fall away with time, as these produce results which are quantitatively and/or qualitatively incompatible with the observations. To simultaneously match the low- and high-luminosity YSO distribution we require the inclusion of a 'swollen-star' pre-main-sequence phase, the length of which is well-described by the Kelvin-Helmholz timescale. Our results suggest that the lifetime of the YSO phase is \sim 10^5yrs, whereas the compact Hii-region phase lasts between \sim 2 - 4 \times 10^5yrs depending on the final mass of the star. Finally, the absolute numbers of YSOs are best matched by a globally averaged star-formation rate for the Galaxy of 1.5-2Msun/yr.

The benefits for the e-MERLIN and EVN arrays of using antennae at the satellite communication station at Goonhilly in Cornwall are discussed. The location of this site - new to astronomy - will provide an almost equal distribution of long baselines in the east-west- and north-south directions, and opens up the possibility to get significantly improved observations of equatorial fields with e-MERLIN. These additional baselines will improve the sensitivity on a set of critical spatial scales and will increase the angular resolution of e-MERLIN by a factor of two. e-MERLIN observations, including many allocated under the e-MERLIN Legacy programme, will benefit from the enhanced angular resolution and imaging capability especially for sources close to or below the celestial equator (where ESO facilities such as ALMA will operate) of including the Goonhilly telescopes. Furthermore, the baselines formed between Goonhilly and the existing stations will close the gap between the baselines of e-MERLIN and those of the European VLBI Network (EVN) and therefore enhance the legacy value of e-MERLIN datasets.

A consortium of universities has recently been formed with the goal of using the decommissioned telecommunications infrastructure at the Goonhilly Earth Station in Cornwall, UK, for astronomical purposes. One particular goal is the introduction of one or more of the ~30-metre parabolic antennas into the existing e-MERLIN radio interferometer. This article introduces this scheme and presents some simulations which quantify the improvements that would be brought to the e-MERLIN system. These include an approximate doubling of the spatial resolution of the array, an increase in its N-S extent with strong implications for imaging the most well-studied equatorial fields, accessible to ESO facilities including ALMA. It also increases the overlap between the e-MERLIN array and the European VLBI Network. We also discuss briefly some niche science areas in which an e-MERLIN array which included a receptor at Goonhilly would be potentially world-leading, in addition to enhancing the existing potential of e-MERLIN in its role as a Square Kilometer Array pathfinder instrument.

We present a determination of the luminosity functions of massive young stellar objects (MYSOs) and compact (C)HII regions within the Milky Way Galaxy using the large, well-selected sample of these sources identified by the Red MSX Source (RMS) survey. The MYSO luminosity function decreases monotonically such that there are few with $L\gtrsim 10^{5}$Lsol, whilst the CHII regions are detected up to ~10$^{6}Lsol. The lifetimes of these phases are also calculated as a function of luminosity by comparison with the luminosity function for local main-sequence OB stars. These indicate that the MYSO phase has a duration ranging from 4x10$^{5}$ yrs for 10$^{4}$Lsol to ~7x10$^{4}$ yrs at 10$^{5}$Lsol, whilst the CHII region phase lasts of order 3x10$^{5}$ yrs or ~3-10% of the exciting star's main-sequence lifetime. MYSOs between 10$^{4} Lsol and ~10$^{5}$ Lsol are massive but do not display the radio continuum or near-IR \HI{} recombination line emission indicative of an HII region, consistent with being swollen due to high ongoing or recent accretion rates. Above ~10$^{5}$ Lsol the MYSO phase lifetime becomes comparable to the main-sequence Kelvin-Helmholtz timescale, at which point the central star can rapidly contract onto the main-sequence even if still accreting, and ionise a CHII region, thus explaining why few highly luminous MYSOs are observed.

We apply mid-infrared spectro-interferometry to the massive young stellar object CRL2136. The observations were performed with the Very Large Telescope Interferometer and the MIDI instrument at a 42m baseline probing angular scales of 50 milli-arcseconds. We model the observed visibilities in parallel with diffraction-limited images at both 24.5micron and in the N-band (with resolutions of 0.6" and 0.3", respectively), as well as the spectral energy distribution. The arcsec-scale spatial information reveals the well-resolved emission from the dusty envelope. By simultaneously modelling the spatial and spectral data, we find that the bulk of the dust emission occurs at several dust sublimation radii (approximately 170 AU). This reproduces the high mid-infrared fluxes and at the same time the low visibilities observed in the MIDI data for wavelengths longward of 8.5micron. However, shortward of this wavelength the visibility data show a sharp up-turn indicative of compact emission. We discuss various potential sources of this emission. We exclude a dust disk being responsible for the observed spectral imprint on the visibilities. A cool supergiant star and an accretion disk are considered and both shown to be viable origins of the compact mid-infrared emission. We propose that CRL2136 is embedded in a dusty envelope, which truncates at several times the dust sublimation radius. A dust torus is manifest in the equatorial region. We find that the spectro-interferometric N-band signal can be reproduced by either a gaseous disk or a bloated central star. If the disk extends to the stellar surface, it accretes at a rate of 3.0 10^(-3) Msun/yr.

The supersonic stellar and disk winds possessed by massive young stellar objects will produce shocks when they collide against the interior of a pre-existing bipolar cavity (resulting from an earlier phase of jet activity). The shock heated gas emits thermal X-rays which may be observable by spaceborne observa- tories such as the Chandra X-ray Observatory. Hydrodynamical models are used to explore the wind-cavity interaction. Radiative transfer calculations are performed on the simulation output to produce synthetic X-ray observations, allowing constraints to be placed on model parameters through comparisons with observations. The model reveals an intricate interplay between the inflowing and outflowing material and is successful in reproducing the observed X-ray count rates from massive young stellar objects.

Context: The Red MSX Source (RMS) survey is returning a large sample of massive young stellar objects (MYSOs) and ultra-compact (UC) \HII{} regions using follow-up observations of colour-selected candidates from the MSX point source catalogue. Aims: To obtain the bolometric fluxes and, using kinematic distance information, the luminosities for young RMS sources with far-infrared fluxes. Methods: We use a model spectral energy distribution (SED) fitter to obtain the bolometric flux for our sources, given flux data from our work and the literature. The inputs to the model fitter were optimised by a series of investigations designed to reveal the effect varying these inputs had on the resulting bolometric flux. Kinematic distances derived from molecular line observations were then used to calculate the luminosity of each source. Results: Bolometric fluxes are obtained for 1173 young RMS sources, of which 1069 have uniquely constrained kinematic distances and good SED fits. A comparison of the bolometric fluxes obtained using SED fitting with trapezium rule integration and two component greybody fits was also undertaken, and showed that both produce considerable scatter compared to the method used here. Conclusions: The bolometric flux results allowed us to obtain the luminosity distributions of YSOs and UC\HII{} regions in the RMS sample, which we find to be different. We also find that there are few MYSOs with L $\geq$ 10$^{5}$\lsol{}, despite finding many MYSOs with 10$^{4}$\lsol{} $\geq$ L $\geq$ 10$^{5}$\lsol{}.

The Red MSX Source (RMS) survey has identified a large sample of massive young stellar objects (MYSOs) and ultra compact (UC) HII regions from a sample of ~2000 MSX and 2MASS colour selected sources. Using a recent catalogue of molecular clouds derived from the Boston University-Five College Radio Astronomy Observatory Galactic Ring Survey (GRS), and by applying a Galactic scaleheight cut off of 120 pc, we solve the distance ambiguity for RMS sources located within 18\degr < |l| > 54\degr. These two steps yield kinematic distances to 291 sources out of a possible 326 located within the GRS longitude range. Combining distances and integrated fluxes derived from spectral energy distributions, we estimate luminosities to these sources and find that > 90% are indicative of the presence of a massive star. We find the completeness limit of our sample is ~10^4 Lsun, which corresponds to a zero age main sequence (ZAMS) star with a mass of ~12 Msun. Selecting only these sources, we construct a complete sample of 196 sources. Comparing the properties of the sample of young massive stars with the general population, we find the RMS-clouds are generally larger, more massive, and more turbulent. We examine the distribution of this sub-sample with respect to the location of the spiral arms and the Galactic bar and find them to be spatially correlated. We identify three significant peaks in the source surface density at Galactocentric radii of approximately 4, 6 and 8 kpc, which correspond to the proposed positions of the Scutum, Sagittarius and Perseus spiral arms, respectively. Fitting a scale height to the data we obtain an average value of ~29+-0.5 pc, which agrees well with other reported values in the literature, however, we note a dependence of the scale height on galactocentric radius with it increases from 30 pc to 45 pc between 2.5 and 8.5 kpc.

We present the results of a programme of scanning and mapping observations of astronomical masers and Jupiter designed to characterise the performance of the Mopra Radio Telescope at frequencies between 16-50 GHz using the 12-mm and 7-mm receivers. We use these observations to determine the telescope beam size, beam shape and overall telescope beam efficiency as a function of frequency. We find that the beam size is well fit by $\lambda$/$D$ over the frequency range with a correlation coefficient of ~90%. We determine the telescope main beam efficiencies are between ~48-64% for the 12-mm receiver and reasonably flat at ~50% for the 7-mm receiver. Beam maps of strong H$_2$O (22 GHz) and SiO masers (43 GHz) provide a means to examine the radial beam pattern of the telescope. At both frequencies the radial beam pattern reveals the presence of three components, a central `core', which is well fit by a Gaussian and constitutes the telescopes main beam, and inner and outer error beams. At both frequencies the inner and outer error beams extend out to approximately 2 and 3.4 times the full-width half maximum of the main beam respectively. Sources with angular sizes a factor of two or more larger than the telescope main beam will couple to the main and error beams, and therefore the power contributed by the error beams needs to be considered. From measurements of the radial beam power pattern we estimate the amount of power contained in the inner and outer error beams is of order one-fifth at 22 GHz rising slightly to one-third at 43 GHz.

We present high resolution (R~50,000) spectroastrometry over the CO 1st overtone bandhead of a sample of seven intermediate/massive young stellar objects. These are primarily drawn from the red MSX source (RMS) survey, a systematic search for young massive stars which has returned a large, well selected sample of such objects. The mean luminosity of the sample is approximately 5 times 10^4 L_\odot, indicating the objects typically have a mass of ~15 solar masses. We fit the observed bandhead profiles with a model of a circumstellar disc, and find good agreement between the models and observations for all but one object. We compare the high angular precision (0.2-0.8 mas) spectroastrometric data to the spatial distribution of the emitting material in the best-fitting models. No spatial signatures of discs are detected, which is entirely consistent with the properties of the best-fitting models. Therefore, the observations suggest that the CO bandhead emission of massive young stellar objects originates in small-scale disks, in agreement with previous work. This provides further evidence that massive stars form via disc accretion, as suggested by recent simulations.

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

Using near infrared UKIDSS Galactic Plane Survey data, we make extinction measurements to individual stars along the same line of sight as molecular clouds. Using an existing 3D extinction map of the inner Galaxy, that provides line of sight specific extinction-distance relationships, we convert the measured extinction of molecular clouds to a corresponding distance. These distances are derived independently from kinematic methods, typically used to derive distances to molecular clouds, and as such they have no near/far ambiguity. The near/far distance ambiguity has been resolved for 27 clouds and distances have been derived to 20 clouds. The results are found to be in good agreement with kinematic measurements to molecular clouds where the ambiguity has already been resolved, using HI self-absorption techniques.

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.

Context: The Red MSX Source (RMS) survey is a multi-wavelength campaign of follow-up observations of a colour-selected sample of candidate massive young stellar objects (MYSOs) in the galactic plane. This survey is returning the largest well-selected sample of MYSOs to date, while identifying other dust contaminant sources with similar mid-infrared colours including a large number of new ultra-compact (UC)HII regions. Aims:To measure the far-infrared (IR) flux, which lies near the peak of the spectral energy distribution (SED) of MYSOs and UCHII regions, so that, together with distance information, the luminosity of these sources can be obtained. Methods:Less than 50% of RMS sources are associated with IRAS point sources with detections at 60 micron and 100 micron, though the vast majority are visible in Spitzer MIPSGAL or IRAS Galaxy Atlas (IGA) images. However, standard aperture photometry is not appropriate for these data due to crowding of sources and strong spatially variable far-IR background emission in the galactic plane. A new technique using a 2-dimensional fit to the background in an annulus around each source is therefore used to obtain far-IR photometry for young RMS sources. Results:Far-IR fluxes are obtained for a total of 1113 RMS candidates identified as young sources. Of these 734 have flux measurements using IGA 60 micron and 100 micron images and 724 using MIPSGAL 70 micron images, with 345 having measurements in both data sets.

The circumstellar structure on 100 AU scales of the massive young stellar object W33A is probed using the VLTI and the MIDI instrument. N-band visibilities on 4 baselines are presented which are inconsistent with a spherically symmetric geometry. The visibility spectra and SED are simultaneously compared to 2D axi-symmetric dust radiative transfer models with a geometry including a rotationally flattened envelope and outflow cavities. We assume an O7.5 ZAMS star as the central source, consistent with the observed bolometric luminosity. The observations are also compared to models with and without (dusty and gaseous) accretion disks. A satisfactory model is constructed which reproduces the visibility spectra for each (u,v) point. It fits the silicate absorption, the mid-IR slope, the far-infrared peak, and the (sub)mm of the SED. It produces a 350 micron morphology consistent with observations. The 10 micron emission on 100 AU scales is dominated by the irradiated walls of the cavity sculpted by the outflow. The visibilities rule out the presence of dust disks with total (gas and dust) masses more than 0.01 Msun. However, optically thick accretion disks, interior to the dust sublimation radius, are allowed to accrete at rates equalling the envelope's mass infall rate (up to 10^(-3) Msun/yr) without substantially affecting the visibilities due to the extinction by the extremely massive envelope of W33A.

The Young Stellar Object (YSO) W33A is one of the best known examples of a massive star still in the process of forming. Here we present Gemini North ALTAIR/NIFS laser-guide star adaptive-optics assisted K-band integral-field spectroscopy of W33A and its inner reflection nebula. In our data we make the first detections of a rotationally-flattened outer envelope and fast bi-polar jet of a massive YSO at near-infrared wavelengths. The predominant spectral features observed are Br-gamma, H_2, and a combination of emission and absorption from CO gas. We perform a 3-D spectro-astrometric analysis of the line emission, the first study of its kind. We find that the object's Br-gamma emission reveals evidence for a fast bi-polar jet on sub-milliarcsecond scales, which is aligned with the larger-scale outflow. The hybrid CO features can be explained as a combination of hot CO emission arising in a disk close to the central star, while cold CO absorption originates in the cooler outer envelope. Kinematic analysis of these features reveals that both structures are rotating, and consistent with being aligned perpendicularly to both the ionised jet and the large-scale outflow. Assuming Keplerian rotation, we find that the circumstellar disk orbits a central mass of >10Msun, while the outer envelope encloses a mass of ~15Msun. Our results suggest a scenario of a central star accreting material from a circumstellar disk at the centre of a cool extended rotating torus, while driving a fast bi-polar wind. These results therefore provide strong supporting evidence for the hypothesis that the formation mechanism for high-mass stars is qualitatively similar to that of low-mass stars.

We determine the slope of the near infrared extinction power law (A$_{\lambda} \propto \lambda^{-\alpha}$) for 8 regions of the Galaxy between l$\sim27^{\circ}$ and $\sim100^{\circ}$. UKIDSS Galactic Plane Survey data are compared, in colour-colour space, with Galactic population synthesis model data reddened using a series of power laws and convolved through the UKIDSS filter profiles. Monte Carlo simulations allow us to determine the best fit value of $\alpha$ and evaluate the uncertainty. All values are consistent with each other giving an average extinction power law of $\alpha$=2.14$^{+0.04}_{-0.05}$. This is much steeper than most laws previously derived in the literature from colour excess ratios, which are typically between 1.6 and 1.8. We show that this discrepancy is due to an inappropriate choice of filter wavelength in conversion from colour excess ratios to $\alpha$ and that effective rather than isophotal wavelengths are more appropriate. In addition, curved reddening tracks, which depend on spectral type and filter system, should be used instead of straight vectors.

Context: The Red MSX Source (RMS) survey has identified a large sample of candidate massive young stellar objects (MYSOs) and ultra compact (UC) HII regions from a sample of ~2000 MSX and 2MASS colour selected sources. Aims: To search for H2O masers towards a large sample of young high mass stars and to investigate the statistical correlation of H2O masers with the earliest stages of massive star formation. Methods: We have used the Mopra Radio telescope to make position-switched observations towards ~500 UCHII regions and MYSOs candidates identified from the RMS survey and located between 190\degr < l < 30\degr. These observations have a 4$\sigma$ sensitivity of ~1 Jy and a velocity resolution of ~0.4 km/s.} Results: We have detected 163 H2O masers, approximately 75% of which were previously unknown. Comparing the maser velocities with the velocities of the RMS sources, determined from 13CO observations, we have identified 135 RMS-H2O maser associations, which corresponds to a detection rate of ~27%. Taking into account the differences in sensitivity and source selection we find our detection rate is in general agreement with previously reported surveys. Conclusions: We find similar detection rates for UCHII regions and MYSOs candidates, suggesting that the conditions needed for maser activity are equally likely in these two stages of the star formation process. Looking at the detection rate as a function of distance from the Galactic centre we find it significantly enhanced within the solar circle, peaking at ~37% between 6-7 kpc, which is consistent with previous surveys of UC HII regions, possibly indicating the presence of a high proportion of more luminous YSOs and HII regions.

In the young stellar object (YSO) phase of their lives, massive stars drive bi-polar molecular outflows. These outflows produce beautiful, often hourglass shaped, cavities. The central star possesses a powerful stellar wind (v ~ 2000 km s^-1), and possibly a dense equatorial disk wind (v ~ 400 km s^-1), which collide with the inner surface of the bi-polar cavity and produces hot (T ~ 10^5 - 10^8 K) shocked plasma. A reverse shock is formed at the point where the ram pressure between the preshock flow balances the thermal pressure of the postshock flow and provides a site for the acceleration of non-thermal particles to relativistic energies. Hydrodynamical models of the wind interaction, coupled with calculations of the non-thermal energy spectrum, are used to explore the observable synchrotron and gamma-ray emission from these objects.

2D axis-symmetric hydrodynamical simulations are presented which explore the interaction of stellar and disk winds with surrounding infalling cloud material. The star, and its accompanying disk, blow winds inside a cavity cleared out by an earlier jet. The collision of the winds with their surroundings generates shock heated plasma which reaches temperatures up to ~10^8 K. Attenuated X-ray spectra are calculated from solving the equation of radiative transfer along lines-of-sight. This process is repeated at various epochs throughout the simulations to examine the evolution of the intrinsic and attenuated flux. We find that the dynamic nature of the wind-cavity interaction fuels intrinsic variability in the observed emission on timescales of several hundred years. This is principally due to variations in the position of the reverse shock which is influenced by changes in the shape of the cavity wall. The collision of the winds with the cavity wall can cause clumps of cloud material to be stripped away. Mixing of these clumps into the winds mass-loads the flow and enhances the X-ray emission measure. The position and shape of the reverse shock plays a key role in determining the strength and hardness of the X-ray emission. In some models the reverse shock is oblique to much of the stellar and disk outflows, whereas in others it is closely normal over a wide range of polar angles. For reasonable stellar and disk wind parameters the integrated count rate and spatial extent of the intensity peak for X-ray emission agree with \textit{Chandra} observations of the deeply embedded MYSOs S106 IRS4, Mon R2 IRS3 A, and AFGL 2591.(abridged)

(Abridged) Context: The Red MSX Source (RMS) survey is an ongoing multi-wavelength observational programme designed to return a large, well-selected sample of massive young stellar objects (MYSOs). We have identified $\sim$2000 MYSO candidates located throughout the Galaxy by comparing the colours of MSX and 2MASS point sources to those of known MYSOs. Aims: To identify the populations of UCHII regions and PNe within the sample and examine their Galactic distribution. Method: We have conducted high resolution radio continuum observations at 6 cm towards 659 MYSO candidates in the northern hemisphere ($10\degr< l < 250\degr$) using the VLA. In addition to these targeted observations we present archival data towards a further 315 RMS sources extracted from a previous VLA survey of the inner Galaxy. Results: We find radio emission towards 272 ($\sim$27% of the observed sample). Using results from other parts of our multi-wavelength survey we separate these RMS-radio associations into two distinct types of objects, classifying 51 as PNe and a further 208 as either compact or UC HII regions. Using this well selected sample of HII regions we estimate their Galactic scale height to be 0.6\degr. Conclusions: Using radio continuum and archival data we have identified 79 PNe and 391 HII regions within the northern RMS catalogue. We estimate the total fraction of contamination by PNe in the RMS sample is of order 10%. The sample of HII regions is probably the best representation to date of the Galactic population of HII regions as a whole.

We discuss VLTI AMBER and MIDI interferometry in addition to single-dish Subaru observations of massive young stellar objects. The observations probe linear size scales between 10 to 1000 AU for the average distance of our sources.

Massive young stellar objects (MYSO) are surrounded by massive dusty envelopes. Our aim is to establish their density structure on scales of ~1000 AU, i.e. a factor 10 increase in angular resolution compared to similar studies performed in the (sub)mm. We have obtained diffraction-limited (0.6") 24.5 micron images of 14 well-known massive star formation regions with Subaru/COMICS. The images reveal the presence of discrete MYSO sources which are resolved on arcsecond scales. For many sources, radiative transfer models are capable of satisfactorily reproducing the observations. They are described by density powerlaw distributions (n(r) ~ r^(-p)) with p = 1.0 +/-0.25. Such distributions are shallower than those found on larger scales probed with single-dish (sub)mm studies. Other sources have density laws that are shallower/steeper than p = 1.0 and there is evidence that these MYSOs are viewed near edge-on or near face-on, respectively. The images also reveal a diffuse component tracing somewhat larger scale structures, particularly visible in the regions S140, AFGL 2136, IRAS 20126+4104, Mon R2, and Cep A. We thus find a flattening of the MYSO envelope density law going from ~10 000 AU down to scales of ~1000 AU. We propose that this may be evidence of rotational support of the envelope (abridged).

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 UKIDSS Galactic Plane Survey (GPS) is one of the five near infrared Public Legacy Surveys that are being undertaken by the UKIDSS consortium, using the Wide Field Camera on the United Kingdom Infrared Telescope. It is surveying 1868 sq.deg. of the northern and equatorial Galactic plane at Galactic latitudes -5<b<5 in the J, H and K filters and a ~200 sq.deg. area of the Taurus-Auriga-Perseus molecular cloud complex in these three filters and the 2.12 um (1-0) H_2 filter. It will provide data on ~2 billion sources. Here we describe the properties of the dataset and provide a user's guide for its exploitation. We also present brief Demonstration Science results from DR2 and from the Science Verification programme. These results illustrate how GPS data will frequently be combined with data taken in other wavebands to produce scientific results. The Demonstration Science includes studies of: (i) the star formation region G28.983-0.603, cross matching with Spitzer-GLIMPSE data to identify YSOs; (ii) the M17 nebula; (iii) H_2 emission in the rho Ophiuchi dark cloud; (iv) X-ray sources in the Galactic Centre; (v) external galaxies in the Zone of Avoidance; (vi) IPHAS-GPS optical-infrared spectrophotometric typing. (abridged).

We present diffraction limited (0.6") 24.5micron Subaru/COMICS images of the red supergiant mu Cep. We report the detection of a circumstellar nebula, that was not detected at shorter wavelengths. It extends to a radius of at least 6" in the thermal infrared. On these angular scales, the nebula is roughly spherical, in contrast, it displays a pronounced asymmetric morphology closer in. We simultaneously model the azimuthally averaged intensity profile of the nebula and the observed spectral energy distribution using spherical dust radiative transfer models. The models indicate a constant mass-loss process over the past 1000 years, for mass-loss rates a few times 10^(-7) Msun/yr. This work supports the idea that at least part of the asymmetries in shells of evolved massive stars and supernovae may be due to the mass-loss process in the red supergiant phase.

We present interferometric and single-dish mid-infrared observations of a sample of massive young stellar objects (BN-type objects), using VLTI-MIDI (10 micron) and Subaru-COMICS (24.5 micron). We discuss the regions S140, Mon R2, M8E-IR, and W33A. The observations probe the inner regions of the dusty envelope at scales of 50 milli arcsecond and 0.6 arcsec (100-1000 AU), respectively. Simultaneous model fits to spectral energy distributions and spatial data are achieved using self-consistent spherical envelope modelling. We conclude that those MYSO envelopes that are best described by a spherical geometry, the commensurate density distribution is a powerlaw with index -1.0. Such a powerlaw is predicted if the envelope is supported by turbulence on the 100-1000AU scales probed with MIDI and COMICS, but the role of rotation at these spatial scales need testing.

Context: The Red MSX Source (RMS) survey is an ongoing multi-wavelength observational programme designed to return a large, high-resolution mid-infrared colour-selected sample of massive young stellar objects. Aims: A critical part of our follow-up programme is to conduct 13CO molecular line observations in order to determine kinematic distances to all of our MYSO candidates. These distances will allow us to identify and remove nearby low-mass YSOs and help in identifying evolved stars which are weak CO emitters. Method: We have used the 15 m James Clerk Maxwell Telescope (JCMT), the 13.7 m telescope of the Purple Mountain Observatory (PMO), the 20 m Onsala telescope and the 22m Mopra telescope to conduct molecular line observations towards 508 MYSOs candidates located in the 1st and 2nd Quadrants. Results: We detect 13CO emission towards 780 RMS sources which corresponds to approximately 84% of those observed (911). A total of 2595 emission components are detected above 3sigma level (typically T^*_{\rm{A}} > 0.3K), with multiple components being observed towards the majority of these sources -- 520 sources (~56%) -- with an average of ~4 molecular clouds detected along each line of sight. We have used archival CS (J=2-1) and maser velocities to resolve the component multiplicity towards 175 sources (~20%) and have derived a criterion which is used to identify the most likely component for a further 191 multiple component sources. Combined with the single component detections we have obtained unambiguous kinematic velocities for 638 of the 780 MYSOs candidates towards which CO is detected (~80% of the detections). Using the Galactic rotation curve we calculate kinematic distances for all detected components.

We present an analysis of wind-blown, parsec-sized, mid-infrared bubbles and associated star-formation using GLIMPSE/IRAC, MIPSGAL/MIPS and MAGPIS/VLA surveys. Three bubbles from the Churchwell et al. (2006) catalog were selected. The relative distribution of the ionized gas (based on 20 cm emission), PAH emission (based on 8 um, 5.8 um and lack of 4.5 um emission) and hot dust (24 um emission) are compared. At the center of each bubble there is a region containing ionized gas and hot dust, surrounded by PAHs. We identify the likely source(s) of the stellar wind and ionizing flux producing each bubble based upon SED fitting to numerical hot stellar photosphere models. Candidate YSOs are also identified using SED fitting, including several sites of possible triggered star formation.

We present K-band polarimetric images of several massive young stellar objects at resolutions $\sim$ 0.1-0.5 arcsec. The polarization vectors around these sources are nearly centro-symmetric, indicating they are dominating the illumination of each field. Three out of the four sources show elongated low-polarization structures passing through the centers, suggesting the presence of polarization disks. These structures and their surrounding reflection nebulae make up bipolar outflow/disk systems, supporting the collapse/accretion scenario as their low-mass siblings. In particular, S140 IRS1 show well defined outflow cavity walls and a polarization disk which matches the direction of previously observed equatorial disk wind, thus confirming the polarization disk is actually the circumstellar disk. To date, a dozen massive protostellar objects show evidence for the existence of disks; our work add additional samples around MYSOs equivalent to early B-type stars.

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.

This chapter reviews progress in the field of massive star formation. It focuses on evidence for accretion and current models that invoke high accretion rates. In particular it is noted that high accretion rates will cause the massive young stellar object to have a radius much larger than its eventual main sequence radius throughout much of the accretion phase. This results in low effective temperatures which may provide the explanation as to why luminous young stellar objects do not ionized their surroundings to form ultra-compact H II regions. The transition to the ultra-compact H II region phase would then be associated with the termination of the high accretion rate phase. Objects thought to be in a transition phase are discussed and diagnostic diagrams to distinguish between massive young stellar objects and ultra-compact H II regions in terms of line widths and radio luminosity are presented.

The paper investigates the milli-arcsecond scale structure of the present-day mass-loss of the post-Red Supergiant IRC+10420. We use three telescopes of the VLT Interferometer in combination with the AMBER near-infrared beam combiner to measure spectrally dispersed correlated fluxes in the K-band around the Br gamma transition. The resulting visibilities are compared to the predicted visibilities of emission structures with various simple models in order to infer the size of the observed emission region. The Br gamma line is resolved by VLTI+AMBER on all three baselines, with the maximum projected baseline extending 69 meter and a P.A. ranging between 10 and 30 degrees. A differential phase between line and continuum is detected on the longest baseline. The Br gamma emission region is found to have a diameter of 3.3 milli-arcseconds (FWHM), when compared to a Gaussian intensity distribution. A uniform disk and a ring-like intensity distribution do not fit the line visibilities. Comparing the AMBER equivalent width of Br gamma with measurements from various epochs, we find that the stellar photosphere contributes about 60% of the total continuum light at 2.2 micron. The remaining 40% continuum emission is found on scales larger than the 66mas AMBER field of view. Using simple arguments, and assuming optically thick line emission, we find that the line emitting region is elongated. We briefly discuss the possibilities whether such a structure is due to a bi-polar flow or a circumstellar disk. (Abridged).

Here we describe the Red MSX Source (RMS) survey which is the largest, systematic, galaxy-wide search for massive young stellar objects (MYSOs) yet undertaken. Mid-IR bright point sources from the MSX satellite survey have been followed-up with ground-based radio, millimetre, and infrared observations to identify the contaminating sources and characterise the MYSOs and UCHII regions. With the initial classification now complete the distribution of sources in the galaxy will be discussed, as well as some programmes being developed to exploit our sample.

We report on resolved interferometric observations with VLTI/MIDI of the massive young stellar object (MYSO) W33A. The MIDI observations deliver spectrally dispersed visibilities with values between 0.03 and 0.06, for a baseline of 45m over the wavelength range 8-13 micron. The visibilities indicate that W33A has a FWHM size of approximately 120AU (0.030'') at 8 micron which increases to 240AU at 13 micron, scales previously unexplored among MYSOs. This observed trend is consistent with the temperature falling off with distance. 1D dust radiative transfer models are simultaneously fit to the visibility spectrum, the strong silicate feature and the shape of the mid infrared spectral energy distribution (SED). For any powerlaw density distribution, we find that the sizes (as implied by the visibilities) and the stellar luminosity are incompatible. A reduction to a third of W33A's previously adopted luminosity is required to match the visibilities; such a reduction is consistent with new high resolution 70 micron data from Spitzer's MIPSGAL survey. We obtain best fits for models with shallow dust density distributions of r^(-0.5) and r^(-1.0) and for increased optical depth in the silicate feature produced by decreasing the ISM ratio of graphite to silicates and using optical grain properties by Ossenkopf et al. (1992).

Abridged abstract: The Red MSX Source (RMS) survey is an ongoing effort to return a large, well-selected sample of massive young stellar objects (MYSOs) within our Galaxy. A series of ground-based follow-up observations are being undertaken in order to remove contaminant objects from our list of 2000 candidates, and to begin characterising these MYSOs. As a part of these follow-up observations, high resolution (~1") mid-IR imaging aids the identification of contaminant objects which are resolved (UCHII regions, PN) as opposed to those which are unresolved (YSOs, evolved stars) as well as identifying YSOs near UCHII regions and other multiple sources. We present 10.4 micron imaging observations for 346 candidate MYSOs in the RMS survey in the Southern Hemisphere, primarily outside the region covered by the GLIMPSE Spitzer Legacy Survey. These were obtained using TIMMI2 on the ESO 3.6m telescope in La Silla, Chile. Our photometric accuracy is of order 0.05Jy, and our astrometric accuracy is 0.8", which is an improvement over the nominal 2" accuracy of the MSX PSC.

Abridged: The Red MSX Source (RMS) survey is an ongoing multi-wavelength observational programme designed to return a large, well-selected sample of massive young stellar objects (MYSOs). Here we present 13CO observations made towards 854 MYSOs candidates located in the 3rd and 4th quadrants. We detected 13CO emission towards a total of 751 of the 854 RMS sources observed (~88%). In total 2185 emission components are detected above 3$\sigma$ level. Multiple emission profiles are observed towards the majority of these sources - 455 sources (~60%) - with an average of ~4 molecular clouds along the line of sight. These multiple emission features make it difficult to assign a kinematic velocity to many of our sample. We have used archival CS (J=2-1) and maser velocities to resolved the component multiplicity towards 82 sources and have derived a criterion which is used to identify the most likely component for a further 202 multiple component sources. Combined with the single component detections we have obtained unambiguous kinematic velocities towards 580 sources (~80% of the detections). The 171 sources for which we have not been able to determine the kinematic velocity will require additional line data. Using the rotation curve of Brand and Blitz (1993) and their radial velocities we calculate kinematic distances for all components detected.

We present high-resolution observations made with the Very Large Array (VLA) in its A configuration at frequencies between 5 and 43 GHz of a sample of five massive young stellar objects (YSOs): LkHa101, NGC2024-IRS2, S106-IR, W75N and S140-IRS1. The resolution varied from 0.04 arcsec (at 43 GHz) to 0.5 arcsec (at 5 GHz), corresponding to a linear resolution as high as 17 AU for our nearest source. A MERLIN observation of S106-IR at 23 GHz with 0.03-arcsec resolution is also presented. S106-IR and S140-IRS1 are elongated at 43 GHz perpendicular to their large scale bipolar outflows. This confirms the equatorial wind picture for these sources seen previously in MERLIN 5 GHz observations. The other sources are marginally resolved at 43 GHz. The spectral indices we derive for the sources in our sample range from +0.2 to +0.8, generally consistent with ionized stellar winds. We have modelled our sources as uniform, isothermal spherical winds, with LkHa101 and NGC2024-IRS2 yielding the best fits. However, in all cases our fits give wind temperatures of only 2000 to 5000 K, much less than the effective temperatures of main-sequence stars of the same luminosity, a result which is likely due to the clumpy nature of the winds.

We have surveyed a ~0.9-square-degree area of the W3 giant molecular cloud and star-forming region in the 850-micron continuum, using the SCUBA bolometer array on the James Clerk Maxwell Telescope. A complete sample of 316 dense clumps was detected with a mass range from around 13 to 2500 Msun. Part of the W3 GMC is subject to an interaction with the HII region and fast stellar winds generated by the nearby W4 OB association. We find that the fraction of total gas mass in dense, 850-micron traced structures is significantly altered by this interaction, being around 5% to 13% in the undisturbed cloud but ~25 - 37% in the feedback-affected region. The mass distribution in the detected clump sample depends somewhat on assumptions of dust temperature and is not a simple, single power law but contains significant structure at intermediate masses. This structure is likely to be due to crowding of sources near or below the spatial resolution of the observations. There is little evidence of any difference between the index of the high-mass end of the clump mass function in the compressed region and in the unaffected cloud. The consequences of these results are discussed in terms of current models of triggered star formation.

The Red MSX Source (RMS) survey is a multi-wavelength programme of follow-up observations designed to distinguish between genuine massive young stellar objects (MYSOs) and other embedded or dusty objects, such as ultra compact (UC) HII regions, evolved stars and planetary nebulae (PNe), from a sample of ~2000 MYSOs candidates. These were identified by comparing the colours of objects from the MSX and 2MASS point source catalogues to those of known MYSOs, in order to develop colour selection criteria which have been used to produced the RMS sample of MYSOs candidates. Our ultimate aim is to produce a large unbiased sample of MYSOs (~500) with complementary multi-wavelength data with which to study their properties. Here we report the results of 826 radio continuum and 346 mid-IR imaging observations carried out using the Australia Telescope Compact Array and TIMMI2 on the ESO 3.6 m telescope respectively. These observations are aimed at identifying and removing contaminating sources.

The Red MSX Source (RMS) survey is a multi-wavelength program of follow-up observations designed to distinguish between genuine massive young stellar objects (MYSOs) and other embedded or dusty objects, such as ultra compact (UC) HII regions, evolved stars and planetary nebulae (PNe). We have identified nearly 2000 massive YSOs candidates by comparing the colours of MSX and 2MASS point sources to those of known MYSOs. Unfortunately, there are several other types of embedded or dust enshrouded objects that have similar colours as MYSOs and contaminate our sample. Two sources of contamination are from UCHII regions and PNe, both of which can be identified from the radio emission emitted by their ionised nebulae. In order to identify UCHII regions and PNe that contaminate our sample we have conducted high resolution radio continuum observations at 3.6 and 6 cm of all southern MYSOs candidates ($235< l < 350$) using the Australia Telescope Compact Array (ATCA). These observations have a spatial resolution of $\sim$1--2\arcsec and typical image rms noise values of $\sim$0.3 mJy -- sensitive enough to detect a B0.5 star at the far side of the Galaxy. Of the 826 RMS sources observed we found 199 to be associated with radio emission, $\sim$25% of the sample. More significant than those sources associated with radio emission are the 627 RMS sources for which no radio emission was detected as these sources are still potential MYSOs. The Galactic distribution, morphologies and spectral indices of the radio sources associated with the RMS sources are consistent with these sources being UCHII regions. In addition to the 826 RMS sources observed we present observations of a further 190 fields. These observations were made towards MSX sources that passed cuts in earlier versions of the survey, but were later excluded.

As part of the Red MSX Source Survey of Massive Young Stellar Objects (MYSOs) we have conducted multi-wavelength follow up observations of the well-known object V645 Cygni. We present our data on this object, whose near-infrared spectrum is exceptional and place these in context with previous observations. Our observations of V645 Cyg included near/mid infrared imaging observations, 13CO 2-1 line observations and high signal-to-noise velocity resolved near-infrared spectroscopy. The spectrum shows P-Cygni hydrogen Brackett emission, consistent with a high velocity stellar wind. A red-shifted emission component to a number of near-IR emission lines was also uncovered. This is associated with a similar component in the H alpha line. V645 Cyg is also found to have variable CO first overtone bandhead emission. The data clearly indicate that the outflow of V645 Cyg is variable. The unidentified feature in a previously published optical spectrum is identified with a receding outflow at 2000 km per second. The nature of this feature, which is found in hydrogen and helium atomic lines and CO molecular lines remains a puzzle.

The discovery of the second equatorial ionized stellar wind from a massive young stellar object is reported. High resolution radio continuum maps of S140 IRS 1 reveal a highly elongated source that is perpendicular to the larger scale bipolar molecular outflow. This picture is confirmed by location of a small scale monopolar near-IR reflection nebula at the base of the blueshifted lobe. A second epoch of observations over a five year baseline show little ordered outward proper motion of clumps as would have been expected for a jet. A third epoch, taken only 50 days after the second, did show significant changes in the radio morphology. These radio properties can all be understood in the context of an equatorial wind driven by radiation pressure from the central star and inner disc acting on the gas in the surface layers of the disc as proposed by Drew et al. (1998). This equatorial wind system is briefly compared with the one in S106IR, and contrasted with other massive young stellar objects that drive ionized jets.

We present numerical radiation-hydrodynamic simulations of cometary HII regions for a number of champagne flow and bowshock models. For the champagne flow models we study smooth density distributions with both steep and shallow gradients. We also consider cases where the ionizing star has a strong stellar wind, and cases in which the star additionally has a proper motion within the ambient density gradient. We present simulated emission-measure maps and long-slit spectra of our results. Our numerical models are not tailored to any particular object but comparison with observations from the literature shows that, in particular, the models combining density gradients and stellar winds are able to account for both the morphology and radial velocity behavior of several observed cometary HII regions, such as the well-studied object G29.96-0.02.

We review the phenomenon of ultra-compact H II regions (UCHIIs) as a key phase in the early lives of massive stars. This most visible manifestation of massive star formation begins when the Lyman continuum output from the massive young stellar object becomes sufficient to ionize the surroundings from which it was born. Knowledge of this environment is gained through an understanding of the morphologies of UCHII regions and we examine the latest developments in deep radio and mid-IR imaging. SPITZER data from the GLIMPSE survey are an important new resource in which PAH emission and the ionizing stars can be seen. We review the role played by strong stellar winds from the central stars in sweeping out central cavities and causing the limb-brightened appearance. A range of evidence from velocity structure, proper motions, the molecular environment and recent hydrodynamical modeling indicates that cometary UCHII regions require a combination of champagne flow and bow shock motion. Finally, we discuss the class of hyper-compact H II regions or broad recombination line objects. They are likely to mark the transition soon after the breakout of the Lyman continuum radiation from the young star. Models for these objects are presented, including photo-evaporating disks and ionized accretion flows that are gravitationally trapped. Evolutionary scenarios tracing young massive stars passage through these ionized phases are discussed.

We present the initial results of a 3-mm spectral line survey towards 83 methanol maser selected massive star-forming regions. Here we report observations of the J=5-4 and 6-5 rotational transitions of methyl cyanide (CH3CN) and the J=1-0 transition of HCO+and H13CO+. CH3CN emission is detected in 58 sources (70 %) of our sample). We estimate the temperature and column density for 37 of these using the rotational diagram method. The temperatures we derive range from 28-166 K, and are lower than previously reported temperatures, derived from higher J transitions. We find that CH3CN is brighter and more commonly detected towards ultra-compact HII (UCHII) regions than towards isolated maser sources. Detection of CH3CN towards isolated maser sources strongly suggests that these objects are internally heated and that CH3CN is excited prior to the UCHII phase of massive star-formation. HCO+ is detected towards 82 sources (99 % of our sample), many of which exhibit asymmetric line profiles compared to H13CO+. Skewed profiles are indicative of inward or outward motions, however, we find approximately equal numbers of red and blue-skewed profiles among all classes. Column densities are derived from an analysis of the HCO+ and H13CO+ line profiles. 80 sources have mid-infrared counterparts: 68 seen in emission and 12 seen in absorption as `dark clouds'. Seven of the twelve dark clouds exhibit asymmetric HCO+ profiles, six of which are skewed to the blue, indicating infalling motions. CH3CN is also common in dark clouds, where it has a 90 % detection rate.

We investigate the use of HI data to resolve the near/far ambiguity in kinematic distances of massive young stellar object (MYSO) candidates. Kinematic distances were obtained from 13CO 1-0 spectral line observations with the Mopra Telescope towards 94 candidates selected from the Red MSX Source (RMS) survey in the fourth Galactic quadrant. HI data from the Southern Galactic Plane Survey (SPGS) was used in conjunction with the HI self-absorption technique to determine the near or far distance. We resolved the kinematic distance ambiguity to 70% of the sources. We can also simultaneously solve for any multiple line-of-sight component sources. We discuss the advantages and disadvantages of this technique in comparison with other methods, and also perform confidence checks on the reliability of using the HI self-absorption technique. We examined the projected location of these sources in both the Galactic plane and longitude-velocity diagrams to ascertain any recognisable spiral arm pattern. Although no obvious spiral pattern was found when compared to that proposed by Cordes & Lazio, far distance sources tended to lie on or near spiral arm loci. Near distance sources, however, had peculiar velocity uncertainties consistent with the separation between the spiral arms themselves. The longitude-velocity plot shows a more consistent picture, with tangent points of the spiral arm loci easily seen. We conclude that using the HI self-absorption technique to determine kinematic distance ambiguities is a quick and reliable method in most cases, with an 80% success rate in determining the correct designation of whether an object is at the near or far distance.

The role of the SKA in high resolution observations of the formation of low and high mass stars in the Galaxy is examined. The large collecting area will have a large impact on dynamical studies at high resolution using spectral lines. The SKA will allow dramatic progress in the investigation of the ionised, atomic and molecular phases over a range of evolutionary stages. Attention is focused here on the potential of radio recombination lines, H I 21 cm line and molecular Zeeman measurements. This combination will unveil the physical processes that drive jets and winds from young stars. Other important areas such as pre-biotic molecules, disc ablation in OB clusters and massive star formation in nearby galaxies are also briefly discussed.

We present near-IR speckle images of 21 massive Young Stellar Objects (YSOs) associated with outflows. The aim of this study is to search for sub-arcsecond reflection nebulae associated with the outflow cavity. We find that 6 of the massive YSOs show a conical nebula which can be interpreted in terms of reflected light from the dusty walls of the outflow cavity. In all cases, the small scale structures seen in our images are compared with outflow indicators found in the literature. No clear correlation is found between the presence of the reflection nebulosity and any property such as degree of embeddedness. We also note that 3 of the sources show close companions, one of them belonging also to the sample with conical nebula.

We have selected red MSX sources (RMS) that have the colours of massive young stellar objects (MYSOs). Our aim is to generate a large, systematically selected sample to address questions such as their luminosity function, lifetimes, clustering and triggering. Other objects such as UCHIIs, PN, PPN and AGB stars have similar IR colours and a large programme of ground-based follow-up observations is underway to identify and eliminate these from the sample of the red MSX sources. These include radio continuum observations, kinematic distances, ground-based mid-IR imaging, near-IR imaging and spectroscopy to distinguish. We report the progress of these campaigns on the 3000 candidates, with initial indications showing that a substantial fraction are indeed massive YSOs.

We use a Monte Carlo code to generate synthetic near-IR reflection nebulae that resemble those (normally associated with a bipolar outflow cavity) seen towards massive young stellar objects (YSOs). The 2D axi-symmetric calculations use an analytic expression for a flattened infalling rotating envelope with a bipolar cavity representing an outflow. We are interested in which aspects of the circumstellar density distribution can be constrained by servations of these reflection nebulae. We therefore keep the line of sight optical depth constant in the model grid, as this is often constrained independently by observations. It is found that envelopes with density distributions corresponding to mass infall rates of ~ 1E-4 M_sun/yr (for an envelope radius of 4700 AU) seen at an inclination angle of ~ 45 deg approximately reproduce the morphology and extension of the sub-arcsecond nebulae observed in massive YSOs. Based on the flux ratio between the approaching and receding lobe of the nebula, we can constrain the system inclination angle. The cavity opening angle is well constrained from the nebula opening angle. Our simulations indicate that to constrain the outflow cavity shape and the degree of flattening in the envelope, near-IR imaging with higher resolution and dynamic range than speckle imaging in 4m-class telescopes is needed. The radiative transfer code is also used to simulate the near-IR sub-arcsecond nebula seen in Mon R2 IRS3. We find indications of a shallower opacity law in this massive YSO than in the interstellar medium, or possibly a sharp drop in the envelope density distribution at distances of ~ 1000 AU from the illuminating source.

We present subarcsecond observations at 2.7 and 1.4 mm of a sample of massive young stellar objects made with the BIMA millimetre array. For most sources the continuum emission on the smallest scales at 2.7 mm is dominated by free-free emission from the stellar wind or jet. Strong emission at 1.4 mm shows the presence of significant dust associated with Cep A and GL 490 but our resolution is not sufficient to resolve any structure. The 2.7-mm emission from GL 490 is resolved but it is not clear whether we are seeing a single circumstellar disk or a secondary companion, although near-infrared data support the disk hypothesis. Estimates of the dust mass yield values of ~1--4 Msun within radii of 150 to 1000 AU.

We present observations and models of the behaviour of the HI and HeI lines between 1.6 and 2.2um in a small sample of compact HII regions. As in our previous papers on planetary nebulae, we find that the `pure' 1.7007um 4^3D-3^3P and 2.16475um 7^(3,1)G-4^(3,1)F HeI recombination lines behave approximately as expected as the effective temperature of the central exciting star(s) increases. However, the 2.058um 2^1P-2^1S HeI line does not behave as the model predicts, or as seen in planetary nebulae. Both models and planetary nebulae showed a decrease in the HeI 2^1P-2^1S/HI Br gamma ratio above an effective temperature of 40000K. The compact HII regions do not show any such decrease. The problem with this line ratio is probably due to the fact that the photoionisation model does not account correctly for the high densities seen in these HII regions, and that we are therefore seeing more collisional excitation of the 2^1P level than the model predicts. It may also reflect some deeper problem in the assumed model stellar atmospheres. In any event, although the normal HeI recombination lines can be used to place constraints on the temperature of the hottest star present, the HeI 2^1P-2^1S/HI Br gamma ratio should not be used for this purpose in either Galactic HII regions or in starburst galaxies, and conclusions from previous work using this ratio should be regarded with extreme caution. We also show that the combination of the near infrared `pure' recombination line ratios with mid-infrared forbidden line data provides a good discriminant of the form of the far ultraviolet spectral energy distribution of the exciting star(s). From this we conclude that CoStar models are a poor match to the available data for our sources, though the more recent WM-basic models are a better fit.

We present a series of JCMT, BIMA and VLA observations of the massive star-forming region associated with G35.2--0.7N. These new observations shed considerable light on the nature of the outflows in this region. The combination of our CO, SiO and radio data suggest that there are perhaps as many as four outflows emanating from the core containing G35.2--0.7N. The most striking result is that the radio jet is ruled out as the driving source of the large-scale CO flow. The jet drives its own highly-collimated CO outflow. Molecular line emission from H13CN and H13CO+ peaks towards the centre of the large-scale CO outflow, but no radio source is detected within. A third flow may be seen in SiO emission, which appears to be centred on another peak in H13CO+. The fourth flow is not seen in CO or SiO, but is inferred from previous L'-band data and the discovery of a radio source at the proposed centre of this flow.

The combination of mid-infrared data from the MSX satellite mission and ground-based near-infrared photometry is used to characterise the properties of the mid-infrared population of the Galactic plane. The colours of the youngest sources still heavily embedded within their natal molecular clouds are in general different from evolved stars shrouded within their own dust shells. Our main motivation is to use MSX for an unbiased search for a large (~1000) sample of massive young stellar objects (MYSOs). A simple analysis shows that the MSX point source catalogue should contain most of the MYSOs within our Galaxy. We develop colour selection criteria using combined near- and mid-infrared data for MYSOs, which produces a list of 3071 objects, excluding the galactic centre region. The programme of follow-up observations already underway to separate the MYSOs from compact H II regions and other remaining objects is briefly described. We also show that these data can be used, just as IRAS data has been previously, to provide a separation between evolved stars with carbon rich and oxygen rich dust. These data may also be used to search for evidence of dust around normal main sequence stars, such as low mass pre-main sequence stars or the Vega-excess class of objects where debris disks are presumed to remain from the planet formation process. We discuss the accuracy and completeness of the MSX point source catalogue, and show that the errors present tend to be of a kind that is not significant for the main stellar populations we discuss in this paper.

This paper continues our study of the behaviour of near infrared helium recombination lines in planetary nebula. We find that the 1.7007um 4^3D-3^3P HeI line is a good measure of the HeI recombination rate, since it varies smoothly with the effective temperature of the central star. We were unable to reproduce the observed data using detailed photoionisation models at both low and high effective temperatures, but plausible explanations for the difference exist for both. We therefore conclude that this line could be used as an indicator of the effective temperature in obscured nebula. We also characterised the nature of the molecular hydrogen emission present in a smaller subset of our sample. The results are consistent with previous data indicating that ultraviolet excitation rather than shocks is the main cause of the molecular hydrogen emission in planetary nebulae.

We have observed a large sample of compact planetary nebulae in the near-infrared to determine how the 2^1P-2^1S HeI line at 2.058um varies as a function of stellar effective temperature, Teff. The ratio of this line with HI Br g at 2.166um has often been used as a measure of the highest Teff present in a stellar cluster, and hence on whether there is a cut-off in the stellar initial mass function at high masses. However, recent photoionisation modelling has revealed that the behaviour of this line is more complex than previously anticipated. Our work shows that in most aspects the photoionisation models are correct. In particular, we confirm the weakening of the 2^1P-2^1S as Teff increases beyond 40000K. However, in many cases the model underpredicts the observed ratio when we consider the detailed physical conditions in the individual planetary nebulae. Furthermore, there is evidence that there is still significant 2^1P-2^1S HeI line emission even in the planetary nebulae with very hot (Teff>100000K) central stars. It is clear from our work that this ratio cannot be considered as a reliable measure of effective temperature on its own.

We present a new observation of the compact HII region, G29.96-0.02, that allows us to compare the velocity structure in the ionised gas and surrounding molecular gas directly. This allows us to remove most of the remaining ambiguity about the nature of this source. In particular, the comparison of the velocity structure present in the 4S-3P HeI lines with that found in the 1-0 S(1) of molecular hydrogern convincingly rules out a bow shock as being important to the kinematics of this source. Our new observation therefore agrees with our previous conclusion, drawn from a velocity resolved HI Br gamma map, that most of the velocity structure in G29.96-0.02 can largely be explained as a result of a champagne flow model. We also find that the best simple model must invoke a powerful stellar wind to evacuate the `head' of the cometary HII region of ionised gas. However, residual differences between model and data tend to indicate that no single simple model can adequately explain all the observed features.

We have mapped the velocity structure of the cometary compact HII region G29.96-0.02 using long-slit echelle spectra of the HI Br gamma line. This technique detects line emission over a much wider area at the necessary spatial resolution compared to radio recombination line observations. Significant structure in both the velocity centroids and the line widths is seen over the entire nebula. Large line widths are seen ahead of the bow and in the tail which may be due to turbulent motions in shocked and interface regions respectively. We construct analytic models of the density and velocity structure in order to attempt to distinguish between the bow shock and champagne flow models which have been put forward to explain the cometary morphology of many compact HII regions. The bow shock model is unable to explain the large velocity gradient that we see right across the tail of the cometary region which can only be explained by the streaming motions towards low density regions in the champagne model. However, our approximation to the champagne model is also not able to fit all of the features of the data. More realistic versions of this model which include the effects of stellar winds and density gradients may be able to provide a better match to these data.