Dust properties inside molecular clouds from coreshine modeling and observations

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

STAR FORMATION LAWS IN BOTH GALACTIC MASSIVE CLUMPS AND EXTERNAL GALAXIES : EXTENSIVE STUDY WITH DUST CONINUUM, HCN (4-3), AND CS (7-6)

We observed 146 Galactic clumps in HCN (4-3) and CS (7-6) with the Atacama Submillimeter Telescope Experiment 10 m telescope. A tight linear relationship between star formation rate and gas mass traced by dust continuum emission was found for both Galactic clumps and the high redshift (z > 1) star forming galaxies (SFGs), indicating a constant gas depletion time of similar to 100 Myr for molecular gas in both Galactic clumps and high z SFGs. However, low z galaxies do not follow this relation and seem to have a longer global gas depletion time. The correlations between total infrared luminosities (L-TIR) and molecular line luminosities (L-mol') of HCN (4-3) and CS (7-6) are tight and sublinear extending down to clumps with L-TIR similar to 10(3) L-circle dot. These correlations become linear when extended to external galaxies. A bimodal behavior in the L-TIR-L-mol' correlations was found for clumps with different dust temperature, luminosity-to-mass ratio, and sigma(line)/sigma(vir). Such bimodal behavior may be due to evolutionary effects. The slopes of L-TIR-L-mol' correlations become more shallow as clumps evolve. We compared our results with lower J transition lines in Wu et al. (2010). The correlations between clump masses and line luminosities are close to linear for low effective excitation density tracers but become sublinear for high effective excitation density tracers for clumps with L-TIR larger than L-TIR similar to 10(4.5) L-circle dot. High effective excitation density tracers cannot linearly trace the total clump masses, leading to a sublinear correlations for both M-clump-L-mol' and L-TIR-L-mol' relations.Peer reviewe

Synthetic observations of dust emission and polarisation of Galactic cold clumps

Context. The Planck Catalogue of Galactic Cold Clumps (PGCC) contains over 13 000 sources that are detected based on their cold dust signature. They are believed to consist of a mixture of quiescent, pre-stellar, and already star-forming objects within molecular clouds. Aims. 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 properties of this large sample of Galactic objects and, conversely, provides valuable tests for numerical simulations of large volumes of the interstellar medium and the evolution towards pre-stellar cores. Methods. We used several magnetohydrodynamical simulation snapshots to define the density field of our model clouds. Sub-millimetre images of the surface brightness and polarised signal were obtained with radiative transfer calculations. We examined the statistics of synthetic cold clump catalogues extracted with methods similar to the PGCC. We also examined the variations of the polarisation fraction p in the clumps. Results. The clump sizes, aspect ratios, and temperatures in the synthetic catalogue are similar to the PGCC. The fluxes and column densities of synthetic clumps 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 model by a factor of two to three, close to N(H-2) = 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, which is 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 similar to 8% to p similar to 7%). Conclusions. Magnetic field geometry reduces the polarisation fraction in the simulated clumps by only Delta p similar to 1% on average. The larger drop seen towards the actual PGCC clumps therefore 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.Peer reviewe

Chemical exploration of Galactic cold cores

Context. A solar-type system starts from an initial molecular core that acquires organic complexity as it evolves. The so-called prestellar cores that can be studied are rare, which has hampered our understanding of how organic chemistry sets in and grows. Aims. We selected the best prestellar core targets from the cold core catalogue (based on Planck and Herschel observations) that represent a diversity in terms of their environment to explore their chemical complexity: 1390 (in the compressed shell of Lambda Ori), 869 (in the MBM12 cloud), and 4149 (in the California nebula). Methods. We obtained a spectral survey with the IRAM 30 m telescope in order to explore the molecular complexity of the cores. We carried out a radiative transfer analysis of the detected transitions in order to place some constraints on the physical conditions of the cores and on the molecular column densities. We also used the molecular ions in the survey to estimate the cosmic-ray ionisation rate and the S/H initial elemental abundance using a gas-phase chemical model to reproduce their abundances. Results. We found large differences in the molecular complexity (deuteration, complex organic molecules, sulphur, carbon chains, and ions) and compared their chemical properties with a cold core and two prestellar cores. The chemical diversity we found in the three cores seems to be correlated with their chemical evolution: two of them are prestellar (1390 and 4149), and one is in an earlier stage (869). Conclusions. The influence of the environment is likely limited because cold cores are strongly shielded from their surroundings. The high extinction prevents interstellar UV radiation from penetrating deeply into the cores. Higher spatial resolution observations of the cores are therefore needed to constrain the physical structure of the cores, as well as a larger-scale distribution of molecular ions to understand the influence of the environment on their molecular complexity.Peer reviewe

Astrochemical Properties of Planck Cold Clumps

