Modeling the 3-D secular planetary three-body problem: Discussion on the outer υ Andromedae planetary system

The three-dimensional secular behavior of a system composed of a central star and two massive planets is modeled semi-analytically in the frame of the general three-body problem. The main dynamical features of the system are presented in geometrical pictures allowing us to investigate a large domain of the phase space of this problem without time-expensive numerical integrations of the equations of motion and without any restriction on the magnitude of the planetary eccentricities, inclinations and mutual distance. Several regimes of motion of the system are observed. With respect to the secular angle Δ?, possible motions are circulations, oscillations (around 0° and 180°), and high-eccentricity/inclination librations in secular resonances. With respect to the arguments of pericenter, ω₁ and ω₂, possible motions are direct circulation and high-inclination libration around ±90° in the Lidov-Kozai resonance. The regions of transition between domains of different regimes of motion are characterized by chaotic behavior. We apply the analysis to the case of the two outer planets of the υ Andromedae system, observed edge-on. The topology of the 3-D phase space of this system is investigated in detail by means of surfaces of section, periodic orbits and dynamical spectra, mapping techniques and numerical simulations. We obtain the general structure of the phase space, and the boundaries of the spatial secular stability. We find that this system is secularly stable in a large domain of eccentricities and inclinations.     

The effects of metallicity and grain growth and settling on the early evolution of gaseous protoplanets

Giant protoplanets formed by gravitational instability in the outer regions of circumstellar disks go through an early phase of quasi-static contraction during which radii are large (∼1 AU) and internal temperatures are low (<2000 K). The main source of opacity in these objects is dust grains. We investigate two problems involving the effect of opacity on the evolution of isolated, non-accreting planets of 3, 5, and 7 M_J. First, we pick three different overall metallicities for the planet and simply scale the opacity accordingly. We show that higher metallicity results in slower contraction as a result of higher opacity. It is found that the pre-collapse time scale is proportional to the metallicity. In this scenario, survival of giant planets formed by gravitational instability is predicted to be more likely around low-metallicity stars, since they evolve to the point of collapse to small size on shorter time scales. But metal-rich planets, as a result of longer contraction times, have the best opportunity to capture planetesimals and form heavy-element cores. Second, we investigate the effects of opacity reduction as a result of grain growth and settling, for the same three planetary masses and for three different values of overall metallicity. When these processes are included, the pre-collapse time scale is found to be of order 1000 years for the three masses, significantly shorter than the time scale calculated without these effects. In this case the time scale is found to be relatively insensitive to planetary mass and composition. However, the effects of planetary rotation and accretion of gas and dust, which could increase the timescale, are not included in the calculation. The short time scale we find would preclude metal enrichment by planetesimal capture, as well as heavy-element core formation, over a large range of planetary masses and metallicities.     

The web of three-planet resonances in the outer Solar System: II. A source of orbital instability for Uranus and Neptune

The motion of the giant planets from Jupiter to Neptune is chaotic with Lyapunov time of approximately 10 Myr. A recent theory explains the presence of this chaos with three-planet mean-motion resonances, i.e. resonances among the orbital periods of at least three planets. We find that the distribution of these resonances with respect to the semi-major axes of all the planets is compatible with orbital instability. In particular, they overlap in a region of 10⁻³ AU with respect to the variation of the semi-major axes of Uranus and Neptune. Fictitious planetary systems with initial conditions in this region can undergo systematic variations of semi-major axes. The true Solar System is marginally in this region, and Uranus and Neptune undergo very slow systematic variations of semi-major axes with speed of order 10⁻⁴ AU/Gyr.     

Formation of gas giant planets: core accretion models with fragmentation and planetary envelope

We have calculated formation of gas giant planets based on the standard core accretion model including effects of fragmentation and planetary envelope. The accretion process is found to proceed as follows. As a result of runaway growth of planetesimals with initial radii of ∼10 km, planetary embryos with a mass of ∼10²⁷ g (∼ Mars mass) are found to form in ∼10⁵ years at Jupiter's position (5.2 AU), assuming a large enough value of the surface density of solid material (25 g/cm²) in the accretion disk at that distance. Strong gravitational perturbations between the runaway planetary embryos and the remaining planetesimals cause the random velocities of the planetesimals to become large enough for collisions between small planetesimals to lead to their catastrophic disruption. This produces a large number of fragments. At the same time, the planetary embryos have envelopes, that reduce energies of fragments by gas drag and capture them. The large radius of the envelope increases the collision rate between them, resulting in rapid growth of the planetary embryos. By the combined effects of fragmentation and planetary envelope, the largest planetary embryo with 21M_⊕ forms at 5.2 AU in 3.8×10⁶ years. The planetary embryo is massive enough to start a rapid gas accretion and forms a gas giant planet.     

