Open Access
Issue
A&A
Volume 631, November 2019
Article Number A88
Number of page(s) 18
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201936089
Published online 25 October 2019

© N. Ysard et al. 2019

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Interstellar grains constantly evolve throughout the interstellar medium (ISM) in response to their local environment. This evolution occurs from diffuse regions to protoplanetary discs via dense molecular clouds, prestellar cores, and young stellar objects (YSOs). The term “evolution” can be understood in different ways, namely in terms of chemical evolution (e.g. reactions on the grain surface, accretion of gas phase elements), size evolution (e.g. growth by accretion and/or coagulation, decay by erosion, collisional fragmentation, and photoprocessing), or evolution of the grain structure (e.g. compact, porous, single grains, and aggregates). All these facets of grain evolution are important to probe the ISM chemical, physical, and dynamical properties but in the following we focus only on the impact of grain size evolution on the dust spectral energy distribution (SED) at long wavelengths.

Grain growth most probably starts by accretion and/or coagulation in dense molecular clouds. Dust growth is evidenced by a decrease in temperature, an increase in the far-infrared (FIR) and submillimetre opacity, an increase in the selective extinction RV = AVE(BV), and an increase in the scattering efficiency from the visible to the near-infrared (NIR; e.g. Whittet et al. 2001; Campeggio et al. 2007; Pagani et al. 2010; Schlafly et al. 2016; Remy et al. 2017, 2018). Along with these three phenomena, the spectral index of the emissivity in the FIR and submillimetre increases (e.g. Juvela et al. 2015). Two explanations have been put forward to explain this increase: (i) the physics of the emission of amorphous solids corroborated by laboratory experiments on ISM dust grain analogues (e.g. Agladze et al. 1996; Mennella et al. 1998; Boudet et al. 2005; Meny et al. 2007; Coupeaud et al. 2011; Demyk et al. 2017); and (ii) a change in the chemical composition of the grain surface (Köhler et al. 2015; Jones et al. 2016; Ysard et al. 2016). According to all these observational facts and current dust models, grain sizes are not expected to increase much more than by one order of magnitude from ~0.1 μm in the diffuse ISM to ~1 μm in molecular clouds (e.g. Köhler et al. 2015; Steinacker et al. 2015) along with a structural change from rather compact single grains to fractal or fluffy aggregates.

Denser and more evolved regions exhibit much shallower dust SEDs in the FIR to millimetre wavelength range with spectral indices lower than in molecular clouds and even lower than in the diffuse ISM. This was observed in the first hydrostatic core candidates (Young et al. 2018) and YSOs (Chiang et al. 2012; Miettinen et al. 2012; Miotello et al. 2014; Choi et al. 2017; Gerin et al. 2017; Sheehan & Eisner 2017) with spectral index values below 1.5. Observations of protoplanetary discs also show very low spectral index values (e.g. Isella et al. 2010; Kwon et al. 2015; Pérez et al. 2015; Ribas et al. 2017; Tripathi et al. 2018), which decrease from the outer to inner disc regions (e.g. Guilloteau et al. 2011; Pérez et al. 2012; Trotta et al. 2013; Tazzari et al. 2016). Spectral index values lower than unity have been attributed to the growth of dust particles to radii of around 1 mm to 1 cm in size (e.g. Natta & Testi 2004; Natta et al. 2007; Draine 2006; Birnstiel et al. 2010; Liu et al. 2017). The scheme is that after an initial fractal growth of grains to micronic aggregates, grains continue to grow andare compacted by grain–grain collisions, ram pressure of the gas, and self-gravity for the largest bodies of ~ 1 km in size (e.g. Blum & Wurm 2008; Blum 2018; Wada et al. 2008; Paszun & Dominik 2009; Dominik et al. 2016) and drift towards disccentres. The grain chemical composition is also expected to vary inside protoplanetary discs from the mid-plane to the surface and from the inner to outer regions. Most grain growth studies have assumed the same grain composition at all radii and heights (e.g. Pérez et al. 2012; Tazzari et al. 2016). Therefore, the general questions arise of up to what sizes can dust grains grow dependent on local conditions within molecular clouds, prestellar cores, YSOs, and protoplanetary discs, and how the derived sizes are affected by changes in the dust grain chemical composition. In order to answer these questions, itis necessary to carry out detailed model calculations of dust emission at multiple wavelengths considering realistic dust properties with a wide range of grain compositions, and to compare these theoretical predictions with observations. Instruments such as the Atacama Large Millimeter/submillimeter Array (ALMA) and, in the future, Square Kilometre Array (SKA) will allow us to measure the distribution of centimetre-sized grains down to AU spatial scales, opening the possibility to test grain growth and any dependence on grain composition in dense interstellar regions.

After exploring the effect of complex grain structures with a single chemical composition in Ysard et al. (2018), the intent of this study isto investigate the effect of complex chemical compositions for isolated spherical grains up to very large sizes. This rather simple approach may not be too unrealistic in the case of very dense environments where grain–grain compaction and gas compression may be efficient processes (e.g. Kataoka et al. 2013). In any case, this study could be used as a basis to compare the influence of composition verses structure complexity on the dust optical property variations.

This paper is organised as follows. Section 2 presents the dust chemical compositions considered in this study and the methods used to first translate them into optical properties and then into dust emission spectra. Section 3 shows the resulting absorption and scattering efficiencies along with the corresponding mass absorption coefficients. Section 4 presents the associated emission spectra in the case of optically thin (Sect. 4.1) and optically thick regions (Sect. 4.2) in addition to the associated FIR to centimetre spectral indices (Sect. 4.3). Finally, Sect. 5 summarises our results.

2 Dust properties and methods

Dust emission depends on grain optical properties and size distribution, and on the strength and intensity of the radiation field. In this section, we describe the variations in dust composition considered in this study and how we translate them into optical properties. The methods used to derive the corresponding dust emission are also outlined.

2.1 Dust optical properties

Our starting point is the optical constants (i.e. complex refractive indices m = n + ik) from the dust modelling framework THEMIS1 (The Heterogeneous dust Evolution Model for Interstellar Solids), which is briefly summarised in Jones et al. (2017)2. Two main sets of optical constants are included in THEMIS: amorphous magnesium-rich silicates with metallic iron and iron sulphide nano-inclusions and amorphous semi-conducting hydrocarbon grains. For the amorphous silicate component, THEMIS consists of two chemical compositions: one similar to that of enstatite and another similar to that of forsterite (Scott & Duley 1996). As the differences are small when moving away from the 10 and 18 μm silicate spectral features, we only consider forsterite normative composition in the following and refer to it as a-Sil. For the carbonaceous component, we consider two extreme cases in terms of hydrogen content: aromatic-rich grains with a hydrogen fraction XH ~ 0.02, referred to as a-C, and aliphatic-rich grains with XH ~ 0.58, referred to as a-C:H. In addition to these three materials, we also consider water ice with the optical constants given by Warren (1984).

