We present initial results from the K-band focal plane array Examinations of Young STellar Object Natal Environments (KEYSTONE) survey, a large project on the 100-m Green Bank Telescope mapping ammonia emission across eleven giant molecular clouds at distances of $0.9-3.0$ kpc (Cygnus X North, Cygnus X South, M16, M17, MonR1, MonR2, NGC2264, NGC7538, Rosette, W3, and W48). This data release includes the NH$_3$ (1,1) and (2,2) maps for each cloud, which are modeled to produce maps of kinetic temperature, centroid velocity, velocity dispersion, and ammonia column density. Median cloud kinetic temperatures range from $11.4\pm2.2$ K in the coldest cloud (MonR1) to $23.0\pm6.5$ K in the warmest cloud (M17). Using dendrograms on the NH$_3$ (1,1) integrated intensity maps, we identify 856 dense gas clumps across the eleven clouds. Depending on the cloud observed, $40-100\%$ of the clumps are aligned spatially with filaments identified in H$_2$ column density maps derived from SED-fitting of dust continuum emission. A virial analysis reveals that 523 of the 835 clumps ($\sim63\%$) with mass estimates are bound by gravity alone. We find no significant difference between the virial parameter distributions for clumps aligned with the dust-continuum filaments and those unaligned with filaments. In some clouds, however, hubs or ridges of dense gas with unusually high mass and low virial parameters are located within a single filament or at the intersection of multiple filaments. These hubs and ridges tend to host water maser emission, multiple 70$\mu$m-detected protostars, and have masses and radii above an empirical threshold for forming massive stars.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

In order to understand the initial conditions and early evolution of star formation in a wide range of Galactic environments, we carried out an investigation of 64 \textit{Planck} Galactic Cold Clumps (PGCCs) in the second quadrant of the Milky Way. Using the $^{13}$CO and C$^{18}$O $J = 1 - 0$ lines, and 850\,$\mu$m continuum observations, we investigated cloud fragmentation and evolution associated with star formation. We extracted 468 clumps and 117 cores from the $^{13}$CO line and 850\,$\mu$m continuum maps, respectively. We make use of the Bayesian Distance Calculator and derived the distances of all 64 PGCCs. We found that in general, the mass-size plane follows a relation of $m\sim r^{1.67}$. At a given scale, the masses of our objects are around 1/10 of that of typical Galactic massive star-forming regions. Analysis of the clump and core masses, virial parameters, densities, and mass-size relation suggests that the PGCCs in our sample have a low core formation efficiency ($\sim$3.0\%), and most PGCCs are likely low-mass star-forming candidates. Statistical study indicates that the 850\,$\mu$m cores are more turbulent, more optically thick, and denser than the $^{13}$CO clumps for star formation candidates, suggesting that the 850\,$\mu$m cores are likely more appropriate future star-formation candidates than the $^{13}$CO clumps.

We have investigated and applied machine-learning algorithms for Infrared Colour Selection of Galactic Wolf-Rayet (WR) candidates. Objects taken from the GLIMPSE catalogue of the infrared objects in the Galactic plane can be classified into different stellar populations based on the colours inferred from their broadband photometric magnitudes ($J$, $H$ and $K_s$ from 2MASS, and the four \textit{Spitzer}/IRAC bands). The algorithms tested in this pilot study are variants of the $k$-Nearest Neighbours ($k$-NN) approach, which is ideal for exploratory studies of classification problems where interrelations between variables and classes are complicated. The aims of this study are (1) to provide an automated tool to select reliable WR candidates and potentially other classes of objects, (2) to measure the efficiency of infrared colour selection at performing these tasks and, (3) to lay the groundwork for statistically inferring the total number of WR stars in our Galaxy. We report the performance results obtained over a set of known objects and selected candidates for which we have carried out follow-up spectroscopic observations, and confirm the discovery of 4 new WR stars.

{Context}. The HIFI instrument on the Herschel Space Observatory performed over 9100 astronomical observations, almost 900 of which were calibration observations in the course of the nearly four-year Herschel mission. The data from each observation had to be converted from raw telemetry into calibrated products and were included in the Herschel Science Archive. {Aims}. The HIFI pipeline was designed to provide robust conversion from raw telemetry into calibrated data throughout all phases of the HIFI missions. Pre-launch laboratory testing was supported as were routine mission operations. {Methods}. A modular software design allowed components to be easily added, removed, amended and/or extended as the understanding of the HIFI data developed during and after mission operations. {Results}. The HIFI pipeline processed data from all HIFI observing modes within the Herschel automated processing environment as well as within an interactive environment. The same software can be used by the general astronomical community to reprocess any standard HIFI observation. The pipeline also recorded the consistency of processing results and provided automated quality reports. Many pipeline modules were in use since the HIFI pre-launch instrument level testing. {Conclusions}. Processing in steps facilitated data analysis to discover and address instrument artefacts and uncertainties. The availability of the same pipeline components from pre-launch throughout the mission made for well-understood, tested, and stable processing. A smooth transition from one phase to the next significantly enhanced processing reliability and robustness.

The Wolf-Rayet (WR) nebula NGC3199 has been suggested to be a bow shock around its central star WR18, presumably a runaway star, because optical images of the nebula show a dominating arc of emission south-west of the star. We present the XMM-Newton detection of extended X-ray emission from NGC3199, unveiling the powerful effect of the fast wind from WR18. The X-ray emission is brighter in the region south-east of the star and analysis of the spectral properties of the X-ray emission reveals abundance variations: i) regions close to the optical arc present nitrogen-rich gas enhanced by the stellar wind from WR18 and ii) gas at the eastern region exhibits abundances close to those reported for nebular abundances derived from optical studies, signature of an efficient mixing of the nebular material with the stellar wind. The dominant plasma temperature and electron density are estimated to be $T\approx1.2\times$10$^{6}$ K and $n_\mathrm{e}$=0.3 cm$^{-3}$ with an X-ray luminosity in the 0.3-3.0 keV energy range of $L_\mathrm{X}$=2.6$\times$10$^{34}$ erg s$^{-1}$. Combined with information derived from Herschel and the recent Gaia first data release, we conclude that WR18 is not a runaway star and the formation, chemical variations, and shape of NGC3199 depend on the initial configuration of the interstellar medium.

Filaments are key for star formation models. As part of the study carried out by the Herschel GCC Programme, here we study the filament properties presented in GCC.VII in context with theoretical models of filament formation and evolution. A conservative sample of filaments at a distance D<500pc was extracted with the Getfilaments algorithm. Their physical structure was quantified according to two main components: the central (Gaussian) region (core component), and the power-law like region dominating the filament column density profile at larger radii (wing component). The properties and behaviour of these components relative to the total linear mass density of the filament and its environmental column density were compared with theoretical models describing the evolution of filaments under gravity-dominated conditions. The feasibility of a transition to supercritical state by accretion is dependent on the combined effect of filament intrinsic properties and environmental conditions. Reasonably self-gravitating (high Mline-core) filaments in dense environments (av\sim3mag) can become supercritical in timescales of t\sim1Myr by accreting mass at constant or decreasing width. The trend of increasing Mline-tot (Mline-core and Mline-wing), and ridge Av with background also indicates that the precursors of star-forming filaments evolve coevally with their environment. The simultaneous increase of environment and filament Av explains the association between dense environments and high Mline-core values, and argues against filaments remaining in constant single-pressure equilibrium states. The simultaneous growth of filament and background in locations with efficient mass assembly, predicted in numerical models of collapsing clouds, presents a suitable scenario for the fulfillment of the combined filament mass-environment criterium that is in quantitative agreement with Herschel observations.