Spin rates of small moonlets embedded in planetary rings: I. Three-body calculations

We investigate the spin rates of moonlets embedded in planetary rings, subject to collisions with surrounding small particles, using three-body integrations including friction and spins. All successive impacts of the particle with the moonlet are followed, including a possible sliding phase after the initial inelastic rebounds. Two methods for treating impacts, (1) as instantaneous velocity changes and (2) using an impact force model, are applied after Salo (1995, Icarus 117, 287). Conducting a series of integrations with various initial summed spin velocity of the moonlet and the particle, we determine the equilibrium spin rate for which the averaged torque vanishes. This equilibrium spin rate corresponds to the final spin rate of the moonlet if the moonlet is much larger than the surrounding particles; it also corresponds to the mean spin rate for a ring composed of identical particles. We find that the equilibrium spin rate is enhanced by sliding orbits as compared with the spin rate determined by considering only the first impacts of the particles with the moonlet. If the random velocities of incident particles are small enough, the resulting equilibrium spin rate of the moonlet can be larger than the synchronous rotation rate, for r_p∼1, where r_p denotes the sum of radii of the colliding pair normalized by their mutual Hill radius. In this special case aggregates without internal strength may become rotationally unstable. However, the equilibrium spin rate decreases with increasing random velocity, and aggregates are always rotationally stable in the more likely case where the relative velocities are comparable to the mutual escape velocity. We also compare our results with the mean spin rates found in previous N-body simulations, and find a good agreement for optically thin rings; however the spin rates for optically thick rings are significantly larger than those predicted by our three-body calculations.     

Orbital evolution and accretion of protoplanets tidally interacting with a gas disk: I. Effects of interaction with planetesimals and other protoplanets

We have performed N-body simulations on the stage of protoplanet formation from planetesimals, taking into account so-called “type-I migration,” and damping of orbital eccentricities and inclinations, as a result of tidal interaction with a gas disk without gap formation. One of the most serious problems in formation of terrestrial planets and jovian planet cores is that the migration time scale predicted by the linear theory is shorter than the disk lifetime (10⁶-10⁷ years). In this paper, we investigate retardation of type-I migration of a protoplanet due to a torque from a planetesimal disk in which a gap is opened up by the protoplanet, and torques from other protoplanets which are formed in inner and outer regions. In the first series of runs, we carried out N-body simulations of the planetesimal disk, which ranges from 0.9 to 1.1 AU, with a protoplanet seed in order to clarify how much retardation can be induced by the planetesimal disk and how long such retardation can last. We simulated six cases with different migration speeds. We found that in all of our simulations, a clear gap is not maintained for more than 10⁵ years in the planetesimal disk. For very fast migration, a gap cannot be created in the planetesimal disk. For migration slower than some critical speed, a gap does form. However, because of the growth of the surrounding planetesimals, gravitational perturbation of the planetesimals eventually becomes so strong that the planetesimals diffuse into the vicinity of the protoplanets, resulting in destruction of the gap. After the gap is destroyed, close encounters with the planetesimals rather accelerate the protoplanet migration. In this way, the migration cannot be retarded by the torque from the planetesimal disk, regardless of the migration speed. In the second series of runs, we simulated accretion of planetesimals in wide range of semimajor axis, 0.5 to 2-5 AU, starting with equal mass planetesimals without a protoplanet seed. Since formation of comparable-mass multiple protoplanets (“oligarchic growth”) is expected, the interactions with other protoplanets have a potential to alter the migration speed. However, inner protoplanets migrate before outer ones are formed, so that the migration and the accretion process of a runaway protoplanet are not affected by the other protoplanets placed inner and outer regions of its orbit. From the results of these two series of simulations, we conclude that the existence of planetesimals and multiple protoplanets do not affect type-I migration and therefore the migration shall proceed as the linear theory has suggested.     

The formation and habitability of terrestrial planets in the presence of close-in giant planets

