Accretion rates of planetesimals by protoplanets embedded in nebular gas

When protoplanets growing by accretion of planetesimals have atmospheres, small planetesimals approaching the protoplanets lose their energy by gas drag from the atmospheres, which leads them to be captured within the Hill sphere of the protoplanets. As a result, growth rates of the protoplanets are enhanced. In order to study the effect of an atmosphere on planetary growth rates, we performed numerical integration of orbits of planetesimals for a wide range of orbital elements and obtained the effective accretion rates of planetesimals onto planets that have atmospheres. Numerical results are obtained as a function of planetesimals’ eccentricity, inclination, planet’s radius, and non-dimensional gas-drag parameters which can be expressed by several physical quantities such as the radius of planetesimals and the mass of the protoplanet. Assuming that the radial distribution of the gas density near the surface can be approximated by a power-law, we performed analytic calculation for the loss of planetesimals’ kinetic energy due to gas drag, and confirmed agreement with numerical results. We confirmed that the above approximation of the power-law density distribution is reasonable for accretion rate of protoplanets with 1-10 Earth masses, unless the size of planetesimals is too small. We also calculated the accretion rates of planetesimals averaged over a Rayleigh distribution of eccentricities and inclinations, and derived a semi-analytical formula of accretion rates, which reproduces the numerical results very well. Using the obtained expression of the accretion rate, we examined the growth of protoplanets in nebular gas. We found that the effect of atmospheric gas drag can enhance the growth rate significantly, depending on the size of planetesimals.     

Impact angle influence in high velocity dust collisions during planetesimal formation

We have examined the influence of impact angle in collisions between small dust aggregates and larger dust targets through laboratory experiments. Targets consisted of μm-sized quartz dust and had a porosity of about 67%; the projectiles, between 1 and 5 mm in diameter, were slightly more compact (64% porosity). The collision velocity was centered at 20 m/s and impact angles range from 0° to 45°. At a given impact angle, the target gained mass for projectiles smaller than a threshold size, which decreases with increasing angle from about 3 mm to 1 mm. The fact that growth is possible up to the largest angles studied supports the idea of planetesimal formation by sweep-up of small dust aggregates.     

Planetesimal formation by turbulent concentration

The formation of 1-1000 km diameter planetesimals from dust grains in a protoplanetary disk is a key step in planet formation. Conventional models for planetesimal formation involve pairwise sticking of dust grains, or the sedimentation of dust grains to a thin layer at the disk midplane followed by gravitational instability. Each of these mechanisms is likely to be frustrated if the disk is turbulent. Particles with stopping times comparable to the turnover time of the smallest eddies in a turbulent disk can become concentrated into dense clumps that may be the precursors of planetesimals. Such particles are roughly millimeter-sized for a typical protoplanetary disk. To survive to become planetesimals, clumps need to form in regions of low vorticity to avoid rotational breakup. In addition, clumps must have sufficient self gravity to avoid break up due to the ram pressure of the surrounding gas. Given these constraints, the rate of planetesimal formation can be estimated using a cascade model for the distribution of particle concentration and vorticity within eddies of various sizes in a turbulent disk. We estimate planetesimal formation rates and planetesimal diameters as a function of distance from a star for a range of protoplanetary disk parameters. For material with a solar composition, the dust-to-gas ratio is too low to allow efficient planetesimal formation, and most solid material will remain in small particles. Enhancement of the dust-to-gas ratio by 1-2 orders of magnitude, either vertically or radially, allows most solid material to be converted into planetesimals within the typical lifetime of a disk. Such dust-to-gas ratios may occur near the disk midplane as a result of vertical settling of short-lived clumps prior to clump breakup. Planetesimal formation rates are sensitive to the assumed size and rotational speed of the largest eddies in the disk, and formation rates increase substantially if the largest eddies rotate more slowly than the disk itself. Planetesimal formation becomes more efficient with increasing distance from the star unless the disk surface density profile has a slope of −1.5 or steeper as a function of distance. Planetesimal formation rates typically increase by an order-of-magnitude or more moving outward across the snow line for a solid surface density increase of a factor of 2. In all cases considered, the modal planetesimal size increases with roughly the square root of distance from the star. Typical modal diameters are 100 km and 400 km in the regions corresponding to the asteroid belt and Kuiper belt in the Solar System, respectively.     

Initial sizes of planetesimals and accretion of the asteroids

