Accretion Disks and Eruptive Phenomena

This paper describes eruptive phenomena in pre-main sequence stars. The eruptions of FU Orionis stars have much in common with outbursts in other accreting systems, such as dwarf novae and some symbiotic stars. These common features are best understood as increases in the rate material flows through an accretion disk. The spectroscopic properties, decay of the light curves, and outflow phenomena suggest that these outbursts arise from thermal instabilities in the disk. Available data and estimates for recurrence times indicate that young stars can accrete much, perhaps all, of their mass in FU Ori accretion events. Future observations can test this notion and place better constraints on the importance of eruptive events in the early life of a low mass star.Comment: to appear in The Physics of Star Formation and Early Stellar Evolution, edited by C. J. Lada and N. Kylafis (30 pages, 11 figures) This version corrects several typos and a mistaken impression left in the derivation of the disk radial temperature profil

Terrestrial planet formation: Dynamical shake-up and the low mass of Mars

We consider a dynamical shake-up model to explain the low mass of Mars and the lack of planets in the asteroid belt. In our scenario, a secular resonance with Jupiter sweeps through the inner solar system as the solar nebula depletes, pitting resonant excitation against collisional damping in the Sun's protoplanetary disk. We report the outcome of extensive numerical calculations of planet formation from planetesimals in the terrestrial zone, with and without dynamical shake-up. If the Sun's gas disk within the terrestrial zone depletes in roughly a million years, then the sweeping resonance inhibits planet formation in the asteroid belt and substantially limits the size of Mars. This phenomenon likely occurs around other stars with long-period massive planets, suggesting that asteroid belt analogs are common.Comment: AJ, in press; 49 pages, 8 figure

Numerical Simulations of Gaseous Disks Generated from Collisional Cascades at the Roche Limits of White Dwarf Stars

We consider the long-term evolution of gaseous disks fed by the vaporization of small particles produced in a collisional cascade inside the Roche limit of a 0.6 Msun white dwarf. Adding solids with radius \r0\ at a constant rate $\dot{M}_0$ into a narrow annulus leads to two distinct types of evolution. When $\dot{M}_0 > \dot{M}_{0,crit}$ = $3 \times 10^4 ~ (r_0 / {\rm 1~km})^{3.92}$~g s$^{-1}$, the cascade generates a fairly steady accretion disk where the mass transfer rate of gas onto the white dwarf is roughly $\dot{M}_0$ and the mass in gas is $M_g \approx 2.3 \times 10^{22} ~ (\dot{M}_0 / 10^{10}~g~s^{-1}) ~ ({\rm 1500~K} / T_0) ~ (10^{-3} / \alpha)$~g, where $T_0$ is the temperature of the gas near the Roche limit and $\alpha$ is the dimensionless viscosity parameter. If $\dot{M}_0 < \dot{M}_{0,crit}$, the system alternates between high states with large mass transfer rates and low states with negligible accretion. Although either mode of evolution adds significant amounts of metals to the white dwarf photosphere, none of our calculations yield a vertically thin ensemble of solids inside the Roche limit. X-ray observations can place limits on the mass transfer rate and test this model for metallic line white dwarfs.Comment: 30 pages and 8 figures, ApJ, accepte

Coagulation Calculations of Icy Planet Formation Around 0.1--0.5~\msun\ Stars: Super-Earths From Large Planetestimals

We investigate formation mechanisms for icy super-Earth mass planets orbiting at 2-20 AU around 0.1-0.5 solar mass stars. A large ensemble of coagulation calculations demonstrates a new formation channel: disks composed of large planetesimals with radii of 30-300 km form super-Earths on time scales of roughly 1 Gyr. In other gas-poor disks, a collisional cascade grinds planetesimals to dust before the largest planets reach super-Earth masses. Once icy Earth-mass planets form, they migrate through the leftover swarm of planetesimals at rates of 0.01-1 AU per Myr. On time scales of 10 Myr to 1 Gyr, many of these planets migrate through the disk of leftover planetesimals from semimajor axes of 5-10 AU to 1-2 AU. A few per cent of super-Earths might migrate to semimajor axes of 0.1-0.2 AU. When the disk has an initial mass comparable with the minimum mass solar nebula scaled to the mass of the central star, the predicted frequency of super-Earths matches the observed frequency.Comment: 32 pages, 16 figures, ApJ accepte

Making Planet Nine: Pebble Accretion at 250--750 AU in a Gravitationally Unstable Ring

We investigate the formation of icy super-Earth mass planets within a gravitationally unstable ring of solids orbiting at 250-750 AU around a 1 solar mass star. Coagulation calculations demonstrate that a system of a few large oligarchs and a swarm of pebbles generates a super-Earth within 100-200 Myr at 250 AU and within 1-2 Gyr at 750 AU. Systems with more than ten oligarchs fail to yield super-Earths over the age of the solar system. As these systems evolve, destructive collisions produce detectable debris disks with luminosities of $10^{-5}$ to $10^{-3}$ relative to the central star.Comment: 24 pages of text, 1 table, 8 figures, ApJ submitted, comments welcom