New high-resolution far-infrared (FIR) observations of both ortho- and para-NH3 transitions toward IRC+10216 were obtained with Herschel, with the goal of determining the ammonia abundance and constraining the distribution of NH3 in the envelope of IRC+10216. We used the Heterodyne Instrument for the Far Infrared (HIFI) on board Herschel to observe all rotational transitions up to the J=3 level (three ortho- and six para-NH3 lines). We conducted non-LTE multilevel radiative transfer modelling, including the effects of near-infrared (NIR) radiative pumping through vibrational transitions. We found that NIR pumping is of key importance for understanding the excitation of rotational levels of NH3. The derived NH3 abundances relative to molecular hydrogen were (2.8+-0.5)x10^{-8} for ortho-NH3 and (3.2^{+0.7}_{-0.6})x10^{-8} for para-NH3, consistent with an ortho/para ratio of 1. These values are in a rough agreement with abundances derived from the inversion transitions, as well as with the total abundance of NH3 inferred from the MIR absorption lines. To explain the observed rotational transitions, ammonia must be formed near to the central star at a radius close to the end of the wind acceleration region, but no larger than about 20 stellar radii (1 sigma confidence level).

The radiative feedback of massive stars on molecular clouds creates pillars, globules and other features at the interface between the HII region and molecular cloud. We present here Herschel observations between 70 and 500 micron of the immediate environment of the Cygnus OB2 association, performed within the HOBYS program. All structures were detected based on their appearance at 70 micron, and have been classified as pillars, globules, evaporating gasous globules (EGGs), proplyd-like objects, and condensations. From the 70 and 160 micron flux maps, we derive the local FUV field on the PDR surfaces. In parallel, we use a census of the O-stars to estimate the overall FUV-field, that is 10^3-10^4 G_0 close to the central OB cluster (within 10 pc) and decreases down to a few tens G_0, in a distance of 50 pc. From a SED fit to the four longest Herschel wavelengths, we determine column density and temperature maps and derive masses, volume densities and surface densities for these structures. We find that the morphological classification corresponds to distinct physical properties. Pillars and globules have the longest estimated photoevaporation lifetimes, a few 10^6 yr, while all other features should survive less than that. These lifetimes are consistent with that found in simulations of turbulent, UV-illuminated clouds. We propose a tentative evolutionary scheme in which pillars can evolve into globules, which in turn then evolve into EGGs, condensations and proplyd-like objects.

The association of filaments with protostellar objects has made these structures a priority target in star formation studies. The datasets of the Herschel Galactic Cold Cores Key Programme allow for a statistical study of filaments with a wide range of intrinsic and environmental characteristics. Characterisation of this sample can be used to identify key physical parameters and quantify the role of environment in the formation of supercritical filaments. Filaments were extracted from fields at D<500pc with the getfilaments algorithm and characterised according to their column density profiles and intrinsic properties. Each profile was fitted with a beam-convolved Plummer-like function and quantified based on the relative contributions from the filament 'core', represented by a Gaussian, and 'wing' component, dominated by the power-law of the Plummer-like function. These parameters were examined for populations associated with different background levels. We find that filaments increase their core (Mcore) and wing (Mwing) contributions while increasing their total linear mass density (Mtot). Both components appear to be linked to the local environment, with filaments in higher backgrounds having systematically more massive Mcore and Mwing. This dependence on the environment supports an accretion-based model for filament evolution in the local neighbourhood (D<500pc). Structures located in the highest backgrounds develop the highest central Av, Mcore, and Mwing as Mtot increases with time, favoured by the local availability of material and the enhanced gravitational potential. Our results indicate that filaments acquiring a significantly massive central region with Mcore>Mcrit/2 may become supercritical and form stars. This translates into a need for filaments to become at least moderately self-gravitating in order to undergo localised star formation or become star-forming filaments.

The total population of Wolf-Rayet (WR) stars in the Galaxy is predicted by models to be as many as $\sim$6000 stars, and yet the number of catalogued WR stars as a result of optical surveys was far lower than this ($\sim$200) at the turn of this century. When beginning our WR searches using infrared techniques it was not clear whether WR number predictions were too optimistic or whether there was more hidden behind interstellar and circumstellar extinction. During the last decade we pioneered a technique of exploiting the near- and mid-infrared continuum colours for individual point sources provided by large-format surveys of the Galaxy, including 2MASS and Spitzer/GLIMPSE, to pierce through the dust and reveal newly discovered WR stars throughout the Galactic Plane. The key item to the colour discrimination is via the characteristic infrared spectral index produced by the strong winds of the WR stars, combined with dust extinction, which place WR stars in a relatively depopulated area of infrared colour-colour diagrams. The use of the Spitzer/GLIMPSE 8$\mu$m and, more recently, WISE 22$\mu$m fluxes together with cross-referencing with X-ray measurements in selected Galactic regions have enabled improved candidate lists that increased our confirmation success rate, achieved via follow-up infrared and optical spectroscopy. To date a total of 102 new WR stars have been found with many more candidates still available for follow-up. This constitutes an addition of $\sim$16\% to the current inventory of 642 Galactic WR stars. In this talk we review our methods and provide some new results and a preliminary review of their stellar and interstellar medium environments. We provide a roadmap for the future of this search, including statistical modeling, and what we can add to star formation and high mass star evolution studies.

The W3 GMC is a prime target for investigating the formation of high-mass stars and clusters. This second study of W3 within the HOBYS Key Program provides a comparative analysis of subfields within W3 to further constrain the processes leading to the observed structures and stellar population. Probability density functions (PDFs) and cumulative mass distributions (CMDs) were created from dust column density maps, quantified as extinction Av. The shape of the PDF, typically represented with a lognormal function at low Av "breaking" to a power-law tail at high Av, is influenced by various processes including turbulence and self-gravity. The breaks can also be identified, often more readily, in the CMDs. The PDF break from lognormal (Av(SF)= 6-10 mag) appears to shift to higher Av by stellar feedback, so that high-mass star-forming regions tend to have higher PDF breaks. A second break at Av > 50 mag traces structures formed or influenced by a dynamic process. Because such a process has been suggested to drive high-mass star formation in W3, this second break might then identify regions with potential for hosting high-mass stars/clusters. Stellar feedback appears to be a major mechanism driving the local evolution and state of regions within W3. A high initial star formation efficiency in a dense medium could result in a self-enhancing process, leading to more compression and favourable star-formation conditions (e.g., colliding flows), a richer stellar content, and massive stars. This scenario would be compatible with the "convergent constructive feedback" model introduced in our previous Herschel study.

CONTEXT: Asymptotic giant branch (AGB) stars are in one of the latest evolutionary stages of low to intermediate-mass stars. Their vigorous mass loss has a significant effect on the stellar evolution, and is a significant source of heavy elements and dust grains for the interstellar medium. The mass-loss rate can be well traced by carbon monoxide (CO) line emission. AIMS: We present new Herschel HIFI and IRAM 30m telescope CO line data for a sample of 53 galactic AGB stars. The lines cover a fairly large range of excitation energy from the $J=1\to0$ line to the $J=9\to8$ line, and even the $J=14\to13$ line in a few cases. We perform radiative transfer modelling for 38 of these sources to estimate their mass-loss rates. METHODS: We used a radiative transfer code based on the Monte Carlo method to model the CO line emission. We assume spherically symmetric circumstellar envelopes that are formed by a constant mass-loss rate through a smoothly accelerating wind. RESULTS: We find models that are consistent across a broad range of CO lines for most of the stars in our sample, i.e., a large number of the circumstellar envelopes can be described with a constant mass-loss rate. We also find that an accelerating wind is required to fit, in particular, the higher-J lines and that a velocity law will have a significant effect on the model line intensities. The results cover a wide range of mass-loss rates ($\sim 10^{-8}$ to $2\times 10^{-5}~\mathrm{M}_\odot~\mathrm{ yr}^{-1}$) and gas expansion velocities (2 to $21.5$ km s$^{-1}$), and include M-, S-, and C-type AGB stars. Our results generally agree with those of earlier studies, although we tend to find slightly lower mass-loss rates by about 40%, on average. We also present "bonus" lines detected during our CO observations.

