Impact of discontinuous grain size distributions on the spectral energy distribution of debris disks

The collisional evolution of debris disks is expected to result in a characteristic wavy pattern of the grain size distributions, i.e., an under/overabundance of particles of specific sizes. This perturbed grain size distribution potentially leaves characteristic patterns in the spectral energy distribution (SED) of the disk system. We aim to quantify and understand the specific influence of discontinuous particle size distributions on the appearance of debris disks. For this purpose, we consider dust emission models based on two different grain size distributions, i.e., once with a single and once with a broken power law. We compare the spectral index $\alpha$ ($F_{\nu}\,\propto \nu^{\rm{\,\alpha}}$) in the case of a continuous grain size distribution with that of a discontinuous grain size distribution. We perform this comparison for central stars with different spectral types and two different disk structures (e.g., slim and broad debris dust rings). Within the considered parameter space, we find a characteristic difference between the spectral slopes of the SED in the different scenarios. More specifically, the overabundance of small grains leads to a steeper slope in the far-infrared/submillimeter regime, while the spectral index in the mm regime is hardly affected. On the other hand, the underabundance of medium-sized grains results in a slight steepening of the far-infrared slope of SED, while its primary effect is on the mm slope of SED, causing it to become shallower. We also find that the impact of an overabundance of small dust particles is more pronounced than that of an underabundance of medium-sized dust particles. Furthermore, we find that the difference between the spectral indices for the two different grain size distributions is largest for debris disks around brighter central stars and broader disks.


Introduction
The lifetime of small dust grains in debris disks is much shorter than the stellar age (Artymowicz & Clampin 1997), indicating that dust cannot be primordial in these systems.Thus, the dust in debris disks must be continuously replenished through disruptive collisions.This process leads to the creation of fragments with a broad size distribution, which can be approximately described by a power law.For example, for a steady-state collisional cascade (Dohnanyi 1969), a power law with γ = −3.5 is valid for the differential size distribution, namely, dn (a) ∝ a γ da, where a is the size of the dust grain.However, due to the finite size distribution and size-dependent collisional strength of the colliding bodies and particles, the real slope of the size distribution may strongly deviate from this value and thus may not accurately reflect the true grain size distribution.For example, the collisional cascade does not extend down to arbitrarily small particles because tiny particles with a high β ratio between radiation pressure and stellar gravity are quickly blown away by the stellar radiation pressure from the disks (Krivov et al. 2006).Furthermore, non-gravitational forces (e.g., Poynting-Robertson drag; Backman & Paresce 1993) acting on grains of particular sizes in the range of tens to hundreds of micrometers may further modify the grain size distribution (Löhne et al. 2017).
In previous studies based on numerical simulations, a characteristic wavy pattern in the grain size distribution resulting from collisional evolution within debris disks has been found (e.g., Thébault & Augereau 2007;Kim et al. 2018).These patterns, which are not solely described by a simple power law but, for example, by broken power laws, indicate the overabundance or underabundance of particles of specific sizes.In particular, the wavy pattern of the discontinuous grain size distribution becomes more pronounced (deeper) and broader in extreme scenarios involving the lower material strength of dust grains (Kim et al. 2018).
Furthermore, several observational findings have pointed toward the need for a more complex model for grain size distribution beyond a simple single power law.For instance, Table 1.Disk and dust model parameters for the simulations and corresponding references.
Table 2. Stellar, disk, and dust model parameters for the simulations and their reference.
observational findings based on the comparison between flux and polarization measurements of protoplanetary disks at submillimeter wavelengths have also raised questions regarding the corresponding grain size distribution (e.g., Brunngräber & Wolf 2021).Bürger et al. (2023) also recently demonstrated an apparent overabundance of grains with radii between 3 and 6 mm in the derived grain size distribution of comet 67P/Churyumov-Gerasimenko. Additionally, Rodigas et al. (2015) found that the grain size distributions derived from the observation of thermal reemission radiation are not consistent with those from scattered light colors and/or scattering functions.Finally, Olofsson et al. (2023) recently demonstrated that the overabundance of small grains in debris disks could also explain an intriguing arc-like structure observed in total intensity.These findings suggest the presence of a discontinuous or more complicated grain size distribution and challenge the traditional single power law model.The grain size distribution directly influences the slope of the spectral energy distribution (SED) in the far-infrared (far-IR) to (sub-)millimeter range of the thermal reemission by the dust (Wyatt & Dent 2002).Consequently, the presence of overabundant or underabundant particle sizes at these wavelengths potentially leaves an observable imprint in the SED.For example, the recent NIKA2 observation of HD 107146 indicates a break at the far-IR and sub-millimeter/millimeter wavelength regimes (Lestrade et al. 2020), which is further supported by supplementary observations, possibly showing a change in the slope of the SED (e.g., Holland et al. 2017;Ricci et al. 2015b).In this particular case, the spatial distribution of dust in this debris disk system cannot be a reason for this break at long wavelength because HD 107146 is observed as a narrow belt system (Marino et al. 2018).Moreover, the influence of the dust temperature on the Planck function is weak in the Rayleigh-Jeans regime, diminishing its importance at long wavelengths.However, given the small range of dust temperatures in a narrow belt system, the impact of the underlying grain size distribution becomes more pronounced.Consequently, a change in the SED slope due to an overabundance/underabundance of grains is expected in the wavelength region where grains of the corresponding size emit most efficiently.
The goal of the current study is to evaluate the impact of the discontinuous size distribution of dust in debris disks on the SED of these systems.In turn, the results of this study are expected to serve as a tool to also address the opposite question, namely whether constraints on the parameters describing the collisions in debris disks can be derived from the analysis of the slope of their thermal reemission SED.The present study is structured as follows: in Sect.2, we describe our model for the simulations of debris disk observables.In Sect.3, we quantify and discuss the impact of the grain size distributions on the observational appearance of debris disks, in particular, their SEDs (e.g., the slope of SEDs).In Sect.4, we discuss an application of our results to observations and possible degeneracies as well as connections to chemical and physical properties.We summarize our findings in Sect. 5.

