We consider a scenario where the small satellites of Pluto and Charon grew within a disk of debris from an impact between Charon and a trans-Neptunian Object (TNO). After Charon's orbital motion boosts the debris into a disk-like structure, rapid orbital damping of meter-size or smaller objects is essential to prevent the subsequent re-accretion or dynamical ejection by the binary. From analytical estimates and simulations of disk evolution, we estimate an impactor radius of 30-100 km; smaller (larger) radii apply to an oblique (direct) impact. Although collisions between large TNOs and Charon are unlikely today, they were relatively common within the first 0.1-1 Gyr of the solar system. Compared to models where the small satellites agglomerate in the debris left over by the giant impact that produced the Pluto-Charon binary planet, satellite formation from a later impact on Charon avoids the destabilizing resonances that sweep past the satellites during the early orbital expansion of the binary.

Current theory considers two options for the formation of the Pluto-Charon binary (Canup 2005, 2011; Desch 2015). In the `hit-and-run' model, a lower mass projectile barely hits the more massive Pluto, kicks up some debris, and remains bound to Pluto (see also Asphaug et al. 2006). In a `graze-and-merge' scenario, the projectile ejects substantial debris as it merges with Pluto (see also Canup 2001). To investigate the graze-and-merge idea in more detail, we consider the growth of Charon-mass objects within a circum-Pluto ring of solids. Numerical calculations demonstrate that Charon analogs form rapidly within a swarm of planetesimals with initial radii of 145-230 km. On time scales of roughly 30-100 days, newly-formed Charon analogs have semimajor axes, a = 5-6 Pluto radii, and orbital eccentricities, e = 0.1-0.3, similar to Charon analogs that remain bound after hit-and-run collisions with Pluto. Although the early growth of Charon analogs generates rings of small particles at a = 50-275 Pluto radii, ejection of several 145-230 km leftovers by the central Pluto-Charon binary removes these small solids in 10-100 yr. Simple estimates suggest small particles might survive the passage of 10-20 km objects ejected by the central binary. Our results indicate that the Pluto-Charon circumbinary satellite system was not formed by a graze-and-merge impact when the formation of Charon within a circum-Pluto disk leads to the ejection of several 100-200 km particles through the orbital plane of the Pluto-Charon binary. If a growing Charon ejects only much smaller particles, however, graze-and-merge impacts are a plausible formation channel for the Pluto-Charon binary and an ensemble of small, circumbinary satellites.

We consider the impact of electromagnetic induction and Ohmic heating on a conducting planetary object that orbits a magnetic star. Power dissipated as heat saps orbital energy. If this heat is trapped by an insulating crust or mantle, interior temperatures increase substantially. We provide a quantitative description of this behavior and discuss the astrophysical scenarios in which it might occur. Magnetic fields around some main-sequence stars and white dwarfs are strong enough to cause the decay of close-in orbits of asteroids and dwarf planets, drawing them through the Roche limit on Myr time scales. We confirm that Ohmic heating around neutron stars is driven by the rotation of the stellar magnetic dipole, not orbital dynamics. In any case, heating can raise interior temperatures of asteroids or dwarf planets on close-in orbits to well above liquidus. Hot material escaping to the surface may lead to volcanic ejections that can obscure the host star (as in the light curve of KIC 8462852) and pollute its atmosphere (as observed with metal-rich white dwarfs). We speculate that mixing of a volatile-rich mantle or crust with material from an induction-heated core may lead to an explosion that could destroy the asteroid prior to tidal break-up.

During 2005-2012, images from Hubble Space Telescope (HST) revealed four moons orbiting Pluto-Charon (Weaver et al 2006, Showalter et al 2011, 2012). Although their orbits and geometric shapes are well-known, the 2$\sigma$ uncertainties in the masses of the two largest satellites - Nix and Hydra - are comparable to their HST masses (Brozovic et al 2015, Showalter & Hamilton 2015, Weaver et al 2016). Remarkably, gravitational $n$-body computer calculations of the long-term system stability on 0.1-1 Gyr time scales place much tighter constraints on the masses of Nix and Hydra, with upper limits $\sim$ 10% larger than the HST mass. Constraints on the mass density using size measurements from New Horizons suggest Nix and Hydra formed in icier material than Pluto and Charon.

We investigate the nature of nearby (10-15 kpc) high-speed stars in the Gaia DR2 archive identified on the basis of parallax, proper motion and radial velocity. Together with a consideration of their kinematic, orbital, and photometric properties, we develop a novel strategy for evaluating whether high speed stars are statistical outliers of the bound population or unbound stars capable of escaping the Galaxy. Out of roughly 1.5 million stars with radial velocities, proper motions, and 5-sigma parallaxes, we identify just over 100 high-speed stars. Of these, only two have a nearly 100% chance of being unbound, with indication that they are not just bound outliers; both are likely hyper-runaway stars. The rest of the high speed stars are likely statistical outliers. We use the sample of high-speed stars to demonstrate that radial velocity alone provides a poor discriminant of nearby, unbound stars. However, nearby, unbound stars are efficiently identified from the tangential velocity, using just parallax and proper motion. Within the full Gaia DR2 archive of stars with 5-sigma parallax and proper motion but no radial velocity, we identify a sample of 19 with speeds significantly larger than the local escape speed of the Milky Way based on tangential motion alone.

A stellar flare can brighten a planet in orbit around its host star, producing a light curve with a faint echo. This echo, and others from subsequent flares, can lead to the planet's discovery, revealing its orbital configuration and physical characteristics. A challenge is that an echo is faint relative to the flare and measurement noise. Here, we use a method, based on autocorrelation function estimation, to extract faint planetary echoes from stellar flare light curves. A key component of our approach is that we compensate for planetary motion, measures of echo strength are then co-added into a strong signal. Using simple flare models in simulations, we explore the feasibility of this method with current technology for detecting planets around nearby M dwarfs. We also illustrate how our method can tightly constrain a planet's orbital elements and the mass of its host star. This technique is most sensitive to giant planets within 0.1 au of active flare stars and offers new opportunities for planet discovery in orientations and configurations that are inaccessible with other planet search methods.

Using a large suite of n-body simulations, we explore the discovery space for new satellites in the Pluto-Charon system. For the adopted masses and orbits of the known satellites, there are few stable prograde or polar orbits with semimajor axes $a \lesssim 1.1~a_H$, where $a_H$ is the semimajor axis of the outermost moon Hydra. Small moons with radii $r \lesssim$ 2 km and $a \lesssim 1.1~a_H$ are ejected on time scales ranging from several yr to more than 10 Myr. Orbits with $a \gtrsim 1.1~a_H$ are stable on time scales exceeding 100 Myr. Near-IR and mid-IR imaging with JWST and ground-based occultation campaigns with 2-3-m class telescopes can detect 1-2 km satellites outside the orbit of Hydra. Searches for these moons enable new constraints on the masses of the known satellites and on theories for circumbinary satellite formation.

We consider how the gravity of the Galactic disk and the Large Magellanic Cloud (LMC) modifies the radial motions of hypervelocity stars (HVSs) ejected from the Galactic Center. For typical HVSs ejected towards low (high) Galactic latitudes, the disk bends trajectories by up to 30 degrees (3-10 deg). For many lines-of-sight through the Galaxy, the LMC produces similar and sometimes larger deflections. Bound HVSs suffer larger deflections than unbound HVSs. Gravitational focusing by the LMC also generates a factor of two overdensity along the line-of-sight towards the LMC. With large enough samples, observations can detect the non-radial orbits and the overdensity of HVSs towards the LMC. For any Galactic potential model, the Galactic rest-frame tangential velocity provides an excellent way to detect unbound and nearly bound HVSs within 10 kpc of the Sun. Similarly, the rest-frame radial velocity isolates unbound HVSs beyond 10-15 kpc from the Sun. Among samples of unbound HVSs, measurements of the radial and tangential velocity serve to distinguish Galactic Center ejections from other types of high velocity stars.

