Mapping of interstellar clouds with infrared light scattered from dust: TMC-1N

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.Comment: abstract shortened and figures reduced for astrop

A Corona Australis cloud filament seen in NIR scattered light. III. Modelling and comparison with Herschel sub-millimetre data

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.Comment: 13 pages, 16 figures, accepted to A&

Evaporation Ages: a New Dating Method for Young Star Clusters

The ages of young star clusters are fundamental clocks to constrain the formation and evolution of pre-main-sequence stars and their protoplanetary disks and exoplanets. However, dating methods for very young clusters often disagree, casting doubts on the accuracy of the derived ages. We propose a new method to derive the kinematic age of star clusters based on the evaporation ages of their stars. The method is validated and calibrated using hundreds of clusters identified in a supernova-driven simulation of the interstellar medium forming stars for approximately 40 Myr within a 250 pc region. We demonstrate that the clusters' evaporation-age uncertainty can be as small as about 10% for clusters with a large enough number of evaporated stars and small but realistic observational errors. We have obtained evaporation ages for a pilot sample of 10 clusters, finding a good agreement with their published isochronal ages. The evaporation ages will provide important constraints for modeling the pre-main-sequence evolution of low-mass stars, as well as to investigate the star-formation and gas-evaporation history of young clusters. These ages can be more accurate than isochronal ages for very young clusters, for which observations and models are more uncertain.Comment: 13 pages, 11 figures, 2 tables, submitted to A&A on Nov 14th, 202

From the CMF to the IMF: Beyond the Core-Collapse Model

Observations have indicated that the prestellar core mass function (CMF) is similar to the stellar initial mass function (IMF), except for an offset towards larger masses. This has led to the idea that there is a one-to-one relation between cores and stars, such that the whole stellar mass reservoir is contained in a gravitationally-bound prestellar core, as postulated by the core-collapse model, and assumed in recent theoretical models of the stellar IMF. We test the validity of this assumption by comparing the final mass of stars with the mass of their progenitor cores in a high-resolution star-formation simulation that generates a realistic IMF under physical conditions characteristic of observed molecular clouds. Using a definition of bound cores similar to previous works we obtain a CMF that converges with increasing numerical resolution. We find that the CMF and the IMF are closely related in a statistical sense only; for any individual star there is only a weak correlation between the progenitor core mass and the final stellar mass. In particular, for high mass stars only a small fraction of the final stellar mass comes from the progenitor core, and even for low mass stars the fraction is highly variable, with a median fraction of only about 50%. We conclude that the core-collapse scenario and related models for the origin of the IMF are incomplete. We also show that competitive accretion is not a viable alternative.Comment: 23 pages, 29 figures. Link to supplementary material and full Table 1: http://www.erda.dk/vgrid/core-mass-function/ . Submitted to MNRA

Multiwavelength study of the high-latitude cloud L1642: chain of star formation

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.Comment: 22 pages, 18 figures, accepted to Astronomy & Astrophysics; abstract shortened and figures reduced for astrop

Grain size limits derived from 3.6 {\mu}m and 4.5 {\mu}m coreshine

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.Comment: 23 pages, 18 figures, accepted for publication in A&

Dust properties inside molecular clouds from coreshine modeling and observations

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

A Corona Australis cloud filament seen in NIR scattered light II: Comparison with sub-millimeter data

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.Comment: 9 pages, 9 figures, accepted to A&

Galactic cold cores: Herschel study of first Planck detections

Context. We present the first results from the project Galactic cold cores, where the cold interstellar clouds detected by the Planck satellite are studied with Herschel photometric observations. The final Planck catalogue is expected to contain several thousand sources. The Herschel observations during the science demonstration phase provided the first glimpse into the nature of these sources. Aims. The main goal of the project is to derive the physical properties of the cold core population revealed by Planck. We examine three fields and confirm the Planck detections with Herschel data, which we also use to establish the evolutionary stage of the identified cores. Methods. We study the morphology and spectral energy distribution of the sources using the combined wavelength coverage of Planck and Herschel. The dust colour temperatures and emissivity indices are determined. The masses of the cores are determined with distance estimates which are taken from the literature and are confirmed by kinematic and extinction information. Results. The observations reveal extended regions of cold dust with dust colour temperatures down to T_(dust) ~ 11 K. The fields represent different evolutionary stages ranging from a quiescent, cold filament inMusca to regions of active star formation in Cepheus. Conclusions. The Herschel observations confirm that the all-sky survey of Planck is capable of making a large number of new cold core detections. Our results suggest that many of the sources may already have left the pre-stellar phase or are at least closely associated with active star formation. High-resolution Herschel observations are needed to establish the true nature of the Planck detections

Galactic cold cores VIII. Filament formation and evolution : Filament properties in context with evolutionary models

Context. The onset of star formation is intimately linked with the presence of massive unstable filamentary structures. These filaments are therefore key for theoretical models that aim to reproduce the observed characteristics of the star formation process in the Galaxy. Aims. As part of the filament study carried out by the Herschel Galactic Cold Cores Key Programme, here we study and discuss the filament properties presented in GCC VII (Paper I) in context with theoretical models of filament formation and evolution. Methods. A conservatively selected sample of filaments located at a distance D <500 pc was extracted from the GCC fields with the getfilaments algorithm. The physical structure of the filaments was quantified according to two main components: the central (Gaussian) region of the filament (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 the column density of its environment were compared with the predictions from theoretical models describing the evolution of filaments under gravity-dominated conditions. Results. The feasibility of a transition from a subcritical to supercritical state by accretion at any given time is dependent on the combined effect of filament intrinsic properties and environmental conditions. Reasonably self-gravitating (high M-line,M-core) filaments in dense environments (Av greater than or similar to 3 mag) can become supercritical on timescales of t similar to 1 Myr by accreting mass at constant or decreasing width. The trend of increasing M-line,M-tot (M-line,M-core and M-line,M-wing) and ridge A(v) with background for the filament population also indicates that the precursors of star-forming filaments evolve coevally with their environment. The simultaneous increase of environment and filament Av explains the observed association between dense environments and high Mlille,co values, and it 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 filaments in collapsing clouds, presents a suitable scenario for the fulfillment of the combined filament mass-environment criterium that is in quantitative agreement with Herschel observations.Peer reviewe