The present size frequency distribution (SFD) of bodies in the asteroid belt appears to have preserved some record of the primordial population, with an excess of bodies of diameter D ∼ 100 km relative to a simple power law. The survival of Vesta’s basaltic crust also implies that the early SFD had a shallow slope in the range ∼10-100 km. (Morbidelli, A., Bottke, W.F., Nesvorny, D., Levison, H.F. [2009]. Icarus 204, 558-573) were unable to produce these features by accretion from an initial population of km-sized planetesimals. They concluded that bodies with sizes in the range ∼100-1000 km and a SFD similar to the current population were produced directly from solid particles of sub-meter scale, without experiencing accretion through intermediate sizes. We present results of new accretion simulations in the primordial asteroid region. The requisite SFD can be produced from an initial population of planetesimals of sizes ?0.1 km, smaller than the usual assumption of km-sized bodies. The bump at D ∼ 100 km is produced by a transition from dispersion-dominated runaway growth to a regime dominated by Keplerian shear, before the formation of large protoplanetary embryos. Thus, accretion of the asteroids from an initial population of small (sub-km) planetesimals cannot be ruled out.     

The distributions and ages of refractory objects in the solar nebula

Refractory objects such as Calcium, Aluminum-rich Inclusions, Amoeboid Olivine Aggregates, and crystalline silicates, are found in primitive bodies throughout our Solar System. It is believed that these objects formed in the hot, inner solar nebula and were redistributed during the mass and angular momentum transport that took place during its early evolution. The ages of these objects thus offer possible clues about the timing and duration of this transport. Here we study how the dynamics of these refractory objects in the evolving solar nebula affected the age distribution of the grains that were available to be incorporated into planetesimals throughout the Solar System. It is found that while the high temperatures and conditions needed to form these refractory objects may have persisted for millions of years, it is those objects that formed in the first 10⁵ years that dominate (make up over 90%) those that survive throughout most of the nebula. This is due to two effects: (1) the largest numbers of refractory grains are formed at this time period, as the disk is rapidly drained of mass during subsequent evolution and (2) the initially rapid spreading of the disk due to angular momentum transport helps preserve this early generation of grains as opposed to later generations. This implies that most refractory objects found in meteorites and comets formed in the first 10⁵ years after the nebula formed. As these objects contained live ²⁶Al, this constrains the time when short-lived radionuclides were introduced to the Solar System to no later than 10⁵ years after the nebula formed. Further, this implies that the t=0 as defined by meteoritic materials represents at most, the instant when the solar nebula finished accreting significant amounts of materials from its parent molecular cloud.     

Radial diffusion rate of planetesimals due to gravitational encounters

We study the rate of radial diffusion of planetesimals due to mutual gravitational encounters under Hill’s approximations in the three-body problem. Planetesimals orbiting a central star radially migrate inward and outward as a result of mutual gravitational encounters and transfer angular momentum. We calculate the viscosity in a disk of equal-sized planetesimals due to their mutual gravitational encounters using three-body orbital integrations, and obtain a semianalytic expression that reproduces the numerical results. We find that the viscosity is independent of the velocity dispersion of planetesimals when the velocity dispersion is so small that Kepler shear dominates planetesimals’ relative velocities. On the other hand, in high-velocity cases where random velocities dominate the relative velocities, the viscosity is a decreasing function of the velocity dispersion, and is found to agree with previous estimates under the two-body approximation neglecting the solar gravity. We also calculate the rate of radial diffusion of planetesimals due to gravitational scattering by a massive protoplanet. Using these results, we discuss a condition for formation of nonuniform radial surface density distribution of planetesimals by gravitational perturbation of an embedded protoplanet.     

Accretion among preplanetary bodies: The many faces of runaway growth

