Multi-epoch infrared photometry from Spitzer is used to monitor circumstellar discs at white dwarfs, which are consistent with disrupted minor planets whose debris is accreted and chemically reflected by their host stars. Widespread infrared variability is found across the population of 37 stars with two or more epochs. Larger flux changes occur on longer time-scales, reaching several tens of per cent over baselines of a few years. The canonical model of a geometrically thin, optically thick disc is thus insufficient, as it cannot give rise to the observed behaviour. Optically thin dust best accounts for the variability, where collisions drive dust production and destruction. Notably, the highest infrared variations are seen in systems that show Ca II emission, supporting planetesimal collisions for all known debris discs, with the most energetic occurring in those with detected gaseous debris. The sample includes the only polluted white dwarf with a circumbinary disc, where the signal of the day-night cycle of its irradiated substellar companion appears diluted by dust emission.

White dwarfs are routinely observed to have polluted atmospheres, and sometimes significant infrared excesses, that indicate ongoing accretion of circumstellar dust and rocky debris. Typically this debris is assumed to be in the form of a (circular) disc, and to originate from asteroids that passed close enough to the white dwarf to be pulled apart by tides. However, theoretical considerations suggest that the circularisation of the debris, which initially occupies highly eccentric orbits, is very slow. We therefore hypothesise that the observations may be readily explained by the debris remaining on highly eccentric orbits, and we explore the properties of such debris. For the generic case of an asteroid originating at several au from the white dwarf, we find that all of the tidal debris is always bound to the white dwarf and that the orbital energy distribution of the debris is narrow enough that it executes similar elliptical orbits with only a narrow spread. Assuming that the tidal field of the white dwarf is sufficient to minimise the effects of self-gravity and collisions within the debris, we estimate the time over which the debris spreads into a single elliptical ring, and we generate toy spectra and lightcurves from the initial disruption to late times when the debris distribution is essentially time steady. Finally we speculate on the connection between these simple considerations and the observed properties of these systems, and on additional physical processes that may change this simple picture.

The magnetorotational instability (MRI) occurs when a weak magnetic field destabilises a rotating, electrically conducting fluid with inwardly increasing angular velocity. The MRI is essential to astrophysical disk theory where the shear is typically Keplerian. Internal shear layers in stars may also be MRI unstable, and they take a wide range of profiles, including near-critical. We show that the fastest growing modes of an ideal magnetofluid are three-dimensional provided the shear rate, $S$, is near the two-dimensional onset value, $S_c$. For a Keplerian shear, three-dimensional modes are unstable above $S\approx0.10S_c$, and dominate the two-dimensional modes until $S\approx2.05S_{c}$. These three-dimensional modes dominate for shear profiles relevant to stars and at magnetic Prandtl numbers relevant to liquid-metal laboratory experiments. Significant numbers of rapidly growing three-dimensional modes remain well past $2.05S_{c}$. These finding are significant in three ways. First, weakly nonlinear theory suggests that the MRI saturates by pushing the shear rate to its critical value. This can happen for systems, like stars and laboratory experiments, that can rearrange their angular velocity profiles. Second, the non-normal character and large transient growth of MRI modes should be important whenever three-dimensionality exists. Finally, three-dimensional growth suggests direct dynamo action driven from the linear instability.

Nine metal-polluted white dwarfs are observed with medium-resolution optical spectroscopy,where photospheric abundances are determined and interpreted through comparison against solar system objects. An improved method of making such comparisons is presented that overcomes potential weaknesses of prior analyses, with the numerous sources of error considered to highlight the limitations on interpretation. The stars are inferred to be accreting rocky, volatile-poor asteroidal materials with origins in differentiated bodies, in line with the consensus model. The most heavily polluted star in the sample has 14 metals detected, and appears to be accreting material from a rocky planetesimal, whose composition is mantle-like with a small Fe-Ni core component. Some unusual abundances are present: one star is strongly depleted in Ca, while two others show Na abundances elevated above bulk Earth, speculated either to reflect diversity in the formation conditions of the source material, or to be traces of past accretion events. Another star shows clear signs that accretion ceased around 5 Myr ago,causing Mg to dominate the photospheric abundances, as it has the longest diffusion time of the observed elements. Observing such post-accretion systems allows constraints to be placed on models of the accretion process.

