Free Access
Issue
A&A
Volume 618, October 2018
Article Number A38
Number of page(s) 14
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/201833061
Published online 11 October 2018

© ESO 2018

1 Introduction

Debris disks around main-sequence stars contain fine dust grains, bur the lifetime of these dust grains is much shorter than the stellar age (Artymowicz & Clampin 1997). Thus, these dust grains are not primordial and should be replenished continuously (Backman & Paresce 1993; Artymowicz & Clampin 1997). Consequently, replenishment occurs by collisions, which are at least partially destructive, in other words, planetesimals in debris disks must have sufficiently high relative velocities to create fragments with a broad size distribution, extending down to μm-sized dust grains.

Spatially resolved images of many debris disks show azimuthal asymmetries of surface brightness (e.g., Kalas et al. 2005; Lagrange et al. 2016; Olofsson et al. 2016; Pan et al. 2016; MacGregor et al. 2017). In (sub)millimeter wavelength observations, which trace large grains and so the dust parent bodies, such disks often appear as narrow eccentric belts with some offset between the belt center and the star. The most likely explanation for the eccentricities and offsets is secular perturbations by as yet undiscovered planets (e.g., Wyatt et al. 1999; Faramaz et al. 2014). Therefore, interpretation of debris disk observations may allow one to constrain parameters of the alleged planets.

Another goal of the analysis of debris disk observations is to reconstruct the hidden planetesimal populations that produce the observed dust by collisions. Unlike planetesimals, dust is subject to stellar radiation forces that make its distribution different from that of the parent planetesimals. To summarize, in order to constrain parameters of the putative planets and dust parent planetesimals, one has to understand combined effects of the dust-producing collisions in the dynamically perturbed belts and stellar radiation forces exerted on dust.

The link between the dust distributions and observable quantities is also important. What is actually observed is radiation scattered or re-emitted by dust, rather than dust itself. Observations of debris disks at a given wavelength are primarily sensitive to dust particles of a comparable size and to certain dust locations – those at which dust emission at this particular wavelength is the strongest. As a result, the observational appearance of debris disks does not directly reflect the underlying dust distributions, and the same disk may reveal dissimilar structure and shape when viewed at different wavelengths. All this further complicates the analysis.

In recent years, many studies have been focused on collisions and radiation pressure effects in debris disks (e.g., Dominik & Decin 2003; Thébault et al. 2003, 2014; Wyatt 2005; Kenyon & Bromley 2005; Thébault & Augereau 2007; Stark & Kuchner 2008, 2009; Kuchner & Stark 2010; Thébault 2012; Nesvold et al. 2013; Kral et al. 2013, 2015; Nesvold & Kuchner 2015a,b; Esposito et al. 2016). For these studies, N-body codes have been used that follow the trajectories of individual particles by numerically integrating their equations of motion. In our study, we use the kinetic approach of statistical physics, which is based on the continuity-like equation for the distribution of dust in an appropriate phase space, and model the evolution of the phase space distributions of the material (Krivov et al. 2006).

The goal of the current study is to answer following questions: What are the impact of selected collisional parameters on the wavelength-dependent observational appearance of debris disks? Is it possible to constrain selected collisional parameters based on specific observational quantities of debris disks? Based on the work by Löhne et al. (2017), we have conducted studies on planetesimal collisions in debris disks and investigate the dependence of the appearance of debris disks on essential collisional parameters, such as the eccentricity of the parent belt, dynamical excitation of the system (i.e., the dispersion of the eccentricities of parent belt bodies), and the critical specific energy for fragmentation (i.e., material strength) of dust particles.

This paper is organized as follows. Section 2 describes the underlying physical process of collisions in the debris disks.In Sect. 3, we briefly introduce the applied tools: the collisional code ACE and the code Debris disks around Main-sequence Stars (DMS) for the simulation of observations. In Sect. 4, we quantify and discuss impacts of the collisional evolution on observables, in particular, their spectral energy distribution (SED) and brightness asymmetries in spatially resolved images. We summarize our findings in Sect. 5.

2 Collisional physics in debris disks

The probability of two spherical objects on Keplerian orbits to collide can be expressed as the product of the overlap of the two orbits, the number density of grains at the desired location in the space of orbital elements, the relative velocity and the collisional cross section (Krivov et al. 2006). In addition, the collisional outcome depends on many other factors, such as the impact angle, dust species, porosities of dust grains, and hardnesses of projectile and target (Blum & Wurm 2008). Furthermore, debris disks with larger typical orbital eccentricities can significantly increase the geometrical probability of collision andthe collisional velocities (Krivov et al. 2006; Queck et al. 2007; Wyatt 2009), resulting in an increased efficiency of the collisional cascade.

2.1 Influence of eccentric belts on collisions

Various previous studies discussed the brightness asymmetry in eccentric debris disks, for example, ϵ Eridani (Greaves et al. 1998, 2005; Hatzes et al. 2000; Poulton et al. 2006; Pan et al. 2016), HR4796 (Wyatt et al. 1999; Schneider et al. 2009; Thalmann et al. 2011; Lagrange et al. 2012), Fomalhaut (Stapelfeldt et al. 2004; Kalas et al. 2005; Pan et al. 2016; MacGregor et al. 2017), HD 181327 (Schneider et al. 2006), HD 15115 (Kalas et al. 2007), HD 202628 (Krist et al. 2012), HD 115600 (Currie et al. 2015), HD 106906 (Kalas et al. 2015), and HD 61005 (Olofsson et al. 2016), as a result of two competing effects: higher temperatures at the pericenter and higher densities at the apocenter. If the increased emission at short wavelengths at the pericenter exceeds the increased flux due to the over-density at the apocenter, the so-called pericenter glow phenomenon is observed. This effect was detected for the first time in Fomalhaut, a 133 AU ring offset by 15 AU implying a forced eccentricity of 0.11 (Kalas et al. 2005).

On the other hand, Pan et al. (2016) calculated the lowest order in the time-averaged linear number density of particles as a function of longitude, f, in a mild eccentric disks (e ≪ 1) ring along the orbit as lsingle(f)1ecosf.\begin{eqnarray*} l_{\textrm{single}}(f)\propto1 - e\cos f.\end{eqnarray*}(1)

Using a different approach, Marsh et al. (2005) found a similar enhancement at the apocenter. In summary, since particles orbit faster at periastron than at apastron, their number density in Keplerian orbits decreases at the pericenter and increases at the apocenter by the same fractional amount.

The ejection condition can be derived from the ratio between the radiation pressure of the stellar radiation-field to the gravitational force, β (≡ FrpFgrav). A general expression for the β threshold, for which a dust particle is ejected from the system, can be expressed as follows (Murray & Dermott 1999): βeject12(1e21+ecosϕ),\begin{eqnarray*} \beta_{\textrm{eject}} \ge \frac{1}{2}\left (\frac{1 - e^2}{1 + e \cos\phi}\right),\vspace*{-1pt}\end{eqnarray*}(2)

where ϕ corresponds to the longitude of the orbit, at which the particle ejected (ϕ = 0 and ϕ = π correspond to periastron and apastron, respectively). Apparently, the β-ratio, and thus the ejection condition and blowout size, depend on the eccentricity of parent belts.

The gravitational impact of large planetesimals can be observed, for instance, as higher dynamical excitation (Kenyon & Bromley 2004). Dynamical excitation, which can be characterized by the width of the typical eccentricity e and inclination i, determines the relative velocities at which orbits cross. Even though the exact level of stirring in debris disks remains unknown, previous studies of debris disk have shown that stirring models are consistent with observations if a dynamical excitation of up to ~0.1 is assumed (e.g., Wyatt et al. 2007, 2012; Thébault & Wu 2008; Löhne et al. 2008, 2012; Kennedy & Wyatt 2010; Müller et al. 2010).

2.2 Influence of material strength on the outcome of collisions

