Over the past three decades, we have witnessed one of the great revolutions in our understanding of the cosmos - the dawn of the Exoplanet Era. Where once we knew of just one planetary system (the Solar system), we now know of thousands, with new systems being announced on a weekly basis. Of the thousands of planetary systems we have found to date, however, there is only one that we can study up-close and personal - the Solar system. In this review, we describe our current understanding of the Solar system for the exoplanetary science community - with a focus on the processes thought to have shaped the system we see today. In section one, we introduce the Solar system as a single well studied example of the many planetary systems now observed. In section two, we describe the Solar system's small body populations as we know them today - from the two hundred and five known planetary satellites to the various populations of small bodies that serve as a reminder of the system's formation and early evolution. In section three, we consider our current knowledge of the Solar system's planets, as physical bodies. In section four, we discuss the research that has been carried out into the Solar system's formation and evolution, with a focus on the information gleaned as a result of detailed studies of the system's small body populations. In section five, we discuss our current knowledge of planetary systems beyond our own - both in terms of the planets they host, and in terms of the debris that we observe orbiting their host stars. As we learn ever more about the diversity and ubiquity of other planetary systems, our Solar system will remain the key touchstone that facilitates our understanding and modelling of those newly found systems, and we finish section five with a discussion of the future surveys that will further expand that knowledge.

Resolved images suggest that asymmetric structures are a common feature of cold debris disks. While planets close to these disks are rarely detected, their hidden presence and gravitational perturbations provide plausible explanations for some of these features. To put constraints on the properties of yet undetected planetary companions, we aim to predict what features such a planet imprints in debris disks undergoing continuous collisional evolution. We discuss the basic equations, analytic approximations and timescales governing collisions, radiation pressure and secular perturbations. In addition, we combine our numerical model of the collisional evolution of the size and spatial distributions in debris disks with the gravitational perturbation by a single planet. We find that the distributions of orbital elements in the disks are strongly dependent on grain sizes. Secular precession is differential with respect to involved semi-major axes and grain sizes. This leads to observable differences between the big grains tracing the parent belt and the small grains in the trailing halo. Observations at different wavelengths can be used to constrain the properties of a possible planet.

HR~8799 is a young F0-type star with four directly imaged giant planets and two debris belts, one located exterior and another one interior to the region occupied by the planetary orbits. Having an architecture similar to that of our Solar System, but also revealing dissimilarities such as high masses of planets, a huge radial extent and a high mass of the outer debris belt, HR~8799 is considered to be a benchmark to test formation and evolution models of planetary systems. Here we focus on the outer debris ring and its relation to the planets. We demonstrate that the models of the outer disc, proposed previously to reproduce Herschel observations, are inconsistent with the ALMA data, and vice versa. In an attempt to find a physically motivated model that would agree with both observational sets, we perform collisional simulations. We show that a narrow planetesimal belt and a radiation pressure induced dust halo cannot account for the observed radial brightness profiles. A single, wide planetesimal disc does not reproduce the data either. Instead, we propose a two-population model, comprising a Kuiper-Belt-like structure of a low-eccentricity planetesimal population ("the classical Kuiper Belt") and a high-eccentricity population of comets ("scattered disc"). We argue that such a structure of the exo-Kuiper belt of HR~8799 could be explained with planet migration scenarios analogous to those proposed for the Kuiper Belt of the Solar System.

One of the defining properties of debris discs compared to protoplanetary discs used to be their lack of gas, yet small amounts of gas have been found around an increasing number of debris discs in recent years. These debris discs found to have gas tend to be both young and bright. In this paper we conduct a deep search for CO gas in the system HD 95086 - a 17 Myr old, known planet host that also has a debris disc with a high fractional luminosity of $1.5\times10^{-3}$. Using the Atacama Large Millimeter/submillimeter Array (ALMA) we search for CO emission lines in bands 3, 6 and 7. By implementing a spectro-spatial filtering technique, we find tentative evidence for CO $J$=2-1 emission in the disc located at a velocity, 8.5$\pm$0.2 km s$^{-1}$, consistent with the radial velocity of the star. The tentative detection suggests that the gas on the East side of the disc is moving towards us. In the same region where continuum emission is detected, we find an integrated line flux of 9.5$\pm$3.6 mJy km s$^{-1}$, corresponding to a CO mass of (1.4-13)$\times10^{-6}$ M$_\oplus$. Our analysis confirms that the level of gas present in the disc is inconsistent with the presence of primordial gas in the system and is consistent with second generation production through the collisional cascade.