We observed 13 Planck cold clumps with the James Clerk Maxwell Telescope/SCUBA-2 and with the Nobeyama 45 m radio telescope. The N2H+ distribution obtained with the Nobeyama telescope is quite similar to SCUBA-2 dust distribution. The 82 GHz HC3N, 82 GHz CCS, and 94 GHz CCS emission are often distributed differently with respect to the N2H+ emission. The CCS emission, which is known to be abundant in starless molecular cloud cores, is often very clumpy in the observed targets. We made deep single-pointing observations in DNC, (HNC)-C-13, N2D+, and cyclic-C3H2 toward nine clumps. The detection rate of N2D+ is 50%. Furthermore, we observed the NH3 emission toward 15 Planck cold clumps to estimate the kinetic temperature, and confirmed that most targets are cold (less than or similar to 20 K). In two of the starless clumps we observed, the CCS emission is distributed as it surrounds the N2H+ core (chemically evolved gas), which resembles the case of L1544, a prestellar core showing collapse. In addition, we detected both DNC and N2D+. These two clumps are most likely on the verge of star formation. We introduce the chemical evolution factor (CEF) for starless cores to describe the chemical evolutionary stage, and analyze the observed Planck cold clumps.Peer reviewe

PLANCK COLD CLUMPS IN THE lambda ORIONIS COMPLEX. I. DISCOVERY OF AN EXTREMELY YOUNG CLASS 0 PROTOSTELLAR OBJECT AND A PROTO-BROWN DWARF CANDIDATE IN THE BRIGHT-RIMMED CLUMP PGCC G192.32-11.88

We are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the lambda Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound brightrimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity (similar to 0.8 L-circle dot) and accretion rate (6.3 x 10(-7) M-circle dot yr(-1)) when compared with other young Class 0 sources (e.g., PACS Bright Red Sources) in the Orion complex. It has slightly larger internal luminosity (0.21 +/- 0.01 L-circle dot) and outflow velocity (similar to 14 km s(-1)) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than an FHSC. Considering its low internal luminosity (0.08 +/- 0.01 L-circle dot) and accretion rate (2.8 x 10(-8) M-circle dot yr(-1)), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency (similar to 0.3%-0.4%) and core formation efficiency (similar to 1%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed owing to stellar feedback.Peer reviewe

The magnetic field in colliding filaments G202.3+2.5

We observe the magnetic field morphology towards a nearby star-forming filamentary cloud, G202.3+2.5, by the JCMT/POL-2 850 {\mu}m thermal dust polarization observation with an angular resolution of 14.4" (~0.053 pc). The average magnetic field orientation is found to be perpendicular to the filaments while showing different behaviors in the four subregions, suggesting various effects from filaments' collision in these subregions. With the kinematics obtained by N2H+ observation by IRAM, we estimate the plane-of-sky (POS) magnetic field strength by two methods, the classical Davis-Chandrasekhar-Fermi (DCF) method and the angular dispersion function (ADF) method, B_{pos,dcf} and B_{pos,adf} are ~90 {\mu}G and ~53 {\mu}G. We study the relative importance between the gravity (G), magnetic field (B) and turbulence (T) in the four subregions, find G > T > B, G >= T > B, G ~ T > B and T > G > B in the north tail, west trunk, south root and east wing, respectively. In addition, we investigate the projection effect on the DCF and ADF methods based on a similar simulation case and find the 3D magnetic field strength may be underestimated by a factor of ~3 if applying the widely-used statistical B_{pos}-to-B_{3D} factor when using DCF or ADF method, which may further underestimate/overestimate related parameters.Comment: Accepted by ApJ. 20 pages, 9 figure

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP) : Density Structure of Centrally Concentrated Prestellar Cores from Multiscale Observations

Starless cores represent the initial stage of evolution toward (proto)star formation, and a subset of them, known as prestellar cores, with high density (similar to 10(6) cm(-3) or higher) and being centrally concentrated are expected to be embryos of (proto)stars. Determining the density profile of prestellar cores therefore provides an important opportunity to gauge the initial conditions of star formation. In this work, we perform rigorous modeling to estimate the density profiles of three nearly spherical prestellar cores among a sample of five highly dense cores detected by our recent observations. We employed multiscale observational data of the (sub)millimeter dust continuum emission, including those obtained by SCUBA-2 on the James Clerk Maxwell Telescope with a resolution of similar to 5600 au and by multiple Atacama Large Millimeter/submillimeter Array observations with a resolution as high as similar to 480 au. We are able to consistently reproduce the observed multiscale dust continuum images of the cores with a simple prescribed density profile, which bears an inner region of flat density and an r (-2) profile toward the outer region. By utilizing the peak density and the size of the inner flat region as a proxy for the dynamical stage of the cores, we find that the three modeled cores are most likely unstable and prone to collapse. The sizes of the inner flat regions, as compact as similar to 500 au, signify them as being the highly evolved prestellar cores rarely found to date.Peer reviewe

ALMA ACA and Nobeyama Observations of Two Orion Cores in Deuterated Molecular Lines

We mapped two molecular cloud cores in the Orion A cloud with the 7 m Array of the Atacama Compact Array (ACA) of the Atacama Large Millimeter/submillimeterArray (ALMA) and with the Nobeyama 45 m radio telescope. These cores have bright N2D+ emission in single-pointing observations with the Nobeyama 45 m radio telescope, have a relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N2H+ and show narrow line widths of 0.41 km s(-1) and 0.45 km s(-1) in N2D+, respectively, with the Nobeyama 45 m telescope. Both cores were detected with the ALMA ACA 7 m Array in the continuum and molecular lines at Band 6. The starless core G211 shows a clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N2D+ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N2D+ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N2D+ peaks.Peer reviewe