The second key parameter which determines the outcome of collisions in debris disks is the catastrophic disruption threshold QD*$Q^{*}_{\textrm{D}}$ that is the impact energy per unit of target mass, which is defined as the normalized energy at which the mass of the largest fragment has half of the original target mass (Benz & Asphaug 1999). While for smaller objects this quantity is determined by the material strength, the gravitational binding energy dominates in the case of larger objects (≥1km). Thus, the catastrophic disruption threshold QD*$Q^{*}_{\textrm{D}}$ is described by the sum of two power laws (Davis et al. 1985; Holsapple 1994; Paolicchi et al. 1996; Durda & Dermott 1997; Durda et al. 1998; Benz & Asphaug 1999; Kenyon & Bromley 2004; Stewart & Leinhardt 2009), where we follow Löhne et al. (2012) and assume a size dependence described by Benz & Asphaug (1999), together with the reformulation by Stewart & Leinhardt (2009), regarding the reduced and the total masses of the colliders: QD*=[ Qs(sij1m)bs+Qg(sij1km)bg ](vimp3km s1)0.5+3G(mi+mj)5sij. \begin{align*} Q^{*}_{\textrm{D}} &= \left[Q_{\textrm{s}} \left(\frac{s_{ij}}{1\,\rm m}\right)^{-b_{\textrm{s}}} + Q_{\textrm{g}} \left(\frac{s_{ij}}{1{\rm\,km}}\right)^{b_{\textrm{g}}} \right] \left(\frac{v_{\textrm{imp}}}{3\,\textrm{km~s}^{-1}}\right)^{0.5} \nonumber \\ & \quad + \frac{3G(m_{i}+m_{j})}{5s_{ij}}.\vspace*{-1pt}\end{align*}(3)

Here, the subscripts s and g stand for the material strength and the gravity regimes, respectively. The quantity Qs is relevant for describing shock disruption in the strength regime, while Qg is the corresponding specific energy in the gravity regime, both scaled by the impact velocity vimp. We denote the equivalent radii of the spheres, si and sj, with sij (si3+sj3)1/3$\equiv (s^3_{i} + s^3_{j})^{1/3}$ as the combined volume of the colliders. The last term of Eq. (3) approximates the specific energy required to overcome self-gravity, which is important only for radii ≥30 km in our simulations (Löhne et al. 2012). We choose Qs = Qg = 5 × 102 J kg−1, bs = 0.37, and bg = 1.38 (Benz & Asphaug 1999).

Depending on the masses and the impact energy of the colliders, following outcomes are expected: If the energy is sufficient to overcome both the material strength and the combined gravitational potential of the colliders, disruption and dispersal occur. Collisions with specific energies below this threshold are depicted as cratering if the target retains at least half of its original mass. In both cases, a cloud of smaller fragments is produced from the remnants (Löhne et al. 2017). If both colliders stay intact, they are assumed to separate again, unless the impact velocities are below 10 m s−1. Only for very low relative velocities (≤1–10 m s−1) sticking and hence growth may occur (Poppe et al. 2000). However, such low collision velocities are very rare.

3 Numerical models

3.1 Dynamical modeling

As discussed in Sect. 1, most of the previous studies focused on purely gravitational models of planet-disk interactions and forces using the N-body simulations. Following the trajectories of individual objects is useful when the dynamics are complex, whereas any collisional event is assumed to simply eliminate both colliders (Lecavelier des Etangs et al. 1996). Even though this method is hardly applicable in the case of a sufficiently large number of objects to cover a broad range of particle masses, a superparticle method is applicable in the case of shorter collisional timescales or more violent collisional events (e.g., Nesvold et al. 2013; Kral et al. 2013).

The kinetic method of statistical physics is more suitable for particle distributions and long-term collisional evolution. The kinetic method introduces a multidimensional phase-space distribution of dust (e.g., a distribution of grain sizes, coordinates, and velocities, and the supply, loss, and transport of dust grains). We apply the collisional code Analysis of Collisional Evolution (ACE; Löhne et al. 2017) which allows the simulation of the dust production and removal in the planetesimal belt as well as the dynamical evolution of the whole disk. ACE derives the velocity distribution in particle ensembles with the Boltzmann equation and the size distribution using the Smoluchowski equation for describing the evolution of the solids, including disruptive and cratering collisions. The code allows the implementation of a three-dimensional kinetic model with masses, semi-major axes, and eccentricities as phase-space variables. Thus, the resulting equation can be solved for the phase-space distribution as a function of time, from which mass, size, and spatial distributions can be calculated (Krivov et al. 2006, 2008; Krivov 2010; Löhne et al. 2017). The applied version of the ACE code (Löhne et al. 2017)allows us to consider the radial and azimuthal variations of the disk density distribution. Collisional outcomes are calculated with size-dependent scaling laws for fragmentation and dispersal in both the strength and the gravity regime as mentioned in Sect. 2.2. A more detailed description of this simulation tool can be found in previous studies (e.g., Krivov et al. 2005, 2006, 2008; Löhne et al. 2008, 2012, 2017; Krivov 2010).

3.2 Simulated observations

The results of the dynamical calculations provide the basis for the simulation of the observational appearance of the considered disk models. For this purpose, we apply our newly developed code DMS which allows us to calculate the SED, wavelength- and inclination-dependent scattered light images, and polarization maps as well as thermal re-emission maps. It is based on the approach outlined by Ertel et al. (2011, 2012). To spatially resolve the local density and temperature distribution sufficiently well, an adaptive mesh refinement method in Cartesian coordinates is used. The optical properties of the dust grains are computed using the tool miex (Wolf & Voshchinnikov 2004) which is based on Mie theory (Mie 1908).

In the current study, the density distribution is provided by ACE in a tabular form. Thermal re-emission and scattered radiation maps are computed at different observing wavelengths for each dust particle location and dust property. A more detailed discussion of the DMS code can be found in Appendix A.

3.3 Model parameters

The basic debris disk model applied in our study is described in Löhne et al. (2017). In the following, the basic characteristics of that model are summarized in brief. To quantify and identify the specific influence of collisions on the appearance of debris disks, we model a fiducial idealized typical debris disk system around a Fomalhaut-like star that is an A3V main-sequence star (see Table 1). The stellar photospheric emission is described by the corresponding PHOENIX/NextGen grid models (Hauschildt et al. 1999).

The particles are assumed to be composed of generic astrosilicates (Draine 2003) and ice (Li & Greenberg 1998) in equal volume fractions, with a bulk density of 2.35 g cm−3. We use the effective medium theory (Bruggeman 1935) to calculate the optical properties. Since we consider eccentric debris disks in our model, the blowout limit varies between two extreme orbital points that is periastron and apastron (see Eq. (2)). For the given photospheric spectrum and an eccentricity eb = 0.4, the blowout limit amounts to ~3 and 6 μm at periastron and apastron, respectively. Asymmetries are most easily identified where narrow belts are resolved, which correspond to narrow volumes in the space of orbital elements (Löhne et al. 2017). The relative radial width Δ qb /qb at apastron qb = ab (1 + eb) can be estimated from Δqbqb=(Δabab)2+(Δeb1+eb)2,\begin{eqnarray*} \frac{\mathrm{\Delta} q_{\textrm{b}}}{q_{\textrm{b}}} = \sqrt{\left(\frac{\mathrm{\Delta} a_{\textrm{b}}}{a_{\textrm{b}}}\right)^2 + \left(\frac{\mathrm{\Delta} e_{\textrm{b}}}{1 + e_{\textrm{b}}}\right)^2},\end{eqnarray*}(4)

where ab and eb are varied independently. Here, the quantities ab and eb are the semi-major axes of a perturbed belt object and average belt eccentricity, respectively. The quantity Δ eb denotes the dispersion of the eccentricities of parent belt bodies (i.e., dynamical excitation of the system). For a more detailed description of the model setup, we refer to Löhne et al. (2017).

We consider four different eccentricities of the parent belt, eb = 0.0 (circular), 0.1, 0.2, and 0.4, respectively. In addition, three different levels of dynamical excitation, Δ eb = 0.05, 0.1, and 0.2 are considered. The dispersion is defined by the range of eccentricity values. It is parameterized by the maximum and minimum orbital eccentricities that planetesimal belts have at the onset of the collisional cascade. Furthermore, we consider five different values of the material strength Qs, which is the strength component of the catastrophic disruption threshold QD*$Q_{\textrm{D}}^{*}$, ranging from “very high” to “very low” (see Table 1). These are related to our reference value for Qs by Qs = α× reference value of Qs, with α = 25, 5, 1, 1/5, and 1/25, respectively. Thus, this parameter is spanning a range of 54 –625 in our simulations. We note that QD*$Q_{\textrm{D}}^{*}$ can be as small as 10−2 –10−4 J kg−1 for meter-sized bodies without self-gravity (Bukhari Syed et al. 2017; Whizin et al. 2017).

For each run, we evolve the system for 107 yr. Although many debris disks are much older, we choose this value as a steady-state is achieved within this time (Thébault & Augereau 2007). For the simulation of observations, we consider all grains up to a maximum grain size ofaround 1 mm. The contribution of larger particles to the observable quantities is negligible. Surveys at sub-mm wavelengths have shown that for most debris disks mass is typically ranging from ~10−9 to several 10−7 M (e.g., Greaves et al. 2005, and references therein). Thus, we consider a collisionally evolved steady-state disk with 0.7 × 10−8 M.