The recent 'Planet Nine' hypothesis has led to many observational and archival searches for this giant planet proposed to orbit the Sun at hundreds of astronomical units. While trans-Neptunian object searches are typically conducted in the optical, models suggest Planet Nine could be self-luminous and potentially bright enough at ~3-5 microns to be detected by the Wide-field Infrared Survey Explorer (WISE). We have previously demonstrated a Planet Nine search methodology based on time-resolved WISE coadds, allowing us to detect moving objects much fainter than would be possible using single-frame extractions. In the present work, we extend our 3.4 micron (W1) search to cover more than three quarters of the sky and incorporate four years of WISE observations spanning a seven year time period. This represents the deepest and widest-area WISE search for Planet Nine to date. We characterize the spatial variation of our survey's sensitivity and rule out the presence of Planet Nine in the parameter space searched at W1 < 16.7 in high Galactic latitude regions (90% completeness).

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.

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.

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.

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.

A distant, as yet unseen ninth planet has been invoked to explain various observations of the outer solar system. While such a 'Planet Nine', if it exists, is most likely to be discovered via reflected light in the optical, it may emit much more strongly at 3$-$5$\mu$m than simple blackbody predictions would suggest, depending on its atmospheric properties (Fortney et al. 2016). As a result, Planet Nine may be detectable at 3.4$\mu$m with WISE, but single exposures are too shallow except at relatively small distances ($d_9 \lesssim 430$ AU). We develop a method to search for Planet Nine far beyond the W1 single-exposure sensitivity, to distances as large as 800 AU, using inertial coadds of W1 exposures binned into $\sim$1 day intervals. We apply our methodology to $\sim$2000 square degrees of sky identified by Holman & Payne (2016) as a potentially likely Planet Nine location, based on the Fienga et al. (2016) Cassini ranging analysis. We do not detect a plausible Planet Nine candidate, but are able to derive a detailed completeness curve, ruling out its presence within the parameter space searched at $W1 < 16.66$ (90% completeness). Our method uses all publicly available W1 imaging, spanning 2010 January to 2015 December, and will become more sensitive with future NEOWISE-Reactivation releases of additional W1 exposures. We anticipate that our method will be applicable to the entire high Galactic latitude sky, and we will extend our search to that full footprint in the near future.

We reconsider the commonly held assumption that warm debris disks are tracers of terrestrial planet formation. The high occurrence rate inferred for Earth-mass planets around mature solar-type stars based on exoplanet surveys (roughly 20%) stands in stark contrast to the low incidence rate (less than 2-3%) of warm dusty debris around solar-type stars during the expected epoch of terrestrial planet assembly (roughly 10 Myr). If Earth-mass planets at AU distances are a common outcome of the planet formation process, this discrepancy suggests that rocky planet formation occurs more quickly and/or is much neater than traditionally believed, leaving behind little in the way of a dust signature. Alternatively, the incidence rate of terrestrial planets has been overestimated or some previously unrecognized physical mechanism removes warm dust efficiently from the terrestrial planet region. A promising removal mechanism is gas drag in a residual gaseous disk with a surface density of roughly or somewhat more than 0.001% of the minimum mass solar nebula.

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.

Correlations in the orbits of several minor planets in the outer solar system suggest the presence of a remote, massive Planet Nine. With at least ten times the mass of the Earth and a perihelion well beyond 100 AU, Planet Nine poses a challenge to planet formation theory. Here we expand on a scenario in which the planet formed closer to the Sun and was gravitationally scattered by Jupiter or Saturn onto a very eccentric orbit in an extended gaseous disk. Dynamical friction with the gas then allowed the planet to settle in the outer solar system. We explore this possibility with a set of numerical simulations. Depending on how the gas disk evolves, scattered super-Earths or small gas giants settle on a range of orbits, with perihelion distances as large as 300 AU. Massive disks that clear from the inside out on million-year time scales yield orbits that allow a super-Earth or gas giant to shepherd the minor planets as observed. A massive planet can achieve a similar orbit in a persistent, low-mass disk over the lifetime of the solar system.

We analyze two new sets of coagulation calculations for solid particles orbiting within the terrestrial zone of a solar-type star. In models of collisional cascades, numerical simulations demonstrate that the total mass, the mass in 1 mm and smaller particles, and the dust luminosity decline with time more rapidly than predicted by analytic models, $\propto t^{-n}$ with $n \approx$ 1.1-1.2 instead of 1. Size distributions derived from the numerical calculations follow analytic predictions at radii less than 0.1 km but are shallower than predicted at larger sizes. In simulations of planet formation, the dust luminosity declines more slowly than in pure collisional cascades, with $n \approx$ 0.5-0.8 instead of 1.1-1.2. Throughout this decline, giant impacts produce large, observable spikes in dust luminosity which last roughly 0.01-0.1 Myr and recur every 1-10 Myr. If most solar-type stars have Earth mass planets with $a \lesssim$ 1-2 AU, observations of debris around 1-100 Myr stars allow interesting tests of theory. Current data preclude theories where terrestrial planets form out of 1000 km or larger planetesimals. Although the observed frequency of debris disks among $\gtrsim$ 30 Myr old stars agrees with our calculations, the observed frequency of warm debris among 5-20 Myr old stars is smaller than predicted.

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.

We consider the formation of satellites around the Pluto-Charon binary. An early collision between the two partners likely produced the binary and a narrow ring of debris, out of which arose the moons Styx, Nix, Kerberos and Hydra. How the satellites emerged from the compact ring is uncertain. Here we show that a particle ring spreads from physical collisions and collective gravitational scattering, similar to migration. Around a binary, these processes take place in the reference frames of "most circular" orbits, akin to circular ones in a Keplerian potential. Ring particles damp to these orbits and avoid destructive collisions. Damping and diffusion also help particles survive dynamical instabilities driven by resonances with the binary. In some situations, particles become trapped near resonances that sweep outward with the tidal evolution of the Pluto-Charon binary. With simple models and numerical experiments, we show how the Pluto-Charon impact ring may have expanded into a broad disk, out of which grew the circumbinary moons. In some scenarios, the ring can spread well beyond the orbit of Hydra, the most distant moon, to form a handful of smaller satellites. If these small moons exist, New Horizons will find them.

We examine characteristics of circumbinary orbits in the context of current planet formation scenarios. Analytical perturbation theory predicts the existence of nested circumbinary orbits that are generalizations of circular paths around a single star. These orbits have forced eccentric motion aligned with the binary as well as higher frequency oscillations, yet they do not cross, even in the presence of massive disks and perturbations from large planets. For this reason, dissipative gas and planetesimals can settle onto these "most circular" orbits, facilitating the growth of protoplanets. Outside a region close to the binary where orbits are generally unstable, circumbinary planets form in much the same way as their cousins around a single star. Here, we review the theory and confirm its predictions with a suite of representative simulations. We then consider the circumbinary planets discovered with NASA's Kepler satellite. These Neptune- and Jupiter-size planets, or their planetesimal precursors, may have migrated inward to reach their observed orbits, since their current positions are outside of unstable zones caused by overlapping resonances. In situ formation without migration seems less likely, only because the surface density of the protoplanetary disks must be implausibly high. Otherwise, the circumbinary environment is friendly to planet formation, and we expect that many Earth-like "Tatooines" will join the growing census of circumbinary planets.

