Embedded star clusters and the formation of the Oort cloud: II. The effect of the primordial solar nebula

This paper deals with Oort cloud formation while the Sun was in an embedded cluster and surrounded by its primordial nebula. This work is a continuation of Brasser et al. [Brasser, R., Duncan, M., Levison, H., 2006. Icarus 184, 59-82], building on the model presented therein, and adding the aerodynamic drag and gravitational potential of the primordial solar nebula. Results are presented of numerical simulations of comets subject to the gravitational influence of the Sun, Jupiter, Saturn, star cluster and primordial solar nebula; some of the simulations included the gravitational influence of Uranus and Neptune as well. The primordial solar nebula was approximated by the minimum-mass Hayashi model [Hayashi, C., Nakozawa, K., Nakagawa, Y., 1985. In: Black, D.C., Matthews, M.S. (Eds.). Protostars and Planets II. Univ. of Arizona Press, Tucson, AZ] whose inner and outer radii have been truncated at various distances from the Sun. A comet size of 1.7 km was used for most of our simulations. In all of our simulations, the density of the primordial solar nebula decayed exponentially with an e-folding time of 2 Myr. It turns out that when the primordial solar nebula extends much beyond Saturn or Neptune, virtually no material will end up in the Oort cloud (OC) during this phase. Instead, the majority of the material will be on circular orbits inside of Jupiter if the inner edge of the disk is well inside Jupiter's orbit. If the disk's inner edge is beyond Jupiter's orbit, most comets end up on orbits in exterior mean-motion resonances with Saturn when Uranus and Neptune are not present. In those cases where the outer edge of the disk is close to Saturn or Neptune, the fraction of material that ends up in the subsequently formed OC is much less than that found in Brasser et al. [Brasser, R., Duncan, M., Levison, H., 2006. Icarus 184, 59-82] for the same cluster densities. This implies that for comets of roughly 2 km in size, the presence of the primordial solar nebula hinders OC formation. A byproduct of some of our simulations are endresults with a substantial fraction of the comets in the Uranus-Neptune scattered disk. A subsequent followup of this material is planned for the near future. In order to determine the effect of the size of the comets on OC formation efficiency, a set of runs with the same initial conditions but different cometary radii have been performed as well, from which it is determined that the threshold comet size to begin producing significant Oort clouds is roughly 20 km. This implies that the presence of the primordial solar nebula acts as a size-sorting mechanism, with large bodies unaffected by the gas drag and ending up in the OC while small bodies remain trapped in the planetary region, in the models studied.     

Reassessing the formation of the inner Oort cloud in an embedded star cluster

We re-examine the formation of the inner Oort comet cloud while the Sun was in its birth cluster with the aid of numerical simulations. This work is a continuation of an earlier study (Brasser, R., Duncan, M.J., Levison, H.F. [2006]. Icarus 184, 59–82) with several substantial modifications. First, the system consisting of stars, planets and comets is treated self-consistently in our N-body simulations, rather than approximating the stellar encounters with the outer Solar System as hyperbolic fly-bys. Second, we have included the expulsion of the cluster gas, a feature that was absent previously. Third, we have used several models for the initial conditions and density profile of the cluster – either a Hernquist or Plummer potential – and chose other parameters based on the latest observations of embedded clusters from the literature. These other parameters result in the stars being on radial orbits and the cluster collapses. Similar to previous studies, in our simulations the inner Oort cloud is formed from comets being scattered by Jupiter and Saturn and having their pericentres decoupled from the planets by perturbations from the cluster gas and other stars. We find that all inner Oort clouds formed in these clusters have an inner edge ranging from 100 AU to a few hundred AU, and an outer edge at over 100,000 AU, with little variation in these values for all clusters. All inner Oort clouds formed are consistent with the existence of (90377) Sedna, an inner Oort cloud dwarf planetoid, at the inner edge of the cloud: Sedna tends to be at the innermost 2% for Plummer models, while it is 5% for Hernquist models. We emphasise that the existence of Sedna is a generic outcome. We define a ‘concentration radius’ for the inner Oort cloud and find that its value increases with increasing number of stars in the cluster, ranging from 600 AU to 1500 AU for Hernquist clusters and from 1500 AU to 4000 AU for Plummer clusters. The increasing trend implies that small star clusters form more compact inner Oort clouds than large clusters. We are unable to constrain the number of stars that resided in the cluster since most clusters yield inner Oort clouds that could be compatible with the current structure of the outer Solar System. The typical formation efficiency of the inner Oort cloud is 1.5%, significantly lower than previous estimates. We attribute this to the more violent dynamics that the Sun experiences as it rushes through the centre of the cluster during the latter’s initial phase of violent relaxation.     

