Free Access
Issue
A&A
Volume 641, September 2020
Article Number A75
Number of page(s) 13
Section Astrophysical processes
DOI https://doi.org/10.1051/0004-6361/202037858
Published online 10 September 2020

© ESO 2020

1. Introduction

The excess emission commonly observed in spectral energy distributions (SEDs) reflects the properties and the evolution of the emitting circumstellar dust. Material from a range of grain sizes, radial distances, grain compositions, etc. contributes to the observed flux density. In turn, the distributions across these grain properties result from the past and ongoing processes that create and remove the dust.

The shape of an observed disc-integrated SED is a convolution of several dust properties and processes. A prominent example for the degeneracies that occur is the small grains at greater distances that can have the same temperatures as bigger grains closer to the star, where big and small are relative to the stellar spectrum. In a similar way, the temperature distribution and the resulting SED can be broadened by both a range of grain sizes and a range of distances. A narrow ring containing small and big grains can mimic a wide ring with a narrow size distribution.

The analysis of long-wavelength emission suffers less from a temperature degeneracy. Grains are inefficient emitters in the Rayleigh regime, that is, at wavelengths much longer than their own scale. At long wavelengths, the contribution from smaller grains to the SED is small. With these warmer grains out of the equation, the range of temperatures involved is reduced; temperatures and grain size become less entangled. This allowed, for example, Draine (2006) to derive an analytic approximation that relates the spectral index in the mm regime, αmm, to the power-law index of the underlying grain size distribution, p. The linear relation that results is:

p = α mm α Pl β s + 3 , $$ \begin{aligned} p = \frac{\alpha _{\rm mm} - \alpha _{\rm Pl}}{\beta _{\rm s}} + 3, \end{aligned} $$(1)

where αPl is the slope of the Planck function at the appropriate temperature and βp is the power-law index of the wavelength-dependent opacity.

Knowledge of the size distribution index provides deep insights into the physical environment of the dust. In protoplanetary discs, the processes of agglomeration and compaction result in net grain growth and shape the size distribution (e. g. Dominik et al. 2007). In debris discs, disruption and erosion steadily transfer mass from bigger to smaller objects in a collisional cascade (e. g. Krivov 2010; Matthews et al. 2014). On the one hand, the resulting dust size distribution is steepened by the material strength decreasing with grain size (O’Brien & Greenberg 2003; Wyatt et al. 2011) and by successive collisional damping as material is moving down the collisional cascade towards smaller sizes (Pan & Schlichting 2012). On the other hand, an interplay of radiation pressure and dynamical excitation (Thébault & Wu 2008) and transport processes such as Poynting–Robertson or stellar wind drag (Reidemeister et al. 2011; Wyatt et al. 2011) flatten the size distribution. The canonical value of p = 3.5 is reached for constant impact velocities and material strength (Dohnanyi 1969). Plausible ranges for average size distribution slopes roughly cover 3 ≲ p ≲ 4.

The small-size cutoff due to radiation pressure blowout induces a wavy deviation where slopes can exceed this range locally (Campo Bagatin et al. 1994). That waviness partially transcends to the observed SEDs (Thébault & Augereau 2007). However, Eq. (1) is derived under certain assumptions, with alternative approximations available for other assumptions. The aim of this work is, therefore, to delimit the parameter ranges in which these approximations apply and when it is best to use a full numerical modelling and fitting of thermal emission, which can be considered the benchmark.

Section 2 provides a review of the basic equations that govern all thermal emission models. Possible approximations are discussed, including those that lead to Eq. (1). Section 3 contains a comparison of numerical and analytic results, revealing and explaining discrepancies that significantly affect the desired direct link between observables and dust physics. The immediate impact on the analysis of literature data is discussed in Section 3.7. Sections 4 and 5 discuss the implications more broadly and summarize the findings.

2. Analytic approximations

An optically thin circumstellar dust ring at a distance d from the observer produces a total thermal flux density,

F ν ( λ ) = s min s max π B ν [ T ( s ) , λ ] Q abs ( T ( s ) , s , λ ) s 2 d 2 N ( s ) d s , $$ \begin{aligned} F_\nu (\lambda ) = \int\limits_{s_{\rm min}}^{s_{\rm max}} \pi B_\nu [T(s), \lambda ] \, Q_{\rm abs}(T(s), s, \lambda ) \, \frac{s^2}{d^2} \, N(s) \, \mathrm{d} s, \end{aligned} $$(2)

with Bν the Planck function, T(s) the grain temperature, λ the wavelength, Qabs the absorption and emission efficiency, N(s) the differential size distribution, and smin and smax the minimum and maximum grain radii. In general, the absorption efficiency depends on wavelength, material, and temperature as well as grain size and shape. In thermal equilibrium, grain temperature T and distance r from the star are related through (Krivov et al. 2008):

r = R 2 C abs ( T , s , λ ) F λ , d λ C abs ( T , s , λ ) B λ ( T , λ ) d λ , $$ \begin{aligned} r = \frac{R_*}{2} \sqrt{\frac{\int C_{\rm abs}(T, s, \lambda ) F_{\lambda ,*} \mathrm{d}\lambda }{\int C_{\rm abs}(T, s, \lambda ) B_\lambda (T, \lambda ) \mathrm{d}\lambda }}, \end{aligned} $$(3)

where Cabs ≡ πs2Qabs is the effective absorption cross-section per grain. Whenever the temperature can be considered constant over the relevant range of grain sizes, Eq. (2) can be simplified to

F ν ( λ ) B ν ( T , λ ) d 2 s min s max C abs ( T , s , λ ) N ( s ) d s . $$ \begin{aligned} F_\nu (\lambda ) \approx \frac{B_\nu (T, \lambda )}{d^2} \, \int\limits_{s_{\rm min}}^{s_{\rm max}} C_{\rm abs}(T, s, \lambda ) \, N(s) \, \mathrm{d} s. \end{aligned} $$(4)

Grain temperature is discussed in greater detail in Sect. 3.5.

Defining further a total cross-section,

C abs , total ( λ ) C abs ( s , λ ) N ( s ) d s $$ \begin{aligned} C_{\rm abs,total}(\lambda ) \equiv \int C_{\rm abs}(s, \lambda ) N(s) \mathrm{d} s \end{aligned} $$(5)

and disc mass

M d = M ( s ) N ( s ) d s , $$ \begin{aligned} M_{\rm d} = \int M(s) N(s) \mathrm{d} s, \end{aligned} $$(6)

where M(s) is the mass of a single grain with radius s, the wavelength-dependent total opacity is given by

κ ( λ ) C abs , total ( λ ) M d · $$ \begin{aligned} \kappa (\lambda ) \equiv \frac{C_{\rm abs,total}(\lambda )}{M_{\rm d}}\cdot \end{aligned} $$(7)

After inserting Eqs. (5)–(7) into (4), the resulting flux density is

F ν ( λ ) B ν ( T , λ ) d 2 κ ( λ ) M d . $$ \begin{aligned} F_\nu (\lambda ) \approx \frac{B_\nu (T, \lambda )}{d^2} \, \kappa (\lambda ) M_{\rm d}. \end{aligned} $$(8)

The measured flux density can thus be decomposed into the underlying blackbody and the opacity. If κ follows a power law,

κ λ β , $$ \begin{aligned} \kappa \propto \lambda ^{-\beta }, \end{aligned} $$(9)

the opacity index β is related to the spectral index via

log F ν B ν = β log λ + const . , $$ \begin{aligned} \log \frac{F_\nu }{B_\nu } = -\beta \log \lambda + \mathrm{const.}, \end{aligned} $$(10)

which is equivalent to

log F ν , 2 / F ν , 1 log λ 2 / λ 1 = α mm log B ν , 2 / B ν , 1 log λ 2 / λ 1 = α Pl = β $$ \begin{aligned} \underbrace{\frac{\log F_{\nu ,2}/F_{\nu ,1}}{\log \lambda _2/\lambda _1}}_{= -\alpha _{\rm mm}} - \underbrace{\frac{\log B_{\nu ,2}/B_{\nu ,1}}{\log \lambda _2/\lambda _1}}_{= -\alpha _{\rm Pl}} = -\beta \end{aligned} $$(11)

for any given pair of wavelengths λ1 and λ2. Assuming that the temperature is known and constant over the range of relevant grain sizes, the blackbody component can thus be removed from observed fluxes. The spectral index of the opacity, and hence the size distribution, can be derived (e.g. Ricci et al. 2012, 2015; MacGregor et al. 2016; Marshall et al. 2017).

While the details of Cabs(s, λ) depend on grain morphology, structure, and material, two asymptotic cases and an intermediate case can be identified for a given wavelength:

  • (s) grains that are geometrically small,

  • (t) grains that are large but optically thin and

  • (o) grains that are optically thick and geometrically large,

