Showing 1–49 of 49 results for author: Pelkonen, V

Protoplanetary Disks from Pre-Main Sequence Bondi-Hoyle Accretion

Authors: Paolo Padoan, Liubin Pan, Veli-Matti Pelkonen, Troels Haugboelle, AAke Nordlund

Abstract: Protoplanetary disks are routinely described as finite mass reservoirs left over by the gravitational collapse of the protostar, an assumption that strongly constrains both disk evolution and planet formation models. We propose a different scenario where protoplanetary disks of pre-main sequence stars are assembled primarily by Bondi-Hoyle accretion from the parent gas cloud. We demonstrate that B… ▽ More Protoplanetary disks are routinely described as finite mass reservoirs left over by the gravitational collapse of the protostar, an assumption that strongly constrains both disk evolution and planet formation models. We propose a different scenario where protoplanetary disks of pre-main sequence stars are assembled primarily by Bondi-Hoyle accretion from the parent gas cloud. We demonstrate that Bondi-Hoyle accretion can supply not only the mass, but also the angular momentum necessary to explain the observed size of protoplanetary disks, and we predict the dependence of the disk specific angular momentum on the stellar mass. Our results are based on an analytical derivation of the scaling of the angular momentum in a turbulent flow, which we also confirm with a numerical simulation of supersonic turbulence. This new scenario for disk formation and evolution may alleviate a number of observational problems as well as compel major revisions of disk and planet formation models. △ Less

Submitted 12 May, 2024; originally announced May 2024.

Comments: Submitted to Nature Astronomy

Origin and Evolution of Angular Momentum of Class II Disks

Authors: Veli-Matti Pelkonen, Paolo Padoan, Mika Juvela, Troels Haugbølle, Åke Nordlund

Abstract: Context. While Class II stars have already accreted most of their mass, the continued inflow of fresh material via Bondi-Hoyle accretion acts as an additional mass reservoir for their circumstellar disks. This may explain the observed accretion rates of pre-main- sequence (PMS) stars, as well as observational inconsistencies in the mass and angular momentum balance of their disks. Aims. Using a ne… ▽ More Context. While Class II stars have already accreted most of their mass, the continued inflow of fresh material via Bondi-Hoyle accretion acts as an additional mass reservoir for their circumstellar disks. This may explain the observed accretion rates of pre-main- sequence (PMS) stars, as well as observational inconsistencies in the mass and angular momentum balance of their disks. Aims. Using a new simulation that reproduces the stellar IMF, we want to quantify the role of Bondi-Hoyle accretion in the formation of Class II disks, as well address the prospect of its observational detection with the James Webb Space Telescope (JWST). Methods. We study the mass and angular momentum of the accreting gas using passively-advected tracer particles in the simulation, and we carry out radiative transfer calculations of near-infrared scattering to generate synthetic JWST observations of Bondi-Hoyle trails of PMS stars. Results. Gas accreting on Class II PMS stars approximately 1 Myr after their formation has enough mass and angular momentum to strongly affect the evolution of the preexisting disks. The accreted angular momentum is large enough to also explain the observed size of Class II disks. The orientation of the angular momentum vector can differ significantly from that of the previously accreted gas, which may result in a significant disk warping or misalignment. We also predict that JWST observations of Class II stars will be able to detect Bondi-Hoyle trails with a 90% success rate with only 2 min exposure time, if stars with accretion rates \dot{M} > 5e-10 Msol/yr and luminosity of L > 0.5 Lsol are selected. △ Less

Submitted 10 May, 2024; originally announced May 2024.

Comments: 12 pages, 16 figures, submitted to A&A

Updated kinematics of the Radcliffe Wave: non-synchronous, dipole-like vertical oscillations

Authors: Zhi-Kai Zhu, Min Fang, Zu-Jia Lu, Junzhi Wang, Guang-Xing Li, Shiyu Zhang, Veli-Matti Pelkonen, Paolo Padoan, En-Wei Liang

Abstract: The kinematic structure of the Radcliffe Wave (RW) is crucial for understanding its origin and evolution. In this work, we present an accurate measurement of the vertical velocity $V_Z$ by where the radial velocity (RV) measures are taken into consideration. This is achieved in two ways. First, the velocities are measured towards Young Stellar Objects (YSOs), using their RV and proper motion measu… ▽ More The kinematic structure of the Radcliffe Wave (RW) is crucial for understanding its origin and evolution. In this work, we present an accurate measurement of the vertical velocity $V_Z$ by where the radial velocity (RV) measures are taken into consideration. This is achieved in two ways. First, the velocities are measured towards Young Stellar Objects (YSOs), using their RV and proper motion measurements from APOGEE-2 and Gaia DR3. Second, we combine RV measurements toward clouds with proper motion measurements of associated YSOs to determine the vertical velocities of the clouds. The results reveal that the oscillations in $V_Z$ are not synchronous with the vertical coordinate. The difference is caused by a combination of the effect of the radial velocity which we include in this paper, and the difference in models. By supplementing our analysis with additional young star samples, we find a consistent dipole pattern in $V_Z$. The fact that no significant amplitude differences are found among the analyzed samples indicates that there is no apparent age gradient within the dipole. We propose that RW evolves at a relatively slow rate. The fact that it will take a much longer time for RW to complete a full period compared to the cloud lifetimes challenges its classification as a traditional "wave". This age discrepancy should explain the phase difference, and non-synchronous oscillation found in kinematic studies. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 13 pages, 9 figures, submitted on 6 Feb 2024

Evaporation Ages: a New Dating Method for Young Star Clusters

Authors: V. -M. Pelkonen, N. Miret-Roig, P. Padoan

Abstract: 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.… ▽ More 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. △ Less

Submitted 16 February, 2024; v1 submitted 14 November, 2023; originally announced November 2023.

Comments: 13 pages, 11 figures, 2 tables, Submitted to A&A on Nov 14th 2023, Accepted 9 January 2024