Starting from the four aforementioned materials, we consider several composition mixtures and grain structures. For the sake of comparison, we first consider compact grains of purely a-Sil, a-C, or a-C:H. Subsequently, according to Köhler et al. (2015), we consider compact grains made of two thirds a-Sil and one third a-C (Mix 1) or one third a-C:H (Mix 2), in terms of volume fractions. These allow reproduction of the mass fractions derived by Jones et al. (2013) for the diffuse ISM. The effect of porosity is tested for the Mix 1 mixture, with a porosity degree of 50% (Mix 1:50). We also evaluate theimpact of the presence of a water ice mantle on compact Mix 1 grains (Mix 1:ice). We further consider two material compositions defined in Pollack et al. (1994) based on depletion measurements: (i) 21% a-Sil and 79% a-C (Mix 3); and (ii) 8% a-Sil, 30% a-C, and 62% water ice (Mix 3:ice). The various grain compositions are summarised in Table 1. For each grain composition, we derive the absorption and scattering efficiencies Qabs and Qsca, respectively, and the asymmetry factor of the phase function g = ⟨cosθ⟩. To allow fast calculations, we make the major assumption that the grains are spherical and compute their optical properties using the Mie theory (Mie 1908; Bohren & Huffman 1983) with the Fortran 90 version of the BHMIE routine given in Bohren & Huffman (1998). For grains consisting of two or three materials, we first derive effective optical constants following the Maxwell Garnett mixing rule (Maxwell Garnett 1904; Bohren & Huffman 1998). Indeed, we assume that in Mix 1 grains, for example, carbon appears as proper inclusions in the silicate matrix rather than assuming a completely random inhomogeneous medium. Mishchenko et al. (2016a,b) performed exhaustive studies of the applicability of the Maxwell Garnett mixing rule to heterogeneous particles. These latter authors showed that this rule can provide accurate estimates of the scattering matrix and absorption cross-section of heterogeneous grains at short wavelengths (typically up to the visible for a 0.1 μm grain and to the mid-infrared(MIR) for a 10 μm grain) if twocriteria are met: both the size parameter of the inclusions and the refractive index contrast between the host material and the inclusions have to be small. Moreover, Mishchenko et al. (2016a) demonstrated that the extinction and asymmetry-parameter errors of the Maxwell Garnett mixing rule are significantly smaller than the scattering-matrix errors, remaining small enough for most typical applications and in particular the kind of applications we perform here. It is however well known that this kind of mixing rule systematically underestimates the absorption efficiency in the FIR to millimetre wavelength range, the implications of which are discussed in Sect. 3.2. We perform our computations with the emc routine of V. Ossenkopf3. For Mix 1 and Mix 2, we assume a matrix of a-Sil with inclusions of a-C or a-C:H, and for Mix 3 a matrix of a-C with inclusions of a-Sil. For grains surrounded by an ice mantle, the optical properties are derived with the core-mantle Mie theory using the BHCOAT routine given in Bohren & Huffman (1998).

Table 1

Grain compositions, structures, and volume densities in g cm−3 as described in Sect. 2.1.

2.2 Dust size distribution

In all our calculations, the grain minimum and maximum sizes are 0.01 μm and 10 cm. Following Köhler et al. (2015), we assume a log-normal size distribution: dndaexp[12(ln(a/a0)σ)2],\begin{equation*}\frac{\textrm{d}n}{\textrm{d}a} \propto \exp\left[ -\frac{1}{2} \left( \frac{\ln(a/a_0)}{\sigma} \right)^2 \right], \end{equation*}(1)

where a is the grain radius. The width σ = 0.7 is fixed but we allow variations inthe centroid a0 from 0.1 μm to 1 cm to mimic the effect of grain growth. The gas-to-dust mass ratio is fixed4 to 100.

The choice of the grain size distribution indeed has a significant impact on both the emission and extinction. However, neither observationsnor laboratory experiments have yet made it possible to conclusively determine its shape. Many studies rely on power-law size distributions (dn∕daap), which are consistent with grains resulting from collisional fragmentation cascades (e.g. Davis & Ryan 1990; Tanaka et al. 1996). This is what is typically expected for debris discs for instance. However, recent studies on grain growth during core collapse and in protoplanetary discs have shown that their size distributions are expected to deviate substantially from power laws (e.g. Birnstiel et al. 2011, 2018). In particular, based on experimental studies of collisional growth (Weidling et al. 2009; Güttler et al. 2009, 2010; Windmark et al. 2012; Gundlach & Blum 2015), the model developed by Lorek et al. (2018) for local growth in protoplanetary discs leads to size distributions close to log-normal functions (see their Figs. 2 and 3). These theoretical and experimental results lead us to choose log-normal distributions rather than the more common power-laws which we do however use for comparison purposes several times in this study.

2.3 Dust emission

Two kinds of calculations are required to arrive at the dust emission spectra, which depend on whether the medium considered is optically thin or thick. In the case of an optically thin medium, the dust emission spectrum is computed with the DustEM5 code described in Compiègne et al. (2011). DustEM is a numerical tool, which computes the dust emission and extinction as a function of the grain size distribution and optical properties for a given radiation field.

In the case of an optically thick medium, we use the 3D Monte Carlo radiative transfer model CRT (Continuum Radiative Transfer, Juvela & Padoan 2003; Juvela 2005), which is coupled to DustEM (Ysard et al. 2012). We first assume a filamentary cloud externally heated by the standard interstellar radiation field (ISRF) as defined by Mathis et al. (1983) scaled by the G0 factor. Unless otherwise stated G0 = 1. The cloud is represented by an infinite right circular cylinder with constant density along its axis and the following radial density distribution: n(r)=nC1+(r/H0)q,\begin{equation*}n(r) \,{=}\, \frac{n_{\textrm{C}}}{1 + (r/H_0)^q}, \end{equation*}(2)

where nC is the central density and H0 the internal flat radius. We define H0 for the cloud to be at equilibrium with its mass per unit length equal to the critical mass defined by Ostriker (1964): Mcrit=2cS2/G$M_{\textrm{crit}} \,{=}\, 2c_{\textrm{S}}^2/G$ with G being the gravitational constant. The sound speed, cS, is assumed to be constant and computed for a gas temperature of 12 K. We fix the steepness of the profile to q = 2 (Arzoumanian et al. 2011). Second, and in order to assess the effect of internal heating, we consider a spherical cloud with the radial density distribution defined in Eq. (2) and a central blackbody radiation source of one solar luminosity and a temperature of 3000 K, which is consistent with theoretical predictions of pre-main sequence evolution (see Fig. 3 in Wuchterl & Tscharnuter 2003). The cloud is also externally illuminated by the ISRF.

These calculations provide dust temperatures and emission spectra. From the dust emission spectra, we further extract the spectral indices β = α − 2, where α is the slope of the flux density obtained with linear regressions on the log-log flux density versus frequency plots over 1 GHz frequency steps.

3 Results: optical properties

For the grain chemical compositions detailed in Table 1, we now present the corresponding optical properties for sizes between 0.01 μm and 10 cm as absorption and scattering efficiencies (Sect. 3.1). To illustrate the importance of the size distribution choice, we also show mass absorption coefficients for various log-normal and power-law distributions (Sect. 3.2).

3.1 Scattering and absorption efficiencies

The absorption efficiencies for spherical grains with different radii and material compositions are shown in Fig. 1. For small grains (a ≲ 1 μm) consisting of a-Sil and a-C:H, features typical of these materials are visible in the near- to mid-IR. The material mixtures including a-Sil also show the silicate features, except for Mix 3, in which the a-Sil content is too low for the 10 and 18 μm features to be seen. For grains with ice mantles, the ice spectral features at ~ 3, 13, and 50 μm are visible above the continuum. For larger grains (1 μm ≲ a ≲ 1 mm), the spectral features gradually disappear as the continuum to feature contrast increases. Indeed, the threshold wavelength – under which the absorption and scattering efficiencies tend asymptotically towards one, and above which they sharply decrease – increases with the grain size, its exact spectral position depending on the material conducting properties as discussed in detail in Ysard et al. (2018) for the THEMIS materials.

For grains larger than 1 mm, both the absorption efficiency Qabs and scattering efficiency Qsca are almost constant and equal to one from UV to centimetre wavelengths. This results in relatively flat absorption and scattering efficiencies, only slightly dependent on wavelength. For smaller grains, at wavelengths much larger than the threshold wavelength, the slope of Qabs does not depend on the size but only on the material and is equal to ~2, 1.3, 2, 1.35, 2, 1.3, 1.35, 1.3, and 1.3 for a-Sil, a-C, a-C:H, Mix 1, Mix 2, Mix 1:50, Mix 1:ice, Mix 3, and Mix 3:ice compositions, respectively. Two interesting points may be noted. Firstly, all grains containing a-C, regardless of the proportion, have almost the same slope as pure a-C grains. Indeed, as a-C grains are low band gap semi-conductors (Eg = 0.1 eV for XH = 0.02), free-carrier absorption is substantial resulting in longer threshold wavelengths than in the case of a-Sil grains. Secondly, neither porosity nor the presence of ice mantles influence the slope steepness.

3.2 Massabsorption coefficient

