Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

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.

We explore the evolution of the mass distribution of dust in collision-dominated debris disks, using the collisional code introduced in our previous paper. We analyze the equilibrium distribution and its dependence on model parameters by evolving over 100 models to 10 Gyr. With our numerical models, we confirm that systems reach collisional equilibrium with a mass distribution that is steeper than the traditional solution by Dohnanyi (1969). Our model yields a quasi steady-state slope of n(m) ~ m^{-1.88} [n(a) ~ a^{-3.65}] as a robust solution for a wide range of possible model parameters. We also show that a simple power-law function can be an appropriate approximation for the mass distribution of particles in certain regimes. The steeper solution has observable effects in the submillimeter and millimeter wavelength regimes of the electromagnetic spectrum. We assemble data for nine debris disks that have been observed at these wavelengths and, using a simplified absorption efficiency model, show that the predicted slope of the particle mass distribution generates SEDs that are in agreement with the observed ones.

Core-collapse supernovae (CC SNe), especially Type II-Plateau ones, are thought to be important contributors to cosmic dust production. SN 2004dj, one of the closest and brightest SN since 1987A, offered a good opportunity to examine dust formation processes. To find signs of newly formed dust, we analyze all available mid-infrared (MIR) archival data from the Spitzer Space Telescope. We re-reduce and analyze data from IRAC, MIPS and IRS instruments obtained between +98 and +1381 days after explosion, generating light curves and spectra for each epoch. Observed spectral energy distributions are fitted with both analytic and numerical models, using the radiative-transfer code MOCASSIN for the latter ones. We also use imaging polarimetric data obtained at +425 days by the Hubble Space Telescope. We present convincing evidence of dust formation in the ejecta of SN 2004dj from MIR light curves and spectra. Significant MIR excess flux is detected in all bands between 3.6 and 24 um. In the optical, a ~0.8% polarization is also detected at a 2-sigma level, which exceeds the interstellar polarization in that direction. Our analysis shows that the freshly-formed dust around SN 2004dj can be modeled assuming a nearly spherical shell containing amorphous carbon grains, cooling from ~700 K to ~400 K between +267 and +1246 days. Persistent excess flux has been found above 10 um, which is explained by a cold (~115 K) dust component. If this cold dust is of circumstellar origin, it is likely to be condensed in a cool, dense shell between the forward and reverse shocks. Pre-existing circumstellar dust is less likely, but cannot be ruled out. An upper limit of ~8x10^{-4} M_sun was derived for the dust mass, which is similar to previously published values for other dust-producing SNe.

We present a model that describes stellar infrared excesses due to heating of the interstellar (IS) dust by a hot star passing through a diffuse IS cloud. This model is applied to six lambda Bootis stars with infrared excesses. Plausible values for the IS medium (ISM) density and relative velocity between the cloud and the star yield fits to the excess emission. This result is consistent with the diffusion/accretion hypothesis that lambda Bootis stars (A- to F-type stars with large underabundances of Fe-peak elements) owe their characteristics to interactions with the ISM. This proposal invokes radiation pressure from the star to repel the IS dust and excavate a paraboloidal dust cavity in the IS cloud, while the metal-poor gas is accreted onto the stellar photosphere. However, the measurements of the infrared excesses can also be fit by planetary debris disk models. A more detailed consideration of the conditions to produce lambda Bootis characteristics indicates that the majority of infrared-excess stars within the Local Bubble probably have debris disks. Nevertheless, more distant stars may often have excesses due to heating of interstellar material such as in our model.

We have discovered a bow shock shaped mid-infrared excess region in front of delta Velorum using 24 micron observations obtained with the Multiband Imaging Photometer for Spitzer (MIPS). The excess has been classified as a debris disk from previous infrared observations. Although the bow shock morphology was only detected in the 24 micron observations, its excess was also resolved at 70 micron. We show that the stellar heating of an ambient interstellar medium (ISM) cloud can produce the measured flux. Since delta Velorum was classified as a debris disk star previously, our discovery may call into question the same classification of other stars. We model the interaction of the star and ISM, producing images that show the same geometry and surface brightness as is observed. The modeled ISM is 15 times overdense relative to the average Local Bubble value, which is surprising considering the close proximity (24 pc) of delta Velorum. The abundance anomalies of lambda Bootis stars have been previously explained as arising from the same type of interaction of stars with the ISM. Low resolution optical spectra of delta Velorum show that it does not belong to this stellar class. The star therefore is an interesting testbed for the ISM accretion theory of the lambda Bootis phenomenon.