Generation of a dynamo magnetic field in a protoplanetary accretion disk

A new computational technique is developed that allows realistic calculations of dynamo magnetic field generation in disk geometries corresponding to protoplanetary and protostellar accretion disks. The approach is of sufficient generality to allow, in the future, a wide class of accretion disk problems to be solved. Here, basic modes of a disk dynamo are calculated. Spatially localized oscillatory states are found to occur in Keplerain disks. A physical interpretation is given that argues that spatially localized fields of the type found in these calculations constitute the basic modes of a Keplerian disk dynamo

What Fraction of Sun-like Stars have Planets?

The radial velocities of ~1800 nearby Sun-like stars are currently being monitored by eight high-sensitivity Doppler exoplanet surveys. Approximately 90 of these stars have been found to host exoplanets massive enough to be detectable. Thus at least ~5% of target stars possess planets. If we limit our analysis to target stars that have been monitored the longest (~15 years), ~11% possess planets. If we limit our analysis to stars monitored the longest and whose low surface activity allow the most precise velocity measurements, ~25% possess planets. By identifying trends of the exoplanet mass and period distributions in a sub-sample of exoplanets less-biased by selection effects, and linearly extrapolating these trends into regions of parameter space that have not yet been completely sampled, we find at least ~9% of Sun-like stars have planets in the mass and orbital period ranges Msin(i) > 0.3 M_Jupiter and P 0.1 M_Jupiter and P < 60 years. Even this larger area of the mass-period plane is less than 20% of the area occupied by our planetary system, suggesting that this estimate is still a lower limit to the true fraction of Sun-like stars with planets, which may be as large as ~100%.Comment: Conforms to version accepted by ApJ. Color version and movie available at http://bat.phys.unsw.edu.au/~charley/download/whatfrac

Planetary migration in evolving planetesimals discs

In the current paper, we further improved the model for the migration of planets introduced in Del Popolo et al. (2001) and extended to time-dependent planetesimal accretion disks in Del Popolo and Eksi (2002). In the current study, the assumption of Del Popolo and Eksi (2002), that the surface density in planetesimals is proportional to that of gas, is released. In order to obtain the evolution of planetesimal density, we use a method developed in Stepinski and Valageas (1997) which is able to simultaneously follow the evolution of gas and solid particles for up to 10^7 yrs. Then, the disk model is coupled to migration model introduced in Del Popolo et al. (2001) in order to obtain the migration rate of the planet in the planetesimal. We find that the properties of solids known to exist in protoplanetary systems, together with reasonable density profiles for the disk, lead to a characteristic radius in the range 0.03-0.2 AU for the final semi-major axis of the giant planet.Comment: IJMP A in prin

Protoplanetary Disk Turbulence Driven by the Streaming Instability: Non-Linear Saturation and Particle Concentration

We present simulations of the non-linear evolution of streaming instabilities in protoplanetary disks. The two components of the disk, gas treated with grid hydrodynamics and solids treated as superparticles, are mutually coupled by drag forces. We find that the initially laminar equilibrium flow spontaneously develops into turbulence in our unstratified local model. Marginally coupled solids (that couple to the gas on a Keplerian time-scale) trigger an upward cascade to large particle clumps with peak overdensities above 100. The clumps evolve dynamically by losing material downstream to the radial drift flow while receiving recycled material from upstream. Smaller, more tightly coupled solids produce weaker turbulence with more transient overdensities on smaller length scales. The net inward radial drift is decreased for marginally coupled particles, whereas the tightly coupled particles migrate faster in the saturated turbulent state. The turbulent diffusion of solid particles, measured by their random walk, depends strongly on their stopping time and on the solids-to-gas ratio of the background state, but diffusion is generally modest, particularly for tightly coupled solids. Angular momentum transport is too weak and of the wrong sign to influence stellar accretion. Self-gravity and collisions will be needed to determine the relevance of particle overdensities for planetesimal formation.Comment: Accepted for publication in ApJ (17 pages). Movies of the simulations can be downloaded at http://www.mpia.de/~johansen/research_en.ph

Toward a Deterministic Model of Planetary Formation IV: Effects of Type-I Migration

In a further development of a deterministic planet-formation model (Ida & Lin 2004), we consider the effect of type-I migration of protoplanetary embryos due to their tidal interaction with their nascent disks. During the early embedded phase of protostellar disks, although embryos rapidly emerge in regions interior to the ice line, uninhibited type-I migration leads to their efficient self-clearing. But, embryos continue to form from residual planetesimals at increasingly large radii, repeatedly migrate inward, and provide a main channel of heavy element accretion onto their host stars. During the advanced stages of disk evolution (a few Myr), the gas surface density declines to values comparable to or smaller than that of the minimum mass nebula model and type-I migration is no longer an effective disruption mechanism for mars-mass embryos. Over wide ranges of initial disk surface densities and type-I migration efficiency, the surviving population of embryos interior to the ice line has a total mass several times that of the Earth. With this reservoir, there is an adequate inventory of residual embryos to subsequently assemble into rocky planets similar to those around the Sun. But, the onset of efficient gas accretion requires the emergence and retention of cores, more massive than a few M_earth, prior to the severe depletion of the disk gas. The formation probability of gas giant planets and hence the predicted mass and semimajor axis distributions of extrasolar gas giants are sensitively determined by the strength of type-I migration. We suggest that the observed fraction of solar-type stars with gas giant planets can be reproduced only if the actual type-I migration time scale is an order of magnitude longer than that deduced from linear theories.Comment: 32 pages, 8 figures, 1 table, accepted for publication in Ap

