The first ten million years in the solar nebula

Abstract— The ²⁶Al/²⁷Al ratio in a large number of calcium-aluminum inclusions (CAIs) is a rather uniform 5 × 10^?5, whereas in chondrules the ratio is either undetectable or has a much lower value; the simplest interpretation of this is that there was an interval of a few million years between the times that these two meteoritic constituents formed stable solids. The present investigation was undertaken as an exploration of the physics of the processes in the solar nebula during and after the accumulation of the Sun. Understanding the time scales of events in this nebular model, to see if this would cast light on this apparent CAI to chondrule time interval, was the major motivation for the exploration. There were four stages in the history of the solar nebula; in stage 1, a fragment of an interstellar molecular cloud collapsed to form the Sun and solar nebula; in stage 2, the nebula was in approximate steady state balance between infall from the cloud and accretion onto the Sun and was in its FU Orionis accumulation stage; in stage 3, the Sun had been mainly accumulated and there was a slow residual mass flow into the Sun while it was in its classical T Tauri stage; and in stage 4, the nebula had finished accreting material onto the Sun (now a weak-lined T Tauri star) and was in a static condition with no significant dissipation or motions, other than removal at the inner edge due to the T Tauri solar wind and photoevaporation beyond 9 astronomical units (AU). It is found that the energy source keeping the nebula warm during stages 3 and 4 is recombination of ionized H in the ionized bipolar jets and the T Tauri coronal expansion solar wind. The parameters of the heating model were adjusted to locate the ice sublimation line at 5.2 AU. In this work, a nebular model is used with a surface density of 4.25 × 10³ gm/cm² at 1 AU and a variation with radial distance as the inverse first power. Under normal conditions in the nebula, there is a negative pressure gradient that provides partial radial support for the gas, which thus circles the Sun more slowly than large solid objects do. Large objects undergo a slow inward spiral due to the gas drag; very small objects move essentially with the gas but have a slow inward drift; and intermediate objects (e.g., 1 m) have a fairly large inward drift velocity that traverses the full radial extent of the nebula in considerably less than the CAI to chondrule time interval. Such objects are thus lost unless they can grow rapidly to larger sizes. Near the inner edge (bow) of the nebula during stage 4, the pressure gradient becomes positive, creating a narrow zone of zero gas drag toward which solids drift from both directions, facilitating planetesimal formation in the inner solar nebula. Recent theoretical and experimental results on sticking probabilities of solids show that icy surfaces have the best sticking properties, but icy interstellar grains can only stick together when subjected to impact velocities of less than 2000 cm/sec. However, if the solid objects are very underdense, then a collision leads to interpenetration and many points at which the small constituent grains can adhere to one another, and thus coagulation becomes possible for such underdense objects. Simulations were made of such coagulation in the outer solar nebula, and it was found that the central plane of the nebula quickly becomes filled with meter-sized and larger bodies that rapidly accumulated near the top of the nebula and rapidly descended; in a few thousand years this quickly leads to gravitational instabilities that can form planetesimals. These processes led to the rapid formation of Jupiter in the nebula (and the slightly less rapid formation of the other giant planets). The early formation of Jupiter opens an annular gap in the nebula, and thus a second region is created in the nebula with zero gas drag. It is concluded that CAIs were formed at the end of stage 2 of the nebula history and moved out into the nebula for long-term storage, and that most chondrules were formed by magnetic reconnection flares in the bow region of the nebula during stage 4, several million years later. Carbonaceous meteorites should be formed on the far side of the Jovian gap, with the chondrules being heated by flares on the early Jupiter irradiating materials in the nearby zone of zero gas drag, and they should have essentially the same ²⁶Al ages as the CAIs (this will be very hard to confirm owing to scarcity of Al mineral phases in these chondrules).     

The effects of nebula surface density profile and giant‐planet eccentricities on planetary accretion in the inner solar system

