Debris discs are second generation dusty discs formed by collisions of planetesimals. Many debris discs have been found and resolved around hot and solar-type stars. However, only a handful have been discovered around M-stars, and the reasons for their paucity remain unclear. Here we check whether the sensitivity and wavelength coverage of present-day telescopes are simply unfavourable for detection of these discs or if they are truly rare. We approach this question by looking at the Herschel/DEBRIS survey that has searched for debris discs including M-type stars. Assuming that these cool-star discs are "similar" to those of the hotter stars in some sense (i.e., in terms of dust location, temperature, fractional luminosity, or mass), we check whether this survey should have found them. With our procedure we can reproduce the $2.1^{+4.5}_{-1.7}$% detection rate of M-star debris discs of the DEBRIS survey, which implies that these discs can indeed be similar to discs around hotter stars and just avoid detection. We then apply this procedure to IRAM NIKA-2 and ALMA bands 3, 6 and 7 to predict possible detection rates and give recommendations for future observations. We do not favour observing with IRAM, since it leads to detection rates lower than for the DEBRIS survey, with 0.6%-4.5% for a 15 min observation. ALMA observations, with detection rates 0.9%-7.4%, do not offer a significant improvement either, and so we conclude that more sensitive far-infrared and single dish sub-millimetre telescopes are necessary to discover the missing population of M-star debris discs.

Debris discs are dusty belts of planetesimals around main-sequence stars, similar to the asteroid and Kuiper belts in our solar system. The planetesimals cannot be observed directly, yet they produce detectable dust in mutual collisions. Observing the dust, we can try to infer properties of invisible planetesimals. Here we address the question of what is the best way to measure the location of outer planetesimal belts that encompass extrasolar planetary systems. A standard method is using resolved images at mm-wavelengths, which reveal dust grains with sizes comparable to the observational wavelength. Smaller grains seen in the infrared (IR) are subject to several non-gravitational forces that drag them away from their birth rings, and so may not closely trace the parent bodies. In this study, we examine whether imaging of debris discs at shorter wavelengths might enable determining the spatial location of the exo-Kuiper belts with sufficient accuracy. We find that around M-type stars the dust best visible in the mid-IR is efficiently displaced inward from their birth location by stellar winds, causing the discs to look more compact in mid-IR images than they actually are. However, around earlier-type stars where the majority of debris discs is found, discs are still the brightest at the birth ring location in the mid-IR regime. Thus, sensitive IR facilities with good angular resolution, such as MIRI on JWST, will enable tracing exo-Kuiper belts in nearby debris disc systems.

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.

When a protoplanetary disc loses gas, it leaves behind planets and one or more planetesimal belts. The belts get dynamically excited, either by planets ('planet stirring') or by embedded big planetesimals ('self-stirring'). Collisions between planetesimals become destructive and start to produce dust, creating an observable debris disc. Following Kenyon & Bromley (2008), it is often assumed that self-stirring starts to operate as soon as the first ~1000 km-sized embedded 'Plutos' have formed. However, state-of-the-art pebble concentration models robustly predict planetesimals between a few km and ~200 km in size to form in protoplanetary discs rapidly, before then slowly growing into Pluto-sized bodies. We show that the timescale, on which these planetesimals excite the disc sufficiently for fragmentation, is shorter than the formation timescale of Plutos. Using an analytic model based on the Ida & Makino (1993) theory, backed up by numerical integrations, we find the excitation timescale to be $T_{\rm excite} \approx 100 x_{\rm m}^{-1} M_*^{-3/2} a^3$ Myr, where $x_{\rm m}$ is the total mass of a protoplanetary disc progenitor in the units of the Minimum-Mass Solar Nebula, $a$ its radius in the units of 100 AU, and $M_*$ is the stellar mass in solar masses. These results are applied to a set of 23 debris discs that have been well resolved with ALMA or SMA. We find that the majority of these discs are consistent with being self-stirred. However, three large discs around young early-type stars do require planets as stirrers. These are 49 Cet, HD 95086, and HR 8799, of which the latter two are already known to have planets.

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.

A warm/hot dust component (at temperature $>$ 300K) has been detected around $\sim$ 20% of stars. This component is called "exozodiacal dust" as it presents similarities with the zodiacal dust detected in our Solar System, even though its physical properties and spatial distribution can be significantly different. Understanding the origin and evolution of this dust is of crucial importance, not only because its presence could hamper future detections of Earth-like planets in their habitable zones, but also because it can provide invaluable information about the inner regions of planetary systems. In this review, we present a detailed overview of the observational techniques used in the detection and characterisation of exozodiacal dust clouds ("exozodis") and the results they have yielded so far, in particular regarding the incidence rate of exozodis as a function of crucial parameters such as stellar type and age, or the presence of an outer cold debris disc. We also present the important constraints that have been obtained, on dust size distribution and spatial location, by using state-of-the-art radiation transfer models on some of these systems. Finally, we investigate the crucial issue of how to explain the presence of exozodiacal dust around so many stars (regardless of their ages) despite the fact that such dust so close to its host star should disappear rapidly due to the coupled effect of collisions and stellar radiation pressure. Several potential mechanisms have been proposed to solve this paradox and are reviewed in detail in this paper. The review finishes by presenting the future of this growing field.

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.

Hot exozodiacal dust emission was detected in recent surveys around two dozen main sequence stars at distances of less than $1\,\text{au}$ using H and K band interferometry. Due to the high contrast as well as the small angular distance between the circumstellar dust and the star, direct observation of this dust component is challenging. An alternative way to explore the hot exozodiacal dust is provided by mid-infrared interferometry. We analyze the L, M and N band interferometric signature of this emission in order to find stronger constraints for the properties and the origin of the hot exozodiacal dust. Considering the parameters of nine debris disc systems derived previously, we model the discs in each of these bands. We find that the M band possesses the best conditions to detect hot dust emission, closely followed by L and N bands. The hot dust in three systems - HD 22484 (10 Tau), HD 102647 ($\beta$ Leo) and HD 177724 ($\zeta$ Aql) - shows a strong signal in the visibility functions which may even allow one to constrain the dust location. In particular, observations in the mid-infrared could help to determine whether the dust piles up at the sublimation radius or is located at radii up to $1\,\text{au}$. In addition, we explore observations of the hot exozodiacal dust with the upcoming mid-infrared interferometer MATISSE at the Very Large Telescope Interferometer (VLTI).

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.

