The infrared dust emission from the white dwarf GD 56 is found to rise and fall by 20% peak-to-peak over 11.2 yr, and is consistent with ongoing dust production and depletion. It is hypothesized that the dust is produced via collisions associated with an evolving dust disk, temporarily increasing the emitting surface of warm debris, and is subsequently destroyed or assimilated within a few years. The variations are consistent with debris that does not change temperature, indicating that dust is produced and depleted within a fixed range of orbital radii. Gas produced in collisions may rapidly re-condense onto grains, or may accrete onto the white dwarf surface on viscous timescales that are considerably longer than Poynting-Robertson drag for micron-sized dust. This potential delay in mass accretion rate change is consistent with multi-epoch spectra of the unchanging Ca II and Mg II absorption features in GD 56 over 15 yr, although the sampling is sparse. Overall these results indicate that collisions are likely to be the source of dust and gas, either inferred or observed, orbiting most or all polluted white dwarfs.

The most heavily polluted white dwarfs often show excess infrared radiation from circumstellar dust disks, which are modeled as a result of tidal disruption of extrasolar minor planets. Interaction of dust, gas, and disintegrating objects can all contribute to the dynamical evolution of these dust disks. Here, we report on two infrared variable dusty white dwarfs, SDSS J1228+1040 and G29-38. For SDSS J1228+1040, compared to the first measurements in 2007, the IRAC [3.6] and [4.5] fluxes decreased by 20% by 2014 to a level also seen in the recent 2018 observations. For G29-38, the infrared flux of the 10 $\mu$m silicate emission feature became 10% stronger between 2004 and 2007, We explore several scenarios that could account for these changes, including tidal disruption events, perturbation from a companion, and runaway accretion. No satisfactory causes are found for the flux drop in SDSS J1228+1040 due to the limited time coverage. Continuous tidal disruption of small planetesimals could increase the mass of small grains and concurrently change the strength of the 10 $\mu$m feature of G29-38. Dust disks around white dwarfs are actively evolving and we speculate that there could be different mechanisms responsible for the temporal changes of these disks.

Context. Kepler observations have revealed a class of short period exoplanets, of which Kepler-1520 b is the prototype, which have comet-like dust tails thought to be the result of small, rocky planets losing mass. The shape and chromaticity of the transits constrain the properties of the dust particles originating from the planet's surface, offering a unique opportunity to probe the composition and geophysics of rocky exoplanets. Aims. We aim to approximate the average Kepler long-cadence light curve of Kepler-1520 b and investigate how the optical thickness and transit cross-section of a general dust tail can affect the observed wavelength dependence and depth of transit light curves. Methods. We developed a new 3D model that ejects sublimating particles from the planet surface to build up a dust tail, assuming it to be optically thin, and used 3D radiative transfer computations that fully treat scattering using the distribution of hollow spheres (DHS) method, to generate transit light curves between 0.45 and 2.5 $\mu$m. Results. We show that the transit depth is wavelength independent for optically thick tails, potentially explaining why only some observations indicate a wavelength dependence. From the 3D nature of our simulated tails, we show that their transit cross-sections are related to the component of particle ejection velocity perpendicular to the planet's orbital plane and use this to derive a minimum ejection velocity of 1.2 kms$^{-1}$. To fit the average transit depth of Kepler-1520 b of 0.87%, we require a high dust mas-loss rate of 7 $-$ 80 M$_\oplus$ Gyr$^{-1}$ which implies planet lifetimes that may be inconsistent with the observed sample. Therefore, these mass-loss rates should be considered to be upper limits.