We investigate pathways for the formation of icy super-Earth mass planets orbiting at 125-250 AU around a 1 solar mass star. An extensive suite of coagulation calculations demonstrates that swarms of 1 cm to 10 m planetesimals can form super-Earth mass planets on time scales of 1-3 Gyr. Collisional damping of 0.01-100 cm particles during oligarchic growth is a highlight of these simulations. In some situations, damping initiates a second runaway growth phase where 100-3000 km protoplanets grow to super-Earth sizes. Our results establish the initial conditions and physical processes required for in situ formation of super-Earth planets at large distances from the host star. For nearby dusty disks in HD 107146, HD 202628, and HD 207129, ongoing super-Earth formation at 80-150 AU could produce gaps and other structures in the debris. In the solar system, forming a putative planet X at a < 300 AU (a > 1000 AU) requires a modest (very massive) protosolar nebula.

As gas giant planets evolve, they may scatter other planets far from their original orbits to produce hot Jupiters or rogue planets that are not gravitationally bound to any star. Here, we consider planets cast out to large orbital distances on eccentric, bound orbits through a gaseous disk. With simple numerical models, we show that super-Earths can interact with the gas through dynamical friction to settle in the remote outer regions of a planetary system. Outcomes depend on planet mass, the initial scattered orbit, and the evolution of the time-dependent disk. Efficient orbital damping by dynamical friction requires planets at least as massive as the Earth. More massive, longer-lived disks damp eccentricities more efficiently than less massive, short-lived ones. Transition disks with an expanding inner cavity can circularize orbits at larger distances than disks that experience a global (homologous) decay in surface density. Thus, orbits of remote planets may reveal the evolutionary history of their primordial gas disks. A remote planet with an orbital distance ~100 AU from the Sun is plausible and might explain correlations in the orbital parameters of several distant trans-Neptunian objects.

We predict the distinctive three dimensional space motions of hypervelocity stars (HVSs) and runaway stars moving in a realistic Galactic potential. For nearby stars with distances less than 10~kpc, unbound stars are rare; proper motions alone rarely isolate bound HVSs and runaways from indigenous halo stars. At large distances of 20-100 kpc, unbound HVSs are much more common than runaways; radial velocities easily distinguish both from indigenous halo stars. Comparisons of the predictions with existing observations are encouraging. Although the models fail to match observations of solar-type HVS candidates from SEGUE, they agree well with data for B-type HVS and runaways from other surveys. Complete samples of g <= 20 stars with GAIA should provide clear tests of formation models for HVSs and runaways and will enable accurate probes of the shape of the Galactic potential.

The Milky Way can act as a large-scale weak gravitational lens of the cosmic microwave background (CMB). We study this effect using a photon ray-tracing code and a Galactic mass distribution with disk, bulge and halo components. For an observer at the Sun's coordinates in the Galaxy, the bending of CMB photon paths is limited to less than one arcsecond, and only for rays that pass within a few degrees of the Galactic Center. However, the entire sky is affected, resulting in global distortions of the CMB on large angular scales. These distortions can cause the low-order multipoles of a spherical harmonic expansion of the CMB sky temperature to leak into higher-order modes. Thus the component of the CMB dipole that results from the Local Group's motion relative to the local cosmic frame of rest contributes to higher-order moments for an observer in the solar system. With our ray-tracing code we show that the phenomenon is not sensitive to the specific choice of Galactic potential. We also quantitatively rule it out as a contributor to CMB anomalies such as power asymmetry or correlated alignment of low-order multipole moments.

We consider the ability of three models - impacts, captures, and collisional cascades - to account for a bright cloud of dust in Fomalhaut b. Our analysis is based on a novel approach to the power-law size distribution of solid particles central to each model. When impacts produce debris with (i) little material in the largest remnant and (ii) a steep size distribution, the debris has enough cross-sectional area to match observations of Fomalhaut b. However, published numerical experiments of impacts between 100 km objects suggest this outcome is unlikely. If collisional processes maintain a steep size distribution over a broad range of particle sizes (300 microns to 10 km), Earth-mass planets can capture enough material over 1-100 Myr to produce a detectable cloud of dust. Otherwise, capture fails. When young planets are surrounded by massive clouds or disks of satellites, a collisional cascade is the simplest mechanism for dust production in Fomalhaut b. Several tests using HST or JWST data - including measuring the expansion/elongation of Fomalhaut b, looking for trails of small particles along Fomalhaut b's orbit, and obtaining low resolution spectroscopy - can discriminate among these models.

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.

Motivated by the New Horizons mission, we consider how Pluto's small satellites -- currently Styx, Nix, Kerberos, and Hydra -- grow in debris from the giant impact that forms the Pluto-Charon binary. After the impact, Pluto and Charon accrete some of the debris and eject the rest from the binary orbit. During the ejection, high velocity collisions among debris particles produce a collisional cascade, leading to the ejection of some debris from the system and enabling the remaining debris particles to find stable orbits around the binary. Our numerical simulations of coagulation and migration show that collisional evolution within a ring or a disk of debris leads to a few small satellites orbiting Pluto-Charon. These simulations are the first to demonstrate migration-induced mergers within a particle disk. The final satellite masses correlate with the initial disk mass. More massive disks tend to produce fewer satellites. For the current properties of the satellites, our results strongly favor initial debris masses of 3-10 x 10^{19} g and current satellite albedos of roughly 0.4-1. We also predict an ensemble of smaller satellites with radii less than roughly 1-3 km, and very small particles, with radii of roughly 1-100 cm and optical depth smaller than 10^-10. These objects should have semimajor axes outside the current orbit of Hydra.

We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet $\beta$ Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve high signal-to-noise, high precision measurements. While $\beta$ Pic b's near-IR colors mimick that of a standard, cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared compared to the field L/T dwarf sequence. Few substellar/planet-mass objects -- i.e. $\kappa$ And b and 1RXJ 1609B -- match $\beta$ Pic b's $JHK_{s}L^\prime$ photometry, and its 3.1 $\mu m$ and 5 $\mu m$ photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large ($\sim$ 60 $\mu m$) dust grains fail to reproduce the $\beta$ Pic b spectrum. However, models incorporating thick clouds similar to those found for HR 8799 bcde but also with small (a few microns) modal particle sizes yield fits consistent with the data within uncertainties. Assuming solar abundance models, thick clouds, and small dust particles ($<a>$ = 4 $\mu m$) we derive atmosphere parameters of log(g) = 3.8 $\pm$ 0.2 and $T_{eff}$ = 1575--1650 $K$, an inferred mass of 7$^{+4}_{-3}$ $M_{J}$, and a luminosity of log(L/L$_{\odot}$) $\sim$ -3.80 $\pm$ 0.02. The best-estimated planet radius, $\approx$ 1.65 $\pm$ 0.06 $R_{J}$, is near the upper end of allowable planet radii for hot-start models given the host star's age and likely reflects challenges with constructing accurate atmospheric models. Alternatively, these radii are comfortably consistent with hot-start model predictions if $\beta$ Pic b is younger than $\approx$ 7 Myr, consistent with a late formation, well after its host star's birth $\sim$ 12$^{+8}_{-4}$ Myr ago.

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.