Debris disk morphology is wavelength dependent due to the wide range of particle sizes and size-dependent dynamics influenced by various forces. Resolved images of nearby debris disks reveal complex disk structures that are difficult to distinguish from their spectral energy distributions. Therefore, multi-wavelength resolved images of nearby debris systems provide an essential foundation to understand the intricate interplay between collisional, gravitational, and radiative forces that govern debris disk structures. We present the SOFIA 35 um resolved disk image of epsilon Eri, the closest debris disk around a star similar to the early Sun. Combining with the Spitzer resolved image at 24 um and 15-38 um excess spectrum, we examine two proposed origins of the inner debris in epsilon Eri: (1) in-situ planetesimal belt(s) and (2) dragged-in grains from the cold outer belt. We find that the presence of in-situ dust-producing planetesmial belt(s) is the most likely source of the excess emission in the inner 25 au region. Although a small amount of dragged-in grains from the cold belt could contribute to the excess emission in the inner region, the resolution of the SOFIA data is high enough to rule out the possibility that the entire inner warm excess results from dragged-in grains, but not enough to distinguish one broad inner disk from two narrow belts.

Many debris discs reveal a two-component structure, with a cold outer and a warm inner component. While the former are likely massive analogues of the Kuiper belt, the origin of the latter is still a matter of debate. In this work we investigate whether the warm dust may be a signature of asteroid belt analogues. In the scenario tested here the current two-belt architecture stems from an originally extended protoplanetary disc, in which planets have opened a gap separating it into the outer and inner discs which, after the gas dispersal, experience a steady-state collisional decay. This idea is explored with an analytic collisional evolution model for a sample of 225 debris discs from a Spitzer/IRS catalogue that are likely to possess a two-component structure. We find that the vast majority of systems (220 out of 225, or 98%) are compatible with this scenario. For their progenitors, original protoplanetary discs, we find an average surface density slope of $-0.93\pm0.06$ and an average initial mass of $\left(3.3^{+0.4}_{-0.3}\right)\times 10^{-3}$ solar masses, both of which are in agreement with the values inferred from submillimetre surveys. However, dust production by short-period comets and - more rarely - inward transport from the outer belts may be viable, and not mutually excluding, alternatives to the asteroid belt scenario. The remaining five discs (2% of the sample: HIP 11486, HIP 23497, HIP 57971, HIP 85790, HIP 89770) harbour inner components that appear inconsistent with dust production in an "asteroid belt." Warm dust in these systems must either be replenished from cometary sources or represent an aftermath of a recent rare event, such as a major collision or planetary system instability.

Recent interferometric surveys of nearby main-sequence stars show a faint but significant near-infrared excess in roughly two dozen systems, i.$\,$e. around $10\,\%$ to $30\,\%$ of stars surveyed. This excess is attributed to dust located in the immediate vicinity of the star, the origin of which is highly debated. We used previously published interferometric observations to constrain the properties and distribution of this hot dust. Considering both scattered radiation and thermal reemission, we modelled the observed excess in nine of these systems. We find that grains have to be sufficiently absorbing to be consistent with the observed excess, while dielectric grains with pure silicate compositions fail to reproduce the observations. The dust should be located within $\sim\,0.01-1\,au$ from the star depending on its luminosity. Furthermore, we find a significant trend for the disc radius to increase with the stellar luminosity. The dust grains are determined to be below $0.2-0.5\,\mu\textrm{m}$, but above $0.02-0.15\,\mu\textrm{m}$ in radius. The dust masses amount to $(0.2-3.5)\times10^{-9}\,M_\oplus$. The near-infrared excess is probably dominated by thermal reemission, though a contribution of scattered light up to $35\,\%$ cannot be completely excluded. The polarisation degree predicted by our models is always below $5\,\%$, and for grains smaller than $\sim0.2\,\mu\textrm{m}$ even below $1\,\%$. We also modelled the observed near-infrared excess of another ten systems with poorer data in the mid-infrared. The basic results for these systems appear qualitatively similar, yet the constraints on the dust location and the grain sizes are weaker.

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.

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 aim of this paper is to provide robust numbers for the incidence of debris discs around FGK stars in the solar neighbourhood. The full sample of 177 FGK stars with d<20 pc proposed for the DUNES survey is presented. Herschel/PACS observations at 100 and 160 micron complemented with data at 70 micron, and at 250, 350 and 500 micron SPIRE photometry, were obtained. The 123 objects observed by the DUNES collaboration were presented in a previous paper. The remaining 54 stars, shared with the DEBRIS consortium and observed by them, and the combined full sample are studied in this paper. The incidence of debris discs per spectral type is analysed and put into context together with other parameters of the sample, like metallicity, rotation and activity, and age. The subsample of 105 stars with d<15 pc containing 23 F, 33 G and 49 K stars, is complete for F stars, almost complete for G stars and contains a substantial number of K stars to draw solid conclusions on objects of this spectral type. The incidence rates of debris discs per spectral type 0.26 (6 objects with excesses out of 23 F stars), 0.21 (7 out of 33 G stars) and 0.20 (10 out of 49 K stars), the fraction for all three spectral types together being 0.22 (23 out of 105 stars). Uncertainties corresponding to a 95% confidence level are given in the text for all these numbers. The medians of the upper limits of L_dust/L_* for each spectral type are 7.8E-7 (F), 1.4E-6 (G) and 2.2E-6 (K); the lowest values being around 4.0E-7. The incidence of debris discs is similar for active (young) and inactive (old) stars. The fractional luminosity tends to drop with increasing age, as expected from collisional erosion of the debris belts.