Table 1

Model parameters for simulation of the dynamical evolution.

4 Results

In the following, we analyze impacts of the eccentricity eb, dynamical excitation Δeb, and the catastrophicdisruption threshold QD*$Q^{*}_{\textrm{D}}$ on the observational appearance of the considered debris disks. For this purpose, we first discuss impacts of these parameters on the resulting particle size distribution (Sect. 4.1). Subsequently, we analyze the corresponding SED (Sect. 4.2) and spatially resolved observations (Sect. 4.3). We focus our studies on the comparison between the periastron and the apastron sides, of the disk (Sect. 4.4).

4.1 Particle size distribution

We first investigate the dust grain size distributions at periastron and apastron of debris disks with different material strengths. As the geometrical cross section of the dust grain is a proxy for the emitted and scattered light, it is shown as a function of the grain size inFigs. 1 and 2.

In Fig. 1, we investigate the dust grain size distributions at the periastron and apastron of debris disks a belt eccentricity of eb = 0.4 (with dynamical excitation Δeb = 0.1) for different material strengths. Since the radiation pressure has a strong influence on the cutoff in the grain size distribution at smallest sizes, the wavy patterns of the grain size distributions start around the smallest bound grains. We note that the blowout size can be increased with particle porosity and increasing stellar temperature (Burns et al. 1979; Kirchschlager & Wolf 2013; Brunngräber et al. 2017). We find different characteristic wavy patterns in the grain size distributions starting with the depletion of grains with radii below the blowout radius abo, which are caused by the lack of β > 0.5 dust grains (Thébault et al. 2003), depending on different collisional evolutions (Figs. 1 and 2). For grains with radii ~3 μm at the apastron side and grains with radii ~6 μm at the periastron side, a strong wavy pattern develops. This depletion leads to an over-density of slightly larger grains (“first peak”) because grains with radii a < abo are depleted and thus can no longer contribute efficiently to the destruction and erosion processes anymore. The overabundance of grains with radii around the first peak (~1.5 abo) in turn induces another depletion of grains with radii around the “second depletion” (~4–40 abo) that is caused by small high-β grains originating inside the disks. Thus, an efficient destruction is responsible for the deep depletion of objects with radii of up to ~4−40 abo. Qualitatively, this first wavy pattern is less pronounced in the size distribution of higher material strengths, where the impact velocities and thus their rate of destructive collisions is significantly lower. A strong depletion for lower material strengths is found for grains with radii ~60 μm, while the depletion in the case of higher material strengths is valid for radii ≤10 μm (Fig. 1). This depletion eventually leads to an over-density of grains with radii around the “second peak”. The overabundance of these grains is shifted from ~100 μm (high material strength) to ~1000 μm (very low material strength).

These successive domino effects (Thébault & Augereau 2007) propagate the local maximum in the grain size distribution toward bigger sizes and leave a unique characteristic wavy size distribution with a pronounced succession of over-densities and depletions (e.g., Campo-Bagatin et al. 1994; Thébault et al. 2003; Krivov et al. 2006; Thébault & Augereau 2007). Finally, the amplitude of the depletion of sub-mm grains reaches its maximum for the lowest material strength. As a result, the wavy pattern of the grain size distribution becomes deeper and broader for low material strengths. These results resemble the finding of Krivov et al. (2006) and Thébault & Augereau (2007).

In Fig. 2, we now investigate the dust grain size distributions at the periastron and apastron of debris disks with different levels of dynamical excitation (Δeb = 0.05, 0.1, and 0.2). Again, we verify that the size distribution converges toward a quasi-steady-state wave pattern, whose characteristics depend on the collisional evolution. However, we note that these size distributions (Fig. 2) do not converge toward the same large grain radii (~0.1 km) for different levels of dynamical excitation Δeb as for varying the material strength (Fig. 1). If planetesimal belts are more dynamically excited, the number density is – on average – decreased. Thus, the ring widths are different and so are the peak densities.

We find that the wavy grain size distribution in highly dynamically excited disks, for instance, in the Δ eb = 0.2 regime, becomes shallower and narrower (see Fig. 2). For the high-excitation case (broad dispersion of the eccentricities of parent belt bodies that is higher Δeb) the depletion of small grains is less pronounced. However, it is propagated to larger grain sizes. Thus, the depletion of sub-mm grains in the case of highly dynamical excited disks is more pronounced which is similar to the case of the low material strength. In contrast, for the low-excitation case (narrow dispersion of the eccentricities of parent belt bodies that is lower Δ eb), the depletion of small grains is even more pronounced. Thus, the depletion of sub-mm grains is less pronounced similar to the case of the high material strength. Smaller size particles are strongly affected by the stellar radiation force, which makes their orbits eccentric and their collisional velocities higher. As a result, the collisional lifetimes of small dust grains only weakly vary with the level of dynamical excitation. Thus, the role of destructive collisions for smaller particles becomes less important (Thébault & Augereau 2007; Thébault & Wu 2008).

We note that this effect results in the second peak in the grain size distribution which is shifted toward larger grains corresponding to a low-excitation case (Thébault & Wu 2008). Thus the second peak is shifted toward larger grain sizes for smaller values of Δ eb (Fig. 2). In contrast, the rate at which large grains are produced by impacts between much larger planetesimal-sized objects strongly depends on the level of dynamical excitation. Thus, a high dynamical excitation leads to a depletion of grains with radii >40 μm and an increased overabundance of smaller grains close to the blowout size (see Fig. 2). Consequently, the final size distribution reflects both effects: a more pronounced depletion of sub-mm grains and a shift of the second peak toward smaller grain sizes with increasing dynamical excitation Δ eb.

The stirring level in the planetesimal belt and material strength of grains determine the degree of “destructiveness” of collisions (Matthews et al. 2014). Thus, these parameters strongly affect the grain size distribution. Decreasing the material strength has a stronger impact on the grain size distribution than an increased dynamical excitation Δ eb (see Figs. 1 and 2).

In addition, at large grain sizes around 0.1 km, a change in the slope (“knee”) of the size distribution is expected (Campo-Bagatin et al. 1994; Durda et al. 1998; O’Brien & Greenberg 2005; Thébault & Augereau 2007). This is due to the transition from the strength-dominated regime, where the resistance of colliders to impacts decreases with increasing grain size, to the gravity-dominated regime, where the resistance quickly increases with increasing grain size (Durda & Dermott 1997). In our simulations, that knee is only pronounced in very destructive circumstances that is for the lowest material strength considered (at a ~ 0.4 km; see Fig. 1). For higher material strengths, the evolution time of 10 Myr was too short to sufficiently converge the size distribution to a steady state.

Since the cratering effect was not included in previous models (e.g., Krivov et al. 2007; Löhne et al. 2008), wavy patterns were present that were much more pronounced compared to the results from this study. Thus, we can deduce that the cratering impact smoothens wavy patterns. Also, a large fraction of the sub-mm grain depletion is due to cratering impacts (see also Thébault et al. 2003; Thébault & Augereau 2007; Müller et al. 2010). Indeed, small grains of a few μm cannot directly break-up objects larger than ~1 mm (around the second peak), even with their increased impact velocities. However, our results show that they can gradually erode much biggergrains efficiently.

thumbnail Fig. 1

Grainsize distribution at the apastron (left panel) and periastron side (right panel) of debris disks with a belt eccentricity of eb = 0.4 for different material strengths.

thumbnail Fig. 2

Grainsize distribution at the apastron (left panel) and periastron (right panel) sides of debris disks with a belt eccentricity of eb = 0.4 for different levels of dynamical excitation (Δeb = 0.05, 0.1, and 0.2).

thumbnail Fig. 3

Spectral energy distribution for systems with different eccentricities (eb = 0.0, 0.1, 0.2, and 0.4) and dynamical excitations (Δeb = 0.05, 0.1, and 0.2). The dashed black line represents the photospheric emission of the central star.

thumbnail Fig. 4

Spectral energy distribution for systems with different material strengths. The dashed black line represents the photospheric emission of the central star.

4.2 Spectral energy distribution (SED)

We investigate the influence of different geometries and collisional factors on the resulting SED. We focus on the infrared to the sub-mm wavelength range because the offset from the stellar photospheric SED is largest in this region.

Dependence on the eccentricity and dynamical excitation