Context: Since circumstellar dust in debris disks is short-lived, dust-replenishing requires the presence of a reservoir of planetesimals. These planetesimals in the parent belt of debris disks orbit their host star and continuously supply the disk with fine dust through their mutual collisions. Aims: We aim to understand effects of different collisional parameters on the observational appearance of eccentric debris disks. Methods: The collisional evolution of selected debris disk configurations was simulated with the numerical code ACE. Subsequently, selected observable quantities are simulated with our newly developed code DMS. The impact of the eccentricity, dynamical excitation, and the material strength is discussed with respect to the grain size distribution, the spectral energy distribution, and spatially resolved images of debris disk systems. Results: The most recognizable features in different collisional evolutions are as follows. First, both the increase of dynamical excitation in the eccentric belt of the debris disk system and the decrease of the material strength of dust particles result in a higher production rate of smaller particles. This reduces the surface brightness differences between the periastron and the apastron sides of the disks. For very low material strengths, the "pericenter glow" phenomenon is reduced and eventually even replaced by the opposite effect, the "apocenter glow". Second, it is possible to constrain the level of collisional activity from the appearance of the disk, for example, the wavelength-dependent apocenter-to-pericenter flux ratio. Within the considered parameter space, the impact of the material strength on the appearance of the disk is stronger than that of dynamical excitation of the system. Finally, we find that the impact of the collisional parameters on the net spectral energy distribution is weak.

Two basic routes for planetesimal formation have been proposed over the last few decades. One is a classical "slow-growth" scenario. Another one is particle concentration models, in which small pebbles are concentrated locally and then collapse gravitationally to form planetesimals. Both types of models make certain predictions for the size spectrum and internal structure of newly-born planetesimals. We use these predictions as input to simulate collisional evolution of debris discs left after the gas dispersal. The debris disc emission as a function of a system's age computed in these simulations is compared with several Spitzer and Herschel debris disc surveys around A-type stars. We confirm that the observed brightness evolution for the majority of discs can be reproduced by classical models. Further, we find that it is equally consistent with the size distribution of planetesimals predicted by particle concentration models - provided the objects are loosely bound "pebble piles" as these models also predict. Regardless of the assumed planetesimal formation mechanism, explaining the brightest debris discs in the samples uncovers a "disc mass problem." To reproduce such discs by collisional simulations, a total mass of planetesimals of up to ~1000 Earth masses is required, which exceeds the total mass of solids available in the protoplanetary progenitors of debris discs. This may indicate that stirring was delayed in some of the bright discs, that giant impacts occurred recently in some of them, that some systems may be younger than previously thought, or that non-collisional processes contribute significantly to the dust production.

Context: High-resolution images of circumstellar debris discs reveal off-centred rings that indicate past or ongoing perturbation, possibly caused by secular gravitational interaction with unseen stellar or substellar companions. The purely dynamical aspects of this departure from radial symmetry are well understood. However, the observed dust is subject to additional forces and effects, most notably collisions and drag. Aims: To complement the studies of dynamics, we therefore aim to understand how new asymmetries are created by the addition of collisional evolution and drag forces, and existing ones strengthened or overridden. Methods: We augmented our existing numerical code "Analysis of Collisional Evolution" (ACE) by an azimuthal dimension, the longitude of periapse. A set of fiducial discs with global eccentricities ranging from 0 to 0.4 is evolved over giga-year timescales. Size distribution and spatial variation of dust are analysed and interpreted. The basic impact of belt eccentricity on spectral energy distributions (SEDs) and images is discussed. Results: We find features imposed on characteristic timescales. First, radiation pressure defines size cutoffs that differ between periapse and apoapse, resulting in an asymmetric halo. The differences in size distribution make the observable asymmetry of the halo depend on wavelength. Second, collisional equilibrium prefers smaller grains on the apastron side of the parent belt, reducing the effect of pericentre glow and the overall asymmetry. Third, Poynting-Robertson drag fills the region interior to an eccentric belt such that the apastron side is more tenuous. Interpretation and prediction of the appearance in scattered light is problematic when spatial and size distribution are coupled.