When preplanetary bodies reach proportions of ∼1 km or larger in size, their accretion rate is enhanced due to gravitational focusing (GF). We have developed a new numerical model to calculate the collisional evolution of the gravitationally-enhanced growth stage. The numerical model is novel as it attempts to preserve the individual particle nature of the bodies (like N-body codes); yet it is statistical in nature since it must incorporate the very large number of planetesimals. We validate our approach against existing N-body and statistical codes. Using the numerical model, we explore the characteristics of the runaway growth and the oligarchic growth accretion phases starting from an initial population of single planetesimal radius R₀. In models where the initial random velocity dispersion (as derived from their eccentricity) starts out below the escape speed of the planetesimal bodies, the system experiences runaway growth. We associate the initial runaway growth phase with increasing GF-factors for the largest body. We find that during the runaway growth phase the size distribution remains continuous but evolves into a power-law at the high-mass end, consistent with previous studies. Furthermore, we find that the largest body accretes from all mass bins; a simple two-component approximation is inapplicable during this stage. However, with growth the runaway body stirs up the random motions of the planetesimal population from which it is accreting. Ultimately, this feedback stops the fast growth and the system passes into oligarchy, where competitor bodies from neighboring zones catch up in terms of mass. We identify the peak of GF with the transition between the runaway growth and oligarchy accretion stages. Compared to previous estimates, we find that the system leaves the runaway growth phase at a somewhat larger radius, especially at the outer disk. Furthermore, we assess the relevance of small, single-size fragments on the growth process. In classical models, where the initial velocity dispersion of bodies is small, these do not play a critical role during the runaway growth; however, in models that are characterized by large initial relative velocities due to external stirring of their random motions, a situation can emerge where fragments dominate the accretion, which could lead to a very fast growth.     

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.     

The origin of water on Mars

This paper considers the origin of water on Mars, in the context of a dynamical model that accounts for most of the Earth's water as a product of collisions between the growing Earth and planet-sized “embryos” from the asteroid belt. Mars' history is found to be different; to explain the present mass of Mars requires that it suffer essentially no giant collisions and the bulk of its growth is through addition of smaller bodies. Asteroids and comets from beyond 2.5 AU provide the source of Mars' water, which totals 6-27% of the Earth's present ocean (1 Earth ocean≡1.5×10²¹ kg), equivalent to 600-2700-m depth on the martian surface. The D/H ratio of this material is 1.2-1.6 times Standard Mean Ocean Water, the smaller value obtaining for the larger amount of water accreted. The upper half of the range of total water accreted, while many times less than that acquired by the Earth, is consistent with geological data on Mars, and the D/H value is that derived for martian magmatic water from SNC meteorites. Both together are consistent with published interpretations of the high D/H in present-day martian atmospheric water in terms of water loss through atmospheric escape.     

Towards initial mass functions for asteroids and Kuiper Belt Objects

Our goal is to understand primary accretion of the first planetesimals. Some examples are seen today in the asteroid belt, providing the parent bodies for the primitive meteorites. The primitive meteorite record suggests that sizeable planetesimals formed over a period longer than a million years, each of which being composed entirely of an unusual, but homogeneous, mixture of millimeter-size particles. We sketch a scenario that might help explain how this occurred, in which primary accretion of 10-100 km size planetesimals proceeds directly, if sporadically, from aerodynamically-sorted millimeter-size particles (generically “chondrules”). These planetesimal sizes are in general agreement with the currently observed asteroid mass peak near 100 km diameter, which has been identified as a “fossil” property of the pre-erosion, pre-depletion population. We extend our primary accretion theory to make predictions for outer Solar System planetesimals, which may also have a preferred size in the 100 km diameter range. We estimate formation rates of planetesimals and explore parameter space to assess the conditions needed to match estimates of both asteroid and Kuiper Belt Object (KBO) formation rates. For parameters that satisfy observed mass accretion rates of Myr-old protoplanetary nebulae, the scenario is roughly consistent with not only the “fossil” sizes of the asteroids, and their estimated production rates, but also with the observed spread in formation ages of chondrules in a given chondrite, and with a tolerably small radial diffusive mixing during this time between formation and accretion. As previously noted, the model naturally helps explain the peculiar size distribution of chondrules within such objects. The optimum range of parameters, however, represents a higher gas density and fractional abundance of solids, and a smaller difference between Keplerian and pressure-supported orbital velocities, than “canonical” models of the solar nebula. We discuss several potential explanations for these differences. The scenario also produces 10-100 km diameter primary KBOs, and also requires an enhanced abundance of solids to match the mass production rate estimates for KBOs (and presumably the planetesimal precursors of the ice giants themselves). We discuss the advantages and plausibility of the scenario, outstanding issues, and future directions of research.     

Thermal conductivity measurements of porous dust aggregates: I. Technique, model and first results

