We have used the ratio of column densities derived independently from the 850-um continuum JCMT Plane Survey (JPS) and the $^{13}$CO/C$^{18}$O J=3-2 Heterodyne Inner Milky Way Plane Survey (CHIMPS) to produce maps of the dense-gas mass fraction (DGMF) in two slices of the Galactic Plane centred at l=30 and l=40. The observed DGMF is a metric for the instantaneous clump-formation efficiency (CFE) in the molecular gas and a two-dimensional power-spectrum analysis of the DGMF reveals a break in slope at the approximate size scale of molecular clouds. We interpret this as the characteristic scale of the amplitude of variations in the CFE and a constraint on the dominant mechanism regulating the CFE and, hence, the star-formation efficiency in CO-traced clouds. By splitting the two fields into velocity components corresponding to the spiral arms that cross them, we find that individual molecular clouds have mean CFE values of 40, 41, and 46 per cent for the Scutum--Centaurus, Sagittarius, and Perseus spiral arms, respectively

The role of ionization feedback on high-mass (>8 Msun) star formation (HMSF) is still highly debated. Questions remain concerning the presence of nearby HII regions changes the properties of early HMSF and whether HII regions promote or inhibit the formation of high-mass stars. To characterize the role of HII regions on the HMSF, we study the properties of a sample of candidates high-mass starless clumps (HMSCs), of which about 90% have masses larger than 100 Msun. These high-mass objects probably represent the earliest stages of HMSF; we search if (and how) their properties are modified by the presence of an HII region. We took advantage of the recently published catalog of HMSC candidates. By cross matching the HMSCs and HII regions, we classified HMSCs into three categories: 1) The HMSCs associated with HII regions both in the position in the projected plane of the sky and in velocity; 2) HMSCs associated in the plane of the sky, but not in velocity; and 3) HMSCs far away from any HII regions in the projected sky plane. We carried out comparisons between associated and nonassociated HMSCs based on statistical analyses of multiwavelength data from infrared to radio. Statistical analyses suggest that HMSCs associated with HII regions are warmer, more luminous, more centrally-peaked and turbulent. We also clearly show, for the first time, that the ratio of bolometric luminosity to envelope mass of HMSCs (L/M) could not be a reliable evolutionary probe for early HMSF due to the external heating effects of the HII regions. More centrally peaked and turbulent properties of HMSCs associated with HII regions may promote the formation of high-mass stars by limiting fragmentation. High resolution interferometric surveys toward HMSCs are crucial to reveal how HII regions impact the star formation process inside HMSCs.

We search for far-infrared (FIR) counterparts of known supernova remnants (SNRs) in the Galactic plane (360 degrees in longitude and b = +/- 1 deg ) at 70 - 500 micron with Herschel. We detect dust signatures in 39 SNRs out of 190, made up of 13 core-collapse supernovae (CCSNe), including 4 Pulsar Wind Nebulae (PWNe), and 2 Type Ia SNe. A further 24 FIR detected SNRs have unknown types. We confirm the FIR detection of ejecta dust within G350.1-0.3, adding to the known sample of ~10 SNRs containing ejecta dust. We discover dust features at the location of a radio core at the centre of G351.2+0.1, indicating FIR emission coincident with a possible Crab-like compact object, with dust temperature and mass of Td = 45.8 K and Md = 0.18 solar mass, similar to the PWN G54.1+0.3. We show that the detection rate is higher among young SNRs. We produce dust temperature maps of 11 SNRs and mass maps of those with distance estimates, finding dust at temperatures 15 < Td < 40 K. If the dust is heated by shock interactions the shocked gas must be relatively cool and/or have a low density to explain the observed low grain temperatures.

We have analysed Herschel observations of M31, using the PPMAP procedure. The resolution of PPMAP images is sufficient (31 pc on M31) that we can analyse far-IR dust emission on the scale of Giant Molecular Clouds. By comparing PPMAP estimates of the far-IR emission optical depth at 300 microns (tau_300), and the near-IR extinction optical depth at 1.1 microns (tau_1.1) obtained from the reddening of RGB stars, we show that the ratio R_OBS.tau = tau_1.1/tau_300 falls in the range 500 to 1500. Such low values are incompatible with many commonly used theoretical dust models, which predict values of R_MODEL.kappa = kappa_1.1/kappa_300 (where kappa is the dust opacity coefficient) in the range 2500 to 4000. That is, unless a large fraction, at least 60%, of the dust emitting at 300 microns is in such compact sources that they are unlikely to intercept the lines of sight to a distributed population like RGB stars. This is not a new result: variants obtained using different observations and/or different wavelengths have already been reported by other studies. We present two analytic arguments for why it is unlikely that at least 60% of the emitting dust is in sufficiently compact sources. Therefore it may be necessary to explore the possibility that the discrepancy between observed values of R_OBS.tau and theoretical values of R_MODEL.kappa is due to limitations in existing dust models. PPMAP also allows us to derive optical-depth weighted mean values for the emissivity index, beta = - dln(kappa_lambda)/dln(lambda), and the dust temperature, T, denoted betabar and Tbar. We show that, in M31, R_OBS.tau is anti-correlated with betabar according to R_OBS.tau = 2042(+/-24)-557(+/-10)betabar. If confirmed, this provides a challenging constraint on the nature of interstellar dust in M31.

We have analysed the Herschel and SCUBA-2 dust continuum observations of the main filament in the Taurus L1495 star forming region, using the Bayesian fitting procedure PPMAP. (i) If we construct an average profile along the whole length of the filament, it has fwhm $\simeq 0.087\pm 0.003\,{\rm pc};\;$, but the closeness to previous estimates is coincidental. (ii) If we analyse small local sections of the filament, the column-density profile approximates well to the form predicted for hydrostatic equilibrium of an isothermal cylinder. (iii) The ability of PPMAP to distinguish dust emitting at different temperatures, and thereby to discriminate between the warm outer layers of the filament and the cold inner layers near the spine, leads to a significant reduction in the surface-density, $\varSigma$, and hence in the line-density, $\mu$. If we adopt the canonical value for the critical line-density at a gas-kinetic temperature of $10\,{\rm K}$, $\mu_{_{\rm CRIT}}\simeq 16\,{\rm M_{_\odot}\,pc^{-1}}$, the filament is on average trans-critical, with ${\bar\mu}\sim \mu_{_{\rm CRIT}};\;$ local sections where $\mu >\mu_{_{\rm CRIT}}$ tend to lie close to pre-stellar cores. (iv) The ability of PPMAP to distinguish different types of dust, i.e. dust characterised by different values of the emissivity index, $\beta$, reveals that the dust in the filament has a lower emissivity index, $\beta\leq1.5$, than the dust outside the filament, $\beta\geq 1.7$, implying that the physical conditions in the filament have effected a change in the properties of the dust.