Long-Period Comets and the Oort Cloud

Long-period (LP) comets, Halley-type (HT) comets, and even some comets of the Jupiter family, probably come from the Oort cloud, a huge reservoir of icy bodies that surrounds the solar system. Therefore, these comets become important probes to learn about the distant Oort cloud population. We review the fundamental dynamical properties of LP comets, and what is our current understanding of the dynamical mechanisms that bring these bodies from the distant Oort cloud region to the inner planetary region. Most new comets have original reciprocal semimajor axes in the range2 × 10^-5 < 1/a_orig < 5 × 10^-5AU^-1. Yet, this cannot be taken to represent the actual space distribution of Oort cloud comets, but only the region in the energy space in which external perturbers have the greatest efficiency in bringing comets to the inner planetary region. The flux of Oort cloud comets in the outer planetary region is found to be at least several tens times greater than the flux in the inner planetary region. The sharp decrease closer to the Sun is due to the powerful gravitational fields of Jupiter and Saturn that prevent most Oort cloud comets from reaching the Earth’s neighborhood (they act as a dynamical barrier). A small fraction of ～10^-2 Oort cloud comets become Halley type (orbital periods P < 200 yr), and some of them can reach short-period orbits with P < 20 yr. We analyze whether we can distinguish the latter, very ‘old” LP comets, from comets of the Jupier family coming from the Edgeworth-Kuiper belt.     

Long-term effects of the Galactic tide on cometary dynamics

We introduce a model for integrating the effects of Galactic tides on Oort cloud comets, which involves two procedures, according to the values of the osculating semi-major axis a and eccentricity e. Ten simulations of the dynamics of 10⁶ comets over 5 Gyr are performed using this model. We thus investigate the long-term effects of the Galactic tide with and without a radial component, the effects of the local density of the Galactic disk, and those of the Oort constants. Most of the results may be understood in terms of the integrability or non-integrability of the system. For an integrable system, which occurs for moderate semi-major axes with or without radial component, the dynamics is explained by periodic variation of the cometary perihelion, inducing the depletion of the outer region of the Oort cloud, a constant flux from the inner region after 500 Myr, and the quick formation of a reservoir of comets with argument of perihelion near 26.6°. When the system is non-integrable, the efficiency of the tide in reducing the cometary perihelion distance is enhanced both by replenishing the Oort cloud domain from which comets are sent toward the planetary system, and by reducing the minimal value that the perihelion distance may reach. No effects of varying the Oort constants were observed, showing that the flat rotation curve is a satisfactory approximation in Oort cloud dynamics.     

The frequency and intensity of comet showers from the Oort cloud

We simulate the Oort comet cloud to study the rate and properties of new comets and the intensity and frequency of comet showers. An ensemble of ∼10⁶ comets is perturbed at random times by a population of main sequence stars and white dwarfs that is described by the Bahcall-Soneira Galaxy model. A cloning procedure allows us to model a large ensemble of comets efficiently, without wasting computer time following a large number of low eccentricity orbits. For comets at semimajor axis a = 20,000 AU, about every 100 myr a star with mass in the range 1M_⊙−2M_⊙ passes within ∼10,000 AU of the Sun and triggers a shower that enhances the flux of new comets by more than a factor of 10. The time-integrated flux is dominated by the showers for comets with semimajor axes less than ∼30,000 AU. For semimajor axes greater than ∼30,000 AU the comet loss rate is roughly constant and strong showers do not occur. In some of our simulations, comets are also perturbed by the Galactic tidal field. The inclusion of tidal effects increases the loss rate of comets with semimajor axes between 10,000 and 20,000 AU by about a factor of 4. Thus the Galactic tide, rather than individual stellar perturbations, is the dominant mechanism which drives the evolution of the Oort cloud.     