Abstract— We describe results of 32 N‐body planetary accretion simulations that investigate the dependence of terrestrial‐planet formation on nebula surface density profile σ and evolution of the eccentricities of Jupiter and Saturn e_j,s. Two surface density profiles are examined: a decaying profile with σ ∝ 1/a, where a is orbital semi‐major axis, and a peaked profile in which σ increases for a < 2 AU and decreases for a > 2 AU. The peaked profiles are generated by models of coagulation in an initially hot nebula. Models with initial e_j,s = 0.05 (the current value) and 0.1 are considered. Simulations using the decaying profile with e_j,s = 0.1 produce systems most like the observed planets in terms of mass‐weighted mean a and the absence of a planet in the asteroid belt. Simulations with doubled σ produce planets roughly twice as massive as the nominal case. Most initial embryos are removed in each simulation via ejection from the solar system or collision with the Sun. The asteroid belt is almost entirely cleared on a timescale of 10–100 Ma that depends sensitively on e_j,s. Most initial mass with a < 2 AU survives, with the degree of mass loss increasing with a. Mass loss from the terrestrial region occurs on a timescale that is long compared to the mass loss time for the asteroid belt. Substantial radial mixing of material occurs in all simulations, but is greater in simulations with initital e_j,s = 0.05. The degree of mixing is equivalent to a feeding zone of half width 1.5 and 0.9 AU for an Earth mass planet at 1 AU for the cases e_j,s = 0.05 and 0.1, respectively. In simulations with e_j,s = 0.05, roughly one‐third and 5–10% of the mass contained in final terrestrial planets originated in the region a > 2.5 AU for the decaying and peaked profiles, respectively. In the case e_j,s = 0.1, the median mass accreted from a > 2.5 AU is zero for both profiles.     

Outline of a theory of planet formation in binary systems

A theory is presented for determining regions where planets may form in binary star systems. Planet formation by accretion is assumed possible if mean planetesimal collision velocities do not exceed a critical value. Collision velocities are increased by perturbations due to the companion star, treated by secular perturbation theory. Collision velocities are damped by aerodynamic drag within the solar nebula, taken as the linear case of Cameron and Pine.A general feature of planetary systems in binary stars is the existence of two zones. The inner zone has enough damping to permit unimpeded growth by accretion; in the outer zone, growth proceeds to a limited diameter, beyond which damping is insufficient. It is shown that the asteroids could not have failed to coalesce due to Jupiter perturbations in the primitive solar nebula. Binary star systems with semimajor axis < 30AU are not likely to have planets; these include Alpha Centauri and 70 Ophiuchi. Systems possibly possessing planets include Eta Cassiopeiae, 40 Eridani, and Σ 2398. Epsilon Eridani is a marginal case.     

Long term dynamical evolution and classification of classical TNOs

Classical trans-Neptunian objects (TNOs) are believed to represent the most dynamically pristine population in the trans-Neptunian belt (TNB) offering unprecedented clues about the formation of our Solar System. The long term dynamical evolution of classical TNOs was investigated using extensive simulations. We followed the evolution of more than 17000 particles with a wide range of initial conditions taking into account the perturbations from the four giant planets for 4 Gyr. The evolution of objects in the classical region is dependent on both their inclination and semimajor axes, with the inner (a<45 AU) and outer regions (a>45 AU) evolving differently. The reason is the influence of overlapping secular resonances with Uranus and Neptune (40–42 AU) and the 5:3 (a∼ ∼42.3 AU), 7:4 (a∼ ∼43.7 AU), 9:5 (a∼ ∼44.5 AU) and 11:6 (a∼ ∼ 45.0 AU) mean motion resonances strongly sculpting the inner region, while in the outer region only the 2:1 mean motion resonance (a∼ ∼47.7 AU) causes important perturbations. In particular, we found: (a) A substantial erosion of low-i bodies (i<10°) in the inner region caused by the secular resonances, except those objects that remained protected inside mean motion resonances which survived for billion of years; (b) An optimal stable region located at 45 AU<a<47 AU, q>40 AU and i>5° free of major perturbations; (c) Better defined boundaries for the classical region: 42–47.5 AU (q>38 AU) for cold classical TNOs and 40–47.5 AU (q>35 AU) for hot ones, with i=4.5° as the best threshold to distinguish between both populations; (d) The high inclination TNOs seen in the 40–42 AU region reflect their initial conditions. Therefore they should be classified as hot classical TNOs. Lastly, we report a good match between our results and observations, indicating that the former can provide explanations and predictions for the orbital structure in the classical region.     

The Leonard Award Address Presented 1998 July 27, Dublin,Ireland: On the difficulties of making Earth-like planets