Many debris discs reveal a two-component structure, with an outer Kuiper-belt analogue and a warm inner component whose origin is still a matter of debate. One possibility is that warm emission stems from an "asteroid belt" closer in to the star. We consider a scenario in which a set of giant planets is formed in an initially extended planetesimal disc. These planets carve a broad gap around their orbits, splitting up the disc into the outer and the inner belts. After the gas dispersal, both belts undergo collisional evolution in a steady-state regime. This scenario is explored with detailed collisional simulations involving realistic physics to describe a long-term collisional depletion of the two-component disc. We find that the inner disc may be able to retain larger amounts of material at older ages than thought before on the basis of simplified analytic models. We show that the proposed scenario is consistent with a suite of thermal emission and scattered light observational data for a bright two-temperature debris disc around a nearby solar-type star q$^1$ Eridani. This implies a Solar System-like architecture of the system, with an outer massive "Kuiper belt", an inner "asteroid belt", and a few Neptune- to Jupiter-mass planets in between.

The spatially resolved AU Mic debris disc is among the most famous and best-studied debris discs. We aim at a comprehensive understanding of the dust production and the dynamics of the disc objects with in depth collisional modelling including stellar radiative and corpuscular forces. Our models are compared to a suite of observational data for thermal and scattered light emission, ranging from the ALMA radial surface brightness profile at 1.3mm to polarisation measurements in the visible. Most of the data can be reproduced with a planetesimal belt having an outer edge at around 40au and subsequent inward transport of dust by stellar winds. A low dynamical excitation of the planetesimals with eccentricities up to 0.03 is preferred. The radial width of the planetesimal belt cannot be constrained tightly. Belts that are 5au and 17au wide, as well as a broad 44au-wide belt are consistent with observations. All models show surface density profiles increasing with distance from the star as inferred from observations. The best model is achieved by assuming a stellar mass loss rate that exceeds the solar one by a factor of 50. While the SED and the shape of the ALMA profile are well reproduced, the models deviate from the scattered light data more strongly. The observations show a bluer disc colour and a lower degree of polarisation for projected distances <40au than predicted by the models. The problem may be mitigated by irregularly-shaped dust grains which have scattering properties different from the Mie spheres used. From tests with a handful of selected dust materials, we derive a preference for mixtures of silicate, carbon, and ice of moderate porosity. We address the origin of the unresolved central excess emission detected by ALMA and show that it cannot stem from an additional inner belt alone. Instead, it should derive, at least partly, from the chromosphere of the central star.

The optical properties of ice in the far infrared are important for models of protoplanetary and debris disks. In this report we derive a new set of data for the absorption (represented by the imaginary part of the refractive index $\kappa$) of crystalline water ice in this spectral range, including a detailed inspection of the temperature dependence, which had not been done in such detail before. We measured the transmission of three ice layers with different thicknesses at temperatures $\vartheta = 10...250$K and present data at wavelengths $\lambda=80...625$ microns. We found a change in the spectral dependence of $\kappa$ at a wavelength of $175 \pm 6$ microns. At shorter wavelengths, $\kappa$ exhibits a constant flat slope and no significant temperature dependence. Long-ward of that wavelength, the slope gets steeper and has a clear, approximately linear temperature dependence. This change in the behaviour is probably caused by a characteristic absorption band of water ice. The measured data were fitted by a power-law model that analytically describes the absorption behaviour at an arbitrary temperature. This model can readily be applied to any object of interest, for instance a protoplanetary or a debris disk. To illustrate how the model works, we simulated the spectral energy distribution (SED) of the resolved, large debris disk around the nearby solar-type star HD 207129. Replacing our ice model by another, commonly used data set for water ice results in a different SED slope at longer wavelengths. This leads to changes in the characteristic model parameters of the disk, such as the inferred particle size distribution, and affects the interpretation of the underlying collisional physics of the disk.