Löhne et al. (2017) found that the eccentricity of parent belts affects their SED only weakly. Similarly, we find only a weak influence of the combination of the geometrical offset and collisions on the overall SED for different eccentricities, assuming the reference material strength (Qs = 5 × 102 J kg−1; see Fig. 3). The most noticeable change of the SED is the slight increase of the scattered flux with increasing eccentricity in the near-IR wavelength range. Around a wavelength of 10 μm, the difference of fluxes between eccentricities eb = 0.0 and eb = 0.4 reaches its maximum because smaller grains can survive at the apastron side in highly eccentric disks (the different blow-out limits on periastron and apastron side; see Eq. (2)). However, the circumstellar dust is hardly detectable through photometric measurements because the stellar photosphere is brighter by a factor of 105 –107. At wavelengths of ~30–50 μm, the thermal re-emission of the dust exceeds the direct stellar radiation, increasing by up to 30% with increasing eccentricity of parent belt.

Figure 3 also depicts the influence of different levels of dynamical excitation Δ eb. The most noticeable change of the SED is the increase of the scattered flux with increasing the level of dynamical excitation in the near to mid-IR wavelength range. However, again it is hardly detectable through photometric measurements in this wavelength range. Increasing the level of dynamical excitation of eccentric parent belts the disk emission increases by up to 56% at 30–50 μm.

These results are easily explained from the dynamical point of view. Higher eccentricities and an increasing level of dynamical excitation increase the mutual collisional velocities. The efficiency of the collisional cascade is increased resulting in a higher dust production rate. Consequently, the disk appears brighter in those regions where the above effect is strongest.

Dependence on the material strength

Figure 4 depicts the influence of the combination of geometrical offset and collisions on the overall SED for different catastrophic disruption threshold QD*$Q_{\textrm{D}}^{*}$ (material strengths Qs = 1.25 × 104, 2.5 × 103, 5 × 102, 1 × 102, and 0.2 × 102 J kg−1, fixed eccentricity eb = 0.4). In comparison to the results shown in Fig. 3, variations of the material strength have a more pronounced effect on the SED than the variation of the eccentricity. The flux at optical to mid-IR wavelengths increases with decreasing material strength. Thus, the more destructive collisional systems show higher fluxes at wavelengths ≤80 μm. In contrast, the flux at sub-mm wavelengths decreases with decreasing material strength. The deeper and broader wavy pattern of the grain size distribution for low material strengths (see Fig. 1) results in the slope of the SED in this range becomes flatter. However, the SEDs in these wavelength ranges have similar spectral indices which makes it hardly possible to constrain the material strength from the analysis of the SED alone.

To analyze this finding, we compare the contributions of the dust within different size ranges to the resulting SED, which is relevant for the optical appearance (Table 2). We find a larger surface fraction of grains with radii <6 μm that is grains which mostly contribute to the SED at optical to mid-IR wavelengths, with decreasing material strength (second column in Table 2). For lower material strengths, the increased production of small grains result in an increase of the SED at these wavelengths. However, we note that surface fraction differences between other collisional regimes are getting smaller because the small grains are more easily affected by stellar radiation. For grains of intermediate size (15 μm < adust < 4200 μm; third column in Table 2) that is grains which mostly contribute to the SED at far-IR to sub-mm wavelengths we find the opposite trend. Less destructive collisional systems, corresponding to the higher material strengths, result in higher fluxes in this wavelengths regime. However, because of the marginal change of the slope of the SED in this range, it is not possible to quantify thematerial strength from pure photometric measurements. The net surface fraction of dust grains with radii >200 μm is sufficiently large to not to be affected by the stellar radiation pressure (fourth column in Table 2). Thus, collisions become more important here.

The second most obvious difference of the SEDs of disks with different material strengths is the change of the location of their maximum. They are not only shifted from 50 μm for very low material strengths toward ~70 μm for very high material strengths (i.e., for the systems with fewer destructive collisions). They also show an increase of the peak flux with decreasing material strength.

Table 2

A total surface area of the dust in three different size bins resulting from the collisional evolution for different material strengths.

4.3 Spatially resolved images

In Sect. 4.2, we discussed the dependence between the SED and the varying collisional parameters, for example, dynamical excitation Δeb and material strength Qs. However, as the SED is a representation of net fluxes integrated over the entire disk, we now investigate spatially resolved images and derived quantities, which potentially provide a wealth of additional information.

4.3.1 Spectral index map

We begin our analysis with the sub-mm spectral index α3mm850μm$\alpha ^{850\,\mu{\textrm{m}}}_{3\,\textrm{mm}}$, based on thefluxes at 850 μm and 3 mm of the observed disks as a function of the material strength for a fixed belt eccentricity eb = 0.4 (see Figs. 5 and 7). In addition, Fig. 6 depicts the difference between the spectral index α3mm850μm$\alpha ^{850\,\mu {\textrm{m}}}_{3\,\textrm{mm}}$ for the very lowand very high material strengths, which indicates the location where large grains (small grains) are more abundant than in the reference case.

We see that the spectral index at larger distances hardly depends on the material strength. Radiation pressure is acting most efficiently on small grains, pushing them to larger distances. Therefore, we find a relatively large fraction of small particles in the outer regions. Consequently, the spectral index is higher at larger distances than around periastron and apastron. Different blowout limits (see Eq. (2)), depending on the different locations under which grains are launched from different sites along the parent belts (e.g., around periastron and apastron), shape the local size distributions and eventually result in the asymmetries seen in the halos (Lee & Chiang 2016; Löhne et al. 2017). In consequence, the halo of the periastron side is much more extended than the halo of the apastron side (see Figs. 5 and 7). However, the halo is typically too faint to allow measuring its (sub)mm spectral index (Löhne et al. 2017).

In contrast, the spectral index in the region around the parent belt strongly depends on the material strength. Figure 7 illustrates the radial cut of the spectral index around apastron (at right ascension Δ α ~ –20′′, declination Δ δ ~ 0′′) and periastron (at Δ α ~ 8′′, Δ δ ~ 0′′). One would expect an increased amount of smaller grains in the case of smaller material strengths. However, Fig. 5 (right plot), Fig. 6 (right plot), and Fig. 7 (red line) show the opposite behavior that is a lower spectral index in the case of a very low material strength. At λobs = 850 μm, the radiation is dominated by grains with a characteristic grain size (Backman & Paresce 1993) acλcritical/2π ≈ 145 μm (β ≈ 0.0146), while grains with ac ≈ 450 μm (β ≈ 0.0047) dominate the emissions at λobs = 3000 μm. Figure 1 shows that grains with radii ~450 μm are overabundant by a factor of 1.5 compared to grains with radii ~145 μm in the case of the very low material strength, eventually causing the observed lower spectral index in the case of the very low material strength.

thumbnail Fig. 5

Spectral index α3mm850μm$\alpha^{850\,\mu {\textrm{m}}}_{3\,\textrm{mm}}$ map for different material strengths (eb = 0.4). We note the low spectral index in the birth ring close to the central star (black square box) in the case of the very low material strength.

thumbnail Fig. 6

Differences between the spectral index α3mm850μm$\alpha ^{850\,\mu {\textrm{m}}}_{3\textrm{mm}}$ maps for the very high material strength (left panel) and the very low material strength (right panel) and the reference case (see Fig. 5 for comparison).

thumbnail Fig. 7

Radial cut of the spectral index α3mm850μm$\alpha ^{850\,\mu {\textrm{m}}}_{3\,\textrm{mm}}$ map for different material strengths (eb = 0.4) along Δδ = 0 (see Fig. 5 for comparison). The left and right sides of this figure indicate apastron and periastron sides of the disks, respectively.

4.3.2 Spatially resolved disk images

In this section, we discuss impacts of collisional parameters, for instance, eccentricities eb, dynamical excitation Δeb, and the material strength Qs on spatially resolved images and their radial surface brightness profiles from near-IR to sub-mm wavelengths.

Dependence on the belt eccentricity and dynamical excitation

The panels in Fig. 8 show simulated observations of spatially resolved disks at different belt eccentricities (eb = 0.2 and 0.4), different levels of dynamical excitation (Δeb = 0.05, 0.1, and 0.2) for the reference material strength Qs at wavelengths λobs = 2.2, 21, 160, 850, and 3000 μm. Both scattered and thermal re-emission radiation are considered. Grains which contribute to the near-IR flux are the smallest within the entire distribution and thus most sensitive to the strong radiation pressure from the central star. Consequently, debris disks are featureless in this wavelength range. In contrast, the larger grains which are dominating the appearance (emission) at sub-mmwavelengths are only weakly affected by radiation pressure allowing the parent ring to be traced in this wavelength range.