‘Hot jupiters,’ giant planets with orbits very close to their parent stars, are thought to form farther away and migrate inward via interactions with a massive gas disk. If a giant planet forms and migrates quickly, the planetesimal population has time to re-generate in the lifetime of the disk and terrestrial planets may form [P.J. Armitage, A reduced efficiency of terrestrial planet formation following giant planet migration, Astrophys. J. 582 (2003) L47-L50]. We present results of simulations of terrestrial planet formation in the presence of hot/warm jupiters, broadly defined as having orbital radii ?0.5 AU. We show that terrestrial planets similar to those in the Solar System can form around stars with hot/warm jupiters, and can have water contents equal to or higher than the Earth's. For small orbital radii of hot jupiters (e.g., 0.15, 0.25 AU) potentially habitable planets can form, but for semi-major axes of 0.5 AU or greater their formation is suppressed. We show that the presence of an outer giant planet such as Jupiter does not enhance the water content of the terrestrial planets, but rather decreases their formation and water delivery timescales. We speculate that asteroid belts may exist interior to the terrestrial planets in systems with close-in giant planets.     

Orbital evolution and accretion of protoplanets tidally interacting with a gas disk: II. Solid surface density evolution with type-I migration

This paper investigates the surface density evolution of a planetesimal disk due to the effect of type-I migration by carrying out N-body simulation and through analytical method, focusing on terrestrial planet formation. The coagulation and the growth of the planetesimals take place in the abundant gas disk except for a final stage. A protoplanet excites density waves in the gas disk, which causes the torque on the protoplanet. The torque imbalance makes the protoplanet suffer radial migration, which is known as type-I migration. Type-I migration time scale derived by the linear theory may be too short for the terrestrial planets to survive, which is one of the major problems in the planet formation scenario. Although the linear theory assumes a protoplanet being in a gas disk alone, Kominami et al. [Kominami, J., Tanaka, H., Ida, S., 2005. Icarus 167, 231-243] showed that the effect of the interaction with the planetesimal disk and the neighboring protoplanets on type-I migration is negligible. The migration becomes pronounced before the planet's mass reaches the isolation mass, and decreases the solid component in the disk. Runaway protoplanets form again in the planetesimal disk with decreased surface density. In this paper, we present the analytical formulas that describe the evolution of the solid surface density of the disk as a function of gas-to-dust ratio, gas depletion time scale and semimajor axis, which agree well with our results of N-body simulations. In general, significant depletion of solid material is likely to take place in inner regions of disks. This might be responsible for the fact that there is no planet inside Mercury's orbit in our Solar System. Our most important result is that the final surface density of solid components (Σd) and mass of surviving planets depend on gas surface density (Σg) and its depletion time scale (τdep) but not on initial Σd; they decrease with increase in Σg and τdep. For a fixed gas-to-dust ratio and τdep, larger initial Σd results in smaller final Σd and smaller surviving planets, because of larger Σg. To retain a specific amount of Σd, the efficient disk condition is not an initially large Σd but the initial Σd as small as the specified final one and a smaller gas-to-dust ratio. To retain Σd comparable to that of the minimum mass solar nebula (MMSN), a disk must have the same Σd and a gas-to-dust ratio that is smaller than that of MMSN by a factor of 1.3×(τdep/1 Myr) at ∼1 AU. (Equivalently, type-I migration speed is slower than that predicted by the linear theory by the same factor.) The surviving planets are Mars-sized ones in this case; in order to form Earth-sized planets, their eccentricities must be pumped up to start orbit crossing and coagulation among them. At ∼5 AU, Σd of MMSN is retained under the same condition, but to form a core massive enough to start runaway gas accretion, a gas-to-dust ratio must be smaller than that of MMSN by a factor of 3×τdep/1 Myr.     

Simulations of dense planetary rings: IV. Spinning self-gravitating particles with size distribution

Previous self-gravitating simulations of dense planetary rings are extended to include particle spins. Both identical particles as well as systems with a modest range of particle sizes are examined. For a ring of identical particles, we find that mutual impact velocity is always close to the escape velocity of the particles, even if the total rms velocity dispersion of the system is much larger, due to collective motions associated to wakes induced by near-gravitational instability or by viscous overstability. As a result, the spin velocity (i.e., the product of the particle radius and the spin frequency) maintained by mutual impacts is also of the order of the escape velocity, provided that friction is significant. For the size distribution case, smaller particles have larger impact velocities and thus larger spin velocities, particularly in optically thick rings, since small particles move rather freely between wakes. Nevertheless, the maximum ratio of spin velocities between the smallest and largest particles, as well as the ratio for translational velocities, stays below about 5 regardless of the width of the size distribution. Particle spin state is one of the important factors affecting the temperature difference between the lit and unlit face of Saturn's rings. Our results suggest that, to good accuracy, the spin frequency is inversely proportional to the particle size. Therefore, the mixing ratio of fast rotators to slow rotators on the scale of the thermal relaxation time increases with the width of the particle size distribution. This will offer means to constrain the particle size distribution with the systematic thermal infrared observations carried by the Cassini probe.     

