Debris disks are exoplanetary systems containing planets, minor bodies (such as asteroids and comets) and debris dust. Unseen planets are presumed to perturb the minor bodies into crossing orbits, generating small dust grains that are detected via remote sensing. Debris disks have been discovered around main sequence stars of a variety of ages (from 10 Myr to several Gyr) and stellar spectral types (from early A-type to M-type stars). As a result, they serve as excellent laboratories for understanding whether the architecture and the evolution of our Solar System is common or rare. This white paper addresses two outstanding questions in debris disk science: (1) Are debris disk minor bodies similar to asteroids and comets in our Solar System? (2) Do planets separate circumstellar material into distinct reservoirs and/or mix material during planet migration? We anticipate that SOFIA/HIRMES, JWST, and WFIRST/CGI will greatly improve our understanding of debris disk composition, enabling the astronomical community to answer these questions. However, we note that despite their observational power, these facilities will not provide large numbers of detections or detailed characterization of cold ices and silicates in the Trans Neptunian zone. Origins Space Telescope is needed to revolutionize our understanding of the bulk composition and mixing in exoplanetary systems.

Thousands of exoplanets have been found with many widely different from the ones in our own system. Despite the success, systems with planets in wide orbits analogous to those of Jupiter and Saturn, in the critical first several hundred million years of evolution, are virtually unexplored. Where are the low-mass planets that are hidden from our exoplanet detection techniques? Is our Solar System's planetary architecture unique? High-fidelity debris disk images offer an effective method to answer these questions. We can use them to study the formation and evolution of low-mass planets from youth to the age of the Solar System, providing snapshots of the complex processes and valuable insights into the formation and migration history of giant planets at wide orbits. This white paper focuses on resolving debris structures in thermal emission that is applicable to a large unbiased sample. We summarize the properties of the known debris disks and assess the feasibility of resolving them within our current and future infrared and millimeter facilities by adopting uniform criteria. JWST and the 9-m Origins Space Telescope are the most promising missions in the coming decades to resolve almost half of the known disks at high fidelity. Resolved debris structures at multiple wavelengths and at all stages of evolution would reveal the properties of unseen planet populations, enabling a unique demographic study of overall planet formation and evolution.

Measuring the masses of protoplanetary disks is crucial for understanding their planet-forming potential. Typically, dust masses are derived from (sub-)millimeter flux density measurements plus assumptions for the opacity, temperature, and optical depth of the dust. Here we use radiative transfer models to quantify the validity of these assumptions with the aim of improving the accuracy of disk dust mass measurements. We first carry out a controlled exploration of disk parameter space. We find that the disk temperature is a strong function of disk size, while the optical depth depends on both disk size and dust mass. The millimeter-wavelength spectral index can be significantly shallower than the naive expectation due to a combination of optical depth and deviations from the Rayleigh-Jeans regime. We fit radiative transfer models to the spectral energy distributions (SEDs) of 132 disks in the Taurus-Auriga region using a Markov chain Monte Carlo approach. We used all available data to produce the most complete SEDs used in any extant modeling study. We perform the fitting twice: first with unconstrained disk sizes and again imposing the disk size--brightness relation inferred for sources in Taurus. This constraint generally forces the disks to be smaller, warmer, and more optically thick. From both sets of fits, we find disks to be $\sim$1--5 times more massive than when derived using (sub-)millimeter measurements and common assumptions. With the uncertainties derived from our model fitting, the previously measured dust mass--stellar mass correlation is present in our study but only significant at the 2$\sigma$ level.

We present ALMA 850 $\mu$m continuum observations of the Orion Nebula Cluster that provide the highest angular resolution ($\sim 0\rlap{.}''1 \approx 40$ AU) and deepest sensitivity ($\sim 0.1$ mJy) of the region to date. We mosaicked a field containing $\sim 225$ optical or near-IR-identified young stars, $\sim 60$ of which are also optically-identified "proplyds". We detect continuum emission at 850 $\mu$m towards $\sim 80$% of the proplyd sample, and $\sim 50$% of the larger sample of previously-identified cluster members. Detected objects have fluxes of $\sim 0.5$-80 mJy. We remove sub-mm flux due to free-free emission in some objects, leaving a sample of sources detected in dust emission. Under standard assumptions of isothermal, optically thin disks, sub-mm fluxes correspond to dust masses of $\sim 0.5$ to 80 Earth masses. We measure the distribution of disk sizes, and find that disks in this region are particularly compact. Such compact disks are likely to be significantly optically thick. The distributions of sub-mm flux and inferred disk size indicate smaller, lower-flux disks than in lower-density star-forming regions of similar age. Measured disk flux is correlated weakly with stellar mass, contrary to studies in other star forming regions that found steeper correlations. We find a correlation between disk flux and distance from the massive star $\theta^1$ Ori C, suggesting that disk properties in this region are influenced strongly by the rich cluster environment.