The cross section of material in debris discs is thought to be dominated by the smallest grains that can still stay in bound orbits despite the repelling action of stellar radiation pressure. Thus the minimum (and typical) grain size $s_\text{min}$ is expected to be close to the radiation pressure blowout size $s_\text{blow}$. Yet a recent analysis of a sample of Herschel-resolved debris discs showed the ratio $s_\text{min}/s_\text{blow}$ to systematically decrease with the stellar luminosity from about ten for solar-type stars to nearly unity in the discs around the most luminous A-type stars. Here we explore this trend in more detail, checking how significant it is and seeking to find possible explanations. We show that the trend is robust to variation of the composition and porosity of dust particles. For any assumed grain properties and stellar parameters, we suggest a recipe of how to estimate the "true" radius of a spatially unresolved debris disc, based solely on its spectral energy distribution. The results of our collisional simulations are qualitatively consistent with the trend, although additional effects may also be at work. In particular, the lack of grains with small $s_\text{min}/s_\text{blow}$ for lower luminosity stars might be caused by the grain surface energy constraint that should limit the size of the smallest collisional fragments. Also, a better agreement between the data and the collisional simulations is achieved when assuming debris discs of more luminous stars to have higher dynamical excitation than those of less luminous primaries. This would imply that protoplanetary discs of more massive young stars are more efficient in forming big planetesimals or planets that act as stirrers in the debris discs at the subsequent evolutionary stage.

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 precise mechanisms that provide the non-radiative energy for heating the chromosphere and the corona of the Sun and those of other stars constitute an active field of research. By studying stellar chromospheres one aims at identifying the relevant physical processes. Defining the permittable extent of the parameter space can also serve as a template for the Sun-as-a-star. Earlier observations with Herschel and APEX have revealed the temperature minimum of alpha Cen, but these were unable to spatially resolve the binary into individual components. With the data reported here, we aim at remedying this shortcoming. Furthermore, these earlier data were limited to the wavelength region between 100 and 870mu. In the present context, we intend to extend the spectral mapping to longer wavelengths, where the contrast between stellar photospheric and chromospheric emission becomes increasingly evident. ALMA is particularly suited to point sources, such as unresolved stars. ALMA provides the means to achieve our objectives with both its high sensitivity of the collecting area for the detection of weak signals and the high spatial resolving power of its adaptable interferometer for imaging close multiple stars. This is the first detection of main-sequence stars at a wavelength of 3mm. Furthermore, the individual components of the binary alpha CenAB are clearly detected and spatially well resolved at all ALMA wavelengths. The high S/N of these data permit accurate determination of their relative flux ratios. The previously obtained flux ratio of 0.44, which was based on measurements in the optical and at 70mu, is consistent with the present ALMA results, albeit with a large error bar. Given the distinct difference in their cyclic activity, the similarity of their submm SEDs appears surprising.

A significant fraction of main-sequence stars are encircled by dusty debris discs, where the short-lived dust particles are replenished through collisions between planetesimals. Most destructive collisions occur when the orbits of smaller bodies are dynamically stirred up, either by the gravitational effect of locally formed Pluto-sized planetesimals (self-stirring scenario), or via secular perturbation caused by an inner giant planet (planetary stirring). The relative importance of these scenarios in debris systems is unknown. Here we present new Herschel Space Observatory imagery of 11 discs selected from the most massive and extended known debris systems. All discs were found to be extended at far-infrared wavelengths, five of them being resolved for the first time. We evaluated the feasibility of the self-stirring scenario by comparing the measured disc sizes with the predictions of the model calculated for the ages of our targets. We concluded that the self-stirring explanation works for seven discs. However, in four cases, the predicted pace of outward propagation of the stirring front, assuming reasonable initial disc masses, was far too low to explain the radial extent of the cold dust. Therefore, for HD 9672, HD 16743, HD 21997, and HD 95086, another explanation is needed. We performed a similar analysis for {\ss} Pic and HR 8799, reaching the same conclusion. We argue that planetary stirring is a promising possibility to explain the disk properties in these systems. In HR 8799 and HD 95086 we may already know the potential perturber, since their known outer giant planets could be responsible for the stirring process. Our study demonstrates that among the largest and most massive debris discs self-stirring may not be the only active scenario, and potentially planetary stirring is responsible for destructive collisions and debris dust production in a number of systems.

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.

(Abridged) The radii of debris disks and the sizes of their dust grains are tracers of the formation mechanisms and physical processes operating in these systems. We use a sample of 34 debris disks spatially resolved in various Herschel programs to constrain them. While we modeled disks with both warm and cold components, we focus our analysis only on the cold outer disks, i.e. Kuiper-belt analogs. The disk radii derived from the resolved images reveal a large dispersion, but no significant trend with the stellar luminosity, which argues against ice lines as a dominant player in setting the debris disk sizes. Fixing the disk radii to those inferred from the resolved images, we model the spectral energy distributions to determine the dust temperatures and the grain size distributions. While the dust temperature systematically increases towards earlier spectral types, its ratio to the blackbody temperature at the disk radius decreases with the stellar luminosity. This is explained by an increase of typical grain sizes towards more luminous stars. The sizes are compared to the radiation pressure blowout limit $s_\text{blow}$ that is proportional to the stellar luminosity-to-mass ratio and thus also increases towards earlier spectral classes. The grain sizes in the disks of G- to A-stars are inferred to be several times $s_\text{blow}$ at all stellar luminosities, in agreement with collisional models of debris disks. The sizes, measured in the units of $s_\text{blow}$, appear to decrease with the luminosity, which may be suggestive of the disk's stirring level increasing towards earlier-type stars.

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.