We present an analysis of the debris disc around the nearby K2 V star HIP 17439. In the context of the Herschel DUNES key programme the disc was observed and spatially resolved in the far-IR with the Herschel PACS and SPIRE instruments. In a first model, Ertel et al. (2014) assumed the size and radial distribution of the circumstellar dust to be independent power laws. There, by exploring a very broad range of possible model parameters several scenarios capable of explaining the observations were suggested. In this paper, we perform a follow-up in-depth collisional modelling of these scenarios trying to further distinguish between them. In our models we consider collisions, direct radiation pressure, and drag forces, i.e. the actual physical processes operating in debris discs. We find that all scenarios discussed in Ertel et al. are physically sensible and can reproduce the observed SED along with the PACS surface brightness profiles reasonably well. In one model, the dust is produced beyond 120au in a narrow planetesimal belt and is transported inwards by Poynting-Robertson and stellar wind drag. A good agreement with the observed radial profiles would require stellar winds by about an order of magnitude stronger than the solar value, which is not supported, although not ruled out, by observations. Another model consists of two spatially separated planetesimal belts, a warm inner and a cold outer one. This scenario would probably imply the presence of planets clearing the gap between the two components. Finally, we show qualitatively that the observations can be explained by assuming the dust is produced in a single, but broad planetesimal disc with a surface density of solids rising outwards, as expected for an extended disc that experiences a natural inside-out collisional depletion. Prospects of discriminating between the competing scenarios by future observations are discussed.

Debris disks are thought to be sculptured by neighboring planets. The same is true for the Edgeworth-Kuiper debris disk, yet no direct observational evidence for signatures of giant planets in the Kuiper belt dust distribution has been found so far. Here we model the dust distribution in the outer solar system to reproduce the dust impact rates onto the dust detector onboard the New Horizons spacecraft measured so far and to predict the rates during the Neptune orbit traverse. To this end, we take a realistic distribution of transneptunian objects to launch a sufficient number of dust grains of different sizes and follow their orbits by including radiation pressure, Poynting-Robertson and stellar wind drag, as well as the perturbations of four giant planets. In a subsequent statistical analysis, we calculate number densities and lifetimes of the dust grains in order to simulate a collisional cascade. In contrast to the previous work, our model not only considers collisional elimination of particles, but also includes production of finer debris. We find that particles captured in the 3:2 resonance with Neptune build clumps that are not removed by collisions, because the depleting effect of collisions is counteracted by production of smaller fragments. Our model successfully reproduces the dust impact rates measured by New Horizons out to ~23AU and predicts an increase of the impact rate of about a factor of two or three around the Neptune orbit crossing. This result is robust with respect to the variation of the vaguely known number of dust-producing scattered disk objects, collisional outcomes, and the dust properties.

(abridged) Infrared excesses associated with debris disk host stars detected so far peak at wavelengths around ~100{\mu}m or shorter. However, six out of 31 excess sources in the Herschel OTKP DUNES have been seen to show significant - and in some cases extended - excess emission at 160{\mu}m, which is larger than the 100{\mu}m excess. This excess emission has been suggested to stem from debris disks colder than those known previously. Using several methods, we re-consider whether some or even all of the candidates may be associated with unrelated galactic or extragalactic emission and conclude that it is highly unlikely that none of the candidates represents a true circumstellar disk. For true disks, both the dust temperatures inferred from the SEDs and the disk radii estimated from the images suggest that the dust is nearly as cold as a blackbody. This requires the grains to be larger than ~100{\mu}m, regardless of their material composition. To explain the dearth of small grains, we explore several conceivable scenarios: transport-dominated disks, disks of low dynamical excitation, and disks of unstirred primordial macroscopic grains. Our qualitative analysis and collisional simulations rule out the first two of these scenarios, but show the feasibility of the third one. We show that such disks can survive for gigayears, largely preserving the primordial size distribution. They should be composed of macroscopic solids larger than millimeters, but smaller than kilometers in size. Thus planetesimal formation, at least in the outer regions of the systems, has stopped before "cometary" or "asteroidal" sizes were reached.