Will ALMA Reveal the True Core Mass Function of Protoclusters?

Authors: Paolo Padoan, Veli-Matti Pelkonen, Mika Juvela, Troels Haugbølle, Åke Nordlund

Abstract: Characterizing prestellar cores in star-forming regions is an important step towards the validation of theoretical models of star formation. Thanks to their sub-arcsecond resolution, ALMA observations can potentially provide samples of prestellar cores up to distances of a few kpc, where regions of massive star formation can be targeted. However, the extraction of real cores from dust-continuum ob… ▽ More Characterizing prestellar cores in star-forming regions is an important step towards the validation of theoretical models of star formation. Thanks to their sub-arcsecond resolution, ALMA observations can potentially provide samples of prestellar cores up to distances of a few kpc, where regions of massive star formation can be targeted. However, the extraction of real cores from dust-continuum observations of turbulent star-forming clouds is affected by complex projection effects. In this work, we study the problem of core extraction both in the idealized case of column-density maps and in the more realistic case of synthetic 1.3\,mm ALMA observations. The analysis is carried out on 12 regions of high column density from our 250 pc simulation. We find that derived core masses are highly unreliable, with only {\em a weak correlation between the masses of cores selected in the synthetic ALMA maps and those of the corresponding three-dimensional cores}. The fraction of real three-dimensional cores detected in the synthetic maps increases monotonically with mass and remains always below 50\%. Above $\sim 1\,M_{\odot}$, the core mass function derived from the column-density maps is steeper than that of the three-dimensional cores, while the core mass function from the synthetic ALMA maps has a slope closer to that of the real three-dimensional cores. Because of the mass uncertainties, proper guidance from realistic simulations is essential if ALMA observations of protoclusters at kpc distances are to be used to test star-formation models. △ Less

Submitted 4 May, 2023; v1 submitted 18 January, 2023; originally announced January 2023.

Comments: 20 pages, 21 figures, MNRAS, in press

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Evidence for a Molecular Jet Launched at an Unprecedented Early Phase of Protostellar evolution

Authors: Somnath Dutta, Chin-Fei Lee, Naomi Hirano, Tie Liu, Doug Johnstone, Sheng-Yuan Liu, Kenichi Tatematsu, Paul F. Goldsmith, Dipen Sahu, Neal J. Evans, Patricio Sanhueza, Woojin Kwon, Sheng-Li Qin, Manash Ranjan Samal, Qizhou Zhang, Kee-Tae Kim, Hsien Shang, Chang Won Lee, Anthony Moraghan, Kai-Syun Jhan, Shanghuo Li, Jeong-Eun Lee, Alessio Traficante, Mika Juvela, Leonardo Bronfman , et al. (9 additional authors not shown)

Abstract: Protostellar outflows and jets play a vital role in star formation as they carry away excess angular momentum from the inner disk surface, allowing the material to be transferred toward the central protostar. Theoretically, low velocity and poorly collimated outflows appear from the beginning of the collapse, at the first hydrostatic core (FHSC) stage. With growing protostellar core mass, high-den… ▽ More Protostellar outflows and jets play a vital role in star formation as they carry away excess angular momentum from the inner disk surface, allowing the material to be transferred toward the central protostar. Theoretically, low velocity and poorly collimated outflows appear from the beginning of the collapse, at the first hydrostatic core (FHSC) stage. With growing protostellar core mass, high-density jets are launched which entrain an outflow from the infalling envelope. Until now, molecular jets have been observed at high velocity ($\gtrsim$ 100 km/s) in early Class\,0 protostars. We, for the first time, detect a dense molecular jet in SiO emission with small-velocity ($\sim$ 4.2 km\,s$^{-1}$, deprojected $\sim$ 24 km\,s$^{-1}$) from source G208.89-20.04Walma (hereafter, G208Walma) using ALMA Band\,6 observations. This object has some characteristics of FHSCs, such as a small outflow/jet velocity, extended 1.3\,mm continuum emission, and N$_2$D$^+$ line emission. Additional characteristics, however, are typical of early protostars: collimated outflow and SiO jet. The full extent of the outflow corresponds to a dynamical time scale of $\sim$ 930$^{+200}_{-100}$ years. The spectral energy distribution also suggests a very young source having an upper limit of T$_{bol}$ $\sim$ 31 K and L$_{bol}$ $\sim$ 0.8 L$_\sun$. We conclude that G208Walma is likely in the transition phase from FHSC to protostar, and the molecular jet has been launched within a few hundred years of initial collapse. Therefore, G208Walma may be the earliest object discovered in the protostellar phase with a molecular jet. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: 16 pages, 7 figures. accepted for publication in ApJ

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): How do dense core properties affect the multiplicity of protostars?

Authors: Qiuyi Luo, Tie Liu, Kenichi Tatematsu, ShengYuan Liu, Pak Shing Li, James di Francesco, Doug Johnstone, Paul F. Goldsmith, Somnath Dutta, Naomi Hirano, ChinFei Lee, Di Li, KeeTae Kim, Chang Won Lee, JeongEun Lee, Xunchuan Liu, Mika Juvela, Jinhua He, ShengLi Qin, HongLi Liu, David Eden, Woojin Kwon, Dipen Sahu, Shanghuo Li, FengWei Xu , et al. (12 additional authors not shown)

