Origin of the peculiar eccentricity distribution of the inner cold Kuiper belt

Dawson and Murray-Clay (2012) pointed out that the inner part of the cold population in the Kuiper belt (that with semi major axis a<43.5 AU) has orbital eccentricities significantly smaller than the limit imposed by stability constraints. Here, we confirm their result by looking at the orbital distribution and stability properties in proper element space. We show that the observed distribution could have been produced by the slow sweeping of the 4/7 mean motion resonance with Neptune that accompanied the end of Neptune's migration process. The orbital distribution of the hot Kuiper belt is not significantly affected in this process, for the reasons discussed in the main text. Therefore, the peculiar eccentricity distribution of the inner cold population can not be unequivocally interpreted as evidence that the cold population formed in-situ and was only moderately excited in eccentricity; it can simply be the signature of Neptune's radial motion, starting from a moderately eccentric orbit. We discuss how this agrees with a scenario of giant planet evolution following a dynamical instability and, possibly, with the radial transport of the cold population.Comment: in press in Icaru

The Correlation Between Metallicity and Debris Disk mass

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.Comment: 13 pages, 7 figures, accepted for publication to Ap

The Collisional Evolution of Debris Disks

We explore the collisional decay of disk mass and infrared emission in debris disks. With models, we show that the rate of the decay varies throughout the evolution of the disks, increasing its rate up to a certain point, which is followed by a leveling off to a slower value. The total disk mass falls off ~ t^-0.35 at its fastest point (where t is time) for our reference model, while the dust mass and its proxy -- the infrared excess emission -- fades significantly faster (~ t^-0.8). These later level off to a decay rate of M_tot(t) ~ t^-0.08 and M_dust(t) or L_ir(t) ~ t^-0.6. This is slower than the ~ t^-1 decay given for all three system parameters by traditional analytic models. We also compile an extensive catalog of Spitzer and Herschel 24, 70, and 100 micron observations. Assuming a log-normal distribution of initial disk masses, we generate model population decay curves for the fraction of debris disk harboring stars observed at 24 micron and also model the distribution of measured excesses at the far-IR wavelengths (70-100 micron) at certain age regimes. We show general agreement at 24 micron between the decay of our numerical collisional population synthesis model and observations up to a Gyr. We associate offsets above a Gyr to stochastic events in a few select systems. We cannot fit the decay in the far infrared convincingly with grain strength properties appropriate for silicates, but those of water ice give fits more consistent with the observations.Comment: 32 pages, 16 figures, emulateapj format, Accepted for publication in Ap

The Decay of Debris Disks around Solar-Type Stars

We present a Spitzer MIPS study of the decay of debris disk excesses at 24 and 70 $\mu$m for 255 stars of types F4 - K2. We have used multiple tests, including consistency between chromospheric and X-ray activity and placement on the HR diagram, to assign accurate stellar ages. Within this spectral type range, at 24 $\mu$m, $13.6 \pm 2.8 \%$ of the stars younger than 5 Gyr have excesses at the 3$\sigma$ level or more, while none of the older stars do, confirming previous work. At 70 $\mu$m, $22.5 \pm 3.6\%$ of the younger stars have excesses at $\ge$ 3 $\sigma$ significance, while only $4.7^{+3.7}_{-2.2}$% of the older stars do. To characterize the far infrared behavior of debris disks more robustly, we double the sample by including stars from the DEBRIS and DUNES surveys. For the F4 - K4 stars in this combined sample, there is only a weak (statistically not significant) trend in the incidence of far infrared excess with spectral type (detected fractions of 21.9$^{+4.8}_{-4.3}\%$, late F; 16.5$^{+3.9}_{-3.3}\%$, G; and 16.9$^{+6.3}_{-5.0}\%$, early K). Taking this spectral type range together, there is a significant decline between 3 and 4.5 Gyr in the incidence of excesses with fractional luminosities just under $10^{-5}$. There is an indication that the timescale for decay of infrared excesses varies roughly inversely with the fractional luminosity. This behavior is consistent with theoretical expectations for passive evolution. However, more excesses are detected around the oldest stars than is expected from passive evolution, suggesting that there is late-phase dynamical activity around these stars.Comment: 46 pages. 7 figures. Accepted to Ap

A Comprehensive Dust Model Applied to the Resolved Beta Pictoris Debris Disk from Optical to Radio Wavelengths

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.Comment: 23 pages, 20 figures, accepted to Ap

The first 40 million years of circumstellar disk evolution: the signature of terrestrial planet formation

We characterize the first 40 Myr of evolution of circumstellar disks through a unified study of the infrared properties of members of young clusters and associations with ages from 2 Myr up to ~ 40 Myr: NGC 1333, NGC 1960, NGC 2232, NGC 2244, NGC 2362, NGC 2547, IC 348, IC 2395, IC 4665, Chamaeleon I, Orion OB1a and OB1b, Taurus, the \b{eta} Pictoris Moving Group, \r{ho} Ophiuchi, and the associations of Argus, Carina, Columba, Scorpius-Centaurus, and Tucana-Horologium. Our work features: 1.) a filtering technique to flag noisy backgrounds, 2.) a method based on the probability distribution of deflections, P(D), to obtain statistically valid photometry for faint sources, and 3.) use of the evolutionary trend of transitional disks to constrain the overall behavior of bright disks. We find that the fraction of disks three or more times brighter than the stellar photospheres at 24 {\mu}m decays relatively slowly initially and then much more rapidly by ~ 10 Myr. However, there is a continuing component until ~ 35 Myr, probably due primarily to massive clouds of debris generated in giant impacts during the oligarchic/chaotic growth phases of terrestrial planets. If the contribution from primordial disks is excluded, the evolution of the incidence of these oligarchic/chaotic debris disks can be described empirically by a log-normal function with the peak at 12 - 20 Myr, including ~ 13 % of the original population, and with a post-peak mean duration of 10 - 20 Myr.Comment: accepted for publication, the Astrophysical Journal (2017