Debris discs are a consequence of the planet formation process and constitute the fingerprints of planetesimal systems. Their solar system's counterparts are the asteroid and Edgeworth-Kuiper belts. The DUNES survey aims at detecting extra-solar analogues to the Edgeworth-Kuiper belt around solar-type stars, putting in this way the solar system into context. The survey allows us to address some questions related to the prevalence and properties of planetesimal systems. We used {\it Herschel}/PACS to observe a sample of nearby FGK stars. Data at 100 and 160 $\mu$m were obtained, complemented in some cases with observations at 70 $\mu$m, and at 250, 350 and 500 $\mu$m using SPIRE. The observing strategy was to integrate as deep as possible at 100 $\mu$m to detect the stellar photosphere. Debris discs have been detected at a fractional luminosity level down to several times that of the Edgeworth-Kuiper belt. The incidence rate of discs around the DUNES stars is increased from a rate of $\sim$ 12.1% $\pm$ 5% before \emph{Herschel} to $\sim$ 20.2% $\pm$ 2%. A significant fraction ($\sim$ 52%) of the discs are resolved, which represents an enormous step ahead from the previously known resolved discs. Some stars are associated with faint far-IR excesses attributed to a new class of cold discs. Although it cannot be excluded that these excesses are produced by coincidental alignment of background galaxies, statistical arguments suggest that at least some of them are true debris discs. Some discs display peculiar SEDs with spectral indexes in the 70-160$\mu$m range steeper than the Rayleigh-Jeans one. An analysis of the debris disc parameters suggests that a decrease might exist of the mean black body radius from the F-type to the K-type stars. In addition, a weak trend is suggested for a correlation of disc sizes and an anticorrelation of disc temperatures with the stellar age.

(Abridged) We access the expected EKB dust disk properties by modeling. We treat the debiased population of the known transneptunian objects (TNOs) as parent bodies and generate the dust with our collisional code. The resulting dust distributions are modified to take into account the influence of gravitational scattering and resonance trapping by planets on migrating dust grains as well as the effect of sublimation. A difficulty is that the amount and distribution of dust are largely determined by sub-kilometer-sized bodies. These are directly unobservable, and their properties cannot be accessed by collisional modeling, because objects larger than 10...60m in the present-day EKB are not in a collisional equilibrium. To place additional constraints, we use in-situ measurements of the New Horizons spacecraft within 20AU. We show that the TNO population has to have a break in the size distribution at s<70km. However, even this still leaves us with several models that all correctly reproduce a nearly constant dust impact rates in the region of giant planet orbits and do not violate the constraints from the non-detection of the EKB dust thermal emission by the COBE spacecraft. The modeled EKB dust disks, which conform to the observational constraints, can either be transport-dominated or intermediate between the transport-dominated and collision-dominated regime. The in-plane optical depth of such disks is tau(r>10AU)~10^-6 and their fractional luminosity is f_d~10^-7. Planets and sublimation are found to have little effect on dust impact fluxes and dust thermal emission. The spectral energy distribution of an EKB analog, as would be seen from 10pc distance, peaks at wavelengths of 40...50\mum at F~0.5mJy, which is less than 1% of the photospheric flux at those wavelengths. Therefore, exact EKB analogs cannot be detected with present-day instruments such as Herschel/PACS.