To infer dust masses from FIR to millimetre observations, observers often use a fixed mass absorption coefficient, κ, at a given wavelength: see for instance Planck Collaboration XXII (2011) and Busquet et al. (2019) using the κ value of Beckwith et al. (1990) to derive masses of cold molecular cores and protostellar discs, respectively. The left panel of Fig. 2 shows the κ values at 1.3 mm for the dust compositions presented in Table 1 and log-normal size distributions centred at 0.1 μm ≤ a0 ≤ 1 cm (see Eq. (1)). As expected, κ is highly dependent on grain composition with variations up to a factor of approximately three, with the exception of grains of single composition and those containing a significant volume fraction of a-C:H (i.e. a-Sil, a-C, a-C:H, and Mix 2). However, even greater variations are observed depending on the grain size distribution for a0 ≥ 510 μm (κ is almost constant otherwise). Excluding grains of single composition or grains containing a-C:H, for a0 = 1 to ~ 100 μm, κ increases by factors of ~ 7 to 19. For size distributions centred at larger sizes, κ decreases again to κ[a0 = 1 cm]~0.1 − 0.25 × κ[a0 = 0.1 μm]. This is simply explained by the fact that log-normal size distributions have a characteristic grain size dominating the total dust mass and by the definition of the mass absorption coefficient which, at zero order, can be approximated by κ = (3∕4ρ) × (Qabsa0). For a0 ≲ 5–10 μm depending on the grain composition, a0 is small compared to the wavelength and the increase in Qabs(λa0) is proportional to a0 resulting in a constant κ. For 5 10 ≲ a0 ≲ 60200 μm, the grain sizes become comparable to the wavelength, resulting in a stronger increase in Qabs and therefore in an increase in κ. For larger grain sizes, at λ = 1.3 mm, Qabs ~ 1 is almost independent of size for all dust compositions (see Fig. 1) resulting in a decreasing κ for increasing sizes (κ ∝ 1∕a0).

For the sake of comparison, the middle and right panels of Fig. 2 again show the κ values at 1.3 mm but in the case of power-law size distributions: dn∕daap. The middlepanel shows the case of a variable power-law exponent − 4≤ p ≤−2 for constant size limits 0.01 μm ≤ a ≤ 10 cm, whereas the right panel shows the case of a constant exponent p= − 3.5 with an increasing maximum grain size 0.01 μm ≤ aamax = 0.1 μm to 10 cm. Indeed, for the most part, previous studies have used such size distributions, usually with p = − 3.5 (e.g. Natta & Testi 2004; Draine 2006; Isella et al. 2012; Liu et al. 2017). For p < −2.25 or amax ≥ 1 μm and whatever the grain composition, the mass absorption coefficient is constant and roughly equal to the mass absorption coefficients found for log-normal size distributions centred on a0 ≤ 5 μm. An increase in κ of about a factor of 1.3 to 3 is found only when p ≥−2.5 or amax < 1 μm, far from the highest values found for the log-normal case at a0 ~ 60 − 100 μm. Tables giving the mass absorption coefficient values κ for wavelengths from 250 μm to 2.8 mm are available in Appendix A.

As already stated in the introduction, several laboratory experiments on ISM grain analogues have also shown that their mass absorption coefficient depends not only on their size and composition but also on their temperature(e.g. Agladze et al. 1996; Mennella et al. 1998; Boudet et al. 2005; Meny et al. 2007; Coupeaud et al. 2011; Demyk et al. 2017). For instance, Demyk et al. (2017) found that for Mg-rich amorphous silicate analogues, κ1mm doubles when the temperature increases from 10 to 300 K, while the spectral index at submillimetre wavelengths decreases from about 2.2 to1.6. Even if the inclusion of such effects is beyond the scope of our study, it should be kept in mind that when observing regions encompassing large temperature gradients, intrinsic variations in the mass absorption coefficient of a given material can introduce significant uncertainties on mass determinations based solely on the long-wavelength FIR to millimetre dust emission.

Further, when calculating the grain optical properties, we make the major assumption that these grains are spherical. However, it is well known that their exact shape significantly alters the scattering and absorption efficiencies and therefore their mass absorption coefficient at all wavelengths. This was shown in studies based on interstellar dust, cometary dust, and terrestrial aerosols. This alteration depends in particular on the size and shape of the monomers composing the aggregate grain (e.g. Liu et al. 2015; Wu et al. 2016; Min et al. 2016; Ysard et al. 2018), on the monomer composition and composition mixing (e.g. Köhler et al. 2012; Min et al. 2016), and on the contact area between the monomers (Köhler et al. 2011). All these studies have focused on different points, but compared to the simple case of the porous or non-porous sphere, a common observation can be made: the influence of the exact shape of the aggregates can increase the FIR to millimetre mass absorption coefficient by a few tens of percent and give a lower spectral index in the same wavelength range6.

To conclude, Fig. 2 and previous laboratory and theoretical studies clearly illustrate how the choice of a mass absorption coefficient κ to infer dust masses from observations is very model dependent. Thus, inter-comparisons between the dust masses derived from different studies should always be taken with caution and, in particular, it should be noted that the need to adopt a dust model automatically introduces at least one order of magnitude uncertainty into any dust mass determination.

thumbnail Fig. 1

Absorption (Qabs, solid lines) and scattering efficiencies (Qsca, dashed lines) for the dust compositions described in Sect. 2.1 and Table 1. Blue lines show the case of grains with a radius of 0.01 μm, orange of 0.1 μm, pink of 1 μm, green of 10 μm, red of 100 μm, black of 1 mm, cyan of 1 cm, and brown of 10 cm.

thumbnail Fig. 2

Left: mass absorption coefficients κ at 1.3 mm for log-normal size distributions centred on a0 and for the dust compositions described in Sect. 2.1 and Table 1. Black, blue, pink, green, brown, orange, red, grey, and magenta lines show the cases of a-Sil, a-C, a-C:H, Mix 1, Mix 2, Mix 1:50, Mix 1:ice, Mix 3, and Mix 3:ice grains, respectively. Mass absorption coefficients of pure a-C:H grains are always lower than 10−4 cm2 g−1 and thus not visible on the figures. The black dashed line shows κ = 1.8 cm2 g−1, a value close to those of the models commonly used in the literature with a gas-to-dust ratio of 100 (e.g. Beckwith et al. 1990; D’Alessio et al. 2001; Andrews et al. 2009; Birnstiel et al. 2018). Middle: same for power-law size distributions d n∕d aap with variable exponent p and 0.01 μm ≤ a ≤ 10 cm. Right: same for power-law size distributions with p= − 3.5, amin = 0.01 μm, and variable amax.

4 Results: emission spectra

The purpose of this section is to provide numerical values comparable to those derived from astronomical observations. From the optical properties and size distributions presented in the previous section, we therefore calculate SEDs for different radiation fields and their associated spectral indices.

4.1 Optically thin medium

Using the DustEM code, the equilibrium temperatures of the grains presented in Table 1 and their SEDs as flux densities have been calculated and are shown in Figs. 3 and 4, respectively. As described in Sect. 2, the temperatures and spectra are computed for the standard ISRF with G0 = 1 for various values of the size distribution centroid: a0 = 0.1 μm, 1 μm, 10 μm, 100 μm, 1 mm, and 1 cm.

First of all, pure a-C:H grains have a far different behaviour from the other grain types. Indeed, Fig. 1 shows that their optical properties are quite different with an absorption efficiency always inferior to the scattering efficiency, regardless of the grain size. This makes these grains poorly emissive which results in high grain equilibrium temperatures for all sizes (Fig. 3) and flux densities peaking at relatively short wavelengths (Figs. 4 and 5). Although this behaviour is interesting from a theoretical point of view, we do not discuss this case further since such grains are improbable in the ISM because we expect mixing with silicate or ice components and/or partial dehydrogenation due to UV photons7.

Figure 3 shows the equilibrium temperature as a function of grain size for the grain compositions described in Table 1. At small sizes, the temperature depends on the material through its heat capacity and the size-dependent variations in Qabs. On the other hand, at larger sizes (a ≳ 100 μm), the temperature is practically independent of the material, which is explained by the fact that all materials, except pure a-C:H, have Qabs ~ Qsca ~ 1. Excluding pure a-C:H, the temperature of the largest grains does not vary by more than 1 K from one material to another. A last interesting point regarding the equilibrium temperature is the fact that for pure a-C grains and grains made partly of a-C, the equilibrium temperature slightly increases by about 0.6 K from a = 100 μm to 1 cm. This is the consequence of the spectral shape of the a-C Qabs, which instead of being almost constant for the largest sizes, decreases slightly up to the threshold wavelength, which makes the largest grains slightly less emissive and thus marginally hotter.