Abstract— Here I discuss the series of events that led to the formation and evolution of our planet to examine why the Earth is unique in the solar system. A multitude of factors are involved: These begin with the initial size and angular momentum of the fragment that separated from a molecular cloud; such random factors are crucial in determining whether a planetary system or a double star develops from the resulting nebula. Another requirement is that there must be an adequate concentration of heavy elements to provide the 2% “rock” and “ice” components of the original nebula. An essential step in forming rocky planets in the inner nebula is the loss of gas and depletion of volatile elements, due to early solar activity that is linked to the mass of the central star. The lifetime of the gaseous nebula controls the formation of gas giants. In our system, fine timing was needed to form the gas giant, Jupiter, before the gas in the nebula was depleted. Although Uranus and Neptune eventually formed cores large enough to capture gas, they missed out and ended as ice giants. The early formation of Jupiter is responsible for the existence of the asteroid belt (and our supply of meteorites) and the small size of Mars, whereas the gas giant now acts as a gravitational shield for the terrestrial planets. The Earth and the other inner planets accreted long after the giant planets, from volatile-depleted planetesimals that were probably already differentiated into metallic cores and silicate mantles in a gas-free, inner nebula. The accumulation of the Earth from such planetesimals was essentially a stochastic process, accounting for the differences among the four rocky inner planets—including the startling contrast between those two apparent twins, Earth and Venus. Impact history and accretion of a few more or less planetesimals were apparently crucial. The origin of the Moon by a single massive impact with a body larger than Mars accounts for the obliquity (and its stability) and spin of the Earth, in addition to explaining the angular momentum, orbital characteristics, and unique composition of the Moon. Plate tectonics (unique among the terrestrial planets) led to the development of the continental crust on the Earth, an essential platform for the evolution of Homo sapiens. Random major impacts have punctuated the geological record, accentuating the directionless course of evolution. Thus a massive asteroidal impact terminated the Cretaceous Period, resulted in the extinction of at least 70% of species living at that time, and led to the rise of mammals. This sequence of events that resulted in the formation and evolution of our planet were thus unique within our system. The individual nature of the eight planets is repeated among the 60-odd satellites—no two appear identical. This survey of our solar system raises the question whether the random sequence of events that led to the formation of the Earth are likely to be repeated in detail elsewhere. Preliminary evidence from the “new planets” is not reassuring. The discovery of other planetary systems has removed the previous belief that they would consist of a central star surrounded by an inner zone of rocky planets and an outer zone of giant planets beyond a few astronomical units (AU). Jupiter-sized bodies in close orbits around other stars probably formed in a similar manner to our giant planets at several astronomical units from their parent star and, subsequently, migrated inwards becoming stranded in close but stable orbits as “hot Jupiters”, when the nebula gas was depleted. Such events would prevent the formation of terrestrial-type planets in such systems.     

The positions of secular resonance surfaces

The surfaces for the three strongest secular resonances have been located as a function of proper semimajor axis, eccentricity, and inclination for semimajor axes between 1.25 and 3.5 AU. The results are presented graphically. The ν₅ resonance only occurs at high inclinations (?23°). The ν₆ resonance passes through both the main belt and Mars-crossing space. The ν₁₆ resonance starts near the inner edge of the belt and, at low inclinations at least, folds around a portion of the Mars-crossing space until it runs nearly parallel with the Earth-crossing boundary.     

Mean Motion Resonances, Gas Drag, and Supersonic Planetesimals in the Solar Nebula

We examine the orbital evolution of planetesimals under the influence of Jupiter's perturbations and nebular gas drag, under the assumption that gas persisted in the asteroid region for some time after Jupiter attained its final mass. Two distinct mechanisms, associated with the 2 : 1 and 3 : 2 mean motion resonances, can excite eccentricities to high values, despite the damping effect of drag. If Jupiter's eccentricity was comparable to its present value, planetesimals can be temporarily trapped in the 2 : 1 resonance. Bodies crossing the 3 : 2 resonance can enter a region of phase space with overlapping high-order resonances. Both mechanisms can produce eccentricities greater than 0.5 for asteroid-sized planetesimals. The combination of resonant perturbations and drag causes secular decay of semimajor axes, resulting in migration of bodies from the outer to inner belt. Inclinations remain low, implying significant collisional evolution during this migration. Velocities of resonant bodies relative to the gas are highly supersonic; these would have been a source of shock waves in the solar nebula.This revised version was published online in October 2005 with corrections to the Cover Date.     