Debris discs are known to exist around many planet-host stars, but no debris dust has been found so far in systems with transiting planets. Using publicly available catalogues, we searched for infrared excesses in such systems. In the recently published Wide-Field Infrared Survey Explorer (WISE) catalogue, we found 52 stars with transiting planets. Two systems with one transiting "hot Jupiter" each, TrES-2 and XO-5, exhibit small excesses both at 12 and 22 microns at a > 3 sigma level. Provided that one or both of these detections are real, the frequency of warm excesses in systems with transiting planets of 2-4 % is comparable to that around solar-type stars probed at similar wavelengths with Spitzer's MIPS and IRS instruments. Modelling suggests that the observed excesses would stem from dust rings with radii of several AU. The inferred amount of dust is close to the maximum expected theoretically from a collisional cascade in asteroid belt analogues. If confirmed, the presence of debris discs in systems with transiting planets may put important constraints onto formation and migration scenarios of hot Jupiters.

Dusty debris discs around main sequence stars are thought to be the result of continuous collisional grinding of planetesimals in the system. The majority of these systems are unresolved and analysis of the dust properties is limited by the lack of information regarding the dust location.vThe Herschel DUNES key program is observing 133 nearby, Sun-like stars (<20 pc, FGK spectral type) in a volume limited survey to constrain the absolute incidence of cold dust around these stars by detection of far infrared excess emission at flux levels comparable to the Edgeworth-Kuiper belt (EKB). We have observed the Sun-like star HD 207129 with Herschel PACS and SPIRE. In all three PACS bands we resolve a ring-like structure consistent with scattered light observations. Using {\alpha} Bo\"otis as a reference point spread function (PSF), we deconvolved the images, clearly resolving the inner gap in the disc at both 70 and 100 {\mu}m. We have resolved the dust-producing planetesimal belt of a debris disc at 100 {\mu}m for the first time. We measure the radial profile and fractional luminosity of the disc, and compare the values to those of discs around stars of similar age and/or spectral type, placing this disc in context of other resolved discs observed by Herschel/DUNES.

(Abridged) The Edgeworth-Kuiper belt with its presumed dusty debris is a natural reference for extrsolar debris disks. We employ a new algorithm to eliminate the inclination and the distance selection effects in the known TNO populations to derive expected parameters of the "true" EKB. Its estimated mass is M_EKB=0.12 M_earth, which is by a factor of \sim 15 larger than the mass of the EKB objects detected so far. About a half of the total EKB mass is in classical and resonant objects and another half is in scattered ones. Treating the debiased populations of EKB objects as dust parent bodies, we then "generate" their dust disk with our collisional code. Apart from accurate handling of collisions and direct radiation pressure, we include the Poynting-Robertson (P-R) drag, which cannot be ignored for the EKB dust disk. Outside the classical EKB, the radial profile of the optical depth approximately follows tau \sim r^-2 which is roughly intermediate between the slope predicted analytically for collision-dominated (r^-1.5) and transport-dominated (r^-2.5) disks. The cross section-dominating grain size still lies just above the blowout size (\sim 1...2 \microm), as it would without the P-R transport. However, if the EKB were by one order of magnitude less massive, the optical depth profile would fall off as tau \sim r^-3, and the cross section-dominating grain size would shift from \sim 1...2\microm to ~100 \microm. These properties are seen if dust is assumed to be generated only by known TNOs. If the solar system were observed from outside, the thermal emission flux from the EKB dust would be about two orders of magnitude lower than for solar-type stars with the brightest known infrared excesses observed from the same distance. Herschel and other new-generation facilities should reveal extrasolar debris disks nearly as tenuous as the EKB disk. The Herschel/PACS instrument should be able to detect disks at a \sim 1...2M_EKB level.

We present the first far-IR observations of the solar-type stars delta Pav, HR 8501, 51 Peg and zeta^2 Ret, taken within the context of the DUNES Herschel Open Time Key Programme (OTKP). This project uses the PACS and SPIRE instruments with the objective of studying infrared excesses due to exo-Kuiper belts around nearby solar-type stars. The observed 100 um fluxes from delta Pav, HR 8501, and 51 Peg agree with the predicted photospheric fluxes, excluding debris disks brighter than Ldust/Lstar ~ 5 x 10^-7 (1 sigma level) around those stars. A flattened, disk-like structure with a semi-major axis of ~ 100 AU in size is detected around zeta^2 Ret. The resolved structure suggests the presence of an eccentric dust ring, which we interpret as an exo-Kuiper belt with Ldust/Lstar ~ 10^-5.