The dependence of temperature on size and material is reflected in the flux density. For grains without a-C (a-Sil, a-C:H, and Mix 2), the wavelength of the SED peak increases with the centroid size a0, whereas for grains containing a-C (i.e. all except a-Sil, a-C:H, and Mix 2), it first increases before decreasing slightly or stagnating when a0 ≳ 100 μm or 1 mm depending on the grain composition (Fig. 5). Figure 4 also shows that the slope of the flux density in the FIR and submillimetre strongly depends on both the grain size and composition with steeper slopes for the smallest a0 values. This is easily visible in Fig. 6 where β = α − 2 is plotted for λ ≥ 200 μm. For the three smallest centroid values, a0 = 0.1, 1, and 10 μm, β(λ ≳ 1 cm) is asymptotically equal to the slope of the material absorption efficiency Qabs (see Table 1). For size distributions centred on sizes a0 larger than 10 μm, the spectral indices are globally lower as long as λ ≲ 10a0. The break in the dependence of the equilibrium temperature on the grain size for a ~ 100 μm (see Fig. 3) directly impacts the spectral index. Compared to the smaller ones, the largest grains have slightly higher temperatures when they contain a-C or are of a comparable temperature otherwise. Their contribution is enhanced when a0 increases which distorts the spectral shape of the flux density peak compared to smaller a0 cases. This results in low β values increasing with wavelength. Based solely on (sub)millimetre/centimetre β measurements, for a0 smaller than 100 μm, it would be very difficult to determine the grain sizes for most of the grain compositions considered in this study. In Sect. 4.2, we investigate to what extent this statement holds in optically thick media in which dust growth is expected to occur.

In addition, it should be noted that the negative spectral index values at short wavelengths are artificial and result from the application of the Rayleigh-Jeans approximation outside its validity domain, that is, too far ahead or too close to the SED peak (see Fig. 5). In order to fully characterise the SEDs, Fig. 6 also shows the spectral indices obtained by fitting the Wien part of the SEDs with a modified blackbody: IνBν(T)νβ, where Bν (T) is the Planck function at temperature T and β is the SED spectral index. The wavelength range corresponding to the Wien part of the SEDs depends on both the dust composition and size distribution and is therefore different for all the cases presented in Fig. 6. Our choice is to perform the fits on wavelength ranges matching the full width at half maximum of each SED. Positive spectral indices are found for the smallest sizes and null for the largest (a0 = 1 mm and 1 cm). It also appears that for a0 ≤ 10 μm the spectral indices of the modified blackbodies are close to the material absorption efficiencies Qabs slopes (see Table 1).

For the sake of comparison, Fig. 4 shows the three extreme cases of power-law size distributions already presented in Sect. 3.2: (i) p =−3.5, 0.01 μm ≤ a ≤ 10 cm (solid); (ii) p = −2, 0.01 μm ≤ a ≤ 10 cm (dashed); (iii) p =−3.5, 0.01 μm ≤ a ≤ 1 μm (dotted). This comparison leads to two main results. Firstly, case (ii), in which a large part of the mass is distributed in the largest grains, gives a SED similar to the case where a0 = 1 cm. Secondly, case (iii), which places most of the mass in the smallest grains, is similar to the case where a0 = 0.1 μm. Both the peak of the SEDs and their long wavelength spectral indices are similar. Case (i) is intermediate where no size dominates the grain size distribution. This results in a very broad SED peaking at about the same wavelength as when a0 = 0.1 μm but its centimetre emission in terms of spectral index is close to that of the case where a0 = 1 mm.

Figure 7 shows the influence of the intensity of the radiation field on the SED for a0 = 0.1 μm and 1 cm. For a0 = 0.1 μm, the peak of the SED shifts to longer wavelengths when G0 decreases. Thus, the spectral index reaches its asymptotic value, equal to the slope of the material absorption efficiency Qabs, at shorter wavelengths for higher G0. For a0 = 1 cm, the peak of the SED also shifts to longer wavelengths when G0 decreases provided that G0 ≥ 0.1. For lower G0, there are no longer enough energetic photons available to warm up the very large grains which make most of the mass in the size distribution and these grains end up being at equilibrium with the cosmic microwave background (CMB) temperature. This break starts appearing for a0 ≳ 1 mm and G0 ≲ 10−2.

thumbnail Fig. 3

Equilibrium temperature as a function of grain size for the dust compositions described in Sect. 2.1 and Table 1, illuminated by the ISRF with G0 = 1. The same colour code is used as in Fig. 2.

thumbnail Fig. 4

Spectral energy distributions for the dust compositions described in Sect. 2.1 and Table 1. The SEDs were obtained considering an optically thin medium with a column density NH = 1020 H cm−2 illuminated by the standard ISRF (see Sect. 2.3 for details). The orange lines show the case where the centroid of the log-normal size distribution is a0 = 0.1 μm, the pink lines where a0 = 1 μm, the green lines where a0 = 10 μm, the red lines where a0 = 100 μm, the black lines where a0 = 1 mm, and the cyanlines where a0 = 1 cm. The thick grey curves on the Mix 1:50 plot present the case of power-law size distributions with amin = 0.01 μm and where the dotted line shows the case of p =−3.5 and amax = 1 μm, the solid lineof p =−3.5 and amax = 10 cm, and the dashed line of p = −2 and amax = 10 cm.

thumbnail Fig. 5

Wavelength of the maxima of the SEDs shown in Fig. 4 as a function of the size distribution centroid a0. The same colour code is used as in Fig. 2.

4.2 Optically thick medium

As described in Sect. 2.3, we compute the emission spectra at the centre of an optically thick cloud (see Eq. (2)), either externally illuminated by the ISRF alone or with an additional central heating source. The resulting SEDs for various size distribution centroids 0.1 μm ≤ a0 ≤ 1 cm and nC = 105 H cm−3 are presented in Fig. 8.

Compared to the optically thin case, the SEDs of the dense clouds illuminated by the ISRF alone, regardless of grain composition, are slightly redshifted (solid lines in Fig. 8). The shift in wavelength decreases with increasing a0 values. As a result, their spectral indices are somewhat lower at short wavelengths but tend asymptotically towards the same value at long wavelengths, that of the grain absorption efficiency Qabs (Figs. 1 and 9). The fact that there is little difference between the optically thin and optically thick cases shows that the SEDs of the latter are dominated by the outer layers of the clouds where the grains are the warmest(e.g. Ysard et al. 2012, and references therein). In prestellar cores, the density distribution can be steeper than what we consider here (e.g. Tafalla et al. 2002, 2004). To test this, we ran calculations with q = 3 and nC = 105 H cm−3 (see Eq. (2)) and found that the spectral index does not vary. The only difference is that the SED is slightly lower than in the q = 2 case due to a lower column density.

As expected, grains heated by an additional internal source are hotter, hence the shift of the SED peak positions towards shorter wavelengths (dashed lines in Fig. 8). The SEDs are also broader than for clouds illuminated only externally by the ISRF. In this latter case, the SEDs are dominated by the outer layers of the clouds, whereas for internally heated clouds, the inner parts are hot enough to contribute significantly to the total emitted flux resulting in a broader dust temperature distribution and thus in wider SEDs. It can be seen in Fig. 8 that the SEDs are composed of two components: a “cold” component that corresponds to the cloud outer layers and very closely resembles the ISRF-only case, and a “hot” component responsible for the shorter wavelength peak corresponding to the inner layers of the heated cloud. For compositions including water ice mantles (Mix 1:ice and Mix 3:ice, see bottom left and right plots in Fig. 8), the SEDs appear distorted at the peak for a0 = 0.1 and 1 μm. This distortion is due to the water ice bands at ~13 and 50 μm that become visible because these small grains are relatively hot. Again, the spectral indices for all size distributions and compositions tend asymptotically towards that of their material absorption efficiencies at long wavelengths (Fig. 9). Yet, since the SEDs peak at shorter wavelengths and are wider, there are two significant differences between the case with internal source compared to the case without: (i) the spectral index asymptotic value is reached at shorter wavelengths – for example λ ~ 1 mm instead of 1 cm for a0 ≤ 1 μm; (ii) for the largest grains, the spectral index is systematically higher with β = α − 2 > −1 for λ > 1 mm. As for the case of clouds without an internal source, we tested the influence of changing the density distribution steepness q by calculating the SEDs for q = 3 and nC = 105 H/cm3 (see Eq. (2)). Since the density decreases faster from the inside out, the grains are globally hotter. The difference in the SED peak position is negligible for a0 ≤ 10 μm and of the order of 100 μm for larger a0. The resulting spectral indices are thus strictly similar to those of the q = 2 case for small a0 values at all wavelengths and at λ ≳ 1 mm for larger a0.