Feedback from high-mass stars shapes the ISM of galaxies and thereby impacts gas that will form future generations of stars. However, due to our inability to track the time evolution of individual molecular clouds, quantifying the exact role of feedback on their star formation history is an observationally challenging task. In the present study, we take advantage of the unique properties of the G316.75-00.00 high-mass star-forming ridge to determine how feedback from O-stars impacts the dynamical stability of filaments. G316.75 is a 13.6pc long ridge containing 18,900Msun of H2 gas which is half IR dark, half IR bright. The IR bright half has already formed 4 O-type stars over the past 2Myr, while the IR dark half is still quiescent. Therefore, by assuming the gas properties of the dark part represents an earlier evolutionary stage of the bright part, we can quantify how feedback impacts these properties by contrasting them. We use archive Herschel and molecular line data, tracing both dense (NH3 and N2H+) and more diffuse (13CO) gas, to measure the ratio of kinetic to gravitational energy per-unit-length, virial-ratio-per-line, across the ridge for a range of gas volume densities. We show that despite the presence of 4 embedded O-stars, feedback cannot unbind the ridge's mass except for some small gas pockets near the O-stars. In fact, the virial-ratio-per-line is almost indistinguishable for both parts of the ridge. These results are at odds with most simulations where O-star-forming clouds are dispersed by feedback within a few cloud free-fall time. We conclude that the discrepancy between such simulations and our observations originates from different cloud morphologies and average densities when the first O-star forms. These results have important implications regarding, for instance, how feedback is implemented within cosmological and galactic-scale simulations.

We present resolution-enhanced images of warm dust at multiple temperatures and opacity index values in the star-forming bubble/HII region, RCW 120. The image set, representing a 4D hypercube of differential column density, was obtained using our Bayesian procedure, ppmap. The cool peripheral material ($\sim16$-22 K) exhibits ragged clumpy structure as noted previously by others. However, at higher temperatures ($\stackrel{>}{_\sim}26$ K) the geometry changes dramatically, showing a bubble boundary which is accurately circular in projection, except for the previously-reported opening in the north. Comparison with Spitzer 8 $\mu$m data suggests that the $\sim26$-30 K dust seen by Herschel resides in the photodissociation region (PDR) surrounding the HII region. Its projected radial profile is consistent with that of a spherical shell, thus arguing against previous suggestions of cylindrical or planar geometry. The inferred geometry is, in fact, consistent with previous interpretations of the HII region as a classical Str\"omgren sphere, except for the fact that the ionising star (CD -38.11636; O8V) is displaced by more than half a radius from its geometric centre. None of the previously published models has satisfactorily accounted for that displacement. It could, however, be explained by proper motion of the O star at $\sim2$-4 km s$^{-1}$ since its formation, possibly due to a cloud-cloud collision. We suggest that the current spherical bubble constitutes the fossilised remnant of the initial expansion of the HII region following the formation of the star, which now continues to flee its formation site.

We search for far-infrared (FIR) counterparts of known supernova remnants (SNRs) in the Galactic plane (10 deg <| l |< 60 deg) at 70-500 micron using the Herschel Infrared Galactic Plane Survey (Hi-GAL). Of 71 sources studied, we find that 29 (41 per cent) SNRs have a clear FIR detection of dust emission associated with the SNR. Dust from 8 of these is in the central region, and 4 indicate pulsar wind nebulae (PWNe) heated ejecta dust. A further 23 have dust emission in the outer shell structures which is potentially related to swept up material. Many Galactic SNe have dust signatures but we are biased towards detecting ejecta dust in young remnants and those with a heating source (shock or PWN). We estimate the dust temperature and mass contained within three PWNe, G11.2-0.3, G21.5-0.9, and G29.7-0.3 using modified blackbody fits. To more rigorously analyse the dust properties at various temperatures and dust emissivity index beta, we use point process mapping (PPMAP). We find significant quantities of cool dust (at 20-40 K) with dust masses of Md = 0.34 +/- 0.14 solar mass, Md = 0.29 +/- 0.08 solar mass, and Md = 0.51 +/- 0.13 solar mass for G11.2-0.3, G21.5-0.9, and G29.7-0.3 respectively. We derive the dust emissivity index for the PWN ejecta dust in G21.5-0.3 to be beta = 1.4 +/- 0.5 compared to dust in the surrounding medium where beta = 1.8 +/- 0.1.

Aims. To constrain the physical processes that lead to the birth of high-mass stars it is mandatory to study the very first stages of their formation. We search for high-mass analogs of low-mass prestellar cores in W43-MM1. Methods. We conducted a 1.3 mm ALMA mosaic of the complete W43-MM1 cloud, which has revealed numerous cores with ~ 2000 au FWHM sizes. We investigated the nature of cores located at the tip of the main filament, where the clustering is minimum. We used the continuum emission to measure the core masses and the $^{13}$CS(5-4) line emission to estimate their turbulence level. We also investigated the prestellar or protostellar nature of these cores by searching for outflow signatures traced by CO(2-1) and SiO(5-4) line emission, and for molecular complexity typical of embedded hot cores. Results. Two high-mass cores of ~ 1300 au diameter and ~ $60~M_\odot$ mass are observed to be turbulent but gravitationally bound. One drives outflows and is associated with a hot core. The other core, W43-MM1\#6, does not yet reveal any star formation activity and thus is an excellent high-mass prestellar core candidate.

We present new Herschel-derived images of warm dust in the Andromeda Galaxy, M31, with unprecedented spatial resolution (~ 30 pc), column density accuracy, and constraints on the three-dimensional distributions of dust temperature and dust opacity index (hence grain size and composition), based on the new PPMAP Bayesian analysis procedure. We confirm the overall radial variation of dust opacity index reported by other recent studies, including the central decrease within ~ 3 kpc of the nucleus. We also investigate the detailed distribution of dust in the nuclear region, a prominent feature of which is a ~ 500 pc bar-like structure seen previously in H{\alpha}. The nature of this feature has been the subject of some debate. Our maps show it to be the site of the warmest dust, with a mean line-of-sight temperature ~ 30 K. A comparison with the stellar distribution, based on 2MASS data, provides strong evidence that it is a gravitationally induced bar. A comparison with radial velocity maps suggests the presence of an inflow towards the nucleus from opposite directions along this bar, fed presumably by the nuclear spiral with which it appears to connect. Such behaviour is common in large-scale bars in spiral galaxies, as is the phenomenon of nested bars whereby a subkiloparsec nuclear bar exists within a large-scale primary bar. We suggest that M31 represents an example of such nesting.

Understanding the processes that determine the stellar Initial Mass Function (IMF) is a critical unsolved problem, with profound implications for many areas of astrophysics. In molecular clouds, stars are formed in cores, gas condensations which are sufficiently dense that gravitational collapse converts a large fraction of their mass into a star or small clutch of stars. In nearby star-formation regions, the core mass function (CMF) is strikingly similar to the IMF, suggesting that the shape of the IMF may simply be inherited from the CMF. Here we present 1.3 mm observations, obtained with ALMA, the world's largest interferometer, of the active star-formation region W43-MM1, which may be more representative of the Galactic-disk regions where most stars form. The unprecedented resolution of these observations reveals, for the first time, a statistically robust CMF at high masses, with a slope that is markedly shallower than the IMF. This seriously challenges our understanding of the origin of the IMF.

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