Practically all known planet hosts will evolve into white dwarfs, and large parts of their planetary systems will survive this transition - the same is true for the solar system beyond the orbit of Mars. Spectroscopy of white dwarfs accreting planetary debris provides the most accurate insight into the bulk composition of exo-planets. Ground-based spectroscopic surveys of ~260, 000 white dwarfs detected with Gaia will identify >1000 evolved planetary systems, and high-throughput high-resolution space-based ultraviolet spectroscopy is essential to measure in detail their abundances. So far, evidence for two planetesimals orbiting closely around white dwarfs has been obtained, and their study provides important constraints on the composition and internal structure of these bodies. Major photometric and spectroscopic efforts will be necessary to assemble a sample of such close-in planetesimals that is sufficiently large to establish their properties as a population, and to deduce the architectures of the outer planetary systems from where they originated. Mid-infrared spectroscopy of the dusty disks will provide detailed mineralogical information of the debris, which, in combination with the elemental abundances measured from the white dwarf spectroscopy, will enable detailed physical modelling of the chemical, thermodynamic, and physical history of the accreted material. Flexible multi-epoch infrared observations are essential to determine the physical nature, and origin of the variability observed in many of the dusty disks. Finally, the direct detection of the outer reservoirs feeding material to the white dwarfs will require sensitive mid- and far-infrared capabilities.

This paper presents combined ${\it Spitzer}$ IRAC and ${\it Hubble}$ COS results for a double-blind survey of 195 single and 22 wide binary white dwarfs for infrared excesses and atmospheric metals. The selection criteria include cooling ages in the range 9 to 300 Myr, and hydrogen-rich atmospheres so that the presence of atmospheric metals can be confidently linked to ongoing accretion from a circumstellar disc. The entire sample has infrared photometry, whereas 168 targets have corresponding ultraviolet spectra. Three stars with infrared excesses due to debris discs are recovered, yielding a nominal frequency of $1.5_{-0.5}^{+1.5}$ per cent, while in stark contrast, the fraction of stars with atmospheric metals is $45\pm4$ per cent. Thus, only one out of 30 polluted white dwarfs exhibits an infrared excess at 3-4 $\mu$m in IRAC photometry, which reinforces the fact that atmospheric metal pollution is the most sensitive tracer of white dwarf planetary systems. The corresponding fraction of infrared excesses around white dwarfs with wide binary companions is consistent with zero, using both the infrared survey data and an independent assessment of potential binarity for well-established dusty and polluted stars. In contrast, the frequency of atmospheric pollution among the targets in wide binaries is indistinct from apparently single stars, and moreover the multiplicity of polluted white dwarfs in a complete and volume-limited sample is the same as for field stars. Therefore, it appears that the delivery of planetesimal material onto white dwarfs is ultimately not driven by stellar companions, but by the dynamics of planetary bodies.

Archival data from the WISE satellite reveals infrared flux variations of tens of per cent around numerous dusty white dwarfs. Data spanning more than seven years reveal more than half of known systems are varying in the 3.4 micron band, while the 4.6 micron data are challenging to interpret due to lower signal-to-noise. The sparsely-sampled data limit interpretation, but the heterogeneous light curves suggest each source may be idiosyncratic, where there maybe competing processes operating on different time-scales. Collisions are likely driving the observed decays in flux, and this finding suggests that dust production is operating more often than indicated by previous observations. The observed variation is at odds with the canonical flat disc model in isolation, and underscores the need for infrared monitoring of these evolved planetary systems to inform the next generation of theoretical models.

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.