We started a multi-scale analysis of G202.3+2.5, an intertwined filamentary region of Monoceros OB1. In Paper I, we examined the distributions of dense cores and protostars and found enhanced star formation (SF) activity in the junction region of the filaments. In this second paper, we aim to unveil the connections between the core and filament evolutions, and between the filament dynamics and the global evolution of the cloud. We characterise the gas dynamics and energy balance using Herschel and WISE observations and molecular tracers observed with the IRAM 30m and TRAO 14m telescopes. The velocity field of the cloud is examined and velocity-coherent structures are put in perspective with the cloud environment. Two main velocity components (VCs) are revealed, well separated in the north and merged around the location of intense N2H+ emission where Paper I found the peak of SF activity. The relative position of the two VCs along the sightline, and the velocity gradient in N2H+ emission imply that the VCs have been undergoing collision for ~10^5 yrs. The dense gas where N2H+ is detected is interpreted as the compressed region between the two filaments, which corresponds to a high mass inflow rate of ~1e-3 Msun/yr and possibly leads to an increase in its SF efficiency. We identify a protostar in the junction region that possibly powers two crossed intermittent outflows. We show that the HII region around the nearby cluster NCG 2264 is still expanding and its role in the collision is examined. However, we cannot rule out the idea that the collision arises mostly from the global collapse of the cloud. The (sub-)filament-scale observables examined in this paper reveal a collision between G202.3+2.5 sub-structures and its probable role in feeding the cores in the junction region. One must now characterise the cloud morphology, its fragmentation, and magnetic field, all at high resolution.

Current theories and models attempt to explain star formation globally, from core scales to giant molecular cloud scales. A multi-scale observational characterisation of an entire molecular complex is necessary to constrain them. We investigate star formation in G202.3+2.5, a ~10x3 pc sub-region of the Monoceros OB1 cloud with a complex morphology harbouring interconnected filamentary structures. We aim to connect the evolution of cores and filaments in G202.3+2.5 with the global evolution of the cloud and to identify the engines of the cloud dynamics. In this first paper, the star formation activity is evaluated by surveying the distributions of dense cores and protostars, and their evolutionary state, as characterised using both infrared observations from the Herschel and WISE telescopes and molecular line observations with the IRAM 30-m telescope. We find ongoing star formation in the whole cloud, with a local peak in star formation activity around the centre of G202.3+2.5 where a chain of massive cores (10-50 Msun) forms a massive ridge (>150 Msun). All evolutionary stages from starless cores to Class II protostars are found in G202.3+2.5, including a possibly starless, large column density (8x10^{22} cm-2), and massive (52 Msun) core. All the core-scale observables examined in this paper point to an enhanced star formation activity centred on the junction between the three main branches of the ramified structure of G202.3+2.5. This suggests that the increased star-formation activity results from the convergence of these branches. To further investigate the origin of this enhancement, it is now necessary to extend the analysis to larger scales, in order to examine the relationship between cores, filaments and their environment. We address these points through the analysis of the dynamics of G202.3+2.5 in a joint paper.

The Planck Catalogue of Galactic Cold Clumps (PGCC) contains over 13000 sources detected based on their cold dust signature. They are believed to consist of a mixture of quiescent, pre-stellar, and already star-forming objects. We extracted PGCC-type objects from cloud simulations and examined their physical and polarisation properties. The comparison with the PGCC catalogue helps to characterise the PGCC sample and provides valuable tests for numerical simulations of interstellar medium. We used several MHD snapshots to define the density field of our models. Sub-millimetre images of the surface brightness and polarisation were obtained with radiative transfer calculations. We examined the statistics of synthetic cold clump catalogues and examined the variations of the clump polarisation fraction p. The clump sizes, aspect ratios, and temperatures in the synthetic catalogue are similar to the PGCC. The fluxes and column densities are smaller by a factor of a few. Rather than with an increased dust opacity, this could be explained by increasing the average column density of the models by a factor of two to three, close to N(H2)= 10^22 cm-2. When the line of sight is parallel to the mean magnetic field, the polarisation fraction tends to increase towards the clump centres, contrary to observations. When the field is perpendicular, the polarisation fraction tends to decrease towards the clumps, but the drop in $p$ is small (e.g. from p~8% to p~7%). Magnetic field geometry reduces the polarisation fraction in the simulated clumps by only \Delta p~1% on average. The larger drop seen towards the actual PGCC clumps suggests some loss of grain alignment in the dense medium, such as predicted by the radiative torque mechanism. The statistical study is not able to quantify dust opacity changes at the scale of the PGCC clumps.

The Vela-Puppis region is known to host the Vela OB2 association as well as several young clusters featuring OB and pre-main sequence stars. Several spatial and kinematic subgroups have been identified in the recent years. By grouping stars based on their positions and velocity we can address the quetion of the dynamical history of the region and the mechanisms that drove stellar formation. The Gaia DR2 astrometry and photometry enables us to characterise the 3D spatial and 3D kinematic distribution of young stars and to estimate the ages of the identified components. We use an unsupervised classification method to group stars based on their proper motions and parallax. We perform a study of the expansion rates of the different identified groups from 3D velocities, and from corrected tangential velocities. We make use of theoretical isochrones to estimate ages. The young stars can be separated into seven main groups of different ages and kinematical distribution. All groups are found to be expanding, although the expansion is mostly not isotropic. The size of the region, the age substructure, and the anistropic expansion rates are compatible with a prolonged period of star formation in a turbulent molecular cloud, and that the dispersion of the stars cannot be explained by gas expulsion alone.

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