Abstract: During the transition phase from a prestellar to a protostellar cloud core, one or several protostars can form within a single gas core. The detailed physical processes of this transition, however, still remain unclear. We present 1.3 mm dust continuum and molecular line observations with the Atacama Large Millimeter/submillimeter Array (ALMA) toward 43 protostellar cores in the Orion Molecular Cl… ▽ More During the transition phase from a prestellar to a protostellar cloud core, one or several protostars can form within a single gas core. The detailed physical processes of this transition, however, still remain unclear. We present 1.3 mm dust continuum and molecular line observations with the Atacama Large Millimeter/submillimeter Array (ALMA) toward 43 protostellar cores in the Orion Molecular Cloud Complex ($λ$ Orionis, Orion B, and Orion A) with an angular resolution of $\sim$ 0.35" ($\sim$ 140 au). In total, we detect 13 binary/multiple systems. We derive an overall multiplicity frequency (MF) of 28$\%$ $\pm$ 4$\%$ and a companion star fraction (CSF) of 51$\%$ $\pm$ 6$\%$, over a separation range of 300-8900 au. The median separation of companions is about 2100 au. The occurrence of stellar multiplicity may depend on the physical characteristics of the dense cores. Notably, those containing binary/multiple systems tend to show higher gas density and Mach number than cores forming single stars. The integral-shaped filament (ISF) of Orion A giant molecular cloud (GMC), which has the highest gas density and hosts high-mass star formation in its central region (the Orion Nebula cluster), shows the highest MF and CSF among the Orion GMCs. In contrast, the $λ$ Orionis Giant Molecular Cloud (GMC) has a lower MF and CSF than the Orion B and Orion A GMCs, indicating that feedback from HII regions may suppress the formation of multiple systems. We also find that the protostars comprising a binary/multiple system are usually at different evolutionary stages. △ Less

Submitted 9 May, 2022; v1 submitted 13 April, 2022; originally announced April 2022.

Comments: This paper was accepted by APJ. 31 pages and 9 figures

The Dynamical State of Massive Clumps

Authors: Zu-Jia Lu, Veli-Matti Pelkonen, Mika Juvela, Paolo Padoan, Troels Haugbølle, Åke Nordlund

Abstract: The dynamical state of massive clumps is key to our understanding of the formation of massive stars. In this work, we study the kinematic properties of massive clumps using synthetic observations. We have previously compiled a very large catalog of synthetic dust-continuum compact sources from our 250 pc, SN-driven, star formation simulation. Here, we compute synthetic $\rm N_{2}H^{+}$ line profil… ▽ More The dynamical state of massive clumps is key to our understanding of the formation of massive stars. In this work, we study the kinematic properties of massive clumps using synthetic observations. We have previously compiled a very large catalog of synthetic dust-continuum compact sources from our 250 pc, SN-driven, star formation simulation. Here, we compute synthetic $\rm N_{2}H^{+}$ line profiles for a subsample of those sources and compare their properties with the observations and with those of the corresponding three-dimensional (3D) clumps in the simulation. We find that the velocity dispersion of the sources estimated from the $\rm N_{2}H^{+}$ line is a good estimate of that of the 3D clumps, although its correlation with the source size is weaker than the velocity-size correlation of the 3D clumps. The relation between the mass of the 3D clumps, $M_{\rm main}$, and that of the corresponding synthetic sources, $M_{\rm SED}$, has a large scatter and a slope of 0.5, $M_{\rm main} \propto M_{\rm SED}^{0.5}$, due to uncertainties arising from the observational band-merging procedure and from projection effects along the line of sight. As a result, the virial parameters of the 3D clumps are not correlated with the clump masses, even if a negative correlation is found for the compact sources, and the virial parameter of the most massive sources may significantly underestimate that of the associated clumps. △ Less

Submitted 16 November, 2021; originally announced November 2021.

Comments: 19 pages, 12 figures, 1 table. Accepted for publication in MNRAS

Physical Properties and Real Nature of Massive Clumps in the Galaxy

Authors: Zu-Jia Lu, Veli-Matti Pelkonen, Mika Juvela, Paolo Padoan, Troels Haugbølle, Åke Nordlund

Abstract: Systematic surveys of massive clumps have been carried out to study the conditions leading to the formation of massive stars. These clumps are typically at large distances and unresolved, so their physical properties cannot be reliably derived from the observations alone. Numerical simulations are needed to interpret the observations. To this end, we generate synthetic Herschel observations using… ▽ More Systematic surveys of massive clumps have been carried out to study the conditions leading to the formation of massive stars. These clumps are typically at large distances and unresolved, so their physical properties cannot be reliably derived from the observations alone. Numerical simulations are needed to interpret the observations. To this end, we generate synthetic Herschel observations using our large-scale star-formation simulation, where massive stars explode as supernovae driving the interstellar-medium turbulence. From the synthetic observations, we compile a catalog of compact sources following the exact same procedure as for the Hi-GAL compact source catalog. We show that the sources from the simulation have observational properties with statistical distributions consistent with the observations. By relating the compact sources from the synthetic observations to their three-dimensional counterparts in the simulation, we find that the synthetic observations overestimate the clump masses by about an order of magnitude on average due to line-of-sight projection, and projection effects are likely to be even worse for Hi-GAL Inner Galaxy sources. We also find that a large fraction of sources classified as protostellar are likely to be starless, and propose a new method to partially discriminate between true and false protostellar sources. △ Less

Submitted 19 August, 2021; originally announced August 2021.

Comments: 17 pages, 15 figures, submitted to MNRAS

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Detection of extremely high density compact structure of prestellar cores and multiple substructures within

Authors: Dipen Sahu, Sheng-Yuan Liu, Tie Liu, Neal J. Evans II, Naomi Hirano, Ken'ichi Tatematsu, Chin-Fei Lee, Kee-Tae Kim, Somnath Dutta, Dana Alina, Leonardo Bronfman, Maria Cunningham, David J. Eden, Guido Garay, Paul F. Goldsmith, Jinhua He, Shih-Ying Hsu, Kai-Syun Jhan, Doug Johnstone, Mika Juvela, Gwanjeong Kim, Yi-Jehng Kuan, Woojin Kwon, Chang Won Lee, Jeong-Eun Lee , et al. (16 additional authors not shown)