Astronomical constraints on nebular temperatures: Implications for planetesimal formation

Abstract— Motivated by recent observations of T-Tauri stars and the interpretation of these observations in terms of the properties of circumstellar disks, we derive internal (midplane) temperatures for disks around mature (age ～1 Ma) T-Tauri stars. The estimates are obtained by combining published results for disk masses, sizes, accretion rates, and surface temperatures. For 26 stars (for which adequate data are available), we derive midplane temperatures at 1 AU primarily in the range 200–800 K, and 100–400 K at 2.5 AU. It is likely that the solar nebula, at the same stage of evolution, contained planetesimals and objects destined to become meteorite parent bodies. Observations of young stellar objects at earlier stages of evolution (age ～0.1 Ma) imply that accretion rates were, on the average, at least two orders of magnitude greater than the 10^?8 M_⊙/year rates typical for mature T-Tauri stars. Such high values would result in midplane temperatures at or near the silicate vaporization temperature in the terrestrial planet region. If cooling of the solar nebula from such a hot epoch was responsible for establishing the pervasive elemental fractionation patterns found in chondritic meteorites, then objects in the asteroid belt must have grown rapidly (within 0.1 Ma) to sizes of ～1 km, a conclusion consistent with current theories of planetesimal formation. However, the fact that primitive meteorite parent bodies escaped being melted by the decay of ²⁶Al then implies that further growth of at least some objects was essentially delayed for 2 Ma or more. Such a diminished growth rate appears to be consistent with simulations of the dynamics of solid bodies in the asteroid belt. Other hypotheses seem less attractive. One might assume that the final cooling occurred only after the decay of ²⁶Al (i.e., more than a million years after calcium-aluminum rich inclusion formation), or that ²⁶Al was not ubiquitous in the early solar system. But the first of these conjectures is incompatible with astronomical observations of T-Tauri systems, and the second appears to be contradicted by the evidence for ²⁶Al in diverse meteoritic components. The remaining alternative would then appear to be that, despite a lack of supporting evidence, chondritic fractionation patterns reflect the net effect of many local heating and cooling events and have nothing to do with global nebular cooling. We conclude that the most plausible hypothesis is that both nebular cooling and coagulation of solids to kilometer-sized objects occurred rapidly and that a substantial number of planetesimals in the asteroid belt remained smaller than a few kilometers in radius for at least 2 Ma.     

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.     

Asteroid sizes and albedos

The radiometric method of determining diameters of asteroids is reviewed, and a synthesis of radiometric and polarimetric measurements of the diameters and geometric albedos of a total of 187 asteroids is presented. All asteroids with diameters greater than 250 km are identified, and statistical studies can be carried out of the size distributions of different albedo classes down to 80-km diameter for the entire main asteroid belt (2.0–3.5 AU). The distribution of albedos is strongly bimodal, with mean albedos for the C and S groups of 0.035 and 0.15, respectively. The C asteroids outnumber the S at all sizes and all values of semi-major axis, increasing from a little over half the population inside 2.5 AU to more than 95% beyond 3.0 AU; for all objects with D > 70 km, the ratio C/(C+S) is 0.88 ± 0.04. More than half of all asteroids in this size range have a > 3.0 AU. The M asteroids constitute about 5% of the population for a < 3.0 AU, but no members of of this class have been identified in the outer belt. There are no significant differences between the distributions of C, S, and M asteroids for the largest asteroids (D > 200 km) and for those of intermediate size (200–270 km). The total mass in the belt, down to 70-km size, but excluding Ceres, is about 2 × 10²⁴ g. Evidence is presented that several large asteroids rotate in a prograde sense, and that a real difference existsbetween the bulk densities of Ceres and Vesta.     

Evolution of the Oort cloud angular momentum: Numerical simulation

This paper considers the evolution of a flat svarm of cometary bodies (under the effect of the passage of stars), initially moving in one direction along the circular orbits with radii 1.4×10⁴<r<2×10⁴ AU and along elliptic orbits with semi-major axes 5×10³<a<1×10⁴ AU and with perihelia within 50<q<100 AU. Numerical simulation shows that the original flat belt of comets is thermalizing. Its root-mean-squarez-coordinate grows withr. A cometary cloud forms with a dense flattened inner core and a rarefied halo (the Oort cloud proper). The value =N _core/N _halo varies within a wide range (up to the order of magnitude) depending on the model used (N _core andN _halo are the numbers of comets in the core and the halo, respectively).The hypothesis of a massive Oort cloud (Marochniket al., 1988) implies that the Oort cloud should have a large angular momentum. This paper employs numerical simulation to calculate Oort cloud models to which the initially flat located at the periphery of the solar nebula rotating cometary swarms is evolving in time. The loss of the initial angular momentum over the time of the Oort cloud evolution is not large.     