The dust sub-millimetre polarisation of star-forming clouds carries information on dust and the role of magnetic fields in cloud evolution. With observations of a dense filamentary cloud G035.39-00.33, we aim to characterise the dust emission properties and the variations of the polarisation fraction. JCMT SCUBA-2/POL-2 data at 850um are combined with Planck 850um (353GHz) data to map polarisation fractions. With previous SCUBA-2 observations (450um and 850um) and Herschel data, the column densities are determined via modified blackbody fits and via radiative transfer modelling. Models are constructed to examine how the polarisation angles and fractions depend on potential magnetic field geometries and grain alignment. POL-2 data show clear changes in the magnetic field orientation. The filament has a peak column density of N(H2)~7 10^22 cm-2, a minimum dust temperature of T~12 K, and a mass of some 4300Msun for the area N(H2)> 5 10^21 cm-2. The estimated average value of the dust opacity spectral index is beta ~ 1.9. The ratio of sub-millimetre and J band optical depths is tau(250 um)/tau(J) ~ 2.5 10^-3, more than four times the typical values for diffuse medium. The polarisation fraction decreases as a function of column density to p ~ 1% in the central filament. Because of noise, the observed decrease of p(N) is significant only at N(H2)>2 10^22 cm-2. The observations suggest that the grain alignment is not constant. Although the data can be explained with a complete loss of alignment at densities above ~ 10^4 cm-3 or using the predictions of radiative torques alignment, the uncertainty of the field geometry and the spatial filtering of the SCUBA-2 data prevent strong conclusions. G035.39-00.33 shows strong signs of dust evolution and the low polarisation fraction is suggestive of a loss of polarised emission from its densest parts.

The Galactic Cold Cores (GCC) project has made Herschel observations of interstellar clouds where Planck detected compact sources of cold dust emission. Our aim is to characterise the structure of the clumps and their parent clouds. We also examine the accuracy to which the structure of dense clumps can be determined from submillimetre data. We use standard statistical methods to characterise the GCC fields. Clumps are extracted using column density thresholding and we construct for each field a three-dimensional radiative transfer (RT) model. These are used to estimate the relative radiation field intensities, clump stability, and the uncertainty of column density estimates. We examine the radial column density profiles of the clumps. In the GCC fields, the structure noise follows the relations previously established at larger scales. The fractal dimension has no significant dependence on column density and the values D = 1.25 +- 0.07 are only slightly lower than in typical molecular clouds. The column density PDFs exhibit large variations, e.g. in the case of externally compressed clouds. At scales r>0.1 pc, the radial column density distributions of the clouds follow an average relation of N~r^{-1}. In spite of a great variety of clump morphology, clumps tend to follow a similar N~r^{-1} relation below r~0.1 pc. RT calculations indicate only factor of 2.5 variation in the local radiation field intensity. The fraction of gravitationally bound clumps increases significantly in regions with A_V > 5 mag but most bound objects appear to be pressure-confined. The GCC host clouds have statistical properties similar to general molecular clouds. The gravitational stability, peak column density, and clump orientation are connected to the cloud background while most other statistics (e.g. D and radial profiles) are insensitive to the environment.

Analysis of all-sky Planck submillimetre observations and the IRAS 100um data has led to the detection of a population of Galactic cold clumps. The clumps can be used to study star formation and dust properties in a wide range of Galactic environments. Our aim is to measure dust spectral energy distribution (SED) variations as a function of the spatial scale and the wavelength. We examine the SEDs at large scales using IRAS, Planck, and Herschel data. At smaller scales, we compare with JCMT/SCUBA-2 850um maps with Herschel data that are filtered using the SCUBA-2 pipeline. Clumps are extracted using the Fellwalker method and their spectra are modelled as modified blackbody functions. According to IRAS and Planck data, most fields have dust colour temperatures T_C ~ 14-18K and opacity spectral index values of beta=1.5-1.9. The clumps/cores identified in SCUBA-2 maps have T~ 13K and similar beta values. There are some indications of the dust emission spectrum becoming flatter at wavelengths longer than 500um. In fits involving Planck data, the significance is limited by the uncertainty of the corrections for CO line contamination. The fits to the SPIRE data give a median beta value slightly above 1.8. In the joint SPIRE and SCUBA-2 850um fits the value decreases to beta ~1.6. Most of the observed T-beta anticorrelation can be explained by noise. The typical submillimetre opacity spectral index beta of cold clumps is found to be ~1.7. This is above the values of diffuse clouds but lower than in some previous studies of dense clumps. There is only tentative evidence of T-beta anticorrelation and beta decreasing at millimetre wavelengths.