We present a non-invasive technique for measuring the thermal conductivity of fragile and sensitive materials. In the context of planet-formation research, the investigation of the thermal conductivity of porous dust aggregates provide important knowledge about the influence of heating processes, like internal heating by radioactive decay of short-lived nuclei, e.g. ²⁶Al, on the evolution and growth of planetesimals. The determination of the thermal conductivity was performed by a combination of laboratory experiments and numerical simulations. An IR camera measured the temperature distribution of the sample surface heated by a well-characterized laser beam. The thermal conductivity as free parameter in the model calculations, exactly emulating the experiment, was varied until the experimental and numerical temperature distributions showed best agreement. Thus, we determined for three types of porous dust samples, consisting of spherical, 1.5 μm-sized SiO₂ particles, with volume filling factors in the range of 15-54%, the thermal conductivity to be 0.002-0.02 W m⁻¹ K⁻¹, respectively. From our results, we can conclude that the thermal conductivity mainly depends on the volume filling factor. Further investigations, which are planned for different materials and varied contact area sizes (produced by sintering), will prove the appropriate dependencies in more detail.     

Planetary growth with collisional fragmentation and gas drag

As planetary embryos grow, gravitational stirring of planetesimals by embryos strongly enhances random velocities of planetesimals and makes collisions between planetesimals destructive. The resulting fragments are ground down by successive collisions. Eventually the smallest fragments are removed by the inward drift due to gas drag. Therefore, the collisional disruption depletes the planetesimal disk and inhibits embryo growth. We provide analytical formulae for the final masses of planetary embryos, taking into account planetesimal depletion due to collisional disruption. Furthermore, we perform the statistical simulations for embryo growth (which excellently reproduce results of direct N-body simulations if disruption is neglected). These analytical formulae are consistent with the outcome of our statistical simulations. Our results indicate that the final embryo mass at several AU in the minimum-mass solar nebula can reach about ∼0.1 Earth mass within 10⁷ years. This brings another difficulty in formation of gas giant planets, which requires cores with ∼10 Earth masses for gas accretion. However, if the nebular disk is 10 times more massive than the minimum-mass solar nebula and the initial planetesimal size is larger than 100 km, as suggested by some models of planetesimal formation, the final embryo mass reaches about 10 Earth masses at 3-4 AU. The enhancement of embryos’ collisional cross sections by their atmosphere could further increase their final mass to form gas giant planets at 5-10 AU in the Solar System.     

Oligarchic growth with migration and fragmentation

In the core-accretion model, giant-planet cores form by oligarchic growth from a population of planetesimals prior to the dispersal of the disk gas. Once a core reaches a critical mass of roughly 10 Earth masses, it begins to accrete a gaseous envelope, forming a giant planet. Collisions between planetesimals cause fragmentation. Planetesimal fragments are more easily captured by cores, speeding up growth, but fragments are also lost by radial drift, reducing the total solid mass in the disk. Interaction with the gas causes cores to undergo inward type-I migration. Migration allows a core to accrete planetesimals from a larger region, but migrating cores may be lost if they reach the star. Thus, migration and fragmentation have both a positive and a negative impact on core formation. Here we describe results of new simulations of oligarchic growth that include fragmentation and/or migration. In the absence of migration, cores grow until they reach their isolation mass, which increases with distance from the star, or until the disk gas disperses. Fragmentation increases the maximum core mass by increasing growth rates in the outer disk, allowing objects to reach their isolation mass during the disk lifetime. When migration is present, cores migrate inwards rapidly when they approach 1 Earth mass. Most migrating cores are lost. Migrating cores gain little extra mass since they are passing through regions that have been depleted by earlier generations of cores. For a disk viscosity parameter alpha=1e−3 and planetesimal radius = 10 km, the maximum core mass is roughly 4 and 0.5 Earth masses with/without fragmentation, respectively, with little dependence on the disk mass. Formation and survival of 10-Earth-mass cores, in the presence of migration, requires large alpha (1e−2) and a massive disk (0.1 solar masses). When alpha is large, type-I migration rates decrease rapidly with time, allowing large, late-forming cores to survive. The addition of a stochastic (random-walk) migration component makes little difference to the outcome, provided that stochastic migration affects only cores larger than 0.01 Earth masses. Stochastic migration becomes increasingly important if it also affects lower-mass objects.     

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.     

N-Body simulations of growth from 1 km planetesimals at 0.4 AU