Nuclear spin temperature of ammonia in Comet 9P/Tempel 1 before and after the Deep Impact event

The Deep Impact mission succeeded in excavating inner materials from the nucleus of Comet 9P/Tempel 1 on 2005 July 04 (at 05:52 UT). Comet 9P/Tempel 1 is one of Jupiter family short period comets, which might originate in the Kuiper belt region in the solar nebula. In order to characterize the comet and to support the mission from the ground-based observatory, optical high-dispersion spectroscopic observations were carried out with the echelle spectrograph (UVES) mounted on the 8-m telescope VLT (UT2) before and after the Deep Impact event. Ortho-to-para abundance ratios (OPRs) of cometary ammonia were determined from the NH₂ emission spectra. The OPRs of ammonia on July 3.996 UT and 4.997 UT were derived to be 1.28±0.07 (nuclear spin temperature: Tspin=24±2 K) and 1.26±0.08 (Tspin=25±2 K), respectively. There is no significant change between before and after the impact. Actually, most materials ejected from the impact site could have moved away from the nucleus on July 4.997 UT, about 17 h after the impact. However, a small fraction of the ejected materials might remain in the slit of UVES instrument at that time because an excess of about 20% in the NH₂ emission flux is observed above the normal activity level was found [Manfroid, J., Hutsemékers, D., Jehin, E., Cochran, A.L., Arpigny, C., Jackson, W.M., Meech, K.J., Schulz, R., Zucconi, J.-M., 2007. Icarus. This issue]. If the excess of NH₂ on July 04.997 UT was produced from icy materials excavated by the Deep Impact, then an upper-limit of the ammonia OPR would be 1.75 (Tspin>17 K) for those materials. On the other hand, the OPR of ammonia produced from the quiescent sources was similar to that of the Oort cloud comets observed so far. This fact may imply that physical conditions where cometary ices formed were similar between Comet 9P/Tempel 1 and the Oort cloud comets.     

A Model for the Common Origin of Jupiter Family and Halley Type Comets

A numerical simulation of the Oort cloud is used to explain the observed orbital distributions and numbers of Jupiter-family (JF) and Halley-type (HT) short-period (SP) comets. Comets are given initial orbits with perihelion distances between 5 and 36 au, and evolve under planetary, stellar and Galactic perturbations for 4.5 Gyr. This process leads to the formation of an Oort cloud (which we define as the region of semimajor axes a > 1,000 au), and to a flux of cometary bodies from the Oort cloud returning to the planetary region at the present epoch. The results are consistent with the dynamical characteristics of SP comets and other observed cometary populations: the near-parabolic flux, Centaurs, and high-eccentricity trans-Neptunian objects. To achieve this consistency with observations, the model requires that the number of comets versus initial perihelion distance is concentrated towards the outer planetary region. Moreover, the mean physical lifetime of observable comets in the inner planetary region (q < 2.5 au) at the present epoch should be an increasing function of the comets’ initial perihelion distances. Virtually all observed HT comets and nearly half of observed JF comets come from the Oort cloud, and initially (4.5 Gyr ago) from orbits concentrated near the outer planetary region. Comets that have been in the Oort cloud also return to the Centaur (5 < q < 28 au, a < 1,000 au) and near-Neptune high-eccentricity regions. Such objects with perihelia near Neptune are hard to discover, but Centaurs with characteristics predicted by the model (e.g. large semimajor axes, above 60 au, or high inclinations, above 40°) are increasingly being found by observers. The model provides a unified picture for the origin of JF and HT comets. It predicts that the mean physical lifetime of all comets in the region q < 1.5 au is less than ～200 revolutions.     