corresponding to cases 1, 2, 3 in Kataoka et al. (2014), respectively. Grains in categories (s) and (t) absorb with all their transparent volume, thus having Cabs  ∝  s3 and Qabs  ∝  s. The opaque grains in category (o) absorb all radiation that is not reflected off of their surface. They absorb (and emit) only in their surface layers, corresponding to Cabs  ∝  s2 and Qabs  ∼  1. Bohren & Huffman (1983) and Kataoka et al. (2014) give the following explicit absorption efficiencies for the three regimes:

Q abs , s 4 x Im ( m 2 1 m 2 + 2 ) = 24 n k x ( n 2 k 2 + 2 ) 2 + ( 2 n k ) 2 s , $$ \begin{aligned} Q_{\rm abs,s}&\approx 4x \mathrm{Im} \left(\frac{m^2 - 1}{m^2 + 2}\right) = \frac{24 nkx}{(n^2 - k^2 + 2)^2 + (2nk)^2} \propto s, \end{aligned} $$(12)

Q abs , t 8 k x 3 n ( n 3 ( n 2 1 ) 3 / 2 ) s , $$ \begin{aligned} Q_{\rm abs,t}&\approx \frac{8kx}{3n} \left(n^3 - (n^2 - 1)^{3/2}\right) \propto s, \end{aligned} $$(13)

Q abs , o 1 , $$ \begin{aligned} Q_{\rm abs,o}&\sim 1, \end{aligned} $$(14)

where m ≡ n + ik is the complex refractive index and x ≡ 2πs/λ the size parameter. The transition from (s) to (t) occurs around grain sizes where 2ns ≈ λ, corresponding to a characteristic size,

s c = λ 2 n · $$ \begin{aligned} s_{\rm c}^\prime = \frac{\lambda }{2n}\cdot \end{aligned} $$(15)

Regime (o) is reached when the grains become opaque, that is, around 8kx ∼ 3n, corresponding to

s c = max ( 3 n λ 16 π k , s c ) . $$ \begin{aligned} s_{\rm c} = \mathrm{max}\left(\frac{3n\lambda }{16\pi k}, s_{\rm c}^\prime \right). \end{aligned} $$(16)

Two extremes are plausible:

  • (a)

    a wide transparent regime (t), that is, s c s c $ s_{\rm c} \gg s_{\rm c}^{\prime} $

  • (b)

    a vanishing transparent regime (t) if Q abs ( s c ) $ Q_{\rm abs}(s_{\rm c}^{\prime}) $ is already close to unity because of a high A, in which case s c s c $ s_{\rm c} \approx s_{\rm c}^{\prime} $

Figure 1 illustrates cases (a) and (b) for Qabs as a function of grain size. Grains becoming opaque while still being small would imply Qabs, s ≳ 1 and hence 2n/λ ≲ 4A/3 or n ≲ 8πk/3, which is only met for very strongly absorbing materials.

thumbnail Fig. 1.

Schematic comparison of the two different cases (a) and (b) for the transition from small grains that absorb by volume to big, opaque grains that absorb by cross-section. Case (a) has an intermediate size interval where the grains are already big compared to the wavelength but still transparent. In case (b) the transition from small to opaque is direct. The values assumed in Eqs. (17), (18), and (28) are βs = 1.6, A0 = 0.63 mm−1, λ0 = 1 mm, n = 3, y = n2.

To simplify the expressions for Qabs, s and Qabs, t, the dependence on wavelength can be subsumed in the volume absorption coefficient A:

Q abs , s = 4 3 s A , $$ \begin{aligned} Q_{\rm abs,s}&= \frac{4}{3}sA,\end{aligned} $$(17)

Q abs , t = 4 3 s A y , $$ \begin{aligned} Q_{\rm abs,t}&= \frac{4}{3}sA { y}, \end{aligned} $$(18)

where y ≡ Qabs, t/Qabs, s. A commonly used alternative symbol for A is αλ. The characteristic size for the transition from (t) to (o) then becomes

s c = max ( 3 Q abs , o 4 A y , s c ) . $$ \begin{aligned} s_{\rm c} = \mathrm{max}\left(\frac{3Q_{\rm abs,o}}{4A { y}}, s_{\rm c}^\prime \right). \end{aligned} $$(19)

As y or A increase, sc approaches s c $ s_{\rm c}^{\prime} $ and the transparent regime vanishes.

For low absorption, k ≪ n, the ratio between the two regimes is given by

y ( n 2 + 2 ) 2 9 n 2 ( n 3 ( n 2 1 ) 3 / 2 ) , $$ \begin{aligned} { y} \approx \frac{(n^2 + 2)^2}{9n^2} \left(n^3 - (n^2 - 1)^{3/2}\right), \end{aligned} $$(20)

which is then independent from k. Over the range 1 <  n <  4, y can be very roughly approximated with 0.75n2. For materials with low refractive index n ≈ 1, the difference between (s) and (t) vanishes and A becomes identical with the linear attenuation coefficient, which is related to the imaginary part of the refractive index (Bohren & Huffman 1983):

A 4 π k λ , $$ \begin{aligned} A \approx \frac{4\pi k}{\lambda }, \end{aligned} $$(21)

which is independent from n.

Assuming a power-law for the grain size distribution,

N ( s ) = N 0 ( s s 0 ) p , $$ \begin{aligned} N(s) = N_0 \left(\frac{s}{s_0}\right)^{-p}, \end{aligned} $$(22)

and inserting Eqs. (17), (18) and (14) into (5) and (6), the dust mass and total absorption cross-section are

M d = 4 π ρ N 0 s 0 4 3 ( 4 p ) [ ( s max s 0 ) 4 p ( s min s 0 ) 4 p ] $$ \begin{aligned} M_{\rm d} = \frac{4\pi \rho N_0 s_0^4}{3(4 - p)} \left[\left(\frac{s_{\rm max}}{s_0}\right)^{4 - p} - \left(\frac{s_{\rm min}}{s_0}\right)^{4 - p}\right] \end{aligned} $$(23)

and

C abs , total ( λ ) = 4 π A N 0 s 0 4 3 ( 4 p ) [ ( s c s 0 ) 4 p ( s min s 0 ) 4 p + y ( s c s 0 ) 4 p y ( s c s 0 ) 4 p ] + π N 0 s 0 3 ( 3 p ) Q abs , o [ ( s max s 0 ) 3 p ( s c s 0 ) 3 p ] , $$ \begin{aligned} C_{\rm abs,total}(\lambda )&= \frac{4\pi A N_0 s_0^4}{3(4 - p)} \left[\left(\frac{s_{\rm c}^\prime}{s_0}\right)^{4 - p} - \left(\frac{s_{\rm min}}{s_0}\right)^{4 - p}\right.\nonumber \\&\quad + \left.{ y}\left(\frac{s_{\rm c}}{s_0}\right)^{4 - p} - { y}\left(\frac{s_{\rm c}^\prime}{s_0}\right)^{4 - p}\right]\nonumber \\&\quad + \frac{\pi N_0 s_0^3}{(3 - p)}Q_{\rm abs,o} \left[\left(\frac{s_{\rm max}}{s_0}\right)^{3 - p} - \left(\frac{s_{\rm c}}{s_0}\right)^{3 - p}\right], \end{aligned} $$(24)

for s min < s c < s c < s max $ s_{\rm min} < s_{\rm c}^{\prime} < s_{\rm c} < s_{\rm max} $. The opacity is then

κ ( λ ) = A ρ { ( 1 y ) ( s c s c ) 4 p ( s min s c ) 4 p + y 4 p p 3 3 Q abs , o 4 A s c [ ( s max s c ) 3 p 1 ] } × [ ( s max s c ) 4 p ( s min s c ) 4 p ] 1 , $$ \begin{aligned} \kappa (\lambda )&= \frac{A}{\rho } \left\{ (1 - { y})\left(\frac{s_{\rm c}^\prime }{s_{\rm c}}\right)^{4 - p} - \left(\frac{s_{\rm min}}{s_{\rm c}}\right)^{4 - p} + { y}\right.\nonumber \\&\quad - \left.\frac{4 - p}{p - 3}\frac{3Q_{\rm abs,o}}{4A s_{\rm c}}\left[\left(\frac{s_{\rm max}}{s_{\rm c}}\right)^{3 - p} - 1\right] \right\} \nonumber \\&\quad \times \left[ \left(\frac{s_{\rm max}}{s_{\rm c}}\right)^{4 - p} - \left(\frac{s_{\rm min}}{s_{\rm c}}\right)^{4 - p} \right]^{-1}, \end{aligned} $$(25)

where s0  =  sc is assumed. For a top-heavy size distribution, where big grains dominate both total mass and total cross-section,