Context. The high intensities and narrow bandwidths exhibited by some astronomical masers make them ideal tools for studying star-forming giant molecular clouds. The water maser transition $J_{K_{-1}K_{1}}=5_{32}-4_{41}$ at 620.701 GHz can only be observed from above Earth's strongly absorbing atmosphere; its emission has recently been detected from space. Aims. We sought to further characterize the star-forming environment of Orion KL by investigating the linear polarisation of a source emitting a narrow 620.701 GHz maser feature with the heterodyne spectrometer HIFI on board the Herschel Space Observatory. Methods. High-resolution spectral datasets were collected over a thirteen month period beginning in 2011 March, to establish not only the linear polarisation but also the temporal variability of the source. Results. Within a $3\sigma$ uncertainty, no polarisation was detected to an upper limit of approximately 2%. These results are compared with coeval linear polarisation measurements of the 22.235 GHz $J_{K_{-1}K_{1}}=6_{16}-5_{23}$ maser line from the Effelsberg 100-m radio telescope, typically a much stronger maser transition. Although strongly polarised emission is observed for one component of the 22.235 GHz maser at 7.2 km s$^{-1}$, a weaker component at the same velocity as the 620.701 GHz maser at 11.7 km s$^{-1}$ is much less polarised.

We present the activities carried out to calibrate and characterise the performance of the elements of attitude control and measurement on board the Herschel spacecraft. The main calibration parameters and the evolution of the indicators of the pointing performance are described, from the initial values derived from the observations carried out in the performance verification phase to those attained in the last year and half of mission, an absolute pointing error around or even below 1 arcsec, a spatial relative pointing error of some 1 arcsec and a pointing stability below 0.2 arsec. The actions carried out at the ground segment to improve the spacecraft pointing measurements are outlined. On-going and future developments towards a final refinement of the Herschel astrometry are also summarised. A brief description of the different components of the attitude control and measurement system (both in the space and in the ground segments) is also given for reference. We stress the importance of the cooperation between the different actors (scientists, flight dynamics and systems engineers, attitude control and measurement hardware designers, star-tracker manufacturers, etc.) to attain the final level of performance.

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

We present the results of Herschel HOBYS photometric mapping combined with BIMA observations and additional archival data, and perform an in-depth study of the evolutionary phases of the star-forming clumps in W 48A and their surroundings. Age estimates for the compact sources were derived from bolometric luminosities and envelope masses, which were obtained from the dust continuum emission, and agree within an order of magnitude with age estimates from molecular line and radio data. The clumps in W 48A are linearly aligned by age (east-old to west-young): we find a ultra compact (UC) HII region, a young stellar object (YSO) with class II methanol maser emission, a YSO with a massive outflow, and finally the NH_2D prestellar cores from Pillai et al. This remarkable positioning reflects the (star) formation history of the region. We find that it is unlikely that the star formation in the W 48A molecular cloud was triggered by the UCHII region and discuss the Aquila supershell expansion as a mayor influence on the evolution of W 48A. We conclude that the combination of Herschel continuum data with interferometric molecular line and radio continuum data is important to derive trustworthy age estimates and interpret the origin of large scale structures through kinematic information.

Ionization feedback should impact the probability distribution function (PDF) of the column density around the ionized gas. We aim to quantify this effect and discuss its potential link to the Core and Initial Mass Function (CMF/IMF). We used in a systematic way Herschel column density maps of several regions observed within the HOBYS key program: M16, the Rosette and Vela C molecular cloud, and the RCW 120 H ii region. We fitted the column density PDFs of all clouds with two lognormal distributions, since they present a double-peak or enlarged shape in the PDF. Our interpretation is that the lowest part of the column density distribution describes the turbulent molecular gas while the second peak corresponds to a compression zone induced by the expansion of the ionized gas into the turbulent molecular cloud. The condensations at the edge of the ionized gas have a steep compressed radial profile, sometimes recognizable in the flattening of the power-law tail. This could lead to an unambiguous criterion able to disentangle triggered from pre-existing star formation. In the context of the gravo-turbulent scenario for the origin of the CMF/IMF, the double peaked/enlarged shape of the PDF may impact the formation of objects at both the low-mass and the high-mass end of the CMF/IMF. In particular a broader PDF is required by the gravo-turbulent scenario to fit properly the IMF with a reasonable initial Mach number for the molecular cloud. Since other physical processes (e.g. the equation of state and the variations among the core properties) have already been suggested to broaden the PDF, the relative importance of the different effects remains an open question.

Pillars and globules are present in many high-mass star-forming regions, such as the Eagle nebula (M16) and the Rosette molecular cloud, and understanding their origin will help characterize triggered star formation. The formation mechanisms of these structures are still being debated. Recent numerical simulations have shown how pillars can arise from the collapse of the shell in on itself and how globules can be formed from the interplay of the turbulent molecular cloud and the ionization from massive stars. The goal here is to test this scenario through recent observations of two massive star-forming regions, M16 and Rosette. The column density structure of the interface between molecular clouds and H ii regions was characterized using column density maps obtained from far-infrared imaging of the Herschel HOBYS key programme. Then, the DisPerSe algorithm was used on these maps to detect the compressed layers around the ionized gas and pillars in different evolutionary states. Finally, their velocity structure was investigated using CO data, and all observational signatures were tested against some distinct diagnostics established from simulations. The column density profiles have revealed the importance of compression at the edge of the ionized gas. The velocity properties of the structures, i.e. pillars and globules, are very close to what we predict from the numerical simulations. We have identified a good candidate of a nascent pillar in the Rosette molecular cloud that presents the velocity pattern of the shell collapsing on itself, induced by a high local curvature. Globules have a bulk velocity dispersion that indicates the importance of the initial turbulence in their formation, as proposed from numerical simulations. Altogether, this study re-enforces the picture of pillar formation by shell collapse and globule formation by the ionization of highly turbulent clouds.

To gain insight into the physical conditions and kinematics of the warm (100-1000 K) gas around the red hyper-giant VY CMa, we performed sensitive high spectral resolution observations of molecular lines in the sub-mm/FIR using the HIFI instrument of the Herschel Space Observatory. We observed CO, H2O, and other molecular species, sampling excitation energies from a few tens to a few thousand K. These observations are part of the Herschel Guaranteed Time Key Program HIFISTARS. We detected the J=6-5, J=10-9, and J=16-15 lines of 12CO and 13CO at about 100, 300, and 750K above the ground state (and the 13CO J=9-8 line). These lines are crucial for improving the modelling of the internal layers of the envelope around VY CMa. We also detected 27 lines of H2O and its isotopomers, and 96 lines of species such as NH3, SiO, SO, SO2 HCN, OH and others, some of them originating from vibrationally excited levels. Three lines were not unambiguously assigned. Our observations confirm that VY CMa's envelope must consist of two or more detached components. The molecular excitation in the outer layers is significantly lower than in the inner ones, resulting in strong self-absorbed profiles in molecular lines that are optically thick in this outer envelope, for instance, low-lying lines of H2O. Except for the most abundant species, CO and H2O, most of the molecular emission detected at these sub-mm/FIR wavelengths arise from the central parts of the envelope. The spectrum of VY CMa is very prominent in vibrationally excited lines, which are caused by the strong IR pumping present in the central regions. Compared with envelopes of other massive evolved stars, VY CMa's emission is particularly strong in these vibrationally excited lines, as well as in the emission from less abundant species such as H13CN, SO, and NH3.

Determining the Galactic distribution and numbers of massive stars, such as Wolf-Rayet stars (WRs), is hampered by intervening Galactic or local circumstellar dust obscuration. In order to probe such regions of the Galaxy we can use infrared observations, which provide a means for finding such hidden populations through the dust. The availability of both 2MASS and Spitzer/GLIMPSE large-scale survey data provides infrared colours from 1.25 to 8$\mu$m for a large fraction of the inner Galactic plane. In 2005 we initiated a pilot study of the combined set of infrared colours for two GLIMPSE fields and showed that WRs typically occupy a sparsely populated region of the colour space. We followed up 42 of our WR candidates spectroscopically in the near-infrared, and with limited additional observations of some of these candidates in the optical. Six new WRs, four late-type WN and two late-type WC stars, were discovered as a result. Of the remaining $\sim$86% of the sample, five appear to be O-type stars. 21 stars are likely of type Be, and 10 stars appear to be of late-type, or possibly young stellar objects, which have "contaminated" the infrared color space. The survey is generally unbiased towards clusters or field stars, and the new WRs found are in both the field and in and around the RCW 49 region (including cluster Westerlund 2). In this work, and in our other recent work, we show that the infrared broad-band colours to be the most efficient means of identifying (particularly, dust-obscured) candidate massive stars, notably WRs.

