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.

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.

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.

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.

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.

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.

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

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.

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.

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.