κ ( λ ) p 3 3 Q abs , o 4 ρ s max 4 p 3 p , $$ \begin{aligned} \kappa (\lambda ) \mathop {\longrightarrow }\limits ^{p \ll 3} \frac{3Q_{\rm abs,o}}{4\rho s_{\rm max}} \frac{4 - p}{3 - p}, \end{aligned} $$(26)

which inherits only the dependence of Qabs, o on wavelength, roughly leading to β → 0 for p ≪ 3. For a bottom-heavy size distribution, where small grains dominate both total mass and total cross-section,

κ ( λ ) p 4 A ρ , $$ \begin{aligned} \kappa (\lambda ) \mathop {\longrightarrow }\limits ^{p \gg 4} \frac{A}{\rho }, \end{aligned} $$(27)

which is independent of p. A power-law approximation to the volume absorption coefficient,

A = A 0 ( λ λ 0 ) β s , $$ \begin{aligned} A = A_0 \left(\frac{\lambda }{\lambda _0}\right)^{-\beta _{\rm s}}, \end{aligned} $$(28)

then leads to β → βs for p ≫ 4.

In the typical intermediate regime, 3 <  p <  4, where small grains dominate total cross-section and big grains the total mass, κ can be simplified to

κ ( λ ) A ρ { y + ( 1 y ) ( s c s c ) 4 p + 4 p p 3 3 Q abs , o 4 A s c } ( s max s c ) p 4 . $$ \begin{aligned} \kappa (\lambda ) \approx \frac{A}{\rho } \left\{ { y} + (1 - { y})\left(\frac{s_{\rm c}^\prime }{s_{\rm c}}\right)^{4 - p} + \frac{4 - p}{p - 3}\frac{3Q_{\rm abs,o}}{4A s_{\rm c}} \right\} \left(\frac{s_{\rm max}}{s_{\rm c}}\right)^{p - 4}. \end{aligned} $$(29)

With Eq. (19), the two cases (a) and (b) imply the opacities

κ a ( λ ) yA ρ { 1 + 4 p p 3 } ( s max s c ) p 4 Q abs , o × s c 3 p Q abs , o × λ β s ( p 3 ) $$ \begin{aligned} \kappa _{\rm a}(\lambda )&\approx \frac{{ y}A}{\rho } \left\{ 1 + \frac{4 - p}{p - 3} \right\} \left(\frac{s_{\rm max}}{s_{\rm c}}\right)^{p - 4}\nonumber \\&\propto Q_{\rm abs,o} \times s_{\rm c}^{3 - p}\, \propto \,Q_{\rm abs,o} \times \lambda ^{-\beta _{\rm s}(p - 3)} \end{aligned} $$(30)

for a wide transparent regime and

κ b ( λ ) A ρ { 1 + 4 p p 3 3 Q abs , o 4 A s c } ( s max s c ) p 4 Q abs , o × ( s c ) 3 p Q abs , o × λ 3 p $$ \begin{aligned} \kappa _{\rm b}(\lambda )&\approx \frac{A}{\rho } \left\{ 1 + \frac{4 - p}{p - 3}\frac{3Q_{\rm abs,o}}{4A s_{\rm c}^\prime } \right\} \left(\frac{s_{\rm max}}{s_{\rm c}^\prime }\right)^{p - 4}\nonumber \\&\propto Q_{\rm abs,o} \times (s_{\rm c}^\prime )^{3 - p}\,\propto \,Q_{\rm abs,o} \times \lambda ^{3 - p} \end{aligned} $$(31)

for a direct transition from small to opaque, respectively. For y ≈ 1 and Qabs, o ≈ const, Eq. (30) is equivalent to Eq. (8) in Draine (2006).

The combination of Eqs. (30) and (11) results in Eq. (1), a relation between observed spectral index αmm and size distribution index p that depends on the material through βs. For case (b), the combination of Eqs. (31) and (11) results in

p = 3 + α mm α Pl , $$ \begin{aligned} p = 3 + \alpha _{\rm mm} - \alpha _{\rm Pl}, \end{aligned} $$(32)

which depends only weakly on the material through αPl and the grain temperature.

Table 1 lists the symbols used throughout this manuscript for quick look-up.

Table 1.

Alphabetic list of recurring symbols.

3. Comparison with numerical results

In this section, the assumptions made for the basic analytic approximations are tested against numerical models for different materials. Absorption and emission are calculated assuming homogeneous, compact spheres. For grains where 2πsn <  104λ, Mie theory as implemented by Wolf & Voshchinnikov (2004) is used. For larger grains the geometrical optics algorithm by Bohren & Huffman (1983) is used.

3.1. Materials and range of parameters

A set of five different materials is assumed: compact astronomical silicate (Draine 2006) with assumed bulk density ρ = 3.3 g cm−3, porous astronomical silicate (mixing rule: Bruggeman 1935) with a filling factor f = 10% and ρ = 0.33 g cm−3, amorphous carbon (Zubko et al. 1996) with ρ = 2 g cm−3, pyroxene (H. Mutschke and P. Mohr, private comm., manuscript in preparation) with ρ = 3.3 g cm−3, and crystalline water ice (Reinert et al. 2015) with ρ = 1.0 g cm−3. Only the data sets for the latter two materials include a dependence on grain temperature. Figure 2 shows the real and imaginary parts of the refractive indexes of these materials. All materials exhibit more or less pronounced absorption features for wavelengths λ ≲ 100 μm and smoother behaviours for λ ≳ 100 μm. Amorphous carbon is exceptional in that absorption is high across all considered wavelengths, without a significant decline in the (sub)mm range.

thumbnail Fig. 2.

Refractive indexes of the five materials discussed: (top) real and (middle, bottom) imaginary parts. The legends are ordered vertically according to the order at long wavelengths. The vertical dashed lines delimit the wavelength range over which the labelled slopes are fitted.

All numerical calculations were done for a grid of parameters. Size distribution indexes are varied from 2.6 to 4.2 in steps of 0.1. Maximum grain radii are smax = 1 km for carbon and the silicates and smax = 100 km for water ice. Minimum grain radii are determined from the radiation pressure blowout limits for main-sequence stars of six different luminosities: 0.0625, 0.25, 1, 4, 16, along with 64 solar luminosities. In between smax and smin, the grid of grain radii is spaced logarithmically with factors of 1.251/3 between neighbouring points. Distances between dust and star are set to 1.5625, 3.125, 6.25, 12.5, 25, 50, 100, and 200 au.

Two computations were performed for each of the 4080 combinations of size distribution index, luminosity, distance, and material: one with absorption and emission directly defined by the optical properties of the materials and another one were actual materials are only used for determining grain temperatures (for s >  1 mm), but emission was then calculated for blackbodies at these temperatures. The latter computations were used to derive αPl as required in Eqs. (1) and (32).

3.2. Absorption by volume and by geometrical cross-section

The resulting absorption efficiencies are plotted in Fig. 3 for a range of wavelengths and grain sizes, illustrating the occurrences of the regimes (s), (t), and (o) as defined in Sect. 2. As expected, the opaque regime (o) with Qabs ∼ 1 is found where k is high and wavelengths λ are short compared to grain radii s. The transparent regime (t) is found for s ≳ λ/(2n) but low absorption. As soon as s ≲ λ/(2n), the grains can be considered small, corresponding to regime (s). The transition radius and the ratio y = Qabs, t/Qabs, s depend, notably, on n.

thumbnail Fig. 3.

Absorption efficiencies for four materials and a range of wavelengths and grain sizes. Wavelength are separated by factors of 10, the longest wavelengths are labelled. Solid lines indicate Mie results, dashed lines the small-grain asymptote given in Eq. (12). The upper boundaries of the light-shaded regions are defined by Eq. (13), the left boundaries by Eq. (15). Black dotted lines follow Eq. (13) multiplied by an additional factor n to roughly trace the actual smoothed Qabs. Values are given for stellar luminosity L = LSun and a distance r = 100 au.

Figure 3 shows that, for silicates and carbon, the transparent regime is missing for λ ≲ 1 mm. Absorption transitions directly from small to opaque regime, corresponding to case (b) and Eq. (32) from Sect. 2. For λ ≳ 1 mm, silicates and carbon show a transparent regime. At these longer wavelengths, a transition to case (a) and Eq. (1) can be expected.

The assumption of homogeneous, perfectly spherical grains leads to artefacts in the Qabs curves. In the transparent regime, the numerical values for Qabs exceed the values predicted by Eq. (13) significantly, showing strong oscillations. This excess amounts on average to an additional factor ≈n. This is due to the monochromatic calculations and the idealized shape resulting in internal resonances, which increase absorption and emission. For more irregular grains and integration over a broader range of wavelengths, the curves are expected to become smoother. In the opaque regime, Qabs of amorphous carbon grains converges to ≈0.4, that is, a value significantly below unity. This is caused by the assumed even surfaces and the resulting Fresnel reflection on the surface of this high-n material. Surface structure on sub-wavelength scales could reduce reflection and increase absorption.

