Browsing by Subject "submillimeter: ISM"

Sort by: Order: Results:

Now showing items 1-20 of 21
  • Juvela, Mika; Neha, Sharma; Mannfors, Emma; Saajasto, Mika; Ysard, Nathalie; Pelkonen, Veli-Matti (2020)
    Context. LDN 1642 is a rare example of a star-forming, high-latitude molecular cloud. The dust emission of LDN 1642 has already been studied extensively in the past, but its location also makes it a good target for studies of light scattering.Aims. We wish to study the near-infrared (NIR) light scattering in LDN 1642, its correlation with the cloud structure, and the ability of dust models to simultaneously explain observations of sub-millimetre dust emission, NIR extinction, and NIR scattering.Methods. We used observations made with the HAWK-I instrument to measure the NIR surface brightness and extinction in LDN 1642. These data were compared with Herschel observations of dust emission and, with the help of radiative transfer modelling, with the predictions calculated for different dust models.Results. We find, for LDN 1642, an optical depth ratio tau (250 mu m)/tau (J) approximate to 10(-3), confirming earlier findings of enhanced sub-millimetre emissivity. The relationships between the column density derived from dust emission and the NIR colour excesses are linear and consistent with the shape of the standard NIR extinction curve. The extinction peaks at A(J) = 2.6 mag, and the NIR surface brightness remains correlated with N(H-2) without saturation. Radiative transfer models are able to fit the sub-millimetre data with any of the tested dust models. However, these predict an NIR extinction that is higher and an NIR surface brightness that is lower than based on NIR observations. If the dust sub-millimetre emissivity is rescaled to the observed value of tau (250 mu m)/tau (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 in the direct extinction measurements, which is also reflected in the shape of the NIR surface brightness spectra.Conclusions. 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, which has not been fully explained by the current dust models yet.
  • Juvela, Mika; Guillet, Vincent; Liu, Tie; Ristorcelli, Isabelle; Pelkonen, Veli-Matti; Alina, Dana; Bronfman, Leonardo; Eden, David J.; Kim, Kee Tae; Koch, Patrick M.; Kwon, Woojin; Lee, Chang Won; Malinen, Johanna; Micelotta, Elisabetta; Montillaud, Julien; Rawlings, Mark G.; Sanhueza, Patricio; Soam, Archana; Traficante, Alessio; Ysard, Nathalie; Zhang, Chuan-Peng (2018)
    Context. The sub-millimetre polarisation of dust emission from star-forming clouds carries information on grain properties and on the effects that magnetic fields have on cloud evolution. Aims. Using 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. Methods. JCMT SCUBA-2/POL-2 observations at 850 mu m were combined with Planck 850 mu m (353 GHz) data to map polarisation fraction at small and large scales. With previous total intensity SCUBA-2 observations (450 and 850 mu m) and Herschel data, the column densities were determined via modified black-body fits and via radiative transfer modelling. Models were constructed to examine how the observed polarisation angles and fractions depend on potential magnetic field geometries and grain alignment processes. Results. POL-2 data show clear changes in the magnetic field orientation. These are not in contradiction with the uniform orientation and almost constant polarisation fraction seen by Planck, because of the difference in the beam sizes and the POL-2 data being affected by spatial filtering. The filament has a peak column density of N(H-2) similar to 7 x 10(22) cm(-2), a minimum dust temperature of T similar to 12 K, and a mass of similar to 4300 M-circle dot for the area N(H-2) > 5 x 10(21) cm(-2). The estimated average value of the dust opacity spectral index is beta similar to 1.9. The ratio of sub-millimetre and J-band optical depths is tau (250 mu m)/tau(J) similar to 2.5 x 10(-3), more than four times the typical values for diffuse medium. The polarisation fraction decreases as a function of column density to p similar to 1% in the central filament. Because of noise, the observed decrease of p(N) is significant only at N(H-2) > 2 x 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 similar to 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. Conclusions. The G035.39-00.33 filament shows strong signs of dust evolution and the low polarisation fraction is suggestive of a loss of polarised emission from its densest parts.
  • Juvela, M.; Malinen, J.; Montillaud, J.; Pelkonen, V.-M.; Ristorcelli, I.; Tóth, L. V. (2018)
    Context. The Galactic Cold Cores (GCC) project has made Herschel photometric observations of interstellar clouds where Planck detected compact sources of cold dust emission. The fields are in different environments and stages of star formation. Aims. Our aim is to characterise the structure of the clumps and their parent clouds, and to study the connections between the environment and the formation of gravitationally bound objects. We also examine the accuracy to which the structure of dense clumps can be determined from sub-millimetre data. Methods. We use standard statistical methods to characterise the GCC fields. Individual clumps are extracted using column density thresholding. Based on sub-millimetre measurements, we construct a three-dimensional radiative transfer (RT) model for each field. These are used to estimate the relative radiation field intensities, to probe the clump stability, and to examine the uncertainty of column density estimates. We examine the structural parameters of the clumps, including their radial column density profiles. Results. In the GCC fields, the structure noise follows the relations previously established at larger scales and in lower-density clouds. The fractal dimension has no significant dependence on column density and the values D-p = 1.25 +/- 0.07 are only slightly lower than in typical molecular clouds. The column density probability density functions (PDFs) exhibit large variations, for example, 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 similar to r(-1). In spite of a great variety of clump morphologies (and a typical aspect ratio of 1.5), clumps tend to follow a similar N similar to r(-1) relation below r similar to 0.1 pc. RT calculations indicate only factor 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. Conclusions. The host clouds of the cold clumps in the GCC sample 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 statistical clump properties (e.g. D-p and radial profiles) are insensitive to the environment. The study of clump morphology should be continued with a comparison with numerical simulations.
  • Rivera-Ingraham, A.; Ristorcelli, I.; Juvela, M.; Montillaud, J.; Men'shchikov, A.; Malinen, J.; Pelkonen, V. -M.; Marston, A.; Martin, P. G.; Pagani, L.; Paladini, R.; Paradis, D.; Ysard, N.; Ward-Thompson, D.; Bernard, J. -P.; Marshall, D. J.; Montier, L.; Toth, L. V. (2017)
    Context. The onset of star formation is intimately linked with the presence of massive unstable filamentary structures. These filaments are therefore key for theoretical models that aim to reproduce the observed characteristics of the star formation process in the Galaxy. Aims. As part of the filament study carried out by the Herschel Galactic Cold Cores Key Programme, here we study and discuss the filament properties presented in GCC VII (Paper I) in context with theoretical models of filament formation and evolution. Methods. A conservatively selected sample of filaments located at a distance D <500 pc was extracted from the GCC fields with the getfilaments algorithm. The physical structure of the filaments was quantified according to two main components: the central (Gaussian) region of the filament (core component), and the power-law-like region dominating the filament column density profile at larger radii (wing component). The properties and behaviour of these components relative to the total linear mass density of the filament and the column density of its environment were compared with the predictions from theoretical models describing the evolution of filaments under gravity-dominated conditions. Results. The feasibility of a transition from a subcritical to supercritical state by accretion at any given time is dependent on the combined effect of filament intrinsic properties and environmental conditions. Reasonably self-gravitating (high M-line,M-core) filaments in dense environments (Av greater than or similar to 3 mag) can become supercritical on timescales of t similar to 1 Myr by accreting mass at constant or decreasing width. The trend of increasing M-line,M-tot (M-line,M-core and M-line,M-wing) and ridge A(v) with background for the filament population also indicates that the precursors of star-forming filaments evolve coevally with their environment. The simultaneous increase of environment and filament Av explains the observed association between dense environments and high Mlille,co values, and it argues against filaments remaining in constant single-pressure equilibrium states. The simultaneous growth of filament and background in locations with efficient mass assembly, predicted in numerical models of filaments in collapsing clouds, presents a suitable scenario for the fulfillment of the combined filament mass-environment criterium that is in quantitative agreement with Herschel observations.
  • Abergel, A.; Ade, P. A. R.; Aghanim, N.; Alves, M. I. R.; Aniano, G.; Armitage-Caplan, C.; Arnaud, M.; Ashdown, M.; Atrio-Barandela, F.; Aumont, J.; Baccigalupi, C.; Banday, A. J.; Barreiro, R. B.; Bartlett, J. G.; Battaner, E.; Benabed, K.; Benoit, A.; Benoit-Levy, A.; Bernard, J. -P.; Bersanelli, M.; Bielewicz, P.; Bobin, J.; Bock, J. J.; Bonaldi, A.; Bond, J. R.; Borrill, J.; Bouchet, F. R.; Boulanger, F.; Bridges, M.; Bucher, M.; Burigana, C.; Butler, R. C.; Cardoso, J. -F.; Catalano, A.; Chamballu, A.; Chary, R. -R.; Chiang, H. C.; Chiang, L. -Y; Christensen, P. R.; Church, S.; Clemens, M.; Juvela, M.; Keihänen, Elina; Keskitalo, R.; Kurki-Suonio, H.; Poutanen, T.; Suur-Uski, A. -S.; Tuovinen, J.; Valiviita, J.; Ysard, N. (2014)
  • Planck Collaboration; Aghanim, N.; Keihanen, E.; Kiiveri, K.; Kurki-Suonio, H.; Lindholm, V.; Savelainen, M.; Suur-Uski, A. -S.; Valiviita, J. (2020)
    Observations of the submillimetre emission from Galactic dust, in both total intensity I and polarization, have received tremendous interest thanks to the Planck full-sky maps. In this paper we make use of such full-sky maps of dust polarized emission produced from the third public release of Planck data. As the basis for expanding on astrophysical studies of the polarized thermal emission from Galactic dust, we present full-sky maps of the dust polarization fraction p, polarization angle psi, and dispersion function of polarization angles ?. The joint distribution (one-point statistics) of p and N-H confirms that the mean and maximum polarization fractions decrease with increasing N-H. The uncertainty on the maximum observed polarization fraction, (max) = 22.0(-1.4)(+3.5) p max = 22 . 0 - 1.4 + 3.5 % at 353 GHz and 80 ' resolution, is dominated by the uncertainty on the Galactic emission zero level in total intensity, in particular towards diffuse lines of sight at high Galactic latitudes. Furthermore, the inverse behaviour between p and ? found earlier is seen to be present at high latitudes. This follows the ?proportional to p(-1) relationship expected from models of the polarized sky (including numerical simulations of magnetohydrodynamical turbulence) that include effects from only the topology of the turbulent magnetic field, but otherwise have uniform alignment and dust properties. Thus, the statistical properties of p, psi, and ? for the most part reflect the structure of the Galactic magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map ?xp, looking for residual trends. While the polarization fraction p decreases by a factor of 3-4 between N-H=10(20) cm(-2) and N-H=2x10(22) cm(-2), out of the Galactic plane, this product ?xp only decreases by about 25%. Because ? is independent of the grain alignment efficiency, this demonstrates that the systematic decrease in p with N-H is determined mostly by the magnetic-field structure and not by a drop in grain alignment. This systematic trend is observed both in the diffuse interstellar medium (ISM) and in molecular clouds of the Gould Belt. Second, we look for a dependence of polarization properties on the dust temperature, as we would expect from the radiative alignment torque (RAT) theory. We find no systematic trend of ?xp with the dust temperature T-d, whether in the diffuse ISM or in the molecular clouds of the Gould Belt. In the diffuse ISM, lines of sight with high polarization fraction p and low polarization angle dispersion ? tend, on the contrary, to have colder dust than lines of sight with low p and high ?. We also compare the Planck thermal dust polarization with starlight polarization data in the visible at high Galactic latitudes. The agreement in polarization angles is remarkable, and is consistent with what we expect from the noise and the observed dispersion of polarization angles in the visible on the scale of the Planck beam. The two polarization emission-to-extinction ratios, R-P/p and R-S/V, which primarily characterize dust optical properties, have only a weak dependence on the column density, and converge towards the values previously determined for translucent lines of sight. We also determine an upper limit for the polarization fraction in extinction, p(V)/E(B-V), of 13% at high Galactic latitude, compatible with the polarization fraction p approximate to 20% observed at 353 GHz. Taken together, these results provide strong constraints for models of Galactic dust in diffuse gas.
  • JCMT Large Program SCOPE Collabora; TRAO Key Sci Program TOP Collabora; Yi, Hee-Weon; Lee, Jeong-Eun; Liu, Tie; Kim, Kee-Tae; Choi, Minho; Eden, David; Evans, Neal J.; Di Francesco, James; Fuller, Gary; Hirano, N.; Juvela, Mika; Kang, Sung-ju; Kim, Gwanjeong; Koch, Patrick M.; Lee, Chang Won; Li, Di; Liu, H-Y B.; Liu, Hong-Li; Liu, Sheng-Yuan; Rawlings, Mark G.; Ristorcelli, I.; Sanhueza, Patrico; Soam, Archana; Tatematsu, Ken'ichi; Thompson, Mark; Toth, L.; Wang, Ke; White, Glenn J.; Wu, Yuefang; Yang, Yao-Lun (2018)
    Based on the 850 mu m dust continuum data from SCUBA-2 at James Clerk Maxwell Telescope (JCMT), we compare overall properties of Planck Galactic Cold Clumps (PGCCs) in the lambda Orionis cloud to those of PGCCs in the Orion A and B clouds. The Orion A and B clouds are well-known active star-forming regions, while the A Orionis cloud has a different environment as a consequence of the interaction with a prominent OB association and a giant H-II region. PGCCs in the lambda Orionis cloud have higher dust temperatures (T-d = 16.13 +/- 0.15 K) and lower values of dust emissivity spectral index (beta = 1.65 +/- 0.02) than PGCCs in the Orion A (T-d = 13.79 +/- 0.21 K, beta = 2.07 +/- 0.03) and Orion B (T-d = 13.82 +/- 0.19 K, beta =1.96 +/- 0.02) clouds. We find 119 substructures within the 40 detected PGCCs and identify them as cores. Out of a total of 119 cores, 15 cores are discovered in the lambda Orionis cloud, while 74 and 30 cores are found in the Orion A and B clouds, respectively. The cores in the lambda Orionis cloud show much lower mean values of size R = 0.08 pc, column density N(H-2) (9.5 +/- 1.2) x 10(22)cm(-2) , number density n(H-2) - (2.9 +/- 0.4) x 10 5 CM -3 , and mass M-core = 1.0 +/- 0.3 M(circle dot)compared to the cores in the Orion A [R = 0.11 pc, N(H-2) = (2.3 +/- 0.3) x 10(23) cm(-2), n(H-2) = (3.8 +/- 0.5) x 10(5)cm(-3) , and M-core = 2.4 +/- 0.3 M-circle dot] and Orion B [R = 0.16 pc, N(H-2) (3.8 +/- 0.4) x 10(23) cm(-2), n(H-2) = (15.6 +/- 1.8) x 10(5) cm(-3) , and M-core = 2.7 +/- 0.3 M-circle dot] clouds. These core properties in the A Orionis cloud can be attributed to the photodissociation and external heating by the nearby H rr region, which may prevent the PGCCs from forming gravitationally bound structures and eventually disperse them. These results support the idea of negative stellar feedback on core formation.
  • Ade, P. A. R.; Juvela, M.; Keihanen, E.; Kurki-Suonio, H.; Suur-Uski, A. -S.; Planck Collaboration (2014)
  • Ade, P. A. R.; Keihanen, E.; Kurki-Suonio, H.; Lähteenmäki, Anne; Suur-Uski, A-S; Valiviita, J.; Ysard, N.; Planck Collaboration (2014)
  • Abergel, A.; Juvela, M.; Kurki-Suonio, H.; Lähteenmäki, Anne; Suur-Uski, A. -S.; Valiviita, J.; Planck Collaboration (2014)
  • Arnaud, M.; Atrio-Barandela, F.; Aumont, J.; Baccigalupi, C.; Banday, A. J.; Barreiro, R. B.; Battaner, E.; Benabed, K.; Benoit-Levy, A.; Bernard, J. T.; Bersanelli, M.; Bielewicz, P.; Bonaldi, A.; Bond, J. R.; Borril, J.; Bouchet, E. R.; Bueini, C. S.; Burigana, C.; Cardoso, J. -F; Casassus, S.; Catalano, A.; Cerrigone, L.; Chamballu, A.; Chiang, H. C.; Colombi, S.; Colombo, L. P. L.; Couchot, F.; Crill, B. P.; Curto, A.; Cuttaia, E.; Davies, R. D.; Davis, R. J.; de Bernardis, P.; de Rosa, A.; de Zotti, G.; Delabrouille, J.; Dickinson, C.; Diego, J. M.; Donzelli, S.; Dore, O.; Dupac, X.; Ensslin, T. A.; Eriksen, H. K.; Finelli, E.; Juvela, M.; Keihanen, E.; Keskitalo, R.; Kurki-Suonio, H.; Suur-Uski, A. -S.; Valiviita, J. (2015)
  • Ade, P. A. R.; Aghanim, N.; Arnaud, M.; Ashdown, M.; Atrio-Barandela, F.; Aumont, J.; Baccigalupi, C.; Banday, A. J.; Barreiro, R. B.; Battaner, E.; Benabed, K.; Benoit-Levy, A.; Bernard, J-P.; Bersanelli, M.; Bielewicz, P.; Bobin, J.; Bonaldi, A.; Bond, J. R.; Bouchet, F. R.; Boulanger, F.; Burigana, C.; Cardoso, J-F.; Catalano, A.; Chamballu, A.; Chiang, H. C.; Christensen, P. R.; Clements, D. L.; Colombi, S.; Colombo, L. P. L.; Combet, C.; Couchot, F.; Crill, B. P.; Cuttaia, F.; Danese, L.; Davies, R. D.; Davis, R. J.; de Bernardis, P.; de Rosa, A.; de Zotti, G.; Delabrouille, J.; Dickinson, C.; Diego, J. M.; Donzelli, S.; Dore, O.; Douspis, M.; Keihänen, E.; Kurki-Suonio, H.; Lähteenmäki, A.; Suur-Uski, A-S.; Valiviita, J. (2015)
    Planck data when combined with ancillary data provide a unique opportunity to separate the diffuse emission components of the inner Galaxy. The purpose of the paper is to elucidate the morphology of the various emission components in the strong star-formation region lying inside the solar radius and to clarify the relationship between the various components. The region of the Galactic plane covered is 1 = 300 degrees -> 0 degrees -> 60 degrees where star-formation is highest and the emission is strong enough to make meaningful component separation. The latitude widths in this longitude range lie between 1 and 2, which correspond to FWHM z-widths of 100-200 pc at a typical distance of 6 kpc. The four emission components studied here are synchrotron, free-free, anomalous microwave emission (AME), and thermal (vibrational) dust emission. These components are identified by constructing spectral energy distributions (SEDs) at positions along the Galactic plane using the wide frequency coverage of Planck (28.4-857 GHz) in combination with low-frequency radio data at 0.408-2.3 GHz plus WMAP data at 23-94 GHz, along with far-infrared (FIR) data from COBE-DIRBE and IRAS. The free-free component is determined from radio recombination line (RRL) data. AME is found to be comparable in brightness to the free-free emission on the Galactic plane in the frequency range 20-40 GHz with a width in latitude similar to that of the thermal dust; it comprises 45 +/- 1% of the total 28.4 GHz emission in the longitude range 1 = 300 degrees -> 0 degrees -> 60 degrees. The free-free component is the narrowest, reflecting the fact that it is produced by current star-formation as traced by the narrow distribution of OB stars. It is the dominant emission on the plane between 60 and 100 GHz. RRLs from this ionized gas are used to assess its distance, leading to a free-free z-width of FWHM approximate to 100 pc. The narrow synchrotron component has a low-frequency brightness spectral index beta(synch) approximate to -2.7 that is similar to the broad synchrotron component indicating that they are both populated by the cosmic ray electrons of the same spectral index. The width of this narrow synchrotron component is significantly larger than that of the other three components, suggesting that it is generated in an assembly of older supernova remnants that have expanded to sizes of order 150 pc in 3 x 10(5) yr; pulsars of a similar age have a similar spread in latitude. The thermal dust is identified in the SEDs with average parameters of T-dust = 20.4 +/- 0.4 K, beta(FIR) = 1.94 +/- 0.03 (>353 GHz), and beta(mm) = 1.67 +/- 0.02 (
  • Ade, P. A. R.; Aghanim, N.; Alina, D.; Alves, M. I. R.; Armitage-Caplan, C.; Amaue, M.; Arzoumanian, D.; Ashdown, M.; Atrio-Barandela, F.; Aumont, J.; Baccigalupi, C.; Banda, A. J.; Barreiro, R. B.; Battaner, E.; Benabed, K.; Benoit-Levy, A.; Bernard, J. -P.; Bersanelli, M.; Bielewicz, P.; Bock, J. J.; Bond, J. R.; Borrill, J.; Bouchet, F. R.; Boulanger, F.; Bracco, A.; Burigana, C.; Butler, R. C.; Cardoso, J. -F.; Catalano, A.; Chamballu, A.; Chary, R. -R.; Chiang, H. C.; Christensen, P. R.; Colombi, S.; Colombo, L. P. L.; Combet, C.; Couchot, F.; Coulais, A.; Crill, B. P.; Curto, A.; Cuttaia, E.; Danese, L.; Davies, R. D.; Davis, R. J.; Juvela, M.; Keihänen, Elina; Keskitalo, R.; Kurki-Suonio, H.; Suur-Uski, A. -S.; Valiviita, J. (2015)
    This paper presents an overview of the polarized sky as seen by Planck HFI at 353GHz, which is the most sensitive Planck channel for dust polarization. We construct and analyse maps of dust polarization fraction and polarization angle at 1 degrees resolution, taking into account noise bias and possible systematic effects. The sensitivity of the Planck HFI polarization measurements allows for the first time a mapping of Galactic dust polarized emission on large scales, including low column density regions. We find that the maximum observed dust polarization fraction is high (p(max) = 19.8%), in particular in some regions of moderate hydrogen column density (N-H <2 x 10(21) cm(-2)). The polarization fraction displays a large scatter at NH below a few 10(21) cm(-2). There is a general decrease in the dust polarization fraction with increasing column density above N-H similar or equal to 1 x 10(21) cm(-2) and in particular a sharp drop above N-H similar or equal to 1.5 x 10(22) cm(-2). We characterize the spatial structure of the polarization angle using the angle dispersion function. We find that the polarization angle is ordered over extended areas of several square degrees, separated by filamentary structures of high angle dispersion function. These appear as interfaces where the sky projection of the magnetic field changes abruptly without variations in the column density. The polarization fraction is found to be anti-correlated with the dispersion of polarization angles. These results suggest that, at the resolution of 1 degrees, depolarization is due mainly to fluctuations in the magnetic field orientation along the line of sight, rather than to the loss of grain alignment in shielded regions. We also compare the polarization of thermal dust emission with that of synchrotron measured with Planck, low-frequency radio data, and Faraday rotation measurements toward extragalactic sources. These components bear resemblance along the Galactic plane and in some regions such as the Fan and North Polar Spur regions. The poor match observed in other regions shows, however, that dust, cosmic-ray electrons, and thermal electrons generally sample different parts of the line of sight.
  • Ade, P. A. R.; Aghanim, N.; Alina, D.; Alves, M. I. R.; Aniano, G.; Annitage-Caplan, C.; Arnaud, M.; Arzoumanian, D.; Ashdown, M.; Atrio-Barandela, F.; Aumont, J.; Baccigalupi, C.; Banday, A. J.; Barreiro, R. B.; Battaner, E.; Benabed, K.; Benoit-Levy, A.; Bernard, J. -P.; Bersanelli, M.; Bielewicz, P.; Bond, J. R.; Borrill, J.; Bouchet, F. R.; Boulanger, F.; Bracco, A.; Burigana, C.; Cardoso, J. -F.; Catalano, A.; Chamballu, A.; Chiang, H. C.; Christensen, P. R.; Colombi, S.; Colombo, L. P. E.; Combet, C.; Couchot, F.; Coulais, A.; Crill, B. P.; Curto, A.; Cuttaia, F.; Danese, L.; Davies, R. D.; Davis, R. J.; de Bernardis, P.; de Rosa, A.; Juvela, M.; Keihanen, E.; Kurki-Suonio, H.; Pelkonen, V. -M.; Suur-Uski, A. -S.; Valiviita, J. (2015)
    Polarized emission observed by Planck HFI at 353GHz towards a sample of nearby fields is presented, focusing on the statistics of polarization fractions p and angles psi. The polarization fractions and column densities in these nearby fields are representative of the range of values obtained over the whole sky. We find that: (i) the largest polarization fractions are reached in the most diffuse fields; (ii) the maximum polarization fraction p(max) decreases with column density N-H in the more opaque fields with N-H > 10(21) 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 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(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 maximum 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 S. However, the dispersion of the polarization angle for a given polarization fraction is found to be larger in the simulations than in the observations, suggesting a shortcoming in the physical content of these numerical models. In summary, we find that the turbulent structure of the magnetic field is able to reproduce the main statistical properties of the dust polarization as observed in a variety of nearby clouds, dense cores excluded, and that the large-scale field orientation with respect to the line of sight plays a major role in the quantitative analysis of these statistical properties.
  • Ade, P. A. R.; Aghanim, N.; Alina, D.; Aniano, G.; Armitage-Caplan, C.; Arnaud, M.; Ashdown, M.; Atrio-Barandela, F.; Aumont, J.; Baccigalupi, C.; Banday, A. J.; Barreiro, R. B.; Battaner, E.; Beichman, C.; Benabed, K.; Benoit-Levy, A.; Bernard, J. -P.; Bersanelli, M.; Bielewicz, P.; Bock, J. J.; Bond, J. R.; Borrill, J.; Bouchet, F. R.; Boulanger, F.; Burigana, C.; Cardoso, J. -F.; Catalano, A.; Chamballu, A.; Chary, R. -R.; Chiang, H. C.; Christensen, P. R.; Colombi, S.; Colombo, L. P. L.; Combet, C.; Couchot, F.; Coulais, A.; Crill, B. P.; Curto, A.; Cuttaia, F.; Danese, L.; Davies, R. D.; Davis, R. J.; de Bernardis, P.; de Rosa, A.; de Zotti, G.; Delabrouille, J.; Desert, F. -X.; Dickinson, C.; Diego, J. M.; Donzelli, S.; Dore, O.; Douspis, M.; Dunkley, J.; Dupac, X.; Efstathiou, G.; Ensslin, T. A.; Eriksen, H. K.; Falgarone, E.; Fanciullo, L.; Finelli, F.; Forni, O.; Frailis, M.; Fraisse, A. A.; Franceschi, E.; Galeotta, S.; Ganga, K.; Ghosh, T.; Giard, M.; Giraud-Heraud, Y.; Gonzalez-Nuevo, J.; Gorski, K. M.; Gregorio, A.; Gruppuso, A.; Guillet, V.; Hansen, F. K.; Harrison, D. L.; Helou, G.; Hernandez-Monteagudo, C.; Hildebrandt, S. R.; Hivon, E.; Hobson, M.; Holmes, W. A.; Hornstrup, A.; Huffenberger, K. M.; Jaffe, A. H.; Jaffe, T. R.; Jones, W. C.; Juvela, M.; Keihanen, E.; Keskitalo, R.; Kisner, T. S.; Kneissl, R.; Knoche, J.; Kunz, M.; Kurki-Suonio, H.; Lagache, G.; Lahteenmaki, A.; Lamarre, J. -M.; Lasenby, A.; Lawrence, C. R.; Leonardi, R.; Levrier, F.; Liguori, M.; Lilje, P. B.; Linden-Vornle, M.; Lopez-Caniego, M.; Lubin, P. M.; Macias-Perez, J. F.; Maffei, B.; Magalhaes, A. M.; Maino, D.; Mandolesi, N.; Maris, M.; Marshall, D. J.; Martin, P. G.; Martinez-Gonzalez, E.; Masi, S.; Matarrese, S.; Mazzotta, P.; Melchiorri, A.; Mendes, L.; Mennella, A.; Migliaccio, M.; Miville-Deschenes, M. -A.; Moneti, A.; Montier, L.; Morgante, G.; Mortlock, D.; Munshi, D.; Murphy, J. A.; Naselsky, P.; Nati, F.; Natoli, P.; Netterfield, C. B.; Noviello, F.; Novikov, D.; Novikov, I.; Oxborrow, C. A.; Pagano, L.; Pajot, F.; Paladini, R.; Paoletti, D.; Pasian, F.; Perdereau, O.; Perotto, L.; Perrotta, F.; Piacentini, F.; Piot, M.; Pietrobon, D.; Plaszczynski, S.; Poidevin, F.; Pointecouteau, E.; Polenta, G.; Popa, L.; Pratt, G. W.; Prunet, S.; Puget, J. -L.; Rachen, J. P.; Reach, W. T.; Rebolo, R.; Reinecke, M.; Remazeilles, M.; Renault, C.; Ricciardi, S.; Riller, T.; Ristorcelli, I.; Rocha, G.; Rosset, C.; Roudier, G.; Rusholme, B.; Sandri, M.; Savini, G.; Scott, D.; Spencer, L. D.; Stolyarov, V.; Stompor, R.; Sudiwala, R.; Sutton, D.; Suur-Uski, A. -S.; Sygnet, J. -F.; Tauber, J. A.; Terenzi, L.; Toffolatti, L.; Tomasi, M.; Tristram, M.; Tucci, M.; Umana, G.; Valenziano, L.; Valiviita, J.; Van Tent, B.; Vielva, P.; Villa, F.; Wade, L. A.; Wandelt, B. D.; Zonca, A. (2015)
    The Planck survey provides unprecedented full-sky coverage of the submillimetre polarized emission from Galactic dust. In addition to the information on the direction of the Galactic magnetic field, this also brings new constraints on the properties of dust. The dust grains that emit the radiation seen by Planck in the submillimetre also extinguish and polarize starlight in the visible. Comparison of the polarization of the emission and of the interstellar polarization on selected lines of sight probed by stars provides unique new diagnostics of the emission and light scattering properties of dust, and therefore of the important dust model parameters, composition, size, and shape. Using ancillary catalogues of interstellar polarization and extinction of starlight, we obtain the degree of polarization, p(V), and the optical depth in the V band to the star, tau(V). Toward these stars we measure the submillimetre polarized intensity, P-S, and total intensity, I-S,I- in the Planck 353 GHz channel. We compare the column density measure in the visible, E(B - V), with that inferred from the Planck product map of the submillimetre dust optical depth and compare the polarization direction (position angle) in the visible with that in the submillimetre. For those lines of sight through the di ff use interstellar medium with comparable values of the estimated column density and polarization directions close to orthogonal, we correlate properties in the submillimetre and visible to find two ratios, R-S/V = (P-S/I-S) = (p(V)/tau(V)) and R-P/p = P-S/p(V), the latter focusing directly on the polarization properties of the aligned grain population alone. We find R-S/V = 4.2, with statistical and systematic uncertainties 0.2 and 0.3, respectively, and R-P/p = 5.4 MJy sr(-1), with uncertainties 0.2 and 0.3 MJy sr(-1), respectively. Our estimate of R-S/V is compatible with predictions based on a range of polarizing dust models that have been developed for the di ff use interstellar medium. This estimate provides new empirical validation of many of the common underlying assumptions of the models, but is not yet very discriminating among them. However, our estimate of R-P/p is not compatible with predictions, which are too low by a factor of about 2.5. This more discriminating diagnostic, R-P/p, indicates that changes to the optical properties in the models of the aligned grain population are required. These new diagnostics, together with the spectral dependence in the submillimetre from Planck, will be important for constraining and understanding the full complexity of the grain models, and for interpreting the Planck thermal dust polarization and refinement of the separation of this contamination of the cosmic microwave background.
  • Ade, P. A. R.; Keihanen, E.; Kurki-Suonio, H.; Suur-Uski, A. -S.; Valiviita, J.; Planck Collaboration (2015)
    Planck has mapped the intensity and polarization of the sky at microwave frequencies with unprecedented sensitivity. We use these data to characterize the frequency dependence of dust emission. We make use of the Planck 353 GHz I, Q, and U Stokes maps as dust templates, and cross-correlate them with the Planck and WMAP data at 12 frequencies from 23 to 353 GHz, over circular patches with 10 degrees radius. The cross-correlation analysis is performed for both intensity and polarization data in a consistent manner. The results are corrected for the chance correlation between the templates and the anisotropies of the cosmic microwave background. We use a mask that focuses our analysis on the diffuse interstellar medium at intermediate Galactic latitudes. We determine the spectral indices of dust emission in intensity and polarization between 100 and 353 GHz, for each sky patch. Both indices are found to be remarkably constant over the sky. The mean values, 1.59 +/- 0.02 for polarization and 1.51 +/- 0.01 for intensity, for a mean dust temperature of 19.6 K, are close, but significantly different (3.6 sigma). We determine the mean spectral energy distribution (SED) of the microwave emission, correlated with the 353 GHz dust templates, by averaging the results of the correlation over all sky patches. We find that the mean SED increases for decreasing frequencies at v <60 GHz for both intensity and polarization. The rise of the polarization SED towards low frequencies may be accounted for by a synchrotron component correlated with dust, with no need for any polarization of the anomalous microwave emission. We use a spectral model to separate the synchrotron and dust polarization and to characterize the spectral dependence of the dust polarization fraction. The polarization fraction (p) of the dust emission decreases by (21 +/- 6)% from 353 to 70 GHz. We discuss this result within the context of existing dust models. The decrease in p could indicate differences in polarization efficiency among components of interstellar dust (e.g., carbon versus silicate grains). Our observational results provide inputs to quantify and optimize the separation between Galactic and cosmological polarization.
  • Paradis, D.; Meny, C.; Juvela, M.; Noriega-Crespo, A.; Ristorcelli, I. (2019)
    Context. Some Galactic molecular clouds show signs of dust evolution as compared to the diffuse interstellar medium, most of the time through indirect evidence such as color ratios, increased dust emissivity, or scattering (coreshine). These signs are not a feature of all Galactic clouds. Moreover, molecular clouds in the Large Magellanic Cloud (LMC) have been analyzed in a previous study based on Spitzer and IRIS data, at 4' angular resolution, with the use of one single dust model, and did not show any signs of dust evolution. Aims. In this present analysis we investigate the dust properties associated with the different gas phases (including the ionized phase this time) of the LMC molecular clouds at 1' angular resolution (four times greater than the previous analysis) and with a larger spectral coverage range thanks to Herschel data. We also ensure the robustness of our results in the framework of various dust models. Methods. We performed a decomposition of the dust emission in the infrared (from 3.6 to 500 mu m) associated with the atomic, molecular, and ionized gas phases in the molecular clouds of the LMC. The resulting spectral energy distributions were fitted with four distinct dust models. We then analyzed the model parameters such as the intensity of the radiation field and the relative dust abundances, as well as the slope of the emission spectra at long wavelengths. Results. This work allows dust models to be compared with infrared data in various environments for the first time, which reveals important differences between the models at short wavelengths in terms of data fitting (mainly in the polycyclic aromatic hydrocarbon bands). In addition, this analysis points out distinct results according to the gas phases, such as dust composition directly affecting the dust temperature and the dust emissivity in the submillimeter and different dust emission in the near-infrared (NIR). Conclusions. We observe direct evidence of dust property evolution from the diffuse to the dense medium in a large sample of molecular clouds in the LMC. In addition, the differences in the dust component abundances between the gas phases could indicate different origins of grain formation. We also point out the presence of a NIR-continuum in all gas phases, with an enhancement in the ionized gas. We favor the hypothesis of an additional dust component as the carrier of this continuum.
  • Juvela, Mika (2019)
    Context. Thermal dust emission carries information on physical conditions and dust properties in many astronomical sources. Because observations represent a sum of emission along the line of sight, their interpretation often requires radiative transfer (RT) modelling. Aims. We describe a new RT program, SOC, for computations of dust emission, and examine its performance in simulations of interstellar clouds with external and internal heating. Methods. SOC implements the Monte Carlo RT method as a parallel program for shared-memory computers. It can be used to study dust extinction, scattering, and emission. We tested SOC with realistic cloud models and examined the convergence and noise of the dust-temperature estimates and of the resulting surface-brightness maps. Results. SOC has been demonstrated to produce accurate estimates for dust scattering and for thermal dust emission. It performs well with both CPUs and GPUs, the latter providing a speed-up of processing time by up to an order of magnitude. In the test cases, accelerated lambda iterations (ALIs) improved the convergence rates but was also sensitive to Monte Carlo noise. Run-time refinement of the hierarchical-grid models did not help in reducing the run times required for a given accuracy of solution. The use of a reference field, without ALI, works more robustly, and also allows the run time to be optimised if the number of photon packages is increased only as the iterations progress. Conclusions. The use of GPUs in RT computations should be investigated further.
  • Park, Geumsook; Kim, Kee-Tae; Johnstone, Doug; Kang, Sung-ju; Liu, Tie; Mairs, Steve; Choi, Minho; Lee, Jeong-Eun; Sanhueza, Patricio; Juvela, Mika; Kang, Miju; Eden, David; Soam, Archana; Montillaud, Julien; Fuller, Gary A.; Koch, Patrick M.; Lee, Chang Won; Stamatellos, Dimitris; Rawlings, Jonathan; Kim, Gwanjeong; Zhang, Chuan-Peng; Kwon, Woojin; Yoo, Hyunju (2019)
    In the early stages of star formation, a protostar is deeply embedded in an optically thick envelope such that it is not directly observable. Variations in the protostellar accretion rate, however, will cause luminosity changes that are reprocessed by the surrounding envelope and are observable at submillimeter wavelengths. We searched for submillimeter flux variability toward 12 Planck Galactic Cold Clumps detected by the James Clerk Maxwell Telescope (JCMT)-SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE) survey. These observations were conducted at 850 mu m using the JCMT/SCUBA-2. Each field was observed three times over about 14 months between 2016 April and 2017 June. We applied a relative flux calibration and achieved a calibration uncertainty of similar to 3.6% on average. We identified 136 clumps across 12 fields and detected four sources with flux variations of similar to 30%. For three of these sources, the variations appear to be primarily due to large-scale contamination, leaving one plausible candidate. The flux change of the candidate may be associated with low- or intermediate-mass star formation assuming a distance of 1.5 kpc, although we cannot completely rule out the possibility that it is a random deviation. Further studies with dedicated monitoring would provide a better understanding of the detailed relationship between submillimeter flux and accretion rate variabilities while enhancing the search for variability in star-forming clumps farther away than the Gould Belt.
  • Juvela, Mika; Padoan, Paolo; Ristorcelli, Isabelle; Pelkonen, Veli-Matti (2019)
    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.