Motion of dust particles under the influence of a latitudinally dependent force

The Kelperian motion of dust particles in the solar system is mainly influenced by the electromagnetic and plasma Poynting-Robertson drag. The first force is isotropic while the second one shows latitudinal variations due to the observed differences of the solar wind parameters in the ecliptic plane and over the solar poles. Close to the Sun other effects become important, e.g. sublimation and sputtering, as well as for submicron particles Lorentz scattering has to be taken into account. These forces are very weak for dust grains of moderate size (10–100 µ) not too close (>0.03 AU) to the Sun and are neglected here. Assuming that the general form of the latidudinally dependent force is a series expansion in Legendre polynomials, we have studied the averaged equations of motion for the classical elements and found the first integral of them. The general character of motion is the same as for the classical Poynting-Robertson drag: particles spiral towards the Sun. The new features in the orbital evolution under the latitudinally dependent force as compared with the isotropic Poynting-Robertson drag are:

1. not only the semimajor axisa and the eccentricity ε but also the argument of the perihelion ω varies with time,
2. the rate of change ofa, ε, ω depends on the inclination.

An example of particle trajectories in the phase space of elements is presented.     

The evidence of an early stellar encounter in Edgeworth-Kuiper belt

We investigate the influence of a stellar fly-by encounter on the Edgeworth-Kuiper belt objects through numerical orbital calculations, in order to explain both mass depletion and high orbital inclinations of the classical Edgeworth-Kuiper belt (CEKB) objects, which have semimajor axis of 42-48 AU and perihelia beyond 35 AU. The observationally inferred total mass of the CEKB is ∼1/10 Earth masses, which is only ∼0.02 of that extrapolated from the minimum-mass solar nebula model. The CEKB consists of bimodal population: “hot population” with inclinations i?0.2-0.6 radians and “cold population” with i?0.1. The observationally suggested difference in size and color of objects between the two populations may imply different origins of the two populations. We find that both the depletion of solid materials in the CEKB and the formation of the hot population are accounted for by a single close stellar encounter with pericenter distance of 80-100 AU and inclination relative to the initial protoplanetary disk ?50°-70°. Such a stellar encounter highly pumps up eccentricities of most objects in the CEKB and then their perihelia migrate within 35 AU. These objects would be removed by Neptune's perturbations after Neptune is formed at or migrates to the current position (30 AU). Less than 10% of the original objects remain in stable orbits with small eccentricities and perihelion distances larger than 35 AU, in the CEKB, which is consistent with the observation. We find that i of the remaining objects are as large as that of the observed hot population. The only problem is how to stop Neptune's migration at ∼30 AU, which is addressed in a separate paper. The depletion by the stellar encounter extends deeply into ∼30-35 AU, which provides the basis of the formation model for the cold population through Neptune's outward migration by Levison and Morbidelli (2003, Nature, 426, 419-421). The combination of our model with Levison and Morbidelli's model could consistently explain the mass depletion, truncation at 50 AU, bimodal distribution in i, and differences in size and color between the hot and the cold populations in the CEKB.     

Solar nebula dispersal and the stability of the planetary system: I. Scanning secular resonance theory

Secular resonances in the early solar system are studied in an effort to establish constraints on the time scale and/or method of solar nebula dispersal. Simplified nebula models and dispersal routines are employed to approximate changes in an assumed axisymmetric nebula potential. These changes, in turn, drive an evolutionary sequence of Laplace-Lagrange solutions for the secular variations of the solar system. A general feature of these sequences is a sweep of one or more giant planet resonances through the inner solar system. Their effect is rate dependent; in the linearized models considered, characteristic dispersal times ≤O(10^4?5 years) are required to avoid the generation of terrestrial eccentricities and inclinations in excess of observed values. These times are short compared to typical estimates of the accretion time scales [i.e., ～O(10^7?9 years)] and may provide an important boundary condition for developing models of nebula dispersal and solar system formation in general.     