We present far-infrared observations of Monoceros R2 (a giant molecular cloud at approximately 830 pc distance, containing several sites of active star formation), as observed at 70 {\mu}m, 160 {\mu}m, 250 {\mu}m, 350 {\mu}m, and 500 {\mu}m by the Photodetector Array Camera and Spectrometer (PACS) and Spectral and Photometric Imaging Receiver (SPIRE) instruments on the Herschel Space Observatory as part of the Herschel imaging survey of OB young stellar objects (HOBYS) Key programme. The Herschel data are complemented by SCUBA-2 data in the submillimetre range, and WISE and Spitzer data in the mid-infrared. In addition, C18O data from the IRAM 30-m Telescope are presented, and used for kinematic information. Sources were extracted from the maps with getsources, and from the fluxes measured, spectral energy distributions were constructed, allowing measurements of source mass and dust temperature. Of 177 Herschel sources robustly detected in the region (a detection with high signal-to-noise and low axis ratio at multiple wavelengths), including protostars and starless cores, 29 are found in a filamentary hub at the centre of the region (a little over 1% of the observed area). These objects are on average smaller, more massive, and more luminous than those in the surrounding regions (which together suggest that they are at a later stage of evolution), a result that cannot be explained entirely by selection effects. These results suggest a picture in which the hub may have begun star formation at a point significantly earlier than the outer regions, possibly forming as a result of feedback from earlier star formation. Furthermore, the hub may be sustaining its star formation by accreting material from the surrounding filaments.

We describe new Hi-GAL based maps of the entire Galactic Plane, obtained using continuum data in the wavelength range 70-500 $\mu$m. These maps are derived with the PPMAP procedure, and therefore represent a significant improvement over those obtained with standard analysis techniques. Specifically they have greatly improved resolution (12 arcsec) and, in addition to more accurate integrated column densities and mean dust temperatures, they give temperature-differential column densities, i.e., separate column density maps in twelve distinct dust temperature intervals, along with the corresponding uncertainty maps. The complete set of maps is available on-line. We briefly describe PPMAP and present some illustrative examples of the results. These include (a) multi-temperature maps of the Galactic HII region W5-E, (b) the temperature decomposition of molecular cloud column-density probability distribution functions, and (c) the global variation of mean dust temperature as a function of Galactocentric distance. Amongst our findings are: (i) a strong localised temperature gradient in W5-E in a direction orthogonal to that towards the ionising star, suggesting an alternative heating source and providing possible guidance for models of the formation of the bubble complex, and (ii) the overall radial profile of dust temperature in the Galaxy shows a monotonic decrease, broadly consistent both with models of the interstellar radiation field and with previous estimates at lower resolution. However, we also find a central temperature plateau within ~ 6 kpc of the Galactic centre, outside of which is a pronounced steepening of the radial profile. This behaviour may reflect the greater proportion of molecular (as opposed to atomic) gas in the central region of the Galaxy.

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

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

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

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

We present a comparison of SCUBA-2 850-$\mu$m and Herschel 70--500-$\mu$m observations of the L1495 filament in the Taurus Molecular Cloud with the goal of characterising the SCUBA-2 Gould Belt Survey (GBS) data set. We identify and characterise starless cores in three data sets: SCUBA-2 850-$\mu$m, Herschel 250-$\mu$m, and Herschel 250-$\mu$m spatially filtered to mimic the SCUBA-2 data. SCUBA-2 detects only the highest-surface-brightness sources, principally detecting protostellar sources and starless cores embedded in filaments, while Herschel is sensitive to most of the cloud structure, including extended low-surface-brightness emission. Herschel detects considerably more sources than SCUBA-2 even after spatial filtering. We investigate which properties of a starless core detected by Herschel determine its detectability by SCUBA-2, and find that they are the core's temperature and column density (for given dust properties). For similar-temperature cores, such as those seen in L1495, the surface brightnesses of the cores are determined by their column densities, with the highest-column-density cores being detected by SCUBA-2. For roughly spherical geometries, column density corresponds to volume density, and so SCUBA-2 selects the densest cores from a population at a given temperature. This selection effect, which we quantify as a function of distance, makes SCUBA-2 ideal for identifying those cores in Herschel catalogues that are closest to forming stars. Our results can now be used by anyone wishing to use the SCUBA-2 GBS data set.

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

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

Image cubes of differential column density as a function of dust temperature are constructed for Galactic Centre molecular cloud G0.253+0.016 ("The Brick") using the recently described PPMAP procedure. The input data consist of continuum images from the Herschel Space Telescope in the wavelength range 70-500 $\mu$m, supplemented by previously published interferometric data at 1.3 mm wavelength. While the bulk of the dust in the molecular cloud is consistent with being heated externally by the local interstellar radiation field, our image cube shows the presence, near one edge of the cloud, of a filamentary structure whose temperature profile suggests internal heating. The structure appears as a cool ($\sim 14$ K) tadpole-like feature, $\sim 6$ pc in length, in which is embedded a thin spine of much hotter ($\sim$ 40-50 K) material. We interpret these findings in terms of a cool filament whose hot central region is undergoing gravitational collapse and fragmentation to form a line of protostars. If confirmed, this would represent the first evidence of widespread star formation having started within this cloud.

We present a catalogue of dense cores in a $\sim 4^\circ\times2^\circ$ field of the Taurus star-forming region, inclusive of the L1495 cloud, derived from Herschel SPIRE and PACS observations in the 70 $\mu$m, 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, and 500 $\mu$m continuum bands. Estimates of mean dust temperature and total mass are derived using modified blackbody fits to the spectral energy distributions. We detect 525 starless cores of which $\sim10$-20% are gravitationally bound and therefore presumably prestellar. Our census of unbound objects is $\sim85$% complete for $M>0.015\,M_\odot$ in low density regions ($A_V\stackrel{<}{_\sim}5$ mag), while the bound (prestellar) subset is $\sim85$% complete for $M>0.1\,M_\odot$ overall. The prestellar core mass function (CMF) is consistent with lognormal form, resembling the stellar system initial mass function, as has been reported previously. All of the inferred prestellar cores lie on filamentary structures whose column densities exceed the expected threshold for filamentary collapse, in agreement with previous reports. Unlike the prestellar CMF, the unbound starless CMF is not lognormal, but instead is consistent with a power-law form below $0.3\,M_\odot$ and shows no evidence for a low-mass turnover. It resembles previously reported mass distributions for CO clumps at low masses ($M\stackrel{<}{_\sim}0.3\,M_\odot$). The volume density PDF, however, is accurately lognormal except at high densities. It is consistent with the effects of self-gravity on magnetized supersonic turbulence. The only significant deviation from lognormality is a high-density tail which can be attributed unambiguously to prestellar cores.

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

We present PPMAP, a Bayesian procedure that uses images of dust continuum emission at multiple wavelengths to produce resolution-enhanced image cubes of differential column-density as a function of dust temperature and position. PPMAP is based on the generic 'point process' formalism, whereby the system of interest (in this case, a dusty astrophysical structure such as a filament or prestellar core) is represented by a collection of points in a suitably defined state space. It can be applied to a variety of observational data, such as Herschel images, provided only that the image intensity is delivered by optically thin dust in thermal equilibrium. PPMAP takes full account of the instrumental point spread functions and does not require all images to be degraded to the same resolution. We present the results of testing using simulated data for a prestellar core and a fractal turbulent cloud, and demonstrate its performance with real data from the Hi-GAL survey. Specifically, we analyse observations of a large filamentary structure in the CMa OB1 giant molecular cloud. Histograms of differential column-density indicate that the warm material (T > 13 K) is distributed log-normally, consistent with turbulence, but the column-densities of the cooler material are distributed as a high density tail, consistent with the effects of self-gravity. The results illustrate the potential of PPMAP to aid in distinguishing between different physical components along the line of sight in star-forming clouds, and aid the interpretation of the associated PDFs of column density.