Abstract: Prestellar cores are self-gravitating dense and cold structures within molecular clouds where future stars are born. They are expected, at the stage of transitioning to the protostellar phase, to harbor centrally concentrated dense (sub)structures that will seed the formation of a new star or the binary/multiple stellar systems. Characterizing this critical stage of evolution is key to our underst… ▽ More Prestellar cores are self-gravitating dense and cold structures within molecular clouds where future stars are born. They are expected, at the stage of transitioning to the protostellar phase, to harbor centrally concentrated dense (sub)structures that will seed the formation of a new star or the binary/multiple stellar systems. Characterizing this critical stage of evolution is key to our understanding of star formation. In this work, we report the detection of high density (sub)structures on the thousand-au scale in a sample of dense prestellar cores. Through our recent ALMA observations towards the Orion molecular cloud, we have found five extremely dense prestellar cores, which have centrally concentrated regions $\sim$ 2000 au in size, and several $10^7$ $cm^{-3}$ in average density. Masses of these centrally dense regions are in the range of 0.30 to 6.89 M$_\odot$. {\it For the first time}, our higher resolution observations (0.8$'' \sim $ 320 au) further reveal that one of the cores shows clear signatures of fragmentation; such individual substructures/fragments have sizes of 800 -1700 au, masses of 0.08 to 0.84 M$_\odot$, densities of $2 - 8\times 10^7$ $cm^{-3}$ and separations of $\sim 1200$ au. The substructures are massive enough ($\gtrsim 0.1~M_\odot$) to form young stellar objects and are likely examples of the earliest stage of stellar embryos which can lead to widely ($\sim$ 1200 au) separated multiple systems. △ Less

Submitted 12 February, 2023; v1 submitted 15 December, 2020; originally announced December 2020.

Comments: Published in The Astrophysical Journal Letters (ApJL)

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP) II. Survey overview: a first look at 1.3 mm continuum maps and molecular outflows

Authors: Somnath Dutta, Chin-Fei Lee, Tie Liu, Naomi Hirano, Sheng-Yuan Liu, Kenichi Tatematsu, Kee-Tae Kim, Hsien Shang, Dipen Sahu, Gwanjeong Kim, Anthony Moraghan, Kai-Syun Jhan, Shih-Ying Hsu, Neal J. Evans, Doug Johnstone, Derek Ward-Thompson, Yi-Jehng Kuan, Chang Won Lee, Jeong-Eun Lee, Alessio Traficante, Mika Juvela, Charlotte Vastel, Qizhou Zhang, Patricio Sanhueza, Archana Soam , et al. (10 additional authors not shown)

Abstract: Planck Galactic Cold Clumps (PGCCs) are contemplated to be the ideal targets to probe the early phases of star formation. We have conducted a survey of 72 young dense cores inside PGCCs in the Orion complex with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.3\,mm (band 6) using three different configurations (resolutions $\sim$ 0$\farcs$35, 1$\farcs$0, and 7$\farcs$0) to statistical… ▽ More Planck Galactic Cold Clumps (PGCCs) are contemplated to be the ideal targets to probe the early phases of star formation. We have conducted a survey of 72 young dense cores inside PGCCs in the Orion complex with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.3\,mm (band 6) using three different configurations (resolutions $\sim$ 0$\farcs$35, 1$\farcs$0, and 7$\farcs$0) to statistically investigate their evolutionary stages and sub-structures. We have obtained images of the 1.3\,mm continuum and molecular line emission ($^{12}$CO, and SiO) at an angular resolution of $\sim$ 0$\farcs$35 ($\sim$ 140\,au) with the combined arrays. We find 70 substructures within 48 detected dense cores with median dust-mass $\sim$ 0.093\,M$_{\sun}$ and deconvolved size $\sim$ 0$\farcs$27. Dense substructures are clearly detected within the central 1000\,au of four candidate prestellar cores. The sizes and masses of the substructures in continuum emission are found to be significantly reduced with protostellar evolution from Class\,0 to Class\,I. We also study the evolutionary change in the outflow characteristics through the course of protostellar mass accretion. A total of 37 sources exhibit CO outflows, and 20 ($>$50\%) show high-velocity jets in SiO. The CO velocity-extents ($Δ$Vs) span from 4 to 110 km/s with outflow cavity opening angle width at 400\,au ranging from $[Θ_{obs}]_{400}$ $\sim$ 0$\farcs$6 to 3$\farcs$9, which corresponds to 33$\fdg$4$-$125$\fdg$7. For the majority of the outflow sources, the $Δ$Vs show a positive correlation with $[Θ_{obs}]_{400}$, suggesting that as protostars undergo gravitational collapse, the cavity opening of a protostellar outflow widens and the protostars possibly generate more energetic outflows. △ Less

Submitted 27 October, 2020; originally announced October 2020.

Comments: 43 pages, 25 figures (14 in the main text, 9 in the appendix), accepted for publication in the Astrophysical Journal Supplement Series (ApJS)

Dust emission, extinction, and scattering in LDN 1642

Authors: Mika Juvela, Sharma Neha, Emma Mannfors, Mika Saajasto, Nathalie Ysard, Veli-Matti Pelkonen

Abstract: We study the near-infrared (NIR) scattering in LDN 1642, its correlation with the cloud structure, and the ability of dust models to simultaneously explain sub-millimetre emission, NIR extinction, and NIR scattering. We use observations from the HAWK-I instrument to measure the NIR surface brightness and extinction. These are compared with Herschel data on dust emission and, with radiative trans… ▽ More We study the near-infrared (NIR) scattering in LDN 1642, its correlation with the cloud structure, and the ability of dust models to simultaneously explain sub-millimetre emission, NIR extinction, and NIR scattering. We use observations from the HAWK-I instrument to measure the NIR surface brightness and extinction. These are compared with Herschel data on dust emission and, with radiative transfer modelling, with predictions calculated for different dust models. We find an optical depth ratio $τ(250\,μ{\rm m})/τ(J)\approx 10^{-3}$, confirming earlier findings of high sub-millimetre emissivity. The relationships between the column density derived from dust emission and the NIR colour excesses is linear and consistent with the standard NIR extinction curve. The extinction peaks at $A_J=2.6\,$mag, the NIR surface brightness remaining correlated with $N({\rm H}_2)$ without saturation. Radiative transfer models can fit the sub-millimetre data with any of the tested dust models. However, these predict a NIR extinction that is higher and a NIR surface brightness that is lower than in observations. If the dust sub-millimetre emissivity is rescaled to the observed value of $τ(250\,μ{\rm m})/τ(J)$, dust models with high NIR albedo can reach the observed level of NIR surface brightness. The NIR extinction of the models tends to be higher than directly measured, which is reflected in the shape of the NIR surface brightness spectra. The combination of emission, extinction, and scattering measurements provides strong constraints on dust models. The observations of LDN 1642 indicate clear dust evolution, including a strong increase in the sub-millimetre emissivity, not yet fully explained by the current dust models. △ Less