Rotation rate and velocity dispersion of planetary ring particles with size distribution: I. Formulation and analytic calculation

We derive an equation for the evolution of rotational energy of Keplerian particles in a dilute disk due to mutual collisions. Three-dimensional Keplerian motion of particles is taken into account precisely, on the basis of Hill's approximation. The Rayleigh distribution of particles' orbital eccentricities and inclinations, and the Gaussian distribution of their rotation rates are also taken into account. Performing appropriate variable transformation, we show that the equation can be expressed with two terms. The first term, which we call collisional stirring term, represents energy exchange between rotation and random motion via collisions. The second term, which we call rotational friction term, tends to equalize the mean rotational energy of particles with different sizes. The equation can describe the evolution of rotational energy of Keplerian particles with an arbitrary size distribution. We analytically evaluate the rates of stirring and friction for the random kinetic energy and rotational energy due to inelastic collisions, for non-gravitating particles in a dilute disk. Using these results, we discuss equilibrium states in a disk of spinning, non-gravitating Keplerian particles.     

The gravitational instability in the dust layer of a protoplanetary disk:: axisymmetric solutions for nonuniform dust density distributions in the direction vertical to the midplane

The gravitational instability in the dust layer of a protoplanetary disk with nonuniform dust density distributions in the direction vertical to the midplane is investigated. The linear analysis of the gravitational instability is performed. The following assumptions are used: (1) One fluid model is adopted, that is, difference of velocities between dust and gas are neglected. (2) The gas is incompressible. (3) Models are axisymmetric with respect to the rotation axis of the disk. Numerical results show that the critical density at the midplane is higher than the one for the uniform dust density distribution by Sekiya (1983, Prog. Theor. Phys. 69, 1116-1130). For the Gaussian dust density distribution, the critical density is 1.3 times higher, although we do not consider this dust density distribution to be realistic because of the shear instability in the dust layer. For the dust density distribution with a constant Richardson number, which is considered to be realized due to the shear instability, the critical density is 2.85 times higher and is independent of the value of the Richardson number. Further, if a constant Richardson number could decrease to the order of 0.001, the gravitational instability would be realized even for the dust to gas surface density ratio with the solar abundance. Our results give a new restriction on planetesimal formation by the gravitational instability.     

High-resolution simulations of the final assembly of Earth-like planets I. Terrestrial accretion and dynamics

The final stage in the formation of terrestrial planets consists of the accumulation of ∼1000-km “planetary embryos” and a swarm of billions of 1-10 km “planetesimals.” During this process, water-rich material is accreted by the terrestrial planets via impacts of water-rich bodies from beyond roughly 2.5 AU. We present results from five high-resolution dynamical simulations. These start from 1000-2000 embryos and planetesimals, roughly 5-10 times more particles than in previous simulations. Each simulation formed 2-4 terrestrial planets with masses between 0.4 and 2.6 Earth masses. The eccentricities of most planets were ∼0.05, lower than in previous simulations, but still higher than for Venus, Earth and Mars. Each planet accreted at least the Earth's current water budget. We demonstrate several new aspects of the accretion process: (1) The feeding zones of terrestrial planets change in time, widening and moving outward. Even in the presence of Jupiter, water-rich material from beyond 2.5 AU is not accreted for several millions of years. (2) Even in the absence of secular resonances, the asteroid belt is cleared of >99% of its original mass by self-scattering of bodies into resonances with Jupiter. (3) If planetary embryos form relatively slowly, then the formation of embryos in the asteroid belt may have been stunted by the presence of Jupiter. (4) Self-interacting planetesimals feel dynamical friction from other small bodies, which has important effects on the eccentricity evolution and outcome of a simulation.     

Exoplanetary systems: The role of an equilibrium at high mutual inclination in shaping the global behavior of the 3-D secular planetary three-body problem

