We present the projected Rayleigh statistic (PRS), a modification of the classic Rayleigh statistic, as a test for non-uniform relative orientation between two sets of pseudo-vector fields. This gives an effective way of investigating whether polarization is preferentially parallel or perpendicular to filaments in the interstellar medium, for example; there are other potential applications in astrophysics, e.g. when comparing small-scale orientations with larger-scale shear patterns. We compare the efficiency of the PRS against histogram binning methods that have previously been used for characterising the relative orientations of gas column density structures with the magnetic field projected on the plane of the sky. We specifically examine data for the Vela C molecular cloud, where the column density is inferred from Herschel submillimetre observations, and the magnetic field from observations by the Balloon-borne Large-Aperture Submillimetre Telescope in the 250-, 350-, and 500-{\mu}m wavelength bands. We find that the projected Rayleigh statistic has greater statistical power than approaches that bin the relative orientation angles, and hence this circular statistic makes more efficient use of the information contained in the orientation data.

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

We statistically evaluate the relative orientation between gas column density structures, inferred from Herschel submillimetre observations, and the magnetic field projected on the plane of sky, inferred from polarized thermal emission of Galactic dust observed by BLASTPol at 250, 350, and 500 micron, towards the Vela C molecular complex. First, we find very good agreement between the polarization orientations in the three wavelength-bands, suggesting that, at the considered common angular resolution of 3.0 arcminutes that corresponds to a physical scale of approximately 0.61 pc, the inferred magnetic field orientation is not significantly affected by temperature or dust grain alignment effects. Second, we find that the relative orientation between gas column density structures and the magnetic field changes progressively with increasing gas column density, from mostly parallel or having no preferred orientation at low column densities to mostly perpendicular at the highest column densities. This observation is in agreement with previous studies by the Planck collaboration towards more nearby molecular clouds. Finally, we find a correspondence between the trends in relative orientation and the shape of the column density probability distribution functions. In the sub-regions of Vela C dominated by one clear filamentary structure, or "ridges", we find a sharp transition from preferentially parallel or having no preferred relative orientation at low column densities to preferentially perpendicular at highest column densities. In the sub-regions of Vela C dominated by several filamentary structures with multiple orientations, or "nests", such a transition is also present, but it is clearly less sharp than in the ridge-like sub-regions. Both of these results suggest that the magnetic field is dynamically important for the formation of density structures in this region.