Kuiper带天体的轨道分布特性

1992年9月,夏威夷大学的D.Jewitt和加利福尼亚大学的J.Lun发现了海王星外绕太阳运行的第一个小天体1992QB1[1],开创了人类对于海王星外天体的实际观测的研究.近10年的接连不断发现,已经证实了海王星轨道外面存在着一个由大量的环绕太阳运动的小天体组成的环带[2].由于G.P.Kuiper曾在1951年的文章中提出过在冥王星的外边可能存在小天体的问题,因此人们一般把这个环带称为Kuiper带,你这些天体为“KuiperBelt　Objects”(KBOs),或从逻辑上称它们为“Trans-NeptunianObjects”(TNOs)[3]     

The Villalbeto de la Peña meteorite fall: II. Determination of atmospheric trajectory and orbit

Abstract— The L6 ordinary chondrite Villalbeto de la Peña fall occurred on January 4, 2004, at 16: 46: 45 ± 2 s UTC. The related daylight fireball was witnessed by thousands of people from Spain, Portugal, and southern France, and was also photographed and videotaped from different locations of León and Palencia provinces in Spain. From accurate astrometric calibrations of these records, we have determined the atmospheric trajectory of the meteoroid. The initial fireball velocity, calculated from measurements of 86 video frames, was 16.9 ± 0.4 km/s. The slope of the trajectory was 29.0 ± 0.6° to the horizontal, the recorded velocity during the main fragmentation at a height of 27.9 ± 0.4 km was 14.2 ± 0.2 km/s, and the fireball terminal height was 22.2 ± 0.2 km. The heliocentric orbit of the meteoroid resided in the ecliptic plane (i = 0.0 ± 0.2°), having a perihelion distance of 0.860 ± 0.007 AU and a semimajor axis of 2.3 ± 0.2 AU. Therefore, the meteorite progenitor body came from the Main Belt, like all previous determined meteorite orbits. The Villalbeto de la Peña fireball analysis has provided the ninth known orbit of a meteorite in the solar system.     

Asteroidal source of L chondrite meteorites

Establishing connections between meteorites and their parent asteroids is an important goal of planetary science. Several links have been proposed in the past, including a spectroscopic match between basaltic meteorites and (4) Vesta, that are helping scientists understand the formation and evolution of the Solar System bodies. Here we show that the shocked L chondrite meteorites, which represent about two thirds of all L chondrite falls, may be fragments of a disrupted asteroid with orbital semimajor axis a=2.8 AU. This breakup left behind thousands of identified 1–15 km asteroid fragments known as the Gefion family. Fossil L chondrite meteorites and iridium enrichment found in an ≈467 Ma old marine limestone quarry in southern Sweden, and perhaps also ∼5 large terrestrial craters with corresponding radiometric ages, may be tracing the immediate aftermath of the family-forming collision when numerous Gefion fragments evolved into the Earth-crossing orbits by the 5:2 resonance with Jupiter. This work has major implications for our understanding of the source regions of ordinary chondrite meteorites because it implies that they can sample more distant asteroid material than was previously thought possible.     

The planetesimal bow shock model for chondrule formation: A more quantitative assessment of the standard (fixed Jupiter) case

Abstract– One transient heating mechanism that can potentially explain the formation of most meteoritic chondrules 1–3 Myr after CAIs is shock waves produced by planetary embryos perturbed into eccentric orbits via resonances with Jupiter following its formation. The mechanism includes both bow shocks upstream of resonant bodies and impact vapor plume shocks produced by high‐velocity collisions of the embryos with small nonresonant planetesimals. Here, we investigate the efficiency of both shock processes using an improved planetesimal accretion and orbital evolution code together with previous simulations of vapor plume expansion in the nebula. Only the standard version of the model (with Jupiter assumed to have its present semimajor axis and eccentricity) is considered. After several hundred thousand years of integration time, about 4–5% of remaining embryos have eccentricities greater than about 0.33 and shock velocities at 3 AU exceeding 6 km s^?1, currently considered to be a minimum for melting submillimeter‐sized silicate precursors in bow shocks. Most embryos perturbed into highly eccentric orbits are relatively large—half as large as the Moon or larger. Bodies of this size could yield chondrule‐cooling rates during bow shock passage compatible with furnace experiment results. The cumulative area of the midplane that would be traversed by highly eccentric embryos and their associated bow shocks over a period of 1–2 Myr is <1% of the total area. However, future simulations that consider a radially migrating Jupiter and alternate initial distributions of the planetesimal swarm may yield higher efficiencies.     