Red supergiant stars (RSGs) and yellow hypergiant stars (YHGs) are believed to be the high-mass counterparts of stars in the AGB and early post-AGB phases. We study the mass-loss in the post main-sequence evolution of massive stars, through the properties of their envelopes in the intermediate and warm gas layers. These are the regions where the acceleration of the gas takes place and the most recent mass-loss episodes can be seen. We used the HIFI instrument on-board the Herschel Space Observatory to observe sub-mm and FIR transitions of CO, water, and their isotopologues in a sample of two RSGs (NML Cyg and Betelgeuse) and two YHGs (IRC+10420 and AFGL 2343) stars. We present an inventory of the detected lines and analyse the information revealed by their spectral profiles. On the basis of the results presented in an earlier study, we model the CO and 13CO emission in IRC+10420 and compare it to a set of lines ranging from the mm, to the FIR. Red supergiants have stronger high-excitation lines than the YHGs, indicating that they harbour dense and hot inner shells contributing to these transitions. Consequently, these high-J lines in RSGs originate from acceleration layers that have not yet reached the circumstellar terminal velocity and have narrower profiles than their flat-topped lower-J counterparts. The YHGs tend to lack this inner component, in line with the picture of detached, hollow envelopes derived from studies at longer wavelengths. NH3 is only detected in two sources (NML Cyg, IRC+10420), which are also observed to be the strongest water-line emitters of the studied sample. In contrast, OH is detected in all sources and does not seem to correlate with the water line intensities. We show that the IRC+10420 model derived solely from mm low-J CO transitions is capable of reproducing the high-J transitions when the temperature in the inner shell is simply lowered by about 30%.

In order to characterise the cloud structures responsible for the formation of high-mass stars, we present Herschel observations of the DR21 environment. Maps of the column density and dust temperature unveil the structure of the DR21 ridge and several connected filaments. The ridge has column densities larger than 1e23/cm^2 over a region of 2.3 pc^2. It shows substructured column density profiles and branching into two major filaments in the north. The masses in the studied filaments range between 130 and 1400 Msun whereas the mass in the ridge is 15000 Msun. The accretion of these filaments onto the DR21 ridge, suggested by a previous molecular line study, could provide a continuous mass inflow to the ridge. In contrast to the striations seen in e.g., the Taurus region, these filaments are gravitationally unstable and form cores and protostars. These cores formed in the filaments potentially fall into the ridge. Both inflow and collisions of cores could be important to drive the observed high-mass star formation. The evolutionary gradient of star formation running from DR21 in the south to the northern branching is traced by decreasing dust temperature. This evolution and the ridge structure can be explained by two main filamentary components of the ridge that merged first in the south.

We present Herschel images from the HOBYS key program of the Eagle Nebula (M16) in the far-infrared and sub-millimetre, using the PACS and SPIRE cameras at 70{\mu}m, 160{\mu}m, 250{\mu}m, 350{\mu}m, 500{\mu}m. M16, home to the Pillars of Creation, is largely under the influence of the nearby NGC6611 high-mass star cluster. The Herschel images reveal a clear dust temperature gradient running away from the centre of the cavity carved by the OB cluster. We investigate the heating effect of NGC6611 on the entire M16 star-forming complex seen by Herschel including the diffuse cloud environment and the dense filamentary structures identified in this region. In addition, we interpret the three-dimensional geometry of M16 with respect to the nebula, its surrounding environment, and the NGC6611 cavity. The dust temperature and column density maps reveal a prominent eastern filament running north-south and away from the high-mass star-forming central region and the NGC6611 cluster, as well as a northern filament which extends around and away from the cluster. The dust temperature in each of these filaments decreases with increasing distance from the NGC6611 cluster, indicating a heating penetration depth of \sim 10 pc in each direction in 3 - 6 \times 10^{22} cm-2 column density filaments. We show that in high-mass star-forming regions OB clusters impact the temperature of future star-forming sites, modifying the initial conditions for collapse and effecting the evolutionary criteria of protostars developed from spectral energy distributions. Possible scenarios for the origin of the morphology seen in this region are discussed, including a western equivalent to the eastern filament, which was destroyed by the creation of the OB cluster and its subsequent winds and radiation.

Molecular globules and pillars are spectacular features, found only in the interface region between a molecular cloud and an HII-region. Impacting Far-ultraviolet (FUV) radiation creates photon dominated regions (PDRs) on their surfaces that can be traced by typical cooling lines. With the GREAT receiver onboard SOFIA we mapped and spectrally resolved the [CII] 158 micron atomic fine-structure line and the highly excited 12CO J=11-10 molecular line from three objects in Cygnus X (a pillar, a globule, and a strong IRAS source). We focus here on the globule and compare our data with existing Spitzer data and recent Herschel Open-Time PACS data. Extended [CII] emission and more compact CO-emission was found in the globule. We ascribe this emission mainly to an internal PDR, created by a possibly embedded star-cluster with at least one early B-star. However, external PDR emission caused by the excitation by the Cyg OB2 association cannot be fully excluded. The velocity-resolved [CII] emission traces the emission of PDR surfaces, possible rotation of the globule, and high-velocity outflowing gas. The globule shows a velocity shift of ~2 km/s with respect to the expanding HII-region, which can be understood as the residual turbulence of the molecular cloud from which the globule arose. This scenario is compatible with recent numerical simulations that emphazise the effect of turbulence. It is remarkable that an isolated globule shows these strong dynamical features traced by the [CII]-line, but it demands more observational studies to verify if there is indeed an embedded cluster of B-stars.

Spectra, taken with the heterodyne instrument, HIFI, aboard the Herschel Space Observatory, of O-rich asymptotic giant branch (AGB) stars which form part of the guaranteed time key program HIFISTARS are presented. The aim of this program is to study the dynamical structure, mass-loss driving mechanism, and chemistry of the outflows from AGB stars as a function of chemical composition and initial mass. We used the HIFI instrument to observe nine AGB stars, mainly in the H2O and high rotational CO lines We investigate the correlation between line luminosity, line ratio and mass-loss rate, line width and excitation energy. A total of nine different molecules, along with some of their isotopologues have been identified, covering a wide range of excitation temperature. Maser emission is detected in both the ortho- and para-H2O molecules. The line luminosities of ground state lines of ortho- and para-H2O, the high-J CO and NH3 lines show a clear correlation with mass-loss rate. The line ratios of H2O and NH3 relative to CO J=6-5 correlate with the mass-loss rate while ratios of higher CO lines to the 6-5 is independent of it. In most cases, the expansion velocity derived from the observed line width of highly excited transitions formed relatively close to the stellar photosphere is lower than that of lower excitation transitions, formed farther out, pointing to an accelerated outflow. In some objects, the vibrationally excited H2O and SiO which probe the acceleration zone suggests the wind reaches its terminal velocity already in the innermost part of the envelope, i.e., the acceleration is rapid. Interestingly, for R Dor we find indications of a deceleration of the outflow in the region where the material has already escaped from the star.

We present the first Herschel PACS and SPIRE results of the Vela C molecular complex in the far-infrared and submillimetre regimes at 70, 160, 250, 350, and 500 um, spanning the peak of emission of cold prestellar or protostellar cores. Column density and multi-resolution analysis (MRA) differentiates the Vela C complex into five distinct sub-regions. Each sub-region displays differences in their column density and temperature probability distribution functions (PDFs), in particular, the PDFs of the `Centre-Ridge' and `South-Nest' sub-regions appear in stark contrast to each other. The Centre-Ridge displays a bimodal temperature PDF representative of hot gas surrounding the HII region RCW 36 and the cold neighbouring filaments, whilst the South-Nest is dominated by cold filamentary structure. The column density PDF of the Centre-Ridge is flatter than the South-Nest, with a high column density tail, consistent with formation through large-scale flows, and regulation by self-gravity. At small to intermediate scales MRA indicates the Centre-Ridge to be twice as concentrated as the South-Nest, whilst on larger scales, a greater portion of the gas in the South-Nest is dominated by turbulence than in the Centre-Ridge. In Vela C, high-mass stars appear to be preferentially forming in ridges, i.e., dominant high column density filaments.