In addition, Fig. 9 also shows the spectral indices obtained by fitting the Wien part of the SEDs with a modified blackbody (see Sect. 4.1 for details). Compared to the optically thin case, the results differ significantly only when an internal source is considered. Due to the broader temperature distribution along the line of sight, the spectral indices are lower for all size distributions and even become negative for the largest log-normal centroids for most dust compositions. For the two compositions with water ice mantles (Mix 1:ice and Mix 3:ice), the distortion of the SED peaks due to the 13 and 50 μm ice bands produces higher spectral indices for the two smallest values of a0 : β(a0 = 0.1 μm) = 5.6 and β(a0 = 1 μm) = 3.6 in the Mix 1:ice case, and β(a0 = 0.1 μm) = 7.0 and β(a0 = 1 μm) = 1.7 in the Mix 3:ice case.

thumbnail Fig. 6

Spectral indices of the SEDs presented in Fig. 4 (optically thin medium), calculated over 1 GHz frequency steps (solid lines). The circles display the spectral indices calculated using a modified blackbody in the Wien part of the SEDs. The same colour code is used as in Fig. 4.

thumbnail Fig. 7

Top: SEDs for the dust composition Mix 1:50 described in Sect. 2.1 and Table 1. The SEDs were obtained considering an optically thin medium with a column density NH = 1020 H cm−2 illuminated by the standard ISRF scaled by G0 = 10−3, 10−2, 10−1, 1, 10, 102, and 103 plotted in orange, magenta, green, red, black, cyan, and brown, respectively. Bottom: corresponding spectral indices calculated over 1 GHz frequency steps. The spectral indices for a0 = 1 cm and G0 = 10−3, 10−2, and 10−1 are overlapping.

thumbnail Fig. 8

Spectral energy distributions for the dust compositions described in Sect. 2.1and Table 1. The SEDs were obtained considering a cloud with central density nC = 105 H cm−3 (Eq. (2)), illuminated either by the standard ISRF alone (solid lines) or in addition to an internal heating source (dashed lines, see Sect. 2.3 for details). The orange lines show the case where the centroid of the log-normal size distribution is a0 = 0.1 μm, the pink lines where a0 = 1 μm, the green lines where a0 = 10 μm, the red lines where a0 = 100 μm, the black lines where a0 = 1 mm, and the cyanlines where a0 = 1 cm.

thumbnail Fig. 9

Spectral indices of the SEDs presented in Fig. 8 (optically thick medium), calculated over 1 GHz frequency steps (lines). The symbols display the spectral indices calculated using a modified blackbody in the Wien part of the SEDs. The first case, clouds illuminated by the standard ISRF alone, is displayed by solid lines and circles, and the second case, clouds illuminated by the ISRF in addition to an internal heating source, by dashed lines and squares. The same colour code is used as in Fig. 8.

4.3 Spectral indices in ALMA and VLA bands

Now that we have a clear picture of the shape of the dust SED and of its variations in optically thick regions, we present the spectral indices calculated between ALMA and VLA bands in Fig. 10 for the same dust and cloud characteristics as in the previous Sect. 4.2. The spectral indices are calculated between ALMA bands 3 and 4 (centred at 2.8 and 1.9 mm, respectively), between ALMA bands 6 and 7 (centred at 1.3 and 0.87 mm, respectively), between ALMA band 3 and VLA K-band (1.3 cm), and between ALMA band 3 and VLA X-band (2.5 cm).

The case of a cloud externally heated by the ISRF only is shown in the left plots in Fig. 10, where the spectral indices calculated between ALMA bands 6 and 7 are always lower than those calculated between bands 3 and 4. Indeed, the former bandpositions roughly match the SED peak position whereas the latter ones are in the Rayleigh domain. In both cases, the spectral indices decrease with increasing size distribution centroid in agreement with Fig. 9. For bands 6 and 7, 0 ≤ β ≤ 1 implies that a0 ≲ 10 μm and β <0 otherwise. For bands 3 and 4, β is greater than one when a0 ≲ 10 μm, 0 ≤ β ≤ 1 for 10 ≲ a0 ≲ 100 μm and negative otherwise. When we consider the spectral indices calculated between ALMA band 3 and VLA bands K and X, we see that, for grains of mixed composition and excluding those containing a-C:H, its value depends mainly on their size. An increase in β around a0 = 100 μm is found as already shown in Fig. 9 (see also Fig. 6 and associated text in Sect. 4.1). In both cases, β > 1 means that a0 < 100 μm, 0 ≤ β ≤ 1 that 100 μm < a0 < 1 mm, and β is negative when a0 > 1 mm.

The case of a cloud with an internal source is shown in the right plots in Fig. 10. As in the previous case, the spectral indices calculated between ALMA bands 6 and 7 are inferior to those calculated between ALMA bands 3 and 4, both being always higher than in the case of clouds only externally heated. For a0 = 0.1 and 1 μm, the increase is around 0.5. A decrease in the spectral index for increasing grain sizes is again found, with one notable exception: size distributions centred at a0 = 10 μm, for which the spectral indices are higher than for smaller and larger size distribution centroids. The magnitude of this increase depends on the grain material composition and is explained by the shape of the SEDs presented in Fig. 10. From these SEDs, one can see that the connection between the contributions to the total emitted flux of the cloud inner and outer layers (see Sect. 4.2) arises between 1 and 3 mm, thus matching the positions of the considered ALMA bands. The exact size distribution centroid at which this effect becomes visible depends mostly on the internal source characteristics. However, despite the simplicity of our toy model, this clearly illustrates the difficulty of assessing dust grain sizes only from the millimetre slope of the SED in YSOs. In the test case shown in Fig. 10, such high spectral index values could be mistaken for grain size distributions centred on smaller grains. For spectral indices calculated between ALMA band 3 and VLA bands K and X, the trends are similar to those of clouds without an internal source. The only difference is slightly higher values for a0 > 100 μm, with β ~ 0 ± 0.2.

Figure 11 shows the influence of the cloud central density on the spectral index measured between the same ALMA and VLA bands: six nC values are considered, nC = 105, 106, 107, 108, 109, and 1010 H cm−3. In the case of clouds only externally illuminated by the ISRF, the cloud density influence is small for all centroid values except of the smallest one, a0 = 0.1 μm. Indeed, only the smallest grains have heat capacities small enough to be sensitive to the small variations in the radiation field intensity at the cloud centre induced by the increase in nC. Moreover, when bands at longer wavelengths are considered, the variations in β decrease. In the case of clouds with internal sources, the differences are more significant in ALMA bands: from 0.12 to 0.27 for spectral indices calculated between bands 3 and 4, and from 0.06 to 0.16 between bands 6 and 7. This results from the decrease in the radiation field intensity at the cloud centres with increasing nC leading to lower grain temperatures and thus slightly shallower dust SEDs. As in the case without internal source, the variations with nC are reduced to almost nothing when one considers the spectral indices calculated between the ALMA band 3 and the VLA bands. This impliesthat in the case of observations for which no information on the density of the region is available, it is preferable to consider thelongest possible wavelengths to determine the grain size.

thumbnail Fig. 10

Spectral indices of the SEDs presented in Fig. 8 obtained considering a cloud with central density nC = 105 H cm−3 (Eq. (2)). From top to bottom: the spectral indices are computed between ALMA bands 6 and 7 (centred at 1.3 and 0.87 mm, respectively), ALMA bands 3 and 4 (centred at 2.8 and 1.9 mm, respectively), ALMA band 3 and VLA K band at 23 GHz (≡1.3 cm), and ALMA band 3 and VLA X band at 12 GHz (≡2.5 cm). Same colour code as in Fig. 3 for the grain composition. Left: cloud externally illuminated by the standard ISRF alone. Right: cloud with an additional internal heating source (see Sect. 2.3 for details).

5 Conclusion