Context. Debris discs are thought to be formed through the collisional grinding of planetesimals, and can be considered as the outcome of planet formation. Understanding the properties of gas and dust in debris discs can help us to comprehend the architecture of extrasolar planetary systems. Herschel Space Observatory far-infrared (IR) photometry and spectroscopy have provided a valuable dataset for the study of debris discs gas and dust composition. This paper is part of a series of papers devoted to the study of Herschel PACS observations of young stellar associations. Aims. This work aims at studying the properties of discs in the Beta Pictoris Moving Group (BPMG) through far-IR PACS observations of dust and gas. Methods. We obtained Herschel-PACS far-IR photometric observations at 70, 100 and 160 microns of 19 BPMG members, together with spectroscopic observations of four of them. Spectroscopic observations were centred at 63.18 microns and 157 microns, aiming to detect [OI] and [CII] emission. We incorporated the new far-IR observations in the SED of BPMG members and fitted modified blackbody models to better characterise the dust content. Results. We have detected far-IR excess emission toward nine BPMG members, including the first detection of an IR excess toward HD 29391.The star HD 172555, shows [OI] emission, while HD 181296, shows [CII] emission, expanding the short list of debris discs with a gas detection. No debris disc in BPMG is detected in both [OI] and [CII]. The discs show dust temperatures in the range 55 to 264 K, with low dust masses (6.6*10^{-5} MEarth to 0.2 MEarth) and radii from blackbody models in the range 3 to 82 AU. All the objects with a gas detection are early spectral type stars with a hot dust component.

Main sequence stars, like the Sun, are often found to be orbited by circumstellar material that can be categorized into two groups, planets and debris. The latter is made up of asteroids and comets, as well as the dust and gas derived from them, which makes debris disks observable in thermal emission or scattered light. These disks may persist over Gyrs through steady-state evolution and/or may also experience sporadic stirring and major collisional breakups, rendering them atypically bright for brief periods of time. Most interestingly, they provide direct evidence that the physical processes (whatever they may be) that act to build large oligarchs from micron-sized dust grains in protoplanetary disks have been successful in a given system, at least to the extent of building up a significant planetesimal population comparable to that seen in the Solar System's asteroid and Kuiper belts. Such systems are prime candidates to host even larger planetary bodies as well. The recent growth in interest in debris disks has been driven by observational work that has provided statistics, resolved images, detection of gas in debris disks, and discoveries of new classes of objects. The interpretation of this vast and expanding dataset has necessitated significant advances in debris disk theory, notably in the physics of dust produced in collisional cascades and in the interaction of debris with planets. Application of this theory has led to the realization that such observations provide a powerful diagnostic that can be used not only to refine our understanding of debris disk physics, but also to challenge our understanding of how planetary systems form and evolve.

[Abridged] Debris discs around main-sequence stars indicate the presence of larger rocky bodies. The components of the nearby binary aCentauri have higher than solar metallicities, which is thought to promote giant planet formation. We aim to determine the level of emission from debris in the aCen system. Having already detected the temperature minimum, Tmin, of aCenA, we here attempt to do so also for the companion aCenB. Using the aCen stars as templates, we study possible effects Tmin may have on the detectability of unresolved dust discs around other stars. We use Herschel and APEX photometry to determine the stellar spectral energy distributions. In addition, we use APEX for spectral line mapping to study the complex background around aCen seen in the photometric images. Models of stellar atmospheres and discs are used to estimate the amount of debris around these stars. For solar-type stars, a fractional dust luminosity fd 2e-7 could account for SEDs that do not exhibit the Tmin-effect. Slight excesses at the 2.5 sigma level are observed at 24 mu for both stars, which, if interpreted to be due to dust, would correspond to fd (1-3)e-5. Dynamical disc modelling leads to rough mass estimates of the putative Zodi belts around the aCen stars, viz. <~4e-6 MMoon of 4 to 1000 mu size grains, distributed according to n a^-3.5. Similarly, for filled-in Tmin emission, corresponding EKBs could account for ~1e-3 MMoon of dust. Light scattered and/or thermally emitted by exo-Zodi discs will have profound implications for future spectroscopic missions designed to search for biomarkers in the atmospheres of Earth-like planets. The F-IR SED of aCenB is marginally consistent with the presence of a minimum temperature region in the upper atmosphere. We also show that an aCenA-like temperature minimum may result in an erroneous apprehension about the presence of dust around other stars.

[abridged] Aims. Our Herschel Open Time Key Programme DUNES aims at detecting and characterizing debris disks around nearby, sun-like stars. In addition to the statistical analysis of the data, the detailed study of single objects through spatially resolving the disk and detailed modeling of the data is a main goal of the project. Methods. We obtained the first observations spatially resolving the debris disk around the sun-like star HIP 17439 (HD23484) using the instruments PACS and SPIRE on board the Herschel Space Observatory. Simultaneous multi-wavelength modeling of these data together with ancillary data from the literature is presented. Results. A standard single component disk model fails to reproduce the major axis radial profiles at 70 um, 100 um, and 160 um simultaneously. Moreover, the best-fit parameters derived from such a model suggest a very broad disk extending from few au up to few hundreds of au from the star with a nearly constant surface density which seems physically unlikely. However, the constraints from both the data and our limited theoretical investigation are not strong enough to completely rule out this model. An alternative, more plausible, and better fitting model of the system consists of two rings of dust at approx. 30 au and 90 au, respectively, while the constraints on the parameters of this model are weak due to its complexity and intrinsic degeneracies. Conclusions. The disk is probably composed of at least two components with different spatial locations (but not necessarily detached), while a single, broad disk is possible, but less likely. The two spatially well-separated rings of dust in our best-fit model suggest the presence of at least one high mass planet or several low-mass planets clearing the region between the two rings from planetesimals and dust.