Persistent evidence of a jovian mass solar companion in the Oort cloud

We present updated dynamical and statistical analyses of outer Oort cloud cometary evidence suggesting that the Sun has a wide-binary jovian mass companion. The results support a conjecture that there exists a companion of mass ≈ orbiting in the innermost region of the outer Oort cloud. Our most restrictive prediction is that the orientation angles of the orbit plane in galactic coordinates are centered on Ω, the galactic longitude of the ascending node = 319° and i, the galactic inclination = 103° (or the opposite direction) with an uncertainty in the orbit normal direction subtending <2% of the sky. Such a companion could also have produced the detached Kuiper Belt object Sedna. If the object exists, the absence of similar evidence in the inner Oort cloud implies that common beliefs about the origin of observed inner Oort cloud comets must be reconsidered. Evidence of the putative companion would have been recorded by the Wide-field Infrared Survey Explorer (WISE) which has completed its primary mission and is continuing on secondary objectives.     

The colors of cometary nuclei—Comparison with other primitive bodies of the Solar System and implications for their origin

We present new color results of cometary nuclei obtained with the Hubble Space Telescope (HST) whose superior resolution enables us to accurately isolate the nucleus signals from the surrounding comae. By combining with scrutinized available data obtained with ground-based telescopes, we accumulated a sample of 51 cometary nuclei, 44 ecliptic comets (ECs) and 7 nearly-isotropic comets (NICs) using the nomenclature of Levison [Levison, H.F., 1996. In: Rettig, T.W., Hahn, J.M. (Eds.), Completing the Inventory of the Solar System. In: ASP Conf. Ser., vol. 107, pp. 173-192]. We analyze color distributions and color-color correlations as well as correlations with other physical parameters. We present our compilation of colors of 232 outer Solar System objects—separately considering the different dynamical populations, classical KBOs in low and high-inclination orbits (respectively CKBO-LI and CKBO-HI), resonant KBOs (practically Plutinos), scattered-disk objects (SDOs) and Centaurs—of 12 candidate dead comets, and of 85 Trojans. We perform a systematic analysis of all color distributions, and conclude by synthesizing the implications of the dynamical evolution and of the colors for the origin of the minor bodies of the Solar System. We find that the color distributions are remarkably consistent with the scenarios of the formation of TNOs by Gomes [Gomes, R.S., 2003. Icarus 161, 404-418] generalized by the “Nice” model [Levison, H.F., Morbidelli, A., VanLaerhoven, Ch., Gomes, R., Tsiganis, L., 2008. Icarus 196, 258-273], and of the Trojans by Morbidelli et al. [Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R., 2005. Nature 435, 462-465]. The color distributions of the Centaurs are globally similar to those of the CKBO-HI, the Plutinos and the SDOs. However the potential bimodality of their distributions allows to possibly distinguish two groups based on their (B−R) index: Centaur I with (B−R)>1.7 and Centaurs II with (B−R)<1.4. Centaurs I could be composed of TNOs (prominently CKBO-LI) and ultra red objects from a yet unstudied family. Centaurs II could consist in a population of evolved objects which have already visited the inner Solar System, and which has been scattered back beyond Jupiter. The diversity of colors of the ECs, in particular the existence of very red objects, is consistent with an origin in the Kuiper belt. Candidate dead comets represent an ultimate state of evolution as they appear more evolved than the Trojans and Centaurs II.     

Long-period comet flux in the planetary region: Dynamical evolution from the Oort cloud