We present and discuss the results of the Herschel Gould Belt survey observations in a ~11 deg^2 area of the Aquila molecular cloud complex at d~260 pc, imaged with the SPIRE/PACS cameras from 70 to 500 micron. We identify a complete sample of starless dense cores and embedded protostars in this region, and analyze their global properties and spatial distributions. We find a total of 651 starless cores, ~60% of which are gravitationally bound prestellar cores, and they will likely form stars in the future. We also detect 58 protostellar cores. The core mass function (CMF) derived for the prestellar cores is very similar in shape to the stellar initial mass function (IMF), supporting the earlier view that there is a close physical link between the IMF and the CMF. The global shift in mass scale observed between the CMF and the IMF is consistent with a typical star formation efficiency of ~40%. By comparing the numbers of starless cores to the number of young stellar objects, we estimate that the lifetime of prestellar cores is ~1 Myr. We find a strong correlation between the spatial distribution of prestellar cores and the densest filaments. About 90% of the Herschel-identified prestellar cores are located above a background column density corresponding to A_V~7, and ~75% of them lie within filamentary structures with supercritical masses per unit length >~16 M_sun/pc. These findings support a picture wherein the cores making up the peak of the CMF (and probably responsible for the base of the IMF) result primarily from the gravitational fragmentation of marginally supercritical filaments. Given that filaments appear to dominate the mass budget of dense gas at A_V>7, our findings also suggest that the physics of prestellar core formation within filaments is responsible for a characteristic "efficiency" SFR/M_dense ~5+-2 x 10^-8 yr^-1 for the star formation process in dense gas.

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

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

We present a sample of brown dwarfs identified with the {\it Wide-field Infrared Survey Explorer} (WISE) for which we have obtained {\it Hubble Space Telescope} ({\it HST}) Wide Field Camera 3 (WFC3) near-infrared grism spectroscopy. The sample (twenty-two in total) was observed with the G141 grism covering 1.10$-$1.70 $\mu$m, while fifteen were also observed with the G102 grism, which covers 0.90$-$1.10 $\mu$m. The additional wavelength coverage provided by the G102 grism allows us to 1) search for spectroscopic features predicted to emerge at low effective temperatures (e.g.\ ammonia bands) and 2) construct a smooth spectral sequence across the T/Y boundary. We find no evidence of absorption due to ammonia in the G102 spectra. Six of these brown dwarfs are new discoveries, three of which are found to have spectral types of T8 or T9. The remaining three, WISE J082507.35$+$280548.5 (Y0.5), WISE J120604.38$+$840110.6 (Y0), and WISE J235402.77$+$024015.0 (Y1) are the nineteenth, twentieth, and twenty-first spectroscopically confirmed Y dwarfs to date. We also present {\it HST} grism spectroscopy and reevaluate the spectral types of five brown dwarfs for which spectral types have been determined previously using other instruments.

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

The AllWISE processing pipeline has measured motions for all objects detected on WISE images taken between 2010 January and 2011 February. In this paper, we discuss new capabilities made to the software pipeline in order to make motion measurements possible, and we characterize the resulting data products for use by future researchers. Using a stringent set of selection criteria, we find 22,445 objects that have significant AllWISE motions, of which 3,525 have motions that can be independently confirmed from earlier 2MASS images yet lack any published motions in SIMBAD. Another 58 sources lack 2MASS counterparts and are presented as motion candidates only. Limited spectroscopic follow-up of this list has already revealed eight new L subdwarfs. These may provide the first hints of a "subdwarf gap" at mid-L types that would indicate the break between the stellar and substellar populations at low metallicities (i.e., old ages). Another object in the motion list -- WISEA J154045.67-510139.3 -- is a bright (J ~ 9 mag) object of type M6; both the spectrophotometric distance and a crude preliminary parallax place it ~6 pc from the Sun. We also compare our list of motion objects to the recently published list of 762 WISE motion objects from Luhman (2014). While these first large motion studies with WISE data have been very successful in revealing previously overlooked nearby dwarfs, both studies missed objects that the other found, demonstrating that many other nearby objects likely await discovery in the AllWISE data products.

The density and temperature structures of dense cores in the L1495 cloud of the Taurus star-forming region are investigated using Herschel SPIRE and PACS images in the 70 $\mu$m, 160 $\mu$m, 250 $\mu$m, 350 $\mu$m and 500 $\mu$m continuum bands. A sample consisting of 20 cores, selected using spectral and spatial criteria, is analysed using a new maximum likelihood technique, COREFIT, which takes full account of the instrumental point spread functions. We obtain central dust temperatures, $T_0$, in the range 6-12 K and find that, in the majority of cases, the radial density falloff at large radial distances is consistent with the $r^{-2}$ variation expected for Bonnor-Ebert spheres. Two of our cores exhibit a significantly steeper falloff, however, and since both appear to be gravitationally unstable, such behaviour may have implications for collapse models. We find a strong negative correlation between $T_0$ and peak column density, as expected if the dust is heated predominantly by the interstellar radiation field. At the temperatures we estimate for the core centres, carbon-bearing molecules freeze out as ice mantles on dust grains, and this behaviour is supported here by the lack of correspondence between our estimated core locations and the previously-published positions of H$^{13}$CO$^+$ peaks. On this basis, our observations suggest a sublimation-zone radius typically $\sim 10^4$ AU. Comparison with previously-published N$_2$H$^+$ data at 8400 AU resolution, however, shows no evidence for N$_2$H$^+$ depletion at that resolution.

Utilizing multi-wavelength dust emission maps acquired with $Herschel$, we reconstruct local volume density and dust temperature profiles for the prestellar cores B68 and L1689B using inverse-Abel transform based technique. We present intrinsic radial dust temperature profiles of starless cores directly from dust continuum emission maps disentangling the effect of temperature variations along the line of sight which was previously limited to the radiative transfer calculations. The reconstructed dust temperature profiles show a significant drop in core center, a flat inner part, and a rising outward trend until the background cloud temperature is reached. The central beam-averaged dust temperatures obtained for B68 and L1689B are 9.3 $\pm$ 0.5 K and 9.8 $\pm$0.5 K, respectively, which are lower than the temperatures of 11.3 K and 11.6 K obtained from direct SED fitting. The best mass estimates derived by integrating the volume density profiles of B68 and L1689B are 1.6 M_sol and 11 M_sol, respectively. Comparing our results for B68 with the near-infrared extinction studies, we find that the dust opacity law adopted by the HGBS project, $\kappa_{\lambda} =0.1(\lambda/300 \mu m)^{-2}$, agrees to within 50% with the dust extinction constraints

In our effort to complete the census of low-mass stars and brown dwarfs in the immediate Solar Neighborhood, we present spectra, photometry, proper motions, and distance estimates for forty-two low-mass star and brown dwarf candidates discovered by the Wide-field Infrared Survey Explorer (WISE). We also present additional follow-up information on twelve candidates selected using WISE data but previously published elsewhere. The new discoveries include fifteen M dwarfs, seventeen L dwarfs, five T dwarfs, and five objects of other type. Among these discoveries is a newly identified "unusually red L dwarf" (WISE J223527.07+451140.9), four peculiar L dwarfs whose spectra are most readily explained as unresolved L+T binary systems, and a T9 dwarf (WISE J124309.61+844547.8). We also show that the recently discovered red L dwarf WISEP J004701.06+680352.1 (Gizis et al. 2012) may be a low-gravity object and hence young and potentially low mass (< 25 MJup).