V-type asteroids in the middle main belt

V-type asteroids are bodies whose surfaces are constituted of basalt. In the Main Asteroid Belt, most of these asteroids are assumed to come from the basaltic crust of Asteroid (4) Vesta. This idea is mainly supported by (i) the fact that almost all the known V-type asteroids are in the same region of the belt as (4) Vesta, i.e., the inner belt (semi-major axis 2.1<a<2.5 AU), (ii) the existence of a dynamical asteroid family associated to (4) Vesta, and (iii) the observational evidence of at least one large craterization event on Vesta's surface. One V-type asteroid that is difficult to fit in this scenario is (1459) Magnya, located in the outer asteroid belt, i.e., too far away from (4) Vesta as to have a real possibility of coming from it. The recent discovery of the first V-type asteroid in the middle belt (2.5<a<2.8 AU), (21238) 1995WV7 [Binzel, R.P., Masi, G., Foglia, S., 2006. Bull. Am. Astron. Soc. 38, 627; Hammergren, M., Gyuk, G., Puckett, A., 2006. ArXiv e-print, astro-ph/0609420], located at ∼2.54 AU, raises the question of whether it came from (4) Vesta or not. In this paper, we present spectroscopic observations indicating the existence of another V-type asteroid at ∼2.53 AU, (40521) 1999RL95, and we investigate the possibility that these two asteroids evolved from the Vesta family to their present orbits by a semi-major axis drift due to the Yarkovsky effect. The main problem with this scenario is that the asteroids need to cross the 3/1 mean motion resonance with Jupiter, which is highly unstable. Combining N-body numerical simulations of the orbital evolution, that include the Yarkovsky effect, with Monte Carlo models, we compute the probability that an asteroid of a given diameter D evolves from the Vesta family and crosses over the 3/1 resonance, reaching a stable orbit in the middle belt. Our results indicate that an asteroid like (21238) 1995WV7 has a low probability (∼1%) of having evolved through this mechanism due to its large size (D∼5 km), because the Yarkovsky effect is not sufficiently efficient for such large asteroids. However, the mechanism might explain the orbits of smaller bodies like (40521) 1999RL95 (D∼3 km) with ∼70-100% probability, provided that we assume that the Vesta family formed ?3.5 Gy ago. We estimate the debiased population of V-type asteroids that might exist in the same region as (21238) and (40521) (2.5<a?2.62 AU) and conclude that about 10 to 30% of the V-type bodies with D>1 km may come from the Vesta family by crossing over the 3/1 resonance. The remaining 70-90% must have a different origin.     

Dust to planetesimals: Settling and coagulation in the solar nebula

The behavior of solid particles in a low-mass solar nebula during settling to the central plane and the formation of planetesimals is examined. Gravitational instability in a dust layer and collisional accretion are considered as possible mechanisms of planetesimal formation. Non-Keplerian rotation of the nebula results in shear between the gas and a dust layer. This shear produces turbulence within the layer which inhibits gravitational instability, unless the mean particle size exceeds a critical value, ～1 cm at 1 AU. The size requirement is less stringent at larger heliocentric distances, suggesting a possible difference in planetesimal formation mechanisms between the inner and outer nebula. Coagulation of grains during settling is expected in the solar nebula environment. Van der Waals forces appear adequate to produce centimeter-sized aggregates. Growth is primarily due to sweepup of small particles by larger ones due to size-dependent settling velocities. A numerical model for computing simultaneous coagulation and settling is described. Relative velocities are determined by gas drag and the non-Keplerian rotation of the nebula. The settling is very nonhomologous. Most of the solid matter reaches the central plane as centimeter-sized aggregates in a few times 10³ revolutions, but some remains suspended in the form of fine dust. Drag-induced relative velocities result in collisions. The growth of bodies in the central plane is initially rapid. After sizes reach ～10³ cm, relative velocities decrease and the growth rate declines. Gas drag rapidly damps the out-of-plane motions of these intermediate-sized bodies. They settle into a thin layer which is subject to gravitational instability. Kilometer-sized planetesimals are formed by this composite process.