The aim of this paper is to understand the influence of the grain composition and size distribution on the shape of their SED in dense interstellar regions such as molecular clouds, prestellar cores, YSOs, and protoplanetary discs. Three main results emerge from this study.

  • 1.

    Modifiedblackbody fits of dust SEDs in combination with the choice of a mass absorption coefficient reference value at a given FIR/submillimetre/millimetre wavelength is often used to determine dust masses and subsequently interstellar object masses. It is well known that the mass absorption coefficient strongly depends on the dust composition, leading to errors on mass estimates up to a factor of four for the compositions considered in our study. We emphasise that the choice of the mass absorption coefficient is also strongly dependent on the dust size distribution, potentially leading to errors greater than one order of magnitude.

  • 2.

    It appears almost impossible to determine the exact composition of dust grains from the FIR/submillimetre/millimetre SED shape alone. Indeed, from one composition to another, the variation in the position of the SED peak does not seem significant enough to be discriminatory unless one has very precise knowledge of the local radiation field and the geometry of the region. Moreover, it can be noted that for all grains containing carbon, regardless of its proportion to silicate, their spectral index is completely determined by this carbon (valid for both a-C and a-C:H). In addition, neither porosity nor the presence of ice mantles influence the steepness of the slope. Combining information at these wavelengths with observations at shorter wavelengths, and in particular with data showing spectral features, characteristics of the materials expected in dense ISM regions may lift the degeneracies (something that could be done with the instruments onboard the James Webb Space Telescope – JWST – and the SPace Infrared telescope for Cosmology and Astrophysics – SPICA). Even if not presented in this paper, the spectral shape of the water ice and silicate bands also depends on the dust grains exact composition (and shape).

  • 3.

    Lastly, as expected, determining the grain size from the SED (sub)millimetre to centimetre spectral index alone is not an easy task. In particular, using the “true” opacity spectral index value directly from the optical properties may not always be a good idea. The spectral index of the SED can indeed be strongly modified by the dust temperature distribution along the line of sight, in particular when observing regions with internal heating sources and/or complex geometries that were beyond the scope of this study (i.e. YSOs and protoplanetary discs). Finally, considering the longest possible wavelengths is better to ensure that observations are done in the Rayleigh domain (far from the SED peak) and in optically thin conditions.

In conclusion, the determination of grain masses, sizes, and compositions from only (sub)millimetre to centimetre continuum data is highly uncertain, a situation that can be improved with complementary observations in the mid- to far-infrared.

thumbnail Fig. 11

Spectral indices of SEDs obtained considering clouds with Mix 1:50 grains and central densities nC = 105, 106, 107, 108, 109, and 1010 H cm−3 (Eq. (2)), plotted in orange, green, red, blue, brown, and magenta, respectively. From top to bottom: the spectral indices are computed between ALMA bands 6 and 7 (centred at 1.3 and 0.87 mm, respectively), ALMA bands 3 and 4 (centred at 2.8 and 1.9 mm, respectively), ALMA band 3 and K band at 23 GHz (≡1.3 cm), and ALMA band 3 and X band at 12 GHz (≡2.5 cm). Left: cloud externally illuminated by the standard ISRF alone. Right: cloud with an additional internal heating source (see Sect. 2.3 for details).

Acknowledgements

We would like to thank A. Dutrey, A. Maury, and S. Guilloteau for fruitful discussions. We also thank M. Juvela and V. Ossenkopf for making their softwares publicly available (for radiative transfer and dust optical properties calculations, respectively). Finally, we thank our anonymous referee whose careful reading and interesting comments helped to clarify and improve the paper. I.J.-S. acknowledges the support from the STFC through an Ernest Rutherford fellowship (grant number ST/L004801), from the FEDER funding under grant ESP2017-86582-C4-1-R, and from the Spanish State Research Agency (AEI) through project number MDM-2017-0737 Unidad de Excelencia “María de Maeztu”– Centro de Astrobiología (INTA-CSIC). A.P.J., L.V., and N.Y. acknowledge the support of the Programme National PCMI of CNRS/INSU with INC/INP co-funded by CEA and CNES.

Appendix A: Mass absorption coefficient

Tables A.1A.9 show the mass absorption coefficient κ values for the dust compositions presented in Table 1 in the case of log-normal and power-law size distributions. The κ values are presented for six wavelengths: 250 and 500 μm, 0.87, 1.3, 1.9, and 2.8 mm. The trends are similar to those described in Sect. 3.2 for all wavelengths and compositions, even if it should be noticed that, in the case of a log-normal size distribution, the peak size of κ depends on the considered wavelength.

Table A.1

Mass absorption coefficient κ for grains made of a-Sil at λ = 250 and 500 μm, and 0.87, 1.3, 1.9, and 2.8 mm.

Table A.2

Same as Table A.1 but for grains made of a-C.

Table A.3

Same as Table A.1 but for grains made of a-C:H.

Table A.4

Same as Table A.1 but for grains made of Mix 1.

Table A.5

Same as Table A.1 but for grains made of Mix 2.

Table A.6

Same as Table A.1 but for grains made of Mix 1:50.

Table A.7

Same as Table A.1 but for grains made of Mix 1:ice.

Table A.8

Same as Table A.1 but for grains made of Mix 3.

Table A.9

Same as Table A.1 but for grains made of Mix 3:ice.