3.3. Opacities

The total mass absorption coefficient or opacity, κ, is an intermediary between absorption efficiencies and SEDs. Figure 4 shows results for four different materials as functions of size distribution indexes and wavelengths. These curves partly inherit the absorption features seen for k, but they are smoother because they tend to be dominated by grains that are rather opaque. The opacities of astronomical silicate and amorphous carbon follow smooth power-laws already for λ ≳ 100 μm. The data sets for pyroxene and water ice have measured changes in slope for λ ≳ 100 μm, which lead to a corresponding significant departure of κ from a simple power law.

thumbnail Fig. 4.

Opacities for four different materials for a range of size distribution slopes and wavelengths. Values are given for stellar luminosity L = L, a distance r = 100 au, and a total disc mass Md = 0.1MEarth. The maximum object radius is set to smax = 100 km for water ice and 1 km for carbon and the silicates. Mie results are plotted with solid lines for five size distribution indexes from (top) p = 4.2 to (bottom) p = 2.6. Sloped dashed lines trace power-law fits to wavelengths from λ1 = 1 mm to the shorter of λ2 = 7 mm and the maximum wavelength in the tabulated optical data. Vertical dashed lines indicate λ1 and λ2.

All four materials show near-constant opacities for shallow size distributions around p = 3 because in this case, big and opaque grains near to smax not only dominate the total mass but also contribute significantly to the total cross-section. The disc then has the same effective cross-section at all wavelength. The remaining slope for carbon at p = 2.6 is caused by high k, and n increasing with λ. As a result, reflectance increases and Qabs, o decreases. For the silicates and water ice, k is lower and n constant.

Towards size distribution slopes p ≥ 4, smaller grains become more and more important. Given the same total mass Md, the total cross-section Cabs, total increases, leading to higher values for κ overall. However, the small grains’ effective cross-section more strongly depends on wavelength, leading to steeper opacity slopes for steeper size distributions, as expected from both Eqs. (30) and (31). Around wavelengths of a few tens of microns, k is comparably high for all materials, and so is κ.

The wavelengths at which the κ curves for different size distribution slopes converge are the wavelengths at which the biggest assumed grains become transparent, that is, the wavelengths at which sc(λ) = smax. At yet longer wavelengths, κ would become insensitive to p. For carbon and the silicates, the extrapolated curves converge at wavelengths of several metres or even kilometres. For water ice at ∼10 K they converge already around λ ∼ 1 dm because of the very low absorption, despite the increased smax for this material.

Figure 5 compares the β(p) dependencies expected from the Mie calculations with the analytic approximations in Eqs. (30) and (31). Neither of the two approximations can reproduce the numerical results over the full range of size distribution slopes for all materials. While the solid curve for amorphous carbon approaches the grey line that follows Eq. (31), the results for pyroxene, water ice, and astronomical silicate are closer to the dashed lines that represent Eq. (30). The slopes of the regression lines for the two silicates are only halfway between that of Eqs. (30) and (31), meaning that the sensitivity of β to differences in p is also halfway between cases (a) and (b). These findings for the four materials are consistent with their refractive indexes in the mm wavelength regime. The high n and k of carbon puts it in category (b), the low absorption of pure water ice leads to case (a), with the two silicates remaining in between.

thumbnail Fig. 5.

Relation between the power-law indices of the grain size distributions and the total opacities between λ1 = 1 mm and λ2 = 7 mm for four materials at r = 100 au around a star with L = L. Solid coloured lines show numerical results. Dotted lines of the same colours show linear fits for 3.2 ≤ p ≤ 3.8, with fit coefficients indicated in the figure keys. Dashed lines show the analytic approximation according to Eq. (30). The solid grey lines trace the approximate Eq. (31).

Taking Astronomical Silicate as an example, Fig. 6 shows how the results shift more towards case (b) for shorter, submm wavelengths. In this wavelength range absorption is higher for all materials, resulting in a narrower transparent regime (t).

thumbnail Fig. 6.

Same as the top-right panel in Fig. 5 (astronomical silicate) but for two combinations of shorter wavelengths, as indicated.

Figure 7 shows that the general shape of the β(p) relation is well-reproduced by Eq. (25). The slopes of the near-linear sections around p = 3.5 are determined by the refractive index, with shallower slopes for higher n. Around p ≈ 3 and p ≈ 4 the curves turn towards the asymptotic Eqs. (26) and (27), respectively. For materials with high n, β can exceed βs for p >  4, while still converging to βs for p → ∞. Figure 8 illustrates further dependencies of Eq. (25) on the material properties and grain size ranges. A lower absorption efficiency Qabs,o, caused by a more reflective surface due to a greater n, reduces the step factor y, narrows the transparent regime and brings the results closer to Eq. (32). A wavelength dependence Qabs,o(λ)∝λξ ≠ const in the mm wavelength range, caused by n(λ)≠const, results in a non-zero β → ξ for p <  3. This effect is also seen in the numerical results for amorphous carbon in the top-left panel of Fig. 5. If the maximum grain radius is comparable to or less than the characteristic grain size sc, the opaque regime vanishes and β(p) becomes shallower. If s max s c $ s_{\rm max} \leq s_{\rm c}^{\prime} $, only the small-grain regime remains and β = βs = const. This is relevant for observed young protoplanetary discs where grain growth may not have reached smax >  sc yet (e.g. Krügel & Siebenmorgen 1994; Natta & Testi 2004; Draine 2006). It is also relevant for modelled discs, where systematic errors are introduced if the assumed smax is too small.

thumbnail Fig. 7.

Slope of the opacity as derived from Eq. (25) for the following parameters: λ = 3 mm, βs = 1.6, A0 = 2/cm, λ0 = 1 mm, ρ = 3.5 g cm−3, smin = 0.5 μm, smax = 100 m, Qabs, o = 1, y = 0.75n2, with (black) n = 8.5, (red) n = 3.4, and (blue) n = 1.14. In the range 3 <  p <  4, this order of colours reflects the vertical order of the corresponding lines, bottom to top. The dashed black and solid grey lines follow Eqs. (30) and (31), respectively.

thumbnail Fig. 8.

Same as Fig. 7, but with (top left) Qabs, o, (top right) smax, (middle) smin, and (bottom) A0 varied. Left panels: results for high n, right panels: for low n. The legends are ordered vertically according to the order at p ≈ 4.

The greater the minimum grain radius smin, the fewer small grains are present and the steeper is the size distribution required for the small grains to dominate. With increasing smin, β(p) flattens (and p(b) steepens) around p ≈ 4 because the small-grain regime βs is approached only for higher p. This dependence of β(p) on smin translates to a dependence on stellar luminosity if the latter is high enough to set a lower size cutoff via radiation pressure blowout of the smallest grains. For materials with high n the effect of different smin is smaller, as long as s min < s c $ s_{\rm min} < s_{\rm c}^{\prime} $. Because of the jump by a great factor y from Qabs, s to Qabs, t, the small grains play a minor role and the exact location of their lower size cut-off is not as important. Only as smin approaches or exceeds s c $ s_{\rm c}^{\prime} $ or sc does the β relation react by converging quickly towards a constant β ≈ 0. This is the case for carbon grains around high-luminosity stars, where radiation pressure blowout can eliminate small and transparent grains completely. For materials with low n, small grains and their lower size cutoff are more important. The lower smin, the more abundant are the small grains and the quicker is the convergence β(p)→βs around p ≈ 4.

A lower volume absorption coefficient, as represented by A0 in the bottom row of Fig. 8, brings β(p) closer to Eq. (1) in the range 3 <  p <  4 for both high-n and low-n materials. For materials with high n, this is due to a lower A causing a greater characteristic grain radius, sc, and a wider transparent regime. This effect does not depend on smin or smax. For materials with low n, where the transparent regime is less important, a lower A brings sc closer to smax and narrows the opaque regime. This effect for low n is similar to that shown in the top-right panel of Fig. 8. It disappears for smax ≫ sc.

3.4. Spectral index

For a direct link between observed spectral index and size distribution index, SEDs are computed for the range of parameters. As long as temperatures do not vary strongly across the relevant grain sizes, SEDs and their slopes simply reflect the opacity index and the slope of the blackbody function at that temperature. In Fig. 9 the analytic predictions for size distribution index p as a function of the observable spectral index in the mm regime, αmm, are compared to the numerical results. As expected from the close relation between β and αmm, the findings for the p(αmm) relation reflect those for β(p) and the curves shown in Fig. 5. The full numerical solution typically falls in between cases (a) and (b), represented by Eqs. (32) and (1), respectively. With αmm (and β) being less sensitive to p than predicted for case (b), p is in turn more sensitive to differences or uncertainties in αmm. All else being equal, the size distribution indexes for systems with different αmm differ more strongly than expected from Eq. (1).