We present Herschel-SPIRE Fourier Transform Spectrometer (FTS) and radio follow-up observations of two Herschel-ATLAS (H-ATLAS) detected strongly lensed distant galaxies. In one of the targeted galaxies H-ATLAS J090311.6+003906 (SDP.81) we detect [OIII] 88\mum and [CII] 158\mum lines at a signal-to-noise ratio of ~5. We do not have any positive line identification in the other fainter target H-ATLAS J091305.0-005343 (SDP.130). Currently SDP.81 is the faintest sub-mm galaxy with positive line detections with the FTS, with continuum flux just below 200 mJy in the 200-600 \mum wavelength range. The derived redshift of SDP.81 from the two detections is z=3.043 +/-0.012, in agreement with ground-based CO measurements. This is the first detection by Herschel of the [OIII] 88\mum line in a galaxy at redshift higher than 0.05. Comparing the observed lines and line ratios with a grid of photo-dissociation region (PDR) models with different physical conditions, we derive the PDR cloud density n ~ 2000 cm^{-3} and the far-UV ionizing radiation field G_0 ~ 200 (in units of the Habing field -- the local Galactic interstellar radiation field of 1.6x10^{-6} W/m^2). Using the CO derived molecular mass and the PDR properties we estimate the effective radius of the emitting region to be 500-700 pc. These characteristics are typical for star-forming, high redshift galaxies. The radio observations indicate that SDP.81 deviates significantly from the local FIR/radio correlation, which hints that some fraction of the radio emission is coming from an AGN. The constraints on the source size from millimiter-wave observations put a very conservative upper limit of the possible AGN contribution to less than 33%. These indications, together with the high [OIII]/FIR ratio and the upper limit of [OI] 63\mum/[CII] 158\mum suggest that some fraction of the ionizing radiation is likely to originate from an AGN.

We report the preliminary results of a survey for water vapor in a sample of eight C stars with large mid-IR continuum fluxes: V384 Per, CIT 6, V Hya, Y CVn, IRAS 15194-5115, V Cyg, S Cep, and IRC+40540. This survey, performed using the HIFI instrument on board the Herschel Space Observatory, entailed observations of the lowest transitions of both ortho- and para-water: the 556.936 GHz 1(10)-1(01) and 1113.343 GHz 1(11)-0(00) transitions, respectively. Water vapor was unequivocally detected in all eight of the target stars. Prior to this survey, IRC+10216 was the only carbon-rich AGB star from which thermal water emissions had been discovered, in that case with the use of the Submillimeter Wave Astronomy Satellite (SWAS). Our results indicate that IRC+10216 is not unusual, except insofar as its proximity to Earth leads to a large line flux that was detectable with SWAS. The water spectral line widths are typically similar to those of CO rotational lines, arguing against the vaporization of a Kuiper belt analog (Ford & Neufeld 2001) being the general explanation for water vapor in carbon-rich AGB stars. There is no apparent correlation between the ratio of the integrated water line fluxes to the 6.3 micron continuum flux - a ratio which measures the water outflow rate - and the total mass-loss rate for the stars in our sample.

In the framework of the HIFISTARS guaranteed time key programme, we measured more than 70 molecular emission lines with high signal-to-noise ratio towards VY CMa using the high-resolution HIFI spectrometer on board the Herschel satellite. The kinematic information obtained from the measured water lines supports the hypothesis of multiple outflow components. The observed high-intensity maser lines give no indication for strong polarisation.

CONTEXT: Water vapour maser emission from evolved oxygen-rich stars remains poorly understood. Additional observations, including polarisation studies and simultaneous observation of different maser transitions may ultimately lead to greater insight. AIMS: We have aimed to elucidate the nature and structure of the VY CMa water vapour masers in part by observationally testing a theoretical prediction of the relative strengths of the 620.701 GHz and the 22.235 GHz maser components of ortho water vapour. METHODS: In its high-resolution mode (HRS) the Herschel Heterodyne Instrument for the Infrared (HIFI) offers a frequency resolution of 0.125 MHz, corresponding to a line-of-sight velocity of 0.06 km/s, which we employed to obtain the strength and linear polarisation of maser spikes in the spectrum of VY CMa at 620.701 GHz. Simultaneous ground based observations of the 22.235 GHz maser with the Max-Planck-Institut f\"ur Radioastronomie 100-meter telescope at Effelsberg, provided a ratio of 620.701 GHz to 22.235 GHz emission. RESULTS:We report the first astronomical detection to date of water vapour maser emission at 620.701 GHz. In VY CMa both the 620.701 and the 22.235 GHz polarisation are weak. At 620.701 GHz the maser peaks are superposed on what appears to be a broad emission component, jointly ejected asymmetrically from the star. We observed the 620.701 GHz emission at two epochs 21 days apart, both to measure the potential direction of linearly polarised maser components and to obtain a measure of the longevity of these components. Although we do not detect significant polarisation levels in the core of the line, they rise up to approximately 6% in its wings.

Circumstellar envelopes (CSEs) of a variety of evolved stars have been found to contain ammonia (NH3) in amounts that exceed predictions from conventional chemical models by many orders of magnitude. The observations reported here were performed in order to better constrain the NH3 abundance in the CSEs of four, quite diverse, oxygen-rich stars using the NH3 ortho J_K = 1_0 - 0_0 ground-state line. We used the Heterodyne Instrument for the Far Infrared aboard Herschel to observe the NH3 J_K = 1_0 - 0_0 transition near 572.5 GHz, simultaneously with the ortho-H2O J_Ka,Kc = 1_1,0 -1_0,1 transition, toward VY CMa, OH 26.5+0.6, IRC+10420, and IK Tau. We conducted non-LTE radiative transfer modeling with the goal to derive the NH3 abundance in these objects' CSEs. For the latter two stars, Very Large Array imaging of NH3 radio-wavelength inversion lines were used to provide further constraints, particularly on the spatial extent of the NH3-emitting regions. Results. We find remarkably strong NH3 emission in all of our objects with the NH3 line intensities rivaling those obtained for the ground state H2O line. The NH3 abundances relative to H2 are very high and range from 2 x 10-7 to 3 x 10-6 for the objects we have studied. Our observations confirm and even deepen the circumstellar NH3 enigma. While our radiative transfer modeling does not yield satisfactory fits to the observed line profiles, it leads to abundance estimates that confirm the very high values found in earlier studies. New ways to tackle this mystery will include further Herschel observations of more NH3 lines and imaging with the Expanded Very Large Array.

A set of new, sensitive, and spectrally resolved, sub-millimeter line observations are used to probe the warm circumstellar gas around the S-type AGB star chi Cyg. The observed lines involve high rotational quantum numbers, which, combined with previously obtained lower-frequency data, make it possible to study in detail the chemical and physical properties of, essentially, the entire circumstellar envelope of chi Cyg. The data were obtained using the HIFI instrument aboard Herschel, whose high spectral resolution provides valuable information about the line profiles. Detailed, non-LTE, radiative transfer modelling, including dust radiative transfer coupled with a dynamical model, has been performed to derive the temperature, density, and velocity structure of the circumstellar envelope. We report the first detection of circumstellar H2O rotational emission lines in an S-star. Using the high-J CO lines to derive the parameters for the circumstellar envelope, we modelled both the ortho- and para-H2O lines. Our modelling results are consistent with the velocity structure expected for a dust-driven wind. The derived total H2O abundance (relative to H2) is (1.1 +/- 0.2)E-5, much lower than that in O-rich stars. The derived ortho-to-para ratio of 2.1 +/- 0.6 is close to the high-temperature equilibrium limit, consistent with H2O being formed in the photosphere.