In Fig. 9, radial surface brightness profiles are shown for different levels of dynamical excitation (Δ eb = 0.05, 0.1, and 0.2) of belt eccentricities eb = 0.4 at wavelengths λobs = 21, 850, and 3000 μm, assuming the reference material strength. As outlined in Sect. 4.3.1, the abundance of the smaller particles on the periastron side is reduced. Therefore, the surface brightness is dominated by the large grains on the periastron side, whereas on the apastron side the smaller grains are the dominant population. Thus, if there are more destructive collisions, for instance, due to higher eccentricities or higher dynamical excitations, the production rate of small particles is increased on the apastron side (see Fig. 2). This situation is depicted in the simulated surface brightness profile at λobs = 21 μm. In the case of a higher dynamical excitation (Δeb = 0.2) a brighter apastron side is observed than at lower dynamical excitation (Δ eb = 0.05).

The abundance of small dust grains is also subject to further forces and effects, including Poynting–Robertson effect and stellar drag forces (e.g., Wyatt et al. 1999; Löhne et al. 2017). For example, disks with a higher dynamical excitation show a higher production rate. Furthermore, smaller particles are most efficiently affected by Poynting–Robertson drag. Consequently, the inner regions of dynamically excited debris disks appear more empty if observed at larger wavelengths as compared to observations at shorter wavelengths (tracing these smaller particles). In addition, Fig. 9 also shows that the higher dynamical excitation of the planetesimal belt come along with the increased overall brightness of the halo.

Dependence on the material strength

The panels in Fig. 10 show simulated observations of spatially resolved disks as a function of the material strength Qs at wavelengths λobs = 2.2, 21, 160, 850, and 3000 μm (eb = 0.4). Both scattered and thermal re-emission radiation are considered. The different appearances are due to the different spatial and size distribution of the dust resulting from simulated collisional evolutions with different material strengths (Fig. 1).

As the periastron side of disks is closer to the star, one could expect that it is always brighter than the apastron side due to the higher dust temperatures and a higher rate of scattered radiation. However, as the particles travel more slowly and hence spend more time at the apocenter, an eccentric disk should be denser in this region (see Eq. (1)). Furthermore, a very low material strength amplifies the increased rate of dust production. Thus, destructive collisions may reverse the situation of “pericenter glow” to an “apocenter glow”. This can indeed be seen in Figs. 1113, where the brightness of the apastron side increases toward low material strengths. At sub-mm wavelengths (λobs ≥ 0.5 mm), far from the wavelength of the maximum emission, the impact of the dust temperature on the brightness distribution of the disk is reduced. Consequently, the effect of the pericenter glow is reduced. Instead, the impact of the dust number density on the surface brightness becomes more important. A complementary effect has been recently studied by Pan et al. (2016) and MacGregor et al. (2017). They found that azimuthal temperature asymmetries, due to disk offset, could be compensated by azimuthal asymmetries in the dust density.

In contrast, in the case of an increase of the material strength, the grain size distribution is critically changed by collisions (see Fig. 1). Thus, for very high material strengths, the surface brightness differences (asymmetries) are increased. Consequently, the pericenter glow phenomenon is even more pronounced for higher material strengths.

At shorter observing wavelengths λobs, different material strengths Qs seem only to result in a change of the brightness distribution at the apocenter (see Fig. 11), while the brightness distribution at the pericenter remains almost unchanged. However, the comparison at large distances that is the halo, different material strengths result in a significant change of the brightness distribution in the halo on the periastron side. At longer observing wavelengths λobs, the impact of collisions becomes less important. Only the halo regions show differences in the brightness distribution. This behavior is in agreement with the findings outlined in Sect. 4.3.2: higher collisional rate results in an increased overall brightness of the halo.

At shorter observing wavelengths (see, e.g., left panel of Fig. 11) we find again that smaller particles are abundant inthe inner regions of the apastron sides due to the Poynting–Robertson drag force. In contrast, at longer observing wavelengths, inner regions appear void. In conclusion, we find that the higher collisional activity is responsible for symmetrizing their halo, changing the asymmetry of the peak of parent belts and even filling inner regions with a large fraction of smaller grains.

Figure 12 shows the wavelength-dependent radial brightness profiles for different collisional evolutions for a belt eccentricity eb = 0.2 at wavelengths λobs = 21, 850, and 3000 μm. We find that the effects observed for the belt eccentricity eb = 0.4 and lower material strength, for instance, the apocenter glow, are still present and more pronounced in the disk with a belt eccentricity eb = 0.2 and lower material strength.

thumbnail Fig. 8

Simulated surface brightness from near-IR to sub-mm wavelengths (λobs = 2.2, 21, 160, 850, and 3000 μm), for different belt eccentricities (eb = 0.2 and 0.4) and dynamical excitations (Δeb = 0.05, 0.1, and 0.2) for the reference material strength. Both scattered and thermal re-emission are considered.

thumbnail Fig. 9

Radial surface brightness profile as a function of dynamical excitations (Δ eb) at wavelengths λobs = 21, 850, and 3000 μm.

thumbnail Fig. 10

Simulated surface brightness from near-IR to sub-mm wavelengths at λobs = 2.2, 21, 160, 850, and 3000 μm (different material strengths Qs, fixed belt eccentricities eb = 0.4). Both scattered and thermal re-emission are considered.

thumbnail Fig. 11

Radial surface brightness profile as a function of the material strength; belt eccentricity eb = 0.4; λobs = 21, 850, and 3000 μm; Δ δ = 0.

thumbnail Fig. 12

Radial surface brightness profiles as a function of material strength; belt eccentricity eb = 0.2; λobs = 21, 850, and 3000 μm; Δ δ = 0.

thumbnail Fig. 13

Wavelength-dependence of the apocenter-to-pericenter flux ratio as a function of eccentricity. Black line: no collisions (t = 0 yr) with eb = 0.4,; green line: collisional evolution (t = 107 yr) eb = 0.1; red line: collisional evolution (t = 107 yr) eb = 0.4.

4.4 Constraining collisional parameters from observational quantities

In this section, we summarize the potential to constrain the collisional parameters from the observational quantities of debris disks.

Belt eccentricity eb

Figure 13 depicts the apocenter-to-pericenter flux ratio as a function of wavelength for systems before and after the collisional evolution. We find that collisions have the potential to reduce the pericenter glow at shorter wavelengths. Furthermore, azimuthal asymmetries in the surface brightness at sub-mm wavelengths, which result from an increased dust density at the apastron side (Pan et al. 2016; MacGregor et al. 2017), are decreased due to collisions.

For two different belt eccentricities (eb = 0.1 and eb = 0.4) but the same level of dynamical excitation (Δeb = 0.1), the situation is more complicated. The apocenter-to-pericenter flux ratio in the case of eb = 0.1 is nearly constant for all wavelengths. In contrast, the apocenter-to-pericenter flux ratio for eb = 0.4 is varying, which is due to the different size distributions at the apocenter and pericenter. For eb = 0.1, the grain size distribution at periastron and apastron is much shallower and broader than for eb = 0.4. This is not only due to the large number of smaller dust grains that survive at the apastron side of debris disks. It is also because more destructive collisions lead to an over-density of smaller particles and a depletion of larger particles. As a result, the pericenter glow phenomenon is reduced at shorter observing wavelengths and increased at longer observing wavelengths in debris disks with higher belt eccentricities (e.g., eb = 0.4).

Dynamical excitation Δeb

Figures 14 and 15 depict the wavelength-dependent apocenter-to-pericenter flux ratio for different levels of dynamical excitation (Δeb = 0.05, 0.1, and 0.2) of different belt eccentricities (eb = 0.2 and 0.4). For the case of eb = 0.2, we find the apocenter glow phenomenon at shorter wavelengths in the case of the reference and higher dynamical excitation level (Δ eb = 0.1 and 0.2). The main reasons for this behavior are the proximity of the star to the apastron side and a higher rate of dust production due to higher relative velocities compared to the case of lower dynamical excitation level. For the case of eb = 0.4, the minimum of the apocenter-to-pericenter flux ratio is smaller than in the case of the disk with eb = 0.2 disks that is the more eccentric belt shows a stronger pericenter glow. The temperature dependence due to the proximity to the star is dominating this behavior over the higher production rate of smaller particles.

At short wavelengths, we find again a decrease of the pericenter glow phenomenon with increasing Δ eb (see also Fig. 9). At long wavelengths, the basic trend is similar but we find that the spread in observed flux ratios is much narrower. In the less dynamically excited case, for instance, Δ eb = 0.05 for eb = 0.2, higher densities on the apastron sides can weaken the pericenter glow phenomenon. In conclusion, collisional parameters can be best constrained through observations of highly eccentric belts at short wavelengths.

thumbnail Fig. 14