This study analyzes the evolution of 2 × 10⁵ orbits with initial parameters corresponding to the orbits of comets of the Oort cloud under the action of planetary, galactic, and stellar perturbations over 2 × 10⁹ years. The dynamical evolution of comets of the outer (orbital semimajor axes a > 10⁴ AU) and inner (5 × 10³ < a (AU) < 10⁴) parts of the comet cloud is analyzed separately. The estimates of the flux of “new” and long-period comets for all perihelion distances q in the planetary region are reported. The flux of comets with a > 10⁴ AU in the interval 15 AU < q < 31 AU is several times higher than the flux of comets in the region q < 15 AU. We point out the increased concentration of the perihelia of orbits of comets from the outer cloud, which have passed several times through the planetary system, in the Saturn-Uranus region. The maxima in the distribution of the perihelia of the orbits of comets of the inner Oort cloud are located in the Uranus-Neptune region. “New” comets moving in orbits with a < 2 × 10⁴ AU and arriving at the outside of the planetary system (q > 25 AU) subsequently have a greater number of returns to the region q < 35 AU. The perihelia of the orbits of these comets gradually drift toward the interior of the Solar System and accumulate beyond the orbit of Saturn. The distribution of the perihelia of long-period comets beyond the orbit of Saturn exhibits a peak. We discuss the problem of replenishing the outer Oort cloud by comets from the inner part and their subsequent dynamical evolution. The annual rate of passages of comets of the inner cloud, which replenish the outer cloud, in the region q < 1 AU in orbits with a > 10⁴ AU (～ 5.0 × 10^?14 yr^?1) is one order of magnitude lower than the rate of passage of comets from the outer Oort cloud (～ 9.1 × 10^?13 yr^?1).     

The formation of the Oort cloud in open cluster environments

We study the influence of an open cluster environment on the formation and current structure of the Oort cloud. To do this, we have run 19 different simulations of the formation of the Oort cloud for 4.5 Gyrs. In each simulation, the Solar System spends its first 100 Myrs in a different open cluster environment before transitioning to its current field environment. We find that, compared to forming in the field environment, the inner Oort cloud is preferentially loaded with comets while the Sun resides in the open cluster and that most of this material remains locked in the interior of the cloud for the next 4.4 Gyrs. In addition, the outer Oort cloud trapping efficiencies we observe in our simulations are lower than previous formation models by about a factor of 2, possibly implying an even more massive early planetesimal disk. Furthermore, some of our simulations reproduce the orbits of observed extended scattered disk objects, which may serve as an observational constraint on the Sun's early environment. Depending on the particular open cluster environment, the properties of the inner Oort cloud and extended scattered disk can vary widely. On the other hand, the outer portions of the Oort cloud in each of our simulations are all similar.     

Injection of Oort Cloud comets: the fundamental role of stellar perturbations

We present Monte Carlo simulations of the dynamical evolution of the Oort cloud over the age of the Solar System, using an initial sample of one million test comets without any cloning. Our model includes perturbations due to the Galactic tide (radial and vertical) and passing stars. We present the first detailed analysis of the injection mechanism into observable orbits by comparing the complete model with separate models for tidal and stellar perturbations alone. We find that a fundamental role for injecting comets from the region outside the loss cone (perihelion distance q > 15 AU) into observable orbits (q < 5 AU) is played by stellar perturbations. These act in synergy with the tide such that the total injection rate is significantly larger than the sum of the two separate rates. This synergy is as important during comet showers as during quiescent periods and concerns comets with both small and large semi-major axes. We propose different dynamical mechanisms to explain the synergies in the inner and outer parts of the Oort Cloud. We find that the filling of the observable part of the loss cone under normal conditions in the present-day Solar System rises from <1% for a < 20 000 AU to about 100% for a ? 100 000 AU.     

The size-frequency distribution of dormant Jupiter family comets