We present N-body simulations of planetary accretion beginning with 1 km radius planetesimals in orbit about a 1 M_⊙ star at 0.4 AU. The initial disk of planetesimals contains too many bodies for any current N-body code to integrate; therefore, we model a sample patch of the disk. Although this greatly reduces the number of bodies, we still track in excess of 10⁵ particles. We consider three initial velocity distributions and monitor the growth of the planetesimals. The masses of some particles increase by more than a factor of 100. Additionally, the escape speed of the largest particle grows considerably faster than the velocity dispersion of the particles, suggesting impending runaway growth, although no particle grows large enough to detach itself from the power law size-frequency distribution. These results are in general agreement with previous statistical and analytical results. We compute rotation rates by assuming conservation of angular momentum around the center of mass at impact and that merged planetesimals relax to spherical shapes. At the end of our simulations, the majority of bodies that have undergone at least one merger are rotating faster than the breakup frequency. This implies that the assumption of completely inelastic collisions (perfect accretion), which is made in most simulations of planetary growth at sizes 1 km and above, is inappropriate. Our simulations reveal that, subsequent to the number of particles in the patch having been decreased by mergers to half its initial value, the presence of larger bodies in neighboring regions of the disk may limit the validity of simulations employing the patch approximation.     

Shock-wave heating model for chondrule formation: Hydrodynamic simulation of molten droplets exposed to gas flows

Millimeter-sized, spherical silicate grains abundant in chondritic meteorites, which are called as chondrules, are considered to be a strong evidence of the melting event of the dust particles in the protoplanetary disk. One of the most plausible scenarios is that the chondrule precursor dust particles are heated and melt in the high-velocity gas flow (shock-wave heating model). We developed the non-linear, time-dependent, and three-dimensional hydrodynamic simulation code for analyzing the dynamics of molten droplets exposed to the gas flow. We confirmed that our simulation results showed a good agreement in a linear regime with the linear solution analytically derived by Sekyia et al. [Sekyia, M., Uesugi, M., Nakamoto, T., 2003. Prog. Theor. Phys. 109, 717-728]. We found that the non-linear terms in the hydrodynamical equations neglected by Sekiya et al. [Sekiya, M., Uesugi, M., Nakamoto, T., 2003. Prog. Theor. Phys. 109, 717-728] can cause the cavitation by producing negative pressure in the droplets. We discussed that the fragmentation through the cavitation is a new mechanism to determine the upper limit of chondrule sizes. We also succeeded to reproduce the fragmentation of droplets when the gas ram pressure is stronger than the effect of the surface tension. Finally, we compared the deformation of droplets in the shock-wave heating with the measured data of chondrules and suggested the importance of other effects to deform droplets, for example, the rotation of droplets. We believe that our new code is a very powerful tool to investigate the hydrodynamics of molten droplets in the framework of the shock-wave heating model and has many potentials to be applied to various problems.     

Origin of three-dimensional shapes of chondrules: I. Hydrodynamics simulations of rotating droplet exposed to high-velocity rarefied gas flow

The origin of three-dimensional shapes of chondrules is an important information to identify their formation mechanism in the early solar nebula. The measurement of their shapes by using X-ray computed topography suggested that they are usually close to perfect spheres, however, some of them have rugby-ball-like (prolate) shapes [Tsuchiyama, A., Shigeyoshi, R., Kawabata, T., Nakano, T., Uesugi, K., Shirono, S., 2003. Lunar Planet. Sci. 34, 1271-1272]. We considered that the prolate shapes reflect the deformations of chondrule precursor dust particles when they are heated and melted in the high velocity gas flow. In order to reveal the origin of chondrule shapes, we carried out the three-dimensional hydrodynamics simulations of a rotating molten chondrule exposed to the gas flow in the framework of the shock-wave heating model for chondrule formation. We adopted the gas ram pressure acting on the chondrule surface of in a typical shock wave. Considering that the chondrule precursor dust particle has an irregular shape before melting, the ram pressure causes a net torque to rotate the particle. The estimated angular velocity is for the precursor radius of r₀=1 mm, though it has a different value depending on the irregularity of the shape. In addition, the rotation axis is likely to be perpendicular to the direction of the gas flow. Our calculations showed that the rotating molten chondrule elongates along the rotation axis, in contrast, shrinks perpendicularly to it. It is a prolate shape. The reason why the molten chondrule is deformed to a prolate shape was clearly discussed. Our study gives a complementary constraint for chondrule formation mechanisms, comparing with conventional chemical analyses and dynamic crystallization experiments that have mainly constrained the thermal evolutions of chondrules.     

Fragment-collision model for compound chondrule formation: Estimation of collision probability

