De-biased populations of Kuiper Belt objects from the Deep Ecliptic Survey
The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters.
Citation Adams, E. R., A. A. S. Gulbis, J. L. Elliot, S. D. Benecchi, M. W. Buie, D. E. Trilling, and L. H. Wasserman. “De-Biased Populations of Kuiper Belt Objects from the Deep Ecliptic Survey.” The Astronomical Journal 148, no. 3 (August 18, 2014): 55. © 2014 The American Astronomical Society
As Published http://dx.doi.org/10.1088/0004-6256/148/3/55
Publisher IOP Publishing
Version Final published version
Citable link http://hdl.handle.net/1721.1/92725
Terms of Use Article is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use. The Astronomical Journal, 148:55 (17pp), 2014 September doi:10.1088/0004-6256/148/3/55 C 2014. The American Astronomical Society. All rights reserved. Printed in the U.S.A.
DE-BIASED POPULATIONS OF KUIPER BELT OBJECTS FROM THE DEEP ECLIPTIC SURVEY
E. R. Adams1,A.A.S.Gulbis2,3, J. L. Elliot3,8, S. D. Benecchi1,4,M.W.Buie5, D. E. Trilling6, and L. H. Wasserman7 1 Planetary Science Institute, 1700 East Fort Lowell, Suite 106, Tucson, AZ 85719, USA 2 The Southern African Large Telescope and South African Astronomical Observatory, Cape Town, 7935, South Africa 3 Department of Earth, Atmospheric, and Planetary Sciences, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA 4 Carnegie Institution of Washington, Department of Terrestrial Magnetism, 5241 Broad Branch Road NW, Washington, DC 20015-1305, USA 5 Southwest Research Institute, 6220 Culebra Road, San Antonio, TX 78238, USA 6 Department of Physics & Astronomy, Northern Arizona University, S San Francisco Street, Flagstaff, AZ 86011, USA 7 Lowell Observatory, 1400 W. Mars Hill Road, Flagstaff, AZ 86001, USA Received 2013 October 11; accepted 2014 June 3; published 2014 August 18
ABSTRACT The Deep Ecliptic Survey (DES) was a survey project that discovered hundreds of Kuiper Belt objects from 1998 to 2005. Extensive follow-up observations of these bodies has yielded 304 objects with well-determined orbits and dynamical classifications into one of several categories: Classical, Scattered, Centaur, or 16 mean-motion resonances with Neptune. The DES search fields are well documented, enabling us to calculate the probability on each frame of detecting an object with its particular orbital parameters and absolute magnitude at a randomized point in its orbit. The detection probabilities range from a maximum of 0.32 for the 3:2 resonant object 2002 GF32 −7 to a minimum of 1.5 × 10 for the faint Scattered object 2001 FU185. By grouping individual objects together by dynamical classes, we can estimate the distributions of four parameters that define each class: semimajor axis, eccentricity, inclination, and object size. The orbital element distributions (a, e, and i) were fit to the largest three classes (Classical, 3:2, and Scattered) using a maximum likelihood fit. Using the absolute magnitude (H magnitude) as a proxy for the object size, we fit a power law to the number of objects versus H magnitude for eight classes with at least five detected members (246 objects). The Classical objects are best fit with a power-law slope of α = 1.02 ± 0.01 (observed from 5 H 7.2). Six other dynamical classes (Scattered plus five resonances) have consistent magnitude distribution slopes with the Classicals, provided that the absolute number of objects is scaled. Scattered objects are somewhat more numerous than Classical objects, while there are only a quarter as many 3:2 objects as Classicals. The exception to the power law relation is the Centaurs, which are non-resonant objects with perihelia closer than Neptune and therefore brighter and detectable at smaller sizes. Centaurs were observed from 7.5 1. INTRODUCTION the solar system. Scattering and migration of the giant planets should leave directly observable signatures in the dynamical Studies of the Kuiper Belt have progressed in the last two structure and numbers of small bodies remaining in the asteroid decades from finding individual objects (Jewitt & Luu 1993) and Kuiper Belts and the outer satellites (Levison et al. 2008; to estimating the total number of objects in the outer solar Morbidelli et al. 2009; Bottke et al. 2012). Some models, such system through observations and statistics (e.g., this work; as the Nice Model (Gomes et al. 2005; Morbidelli et al. 2005; Petit et al. 2011; Gladman et al. 2012). The Kuiper Belt itself Tsiganis et al. 2005) and the “Grand Tack” (Walsh et al. 2012), has been subdivided into many distinct dynamical classes of invoke large-scale, abrupt motion by the giant planets. Other objects, including objects in mean-motion resonances with models assume that giant planet migration proceeded more Neptune, high and low inclination populations (Hot versus smoothly, and predict different ratios of objects captured into Cold; e.g., Morbidelli et al. 2008), and various groupings of mean-motion resonances with Neptune (Malhotra 1993; Hahn Scattered objects that have undergone dynamical interactions & Malhotra 2005). with Neptune. The relative numbers of objects in all of these Most models to date can successfully account for some—but populations offer one of the best direct observational constraints not all—features of the Kuiper Belt, such as reproducing the on different dynamical models of solar system formation and Cold but not the Hot Classical populations (smooth migra- evolution. tion), or leaving particular resonances over- or under-populated A recent renaissance in modeling has resulted in a number relative to the apparent populations (most models). Observa- of mostly successful attempts to reproduce the architecture of tional constraints can be used to distinguish between models and resolve fundamental questions, such as how Neptune mi- 8 Deceased. grated outward: smoothly on a nearly circular orbit (e.g., Hahn 1 The Astronomical Journal, 148:55 (17pp), 2014 September Adams et al. & Malhotra 2005), through dynamic scattering with potential field was within ±6◦.5 of the ecliptic and was selected to have at high eccentricity which was later damped down (e.g., Levison least 35 USNO-B astrometric reference stars on each chip. The et al. 2008; Walsh et al. 2012), or some other mechanism. This fields were also chosen to exclude stars brighter than magnitude fundamental debate is evident in recent meetings, where one V = 9.5 (Millis et al. 2002). Typically two 300 s exposures presentation concluded that chaotic capture models are more would be taken on the same night a few hours apart, with a likely based on colors of resonant Kuiper Belt objects, or KBOs third frame on another night during the same observing run. (Sheppard 2012), while another (Noll et al. 2012) claimed that This resulted in orbital arcs of at least 24 hr (required for the large binary fraction in the Cold population favors smooth a designation by the MPC). Most of the fields imaged were migration, since the mechanism of Levison et al. (2008) to push distinct, but sometimes frames of the same field taken several out the Cold population would not preserve as many binary ob- years apart were counted as new search fields, since by that time jects. However, it is also possible to have models with elements a new set of objects would have moved into the field area. In of both: a Nice-like initial migration that proceeds smoothly 2005, the DES stopped surveying new fields, although recovery into a pre-existing Cold population of objects that formed in efforts to improve orbits and classifications of known objects situ would not be ruled out by the binary fraction (Parker & have continued. Kavelaars 2010). In order to obtain as uniform a sample as possible, for this Until recently, comparisons to observations have been disad- work several selection criteria were applied to both the search vantaged because the known sample (over 1600 objects) is an fields and the objects discovered by the DES, as described below. inherently biased population, due to the difficulty of obtaining observations of faint, distant objects. What is needed is a de- 2.1. Search Frame Selection biased population of objects, that is to say, a relatively large sample of objects discovered by a wide-field, all-sky survey that The discovery phase of the DES ended on 2005 May 11. A have had the sources of bias accounted for and removed. Two few KBOs were discovered during subsequent recovery efforts, surveys with large discovery samples have recently begun to but are not included in our sample since they were found under provide such results. The Canada-France Ecliptic Plane Survey different search conditions. (For objects that were discovered (CFEPS), with a sample of 169 objects with high-quality orbits, during the main survey, we used the latest available orbital reported de-biased population estimates in Petit et al. (2011) information to assign dynamical classes and calculate discovery and Gladman et al. (2012). probabilities.) Here we report the results of the Deep Ecliptic Survey (DES; A total of 62,392 individual frames are in the DES database, Millis et al. 2002; Elliot et al. 2005), the largest single survey to corresponding to one of the eight CCD detectors on the Mosaic date, which has discovered about 500 objects with provisional camera. Of these, 14,440 were not suitable for de-biasing designations by the Minor Planet Center (MPC). Of these ob- because these frames were specifically targeted at the recovery jects, 304 have high-quality orbits suitable for dynamical clas- of particular objects. Targeted recovery frames were not held sifications. The de-biased class populations presented in this to the same criteria for the number and brightness of stars, and have different statistical properties. An additional 1312 frames paper offer a completely independent check on the CFEPS re- ◦ sults, on both the absolute and relative numbers of objects, since with solar elongation angle less than 140 were eliminated, the DES uses a different set of objects and discovery fields and because they have a different search efficiency than the bulk of a new approach to de-biasing. The DES sample also more than the frames, which were taken near opposition. (Over twice as doubles the number of de-biased objects available, particularly many objects were found on high solar elongation frames than in the small-number classes of high-order resonances and Scat- on low-elongation frames.) tered objects. Finally, we hope that the approach described in The 46,640 individual frames left were then paired up with this paper will be useful for obtaining uniformly de-biased re- matching frames. Each field was observed at least three times: sults of the expected thousands or tens of thousands of objects a pair of images taken on the same night a few hours apart, found by future large-scale surveys, such as the Large Synoptic and one image on another night. Sometimes additional frames Survey Telescope (LSST; Ivezic et al. 2008) and Pan-STARRS were taken, for instance if weather interrupted the original (Grav et al. 2011). observations. Only the pair of frames on the same night was In Section 2 we describe the data and briefly review the choice used to search for objects. There were 18,668 pairs of frames that met these criteria. The overlap fraction for each pair is quite of classification scheme. In Section 3 we present the analytical framework for calculating the detection probabilities. Section 4 high, with a median positional offset of 0.026. Each frame is 1065 × 532, and the largest mis-registration was less than 80 contains the numerical implementation of the detection proba- bilities and the methods used to estimate the number of objects in R.A. and less than 90 in declination. in each class. In Section 5 we compare our results to other surveys. In Section 6 we relate our observations to theories of 2.2. Object Selection planetary formation. As of 2012 August, a total of 913 KBOs were listed in 2. OBSERVATIONS the DES internal database. Some objects have been excluded from this analysis because they were found outside the main A full description of the DES can be found in Millis et al. survey, during recovery or other follow-up observations (101 (2002) and Elliot et al. (2005, henceforth E05), but the relevant objects). We also exclude eight objects found on low solar details are summarized here. The DES was a National Optical elongation fields (<140◦). Several additional objects, typically Astronomical Observatory (NOAO) survey. Two 4 m telescopes those that were lost shortly after discovery, are excluded because (the Mayall at Kitt Peak in Arizona and the Blanco at Cerro of missing information: objects lack discovery field coordinates Tololo in Chile) were used, with identical Mosaic cameras with (36 objects), discovery distances (104 objects), or discovery eight SITe CCD chips, each chip measuring 4096 × 2048 pixels. magnitudes (1 object). This leaves 663 objects discovered during Each mosaic image covers about 0.35 deg2 on the sky. Each the main DES survey. 2 The Astronomical Journal, 148:55 (17pp), 2014 September Adams et al. Of these objects, 478 objects had sufficient observations post- An excitation statistic, s = e2 +sin(i)2, is applied to the discovery to receive preliminary designations from the MPC remaining objects, which are sent to two additional classes: (72%). Dynamical classification typically requires additional (4) Scattered Extended objects have s 0.25, and (5) Classical observations to achieve low enough errors, and 316 objects objects have s<0.25. Note that the exact criteria used can (48%) have been classified (see Section 2.3). Almost all (304) have an important impact on class membership: for instance, of these classifications are considered to be secure (quality 3 what might in other works be considered Hot Classical objects as defined in E05). For our analyses, we used the 304 securely (low e, large i) with i>15◦ would be classified as Scattered classified DES objects. Extended under this scheme. We note that there is a potential for recovery bias in the For the purposes of this paper, we have re-grouped a few particular subset of objects that were both discovered and objects in order to be more broadly consistent with the classes classified by the DES. One could imagine that objects in unusual used by CFEPS and elsewhere. First, we restrict the population orbits could be missed more frequently than objects that were of Centaurs to only those with semimajor axes less than closer to their predicted location when follow up observations Neptune’s (resulting in a population of smaller objects than were taken a few months later. However, we do not think that found in any other class, due to the magnitude bias toward this recovery bias introduces a major problem for the current finding closer objects). Second, we group all of the dynamically work for the following reasons. (1) The single biggest reason excited objects with aNep a 80 AU into a single Scattered for failure to recover an object was insufficient telescope time class, which contains objects from the Centaur, Scattered- to re-observe an object within three months of discovery, after Near, and Scattered-Extended classes. Ten objects with a = which observations were typically not attempted; this affected 82–739 AU are not considered in this analysis, because they are objects in all orbits alike. It was more likely that an object too sparse in parameter space to adequately define a dynamical was lost because we did not have enough clear weather nights class; for a brief discussion of such rare objects see Section 5.3. than because we looked in the wrong place and did not find it. (2) When the estimated orbital parameters for the lost objects 3. ANALYTICAL FRAMEWORK were compared to those that were ultimately designated and classified, Elliot et al. (2005) found that the only main difference We now consider how we can use the sample of KBOs was that the lost objects were fainter (and thus would be harder discovered by the DES to learn about the entire population of to recover on an average frame). In this work we find that the KBOs. Dynamical classes for KBOs have been defined with rate of recovery was twice as high for objects with brighter the underlying assumption that the members of each class discovery magnitudes (Mdisc 23) than for fainter objects, and have experienced a common formation process and common we only do analyses using the better-recovered subset of brighter dynamical evolution. Hence, we shall consider the objects in objects. Future work is planned to deal with the fainter objects, each dynamical class separately. Within a dynamical class, we including the improved magnitude calibration, and we will also assume that the DES discoveries give us a sampling of the more fully explore the recovery bias at that time. greater population of the objects, which can be characterized by four distribution functions: one each for the H magnitude, 2.3. Dynamical Classification semimajor axis a, eccentricity e, and inclination i. In Section 3.1, we develop the equations to compute the Over the years, several different schema for dynamically probability of detection by the DES for each object, as a function classifying KBOs have been proposed (e.g., Elliot et al. 2005; of H, a, e, and i. In Section 3.2, we present the distribution Lykawka & Mukai 2007; Gladman et al. 2008). In this paper, we functions that we assume are followed by the general population use the current DES classification scheme, a modified version of the dynamical class for the same four parameters. Finally, in of the scheme laid out in E05 which incorporates the Scattering Section 3.3 we develop a maximum likelihood (ML) model to concept from Gladman et al. (2008). Since differences in the fit the data. From this fit, we determine, with error bars, the precise definitions of classes can affect the membership of parameters of the four distribution functions and the number objects and consequently the parameters derived for the class, of objects in each dynamical class, within the range of orbital we describe our classification scheme in full below. parameters probed by the DES discoveries. We have chosen To test for the dynamical class of an object, we use both to apply the ML technique directly, rather than applying least- the current observed orbital parameters as well as a 10 Myr squares fitting to binned data, since the objects per bin are sparse integration of the orbit forward in time, and two additional in the four-dimensional H, a, e, and i space (see chapter 10 of 10 Myr integrations for the ±3σ values for a, e, and i. An object Bevington & Robinson 2003). is said to be securely classified (quality 3) only if all three integrations agree on the same dynamical class. (In quality 3.1. Probabilities of Detection 2 objects, the nominal classification only agrees with one of the +3σ or −3σ integrations, while in quality 1 it agrees with To compute the probability that an object would have been neither.) detected by the DES at some time during the survey, we extend An object is tested for membership in each of the following the methods described in E05 to include observational biases dynamical classes, in order. (1) Resonant objects occupy a introduced by an object’s H, a, e, and i values. We refer all mean-motion resonance with Neptune (resonances up to ninth inclinations to the Kuiper Belt plane (KBP), as derived in E05 order are identified from libration in the resonant angle). (2) (i = 1◦.74 ± 0◦.23, Ω = 99◦.2 ± 6◦.6), which is consistent Centaurs have osculating perihelia that reach values less than with the invariable plane of the solar system at the 1σ level. the osculating semimajor axis of Neptune. (3) Scattered Near Following the approach used in E05, we consider a set of NO objects have a values that vary by more than an arbitrary amount, objects (j = 1,...,NO ) discovered in a set of NF search fields Δa 0.02, over the 10 Myr integration (where Δa = (amax − (k = 1,...,NF ). For the jth object discovered by the survey, amin)/amean). This is along the same lines as the Scattering class we denote its H magnitude, semimajor axis, eccentricity, and of Gladman et al. (2008), which uses amax − amin > 1.5AU. inclination by the symbols Hj, aj, ej, and ij. We assume that 3 The Astronomical Journal, 148:55 (17pp), 2014 September Adams et al. the H magnitude of an object is not correlated with its orbital To find the amount of time that the jth object spends at elements, and that angles that describe the orientation of the orbit a given heliocentric distance, R, it is possible to directly in three-dimensional space average out over time. (For resonant numerically integrate a two-body orbit and then calculate objects, this assumption is not quite true, and is addressed in a the amount of time spent in each distance bin. However, separate bias factor discussed below.) this approach is far too slow for an equation that is called Each of the NF fields of the survey data set is characterized frequently, and there are no good analytical approximations by (1) the magnitude for which the detection efficiency has that work near perihelion and aphelion, where the proba- dropped to 1/2, which we denote by m1/2,k, (2) the Kuiper Belt bility of finding an object peaks sharply. We thus calcu- latitude and orientation, or tilt, of the frame (βK and θK ), used lated a grid of numerical integrations using e = 0.001 and to calculate which orbits pass through a field, (3) the Kuiper e = 0.01–0.99 (Δe = 0.01). For each eccentricity, we fit a Belt longitude of the frame (λK ), and (4) the solid angle, Ωk, piecewise function composed of five line segments (f1 to f5) subtended by the search field. Below we describe how biases at predefined break points (p1 to p4). The coordinates for the related to all four of these characteristics have been removed. slopes and offsets and break points of each line are stored 1. We assume that all search fields have the same maximum for each eccentricity in a python executable file, and may be quickly called during numeric integrations. The function, efficiency at bright magnitudes, max, and the same char- acteristic range, σ , over which the detection efficiency which has been normalized over the allowed R range, is m scaled to the appropriate a value when called. Tests indi- drops from max to 0; as in E05, we have fixed the values = = cate that the results are within ±15% of the values obtained of σm 0.58 and max 1. For the kth search field, we = employ the functional form used by Trujillo et al. (2001) using numeric integration (for objects with a< 80 AU, and in E05 (their Equation (19)) for the detection efficiency which is the cutoff in this paper for the Scattered class). ⎧ as a function of magnitude, (m, m1/2,k): ⎪0 R