thumbnail Fig. 9.

Relation between the power-law indexes of the grain size distribution and the observed spectrum for four materials. Parameters and line meanings are the same as in Fig. 5. Ranges for the horizontal axes are different for the four panels, indicated above the top panels and below the bottom panels.

The curves for p(αmm) are roughly linear for 3.2 ≤ p ≤ 3.8 and can be approximated well by

p 3 + α mm α Pl β eff + Δ p $$ \begin{aligned} p \approx 3 + \frac{\alpha _{\rm mm} - \alpha _{\rm Pl}}{\beta _{\rm eff}} + \Delta p \end{aligned} $$(33)

in that range, underlining the similarity with Eqs. (32) and (1). The coefficients βeff and Δp for all assumed combinations of materials, distances, and stellar luminosities can be read from Fig. 10. The dependence on stellar luminosity is mostly weak because grains near to smin, as defined by radiation pressure blowout, contribute only slightly to mm emission. One exception is carbon with its high absorption efficiency, which shifts the characteristic grain size close to or below the lower size limit. At L = 64 L, the values for βeff are around 0.4, outside the plotted range. The other exception is the porous silicate, where n and ρ are low, the blowout limit higher, the small-grain regime more important, and the dependence on stellar luminosity more pronounced.

thumbnail Fig. 10.

Regression coefficients for linear fits to the modelled relation between observable spectral slope αmm and the size distribution index p as defined in Eq. (33). Coloured lines indicate values for longer wavelengths, grey lines for short wavelengths. Line width increases with stellar luminosity from L = 0.0625 L to L = 64 L.

The following recipe can then be used for determining the size distribution indexes from observed spectral energy distributions: (1) pick a material (mix) that is appropriate for your system; (2) derive the temperature of the big grains; (3) derive αPl for that temperature and your pair of wavelengths; (4) read βeff and Δp off of Fig. 10; and (5) calculate p from Eq. (33).

Figure 10 also shows a notable increase of βeff for temperatures below ≈30 K for all materials. This dependence on temperature is where the behaviours of αmm(p) and β(p) differ. The following section has a more detailed discussion.

3.5. Grain temperature

Grain temperature influences the relation between size distribution and SED in two ways. First, the optical properties can vary with temperature. The relation between opacity index and size distribution index thus depends on stellar luminosity and disc radius. This variability in βeff is visible in Fig. 10. The temperature-dependent water ice shows stronger variability in βs, βeff, and Δp. The slope βs of pure water ice at low temperatures is so steep that even a very steep SED could correspond to a moderate p. For carbon, the two variants of Astronomical Silicate, and pyroxene, βs does not or only weakly depend on temperature, and βeff is less variable. Second, the Planck function in Eq. (2) and the corresponding power-law index αPl depend on temperature. Figure 11 shows how αPl varies with observing wavelengths and temperature. The colder the grains and the shorter the wavelengths, the stronger is the deviation of αPl from the value of 2 found in the Rayleigh limit. Any uncertainty in the physical grain temperature therefore translates to an uncertainty in the derived size distribution index.

thumbnail Fig. 11.

Blackbody spectral indices as a function of combinations of wavelengths, normalized to temperature-dependent blackbody peak wavelengths λmax.

Figure 12 shows the equilibrium temperatures resulting from Eq. (3). For grains with radii s ≳ 100 μm temperature is nearly constant. These grains are large enough to emit efficiently because their Qabs is close to unity in the far infrared and submm. Their temperature is no longer a function of size. With the exception of very transparent materials, such as pure water ice, this temperature is close to the blackbody value (see, e. g. Pawellek et al. 2014). However, the temperature, which is assumed constant from Eq. (4) onwards, actually depends on grain size. Smaller grains, which cannot re-emit as efficiently, can be significantly hotter. The top-left panel of Fig. 12 shows an extreme case, where the smallest carbon grains at 1000 au are hotter by up to a factor of five. When the size distribution is steep enough, the contribution from small (≲100 μm) grains to the mm emission is significant. The effective blackbody spectral index, determined from the weighted average over all grain sizes, increases for higher p because of the increased abundance of smaller, hotter grains. The αmm(p) relation inherits that dependence and steepens, while p(αmm) flattens, corresponding to a higher βeff. This effects becomes more pronounced for distant discs or less luminous stars, where temperature and spectral index are lower in general and the relative spread in T is higher. Hence, the βeff curves in Fig. 10 rise as temperatures drop below ∼30 K.

thumbnail Fig. 12.

Equilibrium grain temperature as a function of grain size for four materials and four distances around a star of solar luminosity.

For spatially resolved circumstellar rings, the temperature of big grains can be computed from the distance to the star and the stellar luminosity. In the more common case, where the disc is only barely resolved or not at all, distance and temperature can only be assessed through the SED. However, this procedure is partly degenerate because the SED is already a convolute of size and temperature distributions. Determining grain temperature independently from the size distribution can be difficult. Near to the flux density peak, which is around λ ≲ 100 μm for cold debris discs, grains with s <  100 μm or even s ≪ 100 μm contribute or dominate as long as p ≳ 3. These warmer grains shift the emission peak to shorter wavelengths, masking the blackbody curve of the bigger grains that dominate at longer wavelengths. The overestimated temperatures result in overestimated αPl and underestimated p for given αmm. As shown in Fig. 11 for observations at 1 mm and 7 mm, assumed temperatures of 20 K or 50 K instead of 10 K would result in αPl = 1.83 or αPl = 1.93 instead of αPl = 1.64. For βeff ≈ 1.3 and a given αmm, this would correspond to errors in p of −0.15 or −0.22, respectively. The SED analysis can be further complicated by the different temperatures involved for debris or protoplanetary discs that are very extended or composed of two or more separate rings of comparable luminosity.

3.6. Grain porosity

To illustrate the general impact of porosity, ddscat in version 7.3.2 (Draine & Flatau 1994, 2013) is used to model spherical grains of Astronomical Silicate with filling factors f = 1, 0.75, 0.5, and 0.25, corresponding to porosities of 0, 25, 50, and 75 %. The grain radius is fixed at s = 50 μm. The size of the individual dipoles was is to 0.5 μm, much smaller than the wavelengths in the considered range 30 μm <  λ <  3 mm. The resulting Qabs(λ) curves in Fig. 13 show a trend towards lower Qabs values for lower filling factors over a wide range of wavelengths, consistent with lower absorption for less dense material. At long wavelengths, the combined dependence of Qabs on n and k given in Eq. (12) is reflected. At shorter wavelengths, in the opaque regime, Qabs is more affected by n. There, lower filling factors result in lower n, reduced reflection, and hence, increased absorption.

thumbnail Fig. 13.

Absorption efficiencies for porous spheres of Astronomical Silicate with radii s = 50 μm and filling factors f = 1, 0.75, 0.5, 0.25 (top to bottom at λ ≈ 400 μm).

The transition from (s) small and transparent to (o) large and opaque is also affected by porosity. For f = 1, n is high enough for a direct transition from regime (o) to (s). For lower f and lower n the intermediate regime (t) becomes notable. With decreasing f, the transition from (s) to (t) occurs at decreasing wavelengths because λ ∝ n for constant s c $ s_{\rm c}^{\prime} $ (see Eq. (15)). The transition from (t) to (o) occurs at λ ∝ k/n, i. e. λ decreases with decreasing k, despite n levelling off at unity.

For f → 0, as n → 1, regime (t) becomes indistinguishable from regime (s) because refraction becomes unimportant and the off-setting factor y approaches unity. The β(p) relation for porous materials is therefore expected to follow Eq. (30) more closely compared to their compact counterparts, in line with the steeper curves for lower n seen in Fig. 7. As a result, porous grains exhibit higher values of βeff and αmm when compared to compact grains with the same size distribution index p. This effect has been shown by Brunngräber et al. (2017) and is visible in Fig. 10. At the same time, the lower bulk density of porous material increases the blowout size, and hence, strengthens the dependence on disc temperature and stellar luminosity. The coefficient βeff does therefore not converge towards βs, but keeps increasing for increasing porosity (cf. Brunngräber et al. 2017).

3.7. Application to observed debris discs

Figure 14 shows the size distribution indexes resulting for a set of debris discs with observed αmm from Eq. (1) and from full Mie modelling, ordered by stellar luminosity. Table 2 lists the values used for αmm, L, disc radius r, wavelengths λ1 and λ2, together with the references. Temperatures, T, and black-body spectral indexes, αPl, are computed for mm-sized grains of astronomical silicate. This comparison shows three main effects.

thumbnail Fig. 14.