Model parameters of debris disk system
The parameters and characteristics of stars, disks, and dust in our models are based on observations of currently known debris disk systems.Tables 1 and 2 summarize the detailed characteristics of disk and dust parameters for our simulations.
Disk geometry.We considered the inner and outer radii of disks based on the empirical relation between planetesimal belt A148, page 2 of 14 Kim, M., and Wolf, S.: A&A, 683, A148 (2024)  Standard abundance of large grains 1.00 × 10 −7 SL1 dn (a) ∝ a −3.5 da for [1000 µm, 10 cm] radius and stellar luminosity (Matrà et al. 2018) with two different belt fractional widths: a narrower belt (e.g., ∆R/R ∼ 0.17 of debris disks around ε Eridani; Lestrade & Thilliez 2015) comparable to the classical Kuiper belt and a broader belt (e.g., ∆R/R ∼ 1.25 of debris disks around HD 206893; Marino et al. 2020).We note that broader belt systems show a broad range of dust temperatures.The mass of inner belts in a two-belt system is typically only a few percent compared to outer belts (e.g., Ertel et al. 2011;Marino et al. 2018), implying that the contribution of the inner belt to the slope of the SED is negligible.Thus, we considered a single belt system to reduce the complexity of the model.
Minimum/maximum grain size and dust composition.For the smallest grain size, we applied a blowout grain size a bo (Backman & Paresce 1993) depending on the spectral type.For cases without a blowout size (i.e., for central stars' effective temperatures below 5250 K in our model setup; Kirchschlager & Wolf 2013), we considered a minimum grain size of 0.264 µm that has been used in previous studies (e.g., Löhne et al. 2017;Kim et al. 2018), which shows the similarity to the previous findings of debris disks around M-type stars (e.g., Matthews et al. 2015) based on the strong stellar wind (Plavchan et al. 2005), causing an effect similar to that of the radiation pressure force (Augereau & Beust 2006).Furthermore, using Mie theory, silicate or ice grains comparable to this size are hardly affected by the stellar radiation pressure in debris disks around a solartype star (Cataldi 2016).We adopted a maximum grain radius of 10 cm (Lestrade et al. 2020;Pawellek et al. 2021) to analyze the dust continuum up to a wavelength of a few millimeters (e.g., ATCA Survey; Ricci et al. 2015b).We note that the contribution of larger grains to the net flux is negligible when considering a steep grain size distribution.For the grain composition, we assumed the simplest dust grains to be compact and homogenous spheres composed of astronomical silicate (referred to as "Astrosil") with a bulk density of 3.5 g cm −3 and optical constants taken from Draine (2003).
Disk mass.Based on previous surveys at sub-millimeter wavelengths, we considered a dust mass in debris disks of ∼10 −7 M ⊙ for the single power law grain size distribution model Overabundance of small grains Factor 10 (OS10) Factor 50 (OS50) Factor 100 (OS100) Underabundance of medium grains Factor 10 (UM10) Factor 50 (UM50) Factor 100 (UM100) Single power law (reference) Fig. 1.Simulated total cross section of the dust in our debris disk around the A star (e.g., a bo ∼ 1.99 µm) based on the considered grain size distributions.The red solid line represents the classical single power law (reference case; see Sect.2.2.1).The yellow, green, and blue dashed lines represent the perturbed size distribution with an overabundance of grains with radii ranging from the blowout size (see Table 2) to 100 µm.The purple, gray, and navy dashed lines represent the distribution with the underabundance in the size range of 100 µm to 1000 µm, which are mainly motivated by the collisional model presented in Thébault & Augereau (2007) and Kim et al. (2018).
(e.g., dust mass of most debris disks typically range from ∼10 −9 to several 10 −7 M ⊙ ; Greaves et al. 2005).In the case of the broken power law grain size distribution model, the dust mass is varied depending on the scale factor but remains within a similar range (i.e., 1.0 to 1.26 × 10 −7 M ⊙ ; see Table 3).This is because the contribution of small dust grains mainly affects the dust reemission rather than significantly influencing the overall mass budget, which is dominated by larger grains.We considered the same mass of dust M dust for both types of disk structures (i.e., narrow and broad belt-like disks).We note that within this model, the overall mass of a disk does not significantly affect the spectral index as long as the disk remains optically thin; however, a specific range of dust masses with the under-or overabundance of medium and small grains is indeed an important factor.The detailed characteristics of the modeled classical single power law and wavy pattern (i.e., broken power law) grain size distribution are compiled in Fig. 1 and Table 3. A148, page 3 of 14 Kim, M., and Wolf, S.: A&A, 683, A148 (2024) Central star.We considered central stars of the following five spectral types (in brackets, exemplary systems are listed): A star (β Pictoris), F star (HD 139664), G star (HD 107146), K star (HD 53143), and M star (AU mic).Table 2 summarizes the detailed characteristics of the stellar parameters and corresponding disk and dust parameters.

Simulation of SEDs.
For the calculation of model SEDs, we used the software tool DMS (Kim et al. 2018), which is optimized for optically thin emission of debris disks.Based on wavelength-dependent refractive indices, we calculated the required the optical properties with the tool miex on the basis of the Mie scattering theory (Wolf & Voshchinnikov 2004).

Modeling of grain size distributions
To understand how the grain size distribution affects the SED of debris disks, we considered dust emission models that are based on two different grain size distributions: a single power law grain size distribution (Sect.2.2.1) and a broken power law grain size distribution (Sect.2.2.2).Based on these distributions, we studied the shape of the thermal reemission SED in the far-IR to the sub-millimeter-millimeter wavelength range.

Single power law grain size distribution
As a reference case, we considered the classical single power law grain size distribution dn (a) ∝ a −3.5 da (see Fig. 1 and Table 3 top row for details) resulting from a collisional cascade in debris disks at steady-state equilibrium with size-independent collisional strength (Dohnanyi 1969).Since the absorption cross section C abs of dust grains gradually decreases toward longer wavelengths irrespective of the dust composition (Kim et al. 2019), the dust emission model employing the single power law grain size distribution is expected to result in a continuous decrease in the SED at long wavelengths (i.e., the far-IR to millimeter range).

Broken power law grain size distribution
As indicated by previous studies (Thébault & Augereau 2007;Löhne et al. 2017;Kim et al. 2018), collisional evolution leaves an imprint on the grain size distribution.In particular, the depletion of the smallest grains to the blowout size results in a surplus of slightly larger grains due to reduced collisional erosion.This leads to the propagation of an overabundance around small dust grains with radii ranging from several to approximately hundreds of micrometers and an underabundance (i.e., local minimum) toward larger grains with radii of several hundred to thousands of micrometers.For example, an observed peak in the wavy pattern of the grain size distribution occurs for grains that are slightly larger than the blowout size, specifically at about ∼1.5 times this size.Additionally, the first minimum in this distribution is noted at about 100 times the blowout size (Thébault & Augereau 2007).
Furthermore, for the particular case of debris disks around HD107146, referenced in Sect. 1, a local decrease in the long wavelength regime of the SED ranging from 800 µm to 6 mm has been observed (Lestrade et al. 2020).A possible explanation for this finding is to link it to the peculiarities of the grain size distribution, where the deep depletion of the abundance of grains in the size range of hundreds of micrometers coincides with the characteristic grain size for strongly absorbing materials (∼ approximately ranging from 100 to 1000 µm; Backman & Paresce 1993).
Based on the discussion in this section, we considered a general and rather simplistic grain size distribution that shows the key characteristics that potentially allow us to explain discontinuities in the shape of the thermal reemission SED of debris disks.Our discontinuous grain size distribution models follow a broken power law pattern, comprising an overabundance of small grains ranging from the blowout grain size to 100 µm and an underabundance of medium-sized grains ("medium grain" in the following) ranging from 100 to 1000 µm (Backman & Paresce 1993).We incorporated scaling factors of 10, 50, and 100 to represent the relative abundance of these grain populations (see Fig. 1 and Table 3 for details).We did not consider the depletion of the smallest dust grains to the blowout grain size that results in the excess of slightly larger grains (e.g., Kim et al. 2018), as the effect of their absence on the spectral slope of the SED is negligible in the spectral range explored in this present study.To avoid the introduction of additional parameters, we set all slopes described by the exponent of the grain size distribution γ in the broken power law grain size distribution to -3.5 (Dohnanyi 1969) for all size ranges.We note that the location of the subsequent minima and maxima and their overabundance and underabundance may depend not only on the collisional circumstances but also on further disk parameters (e.g., Thébault & Augereau 2007;Kim et al. 2018).Thus, the location and abundance of the local minima and maxima of the size distribution may not be directly indicated by the physical situations (e.g., the collisional state of the disk).
We also note that a uniform overabundance and/or underabundance of grains of certain sizes across all spectral types may not accurately represent the actual conditions within these systems.Thus, we particularly considered grain size distributions with an overabundance and/or underabundance of certain grain sizes as a function of stellar spectral type (i.e., blowout size).For this purpose, we applied a grain size distribution model following the findings of Thébault & Augereau 2007;Kim et al. 2018, in which the first local maximum of the grain size distribution typically occurs at about 1.5 times the blowout size.Furthermore, the first minimum and second maximum are observed at around 100 and 1000 times the blowout size, respectively (see also Fig. 2 and Table 4).In the case of disks around stars that are not luminous enough to reach β = 0.5 (e.g., K and M stars in our model; Kirchschlager & Wolf 2013), we considered 0.264 µm as the minimum grain size (see also Sect.2.1).The results based on the specifically proposed grain size distribution and concerning its impact on the corresponding SEDs in particular can be found in Sect.4.4.

Results
In the following, we systematically compare and evaluate the influence of the different dust emission models based on the two different grain size distributions on the SED for the various considered debris disk systems.Specifically, we focus on quantifying the slope of the SED at long wavelengths in order to assess how the presence of over-and underabundances of the specific size of grains in the distribution impacts the emission characteristics of debris disks.For this purpose, we compute the spectral index α as where F ν represents the radiative flux density and ν denotes the frequency in the dust continuum emission.based on the considered grain size distributions with exemplary cases (i.e., 100 times overabundance of small grains and underabundance of medium grains).The dashed red, orange, yellow, green, and blue color lines represent the classical single power law (reference case) around A, F, G, K, and M stars.The perturbed size distributions with an overabundance of small grains with radii ranging from the blowout size to 100 times each blowout size and one with an underabundance of medium grains with radii ranging from 100 times to 1000 times each blowout size, depending on spectral types are shown (see Table 2).These size distributions are mainly motivated by the collisional model presented in Thébault & Augereau (2007) and Kim et al. (2018).We note that the K star and M star do not have blowout sizes, suggesting the minimum grain sizes are used often in previous studies (Löhne et al. 2017;Kim et al. 2018).Thus, both grain size distributions look the same.See Table 4 for details.
Considering our specific grain size distribution with overabundant small grains and a deep depletion of grains in the size range of hundreds of micrometers, we expected spectral transitions to occur at sub-millimeter to millimeter wavelengths.To characterize the SED in this wavelength range, we calculated the spectral index within two distinct wavelength ranges: 450/850 µm (α 450 µm 850 µm ) and 1000/2000 µm (α 1000 µm 2000 µm ).These spectral indices enabled us to quantify the differences between the two wavelength ranges represented by ∆ α = α 450 µm 850 µmα 1000 µm 2000 µm .Positive and negative values of ∆ α respectively indicate a lower and higher value of α 1000 µm 2000 µm , representing a shallower and steeper slope in the 1000 to 2000 µm range compared to the 450-850 µm range, respectively.

Single power law grain size distribution
We began by examining the spectral index of the SED for the continuous grain size distribution described by the single power law (i.e., γ = −3.5).This reference case allowed us to identify the influence of the over and underabundance of grains of specific sizes.The top-left plots of Figs. 3 and 4 show the resulting SED for the two different disk structures (i.e., thin belt and broad disk) depending on the spectral type of the central stars (see also Tables 5 and 6, column "OS1+UM1+SL1").For this reference case, the SEDs show the "typical" smooth and steady decrease in the far-IR and millimeter wavelength range resulting from the steady decrease of the absorption cross section C abs toward longer wavelengths and the steady decrease of the dust grain temperature with increasing radial distance from the central star.
The slope of the SED slightly steepens when increasing the wavelength (i.e., the value of α 1000 µm 2000 µm is slightly larger than α 450 µm 850 µm ).The absolute amount of this difference between the spectral indices increases toward stars of the late spectral type.This can be explained by the shift of the wavelength of the dust emission maximum toward longer wavelengths with decreasing effective temperature of the central star.Consequently, the spectral index α 450 µm 850 µm decreases, leading to an increase of the absolute value of ∆ α.

Broken power law grain size distribution
Next, we analyzed the spectral index of the SED for models with a discontinuous grain size distribution described by a broken power law.Our grain size distribution model is defined by three distinct groups of grains: an overabundance of small grains, a depletion of medium-sized grains, and a standard abundance of larger grains (Sect.2.2.2).We assessed the individual impact of each group on the resulting SEDs separately.
Overabundance of small grains.Figure 3 shows the modeled SED with a broken power law grain size distribution depending on the overabundance of small grains, disk structures (thin belts and broad disks), and varying stellar sources.The corresponding spectral indices α 450 µm 850 µm , α 1000 µm 2000 µm , and their differences ∆ α are listed in Table 5.
We found that an overabundance of small grains leads to a modest impact on the millimeter slope of the SED, causing it to become slightly steeper (i.e., slightly higher value of α 1000 µm 2000 µm ), but the overabundance predominantly affects and steepens the far-IR slope of the SED (i.e., higher value of α 450 µm 850 µm ).
A148, page 5 of 14  Table 5.Comparison of the spectral index of the considered discontinuous grain size distributions depending on the overabundance of small grains, stellar parameters (e.g., stellar temperature and radius), and disk parameters (e.g., fractional width ∆R/R).6 for details.

SpT
A148, page 7 of 14 Fig. 5. Modeled SED of debris disks with a narrow belt (∆R/R ∼ 0.17) around an A-type star depending on the overabundance of small grains (colored dotted lines).Total emissions (i.e., net SED including the emission from the star and dust) are indicated with solid lines of the same color for the overabundant of small grains (dashed line).See Table 5 for details.
This behavior becomes more pronounced as the overabundance of small grains increases.In particular, for models with an overabundance of small grains of a factor of 50 (OS50) -compared to the ideal collisional cascade model -we found a significant kink in the slope of the SED in the 800 to 900 µm wavelength range (see the bottom-left plot of Fig. 3).Furthermore, we found that the spectral index transition is more pronounced in debris disks around brighter (hotter) central stars and larger disk radial extent.The systems with broader debris rings (i.e., disks that have an inner radius significantly close to the central star and an outer radius located farther out, compared to their thin-ring counterpart) exhibited a larger range of temperatures, with higher dust temperatures closer to the star resulting in stronger emissions and the location of their maximum being shifted toward shorter wavelengths.This leads to a steeper slope of the SED in the 450-850 µm wavelength range (e.g., larger absolute values of the index α 450 µm 850 µm ).However, compared to the impact of different abundances in the different grain size distributions, the impact of the fractional width ∆ R/R on the spectral index is small.
Next, we investigated the contribution of grains in different size regimes to the observed spectral index.For this purpose, we examined their individual contributions to the total flux and thus their impact on the spectral indices (see Fig. 5).Dust particles with sizes comparable to the observing wavelength emit most efficiently in this wavelength regime (Backman & Paresce 1993).We found that small grains, which constitute a significant portion of the emitting surface area in the grain size distribution, primarily affect the spectral index α 450 µm 850 µm , as they efficiently reemit in the far-IR, while their contribution to the long wavelength regime, and thus their impact on α 1000 µm 2000 µm , is negligible.On the other hand, grains larger than the medium size (e.g., SL1), whose abundance is held constant in this model, have a noticeable impact on the SED and the spectral index in the millimeter range.In particular, grains larger than ∼1 mm, which radiate in a manner similar to black bodies at longer wavelengths, cause a shallower slope of the millimeter SED.

Underabundance of medium grains.
Figure 4 shows the SED for the case of the broken power law grain size distribution Fig. 6.Modeled SED of debris disks with a narrow belt (∆R/R ∼ 0.17) around an A-type star depending on the underabundance of medium grains (color dot lines).Total emissions (i.e., net SED including the emission from the star and dust) are indicated with solid lines of the same color for the underabundance of medium grains.See Table 5 for details.
depending on the underabundance of medium grains, disk structures (thin belts and broad disks), and the spectral type of the central stars.The corresponding spectral indices α 450 µm 850 µm and α 1000 µm 2000 µm and their differences ∆ α are listed in Table 6.As an exemplary case, the contribution of medium grains for different levels of their underabundance to the SED is shown in Fig. 6.We found that an underabundance of medium grains leads to a slight steepening of the far-IR slope of the SED (i.e., the slightly higher value of α 450 µm 850 µm ).However, it primarily affects the millimeter slope of the SED, which becomes shallower, corresponding to a lower spectral index α 1000 µm 2000 µm .Therefore, the presence of both an overabundance of small grains and an underabundance of medium grains increases the absolute difference between α 450 µm 850 µm and α 1000 µm 2000 µm , leaving the 800-900 µm kink as a characteristic impact on the reemission SED.
However, we did not observe a significant spectral inversion solely due to the underabundance of medium grains.Thus, we conclude that the impact of the underabundance of medium grains on the spectral transition is weaker compared to the effect of the overabundance of small grains.Additionally, we found that the spectral index transition is more pronounced in debris disks around brighter (hotter) central stars and broader disk systems, which is consistent with the findings regarding the overabundance of small grains.
We note that the impact of the underabundance of medium grains could lead to a higher value of α 450 µm 850 µm if the disk is massive enough (e.g., more mass in the medium-to-large-sized particles).However, this model could also lead to a steeper slope of the SED at a longer wavelength (i.e., a higher value of α 1000 µm 2000 µm ).Thus, the grain size model with the underabundance of medium grains hardly shows a noticeable spectral inversion.2.0 ± 0.18 to 2.0 ± 0.32.Furthermore, MacGregor et al. (2016) found that millimeter spectral indices α mm (i.e., the flux density for the complete sample at a given millimeter wavelength combined with the VLA 9 millimeter or ATCA 7 millimeter fluxes; Ricci et al. (2012Ricci et al. ( , 2015a,b) ,b) have the value of 1.63 to 3.08, with mean values ∼2.55, and inferred for the slope of the size distribution a range from −2.84 (e.g., HD 141569) to −3.64 (e.g., HD 104860), with the weighted mean of these values as −3.36 ± 0.02.These results indicate a shallower slope of the size distribution compared to standard collisional models.Matthews et al. (2015) also indicated a spectral index of approximately two as being desired at sub-millimeter and millimeter wavelengths for AU Mic, which is lower than expected.These findings are in contrast to the steeper spectral indices α 450 µm 850 µm measured at shorter wavelengths, showing 2.8 to 2.9 (Holland et al. 2017).In general, the SONS JCMT/SCUBA2 Legacy Survey (Holland et al. 2017) found that the spectral indices α 450 µm 850 µm in 48 debris disks range from 2.7 to 4.7, indicating some size distributions steeper than one based on the collisional cascade.Further spectral index inversions were observed in the debris disks around q1 Eri (i.e., from 2.90 for α 350 µm 860 µm to 2.39 for α 1260 µm 6800 µm ; Holland et al. 2017, Ricci et al. 2015b), HD 98 800 (i.e., from 3.21 for α 350 µm 870 µm to 1.00 for α 1350 µm 2000 µm ; Nilsson et al. 2010, Sylvester et al. 2001), andHD 95086 (i.e., from 3.18 for α 500 µm 1300 µm to 2.34 for α 1300 µm 6800 µm ; Su et al. 2017, Ricci et al. 2015b).

Possible degeneracies with the grain composition
Variations in dust grain compositions, such as the inclusion of ice or consideration of porosity, affect the absorption and scattering efficiencies (Kim et al. 2019), thereby influencing the overall shape of the SED (Lebreton et al. 2012;Schüppler et al. 2015).To understand the effect of dust compositions on the spectral inversion and potential degeneracy, we examined the implications of altering chemical compositions, specifically by integrating crystalline water ice into icy astrosilicate mixtures (e.g., 50% ice volume fraction F ice = 0.5), and the microstructure of the dust grains when introducing moderate porosity (e.g., 50% porosity P = 0.5) through effective medium approximations (e.g., Maxwell-Garnett rule; Maxwell Garnett 1904) based on the grain size distribution described by the single power law and broken power law.
We found that a single power law grain size distribution and a broken power law grain size distribution with lower over and underabundance of small and medium grains (e.g., OS10+UM1+SL1 and OS1+UM10+SL1; see the first and second columns of Table 7, and Figs.7 and 8) result in a decrease in the difference in their slopes (i.e., smaller ∆α) and might help the situation.However, these modifications not only steepen the slope in the far-IR range but also affect the millimeter range, resulting in relatively small spectral inversions.Furthermore, both models lead to even smaller deviations in the spectral slopes, particularly with an increase in over and underabundance of small and medium grains based on the broken power law grain size distribution model (e.g., OS50+UM1+SL1 and OS1+UM50+SL1; see the third and fourth columns of Table 7, and Figs.7 and  8).Consequently, these modifications to the dust compositions alone are insufficient to reproduce the observed spectral inversion.This indicates that the specific composition is not the sole factor responsible for the spectral inversions seen in the far-IR to millimeter ranges of debris disk SEDs.
Overall, the effect of ice inclusion seems stronger than that of porosity in cases of an increase in the over/underabundance of small/medium grains.On the other hand, for the case of underabundance of medium grains, all of our models with the modification of dust composition fail to reproduce spectral inversion, highlighting the more significant effect of a higher abundance of small grains.Moreover, the fractional width does not seem to play a significant role in these outcomes, suggesting that the spectral characteristics are largely independent of this A148, page 9 of 14   7. SEDs of the case of the overabundance of small grains depending on the dust composition (ice inclusion F ice = 0.5 and porosity P = 0.5).Solid and dashed lines indicate the modeled SED of the dust reemission for the two different disk structures (thin belts and broad disks).The dotted lines represent the direct stellar radiation from the A star.The top-left plot shows the result based on the single power law grain size distribution (OS1+UM1+SL1), while the other three plots exhibit the results based on the broken power law grain size distribution, specifically the overabundance of small grains with the different scaling factors, 10, 50, and 100 (i.e., OS10, OS50, and OS100).See Sect.3.2 and Table 7 for details.aspect of grain size distribution.This finding is consistent with our previous results in Sect.3.

Possible connection between the inferred size distribution and the collisional state of the disk
As MacGregor et al. ( 2016) suggested, the slope of the grain size distribution γ can be derived from the measurement of the mm spectral index α mm based on the dust opacity power law index of small particles (e.g., compared to observing wavelengths) β s and an inferred power-law exponent from the dust temperature for each disk (i.e., spectral index of the Planck function α PL ) given wavelength ranges.Since most astronomical dust compositions have very similar values of β s (e.g., 1.5-2) and the dust temperature only affects the Planck function weakly at long wavelengths (e.g., α PL ∼ 1.9-1.95with the value of 2 as the Rayleigh-Jeans limit), α mm directly corresponds to the slope of the grain size distribution γ.
The inferred range of the parameters γ based on α mm for the overabundance of small grains from our models (e.g., ∼ − 3.5 to −3.6 regardless of stellar properties and dust compositions) is in good agreement with the predictions of steady-state collisional models, both from the classical model and modern models (e.g., Dohnanyi 1969;Marshall et al. 2017), while the one from underabundance of medium grains showing a bit shallower slopes of grain size distributions (e.g., γ ∼ −3.0), which can be explained by the size dependence of the forces removing the smaller dust grains from the system (e.g., PR drag and stellar radiation).
These observational results yielded valuable insights, but they also highlight the intrinsic challenges faced in such measurements, primarily due to the typical uncertainties (e.g., around 5% to 15%) present in observational data, which lead to spectral inversions that differ by approximately 10%.These differences may stem from the varying sensitivities of telescopes as well as from the diverse calibration and reduction techniques employed to mitigate atmospheric conditions, optical flaws, and environmental sensitivity, all of which impact the accuracy and precision of observations.Additionally, our studies show variability in spectral indices, which appears to depend on the specific properties of the considered debris disks (see Sect. 3).The complexity of these factors underscores advanced solutions, which are anticipated to be addressed by the ngVLA (Next Generation Very Large Array;Di Francesco et al. 2019).

Effect of the blow-out size dependent grain size distribution on the SEDs
We also considered the grain size distribution with the overabundance and/or underabundance of different grain sizes as a function of stellar properties (i.e., blowout size; see also A148, page 10 of 14   8. SEDs of the case of the overabundance of small grains depending on the dust composition (ice inclusion F ice = 0.5 and porosity P = 0.5).Solid and dashed lines indicate the modeled SED of the dust reemission for the two different disk structures (thin belts and broad disks).The dotted lines represent the direct stellar radiation from the A star.The top-left plot shows the result based on the single power law grain size distribution (OS1+UM1+SL1), while the other three plots exhibit the results based on the broken power law grain size distribution, specifically the underabundance of medium grains with the different scaling factors, 10, 50, and 100 (i.e., UM10, UM50, and UM100).See Sect.3.2 and Table 7 for details.
Sect. 2.2.2).We found that the variations in the blowout size and spectral type that are responsible for the location of the first and second maximum/minimum in the grain size distribution (see also Fig. 2 and Table 4), in turn, change the spectral slope of the SED at varying wavelength regimes (Figs. 9 and 10).For example, a grain size distribution with an overabundance of small grains in the case of an A-type star displays a change in the slope at significantly longer wavelengths (e.g., of approximately a few millimeters).On the other hand, K-and M-type stars show this impact at even shorter wavelengths (e.g., at about a few hundred micrometers; see Fig. 9).In the case of an F-type star, the results are comparable to those previously discussed in Sect.2.2.2 because of a similar grain size distribution.These results illustrate that the location of the gap in the grain size distribution has a direct impact on the wavelength at which the slope of the SED is affected.This finding highlights the impact of the blowout size on the resulting (even sub-millimeter/millimeter) slope of the SED.Furthermore, we found that the underabundance of medium-sized grains has a rather weak impact on the changing slope of the SED (see Figs. 9 and 10; Tables 8; and 9), which is consistent with our previous findings discussed in Sect. 3.2.The results discussed in this section suggest the need to constrain the slope of the SED in the different wavelength regimes in order to disentangle the impact of the spectral type from the impact of other potential physical mechanisms on the grain size distribution and thus the slope of the SED.However, we note that previous observational studies (e.g., Lestrade et al. 2020;MacGregor et al. 2016; see also Sect.4.1) have not identified a clear dependence of the impact of stellar spectral type on the slope of the SED, showing consistent observed spectral breaks occurring at similar wavelength regions, irrespective of the central star types.For a comprehensive understanding of the grain size distribution intricacies, further investigation into the collisional evolution of debris disks based on spectral types is suggested.

Summary
Motivated by the expected deviation of the grain size distribution from the often assumed continuous distribution described by a single power law, we investigated the impact of the over-and underabundance of grains in specific size regimes on the SED of debris disks.Based on a simple exemplary model, we quantified and discussed the impact of different values of over-and underabundances of submicrometer-to millimeter-sized grains in a debris dust ring with two different fractional widths ∆ R/R around central stars with spectral types A, F, G, K, and M on A148, page 11 of 14 OS 100 R/R 1.25 R/R 0.17 Fig. 9. SEDs of the case of the overabundance of small grains depending on the stellar properties (i.e., blowout sizes).Solid and dashed lines indicate the modeled SED of the dust reemission for the two different disk structures (thin belts and broad disks).The dotted lines represent the direct stellar radiations.The top-left plot shows the result based on the single power law grain size distribution (OS1+UM1+SL1), while the other three plots exhibit the results based on the broken power law grain size distribution, specifically the overabundance of small grains with the different scaling factors, 10, 50, and 100 (i.e., OS10, OS50, and OS100).See Table 4 and Fig. 2 for the proposed grain size distributions depending on the central star and blowout sizes (Thébault & Augereau 2007;Kim et al. 2018).
Table 8.Comparison of the spectral index of the proposed discontinuous grain size distributions (see Table 4 and Fig. 2) depending on the overabundance of small grains, stellar parameters (e.g., stellar temperature and radius), and disk parameters (e.g., fractional width ∆R/R).
SpT    4 and Fig. 2 for the proposed grain size distributions depending on the central star and blowout sizes (Thébault & Augereau 2007;Kim et al. 2018).
Table 9.Comparison of the spectral index of the proposed discontinuous grain size distributions (see Table 4 and Fig. 2) depending on the underabundance of medium grains, stellar parameters (e.g., stellar temperature and radius), and disk parameters (e.g., fractional width ∆R/R).
SpT Notes.For the locations of the first and second maximum/minimum in grain size distribution depending on the central star for the simulations, see Table 4.
A148, page 13 of 14 Kim, M., and Wolf, S.: A&A, 683, A148 (2024) the slope of the far-IR to millimeter SED of these systems.We described this shape and its the spectral indices α 450 µm 850 µm and α 1000 µm 2000 µm , and we discussed their differences ∆ α.Our key findings are the following: -The deviation of a grain size distribution from that described by a single power law can be potentially derived from the analysis of the SED in the far-IR to millimeter wavelength range.More specifically, the characteristics of the slope of the SED via spectral indices and wavelengths at which photometric data can be obtained even from ground-based telescopes appear promising.In the considered model characterized by a chosen grain size distribution, that is, the specific size range of grains for which an over-or underabundance was assumed, we found a characteristic kink in the SED around 800-900 µm.-As the change of the slope and the wavelength range of its occurrence is directly linked to the deviation of the grain size distribution from a single power law distribution, these observations potentially provide the key for constraining the underlying physical processes leading to the specific grain size distribution.In particular, these observations provide the potential to constrain the collisional physics (e.g., material strength) of the otherwise hardly accessible properties of the larger grains (≫cm).-The presence of both the overabundance of small grains and the underabundance of medium grains contributes to a steeper slope in the SED at the far-IR range.However, their effects differ at millimeter wavelengths.Specifically, the overabundance of small grains primarily results in a much higher value α 450 µm 850 µm (i.e., steeper slope), while the α 1000 µm 2000 µm becomes only slightly higher (i.e., a slightly steeper slope).
On the other hand, the underabundance of medium grains leads to a higher value of α 450 µm 850 µm (i.e., slightly steeper slope), but it primarily results in a lower value of α 1000 µm 2000 µm (i.e., shallower slope).The impact of the overabundance of small dust particles is stronger than that of the underabundance of medium-sized dust particles.
-We find that the spectral index transition is more pronounced for debris disks around brighter (i.e., early-type) central stars and broader disks.However, the disk structure, such as the fractional width of the disk system, has only a small effect on the far-IR to millimeter slope of the SED if compared to the impact of the grain size distributions.-We find that dust composition is not the sole physical mechanism responsible for the spectral inversion observed in the far-IR to (sub-)millimeter part of the SEDs of debris disk systems.-If we assume that the location of the first and second maximum and minimum of the grain size distribution depends on the blowout size and thus on the stellar spectral type, the corresponding impact on the sub-millimeter and millimeter SED is wavelength dependent.Thus, further investigations into the collisional evolution of debris disks depending on spectral types are suggested for a more comprehensive understanding of the intricacies of grain size distribution.-The inherent challenges arising from typical uncertainties in observational data, which result in varying spectral slopes, underscore the need for advanced future observations, such as those provided by ngVLA.

Fig. 3 .
Fig.3.SEDs of the case of the overabundance of small grains.Solid and dashed lines indicate the modeled SED of the dust reemission for the two different disk structures (thin belts and broad disks).The dotted lines represent the direct stellar radiation.The top-left plot shows the result based on the single power law grain size distribution (OS1+UM1+SL1), while the other three plots exhibit the results based on the broken power law grain size distribution, specifically the overabundance of small grains with the different scaling factors, 10, 50, and 100 (i.e., OS10, OS50, and OS100).See Sect.3.2 and Table5for details.

Fig. 4 .
Fig. 4. SEDs of the case of the underabundance of medium grains.Solid and dashed lines indicate the modeled SED of the dust reemission for the two different disk structures (thin belts and broad disks).The dotted lines represent the direct stellar radiation.The top-left plot shows the result based on the single power law grain size distribution (OS1+UM1+SL1), while the other three plots exhibit the results based on the broken power law grain size distribution, specifically the underabundance of medium grains with the different scaling factors, 10, 50, and 100 (i.e., UM10, UM50, and UM100).See Sect.3.2 and Table6for details.

Fig.
Fig. 7. SEDs of the case of the overabundance of small grains depending on the dust composition (ice inclusion F ice = 0.5 and porosity P = 0.5).Solid and dashed lines indicate the modeled SED of the dust reemission for the two different disk structures (thin belts and broad disks).The dotted lines represent the direct stellar radiation from the A star.The top-left plot shows the result based on the single power law grain size distribution (OS1+UM1+SL1), while the other three plots exhibit the results based on the broken power law grain size distribution, specifically the overabundance of small grains with the different scaling factors, 10, 50, and 100 (i.e., OS10, OS50, and OS100).See Sect.3.2 and Table7for details.

Fig.
Fig.8.SEDs of the case of the overabundance of small grains depending on the dust composition (ice inclusion F ice = 0.5 and porosity P = 0.5).Solid and dashed lines indicate the modeled SED of the dust reemission for the two different disk structures (thin belts and broad disks).The dotted lines represent the direct stellar radiation from the A star.The top-left plot shows the result based on the single power law grain size distribution (OS1+UM1+SL1), while the other three plots exhibit the results based on the broken power law grain size distribution, specifically the underabundance of medium grains with the different scaling factors, 10, 50, and 100 (i.e., UM10, UM50, and UM100).See Sect.3.2 and Table7for details.

Fig. 10 .
Fig. 10.SEDs of the case of the underabundance of medium grains depending on the stellar properties (i.e., blowout sizes).Solid and dashed lines indicate the modeled SED of the dust reemission for the two different disk structures (thin belts and broad disks).The dotted lines represent the direct stellar radiations.The top-left plot shows the result based on the single power law grain size distribution (OS1+UM1+SL1), while the other three plots exhibit the results based on the broken power law grain size distribution, specifically the underabundance of medium grains with the different scaling factors, 10, 50, and 100 (i.e., UM10, UM50, and UM100).See Table4and Fig.2for the proposed grain size distributions depending on the central star and blowout sizes(Thébault & Augereau 2007;Kim et al. 2018).

Table 3 .
Considered grain size distributions, dust mass, and corresponding model ID.

Table 4 .
Locations of the first and second maximum/minimum in grain size distribution depending on the central star for the simulations.

Table 6 .
Comparison of the spectral index of the considered discontinuous grain size distributions depending on the underabundance of medium grains, stellar parameters (e.g., stellar temperature and radius), and disk parameters (e.g., fractional width ∆R/R).

Table 7 .
Comparison of the spectral index of the considered discontinuous grain size distributions depending on the overabundance of small grains, disk parameters (e.g., fractional width ∆R/R), and ice volume fraction F ice .
Notes.For the definition of model IDs (OS, UM and SL), see Table3.