We report the discovery of 87 new T dwarfs uncovered with the Wide-field Infrared Survey Explorer (WISE) and three brown dwarfs with extremely red near-infrared colors that exhibit characteristics of both L and T dwarfs. Two of the new T dwarfs are likely binaries with L7+/-1 primaries and mid-type T secondaries. In addition, our follow-up program has confirmed 10 previously identified T dwarfs and four photometrically-selected L and T dwarf candidates in the literature. This sample, along with the previous WISE discoveries, triples the number of known brown dwarfs with spectral types later than T5. Using the WISE All-Sky Source Catalog we present updated color-color and color-type diagrams for all the WISE-discovered T and Y dwarfs. Near-infrared spectra of the new discoveries are presented, along with spectral classifications. To accommodate later T dwarfs we have modified the integrated flux method of determining spectral indices to instead use the median flux. Furthermore, a newly defined J-narrow index differentiates the early-type Y dwarfs from late-type T dwarfs based on the J-band continuum slope. The K/J indices for this expanded sample show that 32% of late-type T dwarfs have suppressed K-band flux and are blue relative to the spectral standards, while only 11% are redder than the standards. Comparison of the Y/J and K/J index to models suggests diverse atmospheric conditions and supports the possible re-emergence of clouds after the L/T transition. We also discuss peculiar brown dwarfs and candidates that were found not to be substellar, including two Young Stellar Objects and two Active Galactic Nuclei. The coolest WISE-discovered brown dwarfs are the closest of their type and will remain the only sample of their kind for many years to come.

The Wide-field Infrared Survey Explorer has revealed a T8.5 brown dwarf (WISE J111838.70+312537.9) that exhibits common proper motion with a solar-neighborhood (8 pc) quadruple star system - Xi Ursae Majoris. The angular separation is 8.5 arc-min, and the projected physical separation is about 4000 AU. The sub-solar metallicity and low chromospheric activity of Xi UMa A argue that the system has an age of at least 2 Gyr. The infrared luminosity and color of the brown dwarf suggests the mass of this companion ranges between 14 and 38 Jupiter masses for system ages of 2 and 8 Gyr respectively.

We have monitored the position of the cool Y dwarf WISEPA J182831.08+265037.8 using a combination of ground- and space-based telescopes and have determined its distance to be 11.2$_{-1.0}^{+1.3}$ pc. Its absolute H magnitude, M$_H=22.21^{+0.25}_{-0.22}$ mag, suggests a mass in the range 0.5-20 M$_{Jup}$ for ages of 0.1-10 Gyr with an effective temperature in the range 250-400 K. The broad range in mass is due primarily to the unknown age of the object. Since the high tangential velocity of the object, 51$\pm5$ km s$^{-1}$, is characteristic of an old disk population, a plausible age range of 2-4 Gyr leads to a mass range of 3-6 M$_{Jup}$ based on fits to the (highly uncertain) COND evolutionary models. The range in temperature is due to the fact that no single model adequately represents the 1-5 $\mu$m spectral energy distribution (SED) of the source, failing by factors of up to 5 at either the short or long wavelength portions of the spectral energy distribution. The appearance of this very cold object may be affected by non-equilibrium chemistry or low temperature condensates forming clouds, two atmospheric processes that are known to be important in brown dwarf atmospheres but have proven difficult to model. Finally, we argue that there would have to be a very steep upturn in the number density of late type Y-dwarfs to account for the putative population of objects suggested by by recent microlensing observations. Whether WISE 1828+2650 sits at the low mass end of the brown dwarf population or is the first example of a large number of "free-floating" planets is not yet known.

We present astrometric measurements of eleven nearby ultracool brown dwarfs of spectral types Y and late-T, based on imaging observations from a variety of space-based and ground-based telescopes. These measurements have been used to estimate relative parallaxes and proper motions via maximum likelihood fitting of geometric model curves. To compensate for the modest statistical significance (<~ 7) of our parallax measurements we have employed a novel Bayesian procedure for distance estimation which makes use of an a priori distribution of tangential velocities, Vtan, assumed similar to that implied by previous observations of T dwarfs. Our estimated distances are therefore somewhat dependent on that assumption. Nevertheless, the results have yielded distances for five of our eight Y dwarfs and all three T dwarfs. Estimated distances in all cases are >~ 3 pc. In addition, we have obtained significant estimates of Vtan for two of the Y dwarfs; both are <100 km/s, consistent with membership in the thin disk population. Comparison of absolute magnitudes with model predictions as a function of color shows that the Y dwarfs are significantly redder in J-H than predicted by a cloud-free model

We present a census of the disk population for UKIDSS selected brown dwarfs in the 5-10 Myr old Upper Scorpius OB association. For 116 objects originally identified in UKIDSS, the majority of them not studied in previous publications, we obtain photometry from the WISE database. The resulting colour-magnitude and colour-colour plots clearly show two separate populations of objects, interpreted as brown dwarfs with disks (class II) and without disks (class III). We identify 27 class II brown dwarfs, 14 of them not previously known. This disk fraction (27 out of 116 or 23%) among brown dwarfs was found to be similar to results for K/M stars in Upper Scorpius, suggesting that the lifetimes of disks are independent of the mass of the central object for low-mass stars and brown dwarfs. 5 out of 27 disks (19%) lack excess at 3.4 and 4.6 microns and are potential transition disks (i.e. are in transition from class II to class III). The transition disk fraction is comparable to low-mass stars. We estimate that the timescale for a typical transition from class II to class III is less than 0.4 Myr for brown dwarfs. These results suggest that the evolution of brown dwarf disks mirrors the behaviour of disks around low-mass stars, with disk lifetimes on the order of 5-10 Myr and a disk clearing timescale significantly shorter than 1 Myr.

We present Spitzer 3.6 and 4.5 $\mu$m photometry and positions for a sample of 1510 brown dwarf candidates identified by the WISE all-sky survey. Of these, 166 have been spectroscopically classified as objects with spectral types M(1), L(7), T(146), and Y(12); Sixteen other objects are non-(sub)stellar in nature. The remainder are most likely distant L and T dwarfs lacking spectroscopic verification, other Y dwarf candidates still awaiting follow-up, and assorted other objects whose Spitzer photometry reveals them to be background sources. We present a catalog of Spitzer photometry for all astrophysical sources identified in these fields and use this catalog to identify 7 fainter (4.5 $\mu$m $\sim$ 17.0 mag) brown dwarf candidates, which are possibly wide-field companions to the original WISE sources. To test this hypothesis, we use a sample of 919 Spitzer observations around WISE-selected high-redshift hyper-luminous infrared galaxy (HyLIRG) candidates. For this control sample we find another 6 brown dwarf candidates, suggesting that the 7 companion candidates are not physically associated. In fact, only one of these 7 Spitzer brown dwarf candidates has a photometric distance estimate consistent with being a companion to the WISE brown dwarf candidate. Other than this there is no evidence for any widely separated ($>$ 20 AU) ultra-cool binaries. As an adjunct to this paper, we make available a source catalog of $\sim$ 7.33 $\times 10^5$ objects detected in all of these Spitzer follow-up fields for use by the astronomical community. The complete catalog includes the Spitzer 3.6 and 4.5 $\mu$m photometry, along with positionally matched $B$ and $R$ photometry from USNO-B; $J$, $H$, and $K_s$ photometry from 2MASS; and $W1$, $W2$, $W3$, and $W4$ photometry from the WISE all-sky catalog.