We estimate the total number and the slope of the size-frequency distribution (SFD) of dormant Jupiter family comets (JFCs) by fitting a one-parameter model to the known population. We first select 61 near-Earth objects (NEOs) that are likely to be dormant JFCs because their orbits are dynamically coupled to Jupiter [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J., Levison, H.F., Michel, P., Metcalfe, T.S., 2002a. Icarus 156, 399-433]. Then, from the numerical simulations of Levison and Duncan [1997. Icarus 127, 13-32], we construct an orbit distribution model for JFCs in the NEO orbital element space. We assume an orbit-independent SFD for all JFCs, the slope of which is our unique free parameter. Finally, we compute observational biases for dormant JFCs using a calibrated NEO survey simulator [Jedicke, R., Morbidelli, A., Spahr, T., Petit, J., Bottke, W.F., 2003. Icarus 161, 17-33]. By fitting the biased model to the data, we estimate that there are ∼75 dormant JFCs with H<18 in the NEO region and that the slope of their cumulative SFD is −1.5±0.3. Our slope for the SFD of dormant JFCs is very close to that of active JFCs as determined by Weissman and Lowry [2003. Lunar Planet. Sci. 34. Abstract 2003]. Thus, we argue that when JFCs fade they are likely to become dormant rather than to disrupt and that the fate of faded comets is size-independent. Our results imply that the size distribution of the JFC progenitors—the scattered disk trans-neptunian population—either (i) has a similar and shallow SFD or (i^′) is slightly steeper and physical processes acting on the comets in a size-dependent manner creates the shallower active comet SFD. Our measured slope, typical of collisionally evolved populations with a size-dependent impact strength [Benz, W., Asphaug, E., 1999. Icarus 142, 5-20], suggests that scattered disk bodies reached collisional equilibrium inside the protoplanetary disk prior to their removal from the planetary region.     

The orbit and size distribution of small Solar System objects orbiting the Sun interior to the Earth's orbit

We present the first observational measurement of the orbit and size distribution of small Solar System objects whose orbits are wholly interior to the Earth's (Inner Earth Objects, IEOs, with aphelion <0.983 AU). We show that we are able to model the detections of near-Earth objects (NEO) by the Catalina Sky Survey (CSS) using a detailed parameterization of the CSS survey cadence and detection efficiencies as implemented within the Jedicke et al. [Jedicke, R., Morbidelli, A., Spahr, T., Petit, J.M., Bottke, W.F., 2003. Icarus 161, 17-33] survey simulator and utilizing the Bottke et al. [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H.F., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433] model of the NEO population's size and orbit distribution. We then show that the CSS detections of 4 IEOs are consistent with the Bottke et al. [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H.F., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433] IEO model. Observational selection effects for the IEOs discovered by the CSS were then determined using the survey simulator in order to calculate the corrected number and H distribution of the IEOs. The actual number of IEOs with H<18 (21) is 36±26 (530±240) and the slope of the H magnitude distribution (∝¹⁰αH) for the IEOs is . The slope is consistent with previous measurements for the NEO population of αNEO=0.35±0.02 [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H.F., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433] and αNEO=0.39±0.013 [Stuart, J.S., Binzel, R.P., 2004. Icarus 170, 295-311]. Based on the agreement between the predicted and observed IEO orbit and absolute magnitude distributions there is no indication of any non-gravitational effects (e.g. Yarkovsky, tidal disruption) affecting the known IEO population.     

Perihelion evolution of observed new comets implies the dominance of the galactic tide in making Oort cloud comets discernable

For an Oort cloud comet to be seen as a new comet, its perihelion must be moved from a point exterior to the loss cylinder boundary to a point interior to observable limits in a single orbit. The galactic tide can do this continuously, in a non-impulsive manner. Near-parabolic comets, with specific angular momentum , will most easily be made observable. Therefore, to reduce the perihelion distance H must decrease. Since weakly perturbed comets are, in general, more numerous than strongly perturbed comets, we can anticipate that new comets made observable by a weak tidal torque will more likely be first observed when their slowly changing perihelion distances are approaching their minimum osculating values under the action of the tide, rather than receding from their minimum values. That is, defining ΔH_tide as the vector change due to the galactic tidal torque during the prior orbit, and H_obs as the observed vector, the sign S≡Sign(H_obs·ΔH_tide) will more likely be −1 than +1 if a weak galactic tidal perturbation indeed dominates in making comets observable. Using comet data of the highest quality class (1A) for new comets (a>10,000 AU), we find that 49 comets have S=−1 and 22 have S=+1. The binomial probability that as many or more would exhibit this characteristic if in fact S=?1 were equally likely is only 0.0009. This characteristic also persists in other long-period comet populations, lending support to the notion that they are dominated by comets recently arrived from the outer Oort cloud. The preponderance of S=−1 also correlates with weakly perturbed (i.e., smaller semimajor axis) new comets in a statistically significant manner. This is strong evidence that the data are of sufficiently high quality and sufficiently free of observational selection effects to detect this unique imprint of the tide.     