We propose a new scenario for compound chondrule formation named as “fragment-collision model,” in the framework of the shock-wave heating model. A molten cm-sized dust particle (parent) is disrupted in the high-velocity gas flow. The extracted fragments (ejectors) are scattered behind the parent and the mutual collisions between them will occur. We modeled the disruption event by analytic considerations in order to estimate the probability of the mutual collisions assuming that all ejectors have the same radius. In the typical case, the molten thin () layer of the parent surface will be stripped by the gas flow. The stripped layer is divided into about 200 molten ejectors (assuming that the radius of ejectors is 300 μm) and then they are blown away by the gas flow in a short period of time (). The stripped layer is leaving from the parent with the velocity of depending on the viscosity, and we assumed that the extracted ejectors have a random velocity Δv of the same order of magnitude. Using above values, we can estimate the number density of ejectors behind the parent as . These ejectors occupy ∼9% of the space behind the parent in volume. Considering that the collision rate (number of collisions per unit time experienced by an ejector) is given by Rcoll=σcollneΔv, where σcoll is the cross-section of collision [e.g., Gooding, J.K., Keil, K., 1981. Meteoritics 16, 17-43], we obtain by substituting above values. Since most collisions occur within the short duration () before the ejectors are blown away, we obtain the collision probability of Pcoll∼0.36, which is the probability of collisions experienced by an ejector in one disruption event. The estimated collision probability is about one order of magnitude larger than the observed fraction of compound chondrules. In addition, the model predictions are qualitatively consistent with other observational data (oxygen isotopic composition, textural types, and size ratios of constituents). Based on these results, we concluded that this new model can be one of the strongest candidates for the compound chondrule formation. It should be noted that all collisions do not necessarily lead to the compound chondrule formation. The formation efficiency and the future works which should be investigated in the forthcoming paper are also discussed.     

Breakup of liquids by high velocity flow and size distribution of chondrules

As a new approach to understanding the chondrule formation process, we carried out aerodynamic experiments in which a liquid layer was attached to solid cores, and the breakup of this layer occurred by means of the interaction with a high-velocity gas flow. The size distribution of the dispersed droplets was investigated and compared with the size distributions of chondrules. Both distributions had an exponential form. Using the experimental results, the hydrodynamic pressure to produce the chondrule size distributions was estimated to be ∼ 10⁴ Pa.     

Sedna and the Oort Cloud around a migrating Sun

Recent numerical simulations have demonstrated that the Sun’s dynamical history within the Milky Way may be much more complex than that suggested by its current low peculiar velocity (Sellwood, J.A., Binney, J.J. [2002]. Mon. Not. R. Astron. Soc. 336, 785-796; Roškar, R., Debattista, V.P., Quinn, T.R., Stinson, G.S., Wadsley, J. [2008]. Astrophys. J. 684, L79-L82). In particular, the Sun may have radially migrated through the galactic disk by up to 5-6 kpc (Roškar, R., Debattista, V.P., Quinn, T.R., Stinson, G.S., Wadsley, J. [2008]. Astrophys. J. 684, L79-L82). This has important ramifications for the structure of the Oort Cloud, as it means that the Solar System may have experienced tidal and stellar perturbations that were significantly different from its current local galactic environment. To characterize the effects of solar migration within the Milky Way, we use direct numerical simulations to model the formation of an Oort Cloud around stars that end up on solar-type orbits in a galactic-scale simulation of a Milky Way-like disk formation. Surprisingly, our simulations indicate that Sedna’s orbit may belong to the classical Oort Cloud. Contrary to previous understanding, we show that field star encounters play a pivotal role in setting the Oort Cloud’s extreme inner edge, and due to their stochastic nature this inner edge sometimes extends to Sedna’s orbit. The Sun’s galactic migration heightens the chance of powerful stellar passages, and Sedna production occurs around ∼20-30% of the solar-like stars we study. Considering the entire Oort Cloud, we find its median distance depends on the minimum galactocentric distance attained during the Sun’s orbital history. The inner edge also shows a similar dependence but with increased scatter due to the effects of powerful stellar encounters. Both of these Oort Cloud parameters can vary by an order of magnitude and are usually overestimated by an Oort Cloud formation model that assumes a fixed galactic environment. In addition, the amount of material trapped in outer Oort Cloud orbits (a > 20,000 AU) can be extremely low and may present difficulties for traditional models of Oort Cloud formation and long-period comet production.