Detailed Calculations of the Efficiency of Planetesimal Accretion in the Core-accretion Model
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
We present calculations of the early stages of the formation of Jupiter via core nucleated accretion and gas capture. The core begins as a seed body of about 350 kilometers in radius and orbits in a swarm of planetesimals whose initial radii range from 15 meters to 100 kilometers. We follow the evolution of the swarm by accounting for growth and fragmentation, viscous and gravitational stirring, and for drag-induced migration and velocity damping. Gas capture by the core substantially enhances the cross-section of the planet for accretion of small planetesimals. The dust opacity within the atmosphere surrounding the planetary core is computed self-consistently, accounting for coagulation and sedimentation of dust particles released in the envelope as passing planetesimals are ablated. The calculation is carried out at an orbital semi-major axis of 5.2 AU and an initial solids' surface density of 10/g/cm^2 at that distance. The results give a core mass of 7 Earth masses and an envelope mass of approximately 0.1 Earth mass after 500,000 years, at which point the envelope growth rate surpasses that of the core. The same calculation without the envelope gives a core mass of only 4 Earth masses.
In the core accretion model of giant planet formation, when the core reaches critical mass, hydrostatic equilibrium is no longer possible and gas accretion ensues. If the envelope is radiative, the critical core mass is nearly independent of the boundary conditions and is roughly M(sub crit) ~ 10Mass of the Earth (with weak dependence on the rate of planetesimal accretion M(sub core) and the disk opacity k). Given that such a core may form at the present location of Jupiter in a time comparable to its Type I migration time (10(exp 5) - 10(exp 6) years) provided that the nebula was significantly enhanced in solids with respect to the MMSN and stall at this location in a weakly turbulent (alpha approximately less than 10(exp -4) disk, it may be appropriate to assume that such objects inevitably form and drive the evolution of late-phase T Tauri star disks. Here we investigate the final masses of giant planets in disks with one or more than one such cores. Although the presence of several planets would lead to Type II migration (due to the effective viscosity resulting from the planetary tidal torques), we ignore this complication for now and simply assume that each core has stalled at its location in the disk. Once a core has achieved critical mass, its gaseous accretion is governed by the given Kelvin-Helmholtz timescale.
One of the justifications motivating efforts to detect and characterize young extrasolar giant planets has been to measure atmospheric composition for comparison with that of the primary star. If the enhancement of heavy elements in the atmospheres of extrasolar giant planets, like it is for their solar system analogs, is inversely proportional to mass, then it is likely that these worlds formed by core accretion. However in practice it has been very difficult to constrain metallicity because of the complex effect of clouds. Cloud opacity varies both vertically and, in some cases, horizontally through the atmosphere. Particle size and composition, both of which impact opacity, are difficult challenges both for forward modeling and retrieval studies. In my presentation I will discuss systematic efforts to improve cloud studies to enable more reliable determinations of atmospheric composition. These efforts are relevant both to discerning composition of directly imaged young planets from ground based telescopes and future space based missions, such as WFIRST and LUVOIR.
This project investigates the origin of giant planets, both in the Solar System and around other stars. It is assumed that the planets form by the core accretion process: small solid particles in a disk surrounding a young star gradually coagulate into objects of a few kilometers in size, known as planetesimals, which then accumulate into solid protoplanetary cores. Once the cores have become large enough, they are able to attract gas from the surrounding disk to form the deep gaseous envelope of the giant planet. Our code simulates giant planet growth in a spherical approximation, and it has been quite successful in addressing a number of basic planetary properties. Further improvements to the code have been made to achieve a more realistic understanding of planetary formation. The computations of the models were based on an earlier version of our code and were stopped at the onset of runaway gas accretion. Now, improved boundary conditions have been incorporated into the code to allow for hydrodynamic inflow of gas and to handle the late stages of evolution when the planet evolves at constant mass. These changes were made to the version of the code that uses a constant accretion rate and to the version that uses a self-consistent method for calculating both the solid and gas accretion rates. The equation of state has been updated to incorporate the detailed tables of Saumon, Chabrier, and Van Horn. The opacities were updated to include the results of Alexander and Ferguson. The outer boundary conditions were modified. During the accretion phase when the planet's radius is between the accretion radius and the tidal radius, we set the outer boundary at a 'modified' accretion radius, which is the point where thermal energy is enough to bring gas to the edge of the Hill sphere.
A generalized model which accounts for the effects of initially uniform and slow rotation is defined for the spherical collapse of a singular isothermal sphere such as protosolar and binary nebulae. An initial unstable equilibrium state is described for a sound speed of 0.35 km/sec and a rotation rate of 10 to the -14th/sec for the molecular cloud surrounding the accreting core. The total angular momentum and mass of the inner cloud is set equal to solar system values. The evolution of the collapse is traced by applying a perturbation analysis to the similarity solution for a nonrotating condition, and matched asymptotic expansions solve the hydrodynamic equations. The model is concluded a valid tool for studying star and nebular disk formation.
Evolution of Mars' noble gases through two stages of hydrodynamic escape early in planetary history has been proposed previously by the author. In the first evolutionary stage of this earlier model, beginning at a solar age of approximately 50 m.y., fractionating escape of a H2-rich primordial atmosphere containing CO2, N2, and the noble gases in roughly the proportions found in primitive carbonaceous (CI) chondrites is driven by intense extreme-ultraviolet (EUV) leads to a long (approximately 80 m.y.) period of quiescence, followed by an abrupt degassing of remnant H2, CO2, and N2 from the mantle and of solar-composition noble gases lighter than Xe from the planet's volatile-rich accretional core. Degassed H refuels hydrodynamic loss in a waning but still potent solar EUV flux. Atmospheric Xe, Kr, and Ar remaining at the end of this second escape stage, approximately 4.2 G.y. ago, have evolved to their present-day abundances and compositions. Residual Ne continues to be modified by accretion of solar wind gases throughout the later history of the planet. This model does not address a number of processes that now appear germane to Martian atmospheric history. One, gas loss and fractionation by sputtering, has recently been shown to be relevant. Another, atmospheric erosion, appears increasingly important. In the absence then of a plausible mechanism, the model did not consider the possibility of isotopic evolution of noble gases heavier than Ne after the termination of hydrodynamic escape. Subsequent non-thermal loss of N was assumed, in an unspecified way, to account for the elevation of N from the model value of approximately 250 percent at the end of the second escape stage to approximately 620 percent today. Only qualitative attention was paid to the eroding effects of impact on abundances of all atmophilic species prior to the end of heavy bombardment approximately 3.8 G.y. ago. No attempt was made to include precipitation and recycling of carbonates in tracking the pressure and isotopic history of CO2.
The two X-ray transform telescopes on the Pinhole/Occulter Facility (P/OF) will provide sub-arcsecond angular resolution, high sensitivity imaging in hard X-rays for the first time, as well as improved capabilities for timing studies and spectroscopy. Although the design was driven by solar observing requirements, the result is an instrument of major importance for astronomical as well as solar observations. Fundamental astrophysical measurements could be made on supernova remnants, clusters of galaxies, active galaxy cores, accreting binary systems, and active coronae of nearby stars. Ways in which astronomical observations will benefit solar research are discussed.
There are lines of evidence which relate bounds on the primordial water content of the Earth's mantle to a magma ocean and the accompanying Earth accretion process. We assume initially (before a magma ocean could form) that as the Earth accreted, it grew from volatile- (H2O, CO2, NH3, CH4, SO2, plus noble) gas-rich planetesimals, which accreted to form an initial 'primitive accretion core' (PAC). The PAC retained the initial complement of planetesimal gaseous components. Shock wave experiments in which both solid, and more recently, the gaseous components of materials such as serpentine and the Murchison meteorite have demonstrated that planetesimal infall velocities of less than 0.5 km/sec, induce shock pressures of less than 0.5 GPa and result in virtually complete retention of planetary gases.
The evolution of interstellar molecular hydrogen was studied, with a special interest for the formation and evolution of molecular clouds and star formation within them, by a two-dimensional hydrodynamical simulation performed on a rectangular grid of physical sizes on the order of 100 pc. It is filled with an initial density of approx. 1 cm(exp -3), except for one cell (approx. 1 pc(exp 2)) at the center of the grid where an accretion core of 1-10(exp 3) solar masses is placed. The grid is co-moving with the gridcenter that is on a circular orbit around the Galactic center and that also is the guiding center of epicyclic approximation of orbits of the matter surrounding it. The initial radial velocity is zero; to account for differential rotation the initial tangential velocity (i.e. the movement around the galactic center) is proportional to the radial distance to the grid center. The rate is comparable to the rotation rate at the Local Standard of Rest. The influence of galactic rotation is noticed by spiral or elliptical forms, but on much longer time scales than self gravitation and cooling processes. Density and temperature are kept constant at the boundaries and no inflow is allowed along the tangential boundaries.
All types of coarse-grained components in CM chondrites are surrounded by fine-grained dust coatings, but the origin of these rims is not yet clear. Although a strictly nebular origin seems likely for rims in the relatively unaltered type 3 chondrites, the rims in CM chondrites are dominated by secondary alteration phases. It has been argued that either the coarse-grained cores accreted altered rim materials while still in the nebula or that alteration of primary rim phases occurred on the CM parent body. To constrain the origin of alteration phases in rim material, we have analyzed the textures and mineral associations from 10 CM chondritic falls by optical and scanning electron microscopy. Our results indicate that the secondary phases in CM chondritic rims were produced by parent body fluid-rock interactions which redefined some primary rim textures and may have produced, in some cases, both coarse-grained components and the rims that surround them. Textural features demonstrate the interactive exchange of alteration fluids between rims, matrix, and chondrules on the CM parent body. For example, most matrix-rim contacts are gradational, suggesting the synchronous alteration of both components. Several observations suggest the possibility of in situ rim production. For example, tochilinite and phyllosilicates commonly form rims around matrix carbonates, which are generally believed to have precipitated from alteration fluids on the CM parent body. This suggests that the rims surrounding matrix carbonates may also have been produced by alteration processes. Partially replaced chondrule olivines bear a striking resemblance to many rimmed olivines in the matrix which suggests, by analogy, that site-specific precipitation of S-bearing phases may also be responsible for the occurrence of many tochilinite-rich rims around isolated matrix olivines. Non-silicate rims precipitate around olivines of any composition, but the process is most effective for fayalitic olivines. Most of the remaining olivines in CM chondrites are relatively Mg-rich, which suggests that the precipitation of S-bearing rims on olivines may not have been an important process in the aqueous alteration of CM chondrites. We conclude that: (1) precursor rim materials in CM chondrites were subjected to pervasive aqueous alteration on the CM parent body; and (2) textures and mineral associations observed in CM chondrites also suggest the possibility of in situ rim production.
In order to create a coherent scenario of satellite formation. the source of the solids (rock-metal and ice) that will eventually make up the satellites must be considered. While it is customary to use a solar composition mixture with a gas/solid mass ratio of about 100, at the tail end of the formation of the giant planet (when satellite formation is thought to have taken place) the fraction of solids entrained in the gas (particles with sizes lower than the decoupling size about 1 m for typical nebula parameters) is likely to be significantly lower than cosmic. In particular, in the core accretion model of giant planet formation one expects low dust and rubble content at late times due to particle coagulation leading to a collisional distribution of particle sizes with most of the mas residing in objects 1 km or larger, which are not coupled to the gas and whose dynamics must be followed independently. As a result, flow of gas into circumplanetary orbits is not sufficient to constrain the mas available to form satellites.
The session"Origin of Planetary Systems" presented the following reports:Planetary Formation in the the Gamma Cephei System by Core:Accretion; Accretion and Heating of Particles by Supersonic Planetesimals; Planetesimal Accretion in Close Binary Systems; Collisions, Gas Flow, and the Formation of Planetesimals; Hydrodynamic Escape of a Proto-Atmosphere Just After a Giant Impact; X-Ray Flare Induced Shock Waves and Chondrule Formation in Upper Solar Nebula; and Making Water Worlds: The Role of 26Al.
IRAS l6253-2429. the source of the Wasp-Waist Nebula seen in Spitzer IRAC images, is an isolated very low luminosity ("VeLLO") Class 0 protostar in the nearby rho Ophiuchi cloud. We present VLA ammonia mapping observations of the dense gas envelope feeding the central core accreting system. We find a flattened envelope perpendicular to the outflow axis, and gas cavities that appear to cradle the outflow lobes as though carved out by the flow and associated (apparently precessing) jet. Based on the NH3 (1,1) and (2,2) emission distribution, we derive the mass, velocity fields and temperature distribution for the envelope. We discuss the combined evidence for this source as possibly one of the youngest and lowest mass sources in formation yet known.
We present new 1-1.25 micron (z and J band) Subaru/IRCS and 2 micron (K band) VLT/NaCo data for HR 8799 and a rereduction of the 3-5 micron MMT/Clio data first presented by Hinz et al. Our VLT/NaCo data yield a detection of a fourth planet at a projected separation of approximately 15 AU--"HR 8799e ." We also report new, albeit weak detections of HR 8799b at 1.03 micron and 3.3 micron. Empirical comparisons to field brown dwarfs show that at least HR 8799b and HR 8799c, and possibly HR 8799d, have near-to-mid-IR colors/ magnitudes significantly discrepant from the L/T dwarf sequence. Standard cloud deck atmosphere models appropriate for brown dwarfs provide only (marginally) statistically meaningful fits to HR 8799b and c for unphysically small radii. Models with thicker cloud layers not present in brown dwarfs reproduce the planets' spectral energy distributions far more accurately and without the need for resealing the planets' radii. Our preliminary modeling suggests that HR 8799b has log(g) = 4-4.5, T(sub eff) = 900 K. while HR 8799c, d, and (by inference) e have log(g) = 4-4.5, T(sub eff) = 1000-1200 K. Combining results from planet evolution models and new dynamical stability limits implies that the masses of HR 8799b, c, d, and e are 6-7 M(sub j), 7-10 M(sub j), 7-10 M(sub j), and 7-10 M(sub j). "Patchy" cloud prescriptions may provide even better fits to the data and may lower the estimated surface gravities and masses. Finally, contrary to some recent claims, forming the HR 8799 planets by core accretion is still plausible, although such systems are likely rare.
Based on data in the aeronautics literature, we have derived functions for the drag coefficients of spheres and cubes as a function of Mach number. Experiments have shown that spheres and cubes exhibit an abrupt factor-of-two decrease in the drag coefficient as the object slows through the transonic regime. Irregularly shaped objects such as meteorites likely exhibit a similar trend. These functions are implemented in an otherwise simple projectile motion model, which is applicable to the non-ablative dark flight of meteors (speeds less than .+3 km/s). We demonstrate how these functions may be used as upper and lower limits on the drag coefficient of meteors whose shape is unknown. A Mach-dependent drag coefficient is potentially important in other planetary and astrophysical situations, for instance, in the core accretion scenario for giant planet formation.