About two dozen exo-solar debris systems have been spatially resolved. These debris discs commonly display a variety of structural features such as clumps, rings, belts, eccentric distributions and spiral patterns. In most cases, these features are believed to be formed, shaped and maintained by the dynamical influence of planets orbiting the host stars. In very few cases has the presence of the dynamically important planet(s) been inferred from direct observation. The solar-type star q1 Eri is known to be surrounded by debris, extended on scales of < 30''. The star is known to host at least one planet, albeit on an orbit far too small to make it responsible for structures at distances of tens to hundreds of AU. The aim of the present investigation is twofold: to determine the optical and material properties of the debris and to infer the spatial distribution of the dust, which may hint at the presence of additional planets. The photodetector array camera and spectrometer (PACS) aboard the Herschel Space Observatory allows imaging observations in the far infrared at unprecedented resolution, i.e. at better than 6'' to 12'' over the wavelength range of 60 {\mu}m to 210 {\mu}m. Together with the results from ground-based observations, these spatially resolved data can be modelled to determine the nature of the debris and its evolution more reliably than would be possible from unresolved data alone. For the first time has the q1 Eri disc been resolved at far infrared wavelengths. The PACS observations at 70, 100 and 160 {\mu}m reveal an oval image showing a disc-like structure in all bands, the size of which increases with wavelength. Assuming a circular shape yields the inclination of its equatorial plane with respect to that of the sky, i > 53deg. The results of image de-convolution indicate that i likely is larger than 63deg, where 90deg corresponds to an edge-on disc. {abridged}

It has been argued that the photometric data and images of the archetypical debris disk around Vega may be in contradiction with the standard, steady-state collisional scenario of the disk evolution. Here we perform physical modeling of the Vega disk "from the sources". We assume that dust is maintained by a "Kuiper belt" of parent planetesimals at ~ 100 AU and employ our collisional and radiative transfer codes to consistently model the size and radial distribution of the disk material and then thermal emission of dust. In doing so, we vary a broad set of parameters, including the stellar properties, the exact location, extension, and dynamical excitation of the planetesimal belt, chemical composition of solids, and the collisional prescription. We are able to reproduce the spectral energy distribution in the entire wavelength range from the near-infrared to millimeter, as well as the mid-IR and sub-millimeter radial brightness profiles of the Vega disk. Thus our results suggest that the Vega disk observations are compatible with a steady-state collisional dust production, and we put important constraints on the disk parameters and physical processes that sustain it. The total disk mass in < 100 km-sized bodies is estimated to be ~ 10 Earth masses. Provided that collisional cascade has been operating over much of the Vega age of ~ 350 Myr, the disk must have lost a few Earth masses of solids during that time. We also demonstrate that using an intermediate luminosity of the star between the pole and the equator, as derived from its fast rotation, is required to reproduce the debris disk observations. Finally, we show that including cratering collisions into the model is mandatory.

HR8799 is a nearby A-type star with a debris disk and three planetary candidates recently imaged directly. We undertake a coherent analysis of various portions of observational data on all known components of the system. The goal is to elucidate the architecture and evolutionary status of the system. We try to further constrain the age and orientation of the system, orbits and masses of the companions, as well as the location of dust. From the high luminosity of debris dust and dynamical constraints, we argue for a rather young system's age of <50Myr. The system must be seen nearly, but not exactly, pole-on. Our analysis of the stellar rotational velocity yields an inclination of 13-30deg, whereas i>20deg is needed for the system to be dynamically stable, which suggests a probable inclination range of 20-30deg. The spectral energy distribution is naturally reproduced with two dust rings associated with two planetesimal belts. The inner "asteroid belt" is located at ~10AU inside the orbit of the innermost companion and a "Kuiper belt" at >100AU is just exterior to the orbit of the outermost companion. The dust masses in the inner and outer ring are estimated to be ~1E-05 and 4E-02 M_earth, respectively. We show that all three planetary candidates may be stable in the mass range suggested in the discovery paper by Marois et al. 2008 (between 5 and 13 Jupiter masses), but only for some of all possible orientations. Stable orbits imply a double (4:2:1) mean-motion resonance between all three companions. We finally show that in the cases where the companions themselves are orbitally stable, the dust-producing planetesimal belts are also stable against planetary perturbations.