Numerical Simulations of Collisional Cascades at the Roche Limits of White Dwarf Stars

We consider the long-term collisional and dynamical evolution of solid material orbiting in a narrow annulus near the Roche limit of a white dwarf. With orbital velocities of 300 km/sec, systems of solids with initial eccentricity $e \gtrsim 10^{-3}$ generate a collisional cascade where objects with radii $r \lesssim$ 100--300 km are ground to dust. This process converts 1-100 km asteroids into 1 $\mu$m particles in $10^2 - 10^6$ yr. Throughout this evolution, the swarm maintains an initially large vertical scale height $H$. Adding solids at a rate $\dot{M}$ enables the system to find an equilibrium where the mass in solids is roughly constant. This equilibrium depends on $\dot{M}$ and $r_0$, the radius of the largest solid added to the swarm. When $r_0 \lesssim$ 10 km, this equilibrium is stable. For larger $r_0$, the mass oscillates between high and low states; the fraction of time spent in high states ranges from 100% for large $\dot{M}$ to much less than 1% for small $\dot{M}$. During high states, the stellar luminosity reprocessed by the solids is comparable to the excess infrared emission observed in many metallic line white dwarfs.Comment: 37 pages of text, 12 figures, ApJ, accepte

Rapid Formation of Icy Super-Earths and the Cores of Gas Giant Planets

We describe a coagulation model that leads to the rapid formation of super-Earths and the cores of gas giant planets. Interaction of collision fragments with the gaseous disk is the crucial element of this model. The gas entrains small collision fragments, which rapidly settle to the disk midplane. Protoplanets accrete the fragments and grow to masses of at least 1 Earth mass in roughly 1 Myr. Our model explains the mass distribution of planets in the Solar System and predicts that super-Earths form more frequently than gas giants in low mass disks.Comment: ApJLetters, accepted; 10 pages of text and 2 figure

Variations on Debris Disks IV. An Improved Analytical Model for Collisional Cascades

We derive a new analytical model for the evolution of a collisional cascade in a thin annulus around a single central star. In this model, $r_{max}~$ the size of the largest object declines with time (t); $r_{max} \propto t^{-\gamma}$, with $\gamma$ = 0.1-0.2. Compared to standard models where $r_{max}~$ is constant in time, this evolution results in a more rapid decline of $M_d$, the total mass of solids in the annulus and $L_d$, the luminosity of small particles in the annulus: $M_d \propto t^{-(\gamma + 1)}~$ and $L_d \propto t^{-(\gamma/2 + 1)}~~$. We demonstrate that the analytical model provides an excellent match to a comprehensive suite of numerical coagulation simulations for annuli at 1 AU and at 25 AU. If the evolution of real debris disks follows the predictions of the analytical or numerical models, the observed luminosities for evolved stars require up to a factor of two more mass than predicted by previous analytical models.Comment: ApJ, in press, 22 pages of text and 14 figure

Collisional Cascade Caclulations for Irregular Satellite Swarms in Fomalhaut b

We describe an extensive suite of numerical calculations for the collisional evolution of irregular satellite swarms around 1--300 M-earth planets orbiting at 120 AU in the Fomalhaut system. For 10--100 M-earth planets, swarms with initial masses of roughly 1% of the planet mass have cross-sectional areas comparable to the observed cross-sectional area of Fomalhaut b. Among 30--300 M-earth planets, our calculations yield optically thick swarms of satellites for ages of 1-10 Myr. Observations with HST and ground-based AO instruments can constrain the frequency of these systems around stars in the beta Pic moving group and possibly other nearby associations of young stars.Comment: 46 pages, 30 figures, ApJ, accepte

Migration of small moons in Saturn's rings

The motions of small moons through Saturn's rings provide excellent tests of radial migration models. In theory, torque exchange between these moons and ring particles leads to radial drift. We predict that moons with Hill radii r_H ~ 2-24 km should migrate through the A ring in 1000 yr. In this size range, moons orbiting in an empty gap or in a full ring eventually migrate at the same rate. Smaller moons or moonlets -- such as the propellers (e.g., Tiscareno et al. 2006) -- are trapped by diffusion of disk material into corotating orbits, creating inertial drag. Larger moons -- such as Pan or Atlas -- do not migrate because of their own inertia. Fast migration of 2-24 km moons should eliminate intermediate-size bodies from the A ring and may be responsible for the observed large-radius cutoff of r_H ~ 1-2 km in the size distribution of the A ring's propeller moonlets. Although the presence of Daphnis (r_H ~ 5 km) inside the Keeler gap challenges this scenario, numerical simulations demonstrate that orbital resonances and stirring by distant, larger moons (e.g., Mimas) may be important factors. For Daphnis, stirring by distant moons seems the most promising mechanism to halt fast migration. Alternatively, Daphnis may be a recent addition to the ring that is settling into a low inclination orbit in ~10^3 yr prior to a phase of rapid migration. We provide predictions of observational constraints required to discriminate among possible scenarios for Daphnis.Comment: ApJ, accepted; 47 pages, 14 figure