Peripheral Structures of Planetary Systems

We analyze the conditions for the formation and time evolution of peripheral comet structures of solar-type planetary systems. In the Solar system, these include the Kuiper belt, the Oort cloud, the comet spear, and the Galactic comet ring that marks the Galactic orbit of the Sun. We consider the role of the viscosity of a protoplanetary gas–dust disk, major planets, field stars, globular clusters, giant molecular clouds, and the Galactic gravitational field in the formation of these peripheral structures marked by comets and asteroids. We give a list of the closest past and future passages of neighboring stars through the solar Oort cloud that perturb the motion of its comets and, thus, contribute to the enhancement of its cometary activity, on the one hand, and to the replenishment of the solar comet spear with new members, on the other hand.     

The internal structure of Jupiter family cometary nuclei from Deep Impact observations: The “talps” or “layered pile” model

We consider the hypothesis that the layering observed on the surface of Comet 9P/Tempel 1 from the Deep Impact spacecraft and identified on other comet nuclei imaged by spacecraft (i.e., 19P/Borrelly and 81P/Wild 2) is ubiquitous on Jupiter family cometary nuclei and is an essential element of their internal structure. The observational characteristics of the layers on 9P/Tempel 1 are detailed and considered in the context of current theories of the accumulation and dynamical evolution of cometary nuclei. The works of Donn [Donn, B.D., 1990. Astron. Astrophys. 235, 441-446], Sirono and Greenberg [Sirono, S.-I., Greenberg, J.M., 2000. Icarus 145, 230-238] and the experiments of Wurm et al. [Wurm, G., Paraskov, G., Krauss, O., 2005. Icarus 178, 253-263] on the collision physics of porous aggregate bodies are used as basis for a conceptual model of the formation of layers. Our hypothesis is found to have implications for the place of origin of the JFCs and their subsequent dynamical history. Models of fragmentation and rubble pile building in the Kuiper belt in a period of collisional activity (e.g., [Kenyon, S.J., Luu, J.X., 1998. Astron. J. 115, 2136-2160; 1999a. Astron. J. 118, 1101-1119; 1999b. Astrophys. J. 526, 465-470; Farinella, P., Davis, D.R., Stern, S.A., 2000. In: Mannings, V., Boss, A.P., Russell, S.S. (Eds.), Protostars and Planets IV. Univ. of Arizona Press, Tucson, pp. 1255-1282; Durda, D.D., Stern, S.J., 2000. Icarus 145, 220-229]) following the formation of Neptune appear to be in conflict with the observed properties of the layers and irreconcilable with the hypothesis. Long-term residence in the scattered disk [Duncan, M.J., Levison, H.F., 1997. Science 276, 1670-1672; Duncan, M., Levison, H., Dones, L., 2004. In: Festou, M., Keller, H.U., Weaver, H.A. (Eds.), Comets II. Univ. of Arizona Press, Tucson, pp. 193-204] and/or a change in fragmentation outcome modeling may explain the long-term persistence of primordial layers. In any event, the existence of layers places constraints on the environment seen by the population of objects from which the Jupiter family comets originated. If correct, our hypothesis implies that the nuclei of Jupiter family comets are primordial remnants of the early agglomeration phase and that the physical structure of their interiors, except for the possible effects of compositional phase changes, is largely as it was when they were formed. We propose a new model for the interiors of Jupiter family cometary nuclei, called the talps or “layered pile” model, in which the interior consists of a core overlain by a pile of randomly stacked layers. We discuss how several cometary characteristics—layers, surface texture, indications of flow, compositional inhomogeneity, low bulk density low strength, propensity to split, etc., might be explained in terms of this model. Finally, we make some observational predictions and suggest goals for future space observations of these objects.     

2-Gyr Simulation of the Oort-cloud Formation II. A Close View of the Inner Oort cloud after the First Two Giga-years