Debris disks around main-sequence stars are believed to derive from planetesimal populations that have accreted at early epochs and survived possible planet formation processes. While debris disks must contain solids in a broad range of sizes - from big planetesimals down to tiny dust grains - debris disk observations are only sensitive to the dust end of the size distribution. Collisional models of debris disks are needed to "climb up" the ladder of the collisional cascade, from dust towards parent bodies, representing the main mass reservoir of the disks. We have used our collisional code to generate five disks around a sun-like star, assuming planetesimal belts at 3, 10, 30, 100, and 200 AU with 10 times the Edgeworth-Kuiper-belt mass density, and to evolve them for 10 Gyr. Along with an appropriate scaling rule, this effectively yields a three-parametric set of reference disks (initial mass, location of planetesimal belt, age). For all the disks, we have generated spectral energy distributions (SEDs), assuming homogeneous spherical astrosilicate dust grains. A comparison between generated and actually observed SEDs yields estimates of planetesimal properties (location, total mass etc.). As a test and a first application of this approach, we have selected five disks around sun-like stars with well-known SEDs. In four cases, we have reproduced the data with a linear combination of two disks from the grid (an "asteroid belt" at 3 AU and an outer "Kuiper belt"); in one case a single, outer component was sufficient. The outer components are compatible with "large Kuiper belts" of 0.2-50 earth masses (in the bodies up to 100 km in size) with radii of 100-200 AU.

We simulated the long-term collisional depletion of debris disks around solar-type (G2V) stars with our code. The numerical results were supplemented by, and interpreted through, a new analytic model. A few general scaling rules for the disk evolution are suggested. The timescale of the collisional evolution is inversely proportional to the initial disk mass and scales with radial distance as r^4.3 and with eccentricities of planetesimals as e^-2.3. Further, we show that at actual ages of debris disks between 10 Myr and 10 Gyr, the decay of the dust mass and the total disk mass follow different laws. The reason is that the collisional lifetime of planetesimals is size-dependent. At any moment, there exists a transitional size, which separates larger objects that still have the ``primordial'' size distribution set in the growth phase from small objects whose size distribution is already set by disruptive collisions. The dust mass and its decay rate evolve as that transition affects objects of ever-larger sizes. Under standard assumptions, the dust mass, fractional luminosity, and thermal fluxes all decrease as t^xi with xi = -0.3...-0.4. Specific decay laws of the total disk mass and the dust mass, including the value of xi, largely depend on a few model parameters, such as the critical fragmentation energy as a function of size, the primordial size distribution of largest planetesimals, as well as the characteristic eccentricity and inclination of their orbits. With standard material prescriptions and a distribution of disk masses and extents, a synthetic population of disks generated with our analytic model agrees quite well with the observed Spitzer/MIPS statistics of 24 and 70 micron fluxes and colors versus age.

Observations have revealed a large variety of structures (global asymmetries, warps, belts, rings) and dynamical phenomena ("falling-evaporating bodies" or FEBs, the "Beta Pic dust stream") in the disc of Beta Pictoris, most of which may indicate the presence of one or more planets orbiting the star. Because planets of Beta Pic have not been detected by observations yet, we use dynamical simulations to find "numerical evidence" for a planetary system. We show that already one planet at 12 AU with a mass of 2 to 5 Mjup and an eccentricity smaller or equal 0.1 can probably account for three major features (main warp, two inner belts, FEBs) observed in the Beta Pic disc. The existence of at least two additional planets at about 25 AU and 45 AU from the star seems likely. We find rather strong upper limits of 0.6 Mjup and 0.2 Mjup on the masses of those planets. The same planets could, in principle, also account for the outer rings observed at 500 - 800 AU.

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.