The star KIC 8462852 (Boyajian's Star) displays both fast dips of up to 20% on time scales of days, plus long-term secular fading by up to 19% on time scales from a year to a century. We report on CCD photometry of KIC 8462852 from 2015.75 to 2018.18, with 19,176 images making for 1,866 nightly magnitudes in BVRI. Our light curves show a continuing secular decline (by 0.023 +- 0.003 mags in the B-band) with three superposed dips with duration 120-180 days. This demonstrates that there is a continuum of dip durations from a day to a century, so that the secular fading is seen to be by the same physical mechanism as the short-duration Kepler dips. The BVRI light curves all have the same shape, with the slopes and amplitudes for VRI being systematically smaller than in the B-band by factors of 0.77 +- 0.05, 0.50 +- 0.05, and 0.31 +- 0.05. We rule out any hypothesis involving occultation of the primary star by any star, planet, solid body, or optically thick cloud. But these ratios are the same as that expected for ordinary extinction by dust clouds. This chromatic extinction implies dust particle sizes going down to ~0.1 micron, suggesting that this dust will be rapidly blown away by stellar radiation pressure, so the dust clouds must have formed within months. The modern infrared observations were taken at a time when there was at least 12.4% +- 1.3% dust coverage (as part of the secular dimming), and this is consistent with dimming originating in circumstellar dust.

This paper shows how the dips and secular dimming in the KIC8462852 light curve can originate in circumstellar material distributed around a single elliptical orbit (e.g., exocomets). The expected thermal emission and wavelength dependent dimming is derived for different orbital parameters and geometries, including dust that is optically thick to stellar radiation, and for a size distribution of dust with realistic optical properties. We first consider dust distributed evenly around the orbit, then show how to derive its uneven distribution from the optical light curve and to predict light curves at different wavelengths. The fractional luminosity of an even distribution is approximately the level of dimming times stellar radius divided by distance from the star at transit. Non-detection of dust thermal emission for KIC8462852 thus provides a lower limit on the transit distance to complement the 0.6au upper limit imposed by 0.4day dips. Unless the dust distribution is optically thick, the putative 16% century-long secular dimming must have disappeared before the WISE 12micron measurement in 2010, and subsequent 4.5micron observations require transits at >0.05au. However, self-absorption of thermal emission removes these constraints for opaque dust distributions. The passage of dust clumps through pericentre is predicted to cause infrared brightening lasting 10s of days and dimming during transit, such that total flux received decreases at wavelengths <5micron, but increases to potentially detectable levels at longer wavelengths. We suggest that lower dimming levels than seen for KIC8462852 are more common in the Galactic population and may be detected in future transit surveys.

Several tens of white dwarfs are known to host circumstellar discs of dusty debris, thought to arise from the tidal disruption of rocky bodies originating in the star's remnant planetary system. This paper investigates the evolution of such discs if they are very massive, as may be the case if their progenitor was a terrestrial planet, moon, or dwarf planet. Assuming the discs are physically thin and flat, like Saturn's rings, their evolution is governed by Poynting-Robertson drag or viscous spreading, where the disc's effective viscosity is due to self-gravity wakes. For discs with masses >10^26 g, located in the outer parts of the tidal disruption zone, viscous spreading dominates the evolution, and mass is transported both in- and outwards. When outwards-spreading material flows beyond the Roche limit, it coagulates into new (minor) planets in a process analogous to the ongoing formation of moonlets at the outer edge of Saturn's rings. The newly formed bodies migrate outwards by exchanging angular momentum with the disc and coalesce into larger objects through mutual collisions. Eventually, the disc's Roche-limit overflow recycles tens of percent of the original disc mass; most ends up in a single large body near 2:1 mean-motion resonance with the disc's outer edge. Hence, the recycling of a tidally disrupted super-Earth, for example, could yield an Earth-mass planet on a ~10-h orbit, located in the habitable zone for 2-to-10-Gyr-old white dwarfs. The recycling process also creates a population of smaller bodies just outside the Roche limit, which may explain the minor planets recently postulated to orbit WD 1145+017.

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.

White dwarf WD 1145+017 is orbited by several clouds of dust, possibly emanating from actively disintegrating bodies. These dust clouds reveal themselves through deep, broad, and evolving transits in the star's light curve. Here, we report two epochs of multi-wavelength photometric observations of WD 1145+017, including several filters in the optical, K$_\mathrm{s}$ and 4.5 $\mu$m bands in 2016 and 2017. The observed transit depths are different at these wavelengths. However, after correcting for excess dust emission at K$_\mathrm{s}$ and 4.5 $\mu$m, we find the transit depths for the white dwarf itself are the same at all wavelengths, at least to within the observational uncertainties of $\sim$5%-10%. From this surprising result, and under the assumption of low optical depth dust clouds, we conclude that there is a deficit of small particles (with radii $s \lesssim$ 1.5 $\mu$m) in the transiting material. We propose a model wherein only large particles can survive the high equilibrium temperature environment corresponding to 4.5 hr orbital periods around WD 1145+017, while small particles sublimate rapidly. In addition, we evaluate dust models that are permitted by our measurements of infrared emission.

We discuss a new class of exoplanets that appear to be emitting a tail of dusty effluents. These disintegrating planets are found close to their host stars and have very hot, and likely molten, surfaces. The properties of the dust should provide a direct probe of the constituent material of these rocky bodies.

Infrared excesses around metal polluted white dwarfs have been associated with the accretion of dusty, planetary material. This work analyses the available infrared data for an unbiased sample of white dwarfs and demonstrates that no more than 3.3% can have a wide, flat, opaque dust disc, extending to the Roche radius, with a temperature at the disc inner edge of $T_{in} = 1,400$K, the standard model for the observed excesses. This is in stark contrast to the incidence of pollution of about 30%. We present four potential reasons for the absence of an infrared excess in polluted white dwarfs, depending on their stellar properties and inferred accretion rates: i) their dust discs are opaque, but narrow, thus evading detection if more than 85% of polluted white dwarfs have dust discs narrower than $\delta r< 0.04r$, ii) their dust discs have been fully consumed, which only works for the oldest white dwarfs with sinking timescales longer than hundreds of years, iii) their dust is optically thin, which can supply low accretion rates of $<10^7 $gs$^{-1}$ if dominated by PR-drag, and higher accretion rates, if inwards transport of material is enhanced, for example due to the presence of gas, iv) their accretion is supplied by a pure gas disc, which could result from the sublimation of optically thin dust for T* > 20, 000K. Future observations sensitive to faint infrared excesses or the presence of gas, can test the scenarios presented here, thereby better constraining the nature of the material fuelling accretion in polluted white dwarfs.