We examine whether disrupted binary stars can fuel black hole growth. In this mechanism, tidal disruption produces a single hypervelocity star (HVS) ejected at high velocity and a former companion star bound to the black hole. After a cluster of bound stars forms, orbital diffusion allows the black hole to accrete stars by tidal disruption at a rate comparable to the capture rate. In the Milky Way, HVSs and the S star cluster imply similar rates of 10^{-5}--10^{-3} yr^{-1} for binary disruption. These rates are consistent with estimates for the tidal disruption rate in nearby galaxies and imply significant black hole growth from disrupted binaries on 10 Gyr time scales.

We investigate whether coagulation models of planet formation can explain the observed size distributions of transneptunian objects (TNOs). Analyzing published and new calculations, we demonstrate robust relations between the size of the largest object and the slope of the size distribution for sizes 0.1 km and larger. These relations yield clear, testable predictions for TNOs and other icy objects throughout the solar system. Applying our results to existing observations, we show that a broad range of initial disk masses, planetesimal sizes, and fragmentation parameters can explain the data. Adding dynamical constraints on the initial semimajor axis of `hot' KBOs along with probable TNO formation times of 10-700 Myr restricts the viable models to those with a massive disk composed of relatively small (1-10 km) planetesimals.

If dark matter self-annihilates then it may produce an observable signal when its density is high. The details depend on the intrinsic properties of dark matter and how it clusters in space. For example, the density profile of some dark matter candidates may rise steeply enough toward the Galactic Center that self-annihilation produces detectable gamma-ray emission. Here, we discuss the possibility that an annihilation signal may arise near a compact object (e.g., neutron star or black hole) even when the density of dark matter in the neighborhood of the object is uniform. Gravitational focusing produces a local enhancement of density, with a profile that falls off approximately as the inverse square-root of distance from the compact star. While geometric dilution may overwhelm the annihilation signal from this local enhancement, magnetic fields tied to the compact object can increase the signal's contrast relative to the background.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

Planetary migration poses a serious challenge to theories of planet formation. In gaseous and planetesimal disks, migration can remove planets as quickly as they form. To explore migration in a planetesimal disk, we combine analytic and numerical approaches. After deriving general analytic migration rates for isolated planets, we use N-body simulations to confirm these results for fast and slow migration modes. Migration rates scale as 1/m (for massive planets) and 1/(1+(e_H/3)^3), where m is the mass of a planet and e_H is the eccentricity of the background planetesimals in Hill units. When multiple planets stir the disk, our simulations yield the new result that large-scale migration ceases. Thus, growing planets do not migrate through planetesimal disks. To extend these results to migration in gaseous disks, we compare physical interactions and rates. Although migration through a gaseous disk is an important issue for the formation of gas giants, we conclude that migration has little impact on the formation of terrestrial planets.

We describe an updated version of our hybrid N-body-coagulation code for planet formation. In addition to the features of our 2006-2008 code, our treatment now includes algorithms for the 1D evolution of the viscous disk, the accretion of small particles in planetary atmospheres, gas accretion onto massive cores, and the response of N-bodies to the gravitational potential of the gaseous disk and the swarm of planetesimals. To validate the N-body portion of the algorithm, we use a battery of tests in planetary dynamics. As a first application of the complete code, we consider the evolution of Pluto-mass planetesimals in a swarm of 0.1-1 cm pebbles. In a typical evolution time of 1-3 Myr, our calculations transform 0.01-0.1 solar mass disks of gas and dust into planetary systems containing super-Earths, Saturns, and Jupiters. Low mass planets form more often than massive planets; disks with smaller alpha form more massive planets than disks with larger alpha. For Jupiter-mass planets, masses of solid cores are 10-100 Earth masses.

We describe comprehensive calculations of the formation of icy planets and debris disks at 30-150 AU around 1-3 solar mass stars. Disks composed of large, strong planetesimals produce more massive planets than disks composed of small, weak planetesimals. The maximum radius of icy planets ranges from roughly 1500 km to 11,500 km. The formation rate of 1000 km objects - `Plutos' - is a useful proxy for the efficiency of icy planet formation. Plutos form more efficiently in massive disks, in disks with small planetesimals, and in disks with a range of planetesimal sizes. Although Plutos form throughout massive disks, Pluto production is usually concentrated in the inner disk. Despite the large number of Plutos produced in many calculations, icy planet formation is inefficient. At the end of the main sequence lifetime of the central star, Plutos contain less than 10% of the initial mass in solid material. This conclusion is independent of the initial mass in the disk or the properties of planetesimals. Debris disk formation coincides with the formation of planetary systems containing Plutos. As Plutos form, they stir leftover planetesimals to large velocities. A cascade of collisions then grinds the leftovers to dust, forming an observable debris disk. In disks with small (< 1-10 km) planetesimals, collisional cascades produce luminous debris disks with maximum luminosity roughly 0.01 times the stellar luminosity. Disks with larger planetesimals produce much less luminous debris disks. Observations of debris disks around A-type and G-type stars strongly favor models with small planetesimals. In these models, our predictions for the time evolution and detection frequency of debris disks agree with published observations. We suggest several critical observations that can test key features of our calculations.

Runaway stars ejected from the Galactic disk populate the halo of the Milky Way. To predict the spatial and kinematic properties of runaways, we inject stars into a Galactic potential, compute their trajectories through the Galaxy, and derive simulated catalogs for comparison with observations. Runaways have a flattened spatial distribution, with higher velocity stars at Galactic latitudes less than 30 degrees. Due to their shorter stellar lifetimes, massive runaway stars are more concentrated towards the disk than low mass runaways. Bound (unbound) runaways that reach the halo probably originate from distances of 6--12 kpc (10--15 kpc) from the Galactic center, close to the estimated origin of the unbound runaway star HD 271791. Because runaways are brighter and have smaller velocities than hypervelocity stars (HVSs), radial velocity surveys are unlikely to confuse runaway stars with HVSs. We estimate that at most 1 runaway star contaminates the current sample. We place an upper limit of 2% on the fraction of A-type main sequence stars ejected as runaways.

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.

We study the distribution of angular positions and angular separations of unbound hypervelocity stars (HVSs). HVSs are spatially anisotropic at the 3-sigma level. The spatial anisotropy is significant in Galactic longitude, not in latitude, and the inclusion of lower velocity, possibly bound HVSs reduces the significance of the anisotropy. We discuss how the observed distribution of HVSs may be linked to their origin. In the future, measuring the distribution of HVSs in the southern sky will provide additional constraints on the spatial anisotropy and the origin of HVSs.

We describe calculations for the formation of icy planets and debris disks at 30-150 AU around 1-3 solar mass stars. Debris disk formation coincides with the formation of planetary systems. As protoplanets grow, they stir leftover planetesimals to large velocities. A cascade of collisions then grinds the leftovers to dust, forming an observable debris disk. Stellar lifetimes and the collisional cascade limit the growth of protoplanets. The maximum radius of icy planets, roughly 1750 km, is remarkably independent of initial disk mass, stellar mass, and stellar age. These objects contain no more than 3% to 4% of the initial mass in solid material. Collisional cascades produce debris disks with maximum luminosity of roughly 0.002 times the stellar luminosity. The peak 24 micron excess varies from roughly 1% of the stellar photospheric flux for 1 solar mass stars to roughly 50 times the stellar photospheric flux for 3 solar mass stars. The peak 70-850 micron excesses are roughly 30-100 times the stellar photospheric flux. For all stars, the 24-160 micron excesses rise at stellar ages of 5-20 Myr, peak at 10-50 Myr, and then decline. The decline is roughly a power law, with f propto t^{-n} and n = 0.6-1.0. This predicted evolution agrees with published observations of A-type and solar-type stars. The observed far-IR color evolution of A-type stars also matches model predictions.