Wavelength-dependence of the apocenter-to-pericenter flux ratio as a function of dynamical excitation (Δ eb = 0.05, 0.1, and 0.2) for eccentricities eb = 0.2 (left panel) and eb = 0.4 (right panel).

thumbnail Fig. 15

Apocenter-to-pericenter flux ratio for different observing wavelengths λobs (21, 850, and3000 μm) as a functionof dynamical excitation Δeb (reference material strength). The number in each box indicates the apocenter-to-pericenter flux ratio.

Material strength Qs

Figures 16 and 17 show the wavelength-dependent apocenter-to-pericenter flux ratio for different material strengths Qs and different belt eccentricities (eb = 0.1, 0.2, and 0.4). First of all, we see that the material strength and eccentricity are degenerate to a large degree. For all eccentric parent belts (eb = 0.1, 0.2, and 0.4), we find that a decrease of the pericenter glow phenomenon with the decreasing material strength. As described in Sect. 4.3.2, the pericenter glow phenomenon is replaced by the apocenter glow phenomenon in the case of the lower material strengths.

Furthermore, the influence of the material strength is strongest at short wavelengths. The mid-IR images show the apocenter glow phenomenon as long as the material strength is below its reference value. Thus, it is clear that a smaller material strength results in the higher asymmetry in surface brightness and density. In contrast, the pericenter glow phenomenon is increased for the higher material strengths with decreasing eccentricity of the parent belts. At longer wavelengths (sub-mm/mm), the dominating grains are not significantly affected by the stellar radiation pressure anymore that is the apocenter-to-pericenter flux ratio becomes insensitive to details of the underlying collisional evolution.

thumbnail Fig. 16

Wavelength-dependence of the apocenter-to-pericenter flux ratio as a function of eccentricity eb and materialstrength Qs.

thumbnail Fig. 17

Wavelength-dependence of the apocenter-to-pericenter flux ratio for selected values of the eccentricity eb and materialstrength Qs. The number in each box indicates the apocenter-to-pericenter flux ratio. The case of the apocenter-to-pericenter flux ratio >1 is discussed as the apocenter glow phenomenon in Sect. 4.4.

5 Summary

We studied impacts of collisional parameters, for instance, belt eccentricity eb, dynamical excitation Δeb, and the material strength Qs, on the observational appearance of eccentric debris disks. We have found surface brightness asymmetries that are caused by the combined effects of collisions and radiation pressurers. Our key results are as follows:

  • 1.

    An increasing belt eccentricity of the debris disk system leads to the survival of a large number of smaller dust grains on the apastron side. Thus the “pericenter glow” phenomenon is reduced with decreasing material strength at shorter observing wavelengths. However, the geometrical parameters describing the debris disks system that is the proximity of the star to the apastron side, must be considered as well.

  • 2.

    An increasing level of dynamical excitation in the eccentric parent belt of the debris disk system leads to a higher production rate of smaller particles and depletion of sub-mm grains. Consequently, the “pericenter glow” is reduced.

  • 3.

    A decreasing material strength results in a higher production rate of smaller particles, which reduces the surface brightness differences between periastron and apastron. For very low material strengths, the pericenter glow phenomenon is reduced and eventually even replaced by the opposite effect, the apocenter glow in the near to mid-IR wavelength range. In contrast, an increasing material strength results in an increase of asymmetries in the surface brightness distribution.

  • 4.

    The wavelength-dependent apocenter-to-pericenter flux ratio at short wavelengths is increased in the case of a higher efficiency of destructive collisions, caused by increasing dynamical excitation of the belt or decreasing material strengths. This effect is less pronounced at sub-mm wavelengths. Within the considered parameter space, the impact of the material strength QD*$Q^{*}_{\textrm{D}}$ on the appearance of the disk is stronger than that of dynamical excitation of the belt.

  • 5.

    The SED alone does not provide unique constraints on the collisional parameters considered. Thus, deriving unique constraints for the impact of considered collisional parameters based on the SED alone is hardly possible.

Our study motivates the need to have a closer look at the impact of the specific type of dust on the presented results. For example, the specific roles of crystalline ice, amorphous ice, astrosilicates, PAH, and graphite deserve further investigations.

Acknowledgements

We would like to thank the anonymous referee for helpful and constructive suggestions and comments that greatly contributed to improving the final version of the paper. This work was supported by the Research Unit FOR 2285 “Debris Disks in Planetary Systems” of the Deutsche Forschungsgemeinschaft. M.K., S.W., and F.K. acknowledge the DFG for financial support under contracts WO 857/15-1. F.K. also acknowledges financial support through European Research Grant SNDUST 694520. T.L. and A.V.K. acknowledge the DFG for financial support under contracts LO 1715/2-1 and KR 2164/13-1.

Appendix A DMS description

The software tool DMS has been developed for the simulation of the SED and wavelength and inclination dependent scattered light images and polarization maps as well as thermal re-emission maps of an optically thin dust distribution. The key equations for calculation of scattered light intensity and polarization maps (as well as the corresponding contribution to the SED) are as follows: Wλsca=Lλ,*Qλsca(a)πa24πr2S11(θ)dθ,(Wλsca)pol=Lλ,*Qλsca(a)πa24πr2S12(θ)dθ.\begin{align*} & W^{\textrm{sca}}_{\rm\lambda} = L_{\rm\lambda, *} Q^{\textrm{sca}}_{\rm\lambda} (a) \frac{\pi a^{2}}{4 \pi r^{2}} S_{11} (\theta) \textrm{d}\theta, \\ & (W^{\textrm{sca}}_{\rm\lambda})_{\textrm{pol}} = L_{\rm\lambda, *} Q^{\textrm{sca}}_{\rm\lambda} (a) \frac{\pi a^{2}}{4 \pi r^{2}} S_{12} (\theta) \textrm{d}\theta. \end{align*}

Here Lλ,*, Wλsca$W^{\textrm{sca}}_{\rm\lambda}$, and (Wλsca)pol$(W ^{\textrm{sca}}_{\rm\lambda})_{\textrm{pol}}$ are the monochromatic luminosity, as well as the power of scattered and linearly polarized linearly scattered radiation of a spherical dust grain, respectively. The quantities S11 and S12 are elements of the Müller matrix, describing the scattering of radiation by the particle. The quantity θ denotes the scattering angle that is the angle between the directions of incident and scattered radiation. The dust particles are assumed to be a spherical (radius a). Thermal re-emission maps and contribution to the SED are calculated based on following equations: Wλabs=Lλ,*Qλabs(a)πa24πr2,Wλre0emi=4πa2Qλabs(a)Bλ(Tg).\begin{align*} & W^{\textrm{abs}} _{\rm\lambda} = L_{{\rm\lambda}, *} Q^{\textrm{abs}}_{\rm\lambda} (a) \frac{\pi a^{2}}{4 \pi r^{2}}, \\[5pt] & W^{\textrm{re}\hbox{-}\textrm{emi}}_{\rm\lambda} = 4 \pi a^{2} Q^{\textrm{abs}}_{\rm\lambda} (a) B_{\rm\lambda} (T_{\textrm{g}}). \end{align*}

Here Wλabs$W^{\textrm{abs}}_{\rm\lambda}$ and Wλre0emi$W^{\textrm{re}\hbox{-}\textrm{emi}}_{\rm\lambda}$ are the monochromatic power absorbed and re-emitted by a spherical dust grain with the absorption efficiency Qλabs(a)$Q^{\textrm{abs}}_{\rm\lambda} (a)$ and temperature Tg at a distance r from the star, respectively.

The optical properties of the dust grains are computed using the software tool miex (Wolf & Voshchinnikov 2004), which is included in DMS.