Evaporating rocky exoplanets, such as KIC 12557548b, eject large amounts of dust grains, which can trail the planet in a comet-like tail. When such objects occult their host star, the resulting transit signal contains information about the dust in the tail. We aim to use the detailed shape of the Kepler light curve of KIC 12557548b to constrain the size and composition of the dust grains that make up the tail, as well as the mass loss rate of the planet. Using a self-consistent numerical model of the dust dynamics and sublimation, we calculate the shape of the tail by following dust grains from their ejection from the planet to their destruction due to sublimation. From this dust cloud shape, we generate synthetic light curves (incorporating the effects of extinction and angle-dependent scattering), which are then compared with the phase-folded Kepler light curve. We explore the free-parameter space thoroughly using a Markov chain Monte Carlo method. Our physics-based model is capable of reproducing the observed light curve in detail. Good fits are found for initial grain sizes between 0.2 and 5.6 micron and dust mass loss rates of 0.6 to 15.6 M_earth/Gyr (2-sigma ranges). We find that only certain combinations of material parameters yield the correct tail length. These constraints are consistent with dust made of corundum (Al2O3), but do not agree with a range of carbonaceous, silicate, or iron compositions. Using a detailed, physically motivated model, it is possible to constrain the composition of the dust in the tails of evaporating rocky exoplanets. This provides a unique opportunity to probe to interior composition of the smallest known exoplanets.

Fomalhaut plays an important role in the study of debris disks and small bodies in other planetary systems. The proximity and luminosity of the star make key features of its debris, like the water ice-line, accessible. Here we present ALMA cycle 1, 870 \mu m (345 GHz) observations targeted at the inner part of the Fomalhaut system with a synthesized beam of 0.45"x0.37" (~3 AU linear resolution at the distance of Fomalhaut) and a rms of 26 \mu Jy/beam. The high angular resolution and sensitivity of the ALMA data enable us to place strong constraints on the nature of the warm excess revealed by Spitzer and Herschel observations. We detect a point source at the star position with a total flux consistent with thermal emission from the stellar photosphere. No structures that are brighter than 3\sigma\ are detected in the central 15 AU x 15 AU region. Modeling the spectral energy distribution using parameters expected for a dust-producing planetesimal belt indicates a radial location in the range ~8-15 AU. This is consistent with the location where ice sublimates in Fomalhaut, i.e., an asteroid-belt analog. The 3\sigma\ upper limit for such a belt is <1.3 mJy at 870 \mu m. We also interpret the 2 and 8-13 \mu m interferometric measurements to reveal the structure in the inner 10 AU region as dust naturally connected to this proposed asteroid belt by Poynting-Robertson drag, dust sublimation, and magnetically trapped nano grains.

Recently, two exoplanet candidates have been discovered, KIC 12557548b and KOI-2700b, whose transit profiles show evidence for a comet-like tail of dust trailing the planet, thought to be fed by the evaporation of the planet's surface. We aim to put constraints on the composition of the dust ejected by these objects from the shape of their transit light curves. We derive a semi-analytical expression for the attenuation of dust cross-section in the tail, incorporating the sublimation of dust grains as well as their drift away from the planet. This expression shows that the length of the tail is highly sensitive to the sublimation properties of the dust material. We compute tail lengths for several possible dust compositions, and compare these to observational estimates of the tail lengths of KIC 12557548b and KOI-2700b, inferred from their light curves. The observed tail lengths are consistent with dust grains composed of corundum (Al2O3) or iron-rich silicate minerals (e.g., fayalite, Fe2SiO4). Pure iron and carbonaceous compositions are disfavoured. In addition, we estimate dust mass loss rates of 1.7 +/- 0.5 M_earth/Gyr for KIC 12557548b, and > 0.007 M_earth/Gyr (1-sigma lower limit) for KOI-2700b.