Hypervelocity stars (HVS) traverse the Galaxy from the central black hole to the outer halo. We show that the Galactic potential within 200 pc acts as a high pass filter preventing low velocity HVS from reaching the halo. To trace the orbits of HVS throughout the Galaxy, we construct two forms of the potential which reasonably represent the observations in the range 5--100,000 pc, a simple spherically symmetric model and a bulge-disk-halo model. We use the Hills mechanism (disruption of binaries by the tidal field of the central black hole) to inject HVS into the Galaxy and compute the observable spatial and velocity distributions of HVS with masses in the range 0.6--4 Msun. These distributions reflect the mass function in the Galactic Center, properties of binaries in the Galactic Center, and aspects of stellar evolution and the injection mechanism. For 0.6--4 Msun main sequence stars, the fraction of unbound HVS and the asymmetry of the velocity distribution for their bound counterparts increases with stellar mass. The density profiles for unbound HVS decline with distance from the Galactic Center approximately as r^{-2} (but are steeper for the most massive stars which evolve off the main sequence during their travel time from the Galactic Center); the density profiles for the bound ejecta decline with distance approximately as r^{-3}. In a survey with a limiting visual magnitude V of 23, the detectability of HVS (unbound or bound) increases with stellar mass.

We describe Spitzer/MIPS observations of the double cluster, h and $\chi$ Persei, covering a $\sim$ 0.6 square-degree area surrounding the cores of both clusters. The data are combined with IRAC and 2MASS data to investigate $\sim$ 616 sources from 1.25-24 $\mu m$. We use the long-baseline $K_{s}$-[24] color to identify two populations with IR excess indicative of circumstellar material: Be stars with 24 $\mu m$ excess from optically-thin free free emission and 17 fainter sources (J$\sim$ 14-15) with [24] excess consistent with a circumstellar disk. The frequency of IR excess for the fainter sources increases from 4.5 $\mu m$ through 24 $\mu m$. The IR excess is likely due to debris from the planet formation process. The wavelength-dependent behavior is consistent with an inside-out clearing of circumstellar disks. A comparison of the 24 $\mu m$ excess population in h and $\chi$ Per sources with results for other clusters shows that 24 $\mu m$ emission from debris disks 'rises' from 5 to 10 Myr, peaks at $\sim$ 10-15 Myr, and then 'falls' from $\sim$ 15/20 Myr to 1 Gyr.

We report the discovery of 3 new unbound hypervelocity stars (HVSs), stars traveling with such extreme velocities that dynamical ejection from a massive black hole (MBH) is their only suggested origin. We also detect a population of possibly bound HVSs. The significant asymmetry we observe in the velocity distribution -- we find 26 stars with v_rf > 275 km/s and 1 star with v_rf < -275 km/s -- shows that the HVSs must be short-lived, probably 3 - 4 Msun main sequence stars. Any population of hypervelocity post-main sequence stars should contain stars falling back onto the Galaxy, contrary to the observations. The spatial distribution of HVSs also supports the main sequence interpretation: longer-lived 3 Msun HVSs fill our survey volume; shorter-lived 4 Msun HVSs are missing at faint magnitudes. We infer that there are 96 +- 10 HVSs of mass 3 - 4 Msun within R < 100 kpc, possibly enough HVSs to constrain ejection mechanisms and potential models. Depending on the mass function of HVSs, we predict that SEGUE may find up to 5 - 15 new HVSs. The travel times of our HVSs favor a continuous ejection process, although a ~120 Myr-old burst of HVSs is also allowed.

We analyze 8 sources with strong mid-infrared excesses in the 13 Myr-old double cluster h and chi Persei. New optical spectra and broadband SEDs (0.36-8 mu_m) are consistent with cluster membership. We show that material with T ~ 300-400 K and Ld/Lstar ~ 10^-4-10^-3 produces the excesses in these sources. Optically-thick blackbody disk models - including those with large inner holes - do not match the observed SEDs. The SEDs of optically-thin debris disks produced from terrestrial planet formation calculations match the observations well. Thus, some h and chi Persei stars may have debris from terrestrial zone planet formation.

This chapter summarizes analytic theory and numerical calculations for the formation and collisional evolution of KBOs at 20--150 AU. We describe the main predictions of a baseline self-stirring model and show how dynamical perturbations from a stellar flyby or stirring by a giant planet modify the evolution. Although robust comparisons between observations and theory require better KBO statistics and more comprehensive calculations, the data are broadly consistent with KBO formation in a massive disk followed by substantial collisional grinding and dynamical ejection. However, there are important problems reconciling the results of coagulation and dynamical calculations. Contrasting our current understanding of the evolution of KBOs and asteroids suggests that additional observational constraints, such as the identification of more dynamical families of KBOs (like the 2003 EL61 family), would provide additional information on the relative roles of collisional grinding and dynamical ejection in the Kuiper Belt. The uncertainties also motivate calculations that combine collisional and dynamical evolution, a `unified' calculation that should give us a better picture of KBO formation and evolution.

Hypervelocity stars (HVSs) are stars ejected completely out of the Milky Way by three-body interactions with the massive black hole in the Galactic center. We describe 643 new spectroscopic observations from our targeted survey for HVSs. We find a significant (3.5 sigma) excess of B-type stars with large velocities +275<v_rf<450 km/s and distances d>10 kpc that are most plausibly explained as a new class of HVSs: stars ejected from the Galactic center on bound orbits. If a Galactic center ejection origin is correct, the distribution of HVSs on the sky should be anisotropic for a survey complete to a fixed limiting apparent magnitude. The unbound HVSs in our survey have a marginally anisotropic distribution on the sky, consistent with the Galactic center ejection picture.

We develop a semi-analytic model for planet formation during the pre-main sequence contraction phase of a low mass star. During this evolution, the stellar magnetosphere maintains a fixed ratio between the inner disk radius and the stellar radius. As the star contracts at constant effective temperature, the `snow line', which separates regions of rocky planet formation from regions of icy planet formation, moves inward. This process enables rapid formation of icy protoplanets that collide and merge into super-Earths before the star reaches the main sequence. The masses and orbits of these super-Earths are consistent with super-Earths detected in recent microlensing experiments.

The disruption of binary stars by the tidal field of the black hole in the Galactic Center can produce the hypervelocity stars observed in the halo. We use numerical models to simulate the full spectrum of observable velocities of stars ejected into the halo by this binary disruption process. Our model includes a range of parameters for binaries with 3-4 M_Solar primaries, consideration of radial orbits of the ejected stars through an approximate mass distribution for the Galaxy, and the impact of stellar lifetimes. We calculate the spectrum of ejection velocities and reproduce previous results for the mean ejection velocity at the Galactic center. The model predicts that the full population of ejected stars includes both the hypervelocity stars with velocities large enough to escape from the Galaxy and a comparable number of ejected, but bound, stars of the same stellar type. The predicted median speeds of the population of ejected stars as a function of distance in the halo are consistent with current observations. Combining the model with the data also shows that interesting constraints on the properties of binaries in the Galactic Center and on the mass distribution in the Galaxy can be obtained even with modest samples of ejected stars.