Size distribution indexes as derived from observed spectral indices for the discs presented in Table 2. Grey plusses and crosses mark the results obtained with Eq. (1) for βs = 1.8 and βs = 1.6, respectively. Black squares, filled and open red circles, green upward and blue downward triangles indicate results obtained with numerical Mie modelling for amorphous carbon, compact and porous Astronomical Silicate, pyroxene, and amorphous water ice, respectively. Where symbols are missing, the observed αmm could not be reproduced within the range 2.6 <  p <  4.2.

Table 2.

Observed spectral indexes for a set of debris discs.

First, the spread in size distribution indexes for a set of discs is wider for the more accurate numerical model than for Eq. (1). Because typically 1 <  βeff <  βs, Eq. (1) tends to underestimate the sensitivity of p on αmm, while Eq. (32) overestimates it. The numerical model falls in between. Extreme values near to or beyond the boundaries of the range 3 <  p <  4 are affected even more because there the p(αmm) relation deviates more strongly from linearity. The wider spread not only can make outliers more significant but also any potential trend, such as p(L*) or p(r).

Second, there is a systematic offset towards higher p values in some of the new results. The lines defined by Eq. (1) and the numerical models typically cross around p = 3.1…3.3 for the two compact silicates. With a mean ⟨p⟩ = 3.36 ± 0.03 for Eq. (1) and βs = 1.8 (or ⟨p⟩ = 3.41 ± 0.03 for βs = 1.6), most values are already higher than that. Hence, the increased spread also slightly increases the mean to ⟨p⟩ = 3.42 ± 0.04 for astronomical silicate and ⟨p⟩ = 3.57 ± 0.05 for pyroxene. For the porous silicate and the given set of discs, the mean (3.39) value is very close to that obtained for the simple βeff = βs = 1.6 relation.

Third, uncertainties in the composition of the dust cause uncertainties in both the mean of and the scatter among p values. When pyroxene is assumed instead of astronomical silicate, for example, ⟨p⟩ is higher by 0.15 because the slope of the imaginary part of the refractive index of pyroxene is shallower.

4. Discussion

Most of the results from the previous sections concern the discrepancies between a full numerical model and an analytic approximation based on power laws. The approximations leading to Eq. (1) tend to underestimate the mean slope and the spread of slopes of the grain size distributions for given observed spectral indexes. Assuming compact grains, the true size distribution slopes are somewhat steeper on average, which has implications for the processes that shape the size distribution. Values for p >  3.5 can be explained with details of the collisional cascade, such as specific material strength Q D $ Q_{\rm D}^* $ decreasing and/or collision velocities increasing with grain size (O’Brien & Greenberg 2003; Pan & Schlichting 2012; Gáspár et al. 2012). Values for p <  3.5 require a relative reduction of the lifetime of smaller grains. Poynting–Robertson or stellar wind drag would have that flattening effect on the size distribution (Reidemeister et al. 2011; Wyatt et al. 2011). The same is true for erosion by sublimation or UV sputtering (e.g. Grigorieva et al. 2007; Potapov et al. 2018), and for a stronger dynamical excitation of small grains by radiation pressure (Thébault & Wu 2008). A correction towards higher values of p could thus mean that collisions are more important, the according timescales shorter than that of drag and erosion.

A steeper size distribution would automatically imply a reduced overall disc mass when extrapolating from dust to planetesimals. If p were higher by 0.1, the mass ratio between km-sized objects and mm-sized dust would be reduced by a factor of 4. This could be one of several aspects that help reduce the severity of the problem with the unrealistically high extrapolated disc masses that are derived for most of the more luminous debris discs (e.g. Krivov et al. 2018).

Additional uncertainties arise from a poor understanding of the dust composition and, hence, of its optical properties. On the one hand, the real part of the refractive index, n, determines the slope 1/βeff of the p(αmm) dependence by defining the absorption and emission behaviour around the transition from transparent to opaque. For materials with higher n, such as non-porous amorphous carbon or silicates, absorption of mid-sized grains is higher. The transition from small to opaque is more abrupt, corresponding to the two characteristic radii, s c $ s_{\rm c}^{\prime} $ and sc, being closer to one another. The transparent regime is more distinguished from the small-grain regime, but narrower. The coefficient βeff is closer to unity. For materials with n closer to unity, such as water ice or porous silicate, the transparent intermediate regime is less pronounced, but wider, and hence βeff closer to βs, the power-law index of the volume absorption coefficient A. The fidelity of Eq. (1) thus depends on the material. Closer to the stars, where high-n silicates are expected to dominate, the deviation that comes with approximation (1) should be higher. Far from the stars, where low-n volatile ices can exist in abundance, the deviation should be lower. Figure 10 shows that the parameter region where Eq. (1) (with constants βs and αPl) is fully applicable is, however, very small.

On the other hand, βs is itself a material property and thus depends directly on, for example, the silicate-to-ice ratio. Measurement and usage of better optical data in the mm regime is required. Some of the commonly used data sets have never been designed to be used at such long wavelengths. The ice data set that was derived from a model fit to the β Pictoris disc (Li & Greenberg 1998), for example, inherits the uncertainties in the assumed size distribution. Figure 2 shows that the adopted k ∝ λ−0.5 is flatter than the slope actually measured for ice at a few tens of K, leading to an overestimation of the size distribution index p. Temperature dependence in optical data at mm wavelengths has been discussed for some time (e.g. Mishima et al. 1983) but only recent observations at higher resolution and sensitivity have brought the issue back into focus. Consequently, temperature-dependent data for more and more materials are becoming available (Reinert et al. 2015; Häßner et al. 2018; Mutschke & Mohr 2019).

Strong uncertainties remain even if a full Mie model is used and the grain composition is known. Two of the very basic assumptions, that of solid, spherical grains and that of a power-law size distribution, can be too crude. As shown in Sect. 3.6, the transition from absorption by cross-section to absorption by volume is smoother for filling factors f ≲ 0.5, which is more in line with the assumptions that lead to Eq. (1). The lower bulk density and the lower refractive index, however, cause smaller grains and the lower size cutoff imposed by radiation pressure to become more important. The resulting relation between observed αmm and deduced p flattens (corresponding to a higher βeff) for lower disc temperatures and stellar luminosities. If a constant βeff = βs were imposed for porous grains, the size distribution index would be underestimated for less luminous stars or distant discs and overestimated for more luminous stars or closer discs, causing an artificial trend in disc statistics. If, indeed, porosity is significant, all results based on Mie theory of compact spheres end up suffering from a systematic error. Brunngräber et al. (2017, their Fig. 3) show that the size distribution index for a filling factor of f = 0.5, derived based on the assumption of solid grains, is already too high by 10 to 15%, or, in absolute terms, 0.35 to 0.5 higher than their reference value of 3.5. However, their curves do not converge for even higher porosities (lower filling factors), indicating that there is no simple porosity threshold above which Eq. (1) is applicable and below which it is not.

Dust size distributions may deviate strongly from power laws when radiation pressure blowout introduces a discontinuity in the collisional cascade, resulting in a wavy pattern (Campo Bagatin et al. 1994; Durda & Dermott 1997; Thébault et al. 2003; Krivov et al. 2006). Parts of that waviness are carried over to the total opacity and the observable SED (Thébault & Augereau 2007). The narrower the range of grain sizes that effectively contribute to emission at any given wavelength, the more of the waviness is preserved. This means that a wavier, more expository SED is produced by the more abrupt transition from surface absorption to volume absorption that is seen for materials of higher optical density. For materials or wavelengths with a lower opacity, the convolution over a wider range of grain sizes reduces the waviness. In addition, in dynamically cold debris discs, smaller grains are under-abundant because they suffer from dynamical excitation due to radiation pressure (Thébault & Wu 2008). If the transition from affected small to unaffected big grains occurs in the (sub)mm size range, the average slope derived from the SED analysis will be shallower than the background power-law slope of an infinite, unperturbed collisional cascade.

Models ought to take into account that the characteristic grain size sc can easily exceed λ by orders of magnitude (cf. Draine 2006): sc ≈ λ/(16k) for n ≈ 1 and s c < s c $ s_{\rm c}^{\prime} < s_{\rm c} $ in Eq. (16). Assuming a moderate k ∼ 10−4 for water ice at λ close to 1 cm, sc is already in the order of 600λ, that is, several metres. For the more opaque compact carbon and silicates, where k ≳ 0.01 (Fig. 2), the result is: sc ≲ 6λ. A porous blend of these material can be expected to fall in between these extremes. With the SED being dominated by grains with radii around sc for 3 <  p <  4, the maximum grain sizes smax in numerical and analytic models should always be set such that smax ≫ sc or reflect a physically motivated maximum. If the assumed smax is too small, say 1 mm, a lack of flux results at the longer wavelengths. The resulting synthetic SED will be too steep, requiring a flatter-than-actual size distribution when fit to observational data. On the same note, any dust mass derived from mm or cm observations is a mass in cm or dm-sized grains (if these are present), which sets the starting point from which to extrapolate towards a total disc mass. The often quoted mass in mm grains and the total disc mass would then be lower by a factor of (0.01…0.1)p − 3 ≈ 0.1…0.3 for p = 3.5 (cf. Krivov et al. 2018).