References

  1. Agladze, N. I., Sievers, A. J., Jones, S. A., Burlitch, J. M., & Beckwith, S. V. W. 1996, ApJ, 462, 1026 [NASA ADS] [CrossRef] [Google Scholar]
  2. Andrews, S. M., Wilner, D. J., Hughes, A. M., Qi, C., & Dullemond, C. P. 2009, ApJ, 700, 1502 [NASA ADS] [CrossRef] [Google Scholar]
  3. Arzoumanian, D., André, P., Didelon, P., et al. 2011, A&A, 529, L6 [Google Scholar]
  4. Beckwith, S. V. W., Sargent, A. I., Chini, R. S., & Guesten, R. 1990, AJ, 99, 924 [NASA ADS] [CrossRef] [Google Scholar]
  5. Birnstiel, T., Ricci, L., Trotta, F., et al. 2010, A&A, 516, L14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Birnstiel, T., Ormel, C. W., & Dullemond, C. P. 2011, A&A, 525, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Birnstiel, T., Dullemond, C. P., Zhu, Z., et al. 2018, ApJ, 869, L45 [NASA ADS] [CrossRef] [Google Scholar]
  8. Blum, J. 2018, Space Sci. Rev., 214, 52 [NASA ADS] [CrossRef] [Google Scholar]
  9. Blum, J., & Wurm, G. 2008, ARA&A, 46, 21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Bohren, C., & Huffman, D. 1983, Absorption and Scattering of Light by Small Particles (New York: John Wiley & Sons) [Google Scholar]
  11. Bohren, C. F., & Huffman, D. R. 1998, Absorption and Scattering of Light by Small Particles (New York: John Wiley & Sons), 544 [Google Scholar]
  12. Boudet, N., Mutschke, H., Nayral, C., et al. 2005, ApJ, 633, 272 [NASA ADS] [CrossRef] [Google Scholar]
  13. Busquet, G., Girart, J. M., Estalella, R., et al. 2019, A&A, 623, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Campeggio, L., Strafella, F., Maiolo, B., Elia, D., & Aiello, S. 2007, ApJ, 668, 316 [NASA ADS] [CrossRef] [Google Scholar]
  15. Chiang, H.-F., Looney, L. W., & Tobin, J. J. 2012, ApJ, 756, 168 [NASA ADS] [CrossRef] [Google Scholar]
  16. Choi, M., Kang, M., Lee, J.-E., et al. 2017, ApJS, 232, 24 [NASA ADS] [CrossRef] [Google Scholar]
  17. Compiègne, M., Verstraete, L., Jones, A., et al. 2011, A&A, 525, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Coupeaud, A., Demyk, K., Meny, C., et al. 2011, A&A, 535, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. D’Alessio, P., Calvet, N., & Hartmann, L. 2001, ApJ, 553, 321 [NASA ADS] [CrossRef] [Google Scholar]
  20. Davis, D. R., & Ryan, E. V. 1990, Icarus, 83, 156 [NASA ADS] [CrossRef] [Google Scholar]
  21. Demyk, K., Meny, C., Lu, X.-H., et al. 2017, A&A, 600, A123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Dominik, C., Paszun, D., & Borel, H. 2016, ArXiv e-prints [arXiv:1611.00167] [Google Scholar]
  23. Draine, B. T. 2006, ApJ, 636, 1114 [Google Scholar]
  24. Gerin, M., Pety, J., Commerçon, B., et al. 2017, A&A, 606, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Guilloteau, S., Dutrey, A., Piétu, V., & Boehler, Y. 2011, A&A, 529, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Gundlach, B., & Blum, J. 2015, ApJ, 798, 34 [NASA ADS] [CrossRef] [Google Scholar]
  27. Güttler, C., Krause, M., Geretshauser, R. J., Speith, R., & Blum, J. 2009, ApJ, 701, 130 [NASA ADS] [CrossRef] [Google Scholar]
  28. Güttler, C., Blum, J., Zsom, A., Ormel, C. W., & Dullemond, C. P. 2010, A&A, 513, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Isella, A., Carpenter, J. M., & Sargent, A. I. 2010, ApJ, 714, 1746 [NASA ADS] [CrossRef] [Google Scholar]
  30. Isella, A., Pérez, L. M., & Carpenter, J. M. 2012, ApJ, 747, 136 [NASA ADS] [CrossRef] [Google Scholar]
  31. Jones, A. P. 2012a, A&A, 540, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Jones, A. P. 2012b, A&A, 540, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Jones, A. P. 2012c, A&A, 542, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Jones, A. P., Fanciullo, L., Köhler, M., et al. 2013, A&A, 558, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Jones, A. P., Köhler, M., Ysard, N., et al. 2016, A&A, 588, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Jones, A. P., Köhler, M., Ysard, N., Bocchio, M., & Verstraete, L. 2017, A&A, 602, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Juvela, M. 2005, A&A, 440, 531 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Juvela, M., & Padoan, P. 2003, A&A, 397, 201 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Juvela, M., Demyk, K., Doi, Y., et al. 2015, A&A, 584, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Kataoka, A., Tanaka, H., Okuzumi, S., & Wada, K. 2013, A&A, 557, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Köhler, M., Guillet, V., & Jones, A. 2011, A&A, 528, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Köhler, M., Stepnik, B., Jones, A. P., et al. 2012, A&A, 548, A61 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Köhler, M., Jones, A., & Ysard, N. 2014, A&A, 565, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Köhler, M., Ysard, N., & Jones, A. P. 2015, A&A, 579, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Kwon, W., Looney, L. W., Mundy, L. G., & Welch, W. J. 2015, ApJ, 808, 102 [NASA ADS] [CrossRef] [Google Scholar]
  46. Liu, C., Yin, Y., Hu, F., Jin, H., & Sorensen, C. M. 2015, Aerosol Sci. Technol., 49, 928 [NASA ADS] [CrossRef] [Google Scholar]
  47. Liu, Y., Henning, T., Carrasco-González, C., et al. 2017, A&A, 607, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Lorek, S., Lacerda, P., & Blum, J. 2018, A&A, 611, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Mathis, J. S., Mezger, P. G., & Panagia, N. 1983, A&A, 128, 212 [NASA ADS] [Google Scholar]
  50. Maxwell Garnett, J. 1904, R. Soc. London Philos. Trans., Ser. A, 203 [Google Scholar]
  51. Mennella, V., Brucato, J. R., Colangeli, L., et al. 1998, ApJ, 496, 1058 [NASA ADS] [CrossRef] [Google Scholar]
  52. Meny, C., Gromov, V., Boudet, N., et al. 2007, A&A, 468, 171 [Google Scholar]
  53. Mie, G. 1908, Ann. Phys., 330, 377 [Google Scholar]
  54. Miettinen, O., Harju, J., Haikala, L. K., & Juvela, M. 2012, A&A, 538, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Min, M., Rab, C., Woitke, P., Dominik, C., & Ménard, F. 2016, A&A, 585, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Miotello, A., Testi, L., Lodato, G., et al. 2014, A&A, 567, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Mishchenko, M. I., Dlugach, J. M., & Liu, L. 2016a, J. Quant. Spectr. Rad. Transf., 178, 284 [NASA ADS] [CrossRef] [Google Scholar]
  58. Mishchenko, M. I., Dlugach, J. M., Yurkin, M. A., et al. 2016b, Phys. Rep., 632, 1 [NASA ADS] [CrossRef] [Google Scholar]
  59. Natta, A., & Testi, L. 2004, ASP Conf. Ser., 323, 279 [Google Scholar]
  60. Natta, A., Testi, L., Calvet, N., et al. 2007, Protostars and Planets V (Tucson, AZ: University of Arizona Press), 767 [Google Scholar]
  61. Ostriker, J. 1964, ApJ, 140, 1056 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  62. Pagani, L., Steinacker, J., Bacmann, A., Stutz, A., & Henning, T. 2010, Science, 329, 1622 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  63. Paszun, D., & Dominik, C. 2009, A&A, 507, 1023 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Pérez, L. M., Carpenter, J. M., Chandler, C. J., et al. 2012, ApJ, 760, L17 [NASA ADS] [CrossRef] [Google Scholar]
  65. Pérez, L. M., Chandler, C. J., Isella, A., et al. 2015, ApJ, 813, 41 [NASA ADS] [CrossRef] [Google Scholar]
  66. Planck Collaboration XXII. 2011, A&A, 536, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Pollack, J. B., Hollenbach, D., Beckwith, S., et al. 1994, ApJ, 421, 615 [NASA ADS] [CrossRef] [Google Scholar]
  68. Remy, Q., Grenier, I. A., Marshall, D. J., & Casandjian, J. M. 2017, A&A, 601, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Remy, Q., Grenier, I. A., Marshall, D. J., & Casandjian, J. M. 2018, A&A, 616, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Ribas, Á., Espaillat, C. C., Macías, E., et al. 2017, ApJ, 849, 63 [NASA ADS] [CrossRef] [Google Scholar]
  71. Schlafly, E. F., Meisner, A. M., Stutz, A. M., et al. 2016, ApJ, 821, 78 [NASA ADS] [CrossRef] [Google Scholar]
  72. Scott, A., & Duley, W. W. 1996, ApJS, 105, 401 [NASA ADS] [CrossRef] [Google Scholar]
  73. Sheehan, P. D., & Eisner, J. A. 2017, ApJ, 851, 45 [NASA ADS] [CrossRef] [Google Scholar]
  74. Steinacker, J., Andersen, M., Thi, W.-F., et al. 2015, A&A, 582, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Tafalla, M., Myers, P. C., Caselli, P., Walmsley, C. M., & Comito, C. 2002, ApJ, 569, 815 [NASA ADS] [CrossRef] [Google Scholar]
  76. Tafalla, M., Myers, P. C., Caselli, P., & Walmsley, C. M. 2004, A&A, 416, 191 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Tanaka, H., Inaba, S., & Nakazawa, K. 1996, Icarus, 123, 450 [Google Scholar]
  78. Tazaki, R., & Tanaka, H. 2018, ApJ, 860, 79 [NASA ADS] [CrossRef] [Google Scholar]
  79. Tazzari, M., Testi, L., Ercolano, B., et al. 2016, A&A, 588, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Tripathi, A., Andrews, S. M., Birnstiel, T., et al. 2018, ApJ, 861, 64 [NASA ADS] [CrossRef] [Google Scholar]
  81. Trotta, F., Testi, L., Natta, A., Isella, A., & Ricci, L. 2013, A&A, 558, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  82. Wada, K., Tanaka, H., Suyama, T., Kimura, H., & Yamamoto, T. 2008, ApJ, 677, 1296 [NASA ADS] [CrossRef] [Google Scholar]
  83. Warren, S. G. 1984, Appl. Opt., 23, 1206 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  84. Weidling, R., Güttler, C., Blum, J., & Brauer, F. 2009, ApJ, 696, 2036 [NASA ADS] [CrossRef] [Google Scholar]
  85. Whittet, D. C. B., Gerakines, P. A., Hough, J. H., & Shenoy, S. S. 2001, ApJ, 547, 872 [NASA ADS] [CrossRef] [Google Scholar]
  86. Windmark, F., Birnstiel, T., Güttler, C., et al. 2012, A&A, 540, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  87. Wu, Y., Cheng, T., Zheng, L., & Chen, H. 2016, J. Quant. Spectr. Rad. Transf., 168, 158 [NASA ADS] [CrossRef] [Google Scholar]
  88. Wuchterl, G., & Tscharnuter, W. M. 2003, A&A, 398, 1081 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  89. Young, A. K., Bate, M. R., Mowat, C. F., Hatchell, J., & Harries, T. J. 2018, MNRAS, 474, 800 [NASA ADS] [CrossRef] [Google Scholar]
  90. Ysard, N., Juvela, M., Demyk, K., et al. 2012, A&A, 542, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  91. Ysard, N., Köhler, M., Jones, A., et al. 2016, A&A, 588, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  92. Ysard, N., Jones, A. P., Demyk, K., Boutéraon, T., & Koehler, M. 2018, A&A, 617, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

