No Arabic abstract
The azimuthal substructure observed in some debris disks, as exemplified by epsilon Eridani, is usually attributed to resonances with embedded planets. In a standard scenario, the Poynting-Robertson force, possibly enhanced by the stellar wind drag, is responsible for the delivery of dust from outer regions of the disk to locations of external mean-motion planetary resonances; the captured particles then create characteristic ``clumps. Alternatively, it has been suggested that the observed features in systems like epsilon Eri may stem from populations of planetesimals that have been captured in resonances with the planet, such as Plutinos and Trojans in the solar system. A large fraction of dust produced by these bodies would stay locked in the same resonance, creating the dusty clumps. To investigate both scenarios and their applicability limits for a wide range of stars, planets, disk densities, and planetesimal families we construct simple analytic models for both scenarios. In particular, we show that the first scenario works for disks with the pole-on optical depths below about ~10^{-4}-10^{-5}. Above this optical depth level, the first scenario will generate a narrow resonant ring with a hardly visible azimuthal structure, rather than clumps. The efficiency of the second scenario is proportional to the mass of the resonant planetesimal family, as example, a family with a total mass of ~0.01 to 0.1 Earth masses could be sufficient to account for the clumps of epsilon Eridani.
In this paper, we present results from a multi-stage numerical campaign to begin to explain and determine why extreme debris disk detections are rare, what types of impacts will result in extreme debris disks and what we can learn about the parameters of the collision from the extreme debris disks. We begin by simulating many giant impacts using a smoothed particle hydrodynamical code with tabulated equations of state and track the escaping vapour from the collision. Using an $N$-body code, we simulate the spatial evolution of the vapour generated dust post-impact. We show that impacts release vapour anisotropically not isotropically as has been assumed previously and that the distribution of the resulting generated dust is dependent on the mass ratio and impact angle of the collision. In addition, we show that the anisotropic distribution of post-collision dust can cause the formation or lack of formation of the short-term variation in flux depending on the orientation of the collision with respect to the orbit around the central star. Finally, our results suggest that there is a narrow region of semi-major axis where a vapour generated disk would be observable for any significant amount of time implying that giant impacts where most of the escaping mass is in vapour would not be observed often but this does not mean that the collisions are not occurring.
Observations of debris disks offer a window into the physical and dynamical properties of planetesimals in extrasolar systems through the size distribution of dust grains. In particular, the millimeter spectral index of thermal dust emission encodes information on the grain size distribution. We have made new VLA observations of a sample of seven nearby debris disks at 9 mm, with 3 resolution and $sim5$ $mu$Jy/beam rms. We combine these with archival ATCA observations of eight additional debris disks observed at 7 mm, together with up-to-date observations of all disks at (sub)millimeter wavelengths from the literature to place tight constraints on the millimeter spectral indices and thus grain size distributions. The analysis gives a weighted mean for the slope of the power law grain size distribution, $n(a)propto a^{-q}$, of $langle q rangle = 3.36pm0.02$, with a possible trend of decreasing $q$ for later spectral type stars. We compare our results to a range of theoretical models of collisional cascades, from the standard self-similar, steady-state size distribution ($q=3.5$) to solutions that incorporate more realistic physics such as alternative velocity distributions and material strengths, the possibility of a cutoff at small dust sizes from radiation pressure, as well as results from detailed dynamical calculations of specific disks. Such effects can lead to size distributions consistent with the data, and plausibly the observed scatter in spectral indices. For the AU Mic system, the VLA observations show clear evidence of a highly variable stellar emission component; this stellar activity obviates the need to invoke the presence of an asteroid belt to explain the previously reported compact millimeter source in this system.
We explore the collisional decay of disk mass and infrared emission in debris disks. With models, we show that the rate of the decay varies throughout the evolution of the disks, increasing its rate up to a certain point, which is followed by a leveling off to a slower value. The total disk mass falls off ~ t^-0.35 at its fastest point (where t is time) for our reference model, while the dust mass and its proxy -- the infrared excess emission -- fades significantly faster (~ t^-0.8). These later level off to a decay rate of M_tot(t) ~ t^-0.08 and M_dust(t) or L_ir(t) ~ t^-0.6. This is slower than the ~ t^-1 decay given for all three system parameters by traditional analytic models. We also compile an extensive catalog of Spitzer and Herschel 24, 70, and 100 micron observations. Assuming a log-normal distribution of initial disk masses, we generate model population decay curves for the fraction of debris disk harboring stars observed at 24 micron and also model the distribution of measured excesses at the far-IR wavelengths (70-100 micron) at certain age regimes. We show general agreement at 24 micron between the decay of our numerical collisional population synthesis model and observations up to a Gyr. We associate offsets above a Gyr to stochastic events in a few select systems. We cannot fit the decay in the far infrared convincingly with grain strength properties appropriate for silicates, but those of water ice give fits more consistent with the observations.
The presence of gas in dusty debris disks around main-sequence stars is reviewed. We present new observational results for the most prominent representative of the class, viz. the southern naked-eye star beta Pictoris. The spatial and spectral distribution of observed atomic lines from the disk around the star is reproducable by a Keplerian rotation model to a high degree of accuracy. The expected velocity dispersion due to radiation pressure in resonance lines is not observed. Modeling the motion of different atomic species under the influence of gravity, radiation pressure and gas friction leads to the conclusion that an underlying decelerating component must be present in the disk. This braking agent is most likely hydrogen, with inferred average densities n(H) > 1e6 per cubic centimeter. This could support the observational result of Thi et al. (2001) which indicated the presence of appreciable amounts of H2 around the star beta Pic.
The presence of dusty debris around main sequence stars denotes the existence of planetary systems. Such debris disks are often identified by the presence of excess continuum emission at infrared and (sub-)millimetre wavelengths, with measurements at longer wavelengths tracing larger and cooler dust grains. The exponent of the slope of the disk emission at sub-millimetre wavelengths, `q, defines the size distribution of dust grains in the disk. This size distribution is a function of the rigid strength of the dust producing parent planetesimals. As part of the survey `PLAnetesimals around TYpical Pre-main seqUence Stars (PLATYPUS) we observed six debris disks at 9-mm using the Australian Telescope Compact Array. We obtain marginal (~3-sigma) detections of three targets: HD 105, HD 61005, and HD 131835. Upper limits for the three remaining disks, HD20807, HD109573, and HD109085, provide further constraint of the (sub-)millimetre slope of their spectral energy distributions. The values of q (or their limits) derived from our observations are all smaller than the oft-assumed steady state collisional cascade model (q = 3.5), but lie well within the theoretically expected range for debris disks q ~ 3 to 4. The measured q values for our targets are all < 3.3, consistent with both collisional modelling results and theoretical predictions for parent planetesimal bodies being `rubble piles held together loosely by their self-gravity.