The combination of line and continuum observations can provide vital insights to the formation and fragmentation of filaments and the initial conditions for star formation. We have carried out line observations to map the kinematics of an evolved, actively star forming filament G82.65-2.00. The filament was first The combination of line and continuum observations can provide vital insights to the formation and fragmentation of filaments and the initial conditions for star formation. We have carried out line observations to map the kinematics of an evolved, actively star forming filament G82.65-2.00. The filament was first identified from the Planck data as a region of particularly cold dust emission and was mapped at 100-500 $\mu$m as a part of Herschel key program Galactic Cold Cores. The observations reveal several velocity components in the field, with strongest line emission concentrated to velocity range $\sim [3,5]$ km $\rm s^{-1}$. The column density of molecular hydrogen along the filament varies from 1.0 to 2.3 $\times 10^{22}$ $\rm cm^{2}$. We have examined six cold clumps from the central part of the filament. The clumps have masses in the range $10 - 20$ $M_{\odot}$ ($\sim70M_{\odot}$ in total) and are close to or above the virial mass. Furthermore, the main filament is heavily fragmented and most of the the substructures have a mass lower than or close to the virial mass, suggesting that the filament is dispersing as a whole. Position-velocity maps of $\rm ^{12}CO$ and $\rm ^{13}CO$ lines indicate that at least one of the striations is kinematically connected to two of the clumps, potentially indicating mass accretion from the striation onto the main filament. We tentatively estimate the accretion rate to be $\dot{M}$ = $2.23 \times 10^{-6}$ $ M_{\odot} / \rm year$.

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