2

For the full details of the model see: Jones (2012a,b,c); Jones et al. (2013, 2016); Köhler et al. (2014, 2015); Ysard et al. (2016, 2018).

4

In a realistic dust model, the gas-to-dust mass ratio should decrease in the case where gas species are accreted onto the grains, for instance in the form of water ice mantles, but this is beyond the scope of this study.

6

For more details about the comparison of exact computation methods with common approximate methods, we refer the reader to Min et al. (2016) and also to the exhaustive study of Tazaki & Tanaka (2018).

7

The optical properties of pure a-C:H grains are discussed in details in Jones et al. (2016) and Ysard et al. (2018).

All Tables

Table 1

Grain compositions, structures, and volume densities in g cm−3 as described in Sect. 2.1.

Table A.1

Mass absorption coefficient κ for grains made of a-Sil at λ = 250 and 500 μm, and 0.87, 1.3, 1.9, and 2.8 mm.

Table A.2

Same as Table A.1 but for grains made of a-C.

Table A.3

Same as Table A.1 but for grains made of a-C:H.

Table A.4

Same as Table A.1 but for grains made of Mix 1.

Table A.5

Same as Table A.1 but for grains made of Mix 2.

Table A.6

Same as Table A.1 but for grains made of Mix 1:50.

Table A.7

Same as Table A.1 but for grains made of Mix 1:ice.

Table A.8

Same as Table A.1 but for grains made of Mix 3.

Table A.9

Same as Table A.1 but for grains made of Mix 3:ice.

All Figures

thumbnail Fig. 1

Absorption (Qabs, solid lines) and scattering efficiencies (Qsca, dashed lines) for the dust compositions described in Sect. 2.1 and Table 1. Blue lines show the case of grains with a radius of 0.01 μm, orange of 0.1 μm, pink of 1 μm, green of 10 μm, red of 100 μm, black of 1 mm, cyan of 1 cm, and brown of 10 cm.

In the text
thumbnail Fig. 2

Left: mass absorption coefficients κ at 1.3 mm for log-normal size distributions centred on a0 and for the dust compositions described in Sect. 2.1 and Table 1. Black, blue, pink, green, brown, orange, red, grey, and magenta lines show the cases of a-Sil, a-C, a-C:H, Mix 1, Mix 2, Mix 1:50, Mix 1:ice, Mix 3, and Mix 3:ice grains, respectively. Mass absorption coefficients of pure a-C:H grains are always lower than 10−4 cm2 g−1 and thus not visible on the figures. The black dashed line shows κ = 1.8 cm2 g−1, a value close to those of the models commonly used in the literature with a gas-to-dust ratio of 100 (e.g. Beckwith et al. 1990; D’Alessio et al. 2001; Andrews et al. 2009; Birnstiel et al. 2018). Middle: same for power-law size distributions d n∕d aap with variable exponent p and 0.01 μm ≤ a ≤ 10 cm. Right: same for power-law size distributions with p= − 3.5, amin = 0.01 μm, and variable amax.

In the text
thumbnail Fig. 3

Equilibrium temperature as a function of grain size for the dust compositions described in Sect. 2.1 and Table 1, illuminated by the ISRF with G0 = 1. The same colour code is used as in Fig. 2.

In the text
thumbnail Fig. 4

Spectral energy distributions for the dust compositions described in Sect. 2.1 and Table 1. The SEDs were obtained considering an optically thin medium with a column density NH = 1020 H cm−2 illuminated by the standard ISRF (see Sect. 2.3 for details). The orange lines show the case where the centroid of the log-normal size distribution is a0 = 0.1 μm, the pink lines where a0 = 1 μm, the green lines where a0 = 10 μm, the red lines where a0 = 100 μm, the black lines where a0 = 1 mm, and the cyanlines where a0 = 1 cm. The thick grey curves on the Mix 1:50 plot present the case of power-law size distributions with amin = 0.01 μm and where the dotted line shows the case of p =−3.5 and amax = 1 μm, the solid lineof p =−3.5 and amax = 10 cm, and the dashed line of p = −2 and amax = 10 cm.

In the text
thumbnail Fig. 5

Wavelength of the maxima of the SEDs shown in Fig. 4 as a function of the size distribution centroid a0. The same colour code is used as in Fig. 2.

In the text
thumbnail Fig. 6

Spectral indices of the SEDs presented in Fig. 4 (optically thin medium), calculated over 1 GHz frequency steps (solid lines). The circles display the spectral indices calculated using a modified blackbody in the Wien part of the SEDs. The same colour code is used as in Fig. 4.

In the text
thumbnail Fig. 7

Top: SEDs for the dust composition Mix 1:50 described in Sect. 2.1 and Table 1. The SEDs were obtained considering an optically thin medium with a column density NH = 1020 H cm−2 illuminated by the standard ISRF scaled by G0 = 10−3, 10−2, 10−1, 1, 10, 102, and 103 plotted in orange, magenta, green, red, black, cyan, and brown, respectively. Bottom: corresponding spectral indices calculated over 1 GHz frequency steps. The spectral indices for a0 = 1 cm and G0 = 10−3, 10−2, and 10−1 are overlapping.

In the text
thumbnail Fig. 8

Spectral energy distributions for the dust compositions described in Sect. 2.1and Table 1. The SEDs were obtained considering a cloud with central density nC = 105 H cm−3 (Eq. (2)), illuminated either by the standard ISRF alone (solid lines) or in addition to an internal heating source (dashed lines, see Sect. 2.3 for details). The orange lines show the case where the centroid of the log-normal size distribution is a0 = 0.1 μm, the pink lines where a0 = 1 μm, the green lines where a0 = 10 μm, the red lines where a0 = 100 μm, the black lines where a0 = 1 mm, and the cyanlines where a0 = 1 cm.

In the text
thumbnail Fig. 9

Spectral indices of the SEDs presented in Fig. 8 (optically thick medium), calculated over 1 GHz frequency steps (lines). The symbols display the spectral indices calculated using a modified blackbody in the Wien part of the SEDs. The first case, clouds illuminated by the standard ISRF alone, is displayed by solid lines and circles, and the second case, clouds illuminated by the ISRF in addition to an internal heating source, by dashed lines and squares. The same colour code is used as in Fig. 8.

In the text
thumbnail Fig. 10

Spectral indices of the SEDs presented in Fig. 8 obtained considering a cloud with central density nC = 105 H cm−3 (Eq. (2)). From top to bottom: the spectral indices are computed between ALMA bands 6 and 7 (centred at 1.3 and 0.87 mm, respectively), ALMA bands 3 and 4 (centred at 2.8 and 1.9 mm, respectively), ALMA band 3 and VLA K band at 23 GHz (≡1.3 cm), and ALMA band 3 and VLA X band at 12 GHz (≡2.5 cm). Same colour code as in Fig. 3 for the grain composition. Left: cloud externally illuminated by the standard ISRF alone. Right: cloud with an additional internal heating source (see Sect. 2.3 for details).

In the text
thumbnail Fig. 11

Spectral indices of SEDs obtained considering clouds with Mix 1:50 grains and central densities nC = 105, 106, 107, 108, 109, and 1010 H cm−3 (Eq. (2)), plotted in orange, green, red, blue, brown, and magenta, respectively. From top to bottom: the spectral indices are computed between ALMA bands 6 and 7 (centred at 1.3 and 0.87 mm, respectively), ALMA bands 3 and 4 (centred at 2.8 and 1.9 mm, respectively), ALMA band 3 and K band at 23 GHz (≡1.3 cm), and ALMA band 3 and X band at 12 GHz (≡2.5 cm). Left: cloud externally illuminated by the standard ISRF alone. Right: cloud with an additional internal heating source (see Sect. 2.3 for details).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.