Hot exozodiacal dust is thought to be responsible for excess near-infrared (NIR) emission emanating from the innermost parts of some debris disks. The origin of this dust, however, is still a matter of debate. We test whether hot exozodiacal dust can be supplied from an exterior parent belt by Poynting-Robertson (P-R) drag, paying special attention to the pile-up of dust that occurs due to the interplay of P-R drag and dust sublimation. Specifically, we investigate whether pile-ups still occur when collisions are taken into account, and if they can explain the observed NIR excess. We compute the steady-state distribution of dust in the inner disk by solving the continuity equation. First, we derive an analytical solution under a number of simplifying assumptions. Second, we develop a numerical debris disk model that for the first time treats the complex interaction of collisions, P-R drag, and sublimation in a self-consistent way. From the resulting dust distributions we generate thermal emission spectra and compare these to observed excess NIR fluxes. We confirm that P-R drag always supplies a small amount of dust to the sublimation zone, but find that a fully consistent treatment yields a maximum amount of dust that is about 7 times lower than that given by analytical estimates. The NIR excess due this material is much smaller (<10^-3 for A-type stars with parent belts at >1 AU) than the values derived from interferometric observations (~10^-2). Pile-up of dust still occurs when collisions are considered, but its effect on the NIR flux is insignificant. Finally, the cross-section in the innermost regions is clearly dominated by barely bound grains.

[Abridged] Debris disks are extrasolar analogs to the solar system planetesimal belts. The star Fomalhaut harbors a cold debris belt at 140 AU as well as evidence of a warm dust component, which is suspected of being a bright analog to the solar system's zodiacal dust. Interferometric observations obtained with the VLTI and the KIN have identified near- and mid-infrared excesses attributed to hot and warm exozodiacal dust in the inner few AU of the star. We performed parametric modeling of the exozodiacal disk using the GRaTeR radiative transfer code to reproduce the interferometric data, complemented by mid- to far-infrared measurements. A detailed treatment of sublimation temperatures was introduced to explore the hot population at the sublimation rim. We then used an analytical approach to successively testing several source mechanisms. A good fit to the data is found by two distinct dust populations: (1) very small, hence unbound, hot dust grains confined in a narrow region at the sublimation rim of carbonaceous material; (2) bound grains at 2 AU that are protected from sublimation and have a higher mass despite their fainter flux level. We propose that the hot dust is produced by the release of small carbon grains following the disruption of aggregates that originate from the warm component. A mechanism, such as gas braking, is required to further confine the small grains for a long enough time. In situ dust production could hardly be ensured for the age of the star, so the observed amount of dust must be triggered by intense dynamical activity. Fomalhaut may be representative of exozodis that are currently being surveyed worldwide. We propose a framework for reconciling the hot exozodi phenomenon with theoretical constraints: the hot component of Fomalhaut is likely the tip of the iceberg since it could originate from a warm counterpart residing near the ice line.

The gas transport through non-volatile random porous media is investigated numerically. We extend our previous research of the transport of molecules inside the uppermost layer of a cometary surface (Skorov and Rickmann, 1995; Skorov et al. 2001). We assess the validity of the simplified capillary model and its assumptions to simulate the gas flux trough the porous dust mantle as it has been applied in cometary physics. A new microphysical computational model for molecular transport in random porous media formed by packed spheres is presented. The main transport characteristics such as the mean free path distribution and the permeability are calculated for a wide range of model parameters and compared with those obtained by more idealized models. The focus in this comparison is on limitations inherent in the capillary model. Finally a practical way is suggested to adjust the algebraic Clausing formula taking into consideration the nonlinear dependence of permeability on layer porosity. The retrieved dependence allows us to accurately calculate the permeability of layers whose thickness and porosity vary in the range of values expected for the near-surface regions of a cometary nucleus.

We present high-spectral-resolution, optical spectra of the Herbig Be star MWC 147, in which we spectrally resolve several emission lines, including the [O I] lines at 6300 and 6363\deg. Their highly symmetric, double-peaked line profiles indicate that the emission originates in a rotating circumstellar disk. We deconvolve the Doppler-broadened [O I] emission lines to obtain a measure of emission as a function of distance from the central star. The resulting radial surface brightness profiles are in agreement with a disk structure consisting of a flat, inner, gaseous disk and a flared, outer, dust disk. The transition between these components at 2 to 3 AU corresponds to the estimated dust sublimation radius. The width of the double-peaked Mg II line at 4481\deg suggests that the inner disk extends to at least 0.10 AU, close to the corotation radius.