We performed Herschel/HIFI observations of several CO lines in the far-infrared/sub-mm in the protoplanetary nebula CRL618. The high spectral resolution provided by HIFI allows measurement of the line profiles. Since the dynamics and structure of the nebula is well known from mm-wave interferometric maps, it is possible to identify the contributions of the different nebular components (fast bipolar outflows, double shells, compact slow shell) to the line profiles. The observation of these relatively high-energy transitions allows an accurate study of the excitation conditions in these components, particularly in the warm ones, which cannot be properly studied from the low-energy lines. The 12CO J=16-15, 10-9, and 6-5 lines are easily detected in this source. 13CO J=10-9 and 6-5 are also detected. Wide profiles showing spectacular line wings have been found, particularly in 12CO 16-15. Other lines observed simultaneously with CO are also shown. Our analysis of the CO high-J transitions, when compared with the existing models, confirms the very low expansion velocity of the central, dense component, which probably indicates that the shells ejected during the last AGB phases were driven by radiation pressure under a regime of maximum transfer of momentum. No contribution of the diffuse halo found from mm-wave data is identified in our spectra, because of its low temperature. We find that the fast bipolar outflow is quite hot, much hotter than previously estimated; for instance, gas flowing at 100 km/s must have a temperature higher than ~ 200 K. Probably, this very fast outflow, with a kinematic age < 100 yr, has been accelerated by a shock and has not yet cooled down. The double empty shell found from mm-wave mapping must also be relatively hot, in agreement with the previous estimate.

We report the first detection of the ground-state rotational transition of the methylidyne cation CH+ towards the massive star-forming region DR21 with the HIFI instrument onboard the Herschel satellite. The line profile exhibits a broad emission line, in addition to two deep and broad absorption features associated with the DR21 molecular ridge and foreground gas. These observations allow us to determine a CH+ J=1-0 line frequency of 835137 +/- 3 MHz, in good agreement with a recent experimental determination. We estimate the CH+ column density to be a few 1e13 cm^-2 in the gas seen in emission, and > 1e14 cm^-2 in the components responsible for the absorption, which is indicative of a high line of sight average abundance [CH+]/[H] > 1.2x10^-8. We show that the CH+ column densities agree well with the predictions of state-of-the-art C-shock models in dense UV-illuminated gas for the emission line, and with those of turbulent dissipation models in diffuse gas for the absorption lines.

We report the discovery of water vapour toward the carbon star V Cygni. We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 1(11) - 0(00) para-water transition at 1113.3430 GHz in the upper sideband of the Band 4b receiver. The observed spectral line profile is nearly parabolic, but with a slight asymmetry associated with blueshifted absorption, and the integrated antenna temperature is 1.69 \pm 0.17 K km/s. This detection of thermal water vapour emission, carried out as part of a small survey of water in carbon-rich stars, is only the second such detection toward a carbon-rich AGB star, the first having been obtained by the Submillimeter Wave Astronomy Satellite toward IRC+10216. For an assumed ortho-to-para ratio of 3 for water, the observed line intensity implies a water outflow rate ~ (3 - 6) E-5 Earth masses per year and a water abundance relative to H2 of ~ (2-5) E-6. This value is a factor of at least 1E+4 larger than the expected photospheric abundance in a carbon-rich environment, and - as in IRC+10216 - raises the intriguing possibility that the observed water is produced by the vapourisation of orbiting comets or dwarf planets. However, observations of the single line observed to date do not permit us to place strong constraints upon the spatial distribution or origin of the observed water, but future observations of additional transitions will allow us to determine the inner radius of the H2O-emitting zone, and the H2O ortho-to-para ratio, and thereby to place important constraints upon the origin of the observed water emission.

We present observations of the rotational ortho-water ground transition, the two lowest para-water transitions, and the ground transition of ionised ortho-water in the archetypal starburst galaxy M82, performed with the HIFI instrument on the Herschel Space Observatory. These observations are the first detections of the para-H2O(111-000) (1113\,GHz) and ortho-H2O+(111-000) (1115\,GHz) lines in an extragalactic source. All three water lines show different spectral line profiles, underlining the need for high spectral resolution in interpreting line formation processes. Using the line shape of the para-H2O(111-000) and ortho-H2O+(111-000) absorption profile in conjunction with high spatial resolution CO observations, we show that the (ionised) water absorption arises from a ~2000 pc^2 region within the HIFI beam located about ~50 pc east of the dynamical centre of the galaxy. This region does not coincide with any of the known line emission peaks that have been identified in other molecular tracers, with the exception of HCO. Our data suggest that water and ionised water within this region have high (up to 75%) area-covering factors of the underlying continuum. This indicates that water is not associated with small, dense cores within the ISM of M82 but arises from a more widespread diffuse gas component.

The Rosette molecular cloud is promoted as the archetype of a triggered star-formation site. This is mainly due to its morphology, because the central OB cluster NGC 2244 has blown a circular-shaped cavity into the cloud and the expanding HII-region now interacts with the cloud. Studying the spatial distribution of the different evolutionary states of all star-forming sites in Rosette and investigating possible gradients of the dust temperature will help to test the 'triggered star-formation' scenario in Rosette. We use continuum data obtained with the PACS (70 and 160 micron) and SPIRE instruments (250, 350, 500 micron) of the Herschel telescope during the Science Demonstration Phase of HOBYS. Three-color images of Rosette impressively show how the molecular gas is heated by the radiative impact of the NGC 2244 cluster. A clear negative temperature gradient and a positive density gradient (running from the HII-region/molecular cloud interface into the cloud) are detected. Studying the spatial distribution of the most massive dense cores (size scale 0.05 to 0.3 pc), we find an age-sequence (from more evolved to younger) with increasing distance to the cluster NGC 2244. No clear gradient is found for the clump (size-scale up to 1 pc) distribution. The existence of temperature and density gradients and the observed age-sequence imply that star formation in Rosette may indeed be influenced by the radiative impact of the central NGC 2244 cluster. A more complete overview of the prestellar and protostellar population in Rosette is required to obtain a firmer result.

We present a preliminary analysis of the small-scale structure found in new 70-520 micron continuum maps of the Rosette molecular cloud (RMC), obtained with the SPIRE and PACS instruments of the Herschel Space Observatory. We find 473 clumps within the RMC using a new structure identification algorithm, with sizes up to ~1.0 pc in diameter. A comparison with recent Spitzer maps reveals that 371 clumps are "starless" (without an associated young stellar object), while 102 are "protostellar." Using the respective values of dust temperature, we determine the clumps have masses (M_C) over the range -0.75 <= log (M_C/M_sun) <= 2.50. Linear fits to the high-mass tails of the resulting clump mass spectra (CMS) have slopes that are consistent with those found for high-mass clumps identified in CO emission by other groups.

The Herschel OB young stellar objects survey (HOBYS) has observed the Rosette molecular cloud, providing an unprecedented view of its star formation activity. These new far-infrared data reveal a population of compact young stellar objects whose physical properties we aim to characterise. We compiled a sample of protostars and their spectral energy distributions that covers the near-infrared to submillimetre wavelength range. These were used to constrain key properties in the protostellar evolution, bolometric luminosity, and envelope mass and to build an evolutionary diagram. Several clusters are distinguished including the cloud centre, the embedded clusters in the vicinity of luminous infrared sources, and the interaction region. The analysed protostellar population in Rosette ranges from 0.1 to about 15 Msun with luminosities between 1 and 150 Lsun, which extends the evolutionary diagram from low-mass protostars into the high-mass regime. Some sources lack counterparts at near- to mid-infrared wavelengths, indicating extreme youth. The central cluster and the Phelps & Lada 7 cluster appear less evolved than the remainder of the analysed protostellar population. For the central cluster, we find indications that about 25% of the protostars classified as Class I from near- to mid-infrared data are actually candidate Class 0 objects. As a showcase for protostellar evolution, we analysed four protostars of low- to intermediate-mass in a single dense core, and they represent different evolutionary stages from Class 0 to Class I. Their mid- to far-infrared spectral slopes flatten towards the Class I stage, and the 160 to 70um flux ratio is greatest for the presumed Class 0 source. This shows that the Herschel observations characterise the earliest stages of protostellar evolution in detail.