The nearby cloud L1642 is one of only two known very high latitude (|b| > 30 deg) clouds actively forming stars. It is a rare example of star formation in isolated conditions, and can reveal important details of star formation in general, e.g., of the effect of magnetic fields. We compare Herschel dust emission structures and magnetic field orientation revealed by Planck polarization maps in L1642. The high-resolution ($\sim20"$) Herschel data reveal a complex structure including a dense, compressed central clump, and low density striations. The Planck polarization data (at 10$'$ resolution) reveal an ordered magnetic field pervading the cloud and aligned with the surrounding striations. There is a complex interplay between the cloud structure and large scale magnetic field. This suggests that the magnetic field is closely linked to the formation and evolution of the cloud. CO rotational emission confirms that the striations are connected with the main clumps and likely to contain material either falling into or flowing out of the clumps. There is a clear transition from aligned to perpendicular structures approximately at a column density of $N_{\rm{H}} = 1.6 \times 10^{21}\, {{\rm cm}}^{-2}$. Comparing the Herschel maps with the Planck polarization maps shows the close connection between the magnetic field and cloud structure even in the finest details of the cloud.

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

We observed high S/N, high velocity resolution NH$_3$(1,1) and (2,2) emission on an extended map in TMC-1, a filamentary cloud in a nearby quiescent star forming area. By fitting multiple hyperfine-split line profiles to the NH$_3$(1,1) spectra we derived the velocity distribution of the line components and calculated gas parameters on several positions. Herschel SPIRE continuum observations were reduced and used to calculate the physical parameters of the Planck Galactic Cold Clumps in the region. The Herschel-based column density map of TMC-1 shows a main ridge with two local maxima and a separated peak to the south-west. H$_2$-column densities and dust temperatures are in the range of 0.5-3.3 $\times$ 10$^{22}$ cm$^{-2}$ and 10.5-12 K, respectively. NH$_3$-column densities are 2.8-14.2 $\times$ 10$^{14}$ cm$^{-2}$ and and H$_2$-volume densities are 0.4-2.8 $\times$ 10$^4$ cm$^{-3}$. Kinetic temperatures are typically very low with a minimum of 9 K, and a maximum of 13.7 K was found at the Class I protostar IRAS 04381+2540. The kinetic temperatures vary similarly as the dust temperatures in spite of the fact that densities are lower than the critical density for coupling between the gas and dust phase. The k-means clustering method separated four sub-filaments in TMC-1 in the position-velocity-column density parameter space. They have masses of 32.5, 19.6, 28.9 and 45.9 M$_{\odot}$, low turbulent velocity dispersion (0.13-0.2 kms$^{-1}$) and they are close to gravitational equilibrium. We label them TMC-1F1 through F4. TMC-1F1, TMC-1F2 and TMC-1F4 are very elongated, dense and cold. TMC-1F3 is a little less elongated and somewhat warmer, probably heated by IRAS 04381+2540 that is embedded in it. TMC-1F3 is $\approx$ 0.1 pc behind TMC1-F1. Because of its structure, TMC-1 is a good target to test filament evolution scenarios.

The Galactic Cold Cores project has made Herschel observations of 116 fields where the Planck survey has found signs of cold dust emission. The fields contain sources in different environments and different phases of star formation. The dust opacity spectral index beta and the dust colour temperature T are derived using Herschel and Planck data. The relation between beta and T is examined for the whole sample and inside individual fields. Based on IRAS and Planck data, the fields are characterised by a median colour temperature of 16.1 K and a median opacity spectral index of beta=1.84. We observe a clear T-beta anti-correlation. In Herschel observations, constrained at lower resolution by Planck data, the variations follow the column density structure and beta(FIR) can rise to ~2.2 in individual clumps. The Planck 217 GHz band shows a systematic excess that is consistent with a general flattening of the dust emission spectrum at millimetre wavelengths. When fitted separately below and above 700 um, the median spectral index values are beta(FIR) ~ 1.91 and beta(mm) ~ 1.66. The spectral index changes as a function of column density and wavelength. Beta variations are partly masked by temperature gradients and the changes in the intrinsic grain properties may be even greater.

Recently discovered scattered light from molecular cloud cores in the wavelength range 3-5 {\mu}m (called "coreshine") seems to indicate the presence of grains with sizes above 0.5 {\mu}m. We aim to analyze 3.6 and 4.5 {\mu}m coreshine from molecular cloud cores to probe the largest grains in the size distribution. We analyzed dedicated deep Cycle 9 Spitzer IRAC observations in the 3.6 and 4.5 {\mu}m bands for a sample of 10 low-mass cores. We used a new modeling approach based on a combination of ratios of the two background- and foreground-subtracted surface brightnesses and observed limits of the optical depth. The dust grains were modeled as ice-coated silicate and carbonaceous spheres. We discuss the impact of local radiation fields with a spectral slope differing from what is seen in the DIRBE allsky maps. For the cores L260, ecc806, L1262, L1517A, L1512, and L1544, the model reproduces the data with maximum grain sizes around 0.9, 0.5, 0.65, 1.5, 0.6, and > 1.5 {\mu}m, respectively. The maximum coreshine intensities of L1506C, L1439, and L1498 in the individual bands require smaller maximum grain sizes than derived from the observed distribution of band ratios. Additional isotropic local radiation fields with a spectral shape differing from the DIRBE map shape do not remove this discrepancy. In the case of Rho Oph 9, we were unable to reliably disentangle the coreshine emission from background variations and the strong local PAH emission. Considering surface brightness ratios in the 3.6 and 4.5 {\mu}m bands across a molecular cloud core is an effective method of disentangling the complex interplay of structure and opacities when used in combination with observed limits of the optical depth.

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

Context. he study of prestellar cores is essential to understanding the initial stages of star formation. With $Herschel$ more cold clumps have been detected than ever before. For this study we have selected 21 cold clumps from 20 $Herschel$ fields observed as a follow-up on original $Planck$ detections. We have observed these clumps in $^{13}$CO (1-0), C$^{18}$O (1-0), and N$_2$H$^+$ (1-0) lines. Aims. Our aim is to find out if these cold clumps are prestellar. We have examined to what extent independent analysis of the dust and the molecular lines lead to similar conclusions about the masses of these objects. Methods. We calculate the clump masses and densities from the dust continuum and molecular line observations and compare these to each other and to the virial and Bonnor-Ebert masses calculated for each clump. Finally we examine two of the fields with radiative transfer models to estimate CO abundances. Results. When excitation temperatures could be estimated, the column densities derived from molecular line observations were comparable to those from dust continuum data. The median column density estimates are 4.2$\times 10^{21}$cm$^{-2}$ and 5.5$\times 10^{21}$cm$^{-2}$ for the line and dust emission data, respectively. The calculated abundances, column densities, volume densities, and masses all have large uncertainties and one must be careful when drawing conclusions. Abundance of $^{13}$CO was found in modeling the two clumps in the field G131.65$+$9.75 to be close to the usual value of 10$^{-6}$. The abundance ratio of $^{13}$CO and C$^{18}$O was $\sim$10. Molecular abundances could only be estimated with modeling, relying on dust column density data. Conclusions. The results indicate that most cold clumps, even those with dust color temperatures close to 11 K, are not necessarily prestellar.

We present the Planck Catalogue of Galactic Cold Clumps (PGCC), an all-sky catalogue of Galactic cold clump candidates detected by Planck. This catalogue is the full version of the Early Cold Core (ECC) catalogue, which was made available in 2011 with the Early Release Compact Source Catalogue (ERCSC) and contained 915 high S/N sources. It is based on the Planck 48 months mission data that are currently being released to the astronomical community. The PGCC catalogue is an observational catalogue consisting exclusively of Galactic cold sources. The three highest Planck bands (857, 545, 353 GHz) have been combined with IRAS data at 3 THz to perform a multi-frequency detection of sources colder than their local environment. After rejection of possible extragalactic contaminants, the PGCC catalogue contains 13188 Galactic sources spread across the whole sky, i.e., from the Galactic plane to high latitudes, following the spatial distribution of the main molecular cloud complexes. The median temperature of PGCC sources lies between 13 and 14.5 K, depending on the quality of the flux density measurements, with a temperature ranging from 5.8 to 20 K after removing sources with the 1% largest temperature estimates. Using seven independent methods, reliable distance estimates have been obtained for 5574 sources, which allows us to derive their physical properties such as their mass, physical size, mean density and luminosity. The PGCC sources are located mainly in the solar neighbourhood, up to a distance of 10.5 kpc towards the Galactic centre, and range from low-mass cores to large molecular clouds. Because of this diversity and because the PGCC catalogue contains sources in very different environments, the catalogue is useful to investigate the evolution from molecular clouds to cores. Finally, the catalogue also includes 54 additional sources located in the SMC and LMC.

The project Galactic Cold Cores has made Herschel observations of interstellar clouds where the Planck satellite survey has located cold and compact clumps. The sources range from starless clumps to protostellar cores. We examine 116 Herschel fields to estimate the submillimetre dust opacity and its variations. The submillimetre dust opacity was derived from Herschel data, and near-infrared observations of the reddening of background stars are converted into near-infrared optical depth. We studied the systematic errors affecting these parameters and used modelling to correct for the expected biases. The ratio of 250um and J band opacities is correlated with the cloud location and star formation activity. We find a median ratio of tau(250um)/tau(J)= (1.6+-0.2)*10^-3, which is more than three times the mean value in diffuse medium. Assuming a spectral index beta=1.8 instead of beta=2.0, the value would be lower by ~30%. No significant systematic variation is detected with Galactocentric distance or with Galactic height. The tau(250um)/tau(J) maps reveal six fields with clear increase of submillimetre opacity of up to tau(250um)/tau(J) ~ 4*10^-3. These are all nearby fields with spatially resolved clumps of high column density. We interpret the increase in the far-infrared opacity as a sign of grain growth in the densest and coldest regions of interstellar clouds.

Context. Dust is a good tracer of cold dark clouds but its column density is difficult to quantify. Aims. We want to check whether the far-infrared and submillimeter high-resolution data from Herschel SPIRE and PACS cameras combined with ground-based telescope bolometers allow us to retrieve the whole dust content of cold dark clouds. Methods. We compare far-infrared and submillimeter emission across L183 to the 8 $\mu$m absorption map from Spitzer data and fit modified blackbody functions towards three different positions. Results. We find that none of the Herschel SPIRE channels follow the cold dust profile seen in absorption. Even the ground-based submillimeter telescope observations, although more closely following the absorption profile, cannot help to characterize the cold dust without external information such as the dust column density itself. The difference in dust opacity can reach up to a factor of 3 in prestellar cores of high extinction. Conclusions. In dark clouds, the amount of very cold dust cannot be measured from its emission alone. In particular, studies of dark clouds based only on Herschel data can miss a large fraction of the dust content. This has an impact on core and filament density profiles, masse and stability estimates.

Polarized emission observed by Planck HFI at 353 GHz towards a sample of nearby fields is presented, focusing on the statistics of polarization fractions $p$ and angles $\psi$. The polarization fractions and column densities in these nearby fields are representative of the range of values obtained over the whole sky. We find that: (i) the largest polarization fractions are reached in the most diffuse fields; (ii) the maximum polarization fraction $p_\mathrm{max}$ decreases with column density $N_\mathrm{H}$ in the more opaque fields with $N_\mathrm{H} > 10^{21}\,\mathrm{cm}^{-2}$; and (iii) the polarization fraction along a given line of sight is correlated with the local spatial coherence of the polarization angle. These observations are compared to polarized emission maps computed in simulations of anisotropic magnetohydrodynamical (MHD) turbulence in which we assume a uniform intrinsic polarization fraction of the dust grains. We find that an estimate of this parameter may be recovered from the maximum polarization fraction $p_\mathrm{max}$ in diffuse regions where the magnetic field is ordered on large scales and perpendicular to the line of sight. This emphasizes the impact of anisotropies of the magnetic field on the emerging polarization signal. The decrease of the polarization fraction with column density in nearby molecular clouds is well reproduced in the simulations, indicating that it is essentially due to the turbulent structure of the magnetic field: an accumulation of variously polarized structures along the line of sight leads to such an anti-correlation. In the simulations, polarization fractions are also found to anti-correlate with the angle dispersion function $\mathcal{S}$. [abridged]

Context. Using observations to deduce dust properties, grain size distribution, and physical conditions in molecular clouds is a highly degenerate problem. Aims. The coreshine phenomenon, a scattering process at 3.6 and 4.5 $\mu$m that dominates absorption, has revealed its ability to explore the densest parts of clouds. We want to use this effect to constrain the dust parameters. The goal is to investigate to what extent grain growth (at constant dust mass) inside molecular clouds is able to explain the coreshine observations. We aim to find dust models that can explain a sample of Spitzer coreshine data. We also look at the consistency with near-infrared data we obtained for a few clouds. Methods. We selected four regions with a very high occurrence of coreshine cases: Taurus-Perseus, Cepheus, Chameleon and L183/L134. We built a grid of dust models and investigated the key parameters to reproduce the general trend of surface bright- nesses and intensity ratios of both coreshine and near-infrared observations with the help of a 3D Monte-Carlo radiative transfer code. The grid parameters allow to investigate the effect of coagulation upon spherical grains up to 5 $\mu$m in size derived from the DustEm diffuse interstellar medium grains. Fluffiness (porosity or fractal degree), ices, and a handful of classical grain size distributions were also tested. We used the near- and mostly mid-infrared intensity ratios as strong discriminants between dust models. Results. The determination of the background field intensity at each wavelength is a key issue. In particular, an especially strong background field explains why we do not see coreshine in the Galactic plane at 3.6 and 4.5 $\mu$m. For starless cores, where detected, the observed 4.5 $\mu$m / 3.6 $\mu$m coreshine intensity ratio is always lower than $\sim$0.5 which is also what we find in the models for the Taurus-Perseus and L183 directions. Embedded sources can lead to higher fluxes (up to four times greater than the strongest starless core fluxes) and higher coreshine ratios (from 0.5 to 1.1 in our selected sample). Normal interstellar radiation field conditions are sufficient to find suitable grain models at all wavelengths for starless cores. The standard interstellar grains are not able to reproduce observations and, due to the multi-wavelength approach, only a few grain types meet the criteria set by the data. Porosity does not affect the flux ratios while the fractal dimension helps to explain coreshine ratios but does not seem able to reproduce near-infrared observations without a mix of other grain types. Conclusions. Combined near- and mid-infrared wavelengths confirm the potential to reveal the nature and size distribution of dust grains. Careful assessment of the environmental parameters (interstellar and background fields, embedded or nearby reddened sources) is required to validate this new diagnostic.

L1642 is one of the two high galactic latitude (|b| > 30deg) clouds confirmed to have active star formation. We examine the properties of this cloud, especially the large-scale structure, dust properties, and compact sources in different stages of star formation. We present high-resolution far-infrared and submm observations with the Herschel and AKARI satellites and mm observations with the AzTEC/ASTE telescope, which we combined with archive data from near- and mid-infrared (2MASS, WISE) to mm observations (Planck). The Herschel observations, combined with other data, show a sequence of objects from a cold clump to young stellar objects at different evolutionary stages. Source B-3 (2MASS J04351455-1414468) appears to be a YSO forming inside the L1642 cloud, instead of a foreground brown dwarf, as previously classified. Herschel data reveal striation in the diffuse dust emission around L1642. The western region shows striation towards NE and has a steeper column density gradient on its southern side. The densest central region has a bow-shock like structure showing compression from the west and a filamentary tail extending towards east. The differences suggest that these may be spatially distinct structures, aligned only in projection. We derive values of the dust emission cross-section per H nucleon for different regions of the cloud. Modified black-body fits to the spectral energy distribution of Herschel and Planck data give emissivity spectral index beta values 1.8-2.0 for the different regions. The compact sources have lower beta values and show an anticorrelation between T and beta. Markov chain Monte Carlo calculations demonstrate the strong anticorrelation between beta and T errors and the importance of mm Planck data in constraining the estimates. L1642 reveals a more complex structure and sequence of star formation than previously known.

Mapping of near-infrared (NIR) scattered light is a recent method for the study of interstellar clouds, complementing other, more commonly used methods, like dust emission and extinction. Our goal is to study the usability of this method on larger scale, and compare the properties of a filament using NIR scattering and other methods. We also study the radiation field and differences in grain emissivity between diffuse and dense areas. We have used scattered J, H, and K band surface brightness WFCAM-observations to map filament TMC-1N in Taurus, covering an area of 1dx1d corresponding to ~(2.44 pc)^2. We have converted the data into optical depth and compared the results with NIR extinction and Herschel observations of submm dust emission. We see the filament in scattered light in all three NIR bands. We note that our WFCAM observations in TMC-1N show notably lower intensity than previous results in Corona Australis using the same method. We show that 3D radiative transfer simulations predict similar scattered surface brightness levels as seen in the observations. However, changing the assumptions about the background can change the results of simulations notably. We derive emissivity by using optical depth in the J band as an independent tracer of column density. We obtain opacity sigma(250um) values 1.7-2.4x10^-25 cm^2/H, depending on assumptions of the extinction curve, which can change the results by over 40%. These values are twice as high as obtained for diffuse areas, at the lower limit of earlier results for denser areas. We show that NIR scattering can be a valuable tool in making high resolution maps. We conclude, however, that NIR scattering observations can be complicated, as the data can show relatively low-level artefacts. This suggests caution when planning and interpreting the observations.

With recent Herschel observations, the northern filament of the Corona Australis cloud has now been mapped in a number of bands from 1.2um to 870um. The data set provides a good starting point for the study of the cloud over several orders of magnitude in density. We wish to examine the differences of the column density distributions derived from dust extinction, scattering, and emission, and to determine to what extent the observations are consistent with the standard dust models. From Herschel data, we calculate the column density distribution that is compared to the corresponding data derived in the near-infrared regime from the reddening of the background stars, and from the surface brightness attributed to light scattering. We construct three-dimensional radiative transfer models to describe the emission and the scattering. The scattered light traces low column densities of A_V~1mag better than the dust emission, remaining useful to A_V ~ 10-15 mag. Based on the models, the extinction and the level of dust emission are surprisingly consistent with a sub-millimetre dust emissivity typical of diffuse medium. However, the intensity of the scattered light is very low at the centre of the densest clump and this cannot be explained without a very low grain albedo. Both the scattered light and dust emission indicate an anisotropic radiation field. The modelling of the dust emission suggests that the radiation field intensity is at least three times the value of the normal interstellar radiation field. The inter-comparison between the extinction, light scattering, and dust emission provides very stringent constraints on the cloud structure, the illuminating radiation field, and the grain properties.

We examine the cloud structure around the Planck detections in 71 fields observed with the Herschel SPIRE instrument. We wish to determine the general physical characteristics of the fields and to examine the morphology of the clouds where the cold high column density clumps are found. We derive colour temperature and column density maps of the fields. We examine the infrared spectral energy distributions of the main clumps. The clouds are categorised according to their large scale morphology. With the help of recently released WISE satellite data, we look for signs of enhanced mid-infrared scattering (coreshine), an indication of growth of the dust grains, and examine the star formation activity associated with the cold clumps. The mapped clouds have distances ranging from ~100pc to several kiloparsecs and cover a range of sizes and masses from cores of less than 10 solar masses to clouds with masses in excess of 10000 solar mass. Most fields contain some filamentary structures and in about half of the cases a filament or a few filaments dominate the morphology. In one case out of ten, the clouds show a cometary shape or have sharp boundaries indicative of compression by an external force. The width of the filaments is typically ~0.2-0.3pc. However, there is significant variation from 0.1pc to 1pc and the estimates are sensitive to the methods used and the very definition of a filament. Enhanced mid-infrared scattering, coreshine, was detected in four clouds with six additional tentative detections. The cloud LDN183 is included in our sample and remains the best example of this phenomenon. About half of the fields are associated with active star formation as indicated by the presence of mid-infrared point sources. The mid-infrared sources often coincide with structures whose sub-millimetre spectra are still dominated by the cold dust.

Within the project Galactic Cold Cores we are carrying out Herschel photometric observations of cold interstellar clouds detected with the Planck satellite. The three fields observed as part of the Herschel science demonstration phase (SDP) provided the first glimpse into the nature of these sources. We examine the properties of the dust emission within the fields. We determine the dust sub-millimetre opacity, look for signs of spatial variations in the dust spectral index, and estimate how the apparent variations of the parameters could be affected by different sources of uncertainty. We use the Herschel observations where the zero point of the surface brightness scale is set with the help of the Planck satellite data. We derive the colour temperature and column density maps of the regions and determine the dust opacity by a comparison with extinction measurements. By simultaneously fitting the colour temperature and the dust spectral index values we look for spatial variations in the apparent dust properties. With a simple radiative transfer model we estimate to what extent these can be explained by line-of-sight temperature variations, without changes in the dust grain properties. The analysis of the dust emission reveals cold and dense clouds that coincide with the Planck sources and confirm those detections. The derived dust opacity varies in the range kappa(250um) ~ 0.05-0.2 cm^2/g, higher values being observed preferentially in regions of high column density. The average dust spectral index beta is ~ 1.9-2.2. There are indications that beta increases towards the coldest regions. The spectral index decreases strongly near internal heating sources but, according to radiative transfer models, this can be explained by the line-of-sight temperature variations without a change in the dust properties.

We present the statistical properties of the first version of the Cold Core Catalogue of Planck Objects (C3PO), in terms of their spatial distribution, temperature, distance, mass, and morphology. We also describe the statistics of the Early Cold Core Catalogue (ECC, delivered with the Early Release Compact Source Catalogue, ERCSC) that is the subset of the 915 most reliable detections of the complete catalogue. We have used the CoCoCoDeT algorithm to extract 10783 cold sources. Temperature and dust emission spectral index {\beta} values are derived using the fluxes in the IRAS 100 \mum band and the three highest frequency Planck bands. Temperature spans from 7K to 17K, and peaks around 13K. Data are not consistent with a constant value of {\beta} over the all temperature range. {\beta} ranges from 1.4 to 2.8 with a mean value around 2.1, and several possible scenarios are possible, including {\beta}(T) and the effect of multiple T components folded into the measurements. For one third of the objects the distances are obtained. Most of the detections are within 2 kpc in the Solar neighbourhood, but a few are at distances greater than 4 kpc. The cores are distributed from the deep Galactic plane, despite the confusion, to high latitudes (>30$^{\circle}$). The associated mass estimates range from 1 to $10^5$ solar masses. Using their physical properties these cold sources are shown to be cold clumps, defined as the intermediate cold sub-structures between clouds and cores. These cold clumps are not isolated but mostly organized in filaments associated with molecular clouds. The Cold Core Catalogue of Planck Objects (C3PO) is the first unbiased all-sky catalogue of cold objects. It gives an unprecedented statistical view to the properties of these potential pre-stellar clumps and offers a unique possibility for their classification in terms of their intrinsic properties and environment.

(abridged) We perform a detailed investigation of sources from the Cold Cores Catalogue of Planck Objects (C3PO). Our goal is to probe the reliability of the detections, validate the separation between warm and cold dust emission components, provide the first glimpse at the nature, internal morphology and physical characterictics of the Planck-detected sources. We focus on a sub-sample of ten sources from the C3PO list, selected to sample different environments, from high latitude cirrus to nearby (150pc) and remote (2kpc) molecular complexes. We present Planck surface brightness maps and derive the dust temperature, emissivity spectral index, and column densities of the fields. With the help of higher resolution Herschel and AKARI continuum observations and molecular line data, we investigate the morphology of the sources and the properties of the substructures at scales below the Planck beam size.

We study a northern part of the Corona Australis molecular cloud that consists of a filament and a dense sub-millimetre core inside the filament. Our aim is to measure dust temperature and sub-mm emissivity within the region. We also look for confirmation that near-infrared (NIR) surface brightness can be used to study the structure of even very dense clouds. We extend our previous NIR mapping south of the filament. The dust colour temperatures are estimated using Spitzer 160um and APEX/Laboca 870um maps. The column densities derived based on the reddening of background stars, NIR surface brightness, and thermal sub-mm dust emission are compared. A three dimensional toy model of the filament is used to study the effect of anisotropic illumination on near-infrared surface brightness and the reliability of dust temperature determination. Relative to visual extinction, the estimated emissivity at 870um is kappa(870) = (1.3 +- 0.4) x 10^{-5} 1/mag. This is similar to the values found in diffuse medium. A significant increase in the sub-millimetre emissivity seems to be excluded. In spite of saturation, NIR surface brightness was able to accurately pinpoint, and better than measurements of the colour excesses of background stars, the exact location of the column density maximum. Both near- and far-infrared data show that the intensity of the radiation field is higher south of the filament.

Polarization carries information about the magnetic fields in interstellar clouds. The observations of polarized dust emission are used to study the role of magnetic fields in the evolution of molecular clouds and the initial phases of star-formation. We study the grain alignment with realistic simulations, assuming the radiative torques to be the main mechanism that spins the grains up. The aim is to study the efficiency of the grain alignment as a function of cloud position and to study the observable consequences of these spatial variations. Our results are based on the analysis of model clouds derived from MHD simulations. The continuum radiative transfer problem is solved with Monte Carlo methods to estimate the 3D distribution of dust emission and the radiation field strength affecting the grain alignment. We also examine the effect of grain growth in cores. We are able to reproduce the results of Cho & Lazarian using their assumptions. However, the anisotropy factor even in the 1D case is lower than their assumption of $\gamma = 0.7$, and thus we get less efficient radiative torques. Compared with our previous paper, the polarization degree vs. intensity relation is steeper because of less efficient grain alignment within dense cores. Without grain growth, the magnetic field of the cores is poorly recovered above a few $A_{\rm V}$. If grain size is doubled in the cores, the polarization of dust emission can trace the magnetic field lines possibly up to $A_{\rm V} \sim 10$ magnitudes. However, many of the prestellar cores may be too young for grain coagulation to play a major role. The inclusion of direction dependent radiative torque efficiency weakens the alignment. Even with doubled grain size, we would not expect to probe the magnetic field past a few magnitudes in $A_{\rm V}$.

With current near-infrared (NIR) instruments the near-infrared light scattered from interstellar clouds can be mapped over large areas. The surface brightness carries information on the line-of-sight dust column density. Therefore, scattered light could provide an important tool to study mass distribution in quiescent interstellar clouds at a high, even sub-arcsecond resolution. We wish to confirm the assumption that light scattering dominates the surface brightness in all NIR bands. Furthermore, we want to show that scattered light can be used for an accurate estimation of dust column densities in clouds with Av in the range 1-15mag. We have obtained NIR images of a quiescent filament in the Corona Australis molecular cloud. The observations provide maps of diffuse surface brightness in J, H, and Ks bands. Using the assumption that signal is caused by scattered light we convert surface brightness data into a map of dust column density. The same observations provide colour excesses for a large number of background stars. These data are used to derive an extinction map of the cloud. The two, largely independent tracers of the cloud structure are compared. Results. In regions below Av=15m both diffuse surface brightness and background stars lead to similar column density estimates. The existing differences can be explained as a result of normal observational errors and bias in the sampling of extinctions provided by the background stars. There is no indication that thermal dust emission would have a significant contribution even in the Ks band. The results show that, below Av=15mag, scattered light does provide a reliable way to map cloud structure. Compared with the use of background stars it can also in practice provide a significantly higher spatial resolution.

Our aim is to study the polarization of thermal dust emission to see if the alignment of grain by radiative torques could explain the observed relation between the degree of polarization and the intensity in dense cores. Predictions are made for polarimetry observations with the Planck satellite. Our results are based on model clouds derived from MHD simulations of magnetized turbulent flows, while the continuum radiative transfer problem is solved with Monte Carlo methods in order to estimate the three-dimensional distribution of dust emission and the radiation field strength affecting the grain alignment. The influence of grain alignment efficiency is examined in the calculated polarization maps. We are able to reproduce the P/I-relation with the grain alignment by radiative torques. The decrease in intrinsic polarization and total emission means that sub-mm polarimetry carries only little information about the magnetic fields in dense cores with high visual extinction. The interpretation of the observations will be further complicated by the unknown magnetic field geometry and the fact that what is observed as individual cores may, in fact, be a superposition of several density enhancements. According to our calculations, Planck will be able to map dust polarization reliably when A_V > 2 mag at spatial resolution of 15'.

With current wide-field near-infrared (NIR) instruments the scattered light in the near-infrared can be mapped over large areas. Below A_V ~ 10 the surface brightness is directly proportional to the column density, and at slightly higher column densities the saturation of the intensity values can be corrected using the ratios of the intensity in different NIR bands. NIR scattered light provides a promising new method for the mapping of quiescent interstellar clouds. We develop a method to convert the observed near-infrared surface brightness into estimates of the column density. We study and quantify the effect that different error sources could have on the accuracy of such estimates. We also propose to reduce systematic errors by combining surface brightness data with extinction measurements derived from the near-infrared colour excess of background stars. Our study is based on a set of three-dimensional magnetohydrodynamic turbulence simulations. Maps of near-infrared scattered light are obtained with radiative transfer calculations, and the maps are converted back into column density estimates using the proposed method. The results are compared with the true column densities. Extinction measurements are simulated using the same turbulence simulations, and are used as a complementary column density tracer. We find that NIR intensities can be converted into a reliable estimate of the column density in regions with A_V up to almost 20mag.

We discuss the possibility of mapping interstellar clouds at unprecedentedly high spatial resolution by means of near-IR imaging of their scattered light. We calculate the scattering of the interstellar radiation field by a cloud model obtained from the simulation of a supersonic turbulent flow. Synthetic maps of scattered light are computed in the J, H and K bands and are found to allow an accurate estimate of column density, in the range of visual extinction between 1 and 20 magnitudes. We provide a formalism to convert the intensity of scattered light at these near-IR bands into a total gas column density. We also show that this new method of mapping interstellar clouds is within the capability of existing near-IR facilities, which can achieve a spatial resolution of up to ~ 0.1 arcsec. This opens new perspectives in the study of interstellar dust and gas structure on very small scales. The validity of the method has been recently demonstrated by the extraordinary images of the Perseus region obtained by Foster & Goodman (2005), which motivated this investigation.