Within four nearby (d < 160 pc) molecular clouds, we statistically evaluate the structure of the interstellar magnetic field, projected on the plane of the sky and integrated along the line of sight, as inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz and from the optical and NIR polarization of background starlight. We compare the dispersion of the field orientation directly in vicinities with an area equivalent to that subtended by the Planck effective beam at 353 GHz (10') and using the second-order structure functions of the field orientation angles. We find that the average dispersion of the starlight-inferred field orientations within 10'-diameter vicinities is less than 20 deg, and that at these scales the mean field orientation is on average within 5 deg of that inferred from the submillimetre polarization observations in the considered regions. We also find that the dispersion of starlight polarization orientations and the polarization fractions within these vicinities are well reproduced by a Gaussian model of the turbulent structure of the magnetic field, in agreement with the findings reported by the Planck collaboration at scales greater than 10' and for comparable column densities. At scales greater than 10', we find differences of up to 14.7 deg between the second-order structure functions obtained from starlight and submillimetre polarization observations in the same positions in the plane of the sky, but comparison with a Gaussian model of the turbulent structure of the magnetic field indicates that these differences are small and are consistent with the difference in angular resolution between both techniques.

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

We report the detection of Galactic supernova remnants (SNRs) in the mid-infrared (at 24 and 70 {\mu}m), in the coordinate ranges 10 < l < 65 deg and 285 < l < 350 deg, |b| < 1 deg, using the Multiband Imaging Photometer (MIPS) aboard the Spitzer Space Telescope. We search for infrared counterparts to SNRs in Green's catalog and identify 39 out of 121, i.e., a detection rate of about 32%. Such a relatively low detection fraction is mainly due to confusion with nearby foreground/background sources and diffuse emission. The SNRs in our sample show a linear trend in [F8/F24] versus [F70/F24]. We compare their infrared fluxes with their corresponding radio flux at 1.4 GHz and find that most remnants have ratios of 70 {\mu}m to 1.4 GHz characteristic of SNRs (with the exception of a few which have ratios closer to those of H II regions). Furthermore, we retrieve a slope close to unity when correlating infrared (24 and 70 {\mu}m) with 1.4 GHz emission. Our survey is more successful in detecting remnants with bright X-ray emission, which we find is well correlated with the 24 {\mu}m morphology. Moreover, by comparing the power emitted in X-ray, infrared and radio, we conclude that the energy released in the infrared is comparable to the cooling in the X-ray range.

Observations of HDO are an important complement for studies of water, because they give strong constraints on the formation processes -- grain surfaces versus energetic process in the gas phase, e.g. in shocks. The HIFI observations of multiple transitions of HDO in Sgr~B2(M) presented here allow the determination of the HDO abundance throughout the envelope, which has not been possible before with ground-based observations only. The abundance structure has been modeled with the spherical Monte Carlo radiative transfer code RATRAN, which also takes radiative pumping by continuum emission from dust into account. The modeling reveals that the abundance of HDO rises steeply with temperature from a low abundance ($2.5\times 10^{-11}$) in the outer envelope at temperatures below 100~K through a medium abundance ($1.5\times 10^{-9}$) in the inner envelope/outer core, at temperatures between 100 and 200~K, and finally a high abundance ($3.5\times 10^{-9}$) at temperatures above 200~K in the hot core.

To investigate the accretion and feedback processes in massive star formation, we analyze the shapes of emission lines from hot molecular cores, whose asymmetries trace infall and expansion motions. The high-mass star forming region SgrB2(M) was observed with Herschel/HIFI (HEXOS key project) in various lines of HCN and its isotopologues, complemented by APEX data. The observations are compared to spherically symmetric, centrally heated models with density power-law gradient and different velocity fields (infall or infall+expansion), using the radiative transfer code RATRAN. The HCN line profiles are asymmetric, with the emission peak shifting from blue to red with increasing J and decreasing line opacity (HCN to H$^{13}$CN). This is most evident in the HCN 12--11 line at 1062 GHz. These line shapes are reproduced by a model whose velocity field changes from infall in the outer part to expansion in the inner part. The qualitative reproduction of the HCN lines suggests that infall dominates in the colder, outer regions, but expansion dominates in the warmer, inner regions. We are thus witnessing the onset of feedback in massive star formation, starting to reverse the infall and finally disrupting the whole molecular cloud. To obtain our result, the THz lines uniquely covered by HIFI were critically important.

We have observed CH absorption lines ($J=3/2, N=1 \leftarrow J=1/2, N=1$) against the continuum source Sgr~B2(M) using the \textit{Herschel}/HIFI instrument. With the high spectral resolution and wide velocity coverage provided by HIFI, 31 CH absorption features with different radial velocities and line widths are detected and identified. The narrower line width and lower column density clouds show `spiral arm' cloud characteristics, while the absorption component with the broadest line width and highest column density corresponds to the gas from the Sgr~B2 envelope. The observations show that each `spiral arm' harbors multiple velocity components, indicating that the clouds are not uniform and that they have internal structure. This line-of-sight through almost the entire Galaxy offers unique possibilities to study the basic chemistry of simple molecules in diffuse clouds, as a variety of different cloud classes are sampled simultaneously. We find that the linear relationship between CH and H$_2$ column densities found at lower $A_V$ by UV observations does not continue into the range of higher visual extinction. There, the curve flattens, which probably means that CH is depleted in the denser cores of these clouds.

Mon R2, at a distance of 830 pc, is the only ultracompact HII region (UC HII) where the photon-dominated region (PDR) between the ionized gas and the molecular cloud can be resolved with Herschel. HIFI observations of the abundant compounds 13CO, C18O, o-H2-18O, HCO+, CS, CH, and NH have been used to derive the physical and chemical conditions in the PDR, in particular the water abundance. The 13CO, C18O, o-H2-18O, HCO+ and CS observations are well described assuming that the emission is coming from a dense (n=5E6 cm-3, N(H2)>1E22 cm-2) layer of molecular gas around the UC HII. Based on our o-H2-18O observations, we estimate an o-H2O abundance of ~2E-8. This is the average ortho-water abundance in the PDR. Additional H2-18O and/or water lines are required to derive the water abundance profile. A lower density envelope (n~1E5 cm-3, N(H2)=2-5E22 cm-2) is responsible for the absorption in the NH 1_1-0_2 line. The emission of the CH ground state triplet is coming from both regions with a complex and self-absorbed profile in the main component. The radiative transfer modeling shows that the 13CO and HCO+ line profiles are consistent with an expansion of the molecular gas with a velocity law, v_e =0.5 x (r/Rout)^{-1} km/s, although the expansion velocity is poorly constrained by the observations presented here.

In Herschel images of the Galactic plane and many star forming regions, a major factor limiting our ability to extract faint compact sources is cirrus confusion noise, operationally defined as the "statistical error to be expected in photometric measurements due to confusion in a background of fluctuating surface brightness". The histogram of the flux densities of extracted sources shows a distinctive faint-end cutoff below which the catalog suffers from incompleteness and the flux densities become unreliable. This empirical cutoff should be closely related to the estimated cirrus noise and we show that this is the case. We compute the cirrus noise directly, both on Herschel images from which the bright sources have been removed and on simulated images of cirrus with statistically similar fluctuations. We connect these direct estimates with those from power spectrum analysis, which has been used extensively to predict the cirrus noise and provides insight into how it depends on various statistical properties and photometric operational parameters. We report multi-wavelength power spectra of diffuse Galactic dust emission from Hi-GAL observations at 70 to 500 microns within Galactic plane fields at l= 30 degrees and l= 59 degrees. We find that the exponent of the power spectrum is about -3. At 250 microns, the amplitude of the power spectrum increases roughly as the square of the median brightness of the map and so the expected cirrus noise scales linearly with the median brightness. Generally, the confusion noise will be a worse problem at longer wavelengths, because of the combination of lower angular resolution and the rising power spectrum of cirrus toward lower spatial frequencies, but the photometric signal to noise will also depend on the relative spectral energy distribution of the source compared to the cirrus.

We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158\mu\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2.

Ever since the pioneering study of Spitzer, it has been widely recognized that grains play an important role in the heating and cooling of photo-ionized environments. A detailed code is necessary to model grains in a photo-ionized medium. In this paper we will use the spectral synthesis code Cloudy for this purpose. The grain model of this code has recently undergone a comprehensive upgrade. Part of this is the newly developed hybrid grain charge model, which will be described in detail. This model allows discrete charge states of very small grains to be modelled accurately while simultaneously avoiding the overhead of fully resolving the charge distribution of large grains, thus making the model both accurate and computationally efficient. A comprehensive comparison with the fully resolved charge state models of Weingartner & Draine (2001) shows that the agreement is very satisfactory for realistic size distributions. The effect of the grain size distribution on the line emission from photo-ionized regions is studied by taking standard models for an H II region and a planetary nebula and adding a dust component to the models with varying grain size distributions. A comparison of the models shows that varying the size distribution has a dramatic effect on the emitted spectrum, and affects the ionization balance as well. These results clearly demonstrate that the grain size distribution is an important parameter in photo-ionization models.

In this paper we present results obtained with the new grain code in Cloudy which underline the strong effect of photo-electric heating by grains in photo-ionized regions. We study the effect that the distribution of grain sizes has on the magnitude of the effect, and show that this effect is nothing short of dramatic. This makes the grain size distribution an important parameter in modeling of photo-ionized regions such as H II regions and planetary nebulae.

We present HST/STIS long-slit spectroscopy of NGC 1976. Our goal is to measure the intrinsic line ratio [O III] 4364/5008 and thereby evaluate the electron temperature (T_e) and the fractional mean-square T_e variation (t_A^2) across the nebula. We also measure the intrinsic line ratio [N II] 5756/6585 in order to estimate T_e and t_A^2 in the N^+ region. The interpretation of the [N II] data is not as clear cut as the [O III] data because of a higher sensitivity to knowledge of the electron density as well as a possible contribution to the [N II] 5756 emission by recombination (and cascading). We present results from binning the data along the various slits into tiles that are 0.5" square (matching the slit width). The average [O III] temperature for our four HST/STIS slits varies from 7678 K to 8358 K; t_A^2 varies from 0.00682 to at most 0.0176. For our preferred solution, the average [N II] temperature for each of the four slits varies from 9133 K to 10232 K; t_A^2 varies from 0.00584 to 0.0175. The measurements of T_e reported here are an average along each line of sight. Therefore, despite finding remarkably low t_A^2, we cannot rule out significantly larger temperature fluctuations along the line of sight. The result that the average [N II] T_e exceeds the average [O III] T_e confirms what has been previously found for Orion and what is expected on theoretical grounds. Observations of the proplyd P159-350 indicate: large local extinction associated; ionization stratification consistent with external ionization by theta^1 Ori C; and indirectly, evidence of high electron density.

We present new HST/WFPC2 imagery and STIS long-slit spectroscopy of the planetary nebula NGC 7009. The primary goal was to obtain high spatial resolution of the intrinsic line ratio [O III] 4364/5008 and thereby evaluate the electron temperature (T_e) and the fractional mean-square T_e variation (t_A^2) across the nebula. The WFPC2 T_e map is rather uniform; almost all values are between 9000 - 11,000 K, with the higher T_e's closely coinciding with the inner He^++-zone. The results indicate very small values - <~ 0.01 - for t_A^2 throughout. Our STIS data allow an even more direct determination of T_e and t_A^2, albeit for a much smaller area than with WFPC2. We present results from binning the data along the slit into tiles that are 0.5'' square (matching the slit width). The average [O III] temperature using 45 tiles (excluding the central star and STIS fiducial bars) is 10,139 K; t_A^2 is 0.0035. The measurements of T_e reported here are an average along each line of sight. Therefore, despite finding remarkably low t_A^2, we cannot completely rule out temperature fluctuations along the line of sight as the cause of the large abundance discrepancy between heavy element abundances inferred from collisionally excited emission lines compared to those derived from recombination lines.

Ever since the pioneering study of Spitzer, it has been widely recognized that grains play an important role in the heating and cooling of photo-ionized environments. This includes the diffuse ISM, as well as H II regions, planetary nebulae, and photo-dissociation regions. A detailed code is necessary to model grains in a photo-ionized medium since the interactions of grains with their environment include a host of microphysical processes, and their importance can only be judged by performing a complete simulation. In this paper we will use the spectral synthesis code Cloudy for this purpose. A comprehensive upgrade of the grain model has been recently incorporated in Cloudy, and certain aspects of this upgrade will be discussed. Special emphasis will be on the new grain charge model. We will consider in detail the physics of grains in both ionized and neutral environments, and will present a calculation of photo-electric heating and collisional cooling rates for a range of physical conditions and grain materials and for a range of grain sizes (including a realistic size distribution). We conclude with a brief discussion of the problems currently hampering progress in this field. The new grain model will be used to model the silicate emission in the Ney-Allen nebula, and will help us better understand the nature of the grains in that part of the Orion complex.

This paper describes the most important differences between Cloudy 90 and Cloudy 94, the most recent upgrade of the code which has now been posted on the WWW. It also describes the upgrade of the grain model in Cloudy which is currently underway.

In this paper we study new infrared spectra of the evolved planetary nebula NGC 6445 obtained with ISO. These data show that the thermal emission from the grains is very cool and has a low flux compared to H beta. A model of the ionized region is constructed, using the photo-ionization code CLOUDY 90.05. Based on this model, we show from depletions in the gas phase elements that little grain destruction can have occurred in the ionized region of NGC 6445. We also argue that dust-gas separation in the nebula is not plausible. The most likely conclusion is that grains are residing inside the ionized region of NGC 6445 and that the low temperature and flux of the grain emission are caused by the low luminosity of the central star and the low optical depth of the grains. This implies that the bulk of the silicon-bearing grains in this nebula were able to survive exposure to hard UV photons for at least several thousands of years, contradicting previously published results. A comparison between optical and infrared diagnostic line ratios gives a marginal indication for the presence of a t^2-effect in the nebula. However, the evidence is not convincing and other explanations for the discrepancy are also plausible. The off-source spectrum taken with ISO-LWS clearly shows the presence of a warm cirrus component with a temperature of 24 K as well as a very cold component with a temperature of 7 K. Since our observation encompasses only a small region of the sky, it is not clear how extended the 7 K component is and whether it contributed significantly to the FIRAS spectrum taken by COBE. Because our line of sight is in the galactic plane, the very cold component could be a starless core.

Acquiring data for spectral classification of heavily reddened stars using traditional criteria in the blue-violet region of the spectrum can be prohibitively time consuming using small to medium sized telescopes. One such star is the Vatican Observatory emission-line star VES 735, which we have found excites the HII region KR140. In order to classify VES 735, we have constructed an atlas of stellar spectra of O stars in the yellow-green (4800 - 5420 A). We calibrate spectral type versus the line ratio HeI 4922:HeII 5411, showing that this ratio should be useful for the classification of heavily reddened O stars associated with HII regions. Application to VES 735 shows that the spectral type is O8.5. The absolute magnitude suggests luminosity class V. Comparison of the rateof emission of ionizing photons and the bolometric luminosity of VES 735, inferred from radio and infrared measurements of the KR 140 region, to recent stellar models gives consistent evidence for a main sequence star of mass 25 M_sun and age less than a few million years with a covering factor 0.4 - 0.5 by the nebular material. Spectra taken in the red (6500 - 6700 A) show that the stellar H-alpha emission is double-peaked about the systemic velocity and slightly variable. H-beta is in absorption, so that the emission-line classification is ``(e)''. However, unlike the case of the more well known O(e) star Zeta Oph, the emission from VES 735 appears to be long-lived rather than episodic.