We simulate the formation of the Oort cloud (OC) till the age of 2 Gyr starting from an initial disc of planetesimals made by 10 038 test particles. The results on the outer part of the distant comet reservoir are reported by Neslu?an et al. (this issue). Here we deal with the evolution of the population and structure at 2 Gyr of the complementary inner part of the Oort cloud. The dynamical evolution of the massless test particles was followed via the numerical integration of their orbits. We considered the perturbations produced by four giant planets assuming they have their current orbits and masses, as well as the perturbations caused by the Galactic tide and passing stars. The efficiency of the formation of inner OC is found to be very low: only about 1.1% of all considered particles ended in this part of the OC. At 2 Gyr, the dynamics of the inner cloud is mainly governed by the dominant z-term of the Galactic tide. The number density of the bodies is proportional to the heliocentric distance, r, as r ^?3.53. The directional distribution of orbits is still strongly inhomogeneous. There are large empty regions in the space angles around the Galactic Equator points with the galactic longitude 90 and 270° (non-rotating frame), or there are only few bodies having the ecliptical latitude higher than +60° or lower than 60°. A strong concentration of objects at the Ecliptic is apparent up to ≈1,000 AU, with a possible—but still not proved—extension to ≈1,500 AU. Beyond r ≈ 6,000 AU, bodies directly above and below the Sun, with respect to the Ecliptic, are absent.     

The formation of Uranus and Neptune in solid-rich feeding zones: Connecting chemistry and dynamics

The core accretion theory of planet formation has at least two fundamental problems explaining the origins of Uranus and Neptune: (1) dynamical times in the trans-saturnian solar nebula are so long that core growth can take >15 Myr and (2) the onset of runaway gas accretion that begins when cores reach ∼10M_⊕ necessitates a sudden gas accretion cutoff just as Uranus and Neptune’s cores reach critical mass. Both problems may be resolved by allowing the ice giants to migrate outward after their formation in solid-rich feeding zones with planetesimal surface densities well above the minimum-mass solar nebula. We present new simulations of the formation of Uranus and Neptune in the solid-rich disk of Dodson-Robinson et al. (Dodson-Robinson, S.E., Willacy, K., Bodenheimer, P., Turner, N.J., Beichman, C.A. [2009]. Icarus 200, 672-693) using the initial semimajor axis distribution of the Nice model (Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A. [2005]. Nature 435, 466-469; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R. [2005]. Nature 435, 462-465; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F. [2005]. Nature 435, 459-461), with one ice giant forming at 12 AU and the other at 15 AU. The innermost ice giant reaches its present mass after 3.8-4.0 Myr and the outermost after 5.3-6 Myr, a considerable time decrease from previous one-dimensional simulations (e.g. Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y. [1996]. Icarus 124, 62-85). The core masses stay subcritical, eliminating the need for a sudden gas accretion cutoff.Our calculated carbon mass fractions of 22% are in excellent agreement with the ice giant interior models of Podolak et al. (Podolak, M., Weizman, A., Marley, M. [1995]. Planet. Space Sci. 43, 1517-1522) and Marley et al. (Marley, M.S., Gómez, P., Podolak, M. [1995]. J. Geophys. Res. 100, 23349-23354). Based on the requirement that the ice giant-forming planetesimals contain >10% mass fractions of methane ice, we can reject any Solar System formation model that initially places Uranus and Neptune inside of Saturn’s orbit. We also demonstrate that a large population of planetesimals must be present in both ice giant feeding zones throughout the lifetime of the gaseous nebula. This research marks a substantial step forward in connecting both the dynamical and chemical aspects of planet formation. Although we cannot say that the solid-rich solar nebula model of Dodson-Robinson et al. (Dodson-Robinson, S.E., Willacy, K., Bodenheimer, P., Turner, N.J., Beichman, C.A. [2009]. Icarus 200, 672-693) gives exactly the appropriate initial conditions for planet formation, rigorous chemical and dynamical tests have at least revealed it to be a viable model of the early Solar System.