The architectures of debris disks encode the history of planet formation in these systems. Studies of debris disks via their spectral energy distributions (SEDs) have found infrared excesses arising from cold dust, warm dust, or a combination of the two. The cold outer belts of many systems have been imaged, facilitating their study in great detail. Far less is known about the warm components, including the origin of the dust. The regularity of the disk temperatures indicates an underlying structure that may be linked to the water snow line. If the dust is generated from collisions in an exo-asteroid belt, the dust will likely trace the location of the water snow line in the primordial protoplanetary disk where planetesimal growth was enhanced. If instead the warm dust arises from the inward transport from a reservoir of icy material farther out in the system, the dust location is expected to be set by the current snow line. We analyze the SEDs of a large sample of debris disks with warm components. We find that warm components in single-component systems (those without detectable cold components) follow the primordial snow line rather than the current snow line, so they likely arise from exo-asteroid belts. While the locations of many warm components in two-component systems are also consistent with the primordial snow line, there is more diversity among these systems, suggesting additional effects play a role.

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.

We investigate whether varying the dust composition (described by the optical constants) can solve a persistent problem in debris disk modeling--the inability to fit the thermal emission without over-predicting the scattered light. We model five images of the beta Pictoris disk: two in scattered light from HST/STIS at 0.58 microns and HST/WFC3 at 1.16 microns, and three in thermal emission from Spitzer/MIPS at 24 microns, Herschel/PACS at 70 microns, and ALMA at 870 microns. The WFC3 and MIPS data are published here for the first time. We focus our modeling on the outer part of this disk, consisting of a parent body ring and a halo of small grains. First, we confirm that a model using astronomical silicates cannot simultaneously fit the thermal and scattered light data. Next, we use a simple, generic function for the optical constants to show that varying the dust composition can improve the fit substantially. Finally, we model the dust as a mixture of the most plausible debris constituents: astronomical silicates, water ice, organic refractory material, and vacuum. We achieve a good fit to all datasets with grains composed predominantly of silicates and organics, while ice and vacuum are, at most, present in small amounts. This composition is similar to one derived from previous work on the HR 4796A disk. Our model also fits the thermal SED, scattered light colors, and high-resolution mid-IR data from T-ReCS for this disk. Additionally, we show that sub-blowout grains are a necessary component of the halo.

We find that the initial dust masses in planetary debris disks are correlated with the metallicities of their central stars. We compiled a large sample of systems, including Spitzer, the Herschel DUNES and DEBRIS surveys, and WISE debris disk candidates. We also merged 33 metallicity catalogs to provide homogeneous [Fe/H] and $\sigma_{[Fe/H]}$ values. We analyzed this merged sample, including 222 detected disks (74 warm and 148 cold) around a total of 187 systems (some with multiple components) and 440 disks with only upper limits (125 warm and 315 cold), around a total of 360 systems. The disk dust masses at a common early evolutionary point in time were determined using our numerical disk evolutionary code, evolving a unique model for each of the 662 disks backward to an age of 1 Myr. We find that disk-bearing stars seldom have metallicities less than [Fe/H] = -0.2 and that the distribution of warm component masses lacks examples with large mass around stars of low metallicity ([Fe/H] < -0.085). Previous efforts to find a correlation have been largely unsuccessful; the primary improvements supporting our result are: 1.) basing the study on dust masses, not just infrared excess detections; 2.) including upper limits on dust mass in a quantitative way; 3.) accounting for the evolution of debris disk excesses as systems age; 4.) accounting fully for the range of uncertainties in metallicity measurements; and 5.) having a statistically large enough sample.

A significant fraction of main sequence stars observed interferometrically in the near infrared have slightly extended components that have been attributed to very hot dust. To match the spectrum appears to require the presence of large numbers of very small (< 200 nm in radius) dust grains. However, particularly for the hotter stars, it has been unclear how such grains can be retained close to the star against radiation pressure force. We find that the expected weak stellar magnetic fields are sufficient to trap nm-sized dust grains in epicyclic orbits for a few weeks or longer, sufficient to account for the hot excess emission. Our models provide a natural explanation for the requirement that the hot excess dust grains be smaller than 200 nm. They also suggest that magnetic trapping is more effective for rapidly rotating stars, consistent with the average vsini measurements of stars with hot excesses being larger (at about 2 sigma) than those for stars without such excesses.

