The Astrophysical Journal, 597:1211–1236, 2003 November 10 E # 2003. The American Astronomical Society. All rights reserved. Printed in U.S.A.
RESONANT STRUCTURE IN THE KUIPER DISK: AN ASYMMETRIC PLUTINO DISK Elizabeth K. Holmes National Research Council Resident Research Associate, Jet Propulsion Laboratory, California Institute of Technology, MS 169-506, 4800 Oak Grove Drive, Pasadena, CA 91109; [email protected] Stanley F. Dermott and Bo A˚ . S. Gustafson Department of Astronomy, University of Florida, P.O. Box 112055, 211 Bryant Space Science Center, Gainesville, FL 32611; [email protected]fl.edu, [email protected]fl.edu and Keith Grogan Jet Propulsion Laboratory, California Institute of Technology, MS 169-506, 4800 Oak Grove Drive, Pasadena, CA 91109; [email protected] Received 2002 August 29; accepted 2003 July 11
ABSTRACT In order to develop a dynamical model of the Kuiper disk, we run numerical integrations of particles originating from source bodies trapped in the 3 : 2 external mean motion resonance with Neptune to determine what percentage of particles remain in the resonance for a variety of particle and source body sizes. The dynamical evolution of the particles is followed from source to sink with Poynting-Robertson light drag, solar wind drag, radiation pressure, the Lorentz force, neutral interstellar gas drag, and the effects of plane- tary gravitational perturbations included. We find that the number of particles in the 3 : 2 resonance increases with decreasing (i.e., increasing particle size) for the cases in which the initial source bodies are small ( 10 km in diameter) and that the percentage of particles in resonance is not significantly changed by either the addition of the Lorentz force, as long as the potential of the particles is small ( 5 V), or the effect of neutral interstellar gas drag. The brightness of the entire Kuiper disk is calculated using a model composed of 500 lm diameter particles and fits well with upper limits to the Kuiper disk brightness and previous estimates. A disk with a size-frequency distribution weighted toward large particles, which are more likely to remain in resonance, may have a stronger, more easily identifiable resonant signature than a disk composed of small particles. Subject headings: Kuiper Belt — solar system: general On-line material: color figures
1. INTRODUCTION initiate a collisional cascade that would produce a Kuiper dust disk. In fact, the Kuiper Belt is the region of our solar In 1943 and 1951, respectively, Edgeworth and Kuiper system that is most analogous to the planetary debris disks independently speculated that planetary material may lie in we see around other stars such as Vega, Pic, Fomalhaut, our solar system beyond the orbit of Neptune (Jewitt 1999; and Eri (Backman & Paresce 1993). In this paper we dis- Edgeworth 1943; Kuiper 1951). In 1992, Jewitt and Luu dis- cuss some of the dynamical properties of the Kuiper Belt covered 1992 QB1, the first object in what came to be known and the resonant structure of the Kuiper disk. as the Kuiper Belt or the Edgeworth-Kuiper Belt (Jewitt & A Kuiper dust disk will have a resonant structure, with Luu 1993). As of 2003 March, more than 650 Kuiper Belt two concentrations in brightness along the ecliptic longitude objects (KBOs) have been discovered (Trujillo & Brown (Dermott, Malhotra, & Murray 1988a; Liou & Zook 1999) 2002; Parker 20031), which translates into an estimated arising because an estimated 6% of the Kuiper Belt objects 1 105 KBOs with diameters greater than 100 km (Jewitt are in the 3 : 2 mean motion resonance with Neptune 1999) orbiting with semimajor axes between 40 and 200 AU (Trujillo & Brown 2002). Ultimately we want to develop a (Trujillo & Brown 2002). The Kuiper Belt was formed either complete dynamical model of the Kuiper disk. The first step because planet formation at those heliocentric distances did in this process is to run numerical integrations of particles not occur as a result of the low density of planetesimals and originating from source bodies trapped in the 3 : 2 resonance the long encounter timescales between objects and/or and to determine what percentage of particles remain in the because planetesimals were displaced into the region as a resonance for a variety of particle and source body sizes. result of perturbing encounters with young Uranus and Because the particle properties (shape, composition, etc.) Neptune (Backman & Paresce 1993). In the asteroid belt, are not well known, the variable we use to denote particle collisions between asteroids supply dust particles to the zodiacal cloud, the Sun’s inner dust disk. By comparison, it size is , the magnitude of the ratio of the force of radiation pressure to the force of gravity. The dynamical evolution of has been postulated that collisions between KBOs could the particles is followed from source to sink with Poynting- Robertson light drag (PR drag), solar wind drag, radiation pressure, the Lorentz force, neutral interstellar gas drag, 1 The Distant EKOs Electronic Newsletter is available at and the effects of planetary gravitational perturbations http://www.boulder.swri.edu/ekonews. included. 1211 1212 HOLMES ET AL. Vol. 597
Previous studies of particles in the Kuiper disk (Liou & resonance with Neptune at a semimajor axis of 39.4 AU, Zook 1999; Moro-Martı´n & Malhotra 2002) are different and KBOs that reside in the 2 : 1 mean motion resonance from the study presented in this paper in three ways. Liou & with Neptune (a 47:7 AU). The 4 : 3 and 5 : 3 reson- Zook (1999) considered particles that originate outside of ances are also somewhat populated. We focus here on the 3 : 2 resonance and, by means of PR drag, can later the 3 : 2 resonance. Pluto also resides in this resonance become trapped in resonance. However, it may be more (Malhotra & Williams 1998) and, like the Plutinos, realistic to consider particles that originate from source avoids close encounters with Neptune by being an bodies that are in the 3 : 2 resonance. This is especially aphelion librator. Approximately 15% of the observed important for large particles since they may suffer collisions KBOs are resonant KBOs. However, once the current that will quickly break them down into smaller particles observations are debiased, Plutinos (and KBOs in the before they can evolve into the 3 : 2 resonance. (Moro- other mean motion resonance with Neptune) probably Martı´n & Malhotra [2002] consider a range of initial semi- compose 6% of the total KBO population (Trujillo & major axes: 35 AU a 50 AU.) Secondly, Liou & Zook Brown 2002), or on the order of 6000 objects with (1999) and Moro-Martı´n & Malhotra (2002) have not diameters larger than 100 km. The eccentricities included the Lorentz force and neutral interstellar gas drag (0:1 e 0:34) and inclinations (10 I 20 ) of the in their numerical integrations. Finally, we have included Plutinos bracket the values of Pluto (e ¼ 0:25, I ¼ 17 ). particles with 0 <0:13, while Liou & Zook (1999) con- As KBOs continue to be found, the existence of a Kuiper sidered particles with 0:05 0:4 and Moro-Martı´n& dust disk having an equilibrium grain population orbiting Malhotra (2002) considered particles with 0:01 0:4. at approximately 30–50 AU has been postulated (Backman, Their particle distribution is weighted more toward smaller Dasgupta, & Stencel 1995). As such a disk has not yet been particles, yet they do not include the Lorentz force, which is directly observed, many theoretical questions remain about more important for smaller particles. the size-frequency distribution of the particles in the disk, We find that, for the cases in which the initial source the evolution of the disk, and subsequent implications for bodies are small (i.e., 10 km in diameter), the number of the detectability of the Kuiper disk from existing or future particles in the 3 : 2 mean motion resonance increases with observations. Such a disk could be generated in three ways: decreasing (i.e., increasing particle size). Consequently, a (1) by collisions between KBOs, (2) by impacts of interstel- size-frequency distribution for a Plutino disk must be lar dust grains onto KBOs, and (3) by mutual collisions weighted toward larger particles (i.e., those with smaller between dust particles. Models of collisions between KBOs -values). In addition, as long as the potential, U, of the par- by Stern (1996) estimate a dust production rate in the ticles is small (U 5 V), the Lorentz force does not seem to Kuiper Belt of between 9:5 108 and 3:2 1011 gs 1. These prevent the particles from remaining in resonance. predictions are supported by observational evidence: the Similarly, the addition of the effect of neutral interstellar gas Pioneer 10 and Pioneer 11 spacecrafts have detected the sig- drag does not significantly change the percentage of natures of Kuiper dust grains in the outer solar system particles remaining in resonance. (Landgraf et al. 2002). The dust detectors on board the Pioneer spacecraft registered detections due to particles in the size range of 10–6000 lm, too big to be interstellar 2. KUIPER BELT AND DISK grains. The Pioneer 10 spacecraft detected a nearly constant KBOs are divided into three dynamical classes: classical flux of dust particles outside of Jupiter’s orbit. Inside of KBOs, scattered KBOs (SKBOs), and resonant KBOs. The Saturn’s orbit, the dust particles are attributable to a combi- orbits of the classical objects, which account for approxi- nation of KBOs, short-period Oort cloud comets, and mately two-thirds of the known KBOs, have semimajor axes short-period Jupiter-family comets. However, to match the 41 AU a 47 AU and seem unassociated with resonan- observations, Landgraf et al. (2002) conclude that the flux ces. They have moderate eccentricities (e 0:1) and their of dust particles outside of Saturn’s orbit must be domi- inclinations lie in the range 0 I 32 (Jewitt 1999). nated by grains originating from the Kuiper disk. A Kuiper Assuming that there are 105 KBOs with diameters larger Belt dust production rate of approximately 5 107 gs 1 is than 100 km, there are on the order of 50,000 classical KBOs needed to account for the amount of dust found by the in that size range (Jewitt 1999; Trujillo & Brown 2002). Pioneer spacecrafts outside the orbit of Saturn. Yamamoto Currently, approximately 10% of the observed KBOs are & Mukai (1998) estimate a total dust production rate of dust SKBOs (Trujillo & Brown 2002; Trujillo, Jewitt, & Luu grains due to the impacts of interstellar grains with KBOs of 5 7 1 2000). The first member to be discovered, 1996 TL66, has a 3:7 10 or 3:1 10 gs depending on whether the semimajor axis of 85 AU, an eccentricity of 0.59, and an surfaces of the KBOs are hard and icy or are covered by a inclination of 24 (Luu et al. 1997). As a whole, the SKBOs layer of icy particles, respectively. Finally, Liou, Zook, & are characterized by large eccentricities, large semimajor Dermott (1996) estimate that the average time for two 9 lm axes (>50 AU), and perihelia, q, near 35 AU because they diameter Kuiper grains to have a mutual collision is were scattered by Neptune. The total number of SKBOs 3:11 107 yr. In addition, for Kuiper grains starting out at with diameters larger than 100 km is estimated to be on the 40 AU, Pepin, Palma, & Schlutter (2001) report a mean col- order of 50,000 (Trujillo & Brown 2002). The SKBOs are lisional survival time of between 2 106 and 10 106 yr for thought to have been scattered by Neptune into a scattered grains with diameters ranging from 10 to 10,000 lm. In fact, disk (Trujillo et al. 2000). Finally, it should be noted that their 10 lm diameter grains succeed in penetrating to 1 AU. the observational techniques are biased toward lower (In comparison, the PR drag lifetime for a 10 lm grain at inclinations, so the reported range of inclinations for these 40 AU with density 2 g cm 2 is 4 107 yr for the geom- dynamical classes will be less than the actual range. etrical optics case. The PR drag lifetime for a 10,000 lm The resonant KBOs are primarily composed of the grain is 4 1010 yr.) Since we only follow our numerical Plutinos, objects that reside in the 3 : 2 mean motion integrations for 2:5 105 yr, we neglect mutual collisions. No. 2, 2003 RESONANT STRUCTURE IN KUIPER DISK 1213
3. RESONANT STRUCTURE: MEAN line-of-sight integrator allows the model to be viewed as it MOTION RESONANCES would appear in the 70 lm wave band in a coordinate frame that corotates with Neptune. Note the presence of the two If a Kuiper disk exists, it would most likely have a reso- lobes. A disk of particles trapped in a p þ 1 : p exterior mean nant structure, with two concentrations in brightness along motion resonance will have p lobes, which are located at the the ecliptic longitude (see Fig. 1). The structure arises from particles’ pericenters. In the rotating frame, the lobes are dust particles originating from the Plutinos, the KBOs caused by the fact that the angular velocity of the particles trapped in the 3 : 2 mean motion resonance with Neptune, at pericenter becomes comparable to the angular velocity meaning that the particles orbit the Sun twice for every three of Neptune (Murray & Dermott 1999). The position of times Neptune completes an orbit. We build a simple toy Neptune would lie 90 from the two lobes, since the Plutinos model of a Plutino disk, which is shown in Figure 1. The are aphelion librators. Essentially, Figure 1 represents a particles in the model are all one size: they are assumed to be snapshot of how the resonant cloud would appear, since the large with ¼ jjFrad = Fgrav 0. The particles’ orbits are lobes corotate with Neptune; the positions of the lobes will distributed through a three-dimensional array of cells, and a change with time with respect to a non-Neptunian reference
Fig. 1.—Toy model of a Plutino disk. The particles in the model are assumed to be large ( 300 lm in diameter) with 0. The particles’ orbits are distributed through a three-dimensional array of cells, and a line-of-sight integrator allows the model to be viewed as it would appear in the 70 lm wave band in a coordinate frame that corotates with Neptune. Note the presence of the two lobes. A disk of particles trapped in a p þ 1 : p external mean motion resonance will have p lobes, which are located at the particle’s pericenters. The position of Neptune would lie 90 from the two lobes, since the Plutinos are aphelion librators. The intensity is in arbitrary units of brightness, since the total brightness of the disk is not well determined. [See the electronic edition of the Journal for a color version of this figure.] 1214 HOLMES ET AL. Vol. 597 point. For example, in 2001 January Neptune’s geocentric shown in Figure 1 but with a nonzero libration width. The ecliptic longitude is 126 (USNO and RGO 2001) and the lobes are still visible but have a smoothed appearance. geocentric ecliptic longitudes of the lobes would be 36 and Ultimately, the resonant structure in the Kuiper disk 216 . In 2001 December Neptune’s geocentric ecliptic longi- may be detectable using data from COBE or the Space tude is 128 and, correspondingly, the geocentric ecliptic Infrared Telescope Facility (SIRTF). The main impedi- longitudes of the lobes would be 38 and 218 , demonstrat- ment to detecting the Plutino ring is that, since it is ing that this signal is time varying. Since the Plutinos popu- located at approximately 39.4 AU, it will be very faint. late the densest and closest resonance to the inner edge of Many brighter layers of emission will need to be sub- the Kuiper Belt, dust from these bodies would generate the tracted from the data before the faint residuals can be brightest emission. This Plutino resonant ring is similar in studied for traces of the Plutino resonant structure. More nature to the resonant ring of dust particles corotating with promising may be the idea of detecting resonant struc- the Earth (Jackson & Zook 1992; Dermott et al. 1994b; tures in planetary debris disks. In Figure 3 we take the Jayaraman 1995), except that the lobes are symmetrical with model shown in Figure 1 and show it as it would appear respect to Neptune. There is no trailing cloud because the to SIRTF in the MIPS 70 lm wave band if the disk were drag rate on the particles at 40 AU is very small (Dermott located at 2 pc. Even with a pixel size of 100, the lobes are et al. 1998). Shown in Figure 2 is a similar model to that still visible.
Fig. 2.—Toy model of a Kuiper disk from asteroidal particles (with a distribution of eccentricities and inclinations) in the 3 : 2 mean motion resonance with Neptune with a nonzero libration amplitude. This model is similar to that in Fig. 1 but with a nonzero libration width, which is represented by a Gaussian distribution of width 30 . The two lobes are still visible but have a smoothed appearance. The intensity is in arbitrary units of brightness. [See the electronic edition of the Journal for a color version of this figure.] No. 2, 2003 RESONANT STRUCTURE IN KUIPER DISK 1215
sion. For a more complete discussion we refer the reader to Burns, Lamy, & Soter (1979) and Gustafson (1994). The orbit of a particle or a planet can be described by a set of six osculating elements: the semimajor axis (a), eccen- tricity (e), inclination of the orbit to a reference plane (I ), longitude of ascending node ( ), longitude of pericenter (!~), and the mean longitude ( ). The forces below are calculated for the case of the solar system, but it is trivial to derive expressions for an exosolar system case. In general, forces on a dust particle can be separated into two categories: gravitational and nongravitational forces. The gravitational forces will be discussed in more detail in a subsequent section on mean motion resonances.
4.1. Radiation Pressure Particles are affected by photons streaming out from the Sun. These photons act on the particles, producing a repul- sive force known as the radiation pressure force, Frad.A dimensionless quantity, , is used to describe radiation pressure. It is defined as the ratio of the radiation pressure force to the force of gravity: ¼ jjFrad = Fgrav . The force of 2 gravity can be written as Fgrav ¼ ðGMm=r Þ^r, where G ¼ 6:67 10 11 Nm2 kg 2 is the gravitational constant, M is the mass of the Sun (M ¼ 1:99 1030 kg) or, in gen- Fig. 3.—Toy model of the Plutino disk shown in Fig. 1 as it would eral, the central star, r is the heliocentric distance, and m is appear to SIRTF in the MIPS 70 lm wave band if the disk were located at 2 pc. This image was constructed from 12 footprints. Each pixel is 100 wide. the mass of an individual dust particle. If we combine the The image, which has a linear scale and arbitrary intensity units, was pro- forces of gravity and radiation pressure, the resulting total duced using HiRes, a product of the Jet Propulsion Laboratory’s Center force on the particle is for Long Wavelength Astrophysics. The two lobes are still visible. [See the electronic edition of the Journal for a color version of this figure.] Gð1 ÞMm F ¼ ^r ; ð1Þ r2 The models shown in Figures 1, 2, and 3 are very simplis- which means that when a particle is acted upon by radiation tic and are meant for illustrative purposes only. We have pressure, the gravitational pull on the particle is equivalent included a heuristic distribution of particle eccentricities to that of the Sun if its mass were reduced by a factor of and inclinations in the models and have not included any (1 ) (Burns et al. 1979; Gustafson 1994). nongravitational forces, which need to be considered for particles less than 500 lm in diameter. The goal of the fol- 4.2. Poynting-Robertson Light Drag lowing sections is to calculate, first analytically and then When a particle absorbs incident radiation from the Sun, numerically, the percentage of particles that stay in or it reradiates (or scatters) it isotropically in the particle’s rest become trapped in the 3 : 2 mean motion resonance for a frame. However, since the particle is moving with respect to given set of initial conditions to gain a more complete the Sun, according to special relativity the particle does not understanding of resonant structure in the Kuiper disk. We reradiate (or scatter) the energy flux isotropically in the rest begin by discussing the forces that must be considered. frame of the Sun. The particle reradiates (or scatters) more energy in the direction of its velocity than in the opposite direction, causing a loss of momentum. Since the mass of 4. FORCES ON A DUST PARTICLE the particle is constant, this loss of momentum translates A body, such as an asteroid or dust particle, orbiting a into a deceleration of the particle, resulting in a net drag on star possessing a planetary system is subject to three types the particle. The drag effect is referred to as PR drag of gravitational perturbations due to planets in the system: (Poynting 1903; Robertson 1937; Burns et al. 1979). The long-period or secular perturbations, resonant perturba- loss of momentum causes the particle to move to a lower tions, and gravitational scattering by the planets. In addi- orbit, resulting in the particle slowly spiraling into the Sun tion, for dust particles 500–1000 lm in diameter, over a period of time called the PR drag lifetime (Burns nongravitational forces must also be considered. In this sec- et al. 1979). To the first order in v=c, the total radiation force tion we give a brief overview of the nongravitational forces (radiation pressure force combined with PR drag) acting on on a dust particle that are considered in our numerical simu- a spherical particle is lations, namely, radiation pressure, PR drag, solar wind "# ! drag, the Lorentz force, and neutral interstellar gas drag. 2r_ r _ F 1^r ^r h^ ; ð2Þ Different combinations of these forces are included in the grav c c dynamical models and numerical integrations in this paper. Another drag force, the Yarkovsky effect, typically only where the unit vector h^ is normal to ^r in the orbital plane affects large bodies on the order of 1 m in diameter. Since (Burns et al. 1979; Klacka 1992; Gustafson 1994). The first our studies only encompass small micron-sized particles, we term is the radiation pressure force already discussed above, can omit the Yarkovsky effect from the following discus- and both the second and third terms are PR drag. Since PR 1216 HOLMES ET AL. Vol. 597 drag is a radial force, the inclination and the ascending node This is the classical representation of the Lorentz force. The are not affected and the change in longitude of perihelion is second term in the above equation is independent of the par- negligible for small eccentricities. Given below are the ticle’s motion and can be viewed as resulting from an changes in semimajor axis and eccentricity with respect to induced electric field. Gustafson & Misconi (1979) found time where PR ¼ GM =c (Wyatt & Whipple 1950): that although solar wind parameters are variable and uncer- tain over long timescales, this second term is stable in mag- da ðÞ2 þ 3e2 ¼ PR ; ð3Þ nitude and direction. Gustafson & Misconi (1979) and more dt a ðÞ1 e2 3=2 recently Grogan, Dermott, & Gustafson (1996) modeled the de 5 e magnetic field, utilizing the expanding solar corona model PR : 4 ¼ 2 1=2 ð Þ of the solar magnetic field (Parker 1958). In Parker’s model, dt 2a ðÞ1 e2 the gas flowing outward from the solar corona draws out the solar magnetic field lines as the Sun rotates, so that close 4.3. Solar Wind Corpuscular Drag to the Sun the field is approximately radial, but farther from the Sun the field lines form an Archimedean spiral (Grogan Another nongravitational force acting on dust particles is et al. 1996; Gustafson 1994). In the following equations, B0 the force due to solar wind corpuscular drag, Fsw, which is is the magnetic field strength at some reference distance, caused by solar wind protons streaming out from the Sun b (B0 3 nT), ! is the angular velocity of the Sun and impacting the particles in their paths. Fsw is analogous (! 2:7 10 6), r is the radial distance from the Sun, is to the total radiation force given in equation (2). We can the azimuth angle, and h is the heliocentric colatitude angle; define to be equal to ¼jF j=jF j. Taking v to 1 sw sw sw grav sw the solar wind velocity, vsw, is taken to be 400 km s : be the solar wind speed, we can write the solar wind force as "# ! B0 _ Br ¼ ; ð9Þ 2r_ r ^ r2 Fsw ¼ Fgrav sw 1^r ^r h ; ð5Þ 2 vsw vsw B0!b ðr bÞ sin B ¼ 2 ; ð10Þ vswr which has the same form as equation (2) (Gustafson 1994). We can approximate the solar wind drag force as being B ¼ 0 : ð11Þ 1 roughly 3 the value of the total radiation force given in The above equations are derived using the velocity of the equation (2). gas flowing out from the Sun, v, where vr ¼ vsw, v ! r b sin ,andv 0 (Parker 1958). 4.4. Solar Wind Lorentz Force ¼ ð Þ ¼ Interplanetary dust particles are charged and as a result 4.5. Neutral Interstellar Gas Drag are coupled to the solar wind–driven interplanetary mag- netic field. The particles acquire a charge by a variety of This last force, while small in magnitude, acts from a methods: the photoemission of electrons from the absorp- specific direction and could potentially cause noticeable tion of ultraviolet (UV) solar radiation and the accretion of effects on the orbits of dust particles over a large period solar wind protons tends to positively charge the particles, of time. The solar system inhabits a ‘‘ Local Bubble ’’ in while the accretion of solar wind electrons imparts a the local interstellar medium (LISM) consisting of hot negative charge to the dust grains. The equation for the plasma (106 K) of very low density (5 10 3 cm 3) that equilibrium charge state can be written as spans roughly 100 pc. It is thought that this void was created by a supernova explosion (Gehrels & Chen 1993; eðUÞ fe ¼ pðUÞfp þ y ph fph ; ð6Þ Gru¨n et al. 1994). Within the Local Bubble lies the ‘‘ Local Fluff,’’ consisting of several clouds a few parsecs where , , and are the cross sections of the grains e p ph wide having higher densities and cooler temperatures for electrons, protons, and UV photons, respectively, the than those of the Local Bubble, and within the Local f-values are the corresponding fluxes, and y is the photo- Fluff lies our solar system. The cloud in which the Sun electron yield (Gustafson & Misconi 1979). It is predicted resides has a temperature of 104 K and a density of 0.1 that under normal conditions, photoemission dominates H atoms cm 3 (Bertin et al. 1993; Gru¨n et al. 1994). The over the accretion of solar wind particles, giving the dust interaction of the solar wind with the LISM causes a cav- particles a net positive potential, U,of 5 V (Goertz 1989; ity to form, which extends out to approximately 100 AU Gustafson 1994). This potential is independent of distance and is termed the heliosphere. Inside the heliosphere only since solar wind density follows an approximate 1/r2 neutral gas atoms, primarily hydrogen, reach the plane- decrease in density (Leinert & Gru¨n 1990), and it corre- tary system since interstellar ions are removed from the sponds to a charge q 4 Us on a spherical particle where ¼ 0 solar system by the solar wind (Geiss, Gloeckler, & Mall ¼ 8:859 10 12 CV 1 m 1. The force exerted on the 0 1994; Gloeckler et al. 1993; Gru¨n et al. 1994). particle by the interplanetary magnetic field, B, can be The Ulysses spacecraft has detected a stream of inter- written as stellar helium (Witte et al. 1993; Gru¨n et al. 1994), a
FL ¼ qv µ B ; ð7Þ tracer of hydrogen, originating from an ecliptic longitude of 252 and an ecliptic latitude of +2=5 moving with a where v, the velocity with respect to this field, can be broken speed of 26 km s 1 (Witte et al. 1993; Gru¨n et al. 1994), down into components vg, the heliocentric velocity of the which is in good agreement with the direction of motion particle, and vsw, the solar wind velocity: of the local cloud in which the Sun is embedded with respect to the solar system (Bertin et al. 1993). These µ µ FL ¼ qðvg B vsw BÞ : ð8Þ neutral gas atoms can collide with dust grains in the solar No. 2, 2003 RESONANT STRUCTURE IN KUIPER DISK 1217 system, causing a change in the momentum of the dust distance r, and ðvðrÞÞ is defined by particles. While the effects of one collision are negligible, over a period of time repeated collisions have the 1 mp Mt ; ð17Þ potential to alter the orbital elements of the dust particle ðvðrÞÞ significantly (Scherer 2000). For a spherical dust grain with radius s, mass m,and where Mt is the mass of the target particle (Leinert & Gru¨n 1990). The collisional lifetime of a particle is defined as velocity v, the drag force, FH, on the dust particle due to collisions with neutral interstellar hydrogen particles of 1 mass mH moving with a velocity vH can be written as ¼ ; ð18Þ c Cðm; rÞ FH ¼ NmHuu^ ; ð12Þ the reciprocal of the rate of catastrophic collisions (Leinert 3 where u ¼ v vH, u ¼ jjv vH , u^ is the unit vector in the & Gru¨n 1990). Hence, roughly speaking, nðmp; rÞ r , 0:5 3:5 3:5 u-direction, and N, the average number of collisions of hivðrÞ r ,soCðm; rÞ r making c r . The PR 2 2 hydrogen atoms on a dust particle of collisional cross- drag lifetime only goes as PR a r . Since the colli- sectional area, , per unit time, is defined as sional lifetime at heliocentric distance r is a steeper function of r than the PR drag lifetime, we can make the assumption N ¼ nH uCD ; ð13Þ that the collisional lifetime of particles at r is longer than their PR drag lifetime. We can assume that collisions can be 2 where nH is the number density of hydrogen atoms, ¼ s , neglected for small particles at a distance r from the Sun. where s is the radius of the dust particle, and CD is the free However, we know that collisions do become important for molecular drag coefficient due to hydrogen atoms in a particles larger than 100 lm in diameter (Gustafson et al. Maxwellian velocity distribution of temperature TH impact- 1992) and will need to be considered in future work. ing a spherical dust particle (Gustafson 1994; Scherer 2000), where TH 8000 K. Substituting equation (13) into 4.7. Comparisons of Forces equation (12) gives We calculate the magnitude of various forces (in newtons) F ¼ n s2m jjv v 2C u^ ; ð14Þ on dust particles of different -values (assuming that the H H H H D particles are spheres with diameter d ) and present them in where the numerical values of the constants are as follows: Table 1, assuming that a given particle has a density of 3 3 nH ¼ 0:1 H atoms cm , s is a range of sizes, mH ¼ 1:675 ¼ 2:5gcm , a mass of mdust, a semimajor axis a where 10 24 g, the particle velocity, v, is calculated at each time a ¼ r ¼ 39:4 AU, and an eccentricity e ¼ 0. In this case, the 1 step in numerical integrations, vH is 26 km s ,andCD, the heliocentric velocity of the particle is given by 1=2 1 free molecular drag coefficient, is roughly 1–2 (Gustafson v ¼ ½ GM ð1 Þ=a in km s . The mass of the Sun, M , 30 1994). The value of CD varies with the type of collision that is taken to be 1:99 10 kg. The solar wind velocity is 1 9 1 is assumed to occur and also with the particle shape and taken to be 400 km s , and B0 is 3nT¼ 3 10 Nsm 1 orientation. For instance, a CD of 2 corresponds to the case C . The charge, q, on the dust particle is q ¼ 4 0Ud=2, 12 1 1 in which the total momentum of a proton that hits the dust where 0 ¼ 8:86 10 CV m is the permeability of a vacuum. F is the force of gravity. The magnitude of the particle is transferred to the particle (Gustafson 1994). We grav do not want to limit our calculations by making too many force due to radiation pressure, jjFrad , is simply Fgrav . assumptions and have elected to run simulations with The force due to PR drag, FPR, is represented by the second CD ¼ 1and2. and third terms in equation (2). The Lorentz force, FL,is defined in equation (8), and for the above calculations, we chose the azimuth angle, ,tobe0 and the heliocentric 4.6. Collisions colatitude angle, h,tobe20. Finally, FH, the force due to Collisions between dust particles will also affect some dust neutral interstellar gas particles impacting the dust grains, particles. For small particles (d100 lm) in the asteroid belt, is given by equation (14). Additional constants and collisions are not important since their collisional lifetimes assumptions are given in Holmes (2002). are much longer than their PR drag lifetimes. Leinert & Gru¨n (1990) estimate the average impact speed between two TABLE 1 dust particles as Magnitude of Various Forces on 10 and 100 lm Diameter Spheroidal Dust Particles r 0:5 hivðrÞ ¼ v0 ; ð15Þ Sphere Diameter r0 Parameter 1 10 lm 100 lm where v0 ¼ 20 km s at r0 ¼ 1 AU. The rate Cðm; rÞ of catastrophic collisions, collisions that result in the ¼ jjFrad = Fgrav ...... 0.04868 0.00446 fragmentation of both particles, is given by a ¼ r (AU) ...... 39.4 39.4 12 9 Z m dust (kg)...... 1.3 10 1.3 10 1 18 15 Fgrav (N) ...... 5.0 10 5.0 10 Cðm; rÞ¼ nðmp; rÞhivðrÞ dmp ; ð16Þ 19 17 jjFrad (N)...... 2.4 10 2.2 10 m= ðvðrÞÞ 24 22 jjFPR (N)...... 9.3 10 8.7 10 20 19 jjFL (U ¼ 5 V) (N) ...... 2.9 10 2.9 10 where is the cross-sectional area inside the circle where the 19 18 two particles touch, m is the mass of the projectile, nðm ; rÞ jjFL (U ¼ 20 V) (N) ...... 1.2 10 1.2 10 p p jjF (C ¼ 1) (N)...... 1.2 10 23 1.2 10 21 is the number density of particles of mass m at heliocentric H D 1218 HOLMES ET AL. Vol. 597
TABLE 2 denote those orbital elements of the external orbiting body Ratios of the Magnitude of the Forces on a and the unprimed quantities denote those of the interior 10 ( ¼ 0:04868) and 100 lm ( ¼ 0:00446) Dust Particle body. Physically, the resonant argument can be thought of as the angular distance between the particle’s perihelion and Sphere Diameter the longitude of conjunction of the particle with the planet. We will retain the primes in the subsequent discussion so Parameter 10 lm 100 lm that the orbital elements of the planet, in this case Neptune, 2 3 will be denoted as (aN, eN, IN, N, !~N, N) and those of the Frad=Fgrav ¼ ...... 4.9 10 4.5 10 F =F ...... 1.9 10 6 1.7 10 7 Plutinos and dust particles associated with the Plutinos as PR grav 0 0 0 0 0 0 3 5 (a , e , I , , !~ , ). This distinction is important to make FL=Fgrav (U ¼ 5 V)...... 5.9 10 5.9 10 2 4 FL=Fgrav (U ¼ 20 V)...... 2.3 10 2.3 10 since the disturbing functions are different for interior and 6 7 FH=Fgrav (CD ¼ 1)...... 2.5 10 2.5 10 external resonances, and so the derivation of libration width depends on whether the resonance is interior or exterior. Also note that if p ¼ 1, an extra term must be included in the disturbing function. However, since we are dealing with The ratios of the forces in Table 1 are shown in Table 2. the 3 : 2 resonance, a p ¼ 2 resonance, we need not include In the case of the 10 lm diameter Kuiper disk particle, we that term in our derivation. see that F =F is the same order of magnitude as rad grav FL=Fgrav for U ¼ 20 V and almost the same order of mag- nitude for U ¼ 5 V. However, for the 100 lm particle, 5.1. First-Order External Resonances in the Circular F =F is at least an order of magnitude greater than rad grav Restricted Three-Body Problem FL=Fgrav for both potentials, making the Lorentz force less important for larger particles. For both the 10 and 100 lm For exact first-order external resonances jj¼j 1 ¼ q, 3 size particles, F =F is the same order of magnitude as and j , j , and j are zero, giving PR grav 4 5 6 F =F and in both cases the ratio is 4 orders of H grav 0 0 magnitude less than Frad=Fgrav . ’ ¼ðp þ 1Þ p N !~ ; ð23Þ
if we take j3 to be 1. Now we write the general term in the averaged expansion of the disturbing function, with both 5. THE DYNAMICS OF MEAN secular and resonant terms included. Equation (21) MOTION RESONANCES becomes Consider two large objects moving around a large, central GM 0 N 02 0jjj3 mass in circular, coplanar orbits. For now, we assume that hi¼Rres fs;1ð Þe þ fd ð Þe cos ’ : ð24Þ ¼ 0. The outer object (primed) is said to move in aN resonance with the inner object, if their mean motions, n0 The first term on the right-hand side is the secular term, and n, respectively, can be related as which depends on fs;1ð Þ. The second term is the resonant n p þ q term, and it is a function of f ( ). Both f ð Þ and f ( ) are ¼ ; ð19Þ d s;1 d n0 p functions of Laplace coefficients and are defined in Murray & Dermott (1999). The term containing fs;1ð Þ will be even- where p and q are integers, q denoting the order of the tually neglected. For the 3 : 2 external mean motion with resonance. Since the mean motion can be related to the Neptune, the resonant argument is 2 3 semimajor axis by l ¼ GM ¼ n a , the following relation 0 0 can be written: ’ ¼ 3 2 N !~ : ð25Þ 2=3 p þ q For first-order external resonances, j3 ¼ 1 and the a0 ¼ a : ð20Þ p relationship for the libration width is "s ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi Following the derivations given in Dermott et al. (1988a), 16jjC0 1 jjC0 a0 ¼ r e0 1 þ r Weidenshilling & Jackson (1993), Jayaraman (1995), and max 3n0 27j2e03 n0 Murray & Dermott (1999), if we consider only the resonant 1 # terms, the perturbing force on a particle due to a planet can 2 jjC0 r a0 ; ð26Þ be written as a series expansion in the orbital elements of the 9j e0 n0 planet (a, e, I, , !~, ) and the particle (a0, e0, I0, 0, !~0, 0): 1 X R0 S0 a e I a0 e0 I 0 where res ¼ l ðÞ; ; ; ; ; cos ’; ð21Þ GM where l ¼ GM, G is the gravitational constant, and M is the C0 N f ; 27 r ¼ 02 0 d ð Þ ð Þ mass of the planet. The strength of the resonance, S0, the a aNn strength of an individual term in the series, can be written as where M is the mass of Neptune and G is the solar gravita- S0 S0 ; sin I; sin I 0; e; e0 , where a=a0 . The general N ¼ ð Þ ¼ð Þ tional constant. The full derivation of the libration width form of the resonant argument, ’, can be written as using the pendulum model is given in Murray & Dermott 0 0 0 (1999). Figure 4 shows the maximum libration widths, cal- ’ ¼ j1 þ j2 þ j3!~ þ j4!~ þ j5 þ j6 ; ð22Þ culated using equation (26), in a-e space for the 3 : 2 external where j1 ¼ðp þ qÞ, j2 ¼ q, and the primed quantities mean motion resonance with Neptune. No. 2, 2003 RESONANT STRUCTURE IN KUIPER DISK 1219
Cellino et al. (1999), who parameterized the escape velocity in the following simplified way: d Vesc ¼ vesc 1 ; ð30Þ Rs
where d is the diameter of the fragment, Rs is the radius of the parent body, and vesc is the full escape velocity from the parent body given by sffiffiffiffiffiffiffiffiffiffiffiffiffi
2GMs vesc ¼ ; ð31Þ Rs
where Ms is the mass of the parent body. If we assume that the parent body is spherical and take s to be its density, we 1=2 can rewrite equation (31) as vesc ¼ Rs½ ð8 G sÞ=3 . In the case we are considering, that of micron-sized dust particles, d5 R in equation (30) and therefore V v . Deriving Fig. 4.—Maximum libration widths in a-e space for the 3 : 2 external s esc esc mean motion resonance with Neptune. an effective diameter for Eros of 17.9 km and assuming 3 1 ¼ 2:5gcm , we calculate Vesc ¼ 10:6ms for Eros, which is very close to the average guaranteed escape speed of 10.3 m s 1 derived by Scheeres et al. (2003), Yeomans et al. (2000), and Miller et al. (2002). 6. SOURCE BODIES We take Dvt, Dvr, and Dv? to be the dispersions in the Particles are generated in the Kuiper Belt by the breakup tangential, radial, and perpendicular components of of KBOs. If a parent body has orbital elements (as, es, Is, s, the particles’ velocities upon ejection. If we assume that the !~s), then immediately after breakup the collisional by- collisions disperse particles isotropically, the average com- products of that parent body have a range of orbital ponents of the velocity dispersion are equal to each other elements (Da, De, DI, D , D!~). We can determine if particles (Zappala et al. 1984; Cellino et al. 1999): will stay in a given resonance by comparing Da with the 0 V1 libration width, amax, for that particular resonance. If Dvt ffi Dvr ffi Dv? ffi pffiffiffi : ð32Þ 3 Da <