On the basis of a high-order (order 12) expansion of the perturbative potential in powers of the eccentricities and the inclinations, we analyze the secular interactions of two non-coplanar planets which are not in mean-motion resonance. The model is based on the planetary three-body problem which can be reduced to two degrees of freedom by the well-known elimination of the nodes [Jacobi, C.G.J., 1842. Astron. Nachr. XX, 81-102]. We introduce non-singular canonical variables which bring forward the symmetries of the problem. The main dynamical features depend on the location and stability of the equilibria which are easily found with our analytical model. We find that there exists an equilibrium when both eccentricities are zero. When the mutual inclination is small, this equilibrium is stable, but for larger mutual inclination it becomes unstable, generating a large chaotic zone and, by bifurcation, two regular regions, the so-called Kozai resonances. This analytical study which depends on only two parameters (the ratio of the semi-major axes and the mass ratio of the planets) makes possible a large survey of the problem and enables us to identify and quantify its main dynamical features, periodic orbits, regular and chaotic zones, etc. The results of our analytical model are illustrated and confirmed by numerical integrations.     

The E-ring in the vicinity of Enceladus: II. Probing the moon's interior—The composition of E-ring particles

The population of Saturn's outermost tenuous E-ring is dominated by tiny water ice particles. Active volcanism on the moon Enceladus, embedded in the E-ring, has since late 2005 been known to be a major source of particles replenishing the ring. Therefore particles in the vicinity of Enceladus may provide crucial information about the dynamical and chemical processes occurring below the moon's icy surface. Here we present a statistical evaluation of more than 2000 impact ionisation mass spectra of Saturn's E-ring particles, with sizes predominantly below 1 μm, detected by the Cosmic Dust Analyser onboard the Cassini spacecraft. We focus on the identification of non-water features in spectra otherwise dominated by water ice signatures. Here we specify the categorisation of two different spectrum types, which probably represent two particle populations. Type I spectra imply pure water ice particles, whereas in Type II spectra organic compounds and/or silicate minerals are identified as impurities within the icy particles. This finding supports the hypothesis of a dynamic interaction of Enceladus' rocky core with liquid water.     

Photometric modeling of Saturn's rings: I. Monte Carlo method and the effect of nonzero volume filling factor

The scattering properties of particulate rings with volume filling factors in the interval D=0.001-0.3 are studied, with photometric Monte Carlo ray tracing simulations combining the advantages of direct (photons followed from the source) and indirect methods (brightness as seen from the observing direction). Besides vertically homogeneous models, ranging from monolayers to classical many-particle thick rings, particle distributions obtained from dynamical simulations are studied, possessing a nonuniform vertical profile and a power law distribution of particle sizes. Self-gravity is not included to assure homogeneity in planar directions. Our main goal is to check whether the moderately flattened ring models predicted by dynamical simulations (with central plane D>0.1) are consistent with the basic photometric properties of Saturn's rings seen in ground-based observations, including the brightening near zero phase angle (opposition effect), and the brightening of the B-ring with increasing elevation angle (tilt effect). Our photometric simulations indicate that dense rings are typically brighter in reflected light than those with D→0, due to enhanced single scattering. For a vertically illuminated layer of identical particles this enhancement amounts at intermediate viewing elevations to roughly 1+2D. Increased single scattering is also obtained for low elevation illumination, further augmented at low phase angles α by the opposition brightening when D increases: the simulated opposition effect agrees very well with the Lumme and Bowell (1981, Astron. J. 86, 1694-1704) theoretical formula. For large α the total intensity may also decrease, due to reduced amount of multiple scattering. For the low (α=13°) and high (α=155°) phase angle geometries analyzed in Dones et al. (1993, Icarus 105, 184-215) the brightness change for D=0.1 amounts to 20% and −17%, respectively. In the case of an extended size distribution, dynamical simulations indicate that the smallest particles typically occupy a layer several times thicker than the largest particles. Even if the large particles form a dynamically dense system, a narrow opposition peak can arise due to mutual shadowing among the small particles: for example, a size distribution extending about two decades can account for the observed about 1° wide opposition peak, solely in terms of mutual shadowing. The reduced width of the opposition peak for extended size distribution is in accordance with Hapke's (1986, Icarus 67, 264-280) treatment for semi-infinite layers. Due to vertical profile and particle size distribution, the photometric behavior is sensitive to the viewing elevation: this can account for the tilt-effect of the B-ring, as dense and thus bright central parts of the ring become better visible for larger elevation, whereas in the case of smaller elevation, mainly low volume density upper layers are visible. Since multiple scattering is not involved, the explanation works also for albedo well below unity. Inclusion of nonzero volume density helps also to model some of the Voyager observations. For example, the discrepancy between predicted and observed brightness at large phase angles for much of the A-ring (Dones et al., 1993, Icarus 105, 184-215) is removed when the enhanced low α single scattering and reduced large α multiple scattering is allowed for. Also, a model with vertical thickness increasing with saturnocentric distance offers at least a qualitative explanation for the observed contrast reversal between the inner and outer A-ring in low and high phase Voyager images. Differences in local size distribution and thus on the effective D may also account for the contrast reversal in resonance sites.     