(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.

GASPS is a far-infrared line and continuum survey of protoplanetary and young debris disks using PACS on the Herschel Space Observatory. The survey includes [OI] at 63 microns, as well as 70, 100 and 160um continuum, with the brightest objects also studied in [OI]145um, [CII]157um, H2O and CO. Targets included T Tauri stars and debris disks in 7 nearby young associations, and a sample of isolated Herbig AeBe stars. The aim was to study the global gas and dust content in a wide disk sample, systemically comparing the results with models. In this paper we review the main aims, target selection and observing strategy. We show initial results, including line identifications, sources detected, and a first statistical study. [OI]63um was the brightest line in most objects, by a factor of ~10. Detection rates were 49%, including 100% of HAeBe stars and 43% of T Tauri stars. Comparison with published dust masses show a dust threshold for [OI]63um detection of ~1e-5 M_solar. Normalising to 140pc distance, 32% with mass 1e-6 - 1e-5 M_solar, and a small number with lower mass were also detected. This is consistent with moderate UV excess and disk flaring. In most cases, continuum and line emission is spatially and spectrally unresolved, suggesting disk emission. ~10 objects were resolved, likely from outflows. Detection rates in [OI]145um, [CII]157um and CO J=18-17 were 20-40%, but [CII] was not correlated with disk mass, suggesting it arises instead from a compact envelope. [OI] detection rates in T Tauri associations of ages 0.3-4Myr were ~50%. ~2 stars were detectable in associations of 5-20Myr, with no detections in associations of age >20Myr. Comparing with the total number of young stars, and assuming a ISM-like gas/dust ratio, this indicates that ~18% of stars retain a gas-rich disk of total mass >1M_Jupiter for 1-4Myr, 1-7% keep such disks for 5-10Myr, and none remain beyond 10-20Myr.

HD32297 is a young A-star (~30 Myr) 112 pc away with a bright edge-on debris disk that has been resolved in scattered light. We observed the HD32297 debris disk in the far-infrared and sub-millimeter with the Herschel Space Observatory PACS and SPIRE instruments, populating the spectral energy distribution (SED) from 63 to 500{\mu}m. We aimed to determine the composition of dust grains in the HD32297 disk through SED modeling, using geometrical constraints from the resolved imaging to break degeneracies inherent in SED modeling. We found the best fitting SED model has 2 components: an outer ring centered around 110 AU, seen in the scattered light images, and an inner disk near the habitable zone of the star. The outer disk appears to be composed of grains > 2{\mu}m consisting of silicates, carbonaceous material, and water ice with an abundance ratio of 1:2:3 respectively and 90% porosity. These grains appear consistent with cometary grains, implying the underlying planetesimal population is dominated by comet-like bodies. We also discuss the 3.7{\sigma} detection of [C II] emission at 158{\mu}m with the Herschel PACS Spectrometer, making HD32297 one of only a handful of debris disks with circumstellar gas detected.

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.

Chromospheres and coronae are common phenomena on solar-type stars. Understanding the energy transfer to these heated atmospheric layers requires direct access to the relevant empirical data. Study of these structures has, by and large, been limited to the Sun thus far. The region of the temperature reversal can be directly observed only in the far infrared and submm. We aim at the determination of the characteristics of the atmosphere in the region of the temperature minimum of the solar sister star alpha Cen A. For the nearby binary system alpha Centauri, stellar parameters are known with high accuracy from measurements. For the basic model parameters Teff, log g and [Fe/H], we interpolate in the grid of GAIA/PHOENIX stellar model atmospheres and compute the corresponding model for the G2 V star alpha Cen A. Comparison with photometric measurements shows excellent agreement between observed photospheric data in the optical and infrared. For longer wavelengths, the modelled spectral energy distribution is compared to MIPS, PACS, SPIRE and LABOCA photometry. A specifically tailored Uppsala model based on the MARCS code and extending further in wavelength is used to gauge the emission characteristics of alpha Cen A in the FIR. Similar to the Sun, the FIR emission of alpha Cen A originates in the minimum temperature region above the stellar photosphere in the visible. However, in comparison with the solar case, the FIR photosphere of alpha Cen A appears marginally cooler, Tmin=T160mu=3920+/-375 K. Beyond the minimum near 160mu, the brightness temperatures increase and this radiation likely originates in warmer regions of the chromosphere of alpha Cen A. To the best of our knowledge this is the first time a temperature minimum has been directly measured on a main-sequence star other than the Sun.

We present Herschel PACS photometry of seventeen B- to M-type stars in the 30 Myr-old Tucana-Horologium Association. This work is part of the Herschel Open Time Key Programme "Gas in Protoplanetary Systems" (GASPS). Six of the seventeen targets were found to have infrared excesses significantly greater than the expected stellar IR fluxes, including a previously unknown disk around HD30051. These six debris disks were fitted with single-temperature blackbody models to estimate the temperatures and abundances of the dust in the systems. For the five stars that show excess emission in the Herschel PACS photometry and also have Spitzer IRS spectra, we fit the data with models of optically thin debris disks with realistic grain properties in order to better estimate the disk parameters. The model is determined by a set of six parameters: surface density index, grain size distribution index, minimum and maximum grain sizes, and the inner and outer radii of the disk. The best fitting parameters give us constraints on the geometry of the dust in these systems, as well as lower limits to the total dust masses. The HD105 disk was further constrained by fitting marginally resolved PACS 70 micron imaging.

Aims. We present photometric data of debris disks around HIP 103389 (HD 199260), HIP 107350 (HN Peg, HD206860), and HIP 114948 (HD 219482), obtained in the context of our Herschel Open Time Key Program DUNES (DUst around NEarby Stars). Methods. We used Herschel/PACS to detect the thermal emission of the three debris disks with a 3 sigma sensitivity of a few mJy at 100 um and 160 um. In addition, we obtained Herschel/PACS photometric data at 70 um for HIP 103389. Two different approaches are applied to reduce the Herschel data to investigate the impact of data reduction on the photometry. We fit analytical models to the available spectral energy distribution (SED) data. Results. The SEDs of the three disks potentially exhibit an unusually steep decrease at wavelengths > 70 um. We investigate the significance of the peculiar shape of these SEDs and the impact on models of the disks provided it is real. Our modeling reveals that such a steep decrease of the SEDs in the long wavelength regime is inconsistent with a power-law exponent of the grain size distribution -3.5 expected from a standard equilibrium collisional cascade. In contrast, a very distinct range of grain sizes is implied to dominate the thermal emission of such disks. However, we demonstrate that the understanding of the data of faint sources obtained with Herschel is still incomplete and that the significance of our results depends on the version of the data reduction pipeline used. Conclusions. A new mechanism to produce the dust in the presented debris disks, deviations from the conditions required for a standard equilibrium collisional cascade (grain size exponent of -3.5), and/or significantly different dust properties would be necessary to explain the potentially steep SED shape of the three debris disks presented. (abridged)

(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.

We present Herschel PACS 100 and 160 micron observations of the solar-type stars alpha Men, HD 88230 and HD 210277, which form part of the FGK stars sample of the Herschel Open Time Key Programme (OTKP) DUNES (DUst around NEarby S tars). Our observations show small infrared excesses at 160 micron for all three stars. HD 210277 also shows a small excess at 100 micron, while the 100 micron fluxes of alpha Men and HD 88230 agree with the stellar photospheric predictions. We attribute these infrared excesses to a new class of cold, faint debris discs. alpha Men and HD 88230 are spatially resolved in the PACS 160 micron images, while HD 210277 is point-like at that wavelength. The projected linear sizes of the extended emission lie in the range from ~ 115 to ~ 250 AU. The estimated black body temperatures from the 100 and 160 micron fluxes are $\lesssim$ 22 K, while the fractional luminosity of the cold dust is Ldust/Lstar ~ 10E-6, close to the luminosity of the Solar-System's Kuiper belt. These debris discs are the coldest and faintest discs discovered so far around mature stars and cannot easily be explained by invoking "classical" debris disc models.

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.

Massive planetary cores ($\sim 10$ Earth masses) trigger rapid gas accretion to form gas giant planets \rev{such as} Jupiter and Saturn. We investigate the core growth and the possibilities for cores to reach such a critical core mass. At the late stage, planetary cores grow through collisions with small planetesimals. Collisional fragmentation of planetesimals, which is induced by gravitational interaction with planetary cores, reduces the amount of planetesimals surrounding them, and thus the final core masses. Starting from small planetesimals that the fragmentation rapidly removes, less massive cores are formed. However, planetary cores acquire atmospheres that enlarge their collisional cross section before rapid gas accretion. Once planetary cores exceed about Mars mass, atmospheres significantly accelerate the growth of cores. We show that, taking into account the effects of fragmentation and atmosphere, initially large planetesimals enable formation of sufficiently massive cores. On the other hand, because the growth of cores is slow for large planetesimals, a massive disk is necessary for cores to grow enough within a disk lifetime. If the disk with 100\,km-sized initial planetesimals is 10 times as massive as the minimum mass solar nebula, planetary cores can exceed 10 Earth masses in the Jovian planet region ($>5\,$AU).

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.

Context: The K2V star eps Eri hosts one known inner planet, an outer Kuiper belt analog, and an inner disk of warm dust. Spitzer/IRS measurements indicate that the warm dust is present at distances as close as a few AU from the star. Its origin is puzzling, since an "asteroid belt" that could produce this dust would be unstable because of the known inner planet. Aims: Here we test the hypothesis that the observed warm dust is generated by collisions in the outer belt and is transported inward by Poynting-Robertson (P-R) drag and strong stellar winds. Methods: We simulated a steady-state distribution of dust particles outside 10AU with a collisional code and in the inner region (r<10AU) with single-particle numerical integrations. By assuming homogeneous spherical dust grains composed of water ice and silicate, we calculated the thermal emission of the dust and compared it with observations. We investigated two different orbital configurations for the inner planet inferred from RV measurements, one with a highly eccentric orbit of e=0.7 and another one with a moderate one of e=0.25. We also produced a simulation without a planet. Results: Our models can reproduce the shape and magnitude of the observed SED from mid-IR to sub-mm wavelengths, as well as the Spitzer/MIPS radial brightness profiles. The best-fit dust composition includes both ice and silicates. The results are similar for the two possible planetary orbits and without a planet. Conclusions: The observed warm dust in the system can indeed stem from the outer belt and be transported inward by P-R and stellar wind drag. The inner planet has little effect on the distribution of dust, so that the planetary orbit could not be constrained. Reasonable agreement between the model and observations can only be achieved by relaxing the assumption of purely silicate dust and assuming a mixture of silicate and ice in comparable amounts.

(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.

In an effort to simultaneously study the gas and dust components of the disc surrounding the young Herbig Ae star HD 169142, we present far-IR observations obtained with the PACS instrument onboard the Herschel Space Observatory. This work is part of the Open Time Key Project GASPS, which is aimed at studying the evolution of protoplanetary discs. To constrain the gas properties in the outer disc, we observed the star at several key gas-lines, including [OI] 63.2 and 145.5 micron, [CII] 157.7 micron, CO 72.8 and 90.2 micron, and o-H2O 78.7 and 179.5 micron. We only detect the [OI] 63.2 micron line in our spectra, and derive upper limits for the other lines. We complement our data set with PACS photometry and 12/13CO data obtained with the Submillimeter Array. Furthermore, we derive accurate stellar parameters from optical spectra and UV to mm photometry. We model the dust continuum with the 3D radiative transfer code MCFOST and use this model as an input to analyse the gas lines with the thermo-chemical code ProDiMo. Our dataset is consistent with a simple model in which the gas and dust are well-mixed in a disc with a continuous structure between 20 and 200 AU, but this is not a unique solution. Our modelling effort allows us to constrain the gas-to-dust mass ratio as well as the relative abundance of the PAHs in the disc by simultaneously fitting the lines of several species that originate in different regions. Our results are inconsistent with a gas-poor disc with a large UV excess; a gas mass of 5.0 +/- 2.0 times 10^(-3) Msun is still present in this disc, in agreement with earlier CO observations.

The Herschel GASPS Key Program is a survey of the gas phase of protoplanetary discs, targeting 240 objects which cover a large range of ages, spectral types, and disc properties. To interpret this large quantity of data and initiate self-consistent analyses of the gas and dust properties of protoplanetary discs, we have combined the capabilities of the radiative transfer code MCFOST with the gas thermal balance and chemistry code ProDiMo to compute a grid of 300 000 disc models (DENT). We present a comparison of the first Herschel/GASPS line and continuum data with the predictions from the DENT grid of models. Our objective is to test some of the main trends already identified in the DENT grid, as well as to define better empirical diagnostics to estimate the total gas mass of protoplanetary discs. Photospheric UV radiation appears to be the dominant gas-heating mechanism for Herbig stars, whereas UV excess and/or X-rays emission dominates for T Tauri stars. The DENT grid reveals the complexity in the analysis of far-IR lines and the difficulty to invert these observations into physical quantities. The combination of Herschel line observations with continuum data and/or with rotational lines in the (sub-)millimetre regime, in particular CO lines, is required for a detailed characterisation of the physical and chemical properties of circumstellar discs.

Planets are formed in disks around young stars. With an age of ~10 Myr, TW Hya is one of the nearest T Tauri stars that is still surrounded by a relatively massive disk. In addition a large number of molecules has been found in the TW Hya disk, making TW Hya the perfect test case in a large survey of disks with Herschel-PACS to directly study their gaseous component. We aim to constrain the gas and dust mass of the circumstellar disk around TW Hya. We observed the fine-structure lines of [OI] and [CII] as part of the Open-time large program GASPS. We complement this with continuum data and ground-based 12CO 3-2 and 13CO 3-2 observations. We simultaneously model the continuum and the line fluxes with the 3D Monte-Carlo code MCFOST and the thermo-chemical code ProDiMo to derive the gas and dust masses. We detect the [OI] line at 63 micron. The other lines that were observed, [OI] at 145 micron and [CII] at 157 micron, are not detected. No extended emission has been found. Preliminary modeling of the photometric and line data assuming [12CO]/[13CO]=69 suggests a dust mass for grains with radius < 1 mm of ~1.9 times 10^-4 Msun (total solid mass of 3 times 10^-3 Msun) and a gas mass of (0.5--5) times 10^-3 Msun. The gas-to-dust mass may be lower than the standard interstellar value of 100.

Context - Circumstellar discs are ubiquitous around young stars, but rapidly dissipate their gas and dust on timescales of a few Myr. The Herschel space observatory allows for the study of the warm disc atmosphere, using far-infrared spectroscopy to measure gas content and excitation conditions, and far-IR photometry to constrain the dust distribution. Aims - We aim to detect and characterize the gas content of circumstellar discs in four targets as part of the Herschel science demonstration phase. Methods - We carried out sensitive medium resolution spectroscopy and high sensitivity photometry at lambda ~60-190 micron using the Photodetector Array Camera and Spectrometer instrument on the Herschel space observatory. Results - We detect [OI] 63 micron emission from the young stars HD 169142, TW Hydrae, and RECX 15, but not HD 181327. No other lines, including [CII] 158 and [OI] 145, are significantly detected. All four stars are detected in photometry at 70 and 160 micron. Extensive models are presented in associated papers.

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}

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.

As planetary embryos grow, gravitational stirring of planetesimals by embryos strongly enhances random velocities of planetesimals and makes collisions between planetesimals destructive. The resulting fragments are ground down by successive collisions. Eventually the smallest fragments are removed by the inward drift due to gas drag. Therefore, the collisional disruption depletes the planetesimal disk and inhibits embryo growth. We provide analytical formulae for the final masses of planetary embryos, taking into account planetesimal depletion due to collisional disruption. Furthermore, we perform the statistical simulations for embryo growth (which excellently reproduce results of direct $N$-body simulations if disruption is neglected). These analytical formulae are consistent with the outcome of our statistical simulations. Our results indicate that the final embryo mass at several AU in the minimum-mass solar nebula can reach about $\sim 0.1$ Earth mass within $10^7$ years. This brings another difficulty in formation of gas giant planets, which requires cores with $\sim 10$ Earth masses for gas accretion. However, if the nebular disk is 10 times more massive than the minimum-mass solar nebula and the initial planetesimal size is larger than 100 km, as suggested by some models of planetesimal formation, the final embryo mass reaches about 10 Earth masses at 3-4 AU. The enhancement of embryos' collisional cross sections by their atmosphere could further increase their final mass to form gas giant planets at 5-10 AU in the solar system.

Debris disks are optically thin, almost gas-free dusty disks observed around a significant fraction of main-sequence stars older than about 10 Myr. Since the circumstellar dust is short-lived, the very existence of these disks is considered as evidence that dust-producing planetesimals are still present in mature systems, in which planets have formed - or failed to form - a long time ago. It is inferred that these planetesimals orbit their host stars at asteroid to Kuiper-belt distances and continually supply fresh dust through mutual collisions. This review outlines observational techniques and results on debris disks, summarizes their essential physics and theoretical models, and then places them into the general context of planetary systems, uncovering interrelations between the disks, dust parent bodies, and planets. It is shown that debris disks can serve as tracers of planetesimals and planets and shed light on the planetesimal and planet formation processes that operated in these systems in the past.

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.

Disks around young stars are known to evolve from optically thick, gas-dominated protoplanetary disks to optically thin, almost gas-free debris disks. It is thought that the primordial gas is largely removed at ages of ~10 Myr, but it is difficult to discern the true gas densities from gas observations. This suggests using observations of dust: it has been argued that gas, if present with higher densities, would lead to flatter radial profiles of the dust density and surface brightness than those actually observed. However, here we show that these profiles are surprisingly insensitive to variation of the parameters of a central star, location of the dust-producing planetesimal belt, dustiness of the disk and - most importantly - the parameters of the ambient gas. This result holds for a wide range of gas densities (three orders of magnitude), for different radial distributions of the gas temperature, and different gas compositions. The brightness profile slopes of -3...-4 we find are the same that were theoretically found for gas-free debris disks, and they are the same as actually retrieved from observations of many debris disks. Our specific results for three young (10-30 Myr old), spatially resolved, edge-on debris disks (beta Pic, HD 32297, and AU Mic) show that the observed radial profiles of the surface brightness do not pose any stringent constraints on the gas component of the disk. We cannot exclude that outer parts of the systems may have retained substantial amounts of primordial gas which is not evident in the gas observations (e.g. as much as 50 Earth masses for beta Pic). However, the possibility that gas, most likely secondary, is only present in little to moderate amounts, as deduced from gas detections (e.g. ~0.05 Earth masses in the beta Pic disk), remains open, too.

Analyses of the cosmic microwave background (CMB) radiation maps made by the Wilkinson Microwave Anisotropy Probe (WMAP) have revealed anomalies not predicted by the standard inflationary cosmology. In particular, the power of the quadrupole moment of the CMB fluctuations is remarkably low, and the quadrupole and octopole moments are aligned mutually and with the geometry of the Solar system. It has been suggested in the literature that microwave sky pollution by an unidentified dust cloud in the vicinity of the Solar system may be the cause for these anomalies. In this paper, we simulate the thermal emission by clouds of spherical homogeneous particles of several materials. Spectral constraints from the WMAP multi-wavelength data and earlier infrared observations on the hypothetical dust cloud are used to determine the dust cloud's physical characteristics. In order for its emissivity to demonstrate a flat, CMB-like wavelength dependence over the WMAP wavelengths (3 through 14 mm), and to be invisible in the infrared light, its particles must be macroscopic. Silicate spheres from several millimetres in size and carbonaceous particles an order of magnitude smaller will suffice. According to our estimates of the abundance of such particles in the Zodiacal cloud and trans-neptunian belt, yielding the optical depths of the order of 1E-7 for each cloud, the Solar-system dust can well contribute 10 microKelvin (within an order of magnitude) in the microwaves. This is not only intriguingly close to the magnitude of the anomalies (about 30 microKelvin), but also alarmingly above the presently believed magnitude of systematic biases of the WMAP results (below 5 microKelvin) and, to an even greater degree, of the future missions with higher sensitivities, e.g. PLANCK.

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.

Although known for almost a century, the photophoretic force has only recently been considered in astrophysical context for the first time. In our work, we have examined the effect of photophoresis, acting together with stellar gravity, radiation pressure, and gas drag, on the evolution of solids in transitional circumstellar disks. We have applied our calculations to four different systems: the disks of HR 4796A and HD 141569A, which are several Myr old AB-type stars, and two hypothetical systems that correspond to the solar nebula after disk dispersal has progressed sufficiently for the disk to become optically thin. Our results suggest that solid objects migrate inward or outward, until they reach a certain size-dependent stability distance from the star. The larger the bodies, the closer to the star they tend to accumulate. Photophoresis increases the stability radii, moving objects to larger distances. What is more, photophoresis may cause formation of a belt of objects, but only in a certain range of sizes and only around low-luminosity stars. The effects of photophoresis are noticeable in the size range from several micrometers to several centimeters (for older transitional disks) or even several meters (for younger, more gaseous, ones). We argue that due to gas damping, rotation does not substantially inhibit photophoresis.

We consider a belt of small bodies around a star, captured in one of the external or 1:1 mean-motion resonances with a massive perturber. The objects in the belt collide with each other. Combining methods of celestial mechanics and statistical physics, we calculate mean collisional velocities and collisional rates, averaged over the belt. The results are compared to collisional velocities and rates in a similar, but non-resonant belt, as predicted by the particle-in-a-box method. It is found that the effect of the resonant lock on the velocities is rather small, while on the rates more substantial. The collisional rates between objects in an external resonance are by about a factor of two higher than those in a similar belt of objects not locked in a resonance. For Trojans under the same conditions, the collisional rates may be enhanced by up to an order of magnitude. Our results imply, in particular, shorter collisional lifetimes of resonant Kuiper belt objects in the solar system and higher efficiency of dust production by resonant planetesimals in debris disks around other stars.

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.

Tenuous dust clouds of Jupiter's Galilean moons Io, Europa, Ganymede and Callisto have been detected with the in-situ dust detector on board the Galileo spacecraft. The majority of the dust particles have been sensed at altitudes below five radii of these lunar-sized satellites. We identify the particles in the dust clouds surrounding the moons by their impact direction, impact velocity, and mass distribution. Average particle sizes are 0.5 to $\rm 1 \mu m$, just above the detector threshold, indicating a size distribution with decreasing numbers towards bigger particles. Our results imply that the particles have been kicked up by hypervelocity impacts of micrometeoroids onto the satellites' surfaces. The measured radial dust density profiles are consistent with predictions by dynamical modeling for satellite ejecta produced by interplanetary impactors (Krivov et al., PSS, 2003, 51, 251--269), assuming yield, mass and velocity distributions of the ejecta from laboratory measurements. The dust clouds of the three outer Galilean moons have very similar properties and are in good agreement with the model predictions for solid ice-silicate surfaces. The dust density in the vicinity of Io, however, is more than an order of magnitude lower than expected from theory. This may be due to a softer, fluffier surface of Io (volcanic deposits) as compared to the other moons. The log-log slope of the dust number density in the clouds vs. distance from the satellite center ranges between --1.6 and --2.8. Appreciable variations of number densities obtained from individual flybys with varying geometry, especially at Callisto, might be indicative of leading-trailing asymmetries of the clouds due to the motion of the moons with respect to the field of impactors.

A dust cloud of Ganymede has been detected by in-situ measurements with the dust detector onboard the Galileo spacecraft. The dust grains have been sensed at altitudes below five Ganymede radii (Ganymede radius = $\rm 2,635 km$). Our analysis identifies the particles in the dust cloud surrounding Ganymede by their impact direction, impact velocity, and mass distribution and implies that they have been kicked up by hypervelocity impacts of micrometeoroids onto the satellite's surface. We calculate the radial density profile of the particles ejected from the satellite by interplanetary dust grains. We assume the yields, mass and velocity distributions of the ejecta obtained from laboratory impact experiments onto icy targets and consider the dynamics of the ejected grains in ballistic and escaping trajectories near Ganymede. The spatial dust density profile calculated with interplanetary particles as impactors is consistent with the profile derived from the Galileo measurements. The contribution of interstellar grains as projectiles is negligible. Dust measurements in the vicinities of satellites by spacecraft detectors are suggested as a beneficial tool to obtain more knowledge about the satellite surfaces, as well as dusty planetary rings maintained by satellites through the impact ejecta mechanism.