Submitted 26 September, 2020; originally announced September 2020.

Comments: 16 pages, 16 figures, submitted to A&A

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

Authors: V. -M. Pelkonen, P. Padoan, T. Haugbølle, Å. Nordlund

Abstract: 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,… ▽ More 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. △ Less

Submitted 14 March, 2021; v1 submitted 5 August, 2020; originally announced August 2020.

Comments: 23 pages, 29 figures. Link to supplementary material and full Table 1: http://www.erda.dk/vgrid/core-mass-function/ . Submitted to MNRAS

The Effect of Supernovae on the Turbulence and Dispersal of Molecular Clouds

Authors: Zu-Jia Lu, Veli-Matti Pelkonen, Paolo Padoan, Liubin Pan, Troels Haugbølle, Åke Nordlund

Abstract: While the importance of supernova feedback in galaxies is well established, its role on the scale of molecular clouds is still debated. In this work, we focus on the impact of supernovae on individual clouds, using a high-resolution magneto-hydrodynamic simulation of a region of 250 pc where we resolve the formation of individual massive stars. The supernova feedback is implemented with real super… ▽ More While the importance of supernova feedback in galaxies is well established, its role on the scale of molecular clouds is still debated. In this work, we focus on the impact of supernovae on individual clouds, using a high-resolution magneto-hydrodynamic simulation of a region of 250 pc where we resolve the formation of individual massive stars. The supernova feedback is implemented with real supernovae that are the natural evolution of the resolved massive stars, so their position and timing are self-consistent. We select a large sample of molecular clouds from the simulation to investigate the supernova energy injection and the resulting properties of molecular clouds. We find that molecular clouds have a lifetime of a few dynamical times, less then half of them contract to the point of becoming gravitationally bound, and the dispersal time of bound clouds, of order one dynamical time, is a factor of two shorter than that of unbound clouds. We stress the importance of internal supernovae, that is massive stars that explode inside their parent cloud, in setting the cloud dispersal time, and their huge overdensity compared to models where the supernovae are randomly distributed. We also quantify the energy injection efficiency of supernovae as a function of supernova distance to the clouds. We conclude that intermittent driving by supernovae can maintain molecular-cloud turbulence and may be the main process of cloud dispersal. The role of supernovae in the evolution of molecular clouds cannot be fully accounted for without a self-consistent implementation of their feedback. △ Less

Submitted 18 July, 2020; originally announced July 2020.

Comments: 33 pages, 23 figures, submitted to ApJ

Multi-scale analysis of the Monoceros OB 1 star-forming region: II. Colliding filaments in the Monoceros OB1 molecular cloud