We describe a hybrid algorithm to calculate the formation of planets from an initial ensemble of planetesimals. The algorithm uses a coagulation code to treat the growth of planetesimals into oligarchs and explicit N-body calculations to follow the evolution of oligarchs into planets. To validate the N-body portion of the algorithm, we use a battery of tests in planetary dynamics. Several complete calculations of terrestrial planet formation with the hybrid code yield good agreement with previously published calculations. These results demonstrate that the hybrid code provides an accurate treatment of the evolution of planetesimals into planets.

We use a hybrid, multiannulus, n-body-coagulation code to investigate the growth of km-sized planetesimals at 0.4-2 AU around a solar-type star. After a short runaway growth phase, protoplanets with masses of roughly 10^26 g and larger form throughout the grid. When (i) the mass in these `oligarchs' is roughly comparable to the mass in planetesimals and (ii) the surface density in oligarchs exceeds 2-3 g/sq cm at 1 AU, strong dynamical interactions among oligarchs produce a high merger rate which leads to the formation of several terrestrial planets. In disks with lower surface density, milder interactions produce several lower mass planets. In all disks, the planet formation timescale is roughly 10-100 Myr, similar to estimates derived from the cratering record and radiometric data.

We investigate the prospects for detecting dust from two body collisions during the late stages of planet formation at 1-150 AU. We develop an analytic model to describe the formation of a dusty cloud of debris and use numerical coagulation and n-body calculations to predict observable signals from these events. In a minimum mass solar nebula, collisions of 100-1000 km objects at distances of 3-5 AU or less from the parent star are observable at mid-infrared wavelengths as bright clumps or rings of dust. At 24 microns, the clumps are roughly 0.1-1 mag brighter than emission from dust in the background debris disk. In edge-on systems, dusty clumps produce eclipses with depths of 1.0 mag or less that last for roughly 100 orbital periods. Large-scale surveys for transits from exosolar planets, such as Kepler, can plausibly detect these eclipses and provide important constraints on the terrestrial environment for ages of less than or roughly 100-300 Myr.

The discovery of Sedna places new constraints on the origin and evolution of our solar system. Here we investigate the possibility that a close encounter with another star produced the observed edge of the Kuiper belt, at roughly 50 AU, and the highly elliptical orbit of Sedna. We show that a passing star probably scattered Sedna from the Kuiper Belt into its observed orbit. The likelihood that a planet at 60-80 AU can be scattered into Sedna's orbit is roughly 50%; this estimate depends critically on the geometry of the flyby. Even more interesting, though, is the roughly 10% chance that Sedna was captured from the outer disk of the passing star. Most captures have very high inclination orbits; detection of these objects would confirm the presence of extrasolar planets in our own Solar System.

We describe analytical and numerical collisional evolution calculations for the size distribution of icy bodies in the Kuiper Belt. For a wide range of bulk properties, initial masses, and orbital parameters, our results yield power-law cumulative size distributions, N_C propto r^{-q}, with q_L = 3.5 for large bodies with radii of 10-100 km, and q_s = 2.5-3 for small bodies with radii lesss than 0.1-1 km. The transition between the two power laws occurs at a break radius of 1-30 km. The break radius is more sensitive to the initial mass in the Kuiper Belt and the amount of stirring by Neptune than the bulk properties of individual Kuiper Belt objects (KBOs). Comparisons with observations indicate that most models can explain the observed sky surface density of KBOs for red magnitudes, R = 22-27. For R < 22 and R > 28, the model surface density is sensitive to the amount of stirring by Neptune, suggesting that the size distribution of icy planets in the outer solar system provides independent constraints on the formation of Neptune.

We use a multiannulus accretion code to investigate debris disks in the terrestrial zone, at 0.7-1.3 AU around a 1 solar mass star. Terrestrial planet formation produces a bright dusty ring of debris with a lifetime of at least 1 Myr. The early phases of terrestrial planet formation are observable with current facilities; the late stages require more advanced instruments with adaptive optics.

We use a multiannulus planetesimal accretion code to investigate the growth of icy planets in the outer regions of a planetesimal disk. In a quiescent minimum mass solar nebula, icy planets grow to sizes of 1000--3000 km on a timescale t = 15-20 Myr (a/30 AU)^3 where a is the distance from the central star. Planets form faster in more massive nebulae. Newly-formed planets stir up leftover planetesimals along their orbits and produce a collisional cascade where icy planetesimals are slowly ground to dust. The dusty debris of planet formation has physical characteristics similar to those observed in beta Pic, HR 4796A, and other debris disks. We derive dust masses for small particles, 1 mm and smaller, and large particles, 1 mm and larger, as a function of the initial conditions in the planetesimal disk. The dust luminosities derived from these masses are similar to those observed in Vega, HR 4796A, and other debris disks. The calculations produce bright rings and dark gaps. Bright rings occur where 1000 km and larger planets have recently formed. Dark gaps are regions where planets have cleared out dust or shadows where planets have yet to form.

Many nearby stars are surrounded by a bright ring or disk of cold dust. Our calculations show that these disks and rings of dust are signposts of recent planet formation. Bright rings appear because dust associated with the formation of a planet absorbs and scatters light from the central star. The calculations explain the rings observed so far and predict that all nascent solar systems have dusty rings.

We use a new multiannulus planetesimal accretion code to investigate the evolution of a planetesimal disk following a moderately close encounter with a passing star. The calculations include fragmentation, gas and Poynting-Robertson drag, and velocity evolution from dynamical friction and viscous stirring. We assume that the stellar encounter increases planetesimal velocities to the shattering velocity, initiating a collisional cascade in the disk. During the early stages of our calculations, erosive collisions damp particle velocities and produce substantial amounts of dust. For a wide range of initial conditions and input parameters, the time evolution of the dust luminosity follows a simple relation, L_d/L_{\star} = L_0 / [alpha + (t/t_d)^{beta}]. The maximum dust luminosity L_0 and the damping time t_d depend on the disk mass, with L_0 proportional to M_d and t_d proportional to M_d^{-1}. For disks with dust masses of 1% to 100% of the `minimum mass solar nebula' (1--100 earth masses at 30--150 AU), our calculations yield t_d approx 1--10 Myr, alpha approx 1--2, beta = 1, and dust luminosities similar to the range observed in known `debris disk' systems, L_0 approx 10^{-3} to 10^{-5}. Less massive disks produce smaller dust luminosities and damp on longer timescales. Because encounters with field stars are rare, these results imply that moderately close stellar flybys cannot explain collisional cascades in debris disk systems with stellar ages of 100 Myr or longer.

The radio source Sgr A* in the Galactic center emits a polarized spectrum at millimeter and sub-millimeter wavelengths that is strongly suggestive of relativistic disk accretion onto a massive black hole. We use the well-constrained mass of Sgr A* and a magnetohydrodynamic model of the accretion flow to match both the total flux and polarization from this object. Our results demonstrate explicitly that the shift in the position angle of the polarization vector, seen at wavelengths near the peak of the mm to sub-mm emission from this source, is a signal of relativistic accretion flow in a strong gravitational field. We provide maps of the polarized emission to illustrate how the images of polarized intensity from the vicinity of the black hole would appear in upcoming observations with very long baseline radio interferometers (VLBI). Our results suggest that near-term VLBI observations will be able to directly image the polarized Keplerian portion of the flow near the horizon of the black hole.