References

  1. Acke, B., Min,M., Dominik, C., et al. 2012, A&A, 540, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Artymowicz, P., & Clampin, M. 1997, ApJ, 490, 863 [NASA ADS] [CrossRef] [Google Scholar]
  3. Backman, D. E., & Paresce, F. 1993, in Protostars and Planets III, eds. E. H. Levy, & J. I. Lunine (Tucson, AZ: University of Arizona Press), 1253 [Google Scholar]
  4. Benz, W., & Asphaug, E. 1999, Icarus, 142, 5 [NASA ADS] [CrossRef] [Google Scholar]
  5. Blum, J., & Wurm, G. 2008, ARA&A, 46, 21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Bruggeman, D. A. G. 1935, Ann. Phys., 416, 636 [NASA ADS] [CrossRef] [Google Scholar]
  7. Brunngräber,R., Wolf, S., Kirchschlager, F., & Ertel, S. 2017, MNRAS, 464, 4383 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bukhari Syed, M., Blum, J., Wahlberg Jansson, K., & Johansen, A. 2017, ApJ, 834, 145 [NASA ADS] [CrossRef] [Google Scholar]
  9. Burns, J. A., Lamy, P. L., & Soter, S. 1979, Icarus, 40, 1 [NASA ADS] [CrossRef] [Google Scholar]
  10. Campo-Bagatin, A., Cellino, A., Davis, D. R., Farinella, P., & Paolicchi, P. 1994, Planet. Space Sci., 42, 1079 [NASA ADS] [CrossRef] [Google Scholar]
  11. Currie, T., Lisse, C., & Kuchner, M. 2015, ApJ, 807, L7 [NASA ADS] [CrossRef] [Google Scholar]
  12. Davis, M., Efstathiou, G., Frenk, C. S., & White, S. D. M. 1985, ApJ, 292, 371 [NASA ADS] [CrossRef] [Google Scholar]
  13. Dohnanyi, J. S. 1969, J. Geophys. Res., 74, 2531 [NASA ADS] [CrossRef] [Google Scholar]
  14. Dominik, C., & Decin, G. 2003, ApJ, 598, 626 [NASA ADS] [CrossRef] [Google Scholar]
  15. Draine, B. T. 2003, ApJ, 598, 1017 [NASA ADS] [CrossRef] [Google Scholar]
  16. Durda, D. D., & Dermott, S. 1997, Icarus, 130, 140 [NASA ADS] [CrossRef] [Google Scholar]
  17. Durda, D. D., Greenberg, R. R., & Jedicke, R. 1998, Icarus, 135, 431 [NASA ADS] [CrossRef] [Google Scholar]
  18. Ertel, S., Wolf, S., Metchev, S., et al. 2011, A&A, 533, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Ertel, S., Wolf, S., & Rodmann, J. 2012, A&A, 544, A61 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Esposito, T. M., Fitzgerald, M. P., & Graham, J. R. 2016, ApJ, 152, 85 [Google Scholar]
  21. Faramaz, V., Beust, H., & Thébault, P. 2014, A&A, 563, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Greaves, J. S., Holland, W. S., Moriarty-Schieven, G., et al. 1998, ApJ, 506, 133 [NASA ADS] [CrossRef] [Google Scholar]
  23. Greaves, J. S., Holland, W. S., Wyatt, M. C., et al. 2005, ApJ, 619, L187 [NASA ADS] [CrossRef] [Google Scholar]
  24. Hatzes, A. P., Mkrtichian, D. E., & Kanaan, A. 2000, Variable Stars as Essential Astrophysical Tools, ed. C. Ibanoglu, NATO Sci. Ser. C, 544, 405 [NASA ADS] [Google Scholar]
  25. Hauschildt, P. H., Allard, F., & Baron, E. 1999, ApJ, 512, 337 [Google Scholar]
  26. Holsapple, K. A. 1994, Planet. Space Sci., 42, 1067 [NASA ADS] [CrossRef] [Google Scholar]
  27. Kalas, P., Graham, J. R., & Clampin, M. 2005, Nature, 435, 1067 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  28. Kalas, P., Fitzgerald, M. P., & Graham, J. R. 2007, ApJ, 661, L85 [NASA ADS] [CrossRef] [Google Scholar]
  29. Kalas, P., Rajan, A., Wang, J. J., et al. 2015, ApJ, 814, 32 [NASA ADS] [CrossRef] [Google Scholar]
  30. Kennedy, G. M., & Wyatt, M. C. 2010, MNRAS, 405, 1253 [NASA ADS] [Google Scholar]
  31. Kenyon, S. J., & Bromley, B. C. 2004, ApJ, 602, L133 [NASA ADS] [CrossRef] [Google Scholar]
  32. Kenyon, S. J., & Bromley, B. C. 2005, ApJ, 139, 269 [Google Scholar]
  33. Kirchschlager, K., & Wolf, S. 2013, A&A, 552, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Kral, Q., Thébault, P., & Charnoz, S. 2013, A&A, 558, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Kral, Q., Thébault, P., Augereau, J.-C., Boccaletti, A., & Charnoz, S. 2015, A&A, 573, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Krist, J. E., Stapelfeldt, K. R., Bryden, G., & Plavchan, P. 2012, ApJ, 134, 45 [Google Scholar]
  37. Krivov, A. V. 2010, Res. Astron. Astrophys., 10, 383 [NASA ADS] [CrossRef] [Google Scholar]
  38. Krivov, A. V., Sremčevć, M., & Spahn, F. 2005, Icarus, 174, 105 [NASA ADS] [CrossRef] [Google Scholar]
  39. Krivov, A. V., Löhne, T., & Sremčevć, M. 2006, A&A, 455, 509 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Krivov, A. V., Queck, M., Löhne, T., & Sremčevć, M. 2007, A&A, 462, 199 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Krivov, A. V., Müller, S., Löhne, T., & Mutschke, H. 2008, ApJ, 687, 608 [NASA ADS] [CrossRef] [Google Scholar]
  42. Kuchner, M. J., & Stark, C. C. 2010, AJ, 140, 1007 [NASA ADS] [CrossRef] [Google Scholar]
  43. Lagrange, A.-M., Boccaletti, A., Milli, J., et al. 2012, A&A, 546, 38 [Google Scholar]
  44. Lagrange, A.-M., Mathias, P., Absil, O., et al. 2016, A&A, 586, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Lebreton, J., Augereau, J.-C., Thi, W.-F., et al. 2012, A&A, 539, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Lecavelier des Etangs, A., Scholl, H., Roques, F., Sicardy, B., & Vidal-Madjar, A. 1996, A&A, 307, 542 [NASA ADS] [Google Scholar]
  47. Lee, E. J., & Chiang, E. 2016, ApJ, 827, 125 [NASA ADS] [CrossRef] [Google Scholar]
  48. Li, A., & Greenberg, J. M. 1998, A&A, 331, 291 [NASA ADS] [Google Scholar]
  49. Löhne, T., Krivov, A. V., & Rodmann, J. 2008, ApJ, 673, 1123 [NASA ADS] [CrossRef] [Google Scholar]
  50. Löhne, T., Augereau, J.-C., Ertel, S., et al. 2012, ApJ, 537, A110 [NASA ADS] [CrossRef] [Google Scholar]
  51. Löhne, T., Krivov, A. V., Kirchschlager, F., Sende, J. A., & Wolf, S. 2017, A&A, 605, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. MacGregor, M. A., Kalas, P., Wyatt, M. C., et al. 2017, ApJ, 842, 1 [NASA ADS] [CrossRef] [Google Scholar]
  53. Mamajek, E. E. 2012, ApJ, 754, L20 [NASA ADS] [CrossRef] [Google Scholar]
  54. Marsh, K. A., Velusamy, T., Dowell, C. D., Grogan, K., & Beichman, C. A. 2005, ApJ, 620, L47 [NASA ADS] [CrossRef] [Google Scholar]
  55. Matthews, B. C., Krivov, A. V., Wyatt, M. C., Bryden, G., & Eiroa, C. 2014, in Protostars and Planets VI, eds. H. Beuther, R. Klessen, C. Dullemond, & T. Henning (Tucson, AZ: University of Arizona Press), 521 [Google Scholar]
  56. Mie, G. 1908, Ann. Phys., 330, 377 [Google Scholar]
  57. Müller, S.T.L., & Krivov, A. V. 2010, ApJ, 708, 1728 [NASA ADS] [CrossRef] [Google Scholar]
  58. Murray, C. D., & Dermott, S. F. 1999, Solar System Dynamics (Cambridge: Cambridge University Press), 592 [Google Scholar]
  59. Mustill, A. J., & Wyatt, M. C. 2009, MNRAS, 399, 1403 [NASA ADS] [CrossRef] [Google Scholar]
  60. Nesvold, E. R., & Kuchner, M. J. 2015a, ApJ, 815, 61 [NASA ADS] [CrossRef] [Google Scholar]
  61. Nesvold, E. R., & Kuchner, M. J. 2015b, ApJ, 798, 83 [NASA ADS] [CrossRef] [Google Scholar]
  62. Nesvold, E. R., Kuchner, M. J., Rein, H., & Pan, M. 2013, ApJ, 777, 1444 [Google Scholar]
  63. O’Brien, D. P., & Greenberg, R. 2005, Icarus, 178, 179 [NASA ADS] [CrossRef] [Google Scholar]
  64. Olofsson, J., Samland, M., Avenhaus, H., et al. 2016, A&A, 591, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Pan, M., Nesvold, E. R., & Kuchner, M. J. 2016, ApJ, 832, 81 [NASA ADS] [CrossRef] [Google Scholar]
  66. Paolicchi, P., Verlicchi, A., & Cellino, A. 1996, Icarus, 121, 126 [NASA ADS] [CrossRef] [Google Scholar]
  67. Poppe, T., Blum, J., & Henning, T. 2000, ApJ, 533, 454 [NASA ADS] [CrossRef] [Google Scholar]
  68. Poulton, C. J., Greaves, J. S., & Collier Cameron A. 2006, MNRAS, 372, 53 [NASA ADS] [CrossRef] [Google Scholar]
  69. Queck, M., Krivov, A. V., Sremčevć, M., & Thébault, P., 2007, Celest. Mech. Dyn. Astron., 99, 169 [NASA ADS] [CrossRef] [Google Scholar]
  70. Schneider, G., Silverstone, M. D., Hines, D. C., et al. 2006, ApJ, 650, 414 [NASA ADS] [CrossRef] [Google Scholar]
  71. Schneider, G., Weinberger, A. J., Becklin, E. E., Debes, J. H., & Smith, B. A. 2009, AJ, 137, 53 [NASA ADS] [CrossRef] [Google Scholar]
  72. Stapelfeldt, K. R., Holmes, E. K., Chen, C., et al. 2004, ApJS, 154, 458 [NASA ADS] [CrossRef] [Google Scholar]
  73. Stark, C. C., & Kuchner, M. J. 2008, ApJ, 686, 637 [NASA ADS] [CrossRef] [Google Scholar]
  74. Stark, C. C., & Kuchner, M. J. 2009, ApJ, 707, 543 [NASA ADS] [CrossRef] [Google Scholar]
  75. Stewart, S. T., & Leinhardt, Z. M. 2009, ApJ, 691, L133 [NASA ADS] [CrossRef] [Google Scholar]
  76. Strubbe, L. E., & Chiang, E. I. 2006, ApJ, 638, 652 [NASA ADS] [CrossRef] [Google Scholar]
  77. Thalmann, C., Janson, M., Buenzli, E., et al. 2011, ApJ, 732, L34 [NASA ADS] [CrossRef] [Google Scholar]
  78. Thébault, P. 2012, A&A, 537, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  79. Thébault, P.,& Augereau, J.-C. 2007, A&A, 472, 169 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Thébault, P.,& Wu, Y. 2008, A&A, 481, 713 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  81. Thébault, P., Augereau, J.-C., & Beust, H. 2003, A&A, 408, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  82. Thébault, P., Kral, Q., & Augereau, J.-C. 2014, A&A, 561, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  83. Weissman, P. R. 1984, Science, 224, 987 [NASA ADS] [CrossRef] [Google Scholar]
  84. Whizin, A. D., Blum, J., & Colwell, J. E. 2017, ApJ, 836, 94 [NASA ADS] [CrossRef] [Google Scholar]
  85. Wolf, S., & Voshchinnikov, N. V. 2004, Comput. Phys. Commun., 162, 113 [NASA ADS] [CrossRef] [Google Scholar]
  86. Wyatt, M. C. 2005, A&A, 433, 1007 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  87. Wyatt, M. C. 2008, ARA&A, 46, 339 [NASA ADS] [CrossRef] [Google Scholar]
  88. Wyatt, M. C. 2009, Lect. Notes Phys., 758, 37 [Google Scholar]
  89. Wyatt, M. C., Dermott, S. F., & Telesco, C. M. 1999, ApJ, 527, 918 [NASA ADS] [CrossRef] [Google Scholar]
  90. Wyatt, M. C., Smith, R., Greaves, J. S., et al. 2007, ApJ, 658, 569 [NASA ADS] [CrossRef] [Google Scholar]
  91. Wyatt, M. C., Clarke, C. J., & Booth, M. 2011, Celest. Mech. Dyn. Astron., 111, 1 [NASA ADS] [CrossRef] [Google Scholar]
  92. Wyatt, M. C., Kennedy, G., Sibthorpe, B., et al. 2012, MNRAS, 424, 1206 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Model parameters for simulation of the dynamical evolution.