Authors: Julien Montillaud, Mika Juvela, Charlotte Vastel, J. H. He, Tie Liu, Isabelle Ristorcelli, David J. Eden, Sung-ju Kang, Kee-Tae Kim, Patrick M. Koch, Chang Won Lee, Mark G. Rawlings, Mika Saajasto, Patricio Sanhueza, Archana Soam, Sarolta Zahorecz, Dana Alina, Rebeka Bögner, David Cornu, Yasuo Doi, Johanna Malinen, Douglas Marshall, E. R. Micelotta, V. M. Pelkonen, L. V. Tóth , et al. (2 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 13 September, 2019; originally announced September 2019.

Comments: A&A, in press; IRAM and TRAO data cubes will be available at CDS

Journal ref: A&A 631, A3 (2019)

Multi-scale analysis of the Monoceros OB 1 star-forming region: I. The dense core population

Authors: Julien Montillaud, Mika Juvela, Charlotte Vastel, J. H. He, Tie Liu, Isabelle Ristorcelli, David J. Eden, Sung-ju Kang, Kee-Tae Kim, Patrick M. Koch, Chang Won Lee, Mark G. Rawlings, Mika Saajasto, Patricio Sanhueza, Archana Soam, Sarolta Zahorecz, Dana Alina, Rebeka Bögner, David Cornu, Yasuo Doi, Johanna Malinen, Douglas Marshall, E. R. Micelotta, V. M. Pelkonen, L. V. Tóth , et al. (2 additional authors not shown)

Abstract: 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.… ▽ More 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. △ Less

Submitted 9 September, 2019; originally announced September 2019.

Journal ref: A&A 631, L1 (2019)

Synthetic observations of dust emission and polarisation of Galactic cold clumps

Authors: Mika Juvela, Paolo Padoan, Isabelle Ristorcelli, Veli-Matti Pelkonen

Abstract: 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… ▽ More 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 Δ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. △ Less

Submitted 9 August, 2019; originally announced August 2019.

Comments: Accepted to A&A (20 pages, 26 figures)

SCOPE: SCUBA-2 Continuum Observations of Pre-protostellar Evolution - Survey Description and Compact Source Catalogue

Authors: D. J. Eden, Tie Liu, Kee-Tae Kim, S. -Y. Liu, K. Tatematsu, J. Di Francesco, K. Wang, Y. Wu, M. A. Thompson, G. A. Fuller, Di Li, I. Ristorcelli, Sung-ju Kang, N. Hirano, D. Johnstone, Y. Lin, J. H. He, P. M. Koch, Patricio Sanhueza, S. -L. Qin, Q. Zhang, P. F. Goldsmith, N. J. Evans II, J. Yuan, C. -P. Zhang , et al. (136 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 26 February, 2019; originally announced February 2019.

Comments: 16 pages, 13 figures, 1 table. Accepted for publication in MNRAS

Expanding associations in the Vela-Puppis region: 3D structure and kinematics of the young population

Authors: T. Cantat-Gaudin, C. Jordi, N. J. Wright, J. J. Armstrong, A. Vallenari, L. Balaguer-Núñez, P. Ramos, D. Bossini, P. Padoan, V. M. Pelkonen, M. Mapelli, R. D. Jeffries

Abstract: 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 Ga… ▽ More 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. △ Less

Submitted 7 May, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: accepted A&A

Journal ref: A&A 626, A17 (2019)

Dust spectrum and polarisation at 850 um in the massive IRDC G035.39-00.33

Authors: Mika Juvela, Vincent Guillet, Tie Liu, Isabelle Ristorcelli, Veli-Matti Pelkonen, Dana Alina, Leonardo Bronfman, David J. Eden, Kee Tae Kim, Patrick M. Koch, Woojin Kwon, Chang Won Lee, Johanna Malinen, Elisabetta Micelotta, Julien Montillaud, Mark G. Rawlings, Patricio Sanhueza, Archana Soam, Alessio Traficante, Nathalie Ysard, Chuan-Peng Zhang

Abstract: 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… ▽ More 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. △ Less

Submitted 6 September, 2018; v1 submitted 4 September, 2018; originally announced September 2018.

Comments: accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 620, A26 (2018)

Galactic cold cores IX. Column density structures and radiative transfer modelling

Authors: M. Juvela, J. Malinen, J. Montillaud, V. -M. Pelkonen, I. Ristorcelli, L. V. Toth

Abstract: 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… ▽ More 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. △ Less

Submitted 15 January, 2018; v1 submitted 8 January, 2018; originally announced January 2018.

Comments: accepted to A&A, 21 pages (37 pages with appendices, this version including language corrections)

Journal ref: A&A 614, A83 (2018)

Herschel and SCUBA-2 observations of dust emission in a sample of Planck cold clumps

Authors: Mika Juvela, Jinhua He, Katherine Pattle, Tie Liu, George Bendo, David J. Eden, Orsolya Feher, Michel Fich, Gary Fuller, Naomi Hirano, Kee-Tae Kim, Di Li, Sheng-Yuan Liu, Johanna Malinen, Douglas J. Marshall, Deborah Paradis, Harriet Parsons, Veli-Matti Pelkonen, Mark G. Rawlings, Isabelle Ristorcelli, Manash R. Samal, Ken'ichi Tatematsu, Mark Thompson, Alessio Traficante, Ke Wang , et al. (4 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 27 December, 2017; v1 submitted 26 November, 2017; originally announced November 2017.

Comments: 37 pages, accepted to A&A, final version with language corrections

Journal ref: A&A 612, A71 (2018)

Correlation of gas dynamics and dust in the evolved filament G82.65-02.00

Authors: M. Saajasto, M. Juvela, K. Dobashi, T. Shimoikura, I. Ristorcelli, J. Montillaud, D. J. Marshall, J. Malinen, V. -M. Pelkonen, O. Fehér, A. Rivera-Ingraham, L. V. Toth, L. Montier, J-Ph. Bernard, T. Onishi

Abstract: 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… ▽ More 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 $μ$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$. △ Less

Submitted 3 July, 2017; originally announced July 2017.

Comments: Accepted for publication in A&A

Journal ref: A&A 608, A21 (2017)

Galactic Cold Cores. VIII. Filament formation and evolution: Filament properties in context with evolutionary models

Authors: A. Rivera-Ingraham, I. Ristorcelli, M. Juvela, J. Montillaud, A. Men'shchikov, J. Malinen, V. -M. Pelkonen, A. Marston, P. G. Martin, L. Pagani, R. Paladini, D. Paradis, N. Ysard, D. Ward-Thompson, J. -P. Bernard, D. J. Marshall, L. Montier, L. V. Toth

Abstract: 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… ▽ More 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. △ Less

Submitted 24 February, 2017; originally announced February 2017.

Comments: Accepted in A&A

Journal ref: A&A 601, A94 (2017)

Galactic Cold Cores VII: Filament Formation and Evolution - Methods & Observational Constraints

Authors: A. Rivera-Ingraham, I. Ristorcelli, M. Juvela, J. Montillaud, A. Men'shchikov, J. Malinen, V. -M. Pelkonen, A. Marston, P. G. Martin, L. Pagani, R. Paladini, D. Paradis, N. Ysard, D. Ward-Thompson, J. -P. Bernard, D. J. Marshall, L. Montier, V. Tóth

Abstract: 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… ▽ More 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. △ Less

Submitted 30 March, 2016; originally announced March 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 591, A90 (2016)

Structure and stability in TMC-1: analysis of NH$_3$ molecular line and Herschel continuum data

Authors: O. Fehér, L. V. Tóth, D. Ward-Thompson, J. Kirk, A. Kraus, V. -M. Pelkonen, S. Pintér, S. Zahorecz

Abstract: 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… ▽ More 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. △ Less

Submitted 18 March, 2016; originally announced March 2016.

Comments: 17 pages, 9 figures; accepted for publication in A&A

Journal ref: A&A 590, A75 (2016)

Matching dust emission structures and magnetic field in high-latitude cloud L1642: comparing Herschel and Planck maps

Authors: J. Malinen, L. Montier, J. Montillaud, M. Juvela, I. Ristorcelli, S. E. Clark, O. Berné, J. -Ph. Bernard, V. -M. Pelkonen, D. C. Collins

Abstract: 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… ▽ More 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. △ Less

Submitted 3 May, 2016; v1 submitted 11 December, 2015; originally announced December 2015.

Comments: accepted to MNRAS; 14 pages; changes: moderate text and figure editing

Galactic cold cores VI. Dust opacity spectral index

Authors: M. Juvela, K. Demyk, Y. Doi, A. Hughes, C. Lefevre, D. J. Marshall, C. Meny, J. Montillaud, L. Pagani, D. Paradis, I. Ristorcelli, J. Malinen, L. A. Montier, R. Paladini, V. -M. Pelkonen, A. Rivera-Ingraham

Abstract: 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… ▽ More 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. △ Less

Submitted 26 September, 2015; originally announced September 2015.

Comments: Accepted to A&A, 138 pages

Journal ref: A&A 584, A94 (2015)

Grain size limits derived from 3.6 μm and 4.5 μm coreshine

Authors: J. Steinacker, M. Andersen, W. -F. Thi, R. Paladini, M. Juvela, A. Bacmann, V. -M. Pelkonen, L. Pagani, C. Lefèvre, Th. Henning, A. Noriega-Crespo

Abstract: Recently discovered scattered light from molecular cloud cores in the wavelength range 3-5 μm (called "coreshine") seems to indicate the presence of grains with sizes above 0.5 μm. We aim to analyze 3.6 and 4.5 μ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 μm bands for a… ▽ More Recently discovered scattered light from molecular cloud cores in the wavelength range 3-5 μm (called "coreshine") seems to indicate the presence of grains with sizes above 0.5 μm. We aim to analyze 3.6 and 4.5 μ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 μ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 μ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 μ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. △ Less

Submitted 17 August, 2015; originally announced August 2015.

Comments: 23 pages, 18 figures, accepted for publication in A&A

The physical state of selected cold clumps

Authors: A. Parikka, M. Juvela, V. -M. Pelkonen, J. Malinen, J. Harju

Abstract: 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.… ▽ More 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. △ Less

Submitted 20 March, 2015; originally announced March 2015.

Comments: Accepted to A&A

Journal ref: A&A 577, A69 (2015)

Planck 2015 results. XXVIII. The Planck Catalogue of Galactic Cold Clumps

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, N. Bartolo, E. Battaner, K. Benabed, A. Benoît, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana , et al. (198 additional authors not shown)

Abstract: 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 cur… ▽ More 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. △ Less

Submitted 5 February, 2015; originally announced February 2015.

Planck 2015 results. I. Overview of products and scientific results

Authors: Planck Collaboration, R. Adam, P. A. R. Ade, N. Aghanim, Y. Akrami, M. I. R. Alves, M. Arnaud, F. Arroja, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, J. G. Bartlett, N. Bartolo, S. Basak, P. Battaglia, E. Battaner, R. Battye, K. Benabed, A. Benoît, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, B. Bertincourt , et al. (330 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 9 August, 2015; v1 submitted 5 February, 2015; originally announced February 2015.

Comments: 40 pages, 24 figures

Galactic cold cores V. Dust opacity

Authors: M. Juvela, I. Ristorcelli, D. J. Marshall, J. Montillaud, V. -M. Pelkonen, N. Ysard, P. McGehee, R. Paladini, L. Pagani, J. Malinen, A., Rivera-Ingraham, C. Lefevre, L. V. Toth, L. A. Montier, J. -P. Bernard, P. Martin

Abstract: 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 observ… ▽ More 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. △ Less

Submitted 28 January, 2015; originally announced January 2015.

Comments: 28 pages, A&A, in press

Can we trace very cold dust from its emission alone ?

Authors: Laurent Pagani, Charlène Lefèvre, Mika Juvela, Veli-Matti Pelkonen, Frédéric Schuller

Abstract: 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 acr… ▽ More 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 $μ$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. △ Less

Submitted 5 January, 2015; originally announced January 2015.

Comments: Letter to A&A (accepted for publication). must be viewed with ACROBAT READER for full enhancement. Otherwise, check images in Appendix A

Journal ref: A&A 574, L5 (2015)

Dust properties inside molecular clouds from coreshine modeling and observations

Authors: Charlène Lefèvre, Laurent Pagani, Mika Juvela, Roberta Paladini, Rosine Lallement, D. J. Marshall, Morten Andersen, Aurore Bacmann, Peregrine M. McGehee, Ludovic Montier, Alberto Noriega-Crespo, V. -M. Pelkonen, Isabelle Ristorcelli, Jürgen Steinacker

Abstract: 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 $μ$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 inv… ▽ More 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 $μ$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 $μ$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 $μ$m. For starless cores, where detected, the observed 4.5 $μ$m / 3.6 $μ$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. △ Less

Submitted 22 July, 2014; originally announced July 2014.

Journal ref: A&A 572, A20 (2014)

Planck intermediate results. XX. Comparison of polarized thermal emission from Galactic dust with simulations of MHD turbulence

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, D. Alina, M. I. R. Alves, G. Aniano, C. Armitage-Caplan, M. Arnaud, D. Arzoumanian, M. Ashdown, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, E. Battaner, K. Benabed, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger , et al. (164 additional authors not shown)

Abstract: 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 $ψ$. 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… ▽ More 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 $ψ$. 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] △ Less

Submitted 25 September, 2014; v1 submitted 5 May, 2014; originally announced May 2014.

Comments: Part of a series of four papers, accepted for publication in A&A

Journal ref: A&A 576, A105 (2015)

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

Authors: J. Malinen, M. Juvela, V. -M. Pelkonen, M. G. Rawlings

Abstract: 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 emiss… ▽ More 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. △ Less

Submitted 9 April, 2014; originally announced April 2014.

Comments: abstract shortened and figures reduced for astroph

Journal ref: A&A 558, A44 (2013)

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

Authors: J. Malinen, M. Juvela, S. Zahorecz, A. Rivera-Ingraham, J. Montillaud, K. Arimatsu, J. -Ph. Bernard, Y. Doi, L. Haikala, R. Kawabe, G. Marton, P. McGehee, V. -M. Pelkonen, I. Ristorcelli, Y. Shimajiri, S. Takita, L. V. Toth, T. Tsukagoshi, N. Ysard

Abstract: 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… ▽ More 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. △ Less

Submitted 9 April, 2014; v1 submitted 11 February, 2014; originally announced February 2014.

Comments: 22 pages, 18 figures, accepted to Astronomy & Astrophysics; abstract shortened and figures reduced for astroph

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

Authors: M. Juvela, V. -M. Pelkonen, G. J. White, V. Konyves, J. Kirk, P. Andre

Abstract: 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 d… ▽ More 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. △ Less

Submitted 23 June, 2012; originally announced June 2012.

Comments: 13 pages, 16 figures, accepted to A&A

Galactic cold cores III. General cloud properties

Authors: M. Juvela, I. Ristorcelli, L. Pagani, Y. Doi, V. -M. Pelkonen, D. J. Marshall, J. -P. Bernard, E. Falgarone, J. Malinen, G. Marton, P. McGehee, L. A. Montier, F. Motte, R. Paladini, L. V. Toth, N. Ysard, S. Zahorecz, A. Zavagno

Abstract: 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 distr… ▽ More 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. △ Less

Submitted 8 February, 2012; originally announced February 2012.

Comments: Accepted to A&A. The pdf with high resolution figures and with the figures of the appendices is available at https://wiki.helsinki.fi/display/PlanckHerschel/Publications

Galactic cold cores II. Herschel study of the extended dust emission around the first Planck detections

Authors: M. Juvela, I. Ristorcelli, V. -M. Pelkonen, D. J. Marshall, L. A. Montier, J. -P. Bernard, R. Paladini, T. Lunttila, A. Abergel, Ph. Andre, C. Dickinson, X. Dupac, J. Malinen, P. Martin, P. McGehee, L. Pagani, N. Ysard, A. Zavagno

Abstract: 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-millim… ▽ More 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. △ Less

Submitted 15 January, 2011; originally announced January 2011.

Comments: Accepted to A&A

Planck Early Results: The Galactic Cold Core Population revealed by the first all-sky survey

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. Balbi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoît, J. -P. Bernard, M. Bersanelli, R. Bhatia, J. J. Bock, A. Bonaldi, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana , et al. (180 additional authors not shown)

Abstract: 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 comp… ▽ More 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 β 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 β over the all temperature range. β ranges from 1.4 to 2.8 with a mean value around 2.1, and several possible scenarios are possible, including β(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. △ Less

Submitted 11 January, 2011; originally announced January 2011.

Comments: This paper is part of package of Planck early papers, labelled in their title by "Planck Early results:..." Corresponding Author: L. Montier (Ludovic.Montier@cesr.fr)

Planck Early Results: The submillimetre properties of a sample of Galactic cold clumps

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. Balbi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoît, J. -P. Bernard, M. Bersanelli, R. Bhatia, J. J. Bock, A. Bonaldi, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana , et al. (179 additional authors not shown)

Abstract: (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 sour… ▽ More (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. △ Less

Submitted 11 January, 2011; originally announced January 2011.

Comments: This paper is part of a package of papers describing first results of the Planck mission, to be published together in a special issue of A&A. Corresponding author: Isabelle.Ristorcelli@cesr.fr. Full resolution version available from http://www.sciops.esa.int/index.php?project=PLANCK&page=Planck_Published_Papers

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

Authors: Mika Juvela, Veli-Matti Pelkonen, Sebastian Porceddu

Abstract: 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 s… ▽ More 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. △ Less

Submitted 27 July, 2009; originally announced July 2009.

Comments: 9 pages, 9 figures, accepted to A&A

Predictions of polarized dust emission from interstellar clouds: spatial variations in the efficiency of radiative torque alignment

Authors: V. -M. Pelkonen, M. Juvela, P. Padoan

Abstract: 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… ▽ More 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 $γ= 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}$. △ Less

Submitted 18 June, 2009; v1 submitted 7 January, 2009; originally announced January 2009.

Comments: 12 pages, 15 figures, submitted to A&A 19.12.2008; 09.01.2009: Corrected the name of Juvela; 24.04.2009: revised, added content, 13 pages, 16 figures; 18.06.2009: Language edited, print version

A Corona Australis cloud filament seen in NIR scattered light I. Comparison with extinction of background stars

Authors: M. Juvela, V. -M. Pelkonen, P. Padoan, K. Mattila

Abstract: 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 confi… ▽ More 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. △ Less

Submitted 9 January, 2008; v1 submitted 12 August, 2007; originally announced August 2007.

Comments: 14 pages, 15 figures, accepted to A&A, the version includes small changes in the text and an added appendix

Simulations of polarized dust emission

Authors: V. -M. Pelkonen, M. Juvela, P. Padoan

Abstract: 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, wh… ▽ More 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'. △ Less

Submitted 7 September, 2006; v1 submitted 6 September, 2006; originally announced September 2006.

Comments: Accepted to A&A 31/08/2006. Replaced 07/09/2006 to make minor changes to make the author and affiliation list to conform to standards and tradition

Journal ref: Astron.Astrophys. 461 (2007) 551-564

High-resolution mapping of interstellar clouds with near-infrared scattered light

Authors: M. Juvela, V. -M. Pelkonen, P. Padoan, K. Mattila

Abstract: 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… ▽ More 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. △ Less

Submitted 5 July, 2006; v1 submitted 14 March, 2006; originally announced March 2006.

Comments: 20 pages, accepted to A&A; this version contains some structural changes and two new figures

Journal ref: Astron.Astrophys. 457 (2006) 877-889

High Resolution Mapping of Interstellar Clouds by Near--IR Scattering

Authors: Paolo Padoan, Mika Juvela, Veli-Matti Pelkonen

Abstract: 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 accu… ▽ More 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. △ Less

Submitted 20 October, 2005; originally announced October 2005.

Comments: Submitted to ApJL

Journal ref: Astrophys.J.636:L101-L104,2006