The polarized mm/sub-mm radiation from Sgr A* is apparently produced by a Keplerian structure whose peak emission occurs within several Schwarzschild radii (r_S=2GM/c^2) of the black hole. The Chandra X-ray counterpart, if confirmed, is presumably the self-Comptonized component from this region. In this paper, we suggest that sub-mm timing observations could yield a signal corresponding to the period P_0 of the marginally stable orbit, and therefore point directly to the black hole's spin a. Sgr A*'s mass is now known to be (2.6\pm 0.2)\times 10^6 M_\odot (an unusually accurate value for supermassive black hole candidates), for which 2.7 min<P_0<36 min, depending on the value of a and whether the Keplerian flow is prograde or retrograde. A Schwarzschild black hole (a=0) should have P_0 ~ 20 min. The identification of the orbital frequency with the innermost stable circular orbit is made feasible by the transition from optically thick to thin emission at sub-mm wavelengths. With stratification in the emitter, the peak of the sub-mm bump in Sgr A*'s spectrum is thus produced at the smallest radius. We caution, however, that theoretical uncertainties in the structure of the emission region may still produce some ambiguity in the timing signal. Given that Sgr A*'s flux at $\nu\sim 1$ mm is several Jy, these periods should lie within the temporal-resolving capability of sub-mm telescopes using bolometric detectors. A determination of P_0 should provide not only a value of a, but it should also define the angular momentum vector of the orbiting gas in relation to the black hole's spin axis. In addition, since the X-ray flux detected by Chandra appears to be the self-Comptonized mm to sub-mm component, these temporal fluctuations may also be evident in the X-ray signal.

We use optical rotation curves to investigate the R-band Tully-Fisher properties of a sample of 90 spiral galaxies in close pairs. The galaxies follow the Tully-Fisher relation remarkably well, with the exception of eight distinct 3-sigma outliers. Although most of the outliers show signs of recent star formation, gasdynamical effects are probably the dominant cause of their anomalous Tully-Fisher properties. Four outliers with small emission line widths have very centrally concentrated line emission and truncated rotation curves; the central emission indicates recent gas infall after a close galaxy-galaxy pass. These four galaxies may be local counterparts to compact, blue galaxies at intermediate redshift. The remaining galaxies have a negligible offset from the reference Tully-Fisher relation, but a shallower slope (2.6-sigma significance) and a 25% larger scatter. We characterize the non-outlier sample with measures of distortion and star formation to search for third parameter dependence in the residuals of the TF relation. Severe kinematic distortion is the only significant predictor of TF residuals; this distortion is not, however, responsible for the slope difference from the reference distribution. Because the outliers are easily removed by sigma clipping, we conclude that even in the presence of some tidal distortion, detection of moderate luminosity evolution should be possible with high-redshift samples the size of this 90-galaxy study. (Abridged.)

The observed iron K-alpha fluorescence lines in Seyfert-1 galaxies provide strong evidence for an accretion disk near a supermassive black hole as a source of the emission. Here we present an analysis of the geometrical and kinematic properties of the disk based on the extreme frequency shifts of a line profile as determined by measurable flux in both the red and blue wings. The edges of the line are insensitive to the distribution of the X-ray flux over the disk, and hence provide a robust alternative to profile fitting of disk parameters. Our approach yields new, strong bounds on the inclination angle of the disk and the location of the emitting region. We apply our method to interpret observational data from MCG-6-30-15 and find that the commonly assumed inclination 30 deg for the accretion disk in MCG-6-30-15 is inconsistent with the position of the blue edge of the line at a 3 sigma level. A thick turbulent disk model or the presence of highly ionized iron may reconcile the bounds on inclination from the line edges with the full line profile fits based on simple, geometrically thin disk models. The bounds on the innermost radius of disk emission indicate that the black hole in MCG-6-30-15 is rotating faster than 30 % of theoretical maximum. When applied to data from NGC 4151, our method gives bounds on the inclination angle of the X-ray emitting inner disk of 50\pm 10 deg, consistent with the presence of an ionization cone grazing the disk as proposed by Pedlar et al. (1993). The frequency extrema analysis also provides limits to the innermost disk radius in another Seyfert-1 galaxy, NGC 3516.

We describe gravitational stirring models of planetary debris disks using a new multi-annulus planetesimal evolution code. The current code includes gravitational stirring and dynamical friction; future studies will include coagulation, fragmentation, Poynting-Robertson drag, and other physical processes. We use the results of our calculations to investigate the physical conditions required for small bodies in a planetesimal disk to reach the shattering velocity and begin a collisional cascade. Our results demonstrate that disks composed primarily of bodies with a single size will not undergo a collisional cascade which produces small dust grains at 30-150 AU on timescales of 1 Gyr or smaller. Disks with a size distribution of bodies reach conditions necessary for a collisional cascade in 10 Myr to 1 Gyr if the disk is at least as massive as a minimum mass solar nebula and if the disk contains objects with radii of 500 km or larger. The estimated 500 Myr survival time for these disks is close to the median age of roughly 400 Myr derived for nearby stars with dusty disks.

Two recent papers have claimed detection of non-Gaussian features in the COBE DMR sky maps of the cosmic microwave background. We confirm these results, but argue that Gaussianity is still not convincingly ruled out. Since a score of non-Gaussianity tests have now been published, one might expect some mildly significant results even by chance. Moreover, in the case of one measure which yields a detection, a bispectrum statistic, we find that if the non-Gaussian feature is real, it may well be due to detector noise rather than a non-Gaussian sky signal, since a signal-to-noise analysis localizes it to angular scales smaller than the beam. We study its spatial origin in case it is nonetheless due to a sky signal (eg, a cosmic string wake or flat-spectrum foreground contaminant). It appears highly localized in the direction b=39.5, l=257, since removing a mere 5 pixels inside a single COBE beam area centered there makes the effect statistically insignificant. We also test Guassianity with an eigenmode analysis which allows a sky map to be treated as a random number generator. A battery of tests of this generator all yield results consistent with Gaussianity.

We show that the galaxy density in the Las Campanas Redshift Survey (LCRS) cannot be perfectly correlated with the underlying mass distribution since various galaxy subpopulations are not perfectly correlated with each other, even taking shot noise into account. This rules out the hypothesis of simple linear biasing, and suggests that the recently proposed stochastic biasing framework is necessary for modeling actual data.

We use self-consistent N-body models, in conjunction with models of test particles moving in galaxy potentials, to explore the initial effects of interactions on the rotation curves of spiral galaxies. Using nearly self-consistent disk/bulge/halo galaxy models (Kuijken & Dubinski 1995), we simulate the first pass of galaxies on nearly parabolic orbits; we vary orbit inclinations, galaxy halo masses and impact parameters. For each simulation, we mimic observed rotation curves of the model galaxies. Transient interaction-induced features of the curves include distinctly rising or falling profiles at large radii and pronounced bumps in the central regions. Remarkably similar features occur in our statistical sample of optical emission-line rotation curves of spiral galaxies in tight pairs and n-tuples.

The observed iron K-alpha fluorescence lines in Seyfert-1 galaxies provide strong evidence for an accretion disk near a supermassive black hole as a source of the line emission. These lines serve as powerful probes for examining the structure of inner regions of accretion disks. Previous studies of line emission have considered geometrically thin disks only, where the gas moves along geodesics in the equatorial plane of a black hole. Here we extend this work to consider effects on line profiles from finite disk thickness, radial accretion flow and turbulence. We adopt the Novikov and Thorne (1973) solution, and find that within this framework, turbulent broadening is the dominant new effect. The most prominent change in the skewed, double-horned line profiles is a substantial reduction in the maximum flux at both red and blue peaks. The effect is most pronounced when the inclination angle is large, and when the accretion rate is high. Thus, the effects discussed here may be important for future detailed modeling of high quality observational data.