We present a general method for identifying the pre-main-sequence population of any star-forming region, unbiased with respect to the presence or absence of disks, in contrast to samples selected primarily via their mid-infrared emission from Spitzer surveys. We have applied this technique to a new, deep, wide-field, near-infrared imaging survey of the Rho Ophiuchi cloud core to search for candidate low mass members. In conjunction with published Spitzer IRAC photometry, and least squares fits of model spectra (COND, DUSTY, NextGen, and blackbody) to the observed spectral energy distributions, we have identified 948 candidate cloud members within our 90% completeness limits of J=20.0, H=20.0, and K_S=18.50. This population represents a factor of ~3 increase in the number of known young stellar objects in the Rho Ophiuchi cloud. A large fraction of the candidate cluster members (81% +/- 3%) exhibit infrared excess emission consistent with the presence of disks, thus strengthening the possibility of their being bona fide cloud members. Spectroscopic follow-up will confirm the nature of individual objects, better constrain their parameters, and allow an initial mass function to be derived.

We present the discovery of another seven Y dwarfs from the Wide-field Infrared Survey Explorer (WISE). Using these objects, as well as the first six WISE Y dwarf discoveries from Cushing et al., we further explore the transition between spectral types T and Y. We find that the T/Y boundary roughly coincides with the spot where the J-H colors of brown dwarfs, as predicted by models, turn back to the red. Moreover, we use preliminary trigonometric parallax measurements to show that the T/Y boundary may also correspond to the point at which the absolute H (1.6 um) and W2 (4.6 um) magnitudes plummet. We use these discoveries and their preliminary distances to place them in the larger context of the Solar Neighborhood. We present a table that updates the entire stellar and substellar constituency within 8 parsecs of the Sun, and we show that the current census has hydrogen-burning stars outnumbering brown dwarfs by roughly a factor of six. This factor will decrease with time as more brown dwarfs are identified within this volume, but unless there is a vast reservoir of cold brown dwarfs invisible to WISE, the final space density of brown dwarfs is still expected to fall well below that of stars. We also use these new Y dwarf discoveries, along with newly discovered T dwarfs from WISE, to investigate the field substellar mass function. We find that the overall space density of late-T and early-Y dwarfs matches that from simulations describing the mass function as a power law with slope -0.5 < alpha < 0.0; however, a power-law may provide a poor fit to the observed object counts as a function of spectral type because there are tantalizing hints that the number of brown dwarfs continues to rise from late-T to early-Y. More detailed monitoring and characterization of these Y dwarfs, along with dedicated searches aimed at identifying more examples, are certainly required.

The Wide-field Infrared Survey Explorer Preliminary Data Release Source Catalog contains over 257 million objects. We describe the method used to flag variable source candidates in the Catalog. Using a method based on the chi- square of single-exposure flux measurements, we generated a variability flag for each object, and have identified almost 460,000 candidates sources that exhibit significant flux variability with greater than \sim 7{\sigma} confidence. We discuss the flagging method in detail and describe its benefits and limitations. We also present results from the flagging method, including example light curves of several types of variable sources including Algol-type eclipsing binaries, RR Lyr, W UMa, and a blazar candidate.

We discuss the optimal detection of point sources from multiwavelength imaging data using an approach, referred to as MDET, which requires no prior knowledge of the source spectrum. MDET may be regarded as a somewhat more general version of the so-called "chi squared" technique. We describe the theoretical basis of the technique, and show examples of its performance with four-channel infrared broad-band imaging data from the WISE mission. We also discuss the potential benefits of applying it to the multifrequency data cubes of the ASKAP surveys, and suggest that it could increase the detection sensitivity of searches for neutral hydrogen emission at moderately high redshifts.

We present ground-based spectroscopic verification of six Y dwarfs (see Cushing et al), eighty-nine T dwarfs, eight L dwarfs, and one M dwarf identified by the Wide-field Infrared Survey Explorer (WISE). Eighty of these are cold brown dwarfs with spectral types greater than or equal to T6, six of which have been announced earlier in Mainzer et al and Burgasser et al. We present color-color and color-type diagrams showing the locus of M, L, T, and Y dwarfs in WISE color space. Near-infrared classifications as late as early Y are presented and objects with peculiar spectra are discussed. After deriving an absolute WISE 4.6 um (W2) magnitude vs. spectral type relation, we estimate spectrophotometric distances to our discoveries. We also use available astrometric measurements to provide preliminary trigonometric parallaxes to four our discoveries, which have types of L9 pec (red), T8, T9, and Y0; all of these lie within 10 pc of the Sun. The Y0 dwarf, WISE 1541-2250, is the closest at 2.8 (+1.3,-0.6) pc; if this 2.8 pc value persists after continued monitoring, WISE 1541-2250 will become the seventh closest stellar system to the Sun. Another ten objects, with types between T6 and >Y0, have spectrophotometric distance estimates also placing them within 10 pc. The closest of these, the T6 dwarf WISE 1506+7027, is believed to fall at a distance of roughly 4.9 pc. WISE multi-epoch positions supplemented with positional info primarily from Spitzer/IRAC allow us to calculate proper motions and tangential velocities for roughly one half of the new discoveries. This work represents the first step by WISE to complete a full-sky, volume-limited census of late-T and Y dwarfs. Using early results from this census, we present preliminary, lower limits to the space density of these objects and discuss constraints on both the functional form of the mass function and the low-mass limit of star formation.

We present the discovery of seven ultracool brown dwarfs identified with the Wide-field Infrared Survey Explorer (WISE). Near-infrared spectroscopy reveals deep absorption bands of H_2O and CH_4 that indicate all seven of the brown dwarfs have spectral types later than UGPS J072227.51-054031.2, the latest type T dwarf currently known. The spectrum of WISEP J182831.08+265037.8 is distinct in that the heights of the J- and H-band peaks are approximately equal in units of f_lambda, so we identify it as the archetypal member of the Y spectral class. The spectra of at least two of the other brown dwarfs exhibit absorption on the blue wing of the H-band peak that we tentatively ascribe to NH_3. These spectral morphological changes provide a clear transition between the T dwarfs and the Y dwarfs. In order to produce a smooth near-infrared spectral sequence across the T/Y dwarf transition, we have reclassified UGPS J0722-0540 as the T9 spectral standard and tentatively assign WISEP J173835.52+273258.9 as the Y0 spectral standard. In total, six of the seven new brown dwarfs are classified as Y dwarfs: four are classified as Y0, one is classified as Y0 (pec?), and WISEP J1828+2650 is classified as >Y0. We have also compared the spectra to the model atmospheres of Marley and Saumon and infer that the brown dwarfs have effective temperatures ranging from 300 K to 500 K, making them the coldest spectroscopically confirmed brown dwarfs known to date.

The multiplicity properties of brown dwarfs are critical empirical constraints for formation theories, while multiples themselves provide unique opportunities to test evolutionary and atmospheric models and examine empirical trends. Studies using high-resolution imaging can not only uncover faint companions, but they can also be used to determine dynamical masses through long-term monitoring of binary systems. We have begun a search for the coolest brown dwarfs using preliminary processing of data from the Wide-field Infrared Survey Explorer (WISE) and have confirmed many of the candidates as late-type T dwarfs. In order to search for companions to these objects, we are conducting observations using the Laser Guide Star Adaptive Optics system on Keck II. Here we present the first results of that search, including a T5 binary with nearly equal mass components and a faint companion to a T8.5 dwarf with an estimated spectral type of T9.