Computer-generated global map of valley networks on Mars

The presence of valley networks (VN) on Mars suggests that early Mars was warmer and wetter than present. However, detailed geomorphic analyses of individual networks have not led to a consensus regarding their origin. An additional line of evidence can be provided by the global pattern of dissection on Mars, but the currently available global map of VN, compiled from Viking images, is incomplete and outdated. We created an updated map of VN by using a computer algorithm that parses topographic data and recognizes valleys by their morphologic signature. This computer-generated map was visually inspected and edited to produce the final updated map of VN. The new map shows an increase in total VN length by a factor of 2.3. A global map of dissection density, D, derived from the new VN map, shows that the most highly dissected region forms a belt located between the equator and mid-southern latitudes. The most prominent regions of high values of D are the northern Terra Cimmeria and the Margaritifer Terra where D reaches the value of 0.12 km−1 over extended areas. The average value of D is 0.062 km−1, only 2.6 times lower than the terrestrial value of D as measured in the same fashion. These relatively high values of dissection density over extensive regions of the planet point toward precipitation-fed runoff erosion as the primary mechanism of valley formation. Assuming a warm and wet early Mars, peculiarity of the global pattern of dissection is interpreted in the terms of climate controlling factors influenced by the topographic dichotomy

The Evolution of Protoplanetary Disks Around Millisecond Pulsars: The PSR 1257 +12 System

We model the evolution of protoplanetary disks surrounding millisecond pulsars, using PSR 1257+12 as a test case. Initial conditions were chosen to correspond to initial angular momenta expected for supernova-fallback disks and disks formed from the tidal disruption of a companion star. Models were run under two models for the viscous evolution of disks: fully viscous and layered accretion disk models. Supernova-fallback disks result in a distribution of solids confined to within 1-2 AU and produce the requisite material to form the three known planets surrounding PSR 1257+12. Tidal disruption disks tend to slightly underproduce solids interior to 1 AU, required for forming the pulsar planets, while overproducing the amount of solids where no body, lunar mass or greater, exists. Disks evolving under 'layered' accretion spread somewhat less and deposit a higher column density of solids into the disk. In all cases, circumpulsar gas dissipates on $\lesssim 10^{5}$ year timescales, making formation of gas giant planets highly unlikely.Comment: 16 pages, 17 figures, Accepted for publication in The Astrophysical Journal (September 20, 2007 issue

Growth and migration of solids in evolving protostellar disks I: Methods and Analytical tests

This series of papers investigates the early stages of planet formation by modeling the evolution of the gas and solid content of protostellar disks from the early T Tauri phase until complete dispersal of the gas. In this first paper, I present a new set of simplified equations modeling the growth and migration of various species of grains in a gaseous protostellar disk evolving as a result of the combined effects of viscous accretion and photo-evaporation from the central star. Using the assumption that the grain size distribution function always maintains a power-law structure approximating the average outcome of the exact coagulation/shattering equation, the model focuses on the calculation of the growth rate of the largest grains only. The coupled evolution equations for the maximum grain size, the surface density of the gas and the surface density of solids are then presented and solved self-consistently using a standard 1+1 dimensional formalism. I show that the global evolution of solids is controlled by a leaky reservoir of small grains at large radii, and propose an empirically derived evolution equation for the total mass of solids, which can be used to estimate the total heavy element retention efficiency in the planet formation paradigm. Consistency with observation of the total mass of solids in the Minimum Solar Nebula augmented with the mass of the Oort cloud sets strong upper limit on the initial grain size distribution, as well as on the turbulent parameter \alphat. Detailed comparisons with SED observations are presented in a following paper.Comment: Submitted to ApJ. 23 pages and 13 figure

Dust Size Growth and Settling in a Protoplanetary Disk

We have studied dust evolution in a quiescent or turbulent protoplanetary disk by numerically solving coagulation equation for settling dust particles, using the minimum mass solar nebular model. As a result, if we assume an ideally quiescent disk, the dust particles settle toward the disk midplane to form a gravitationally unstable layer within 2x10^3 - 4x10^4 yr at 1 - 30 AU, which is in good agreement with an analytic calculation by Nakagawa, Sekiya, & Hayashi (1986) although they did not take into account the particle size distribution explicitly. In an opposite extreme case of a globally turbulent disk, on the other hand, the dust particles fluctuate owing to turbulent motion of the gas and most particles become large enough to move inward very rapidly within 70 - 3x10^4 yr at 1 - 30 AU, depending on the strength of turbulence. Our result suggests that global turbulent motion should cease for the planetesimal formation in protoplanetary disks.Comment: 27 pages, 8 figures, accepted for publication in the Ap