The Leftovers of Planet Formation:
Small Body Populations of Our Solar System and Exoplanet Systems
by
Samantha Lawler
B.S., California Institute of Technology, 2005 M.A., Wesleyan University, 2009
A THESIS SUBMITTED IN PARTIAL FULFILLMENT OFTHEREQUIREMENTSFORTHEDEGREEOF
Doctor of Philosophy
in
THE FACULTY OF GRADUATE STUDIES (Astronomy)
The University Of British Columbia (Vancouver)
July 2013
c Samantha Lawler, 2013 Abstract
The small body populations within a planetary system give information about the planet formation and migration history of the system. In our Solar System, we study these bodies (asteroids, comets, and trans-Neptunian objects), by directly observing them in reflected light. In other solar systems, dust traces the position of the planetesimal belts that produce it, and is observed as an excess above the stellar flux in the infrared. The dust is visible and not the planetesimals because of the much greater cross-sectional surface area of a swarm of dust particles. In this thesis, both leftover large planetesimals in our Solar System and dust around other stars are investigated. Data from the Canada-France Ecliptic Plane Survey (CFEPS) are used to mea- sure the absolute populations of trans-Neptunian objects (TNOs) in mean-motion resonances with Neptune, as well as constrain the internal orbital element dis- tributions. Detection biases play a critical role because phase relationships with Neptune make object discovery more likely at certain longitudes. The plutinos (objects in the 3:2 resonance) are given particular attention because the presence of the secular Kozai resonance within the mean-motion resonance causes different detection biases that need to be accounted for to properly debias surveys that in- clude detections of plutinos. Because the TNOs that are trapped in mean-motion resonances with Neptune were likely emplaced there during planet migration late in the giant planet formation process, the structure within and relative populations of the resonances should be a diagnostic of the timescale and method of giant planet migration.
ii Exoplanet systems that host several rocky planets are those that did not expe- rience giant planet migration, and thus are likely to host planetesimal belts which should be detectable as debris disks. The Kepler Mission has detected a host of such systems, and we use data from the Wide-field Infrared Survey Explorer (WISE) Mission to search for debris disks around these stars. Though we tenta- tively detect more excesses toward these stars than would be expected, contami- nation from warm dust in the Milky Way Galaxy makes detection unreliable for these systems, and will have to await future infrared space telescopes.
iii Preface
The text of this dissertation includes modified reprints of previously published material as listed below.
Chapter 2 (published):
B. Gladman, S. M. Lawler, J.-M. Petit, J. Kavelaars, R. L. Jones, J. Wm. • Parker, C. Van Laerhoven, P. Nicholson, P. Rosselot, A. Bieryla, and M. L. N. Ashby, The Resonant Trans-Neptunian Populations, AJ 144, 23 (2012).
This paper presents an analysis of resonant TNOs detected in a large (five- year) observational program by Gladman et al. I contributed to this project by performing all of the modelling that was required to debias the observations and regain the orbital element distributions and true populations for the TNOs that were resonant. This involved planning and running many hundreds of simulations with varying parameters, comparing the resulting models to observed resonant objects, plotting the results, and doing statistical tests. The telescope time was originally proposed for by B. Gladman, J.-M. Petit, and J. Kavelaars. The observations for the survey were performed at the Canada-France-Hawaii Tele- scope during 2003-2008 by B. Gladman, J.-M. Petit, J. Kavelaars, R. L. Jones, J. Wm. Parker, and P. Nicholson. These observations were analyzed and orbits were fit by B. Gladman, R. L. Jones, and C. Van Laerhoven. In addition to writing my own code to analyze and display the results, I utilized the Survey Simulator code originally written by J.-M. Petit, J. Kavelaars, and B. Gladman. To test how
iv well the models fit the data, I modified statistical analysis code originally written by R. L. Jones. The bulk of the paper was written by B. Gladman, with editing by myself, J.-M. Petit, and J. Kavelaars. Section 2.3 and all of the figure captions were written by me, and all of the figures were made by me.
Chapter 3 (published):
S. M. Lawler and B. Gladman, Plutino Detection Biases, Including the • Kozai Resonance, AJ 146, 6 (2013).
I wrote the entire text, performed all of the modelling and simulated surveys, and made all of the figures. B. Gladman provided commentary, including the initial idea to separate this project from the larger resonant paper (Chapter 2). B. Gladman also helped my understanding of the complicated Kozai effect and dynamical disturbing function, as well as proofreading and editing.
Chapter 4 (published):
S. M. Lawler and B. Gladman, Debris Disks in Kepler Exoplanet Systems, • ApJ 752, 53 (2012) .
I had the idea for and formulated the project, wrote the entire text, wrote all of the code for analysis, and made all of the figures. B. Gladman provided understanding of the dynamical background, timescales, and interpretation of results, as well as providing feedback and proofreading.
Chapter 5 (submitted):
S. M. Lawler and B. Gladman, Galactic Dust Limits Searches for Asteroid • Belts in Kepler Exoplanet Systems, submitted July 2013.
This was a follow-up to Chapter 4, using the slightly modified code from the first project with an expanded WISE dataset and more Kepler exoplanet hosts, as well as dealing with the fact that another paper had been published calling in
v to question all of the debris disks claimed in Chapter 4. B. Gladman provided extensive guidance and advice on how to proceed, as well as many rounds of editing and proofreading.
vi Table of Contents
Abstract...... ii
Preface ...... iv
TableofContents ...... vii
ListofTables...... xii
ListofFigures ...... xiii
Glossary ...... xvii
Acknowledgments ...... xxi
Dedication ...... xxii
1 Introduction ...... 1 1.1 ABriefOverviewofPlanetFormation ...... 3 1.2 The Kuiper Belt: Our Solar System’s Outer Planetesimal Belt... 7 1.2.1 DiscoveryHistory...... 8 1.2.2 TNOClassification ...... 10 1.2.3 GiantPlanetMigrationinourSolarSystem ...... 14 1.3 Debris Disks: Indirect Measurement of Extrasolar Planetesimal Belts...... 16
vii 1.3.1 DebrisDiskObservations...... 16 1.3.2 The Evolution of a Protoplanetary Disk into a Debris Disk 20 1.3.3 RemovalandReplenishmentofDust ...... 22 1.3.4 OurSolarSystemasaDebrisDisk ...... 26 1.4 ARelationshipBetweenPlanetsandDebrisDisks? ...... 30 1.4.1 DebrisDisksas aMarkerforDynamicalStability? . . . . 32 1.5 DatasetsUsedinthisThesis...... 33 1.5.1 SolarSystemObservations ...... 33 1.5.2 ExtrasolarObservations...... 39 1.6 ThesisOutline...... 40
2 The Resonant Trans-Neptunian Populations ...... 42 2.1 Introduction ...... 42 2.1.1 ResonanceDynamics ...... 44 2.2 ResonantCFEPSObjects ...... 48 2.3 CFEPS SurveySimulationofaResonant Population ...... 54 2.3.1 Building the Simulated Populations for Each Resonance . 55 2.3.2 Simulatingthe5:2,n:3,andn:4Populations ...... 56 2.3.3 SimulatingthePlutinoPopulation ...... 58 2.3.4 Simulatingthen:1Populations ...... 62 2.4 ThePlutinos(3:2Resonators)...... 64 2.4.1 ThePlutinoInclinationDistribution ...... 66 2.4.2 ThePlutinoKozaiSubcomponent ...... 69 2.4.3 ThePlutino Size and Eccentricity Distribution ...... 70 2.4.4 PlutinoLibrationAmplitudes...... 75 2.5 ThePopulationofPlutinos ...... 76 2.6 The5:2Resonance ...... 77 2.7 Then:3Resonances...... 81 2.7.1 The4:3Resonance ...... 82 2.7.2 The5:3Resonance ...... 83 2.7.3 The7:3Resonance ...... 84
viii 2.8 Then:4Resonances...... 85 2.8.1 The5:4Resonance ...... 85 2.8.2 The7:4Resonance ...... 85 2.9 Then:1Resonances...... 86 2.9.1 TheTwotinos ...... 87 2.9.2 The3:1Resonance ...... 91 2.9.3 The5:1Resonance ...... 92 2.10 ResonanceswithNoCFEPSDetections ...... 93 2.10.1 TheNeptuneTrojans ...... 93 2.10.2 The4:1Resonance ...... 95 2.11 CFEPSComparisontoaCosmogonicModel...... 96 2.12 ResonantPopulations ...... 99 2.13Discussion...... 105
3 Detection Biases for the Plutinos, Including the Kozai Resonance . . 110 3.1 Introduction ...... 110 3.2 ResonantDynamics ...... 114 3.3 Non-KozaiPlutinos ...... 116 3.3.1 0◦ LibrationAmplitudeToyModel ...... 116 3.3.2 95◦ LibrationAmplitudeToyModel ...... 116 3.3.3 ARealisticPlutinoDistribution ...... 118 3.4 KozaiPlutinoDynamics ...... 121 3.4.1 On-SkyDetectionBiasesforKozaiPlutinos ...... 124 3.5 SimulatingthePlutinos ...... 127 3.5.1 Non-KozaiPlutinos ...... 129 3.5.2 KozaiPlutinos...... 130 3.6 On-SkyBiases...... 131 3.6.1 EclipticLatitudeDistributionofDetections ...... 133
3.6.2 fkoz On-Sky...... 134 3.6.3 SizeDistributionEffects ...... 135 3.6.4 ExampleSimulatedSurveys ...... 135
ix 3.7 ComparisonwithPreviousLiterature ...... 141
3.7.1 The Kozai Fraction fkoz ...... 141 3.7.2 DistributionofKozaiParameters ...... 145 3.8 Conclusion ...... 148
4 DebrisDisksinKeplerExoplanetSystems...... 149 4.1 Introduction ...... 149 4.2 WISEData ...... 151 4.2.1 DefiningOurSample ...... 151 4.2.2 FindingExcessesintheWISEData ...... 152 4.3 CandidateDebrisDiskDetections ...... 155 4.3.1 Warm Dust Systems ( 300–500K) ...... 164 ∼ 4.3.2 Cool Dust Systems ( 100K)...... 168 ∼ 4.3.3 HotDust( 1000K)inaMultiplanetSystem ...... 170 ∼ 4.4 DustProduction...... 171 4.5 Discussion...... 174
5 Galactic Dust Limits Detection of Asteroid Belts in Kepler Exo- planetSystems ...... 176 5.1 Introduction ...... 176 5.2 OurSample ...... 177 5.2.1 WISEData ...... 178 5.2.2 FindingExcessesintheWISEData ...... 178 5.2.3 SystemswithApparentExcess ...... 180 5.2.4 ComparisonwithPreviousPapers ...... 184 5.3 StrawmanDustModels ...... 185 5.3.1 PropertiesofSystemswithIRExcess ...... 187 5.4 ContaminationbyBackgroundGalacticCirrus ...... 188 5.4.1 GalacticDustintheWISEandIRASImages ...... 189 5.4.2 IndependentMeasurementsofWISEPhotometry . . . . . 193 5.4.3 CalculatedDustTemperatures ...... 197
x 5.5 Discussion...... 200 5.5.1 ComparisonwithOtherSurveys ...... 200 5.5.2 FutureObservations...... 202
6 FutureWork ...... 203 6.1 SolarSystem ...... 203 6.2 ExoplanetsandDebrisDisks ...... 207
Bibliography ...... 211
xi List of Tables
Table2.1 CFEPS3:2(Plutinos)and5:2Resonators...... 50 Table2.2 CFEPS TNOsinResonancesOtherthan3:2and5:2 ...... 52 Table2.3 ResonancePopulations...... 100
Table3.1 SimulationsoftheObservedKozaiFraction ...... 140 Table 3.2 Kozai Fraction Measurements in the Literature ...... 142
Table 4.1 FGK Stars with No Significant (>5 δ)Excess...... 157 Table4.2 DebrisDiskCandidates...... 160 Table4.3 WISEDataforDebrisDiskCandidates...... 161
Table 5.1 Data for High Excess Tier Systems and Strawman Dust Models 183
xii List of Figures
Figure 1.1 Orbital Elements and Sizes of Catalogued Minor Planets ... 4 Figure1.2 OrbitalElementsofaSolarSystemBody ...... 11 Figure1.3 NomenclatureoftheOuterSolarSystem ...... 13 Figure 1.4 An Assortment of Directly Imaged Debris Disks ...... 18 Figure 1.5 Spectral Energy Distribution for Different ExampleDebrisDisks 19 Figure 1.6 The Fraction of Stars in Young Named Clusters with Disks . . 21 Figure1.7 PRDrag...... 23 Figure 1.8 Fractional Luminosityof Known Debris Disks ...... 25 Figure 1.9 Predictions for the Kuiper belt SED Based on a Dynamical Model...... 27 Figure 1.10 Predictions for the Kuiper belt SED using the MPC Database . 28 Figure 1.11 SED of τ Ceti, a Sun-Like Star with a Massive Kuiper Belt Analogue ...... 29 Figure 1.12 CFEPS Observing Blocks on the Celestial Sphere ...... 34 Figure 1.13 Orbital Elements (a, e, and i) for 176 Objects Detected by CFEPS ...... 35 Figure 1.14 Orbital Elements (a, e, and i) for the CFEPS L7 Debiased KuiperBeltModel ...... 36 Figure 1.15 The CFEPS L7 Debiased Kuiper Belt Model Shown in the ContextofOurSolarSystem,asSeenfromAbove...... 37 Figure 1.16 The CFEPS L7 Debiased Kuiper Belt Model Shown in the ContextofOurSolar System,asSeenfromEdge-On...... 38
xiii Figure 2.1 Toy Models Giving Ecliptic Projections of ResonantTNOs . . 46 Figure 2.2 Model Predictions for the Detection Distance Distribution of Non-KozaiandKozaiPlutinos ...... 60 Figure 2.3 The Hamiltonian Phase Space for the Set of Kozai Librators . 61 Figure 2.4 The Range of Libration Amplitudes, Libration Centers, and Eccentricities Chosen in Our Model for the Symmetric and AsymmetricIslandsinthe2:1Resonances ...... 64 Figure2.5 EclipticProjectionofthePlutinos ...... 65 Figure 2.6 Cumulative Distributions of the Five Variables on which we PerformStatisticalAnalysis,forthePlutinos...... 68 Figure 2.7 Confidence Regions in the Plutinos Orbital Model Parameter Space ...... 72 Figure 2.8 The CumulativePlutino Distance at Detection Distribution for a Model with Size Distribution Exponent α =0.55 ...... 74 Figure 2.9 Confidence Regions for the 5:2 Parameter Space, for α =0.9 . 80 Figure 2.10 Comparison Between CFEPS Plutino Detections and Simu- latedDetectionsfromtheNiceModel ...... 98 Figure 2.11 The Apparent versus Debiased Resonant TNOs in a, e, and i . 101 Figure 2.12 Population Estimates for the 3:2, 5:2, and 2:1 Resonances . . . 102
Figure 3.1 A Simplified Orbital Diagram of the Kozai Effect ...... 112 Figure 3.2 Orbital Integrations of the Kozai Plutino 28978 Ixion, the
CFEPS-Discovered Non-Kozai Plutino L4h15 (2004 HB79),
and a Non-Resonant TNO (2004 PA112) ...... 115 Figure3.3 0◦ LibrationAmplitudeToyModel ...... 117 Figure3.4 95◦ LibrationAmplitudePlutinoToyModel ...... 119 Figure3.5 ARealisticNon-KozaiPlutinoModel ...... 120 Figure 3.6 The Detection Density on the Sky using a Realistic Distribu- tionofNon-KozaiPlutinos ...... 121 Figure 3.7 3 Different Hamiltonian Level Surfaces for Kozai Plutinos . . 122 ◦ Figure3.8 95 Kozai Libration Amplitude (Aω)ToyModel...... 125
xiv Figure 3.9 Distance at Detection for Kozai and Non-Kozai Plutinos in a Simulated,Flux-Limited,All-SkySurvey ...... 126 Figure 3.10 Averaged Distance at Detection for Kozai and Non-Kozai Plutinos in a Simulated, Flux-Limited, All-Sky Survey . . . . 127 Figure 3.11 Distribution of Simulated Plutinos in Semimajor Axis, Eccen- tricity,andInclination...... 128 Figure 3.12 Relative Density of Detections on the Sky in an All-Sky Sur- veyusingourPlutinoModel ...... 132 Figure 3.13 Relative Density of Detections on the Sky in an All-Sky Sur- veyonlyusingtheKozaiComponentofourModel ...... 133 Figure 3.14 Stacked Histograms of the Ecliptic Latitude Distribution of DetectedPlutinos ...... 134 obs Figure 3.15 The Observed Kozai Fraction fkoz ...... 136 Figure 3.16 Relative Density of Plutino Detections using Three Different Values of α ...... 137 obs true Figure 3.17 fkoz at Different Points on the Sky, using fkoz = 10%, and Three Different Values of α ...... 138
Figure 3.18 Imax and Aω Values for Real and Simulated Kozai Plutinos . . 145
Figure 3.19 Imax and Aω Cumulative Distributions for Real and Simulated KozaiPlutinos...... 146
Figure4.1 TheExcessSignificancesforEachBand ...... 153 Figure 4.2 Histograms of Each of the Four Stellar Parameters Reported by the Kepler Team...... 154 Figure 4.3 Eight Stars with WISE Data that are Consistent with the Stel- larPhotosphere ...... 158 Figure 4.4 Eight Stars with WISE Data that are Consistent with the Stel- larPhotosphere(continued)...... 159 Figure 4.5 The Derived Masses, Fractional Luminosities, and Tempera- tures of Dust Rings in our Sample Compared to Other Known Systems...... 163
xv Figure4.6 SEDforKOI1099 ...... 165 Figure4.7 SEDsforStarswithWarmExcesses ...... 166 Figure4.8 SEDsforStarswithCoolExcesses ...... 168 Figure4.9 SEDforKOI904 ...... 169
Figure5.1 ExcessSignificancesinEachBand ...... 181 Figure5.2 SEDsfortheHighExcessTierSystems ...... 186 Figure 5.3 SEDs for the High Excess Tier Systems (continued) ...... 187 Figure 5.4 Locations of Kepler Targets Analyzed in this chapter within the Kepler FieldofView ...... 190 Figure5.5 WISEandIRASImagesofExcessHosts ...... 192 Figure 5.6 Photometric Curve of Growth the Six High Excess Tier Stars . 194 Figure5.7 10′-wide Wide-field Infrared Survey Explorer (WISE) images centeredonKOI2130...... 196 Figure 5.8 The W3/W4 Flux Ratio for Different Blackbody Temperatures 198 Figure5.9 ModelSpectrumoftheMilkyWayGalaxy...... 199
Figure 5.10 The derived Strawman Masses, Ldust/L∗, and Temperatures ofDustRingsComparedtoOtherKnownSystems ...... 201
xvi Glossary
ALMA Atacama Large Millimeter Array
CFEPS Canada-France Ecliptic Plane Survey
CFHT Canada-France-Hawaii Telescope
COBE Cosmic Background Explorer
DES Deep Ecliptic Survey
DEBRIS Disk Emission via a Bias-free Reconnaissance in the Infrared/Submil- limetre
E-ELT European Extremely Large Telescope
GMT Giant Magellan Telescope
IAU International Astronomical Union
IR Infrared
IRAF Image Reduction and Analysis Facility
IRAS Infrared Astronomical Satellite
IPAC Infrared Processing and Analysis Center
JCMT James Clerk Maxwell Telescope
xvii JFC Jupiter Family Comet
JWST James Webb Space Telescope
KOI Kepler Object of Interest
LSST Large Synoptic Survey Telescope
MPC Minor Planet Center
NASA National Aeronautics and Space Administration
PAH Polycyclic Aromatic Hydrocarbon
Pan-STARRS Panoramic Survey Telescope & Rapid Response System
PR Poynting-Robertson
SDSS Sloan Digital Sky Survey
SED Spectral Energy Distribution
SNR Signal-to-Noise Ratio
TMT Thirty Meter Telescope
TNO Trans-Neptunian Object
2MASS Two Micron All Sky Survey
VLT Very Large Telescope
WISE Wide-field Infrared Survey Explorer
xviii AU - Astronomical Unit. Roughly the mean Earth-Sun distance. Defined as 149,597,870,700 metres, exactly. pc - Parsec. The distance at which one AU subtends one arcsecond. 206,264.8 AU or 3.0856 1016 m. × Jy - Jansky. A unit of flux density commonly used in radio, submillimeter, in- −26 W frared, and x-ray astronomy. 1Jy =10 m2Hz .
Ldust/L∗ - Fractional luminosity, commonly used to measure the brightness of debris disks by comparing the luminosity of the dust disk Ldust to the luminosity of the star L∗.