The large sc and smax imply that observations in the mm-cm wavelength range are sensitive to much bigger, dm-sized grains. Any derived size distribution index would also cover that size range, including potential changes in slope or cutoffs. Models for planetesimal formation through pebble aggregation, such as the streaming instability scenario (e. g. Johansen et al. 2015), predict a transition from stronger individual pebbles to weaker and, hence, less abundant pebble piles around cm sizes (cf. Krivov et al. 2018). This transition would translate to a steepening of the SED in the wavelength range of interest. Below and above the transition size from pebbles to piles, where collision physics are assumed continuous, the size distribution would be flatter. At earlier phases, when the pebbles are still forming from smaller grains, smax is a measure of the actual maximum grain size. The opacity curves and observed SEDs are steeper early on, flattening as smax grows with time (e.g. Krügel & Siebenmorgen 1994; Natta & Testi 2004; Draine 2006). At wavelengths for which optical depths ≳1 are reached, Eq. (2) is not applicable and, thus, the relations for p(αmm) derived in Sects. 2 and 3 are invalidated.

5. Summary

When analysing SEDs of individual systems, a full numerical approach is the most flexible in terms of materials and assumed size distributions. The computational cost of Mie calculations is low. If only a quick estimate of the size distribution index and its uncertainties is sought, Eq. (33) can be used in combination with the coefficients provided in Fig. 10. These coefficients are fitted to the numerical results for moderate size distributions with 3.2 <  p <  3.8 for five materials, two wavelength combinations, and a range of stellar luminosities and disc temperatures.

With the exception of the unrealistic pure water ice, the numerical fits to observed systems result in steeper size distribution slopes, on average, tipping the scales more towards size-dependent material strength or impact velocities and, then, away from drag forces. The system-to-system variations are somewhat more pronounced than previously thought. However, these findings are accompanied by strong uncertainties related to the material composition and porosity. Further work is needed to examine the impact of realistic material blends, grain structures, and temperature-dependent optical properties at long wavelengths.

Acknowledgments

TL is grateful to H. Mutschke for fruitful discussions on dust optical properties. Parts of this work have been supported by the Deutsche Forschungsgemeinschaft (grants LO 1715/2-1 and MU 1164/9-1 within Research Unit FOR 2285 “Debris Disks in Planetary Systems”).

References

  1. Bohren, C. F., & Huffman, D. R. 1983, Absorption and Scattering of Light by Small Particles (New York: Wiley) [Google Scholar]
  2. Booth, M., Dent, W. R. F., Jordán, A., et al. 2017, MNRAS, 469, 3200 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bruggeman, D. A. G. 1935, Ann. Phys., 416, 636 [NASA ADS] [CrossRef] [Google Scholar]
  4. Brunngräber, R., Wolf, S., Kirchschlager, F., & Ertel, S. 2017, MNRAS, 464, 4383 [NASA ADS] [CrossRef] [Google Scholar]
  5. Campo Bagatin, A., Cellino, A., Davis, D. R., Farinella, P., & Paolicchi, P. 1994, Planet. Space Sci., 42, 1079 [NASA ADS] [CrossRef] [Google Scholar]
  6. Dohnanyi, J. S. 1969, J. Geophys. Res., 74, 2531 [NASA ADS] [CrossRef] [Google Scholar]
  7. Dominik, C., Blum, J., Cuzzi, J. N., & Wurm, G. 2007, Protostars and Planets V, eds. B. Reipurth, D. Jewitt, & K. Keil, 783 [Google Scholar]
  8. Draine, B. T. 2006, ApJ, 636, 1114 [Google Scholar]
  9. Draine, B. T., & Flatau, P. J. 1994, J. Opt. Soc. Am. A, 11, 1491 [NASA ADS] [CrossRef] [Google Scholar]
  10. Draine, B. T., & Flatau, P. J. 2013, ArXiv e-prints [arXiv:1305.6497] [Google Scholar]
  11. Durda, D. D., & Dermott, S. F. 1997, Icarus, 130, 140 [NASA ADS] [CrossRef] [Google Scholar]
  12. Gaia Collaboration (Prusti, T., et al.) 2016, A&A, 595, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Gáspár, A., Psaltis, D., Rieke, G. H., & Özel, F. 2012, ApJ, 754, 74 [NASA ADS] [CrossRef] [Google Scholar]
  15. Gerbaldi, M., Faraggiana, R., Burnage, R., et al. 1999, A&AS, 137, 273 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Grigorieva, A., Thébault, P., Artymowicz, P., & Brandeker, A. 2007, A&A, 475, 755 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Häßner, D., Mutschke, H., Blum, J., Zeller, T., & Gundlach, B. 2018, MNRAS [Google Scholar]
  18. Hughes, A. M., Lieman-Sifry, J., Flaherty, K. M., et al. 2017, ApJ, 839, 86 [NASA ADS] [CrossRef] [Google Scholar]
  19. Hung, L.-W., Duchêne, G., Arriaga, P., et al. 2015a, ApJ, 815, L14 [NASA ADS] [CrossRef] [Google Scholar]
  20. Hung, L.-W., Fitzgerald, M. P., Chen, C. H., et al. 2015b, ApJ, 802, 138 [NASA ADS] [CrossRef] [Google Scholar]
  21. Johansen, A., Mac Low, M.-M., Lacerda, P., & Bizzarro, M. 2015, Sci. Adv., 1, 1500109 [NASA ADS] [CrossRef] [Google Scholar]
  22. Kataoka, A., Okuzumi, S., Tanaka, H., & Nomura, H. 2014, A&A, 568, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Kennedy, G. M., Marino, S., Matrà, L., et al. 2018, MNRAS, 475, 4924 [Google Scholar]
  24. Krivov, A. V. 2010, Res. Astron. Astrophys., 10, 383 [NASA ADS] [CrossRef] [Google Scholar]
  25. Krivov, A. V., Löhne, T., & Sremčević, M. 2006, A&A, 455, 509 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Krivov, A. V., Müller, S., Löhne, T., & Mutschke, H. 2008, ApJ, 687, 608 [NASA ADS] [CrossRef] [Google Scholar]
  27. Krivov, A. V., Ide, A., Löhne, T., Johansen, A., & Blum, J. 2018, MNRAS, 474, 2564 [NASA ADS] [CrossRef] [Google Scholar]
  28. Krügel, E., & Siebenmorgen, R. 1994, A&A, 282, 407 [NASA ADS] [Google Scholar]
  29. Li, A., & Greenberg, J. M. 1998, A&A, 331, 291 [NASA ADS] [Google Scholar]
  30. Liseau, R., Eiroa, C., Fedele, D., et al. 2010, A&A, 518, L132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. MacGregor, M. A., Wilner, D. J., Andrews, S. M., Lestrade, J.-F., & Maddison, S. 2015, ApJ, 809, 47 [NASA ADS] [CrossRef] [Google Scholar]
  32. MacGregor, M. A., Wilner, D. J., Chandler, C., et al. 2016, ApJ, 823, 79 [Google Scholar]
  33. MacGregor, M. A., Matrà, L., Kalas, P., et al. 2017, ApJ, 842, 8 [NASA ADS] [CrossRef] [Google Scholar]
  34. MacGregor, M. A., Weinberger, A. J., Hughes, A. M., et al. 2018, ApJ, 869, 75 [NASA ADS] [CrossRef] [Google Scholar]
  35. MacGregor, M. A., Weinberger, A. J., Nesvold, E. R., et al. 2019, ApJ, 877, L32 [CrossRef] [Google Scholar]
  36. Marino, S., Matrà, L., Stark, C., et al. 2016, MNRAS, 460, 2933 [NASA ADS] [CrossRef] [Google Scholar]
  37. Marino, S., Wyatt, M. C., Panić, O., et al. 2017, MNRAS, 465, 2595 [NASA ADS] [CrossRef] [Google Scholar]
  38. Marino, S., Carpenter, J., Wyatt, M. C., et al. 2018, MNRAS, 479, 5423 [NASA ADS] [CrossRef] [Google Scholar]
  39. Marshall, J. P., Maddison, S. T., Thilliez, E., et al. 2017, MNRAS, 468, 2719 [CrossRef] [Google Scholar]
  40. Marshall, J. P., Milli, J., Choquet, É., et al. 2018, ApJ, 869, 10 [CrossRef] [Google Scholar]
  41. Matrà, L., Wyatt, M. C., Wilner, D. J., et al. 2019, AJ, 157, 135 [NASA ADS] [CrossRef] [Google Scholar]
  42. Matthews, B., Kennedy, G., Sibthorpe, B., et al. 2014, ApJ, 780, 97 [NASA ADS] [CrossRef] [Google Scholar]
  43. Matthews, B. C., Kennedy, G., Sibthorpe, B., et al. 2015, ApJ, 811, 100 [CrossRef] [Google Scholar]
  44. Mishima, O., Klug, D. D., & Whalley, E. 1983, J. Chem. Phys., 78, 6399 [NASA ADS] [CrossRef] [Google Scholar]
  45. Mutschke, H., & Mohr, P. 2019, A&A, 625, A61 [CrossRef] [EDP Sciences] [Google Scholar]
  46. Natta, A., & Testi, L. 2004, Star Formation in the Interstellar Medium. In Honor of David Hollenbach, eds. D. Johnstone, F. C. Adams, D. N. C. Lin, D. A. Neufeeld, & E. C. Ostriker, ASP Conf. Ser., 323, 279 [Google Scholar]
  47. O’Brien, D. P., & Greenberg, R. 2003, Icarus, 164, 334 [CrossRef] [Google Scholar]
  48. Pan, M., & Schlichting, H. E. 2012, ApJ, 747, 113 [NASA ADS] [CrossRef] [Google Scholar]
  49. Pawellek, N., Krivov, A. V., Marshall, J. P., et al. 2014, ApJ, 792, 65 [NASA ADS] [CrossRef] [Google Scholar]
  50. Potapov, A., Jäger, C., & Henning, T. 2018, ApJ, 865, 58 [NASA ADS] [CrossRef] [Google Scholar]
  51. Reidemeister, M., Krivov, A. V., Stark, C. C., et al. 2011, A&A, 527, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Reinert, C., Mutschke, H., Krivov, A. V., Löhne, T., & Mohr, P. 2015, A&A, 573, A29 [CrossRef] [EDP Sciences] [Google Scholar]
  53. Ricci, L., Testi, L., Maddison, S. T., & Wilner, D. J. 2012, A&A, 539, L6 [CrossRef] [EDP Sciences] [Google Scholar]
  54. Ricci, L., Maddison, S. T., Wilner, D., et al. 2015, ApJ, 813, 138 [CrossRef] [Google Scholar]
  55. Steele, A., Hughes, A. M., Carpenter, J., et al. 2016, ApJ, 816, 27 [NASA ADS] [CrossRef] [Google Scholar]
  56. Su, K. Y. L., MacGregor, M. A., Booth, M., et al. 2017, AJ, 154, 225 [NASA ADS] [CrossRef] [Google Scholar]
  57. Thébault, P., & Augereau, J.-C. 2007, A&A, 472, 169 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Thébault, P., & Wu, Y. 2008, A&A, 481, 713 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Thébault, P., Augereau, J.-C., & Beust, H. 2003, A&A, 408, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Wolf, S., & Voshchinnikov, N. V. 2004, Comput. Phys. Commun., 162, 113 [NASA ADS] [CrossRef] [Google Scholar]
  61. Wyatt, M. C., Clarke, C. J., & Booth, M. 2011, Celestial Mech. Dyn. Astron., 111, 1 [NASA ADS] [CrossRef] [Google Scholar]
  62. Zubko, V. G., Mennella, V., Colangeli, L., & Bussoletti, E. 1996, MNRAS, 282, 1321 [Google Scholar]