Galaxies in the Las Campanas Redshift Survey are classified according to their spectra, and the resulting spectral types are analyzed to determine if local environment affects their properties. We find that the luminosity function of early-type objects varies as a function of local density. Our results suggest that early-type galaxies (presumably ellipticals and S0's) are, on average, fainter when they are located in high-density regions of the Universe. The same effect may operate for some, but not all, late-type galaxies. We discuss the implications of this result for theories of galaxy formation and evolution.

We construct a spectral classification scheme for the galaxies of the Las Campanas Redshift Survey (LCRS) based on a principal component analysis of the measured galaxy spectra. We interpret the physical significance of our six spectral types and conclude that they are sensitive to morphological type and the amount of active star formation. In this first analysis of the LCRS to include spectral classification, we estimate the general luminosity function, expressed as a weighted sum of the type-specific luminosity functions. In the R-band magnitude range of -23 < M <= -16.5, this function exhibits a broad shoulder centered near M = -20, and an increasing faint-end slope which formally converges on an alpha value of about -1.8 in the faint limit. The Schechter parameterization does not provide a good representation in this case, a fact which may partly explain the reported discrepancy between the luminosity functions of the LCRS and other redshift catalogs such as the Century Survey (Geller et al. 1997). The discrepancy may also arise from environmental effects such as the density-morphology relationship for which we see strong evidence in the LCRS galaxies. However, the Schechter parameterization is more effective for the luminosity functions of the individual spectral types. The data show a significant, progressive steepening of the faint-end slope, from alpha = +0.5 for early-type objects, to alpha = -1.8 for the extreme late-type galaxies. The extreme late-type population has a sufficiently high space density that its contribution to the general luminosity function is expected to dominate fainter than M = -16. We conclude that an evaluation of type-dependence is essential to any assessment of the general luminosity function.

The observed iron K-alpha fluorescence lines in Seyfert galaxies provide strong evidence for an accretion disk near a supermassive black hole as a source of the line emission. Previous studies of line emission have considered only geometrically thin disks, where the gas moves along geodesics in the equatorial plane of a black hole. Here we extend this work to include effects on line profiles from finite disk thickness, radial accretion flow and turbulence. We adopt the Novikov-Thorne solution, and find that within this framework, turbulent broadening is the most significant effect. The most prominent changes in the skewed, double-horned line profiles is a substantial reduction in the maximum flux at both red and blue peaks. We show that at the present level of signal-to-noise in X-ray spectra, proper treatment of the actual structure of the accretion disk can change estimates of the inclination angle of the disk. Thus these effects will be important for future detailed modeling of high quality observational data.

The observed profile of spectral lines from a relativistic accretion disk can constrain parameters such as the disk geometry and the rotation of the central black hole. The formation of the spectral line in a disk generally has been modeled with simple assumptions such as local isotropy of emission. Here we consider line formation in the presence of velocity gradients induced by the differential flow in the disk. In this case the emission can have anisotropy in the form of an azimuthal dependence relative to the local principle axes of shear. Since the physical conditions in a disk are uncertain in detail, we investigate this effect with simple parameterized models based on Sobolev theory to highlight the overall character of the changes in the line profile. We find that velocity gradients generally cause a relative increase of flux in the red wing, hence the inner radius of the disk would be underestimated if the effect were not taken into consideration. If the inner radius is used as a signature of black hole rotation, as when the disk is not emissive within the marginally stable circular orbit, then the inferred rotation would be overestimated in cases where the emissivity of the disk has fairly shallow fall-off with radius. If the disk were emissive even within the marginally stable orbit, then the local azimuthal anisotropy of emission produces features in the line profile which distinguish rotating from nonrotating black holes.

Line emission from an accretion disk and a corotating hot spot about a rotating black hole are considered for possible signatures of the frame-dragging effect. We explicitly compare integrated line profiles from a geometrically thin disk about a Schwarzschild and an extreme Kerr black hole, and show that the line profile differences are small if the inner radius of the disk is near or above the Schwarzschild stable-orbit limit of radius 6GM/c^2. However, if the inner disk radius extends below this limit, as is possible in the extreme Kerr spacetime, then differences can become significant, especially if the disk emissivity is stronger near the inner regions. We demonstrate that the first three moments of a line profile define a three-dimensional space in which the presence of material at small radii becomes quantitatively evident in broad classes of disk models. In the context of the simple, thin disk paradigm, this moment-mapping scheme suggests formally that the iron line detected by the Advanced Satellite for Cosmology and Astrophysics mission from MCG-6-30-15 (Tanaka et al. 1995) is 3 times more likely to originate from a disk about a rotating black hole than from a Schwarzschild system. A statistically significant detection of black hole rotation in this way may be achieved after only modest improvements in the quality of data. We also consider light curves and frequency shifts in line emission as a function of time for corotating hot spots in extreme Kerr and Schwarzschild geometries. Both the frequency-shift profile and the light curve from a hot spot are valuable measures of orbital parameters and might possibly be used to detect frame dragging even at radii approaching 6GM/c^2 if the inclination angle of the orbital plane is large.

We propose a method to determine the cosmic mass density Omega from redshift-space distortions induced by large-scale flows in the presence of nonlinear clustering. Nonlinear structures in redshift space such as fingers of God can contaminate distortions from linear flows on scales as large as several times the small-scale pairwise velocity dispersion sigma_v. Following Peacock & Dodds (1994), we work in the Fourier domain and propose a model to describe the anisotropy in the redshift-space power spectrum; tests with high-resolution numerical data demonstrate that the model is robust for both mass and biased galaxy halos on translinear scales and above. On the basis of this model, we propose an estimator of the linear growth parameter beta = Omega^0.6/b, where b measures bias, derived from sampling functions which are tuned to eliminate distortions from nonlinear clustering. The measure is tested on the numerical data and found to recover the true value of beta to within ~10%. An analysis of the IRAS 1.2Jy galaxies yields beta = 0.8+0.4/-0.3 at a scale of 1,000 km/s which is close to optimal given the shot noise and the finite survey volume. This measurement is consistent with dynamical estimates of beta derived from both real-space and redshift-space information. The importance of the method presented here is that nonlinear clustering effects are removed to enable linear correlation anisotropy measurements on scales approaching the translinear regime. We discuss implications for analyses of forthcoming optical redshift surveys in which the dispersion is more than a factor of two greater than in the IRAS data.

We present a new method for reconstructing the cosmological density, peculiar velocity and peculiar gravitational potential on large scales from redshift data. We remove the distorting effects of line-of-sight peculiar motions by using the linear theory of gravitational instability, in which the potential is the solution to a linear partial differential equation first derived by Kaiser. We solve this equation by deriving its Green function; the Green functions for the peculiar velocity and the real-space density follow directly. Reconstruction of the cosmic fields is thus reduced to integration over redshift space with the appropriate Green function kernel. Our algorithm has a single input parameter, beta=Omega^0.6/b, where b is the linear bias factor and Omega is the cosmological density parameter. We discuss the virtues of this method for error control, for estimating beta, and for constraining the bias mechanism.