Table 2

A total surface area of the dust in three different size bins resulting from the collisional evolution for different material strengths.

All Figures

thumbnail Fig. 1

Grainsize distribution at the apastron (left panel) and periastron side (right panel) of debris disks with a belt eccentricity of eb = 0.4 for different material strengths.

In the text
thumbnail Fig. 2

Grainsize distribution at the apastron (left panel) and periastron (right panel) sides of debris disks with a belt eccentricity of eb = 0.4 for different levels of dynamical excitation (Δeb = 0.05, 0.1, and 0.2).

In the text
thumbnail Fig. 3

Spectral energy distribution for systems with different eccentricities (eb = 0.0, 0.1, 0.2, and 0.4) and dynamical excitations (Δeb = 0.05, 0.1, and 0.2). The dashed black line represents the photospheric emission of the central star.

In the text
thumbnail Fig. 4

Spectral energy distribution for systems with different material strengths. The dashed black line represents the photospheric emission of the central star.

In the text
thumbnail Fig. 5

Spectral index α3mm850μm$\alpha^{850\,\mu {\textrm{m}}}_{3\,\textrm{mm}}$ map for different material strengths (eb = 0.4). We note the low spectral index in the birth ring close to the central star (black square box) in the case of the very low material strength.

In the text
thumbnail Fig. 6

Differences between the spectral index α3mm850μm$\alpha ^{850\,\mu {\textrm{m}}}_{3\textrm{mm}}$ maps for the very high material strength (left panel) and the very low material strength (right panel) and the reference case (see Fig. 5 for comparison).

In the text
thumbnail Fig. 7

Radial cut of the spectral index α3mm850μm$\alpha ^{850\,\mu {\textrm{m}}}_{3\,\textrm{mm}}$ map for different material strengths (eb = 0.4) along Δδ = 0 (see Fig. 5 for comparison). The left and right sides of this figure indicate apastron and periastron sides of the disks, respectively.

In the text
thumbnail Fig. 8

Simulated surface brightness from near-IR to sub-mm wavelengths (λobs = 2.2, 21, 160, 850, and 3000 μm), for different belt eccentricities (eb = 0.2 and 0.4) and dynamical excitations (Δeb = 0.05, 0.1, and 0.2) for the reference material strength. Both scattered and thermal re-emission are considered.

In the text
thumbnail Fig. 9

Radial surface brightness profile as a function of dynamical excitations (Δ eb) at wavelengths λobs = 21, 850, and 3000 μm.

In the text
thumbnail Fig. 10

Simulated surface brightness from near-IR to sub-mm wavelengths at λobs = 2.2, 21, 160, 850, and 3000 μm (different material strengths Qs, fixed belt eccentricities eb = 0.4). Both scattered and thermal re-emission are considered.

In the text
thumbnail Fig. 11

Radial surface brightness profile as a function of the material strength; belt eccentricity eb = 0.4; λobs = 21, 850, and 3000 μm; Δ δ = 0.

In the text
thumbnail Fig. 12

Radial surface brightness profiles as a function of material strength; belt eccentricity eb = 0.2; λobs = 21, 850, and 3000 μm; Δ δ = 0.

In the text
thumbnail Fig. 13

Wavelength-dependence of the apocenter-to-pericenter flux ratio as a function of eccentricity. Black line: no collisions (t = 0 yr) with eb = 0.4,; green line: collisional evolution (t = 107 yr) eb = 0.1; red line: collisional evolution (t = 107 yr) eb = 0.4.

In the text
thumbnail Fig. 14

Wavelength-dependence of the apocenter-to-pericenter flux ratio as a function of dynamical excitation (Δ eb = 0.05, 0.1, and 0.2) for eccentricities eb = 0.2 (left panel) and eb = 0.4 (right panel).

In the text
thumbnail Fig. 15

Apocenter-to-pericenter flux ratio for different observing wavelengths λobs (21, 850, and3000 μm) as a functionof dynamical excitation Δeb (reference material strength). The number in each box indicates the apocenter-to-pericenter flux ratio.

In the text
thumbnail Fig. 16

Wavelength-dependence of the apocenter-to-pericenter flux ratio as a function of eccentricity eb and materialstrength Qs.

In the text
thumbnail Fig. 17

Wavelength-dependence of the apocenter-to-pericenter flux ratio for selected values of the eccentricity eb and materialstrength Qs. The number in each box indicates the apocenter-to-pericenter flux ratio. The case of the apocenter-to-pericenter flux ratio >1 is discussed as the apocenter glow phenomenon in Sect. 4.4.

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.