We present the discovery of five late-type T dwarfs identified with the Wide-field Infrared Survey Explorer (WISE). Low-resolution near-infrared spectroscopy obtained with the Magellan Folded-port InfraRed Echellette (FIRE) reveal strong water and methane absorption in all five sources, and spectral indices and comparison to spectral templates indicate classifications ranging from T5.5 to T8.5:. The spectrum of the latest-type source, WISE J1812+2721, is an excellent match to that of the T8.5 companion brown dwarf Wolf 940B. WISE-based spectrophotometric distance estimates place these T dwarfs at 12-13 pc from the Sun, assuming they are single. Preliminary fits of the spectral data to the atmosphere models of Saumon & Marley indicate effective temperatures ranging from 600 K to 930 K, both cloudy and cloud-free atmospheres, and a broad range of ages and masses. In particular, two sources show evidence of both low surface gravity and cloudy atmospheres, tentatively supporting a trend noted in other young brown dwarfs and exoplanets. In contrast, the high proper motion T dwarf WISE J2018-7423 exhibits a suppressed K-band peak and blue spectrophotometric J-K colors indicative of an old, massive brown dwarf; however, it lacks the broadened Y-band peak seen in metal-poor counterparts. These results illustrate the broad diversity of low-temperature brown dwarfs that will be uncovered with WISE.

The \emph{Wide-field Infrared Survey Explorer} has surveyed the entire sky at four infrared wavelengths with greatly improved sensitivity and spatial resolution compared to its predecessors, the \emph{Infrared Astronomical Satellite} and the \emph{Cosmic Background Explorer}. NASA's Planetary Science Division has funded an enhancement to the \WISE\ data processing system called "NEOWISE" that allows detection and archiving of moving objects found in the \WISE\ data. NEOWISE has mined the \WISE\ images for a wide array of small bodies in our Solar System, including Near-Earth Objects (NEOs), Main Belt asteroids, comets, Trojans, and Centaurs. By the end of survey operations in February 2011, NEOWISE identified over 157,000 asteroids, including more than 500 NEOs and $\sim$120 comets. The NEOWISE dataset will enable a panoply of new scientific investigations.

We report the discovery of the first new ultra-cool brown dwarf found with the Wide-field Infrared Survey Explorer (WISE). The object's preliminary designation is WISEPC J045853.90+643451.9. Follow-up spectroscopy with the LUCIFER instrument on the Large Binocular Telescope indicates that it is a very late-type T dwarf with a spectral type approximately equal to T9. Fits to an IRTF/SpeX 0.8-2.5 micron spectrum to the model atmospheres of Marley and Saumon indicate an effective temperature of approximately 600 K as well as the presence of vertical mixing in its atmosphere. The new brown dwarf is easily detected by WISE, with a signal-to-noise ratio of ~36 at 4.6 microns. Current estimates place it at a distance of 6 to 10 pc. This object represents the first in what will likely be hundreds of nearby brown dwarfs found by WISE that will be suitable for follow up observations, including those with the James Webb Space Telescope. One of the two primary scientific goals of the WISE mission is to find the coolest, closest stars to our Sun; the discovery of this new brown dwarf proves that WISE is capable of fulfilling this objective.

A search for young substellar objects in the rho Oph cloud core region has been made using the deep-integration Combined Calibration Scan images of the 2MASS extended mission in J, H and Ks bands, and Spitzer IRAC images at 3.6, 4.5, 5.8 and 8.0 microns. The field of view of the combined observations was 1 deg x 9.3 arcmin, and the 5 sigma limiting magnitude at J was 20.5. Comparison of the observed SEDs with the predictions of the COND and DUSTY models, for an assumed age of 1 Myr, supports the identification of many of the sources with brown dwarfs, and enables the estimation of effective temperature, Teff. The cluster members are then readily distinguishable from background stars by their locations on a plot of flux density versus Teff. The range of estimated Teff extends down to ~ 750 K, suggesting the presence of objects of sub-Jupiter mass. The results also suggest that the mass function for the rho Oph cloud resembles that of the sigma Orionis cluster based on a recent study, with both rising towards lower masses. The other main result from our study is the apparent presence of a progressive blueward skew in the distribution of J-H and H-Ks colors, such that the blue end of the range becomes increasingly bluer with increasing magnitude. We suggest that this behavior might be understood in terms of the 'ejected stellar embryo' hypothesis, whereby some of the lowest-mass brown dwarfs could escape to locations close to the front edge of the cloud, and thereby be seen with less extinction.

We report the discovery of a young planetary-mass brown dwarf in the rho Oph cloud core. The object was identified as such with the aid of a 1.5-2.4 micron low-resolution spectrum obtained using the NIRC instrument on the Keck I telescope. Based on the COND model, the observed spectrum is consistent with a reddened (Av ~ 15-16) brown dwarf whose effective temperature is in the range 1200-1800 K. For an assumed age of 1 Myr, comparison with isochrones further constrains the temperature to ~ 1400 K and suggests a mass of ~ 2-3 Jupiter masses. The inferred temperature is suggestive of an early T spectral type, which is supported by spectral morphology consistent with weak methane absorption. Based on its inferred distance (~ 100 pc) and the presence of overlying visual absorption, it is very likely to be a rho Oph cluster member. In addition, given the estimated spectral type, it may be the youngest and least massive T dwarf found so far. Its existence suggests that the initial mass function for the rho Oph star-forming region extends well into the planetary-mass regime.

We present new details of the structure and morphology of the jets and outflows in HH46/47 as seen in Spitzer infrared images from IRAC and MIPS, reprocessed using the ``HiRes'' deconvolution technique. HiRes improves the visualization of spatial morphology by enhancing resolution (to sub-arcsec levels in IRAC bands) and removing the contaminating side lobes from bright sources. In addition to sharper views of previously reported bow shocks, we have detected: (i) the sharply-delineated cavity walls of the wide-angle biconical outflow, seen in scattered light on both sides of the protostar, (ii) several very narrow jet features at distances 400 AU to 0.1 pc from the star, and, (iii) compact emissions at MIPS 24 micron coincident with the jet heads, tracing the hottest atomic/ionic gas in the bow shocks.

The capability of the Terrestrial Planet Finder Interferometer (TPF-I) for planetary signal extraction, including both detection and spectral characterization, can be optimized by taking proper account of instrumental characteristics and astrophysical prior information. We have developed the Point Process Algorithm (PPA), a Bayesian technique for extracting planetary signals using the sine-chopped outputs of a dual nulling interferometer. It is so-called because it represents the system being observed as a set of points in a suitably-defined state space, thus providing a natural way of incorporating our prior knowledge of the compact nature of the targets of interest. It can also incorporate the spatial covariance of the exozodi as prior information which could help mitigate against false detections. Data at multiple wavelengths are used simultaneously, taking into account possible spectral variations of the planetary signals. Input parameters include the RMS measurement noise and the a priori probability of the presence of a planet. The output can be represented as an image of the intensity distribution on the sky, optimized for the detection of point sources. Previous approaches by others to the problem of planet detection for TPF-I have relied on the potentially non-robust identification of peaks in a "dirty" image, usually a correlation map. Tests with synthetic data suggest that the PPA provides greater sensitivity to faint sources than does the standard approach (correlation map + CLEAN), and will be a useful tool for optimizing the design of TPF-I.