We present high-resolution infrared (2--18 micron) images of the archetypal periodic dust-making Wolf-Rayet binary system WR140 (HD 193793) taken between 2001 and 2005, and multi-colour (J -- [19.5]) photometry observed between 1989 and 2001. The images resolve the dust cloud formed by WR140 in 2001, allowing us to track its expansion and cooling, while the photometry allows tracking the average temperature and total mass of the dust. The combination of the two datasets constrains the optical properties of the dust. The most persistent dust features, two concentrations at the ends of a `bar' of emission to the south of the star, were observed to move with constant proper motions of 324+/-8 and 243+/-7 mas/y. Longer wavelength (4.68-micron and 12.5-micron) images shows dust emission from the corresponding features from the previous (1993) periastron passage and dust-formation episode. A third persistent dust concentration to the east of the binary (the `arm') was found to have a proper motion ~ 320 mas/y. Extrapolation of the motions of the concentrations back to the binary suggests that the eastern `arm' began expansion 4--5 months earlier than those in the southern `bar', consistent with the projected rotation of the binary axis and wind-collision region (WCR) on the sky. Comparison of model dust images and the observations constrain the intervals when the WCR was producing sufficiently compressed wind for dust nucleation in the WCR, and suggests that the distribution of this material was not uniform about the axis of the WCR, but more abundant in the following edge in the orbital plane.

We report the discovery of fifteen previously unknown Wolf-Rayet (WR) stars found as part of an infrared broad-band study of candidate WR stars in the Galaxy. We have derived an empirically-based selection algorithm which has selected ~5000 WR candidate stars located within the Galactic Plane drawn from the GLIMPSE (mid-infrared) and 2MASS (near-infrared) catalogues. Spectroscopic follow-up of 184 of these reveals eleven WN and four WC-type WR stars. Early WC subtypes are absent from our sample and none show evidence for circumstellar dust emission. Of the candidates which are not WR stars, ~120 displayed hydrogen emission line features in their spectra. Spectral features suggest that the majority of these are in fact B supergiants/hypergiants, ~40 of these are identified Be/B[e] candidates. Here, we present the optical spectra for six of the newly-detected WR stars, and the near-infrared spectra for the remaining nine of our sample. With a WR yield rate of ~7% and a massive star detection rate of ~65%, initial results suggest that this method is one of the most successful means for locating evolved, massive stars in the Galaxy.

Spitzer IRAC 3.6-8 micron photometry obtained as part of the GLIMPSE survey has revealed mid-infrared excesses for 33 field stars with known spectral types in a 1.2 sq. degree field centered on the southern Galactic HII region RCW49. These stars comprise a subset of 184 stars with known spectral classification, most of which were pre-selected to have unusually red IR colors. We propose that the mid-IR excesses are caused by circumstellar dust disks that are either very late remnants of stellar formation or debris disks generated by planet formation. Of these 33 stars, 29 appear to be main-sequence stars based on optical spectral classifications. Five of the 29 main-sequence stars are O or B stars with excesses that can be plausibly explained by thermal bremsstrahlung emission, and four are post main-sequence stars. The lone O star is an O4V((f)) at a spectrophotometric distance of 3233+ 540- 535 pc and may be the earliest member of the Westerlund 2 cluster. Of the remaining 24 main-sequence stars, 18 have SEDs that are consistent with hot dusty debris disks, a possible signature of planet formation. Modeling the excesses as blackbodies demonstrates that the blackbody components have fractional bolometric disk-to-star luminosity ratios, LIR/L*, ranging from 10^-3 to 10^-2 with temperatures ranging from 220 to 820 K. The inferred temperatures are more consistent with asteroid belts rather than the cooler temperatures expected for Kuiper belts. Mid-IR excesses are found in all spectral types from late B to early K.

Observations of the 1.083-micron He I line in WR140 (HD 193793) show P-Cygni profiles which varied as the binary system went through periastron passage. A sub-peak appeared on the normally flat-topped emission component and then moved across the profile consistent with its formation in the wind-collision region. Variation of the absorption component provided constraints on the opening angle (theta) of the wind-collision region. Infrared (2-10 micron) images observed with a variety of instruments in 2001-04 resolve the dust cloud formed in 2001, and show it to be expanding at a constant rate. Owing to the high eccentricity of the binary orbit, the dust is spread around the orbital plane in a 'splash' and we compare the dust images with the orientation of the orbit.

We determine and tabulate A(lambda)/A(K), the wavelength dependence of interstellar extinction, in the Galactic plane for 1.25um<lambda<8.0um along two lines of sight: l~42 and l~284. The first is a relatively quiescent and unremarkable region; the second contains the giant HII region RCW 49 as well as a ``field'' region unrelated to the cluster and nebulosity. Areas near these Galactic longitudes were imaged at J, H, and K bands by 2MASS and at 3--8um by Spitzer for the GLIMPSE Legacy program. We measure the mean values of the color excess ratios (A(lambda)-A(K))/(A(J)-A(K)) directly from the color distributions of observed stars. The extinction ratio between two of the filters, e.g. A(J)/A(K), is required to calculate A(lambda)/A(K) from those measured ratios. We use the apparent JHK magnitudes of giant stars along our two sightlines, and fit the reddening as a function of magnitude (distance) to determine A(J)/kpc, A(K)/kpc, and A(J)/A(K). Our values of A(lambda)/A(K) show a flattening across the 3--8um wavelength range, roughly consistent with the Lutz 1996 extinction measurements derived for the sightline toward the Galactic center.

Spitzer Space Telescope imaging as part of the Galactic Legacy Mid-Plane Survey Extraordinaire (GLIMPSE) reveals a previously unidentified low-latitude rich star cluster near l=31.3 degrees, b=-0.1 degrees. Near-infrared JHK' photometry from the Wyoming Infrared Observatory indicates an extinction of A_V ~ 15+/-3 mag for cluster members. Analysis of 13CO features along the same sightline suggests a probable kinematic distance of 3.1 -- 5.2 kpc. The new cluster has an angular diameter of ~1-2 pc, a total magnitude m_{K_0}=2.1 corrected for extinction, and a luminosity of M_K ~ -10.3 at 3.1 kpc. In contrast to young massive Galactic clusters with ages <100 Myr, the new cluster has no significant radio emission. Comparison to theoretical K-band luminosity functions indicates an age of at least several Gyr and a mass of at least 10^5 solar masses. Unlike known old open clusters, this new cluster lies in the inner Galaxy at R_{GC} ~ 6.1 kpc. We designate this object ``GLIMPSE-C01'' and classify it as a Milky Way globular cluster passing through the Galactic disk. We also identify a region of star formation and fan-shaped outflows from young stellar objects in the same field as the cluster. The cluster's passage through the Galactic molecular layer may have triggered this star formation activity.

The optically-dark globule IC 1396A is revealed using Spitzer images at 3.6, 4.5, 5.8, 8, and 24 microns to be infrared-bright and to contain a set of previously unknown protostars. The mid-infrared colors of the 24 microns detected sources indicate several very young (Class I or 0) protostars and a dozen Class II stars. Three of the new sources (IC 1396A: gamma, delta, and epsilon) emit over 90% of their bolometric luminosities at wavelengths greater than 3 microns, and they are located within ~0.02 pc of the ionization front at the edge of the globule. Many of the sources have spectra that are still rising at 24 microns. The two previously-known young stars LkHa 349 a and c are both detected, with component c harboring a massive disk and component a being bare. Of order 5% of the mass of material in the globule is presently in the form of protostars in the 10^5 to 10^6 yr age range. This high star formation rate was likely triggered by radiation from a nearby O star.

GLIMPSE imaging using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope indicates that star formation is ongoing in the RCW 49 giant H II region. A photometric comparison of the sources in RCW 49 to a similar area to its north finds that at least 300 stars brighter than 13th magnitude in band [3.6] have infrared excesses inconsistent with reddening due to foreground extinction. These are likely young stellar objects (YSOs) more massive than 2.5 M_sun suggesting that thousands more low-mass stars are forming in this cloud. Some of the YSOs are massive (B stars) and therefore very young, suggesting that a new generation of star formation is occurring, possibly triggered by stellar winds and shocks generated by the older (2-3 Myr) central massive cluster. The Spitzer IRAC camera has proven to be ideally suited for distinguishing young stars from field stars, and the GLIMPSE survey of the Galactic Plane will likely find thousands of new star formation regions.

GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordinaire), a SIRTF Legacy Science Program, will be a fully sampled, confusion-limited infrared survey of the inner two-thirds of the Galactic disk with a pixel resolution of \~1.2" using the Infrared Array Camera (IRAC) at 3.6, 4.5, 5.8, and 8.0 microns. The survey will cover Galactic latitudes |b| <1 degree and longitudes |l|=10 to 65 degrees (both sides of the Galactic center). The survey area contains the outer ends of the Galactic bar, the Galactic molecular ring, and the inner spiral arms. The GLIMPSE team will process these data to produce a point source catalog, a point source data archive, and a set of mosaicked images. We summarize our observing strategy, give details of our data products, and summarize some of the principal science questions that will be addressed using GLIMPSE data. Up-to-date documentation, survey progress, and information on complementary datasets are available on the GLIMPSE web site: www.astro.wisc.edu/glimpse.

We present optical spectroscopic observations of 93 galaxies taken from the infra-red selected COLA (Compact Objects in Low Power AGN) sample. The sample spans the range of far-IR luminosities from normal galaxies to LIRGs. Of the galaxies observed, 78 (84%) exhibit emission lines. Using a theoretically-based optical emission-line scheme we classify 15% of the emission-line galaxies as Seyferts, 77% as starbursts, and the rest are either borderline AGN/starburst or show ambiguous characteristics. We find little evidence for an increase in the fraction of AGN in the sample as a function of far-IR luminosity but our sample covers only a small range in infrared luminosity and thus a weak trend may be masked. As a whole the Seyfert galaxies exhibit a small, but significant, radio excess on the radio-FIR correlation compared to the galaxies classified as starbursts. Compact (<0.05'') radio cores are detected in 55% of the Seyfert galaxies, and these galaxies exhibit a significantly larger radio excess than the Seyfert galaxies in which cores were not detected. Our results indicate that there may be two distinct populations of Seyferts, ``radio-excess'' Seyferts, which exhibit extended radio structures and compact radio cores, and ``radio-quiet'' Seyferts, in which the majority of the radio emission can be attributed to star-formation in the host galaxy. No significant difference is seen between the IR and optical spectroscopic properties of Seyferts with and without radio cores. (Abridged)

We present the first detection of 12 CO J=2->1 and 12 CO J=1->0 emission from the LBV AG Carinae. AG Carinae resides in a region which is very rich in molecular gas with complex motions. We find evidence of a slow outflow of molecular gas, expanding at ~ 7 km/s. This emission appears spatially unresolved. We argue that it is spatially localised, rather than extended, and possibly associated with the immediate circumstellar region of AG Carinae. Does it originate from a circumstellar envelope, similar to carbon stars, or from a circumstellar disk? The option of the circumstellar disk is preferable because it is consistent with additional independent indications for the existence of wind asymmetries in close proximity to the central star, found from spectropolarimetry and analysis of the UV and optical line profiles, and it provides the conditions of density and shielding necessary for the survival of the CO molecules in proximity to such a hot star (Teff ~ 14000 K - 20000 K). In the assumption that the CO emission originated when AG Carinae was in an evolved state, we derive a lower limit to the mass of molecular gas of 2.8 solar masses. This is smaller, but still comparable with the mass of ionized gas present in the circumstellar environment (4.2 solar masses), with the implication that the molecular gas fraction can contribute significantly to the overall mass lost from the central star in its post main sequence evolution.

We present the first results from the COLA (Compact Objects in Low-power AGN) project which aims to determine the relationship between one facet of AGN activity, the compact radio core, with star formation in the circumnuclear region of the host galaxy. This will be accomplished by the comparison of the multi-wavelength properties of a sample of AGN with compact radio cores to those of a sample of AGN without compact cores and a matched sample of galaxies without AGN. In this paper we discuss the selection criteria for our galaxy samples and present the initial radio observations of the 107 Southern galaxies in our sample. Low-resolution ATCA observations at 4.8, 2.5 and 1.4 GHz and high resolution, single baseline snapshots at 2.3 GHz with the Australian LBA are presented. We find that for the majority of the galaxies in our sample, the radio luminosity is correlated with the FIR luminosity. Compact radio cores are detected in 9 galaxies. The majority (8/9) of these galaxies exhibit a significant radio excess and 50% (7/14) of the galaxies which lie above the radio-FIR correlation by more than 1 sigma have compact radio cores. The emission from the cores is too weak to account for this radio excess and there is no evidence that the radio luminosity of the compact cores is correlated with the FIR galaxy luminosity. The galaxies with compact cores tend to be classified optically as AGN, with two thirds (6/9) exhibiting Seyfert-like optical emission line ratios, and the remaining galaxies classified either as composite objects (2/9) or starburst (1/9). (Abridged)

This is the first of two papers in which a study is made of a sample of 12 southern radio-luminous Seyfert galaxies. Our aim is to investigate possible correlations between radio morphology and nuclear/circumnuclear emission-line properties. In this paper we present radio images at 13, 6, and 3 cm taken with the Australia Telescope Compact Array (ATCA), global far-infrared (FIR) properties for the whole sample, and optical and near-infrared (NIR) spectroscopy of an interesting subset. We find a mixture of radio morphologies, including linear, diffuse and compact sources. When the FIR colors of the galaxies are considered there is an indication that the compact radio sources have warmer FIR colors than the diffuse sources, whereas the linear sources span a wide range of FIR colors. There is a wide variation in radio spectral-indices, suggesting that free-free absorption is significant in some systems, particularly IRAS 11249-2859, NGC 4507, and NGC 7213. Detailed emission-line studies are presented of 4 galaxies IC 3639, NGC 5135, NGC 3393 & IRAS 11249-2859. In IC 3639 we present evidence of vigorous, compact star formation enclosed by very extended [OI]6300 emission, suggestive of the boundary between a diffuse outflow and the surrounding ISM. In another galaxy, IC 5063, we see evidence for the possible interaction of a highly collimated outflow and the surrounding rotating inner disk. Of the 5 galaxies which show compact radio emission, 4 have radio/FIR flux ratios consistent with an energetically dominant AGN, whereas IC 4995 exhibits evidence for a very compact starburst.

We present the results of a complete Einstein Imaging Proportional Counter (IPC) X-ray survey of optically-selected galaxies from Shapley-Ames Catalog, Uppsala General Catalog and the European Southern Observatory Catalog. X-ray fluxes for 1018 galaxies have been surveyed, 21% of which are detected. 827 galaxies have either independent distance estimates or radial velocities. Associated optical, redshift and distance data has been assembled. The five large data tables in the paper are in easily-accessible ascii formats with an ADC-style README file, for easy electronic data handling. The accuracy of the X-ray fluxes has been checked in three different ways; all are consistent with the derived X-ray fluxes being of $\le$0.1 dex accuracy. In particular, there is agreement with previous published X-ray fluxes for galaxies in common with Roberts et al. (1991) and Fabbiano et al. (1992). The data presented here will be used in further studies to characterize the X-ray output of galaxies of various morphological types and thus to enable the determination of the major sources contributing to the X-ray emission from galaxies.

We present images of the peculiar galaxy Arp~10 which reveal two rings of concurrent star formation. Apart from a bright ring visible on early photographs of the system, an even brighter inner ring of \ion{H}{2} regions is found within the nuclear bulge of the galaxy. A very faint ring-arc of \ion{H}{2} regions is also seen associated with a third outer ring or shell. An investigation of the H$\alpha$ fluxes in the rings reveals an increase in the emission where the ring surface density in the R-band light exceeds 21.3$\pm$0.2 mag arcsec$^{-2}$. If the R-band light is dominated by old stars in the underlying density wave, then the results suggest evidence for a Star Formation law which exhibits a threshold dependence on the strength of the density wave in the rings. Even if the R-band continuum in the ring is heavily contaminated with red light from the underlying young stars then a smaller, but still significant non-linear enhancement in the star formation rates (SFR) in one segment of the second ring is required to explain the results. In either case, strong trends in the SFR with azimuth around Ring~2 are in good agreement with off-center collisional ring galaxy models.