From planetesimals to terrestrial planets: N-body simulations including the effects of nebular gas and giant planets

We present results from a suite of N-body simulations that follow the formation and accretion history of the terrestrial planets using a new parallel treecode that we have developed. We initially place 2000 equal size planetesimals between 0.5 and 4.0 AU and the collisional growth is followed until the completion of planetary accretion (>100 Myr). A total of 64 simulations were carried out to explore sensitivity to the key parameters and initial conditions. All the important effect of gas in laminar disks are taken into account: the aerodynamic gas drag, the disk-planet interaction including Type I migration, and the global disk potential which causes inward migration of secular resonances as the gas dissipates. We vary the initial total mass and spatial distribution of the planetesimals, the time scale of dissipation of nebular gas (which dissipates uniformly in space and exponentially in time), and orbits of Jupiter and Saturn. We end up with 1-5 planets in the terrestrial region. In order to maintain sufficient mass in this region in the presence of Type I migration, the time scale of gas dissipation needs to be 1-2 Myr. The final configurations and collisional histories strongly depend on the orbital eccentricity of Jupiter. If today’s eccentricity of Jupiter is used, then most of bodies in the asteroidal region are swept up within the terrestrial region owing to the inward migration of the secular resonance, and giant impacts between protoplanets occur most commonly around 10 Myr. If the orbital eccentricity of Jupiter is close to zero, as suggested in the Nice model, the effect of the secular resonance is negligible and a large amount of mass stays for a long period of time in the asteroidal region. With a circular orbit for Jupiter, giant impacts usually occur around 100 Myr, consistent with the accretion time scale indicated from isotope records. However, we inevitably have an Earth size planet at around 2 AU in this case. It is very difficult to obtain spatially concentrated terrestrial planets together with very late giant impacts, as long as we include all the above effects of gas and assume initial disks similar to the minimum mass solar nebular.     

On the Light-Absorbing Surface Layer of Cometary Nuclei: I. Radiative Transfer

A small but increasing volume of observations of cometary nuclei has accumulated during the past two decades. This development is accelerating with upcoming space missions such as Stardust, Contour, and Rosetta. In response to the growing need for a theoretical understanding of optical properties of cometary nuclei, we have calculated synthetic reflectance spectra in the wavelength region 0.2-2.0 μm, photometric colors in the Johnson-Kron-Cousins UBVRI system, and visual geometric albedos for a large number of porous ice-dust mixtures with differing composition, regolith grain sizes, and grain morphologies, such as core-mantle grains, dense clusters of such grains, and large irregular particles with internal scatterers. The calculations are based on Mie theory, the discrete dipole approximation, Hapke theory, and a numerical solution to the equation of radiative transfer in particulate media. In addition, wavelength-integrated directional-hemispherical albedos and flux attenuation profiles in the regolith as functions of depth have been calculated in order to improve the energy budget and treatment of energy boundary conditions in thermal models of cometary nuclei.Our results are compared with spectra and colors of observed cometary nuclei. Our main conclusions are that only regolith consisting of relatively large core-mantle grains, or clusters of smaller core-mantle grains, is capable of reproducing the red colors seen in comets; that ice-dust mixtures actually can be darker than ice-free regolith in certain circumstances; and that solar radiation sometimes penetrates to a depth that is comparable to the region in which diurnal temperature variations occur.     

Corona-like atmospheric escape from KBOs: I. Gas dynamics

We show that for low temperatures (T∼30 K) and small, but non-negligible, gravitational fields the hydrodynamic escape of gas can be treated by Parker's theory of coronal expansion [Parker, E.N., 1963. Interplanetary Dynamical Processes. Interscience Publishers, New York]. We apply this theory to gas escape from Kuiper belt objects. We derive limits on the density and radius of the bodies for which this theory is applicable, and show how the flow depends on the mean molecular weight and internal degrees of freedom of the gas molecules. We use these results to explain the CH₄ dichotomy seen on KBOs [Schaller, E.L., Brown, M.E., 2007. Astrophys. J., 659, L61-L64].