We have used the SHARC II camera at Caltech Submillimeter Observatory to make 350 and 450 micron images of the Vega dust disk at spatial resolutions (FWHM) of 9.7" and 11.1", respectively. The images show a ring-like morphology (radius \~ 100 AU) with inhomogeneous structure that is qualitatively different from that previously reported at 850 microns and longer wavelengths. We attribute the 350/450 micron emission to a grain population whose characteristic size (~ 1 mm) is intermediate between that of the cm-sized grains responsible for emission longward of 850 microns and the much smaller grains (< 18 microns) in the extensive halo, visible at 70 microns, discussed by Su et al. (2005). We have combined our submillimeter images with Spitzer data at 70 microns to produce 2-d maps of line-of-sight optical depth (relative column density). These "tau maps" suggest that the mm-sized grains are located preferentially in three symmetrically-located concentrations. If so, then this structure could be understood in terms of the Wyatt (2003) model in which planetesimals are trapped in the mean motion resonances of a Neptune-mass planet at 65 AU, provided allowance is made for the spatial distribution of dust grains to differ from that of the parent planetesimals. The peaks of the tau maps are, in fact, located near the expected positions corresponding to the 4:3 resonance. If this identification is confirmed by future observations, it would resolve an ambiguity with regard to the location of the planet.

Results of a comprehensive, new, ground-based mid-infrared imaging survey of the young stellar population of the Rho Ophiuchi cloud are presented. Data were acquired at the Palomar 5-m and at the Keck 10-m telescopes with the MIRLIN and LWS instruments, at 0.25 arcsec and 0.25 arcsec resolutions, respectively. Of 172 survey objects, 85 were detected. Among the 22 multiple systems observed, 15 were resolved and their individual component fluxes determined. A plot of the frequency distribution of the detected objects with SED spectral slope shows that YSOs spend ~400,000 yr in the Flat Spectrum phase, clearing out their remnant infall envelopes. Mid-infrared variability is found among a significant fraction of the surveyed objects, and is found to occur for all SED classes with optically thick disks. Large-amplitude near-infrared variability, also found for all SED classes with optically thick disks, seems to occur with somewhat higher frequency at the earlier evolutionary stages. Although a general trend of mid-infrared excess and NIR veiling exists proceeding through SED classes, with Class I objects generally exhibiting K-veilings > 1, Flat Spectrum objects with K-veilings > 0.58, and Class III objects with K-veilings =0, Class II objects exhibit the widest range of K-band veiling values, 0-4.5. However, the highly variable value of veiling that a single source can exhibit in any of the SED classes in which active disk accretion can take place is striking, and is direct observational evidence for highly time-variable accretion activity in disks. Finally, by comparing mid-infrared vs. near-infrared excesses in a subsample with well-determined effective temperatures and extinction values, disk clearing mechanisms are explored. The results are consistent with disk clearing proceeding from the inside-out.

We have imaged the circumstellar disk of Fomalhaut at 350 microns wavelength, using SHARC II at the Caltech Submillimeter Observatory. The spatial resolution of the raw images (9 arcsec) has been enhanced by a factor of three using the HiRes deconvolution procedure. We find that at this wavelength and signal to noise ratio (~ 12), the observed morphology is that of a simple inclined ring (i ~ 70 deg), with little or no other apparent structure; this is the first observation that shows clearly the ring morphology of the disk. We have combined our 350 micron data with Spitzer Space Telescope images at 24, 70, and 160 microns in order to estimate the 2-dimensional spatial variation of relative column density ("tau map") using our DISKFIT procedure. The tau map is based on the following physical assumptions: (1) the wavelength variation of opacity is the same throughout the disk, (2) the radial variation of dust temperature is dictated by the energy balance of individual grains in the stellar radiation field, and (3) the vertical scale height of the disk follows a power-law radial variation. The results confirm the ring-like morphology, but also show that the geometric center is displaced from the star by about 8 AU and that the ring has an apocentric enhancement of approximately 14% in integrated column density. If we interpret the displacement in terms of elliptical orbital motion due to gravitational perturbation by an unseen planet, then the implied forced eccentricity is ~ 0.06; dynamical modeling then predicts an apocentric density enhancement consistent with that inferred from the tau map.

Images at 12.5 microns of nuclei of three nearby Seyfert galaxies -- NGC 1275, NGC 4151 and NGC 7469 -- have been obtained with the Keck 10-m Telescope. NGC 7469 is resolved and deconvolution delineates a structure <0.04"x0.08" or <13x26 pc at a position angle of 135deg. From a comparison with structure seen at millimeter wavelengths, this structure is interpreted as a disk aligned with the molecular gas in the central few hundred parsecs of the galaxy. NGC 1275 and NGC 4151 are not resolved; limits on the sizes of these nuclei are 0.08" and 0.16", corresponding to physical spatial scales of 28 and 10 pc. The lower limits to the brightness temperatures implied by these size limits and the measured flux densities are within ~50K of the 12 micron to 25 micron color temperatures of these systems as inferred from IRAS observations. The angular size limits are within a factor of 2.5 of the sizes required to spatially resolve thermal emission from dust heated by a central luminosity source. These sizes preclude significant contributions to the nuclear infrared emission from star forming regions.

Mid-infrared observations of the central source of NGC 1068 have been obtained with a spatial resolution in the deconvolved image of 0.1" (~ 7pc). The central source is extended by ~1" in the north-south direction but appears unresolved in the east-west direction over most of its length. About two-thirds of its flux can be ascribed to a core structure which is itself elongated north-south and does not show a distinct unresolved compact source. The source is strongly asymmetric, extending significantly further to the north than to the south. The morphology of the mid-infrared emission appears similar to that of the radio jet, and has features which correlate with the images in [O III]. Its 12.5-24.5 micron color temperature ranges from 215 to 260 K and does not decrease smoothly with distance from the core. Silicate absorption is strongest in the core and to the south and is small in the north. The core, apparently containing two-thirds of the bolometric luminosity of the inner 4" diameter area, may be explained by a thick, dusty torus near the central AGN viewed at an angle of ~65 deg to its plane. There are, however, detailed difficulties with existing models, especially the narrow east-west width of the thin extended mid-infrared "tongue" to the north of the core. We interpret the tongue as re-processed visual and ultraviolet radiation that is strongly beamed and that originates in the AGN.

We have obtained mid-infrared images of the nucleus of NGC 1068 from the Hale 5 m telescope at Mt. Palomar with diffraction-limited resolution and high sensitivity at 8.8, 10.3, and 12.5 microns. Deconvolved images show that the infrared emission extends north to south in the inner 2", consisting of a central peak, a component extending 1" north of the central source, a component extending 1" south of the central source, and several smaller structures located 1" to the northeast. The central peak is extended 0.4" N-S and unresolved (<= 0.2") E-W. We find that 50 +/- 5% of the flux emerges from the central 0.4" and that a single unresolved point source can account for only 27 +/- 5% of the total flux. However, if the central peak arises from optically thick emission, we estimate that the emitting region has a projected area >= 2 pc^2, and may thus contain a compact source such as a parsec-scale torus. We observe a correspondence between the northern extension and the northeastern sources appearing on the mid-infrared images and the [OIII] clouds A-C & E. We interpret the faint optical counterpart to the mid-infrared southern extension as being due to partial obscuration by the intervening disk of the host galaxy. The N-S extension of the mid-infrared emission coincides with one wall of the conical narrow line region and aligns with the N-S orientation of the radio jet close to the nucleus. We interpret the infrared emission as arising from optically thick dust lining the walls of the low density cavity formed by the radio jet and heated by radiation from the central source.