All Tables

Table 1.

Alphabetic list of recurring symbols.

Table 2.

Observed spectral indexes for a set of debris discs.

All Figures

thumbnail Fig. 1.

Schematic comparison of the two different cases (a) and (b) for the transition from small grains that absorb by volume to big, opaque grains that absorb by cross-section. Case (a) has an intermediate size interval where the grains are already big compared to the wavelength but still transparent. In case (b) the transition from small to opaque is direct. The values assumed in Eqs. (17), (18), and (28) are βs = 1.6, A0 = 0.63 mm−1, λ0 = 1 mm, n = 3, y = n2.

In the text
thumbnail Fig. 2.

Refractive indexes of the five materials discussed: (top) real and (middle, bottom) imaginary parts. The legends are ordered vertically according to the order at long wavelengths. The vertical dashed lines delimit the wavelength range over which the labelled slopes are fitted.

In the text
thumbnail Fig. 3.

Absorption efficiencies for four materials and a range of wavelengths and grain sizes. Wavelength are separated by factors of 10, the longest wavelengths are labelled. Solid lines indicate Mie results, dashed lines the small-grain asymptote given in Eq. (12). The upper boundaries of the light-shaded regions are defined by Eq. (13), the left boundaries by Eq. (15). Black dotted lines follow Eq. (13) multiplied by an additional factor n to roughly trace the actual smoothed Qabs. Values are given for stellar luminosity L = LSun and a distance r = 100 au.

In the text
thumbnail Fig. 4.

Opacities for four different materials for a range of size distribution slopes and wavelengths. Values are given for stellar luminosity L = L, a distance r = 100 au, and a total disc mass Md = 0.1MEarth. The maximum object radius is set to smax = 100 km for water ice and 1 km for carbon and the silicates. Mie results are plotted with solid lines for five size distribution indexes from (top) p = 4.2 to (bottom) p = 2.6. Sloped dashed lines trace power-law fits to wavelengths from λ1 = 1 mm to the shorter of λ2 = 7 mm and the maximum wavelength in the tabulated optical data. Vertical dashed lines indicate λ1 and λ2.

In the text
thumbnail Fig. 5.

Relation between the power-law indices of the grain size distributions and the total opacities between λ1 = 1 mm and λ2 = 7 mm for four materials at r = 100 au around a star with L = L. Solid coloured lines show numerical results. Dotted lines of the same colours show linear fits for 3.2 ≤ p ≤ 3.8, with fit coefficients indicated in the figure keys. Dashed lines show the analytic approximation according to Eq. (30). The solid grey lines trace the approximate Eq. (31).

In the text
thumbnail Fig. 6.

Same as the top-right panel in Fig. 5 (astronomical silicate) but for two combinations of shorter wavelengths, as indicated.

In the text
thumbnail Fig. 7.

Slope of the opacity as derived from Eq. (25) for the following parameters: λ = 3 mm, βs = 1.6, A0 = 2/cm, λ0 = 1 mm, ρ = 3.5 g cm−3, smin = 0.5 μm, smax = 100 m, Qabs, o = 1, y = 0.75n2, with (black) n = 8.5, (red) n = 3.4, and (blue) n = 1.14. In the range 3 <  p <  4, this order of colours reflects the vertical order of the corresponding lines, bottom to top. The dashed black and solid grey lines follow Eqs. (30) and (31), respectively.

In the text
thumbnail Fig. 8.

Same as Fig. 7, but with (top left) Qabs, o, (top right) smax, (middle) smin, and (bottom) A0 varied. Left panels: results for high n, right panels: for low n. The legends are ordered vertically according to the order at p ≈ 4.

In the text
thumbnail Fig. 9.

Relation between the power-law indexes of the grain size distribution and the observed spectrum for four materials. Parameters and line meanings are the same as in Fig. 5. Ranges for the horizontal axes are different for the four panels, indicated above the top panels and below the bottom panels.

In the text
thumbnail Fig. 10.

Regression coefficients for linear fits to the modelled relation between observable spectral slope αmm and the size distribution index p as defined in Eq. (33). Coloured lines indicate values for longer wavelengths, grey lines for short wavelengths. Line width increases with stellar luminosity from L = 0.0625 L to L = 64 L.

In the text
thumbnail Fig. 11.

Blackbody spectral indices as a function of combinations of wavelengths, normalized to temperature-dependent blackbody peak wavelengths λmax.

In the text
thumbnail Fig. 12.

Equilibrium grain temperature as a function of grain size for four materials and four distances around a star of solar luminosity.

In the text
thumbnail Fig. 13.

Absorption efficiencies for porous spheres of Astronomical Silicate with radii s = 50 μm and filling factors f = 1, 0.75, 0.5, 0.25 (top to bottom at λ ≈ 400 μm).

In the text
thumbnail Fig. 14.

Size distribution indexes as derived from observed spectral indices for the discs presented in Table 2. Grey plusses and crosses mark the results obtained with Eq. (1) for βs = 1.8 and βs = 1.6, respectively. Black squares, filled and open red circles, green upward and blue downward triangles indicate results obtained with numerical Mie modelling for amorphous carbon, compact and porous Astronomical Silicate, pyroxene, and amorphous water ice, respectively. Where symbols are missing, the observed αmm could not be reproduced within the range 2.6 <  p <  4.2.

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.