Observations of debris disks allow for the study of planetary systems, even where planets have not been detected. However, debris disks are often only characterized by unresolved infrared excesses that resemble featureless blackbodies, and the location of the emitting dust is uncertain due to a degeneracy with the dust grain properties. Here we characterize the Spitzer IRS spectra of 22 debris disks exhibiting 10 micron silicate emission features. Such features arise from small warm dust grains, and their presence can significantly constrain the orbital location of the emitting debris. We find that these features can be explained by the presence of an additional dust component in the terrestrial zones of the planetary systems, i.e. an exozodiacal belt. Aside from possessing exozodiacal dust, these debris disks are not particularly unique; their minimum grain sizes are consistent with the blowout sizes of their systems, and their brightnesses are comparable to those of featureless warm debris disks. These disks are in systems with a range of ages, although the older systems with features are found only around A-type stars. The features in young systems may be signatures of terrestrial planet formation. Analyzing the spectra of unresolved debris disks with emission features may be one of the simplest and most accessible ways to study the terrestrial regions of planetary systems.

Cold debris disks trace the limits of planet formation or migration in the outer regions of planetary systems, and thus have the potential to answer many of the outstanding questions in wide-orbit planet formation and evolution. We characterized the infrared excess spectral energy distributions of 174 cold debris disks around 546 main-sequence stars observed by both Spitzer IRS and MIPS. We found a trend between the temperature of the inner edges of cold debris disks and the stellar type of the stars they orbit. This argues against the importance of strictly temperature-dependent processes (e.g. non-water ice lines) in setting the dimensions of cold debris disks. Also, we found no evidence that delayed stirring causes the trend. The trend may result from outward planet migration that traces the extent of the primordial protoplanetary disk, or it may result from planet formation that halts at an orbital radius limited by the efficiency of core accretion.

The Galactic Arecibo L-band Feed Array HI (GALFA-HI) survey is mapping the entire Arecibo sky at 21-cm, over a velocity range of -700 to +700 km/s (LSR), at a velocity resolution of 0.18 km/s and an angular resolution of 3.5 arcmin. The unprecedented resolution and sensitivity of the GALFA-HI survey have resulted in the detection of many isolated, very compact HI clouds at low Galactic velocities which are distinctly separated from the HI disk emission. In the limited area of ~4600 deg$^2$ searched so far, we have detected 96 such compact clouds. The detected clouds are cold with kinetic temperature less than 300 K. Moreover, they are quite compact and faint, with median values of 5 arcmin in angular size, 0.75 K in peak brightness temperature, and $5\times10^{18}$ cm$^{-2}$ in HI column density. From the modeling of spatial and velocity distributions of the whole compact cloud population, we find that the bulk of clouds are related to the Galactic disk, and are within a few kpc distance. We present properties of the compact clouds sample and discuss various possible scenarios for the origin of this clouds population and its role in the Galactic interstellar medium studies.

The Galactic Arecibo L-band Feed Array HI (GALFA-HI) survey is mapping the entire Arecibo sky at 21-cm, over a velocity range of -700 to +700 km/s (LSR), at a velocity resolution of 0.18 km/s and a spatial resolution of 3.5 arcmin. The unprecedented resolution and sensitivity of the GALFA-HI survey have resulted in the detection of numerous isolated, very compact HI clouds at low Galactic velocities, which are distinctly separated from the HI disk emission. In the limited area of ~4600 deg$^2$ surveyed so far, we have detected 96 of such compact clouds. The detected clouds are cold with a median T$_{k,max}$ (the kinetic temperature in the case in which there is no non-thermal broadening) of 300 K. Moreover, these clouds are quite compact and faint, with median values of 5 arcmin in angular size, 0.75 K in peak brightness temperature, and $5 \times 10^{18}$ cm$^{-2}$ in HI column density. Most of the clouds deviate from Galactic rotation at the 20-30 km/s level, and a significant fraction show evidence for a multiphase medium and velocity gradients. No counterparts for these clouds were found in other wavebands. From the modeling of spatial and velocity distributions of the whole compact cloud population, we find that the bulk of the compact clouds are related to the Galactic disk, and their distances are likely to be in the range of 0.1 to a few kpc. We discuss various possible scenarios for the formation and maintenance of this cloud population and its significance for Galactic ISM studies.