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.

We present the final sample of 98 detached double white dwarf (WD) binaries found in the Extremely Low Mass (ELM) Survey, a spectroscopic survey targeting <0.3 Msun He-core WDs completed in the Sloan Digital Sky Survey footprint. Over the course of the survey we observed ancillary low mass WD candidates like GD278, which we show is a P=0.19 d double WD binary, as well as candidates that turn out to be field blue straggler/subdwarf A-type stars with luminosities too large to be WDs given their Gaia parallaxes. Here, we define a clean sample of ELM WDs that is complete within our target selection and magnitude range 15<g_0<20 mag. The measurements are consistent with 100% of ELM WDs being 0.0089 < P < 1.5 d double WD binaries, 35% of which belong to the Galactic halo. We infer these are mostly He+CO WD binaries given the measurement constraints. The merger rate of the observed He+CO WD binaries exceeds the formation rate of stable mass transfer AM CVn binaries by a factor of 25, and so the majority of He+CO WD binaries must experience unstable mass transfer and merge. The shortest-period systems like J0651+2844 are signature {\it LISA} verification binaries that can be studied with gravitational waves and light.

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.

The Disk Detective citizen science project aims to find new stars with excess 22-$\mu$m emission from circumstellar dust in the AllWISE data release from the Wide-field Infrared Survey Explorer (WISE). We evaluated 261 Disk Detective objects of interest with imaging with the Robo-AO adaptive optics instrument on the 1.5m telescope at Palomar Observatory and with RetroCam on the 2.5m du Pont telescope at Las Campanas Observatory to search for background objects at 0.15''-12'' separations from each target. Our analysis of these data lead us to reject 7% of targets. Combining this result with statistics from our online image classification efforts implies that at most $7.9\% \pm 0.2\%$ of AllWISE-selected infrared excesses are good disk candidates. Applying our false positive rates to other surveys, we find that the infrared excess searches of McDonald et al. (2012), McDonald et al. (2017), and Marton et al. (2016) all have false positive rates $>70\%$. Moreover, we find that all thirteen disk candidates in Theissen & West (2014) with W4 signal-to-noise >3 are false positives. We present 244 disk candidates that have survived vetting by follow-up imaging. Of these, 213 are newly-identified disk systems. Twelve of these are candidate members of comoving pairs based on \textit{Gaia} astrometry, supporting the hypothesis that warm dust is associated with binary systems. We also note the discovery of 22 $\mu$m excess around two known members of the Scorpius-Centaurus association, and identify known disk host WISEA J164540.79-310226.6 as a likely Sco-Cen member. Thirty-one of these disk candidates are closer than $\sim 125$ pc (including 27 debris disks), making them good targets for direct imaging exoplanet searches.

The infrared dust emission from the white dwarf GD 56 is found to rise and fall by 20% peak-to-peak over 11.2 yr, and is consistent with ongoing dust production and depletion. It is hypothesized that the dust is produced via collisions associated with an evolving dust disk, temporarily increasing the emitting surface of warm debris, and is subsequently destroyed or assimilated within a few years. The variations are consistent with debris that does not change temperature, indicating that dust is produced and depleted within a fixed range of orbital radii. Gas produced in collisions may rapidly re-condense onto grains, or may accrete onto the white dwarf surface on viscous timescales that are considerably longer than Poynting-Robertson drag for micron-sized dust. This potential delay in mass accretion rate change is consistent with multi-epoch spectra of the unchanging Ca II and Mg II absorption features in GD 56 over 15 yr, although the sampling is sparse. Overall these results indicate that collisions are likely to be the source of dust and gas, either inferred or observed, orbiting most or all polluted white dwarfs.

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.

We use new Gaia measurements to explore the origin of the highest velocity stars in the Hypervelocity Star Survey. The measurements reveal a clear pattern in the B-type stars. Halo stars dominate the sample at speeds about 100 km/s below Galactic escape velocity. Disk runaway stars have speeds up to 100 km/s above Galactic escape velocity, but most disk runaways are bound. Stars with speeds about 100 km/s above Galactic escape velocity originate from the Galactic center. Two bound stars may also originate from the Galactic center. Future Gaia measurements will enable a large, clean sample of Galactic center ejections for measuring the massive black hole ejection rate of hypervelocity stars, and for constraining the mass distribution of the Milky Way dark matter halo.

We describe the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept that will provide an unprecedented view of the X-ray sky, performing timing and spectroscopy over both a broad energy band (0.2-30 keV) and a wide range of timescales from microseconds to years. STROBE-X comprises two narrow-field instruments and a wide field monitor. The soft or low-energy band (0.2-12 keV) is covered by an array of lightweight optics (3-m focal length) that concentrate incident photons onto small solid-state detectors with CCD-level (85-175 eV) energy resolution, 100 ns time resolution, and low background rates. This technology has been fully developed for NICER and will be scaled up to take advantage of the longer focal length of STROBE-X. The higher-energy band (2-30 keV) is covered by large-area, collimated silicon drift detectors that were developed for the European LOFT mission concept. Each instrument will provide an order of magnitude improvement in effective area over its predecessor (NICER in the soft band and RXTE in the hard band). Finally, STROBE-X offers a sensitive wide-field monitor (WFM), both to act as a trigger for pointed observations of X-ray transients and also to provide high duty-cycle, high time-resolution, and high spectral-resolution monitoring of the variable X-ray sky. The WFM will boast approximately 20 times the sensitivity of the RXTE All-Sky Monitor, enabling multi-wavelength and multi-messenger investigations with a large instantaneous field of view. This mission concept will be presented to the 2020 Decadal Survey for consideration.

We analyze $\sim$32 years of optical spectra and photometry for the symbiotic binary PU Vul. Light curves for the He I $\lambda$4471, He II $\lambda$4686 and H$\beta$ $\lambda$4861 emission lines reveal an illumination effect, where the hot white dwarf ionizes the outflowing wind of the red giant, and evidence for an eccentric orbit with e $\approx$ 0.16. Along with the gradual appearance of high ionization emission from [Fe VII] and O VI, the relative fluxes of these lines suggest an increase in the effective temperature of the hot component, from roughly $10^5$ K on JD 2448000 (1990) to roughly $2\times 10^5$ K on JD 2455000 (2009). During this period, the luminosity of the hot component dropped by a factor of 4$-$6 to a current value of roughly 1000 L$_{\odot}$.

We report on a spectral-timing analysis of the neutron star low-mass X-ray binary Aql~X-1 with the Neutron Star Interior Composition Explorer (NICER) on the International Space Station. Aql~X-1 was observed with NICER during a dim outburst in 2017 July, collecting approximately $50$ ks of good exposure. The spectral and timing properties of the source correspond to that of an (hard) extreme island state in the atoll classification. We find that the fractional amplitude of the low frequency ($<0.3$ Hz) band-limited noise shows a dramatic turnover as a function of energy: it peaks at 0.5 keV with nearly 5\% rms, drops to $12\%$ rms at 2 keV, and rises to $15\%$ rms at 10 keV. Through the analysis of covariance spectra, we demonstrate that band-limited noise exists in both the soft thermal emission and the power-law emission. Additionally, we measure hard time lags, indicating the thermal emission at $0.5$ keV leads the power-law emission at 10 keV on a timescale of $\sim100$ ms at $0.3$ Hz to $\sim10$ ms at $3$ Hz. Our results demonstrate that the thermal emission in the hard state is intrinsically variable, and driving the modulation of the higher energy power-law. Interpreting the thermal spectrum as disk emission, we find our results are consistent with the disk propagation model proposed for accretion onto black holes.

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).

HD 100546 is a young, early-type star and key laboratory for studying gas giant planet formation. GPI data taken in 2015 and reported by Currie et al. (2015) recover the previously-identified protoplanet candidate HD 100546 b and identify a second emission source at ~13--14 au: either a disk hot spot or a second protoplanetary candidate (HD 100546 "c"). In this short research note, we update the status of HD 100546 as observed by the Gemini Planet Imager by rereducing our original data using a different PSF subtraction method (KLIP instead of A-LOCI), rereducing recently public GPI Campaign Team (GPIES) data, and comparing the quality of the two data sets. Our results support the original findings in Currie et al. (2015).

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 examine the possibility that icy super-Earth mass planets, formed over long time scales (0.1--1~Gyr) at large distances ($\sim$ 200--1000~AU) from their host stars, will develop massive H-rich atmospheres. Within the interior of these planets, high pressure converts CH$_4$ into ethane, butane, or diamond and releases H$_2$. Using simplified models which capture the basic physics of the internal structure, we show that the physical properties of the atmosphere depend on the outflux of H$_2$ from the mantle. When this outflux is $\lesssim 10^{10}$ [molec cm$^{-2}$ s$^{-1}$], the outgassed atmosphere has base pressure $\lesssim$ 1 bar. Larger outflows result in a substantial atmosphere where the base pressure may approach $10^3 - 10^4$ bar. For any pressure, the mean density of these planets, 2.4--3 [g cm$^{-3}$], is much larger than the mean density of Uranus and Neptune, 1.3--1.6 [g cm$^{-3}$]. Thus, observations can distinguish between a Planet Nine with a primordial H/He-rich atmosphere accreted from the protosolar nebula and one with an atmosphere outgassed from the core.

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.

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.

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.

The Disk Detective citizen science project aims to find new stars with 22 micron excess emission from circumstellar dust using data from NASA's WISE mission. Initial cuts on the AllWISE catalog provide an input catalog of 277,686 sources. Volunteers then view images of each source online in 10 different bands to identify false-positives (galaxies, background stars, interstellar matter, image artifacts, etc.). Sources that survive this online vetting are followed up with spectroscopy on the FLWO Tillinghast telescope. This approach should allow us to unleash the full potential of WISE for finding new debris disks and protoplanetary disks. We announce a first list of 37 new disk candidates discovered by the project, and we describe our vetting and follow-up process. One of these systems appears to contain the first debris disk discovered around a star with a white dwarf companion: HD 74389. We also report four newly discovered classical Be stars (HD 6612, HD 7406, HD 164137, and HD 218546) and a new detection of 22 micron excess around a previously known debris disk host star, HD 22128.

We analyze a roughly 20 yr set of spectroscopic observations for the symbiotic binary EG And. Radial velocities derived from echelle spectra are best-fit with a circular orbit having orbital period P = 483.3 +- 1.6 days and semi-amplitude K = 7.34 +- 0.07 km/sec. Combined with previous data, these observations rule out an elliptical orbit at the 10-sigma level. Equivalent widths of H~I Balmer emission lines and various absorption features vary in phase with the orbital period. Relative to the radius of the red giant primary, the apparent size of the H II region is consistent with a model where a hot secondary star with effective temperature of roughly 75,000 K ionizes the wind from the red giant.

We estimate the merger rate of double degenerate binaries containing extremely low mass (ELM) <0.3 Msun white dwarfs in the Galaxy. Such white dwarfs are detectable for timescales of 0.1 Gyr -- 1 Gyr in the ELM Survey; the binaries they reside in have gravitational wave merger times of 0.001 Gyr -- 100 Gyr. To explain the observed distribution requires that most ELM white dwarf binary progenitors detach from the common envelope phase with <1 hr orbital periods. We calculate the local space density of ELM white dwarf binaries and estimate a merger rate of 3e-3/yr over the entire disk of the Milky Way; the merger rate in the halo is 10 times smaller. The ELM white dwarf binary merger rate exceeds by a factor of 40 the formation rate of stable mass transfer AM CVn binaries, marginally exceeds the rate of underluminous supernovae, and is identical to the formation rate of R CrB stars. On this basis, we conclude that ELM white dwarf binaries can be the progenitors of all observed AM CVn and possibly underluminous supernovae, however the majority of He+CO white dwarf binaries go through unstable mass transfer and merge, e.g. into single massive ~1 Msun white dwarfs.

We present the discovery of 15 extremely low mass (5 < log{g} < 7) white dwarf candidates, 9 of which are in ultra-compact double-degenerate binaries. Our targeted ELM Survey sample now includes 76 binaries. The sample has a lognormal distribution of orbital periods with a median period of 5.4 hr. The velocity amplitudes imply that the binary companions have a normal distribution of mass with 0.76 Msun mean and 0.25 Msun dispersion. Thus extremely low mass white dwarfs are found in binaries with a typical mass ratio of 1:4. Statistically speaking, 95% of the white dwarf binaries have a total mass below the Chandrasekhar mass and thus are not Type Ia supernova progenitors. Yet half of the observed binaries will merge in less than 6 Gyr due to gravitational wave radiation; probable outcomes include single massive white dwarfs and stable mass transfer AM CVn binaries.

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 report Gemini Planet Imager H band high-contrast imaging/integral field spectroscopy and polarimetry of the HD 100546, a 10 $Myr$-old early-type star recently confirmed to host a thermal infrared bright (super)jovian protoplanet at wide separation, HD 100546 b. We resolve the inner disk cavity in polarized light, recover the thermal-infrared (IR) bright arm, and identify one additional spiral arm. We easily recover HD 100546 b and show that much of its emission originates an unresolved, point source. HD 100546 b likely has extremely red infrared colors compared to field brown dwarfs, qualitatively similar to young cloudy superjovian planets, however, these colors may instead indicate that HD 100546 b is still accreting material from a circumplanetary disk. Additionally, we identify a second point source-like peak at $r_{proj}$ $\sim$ 14 AU, located just interior to or at inner disk wall consistent with being a 10--20 $M_{J}$ candidate second protoplanet-- "HD 100546 c" -- and lying within a weakly polarized region of the disk but along an extension of the thermal IR bright spiral arm. Alternatively, it is equally plausible that this feature is a weakly polarized but locally bright region of the inner disk wall. Astrometric monitoring of this feature over the next 2 years and emission line measurements could confirm its status as a protoplanet, rotating disk hot spot that is possibly a signpost of a protoplanet, or a stationary emission source from within the disk.

We present the discovery of 11 new double degenerate systems containing extremely low-mass white dwarfs (ELM WDs). Our radial velocity observations confirm that all of the targets have orbital periods $\leq$ 1 day. We perform spectroscopic fits and provide a complete set of physical and binary parameters. We review and compare recent evolutionary calculations and estimate that the systematic uncertainty in our mass determinations due to differences in the evolutionary models is small ($\approx$ 0.01 M$_{\odot}$). Five of the new systems will merge due to gravitational wave radiation within a Hubble time, bringing the total number of merger systems found in the ELM Survey to 38. We examine the ensemble properties of the current sample of ELM WD binaries, including the period distribution as a function of effective temperature, and the implications for the future evolution of these systems. We also revisit the empirical boundaries of instability strip of ELM WDs and identify new pulsating ELM WD candidates. Finally, we consider the kinematic properties of our sample of ELM WDs and estimate that a significant fraction of the WDs from the ELM Survey are members of the Galactic halo.

We measure the three components of velocity dispersion, $\sigma_{R},\sigma_{\theta},\sigma_{\phi}$, for stars within 6 < R < 30 kpc of the Milky Way using a new radial velocity sample from the MMT telescope. We combine our measurements with previously published data so that we can more finely sample the stellar halo. We use a maximum likelihood statistical method for estimating mean velocities, dispersions, and covariances assuming only that velocities are normally distributed. The alignment of the velocity ellipsoid is consistent with a spherically symmetric gravitational potential. From the spherical Jeans equation, the mass of the Milky Way is M(<12 kpc) = $1.3\times10^{11}$ M$_{\odot}$ with an uncertainty of 40%. We also find a region of discontinuity, 15 < R < 25 kpc, where the estimated velocity dispersions and anisotropies diverge from their anticipated values, confirming at high significance the break observed by others. We argue that this break in anisotropy is physically explained by coherent stellar velocity structure in the halo, such as the Sgr stream. To significantly improve our understanding of halo kinematics will require combining radial velocities with future Gaia proper motions.

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 describe the discovery of a bright, young Kuiper belt-like debris disk around HD 115600, a $\sim$ 1.4--1.5 M$_\mathrm{\odot}$, $\sim$ 15 Myr old member of the Sco-Cen OB Association. Our H-band coronagraphy/integral field spectroscopy from the \textit{Gemini Planet Imager} shows the ring has a (luminosity scaled) semi major axis of ($\sim$ 22 AU) $\sim$ 48 AU, similar to the current Kuiper belt. The disk appears to have neutral scattering dust, is eccentric (e $\sim$ 0.1--0.2), and could be sculpted by analogues to the outer solar system planets. Spectroscopy of the disk ansae reveal a slightly blue to gray disk color, consistent with major Kuiper belt chemical constituents, where water-ice is a very plausible dominant constituent. Besides being the first object discovered with the next generation of extreme adaptive optics systems (i.e. SCExAO, GPI, SPHERE), HD 115600's debris ring and planetary system provides a key reference point for the early evolution of the solar system, the structure and composition of the Kuiper belt, and the interaction between debris disks and 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.

We measure proper motions with the Hubble Space Telescope for 16 extreme radial velocity stars, mostly unbound B stars in the Milky Way halo. Twelve of these stars have proper motions statistically consistent with zero, and thus have radial trajectories statistically consistent with a Galactic center "hypervelocity star" origin. The trajectory of HE 0437-5439 is consistent with both Milky Way and Large Magellanic Cloud origins. A Galactic center origin is excluded at 3-sigma confidence for two of the lowest radial velocity stars in our sample, however. These two stars are probable disk runaways and provide evidence for ~500 km/s ejections from the disk. We also measure a significant proper motion for the unbound sdO star US 708. Its 1,000 km/s motion is in some tension with proposed supernova ejection models, but can be explained if US 708 was ejected from the stellar halo. In the future, we expect Gaia will better constrain the origin of these remarkable unbound stars.

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.

The high rate of planet detection among solar-type stars argues that planet formation is common. It is also generally assumed that planets form in protoplanetary discs like those observed in nearby star forming regions. On what timescale does the transformation from discs to planets occur? Here we show that current inventories of planets and protoplanetary discs are sensitive enough to place basic constraints on the timescale and efficiency of the planet formation process. A comparison of planet detection statistics and the measured solid reservoirs in T Tauri discs suggests that planet formation is likely already underway at the few Myr age of the discs in Taurus-Auriga, with a large fraction of solids having been converted into large objects with low millimeter opacity and/or sequestered at small disc radii where they are difficult to detect at millimeter wavelengths.

We present new $L^\prime$ (3.8 $\mu m$) and Br-$\alpha$ (4.05 $\mu m$) data and reprocessed archival $L^\prime$ data for the young, planet-hosting star HR 8799 obtained with Keck/NIRC2, VLT/NaCo and Subaru/IRCS. We detect all four HR 8799 planets in each dataset at a moderate to high signal-to-noise (SNR $\gtrsim$ 6-15). We fail to identify a fifth planet, "HR 8799 f", at $r$ $<$ 15 $AU$ at a 5-$\sigma$ confidence level: one suggestive, marginally significant residual at 0.2" is most likely a PSF artifact. Assuming companion ages of 30 $Myr$ and the Baraffe (Spiegel \& Burrows) planet cooling models, we rule out an HR 8799 f with mass of 5 $M_{J}$ (7 $M_{J}$), 7 $M_{J}$ (10 $M_{J}$), and 12 $M_{J}$ (13 $M_{J}$) at $r_{proj}$ $\sim$ 12 $AU$, 9 $AU$, and 5 $AU$, respectively. All four HR 8799 planets have red early T dwarf-like $L^\prime$ - [4.05] colors, suggesting that their SEDs peak in between the $L^\prime$ and $M^\prime$ broadband filters. We find no statistically significant difference in HR 8799 cde's colors. Atmosphere models assuming thick, patchy clouds appear to better match HR 8799 bcde's photometry than models assuming a uniform cloud layer. While non-equilibrium carbon chemistry is required to explain HR 8799 bc's photometry/spectra, evidence for it from HR 8799 de's photometry is weaker. Future, deep IR spectroscopy/spectrophotometry with the Gemini Planet Imager, SCExAO/CHARIS, and other facilities may clarify whether the planets are chemically similar or heterogeneous.

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.

We analyze very deep IRAC and MIPS photometry of $\sim$ 12,500 members of the 14 Myr old Double Cluster, h and $\chi$ Persei, building upon on our earlier, shallower Spitzer studies (Currie et al. 2007a, 2008a). Numerous likely members show infrared (IR) excesses at 8 {\mu}m and 24 $\mu$m indicative of circumstellar dust. The frequency of stars with 8 $\mu$m excess is at least 2% for our entire sample, slightly lower (higher) for B/A stars (later type, lower-mass stars). Optical spectroscopy also identifies gas in about 2% of systems but with no clear trend between the presence of dust and gas. Spectral energy distribution (SED) modeling of 18 sources with detections at optical wavelengths through MIPS 24 $\mu m$ reveals a diverse set of disk evolutionary states, including a high fraction of transitional disks, although similar data for all disk-bearing members would provide better constraints. We combine our results with those for other young clusters to study the global evolution of dust/gas disks. For nominal cluster ages, the e-folding times ($\tau_{o}$) for the frequency of warm dust and gas are 2.75 Myr and 1.75 Myr respectively. Assuming a revised set of ages for some clusters (e.g. Bell et al. 2013), these timescales increase to 5.75 and 3.75 Myr, respectively, implying a significantly longer typical protoplanetary disk lifetime. The transitional disk duration, averaged over multiple evolutionary pathways, is $\sim$ 1 Myr. Finally, 24 $\mu$m excess frequencies for 4-6 M$_{\odot}$ stars appear lower than for 1-2.5 M$_{\odot}$ stars in other 10-30 Myr old clusters.

In most theories of planet formation, the snow-line represents a boundary between the emergence of the interior rocky planets and the exterior ice giants. The wide separation of the snow-line makes the discovery of transiting worlds challenging, yet transits would allow for detailed subsequent characterization. We present the discovery of Kepler-421b, a Uranus-sized exoplanet transiting a G9/K0 dwarf once every 704.2 days in a near-circular orbit. Using public Kepler photometry, we demonstrate that the two observed transits can be uniquely attributed to the 704.2 day period. Detailed light curve analysis with BLENDER validates the planetary nature of Kepler-421b to >4 sigmas confidence. Kepler-421b receives the same insolation as a body at ~2AU in the Solar System and for a Uranian albedo would have an effective temperature of ~180K. Using a time-dependent model for the protoplanetary disk, we estimate that Kepler-421b's present semi-major axis was beyond the snow-line after ~3Myr, indicating that Kepler-421b may have formed at its observed location.

We report the discovery of a detached 20 min orbital period binary white dwarf. WD0931+444 (SDSS J093506.93+441106.9) was previously classified as a WD + M dwarf system based on its optical spectrum. Our time-resolved optical spectroscopy observations obtained at the 8m Gemini and 6.5m MMT reveal peak-to-peak radial velocity variations of 400 km/s every 20 min for the WD, but no velocity variations for the M dwarf. In addition, high-speed photometry from the McDonald 2.1m telescope shows no evidence of variability nor evidence of a reflection effect. An M dwarf companion is physically too large to fit into a 20 min orbit. Thus, the orbital motion of the WD is almost certainly due to an invisible WD companion. The M dwarf must be either an unrelated background object or the tertiary component of a hiearchical triple system. WD0931+444 contains a pair of WDs, a 0.32 Msol primary and a >0.14 Msol secondary, at a separation of >0.19 Rsol. After J0651+2844, WD0931+444 becomes the second-shortest period detached binary WD currently known. The two WDs will lose angular momentum through gravitational wave radiation and merge in <9 Myr. The log h ~ -22 gravitational wave strain from WD0931+444 is strong enough to make it a verification source for gravitational wave missions in the milli-Hertz frequency range, e.g. the evolved Laser Interferometer Space Antenna (eLISA), bringing the total number of known eLISA verification sources to nine.

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 describe our completed spectroscopic survey for unbound hypervelocity stars (HVSs) ejected from the Milky Way. Three new discoveries bring the total number of unbound HVSs to 21. We place new constraints on the nature of HVSs and on their distances using moderate resolution MMT spectroscopy. Half of the HVSs are fast rotators; they are certain 2.5-4 Msun main sequence stars at 50 - 120 kpc distances. Correcting for stellar lifetime, our survey implies that unbound 2.5-4 Msun stars are ejected from the Milky Way at a rate of 1.5e-6 /yr. The observed HVSs are likely ejected continuously over the past 200 Myr and do not share a common flight time. The anisotropic spatial distribution of HVSs on the sky remains puzzling. Southern hemisphere surveys like SkyMapper will soon allow us to map the all-sky distribution of HVSs. Future proper motion measurements with Hubble Space Telescope and Gaia will provide strong constraints on origin. All existing observations are consistent with HVS ejections from encounters with the massive black hole in the Galactic center.

We present the discovery of an unusual, tidally-distorted extremely low mass white dwarf (WD) with nearly solar metallicity. Radial velocity measurements confirm that this is a compact binary with an orbital period of 2.6975 hrs and a velocity semi-amplitude of K = 108.7 km/s. Analysis of the hydrogen Balmer lines yields an effective temperature of Teff = 8380 K and a surface gravity of log g = 6.21 that in turn indicate a mass of M = 0.16 Msol and a cooling age of 4.2 Gyr. In addition, a detailed analysis of the observed metal lines yields abundances of log Mg/H = -3.90, log Ca/H = -5.80, log Ti/H = -6.10, log Cr/H = -5.60, and log Fe/H = -4.50, similar to the sun. We see no evidence of a debris disk from which these metals would be accreted though the possibility cannot entirely be ruled out. Other potential mechanisms to explain the presence of heavy elements are discussed. Finally, we expect this system to ultimately undergo unstable mass transfer and merge to form a ~0.3-0.6 Msol WD in a few Gyr.

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 report the discovery of two new pulsating extremely low-mass (ELM) white dwarfs (WDs), SDSS J161431.28+191219.4 (hereafter J1614) and SDSS J222859.93+362359.6 (hereafter J2228). Both WDs have masses <0.25 Msun and thus likely harbor helium cores. Spectral fits indicate these are the two coolest pulsating WDs ever found. J1614 has Teff = 8880 +/- 170 K and log g = 6.66 +/- 0.14, which corresponds to a roughly 0.19 Msun WD. J2228 is considerably cooler, with a Teff = 7870 +/- 120 K and log g = 6.03 +/- 0.08, which corresponds to a roughly 0.16 Msun WD, making it the coolest and lowest-mass pulsating WD known. There are multiple ELM WDs with effective temperatures between the warmest and coolest known ELM pulsators that do not pulsate to observable amplitudes, which questions the purity of the instability strip for low-mass WDs. In contrast to the CO-core ZZ Ceti stars, which are believed to represent a stage in the evolution of all such WDs, ELM WDs may not all evolve as a simple cooling sequence through an instability strip. Both stars exhibit long-period variability (1184-6235 s) consistent with non-radial g-mode pulsations. Although ELM WDs are preferentially found in close binary systems, both J1614 and J2228 do not exhibit significant radial-velocity variability, and are perhaps in low-inclination systems or have low-mass companions. These are the fourth and fifth pulsating ELM WDs known, all of which have hydrogen-dominated atmospheres, establishing these objects as a new class of pulsating WD.

We present deep Keck/NIRC2 1.6 and 3.8 $\mu m$ imaging of Fomalhaut to constrain the near-infrared (IR) brightness of Fomalhaut b, recently confirmed as a likely planet (Currie et al. 2012a), and search for additional planets at r_proj = 15--150 AU. Using advanced/novel PSF subtraction techniques, we identify seven candidate substellar companions at projected separations comparable to Fomalhaut b. However, multi-epoch data shows them to be background objects. We set a new 3-$\sigma$ upper limit for Fomalhaut b's H-band brightness of m(H) $\sim$ 23.15 or 1.5 to 4.5 M_J, depending on the system age and planet cooling models used. We do not recover the possible point source reported from Spitzer/IRAC data: at its location detection limits are similar to those for Fomalhaut b. Our data when combined with other recent work rule out planets with masses and projected separations comparable to HR 8799 bcde and planets with masses M $>$ 3 M_{J} at r_proj $>$ 45 AU. The James Webb Space Telescope will likely be required to shed substantial further light on Fomalhaut's planetary system.

We analyze Keck ESI spectroscopy of HVS17, a B-type star traveling with a Galactic rest frame radial velocity of +445 km/s in the outer halo of the Milky Way. HVS17 has the projected rotation of a main sequence B star and is chemically peculiar, with solar iron abundance and sub-solar alpha abundance. Comparing measured T_eff and logg with stellar evolution tracks implies that HVS17 is a 3.91 +-0.09 Msun, 153 +-9 Myr old star at a Galactocentric distance of r=48.5 +-4.6 kpc. The time between its formation and ejection significantly exceeds 10 Myr and thus is difficult to reconcile with any Galactic disk runaway scenario involving massive stars. The observations are consistent, on the other hand, with a hypervelocity star ejection from the Galactic center. We show that Gaia proper motion measurements will easily discriminate between a disk and Galactic center origin, thus allowing us to use HVS17 as a test particle to probe the shape of the Milky Way's dark matter halo.

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.

We present optical spectroscopy, astrometry, radio, and X-ray observations of the runaway binary LP 400-22. We refine the orbital parameters of the system based on our new radial velocity observations. Our parallax data indicate that LP 400-22 is significantly more distant (3 sigma lower limit of 840 pc) than initially predicted. LP 400-22 has a tangential velocity in excess of 830 km/s; it is unbound to the Galaxy. Our radio and X-ray observations fail to detect a recycled millisecond pulsar companion, indicating that LP 400-22 is a double white dwarf system. This essentially rules out a supernova runaway ejection mechanism. Based on its orbit, a Galactic center origin is also unlikely. However, its orbit intersects the locations of several globular clusters; dynamical interactions between LP 400-22 and other binary stars or a central black hole in a dense cluster could explain the origin of this unusual binary.

We present the discovery of 17 low mass white dwarfs (WDs) in short-period P<1 day binaries. Our sample includes four objects with remarkable log(g)~5 surface gravities and orbital solutions that require them to be double degenerate binaries. All of the lowest surface gravity WDs have metal lines in their spectra implying long gravitational settling times or on-going accretion. Notably, six of the WDs in our sample have binary merger times <10 Gyr. Four have >=0.9 Msun companions. If the companions are massive WDs, these four binaries will evolve into stable mass transfer AM CVn systems and possibly explode as underluminous supernovae. If the companions are neutron stars, then these may be milli-second pulsar binaries. These discoveries increase the number of detached, double degenerate binaries in the ELM Survey to 54; 31 of these binaries will merge within a Hubble time.

We report the discovery of the second and third pulsating extremely low mass white dwarfs (WDs), SDSS J111215.82+111745.0 (hereafter J1112) and SDSS J151826.68+065813.2 (hereafter J1518). Both have masses < 0.25 Msun and effective temperatures below 10,000 K, establishing these putatively He-core WDs as a cooler class of pulsating hydrogen-atmosphere WDs (DAVs, or ZZ Ceti stars). The short-period pulsations evidenced in the light curve of J1112 may also represent the first observation of acoustic (p-mode) pulsations in any WD, which provide an exciting opportunity to probe this WD in a complimentary way compared to the long-period g-modes also present. J1112 is a Teff = 9590 +/- 140 K and log(g) = 6.36 +/- 0.06 WD. The star displays sinusoidal variability at five distinct periodicities between 1792-2855 s. In this star we also see short-period variability, strongest at 134.3 s, well short of expected g-modes for such a low-mass WD. The other new pulsating WD, J1518, is a Teff = 9900 +/- 140 K and log(g) = 6.80 +/- 0.05 WD. The light curve of J1518 is highly non-sinusoidal, with at least seven significant periods between 1335-3848 s. Consistent with the expectation that ELM WDs must be formed in binaries, these two new pulsating He-core WDs, in addition to the prototype SDSS J184037.78+642312.3, have close companions. However, the observed variability is inconsistent with tidally induced pulsations and is so far best explained by the same hydrogen partial-ionization driving mechanism at work in classic C/O-core ZZ Ceti stars.

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.

Mass transfer from an evolved donor star to its binary companion is a standard feature of stellar evolution in binaries. In wide binaries, the companion star captures some of the mass ejected in a wind by the primary star. The captured material forms an accretion disk. Here, we study the evolution of wind-accretion disks, using a numerical approach which allows us to follow the long term evolution. For a broad range of initial conditions, we derive the radial density and temperature profiles of the disk. In most cases, wind-accretion leads to long-lived stable disks over the lifetime of the AGB donor star. The disks have masses of a few times 10^{-5}-10^{-3} M_sun, with surface density and temperature profiles that follow broken power-laws. The total mass in the disk scales approximately linearly with the viscosity parameter used. Roughly 50% to 80% of the mass falling into the disk accretes onto the central star; the rest flows out through the outer edge of the disk into the stellar wind of the primary. For systems with large accretion rates, the secondary accretes as much as 0.1 M_sun. When the secondary is a white dwarf, accretion naturally leads to nova and supernova eruptions. For all types of secondary star, the surface density and temperature profiles of massive disks resemble structures observed in protoplanetary disks, suggesting that coordinated observational programs might improve our understanding of uncertain disk physics.

We present Subaru/IRCS J band data for Fomalhaut and a (re)reduction of archival 2004--2006 HST/ACS data first presented by Kalas et al. (2008). We confirm the existence of a candidate exoplanet, Fomalhaut b, in both the 2004 and 2006 F606W data sets at a high signal-to-noise. Additionally, we confirm the detection at F814W and present a new detection in F435W. Fomalhaut b's space motion may be consistent with it being in an apsidally-aligned, non debris ring-crossing orbit, although new astrometry is required for firmer conclusions. We cannot confirm that Fomalhaut b exhibits 0.7-0.8 mag variability cited as evidence for planet accretion or a semi-transient dust cloud. The new, combined optical SED and IR upper limits confirm that emission identifying Fomalhaut b originates from starlight scattered by small dust, but this dust is most likely associated with a massive body. The Subaru and IRAC/4.5 micron upper limits imply M < 2 Mj, still consistent with the range of Fomalhaut b masses needed to sculpt the disk. Fomalhaut b is very plausibly "a planet identified from direct imaging" even if current images of it do not, strictly speaking, show thermal emission from a directly imaged planet.

We report the detection of orbital decay in the 12.75-min, detached binary white dwarf (WD) SDSS J065133.338+284423.37 (hereafter J0651). Our photometric observations over a 13-month baseline constrain the orbital period to 765.206543(55) s and indicate the orbit is decreasing as a rate of (-9.8 +/- 2.8) x 10^(-12) s/s (or -0.31 +/- 0.09 ms/yr). We revise the system parameters based on our new photometric and spectroscopic observations: J0651 contains two WDs with M1 = 0.26 +/- 0.04 Msun and M2 = 0.50 +/- 0.04 Msun. General relativity predicts orbital decay due to gravitational wave radiation of (-8.2 +/- 1.7) x 10^(-12) s/s (or -0.26 +/- 0.05 ms/yr). Our observed rate of orbital decay is consistent with this expectation. J0651 is currently the second-loudest gravitational wave source known in the milli-Hertz range and the loudest non-interacting binary, which makes it an excellent verification source for future missions aimed at directly detecting gravitational waves. Our work establishes the feasibility of monitoring this system's orbital period decay at optical wavelengths.

We obtain Keck HIRES spectroscopy of HVS5, one of the fastest unbound stars in the Milky Way halo. We show that HVS5 is a 3.62 +- 0.11 Msun main sequence B star at a distance of 50 +- 5 kpc. The difference between its age and its flight time from the Galactic center is 105 +-18(stat)+-30(sys) Myr; flight times from locations elsewhere in the Galactic disk are similar. This 10^8 yr `arrival time' between formation and ejection is difficult to reconcile with any ejection scenario involving massive stars that live for only 10^7 yr. For comparison, we derive arrival times of 10^7 yr for two unbound runaway B stars, consistent with their disk origin where ejection results from a supernova in a binary system or dynamical interactions between massive stars in a dense star cluster. For HVS5, ejection during the first 10^7 yr of its lifetime is ruled out at the 3-sigma level. Together with the 10^8 yr arrival times inferred for three other well-studied hypervelocity stars (HVSs), these results are consistent with a Galactic center origin for the HVSs. If the HVSs were indeed ejected by the central black hole, then the Galactic center was forming stars ~200 Myr ago, and the progenitors of the HVSs took ~100 Myr to enter the black hole's loss cone.

This pedagogical chapter covers the theory of planet formation, with an emphasis on the physical processes relevant to current research. After summarizing empirical constraints from astronomical and geophysical data, we describe the structure and evolution of protoplanetary disks. We consider the growth of planetesimals and of larger solid protoplanets, followed by the accretion of planetary atmospheres, including the core accretion instability. We also examine the possibility that gas disks fragment directly into giant planets and/or brown dwarfs. We defer a detailed description of planet migration and dynamical evolution to other work, such as the complementary chapter in this series by Morbidelli (available at arXiv:1106.4114).

The Pluto system provides a unique local laboratory for the study of binaries with multiple low mass companions. In this paper, we study the orbital stability of P4, the most recently discovered moon in the Pluto system. This newfound companion orbits near the plane of the Pluto-Charon binary, roughly halfway between the two minor moons Nix and Hydra. We use a suite of few body integrations to constrain the masses of Nix and Hydra, and the orbital parameters of P4. For the system to remain stable over the age of the Solar System, the masses of Nix and Hydra likely do not exceed 5e16 kg and 9e16 kg, respectively. These upper limits assume a fixed mass ratio between Nix and Hydra at the value implied by their median optical brightness. Our study finds that stability is more sensitive to their total mass and that a downward revision of Charon's eccentricity (from our adopted value of 0.0035) is unlikely to significantly affect our conclusions. Our upper limits are an order of magnitude below existing astrometric limits on the masses of Nix and Hydra. For a density at least that of ice, the albedos of Nix and Hydra would exceed 0.3. This constraint implies they are icy, as predicted by giant impact models. Even with these low masses, P4 only remains stable if its eccentricity e < 0.02. The 5:1 commensurability with Charon is particularly unstable, Combining stability constraints with the observed mean motion places the preferred orbit for P4 just exterior to the 5:1 resonance. These predictions will be tested when the New Horizons satellite visits Pluto. Based on the results for the Pluto-Charon system, we expect that circumbinary, multi-planet systems will be more widely spaced than their singleton counterparts. Further, circumbinary exoplanets close to the three-body stability boundary, such as those found by Kepler, are less likely to have other companions nearby.

We report the discovery of the first pulsating extremely low mass (ELM) white dwarf (WD), SDSS J184037.78+642312.3 (hereafter J1840). This DA (hydrogen-atmosphere) WD is by far the coolest and the lowest-mass pulsating WD, with Teff = 9100 \pm 170 K and log g = 6.22 \pm 0.06, which corresponds to a mass ~ 0.17 Msun. This low-mass pulsating WD greatly extends the DAV (or ZZ Ceti) instability strip, effectively bridging the log g gap between WDs and main sequence stars. We detect high-amplitude variability in J1840 on timescales exceeding 4000 s, with a non-sinusoidal pulse shape. Our observations also suggest that the variability is multi-periodic. The star is in a 4.6 hr binary with another compact object, most likely another WD. Future, more extensive time-series photometry of this ELM WD offers the first opportunity to probe the interior of a low-mass, presumably He-core WD using the tools of asteroseismology.

We present new radial velocity and X-ray observations of extremely low-mass (ELM, 0.2 Msol) white dwarf candidates in the Sloan Digital Sky Survey (SDSS) Data Release 7 area. We identify seven new binary systems with 1-18 h orbital periods. Five of the systems will merge due to gravitational wave radiation within 10 Gyr, bringing the total number of merger systems found in the ELM Survey to 24. The ELM Survey has now quintupled the known merger white dwarf population. It has also discovered the eight shortest period detached binary white dwarf systems currently known. We discuss the characteristics of the merger and non-merger systems observed in the ELM Survey, including their future evolution. About half of the systems have extreme mass ratios. These are the progenitors of the AM Canum Venaticorum systems and supernovae .Ia. The remaining targets will lead to the formation of extreme helium stars, subdwarfs, or massive white dwarfs. We identify three targets that are excellent gravitational wave sources. These should be detected by the Laser Interferometer Space Antenna (LISA)-like missions within the first year of operation. The remaining targets are important indicators of what the Galactic foreground may look like for gravitational wave observatories.

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 present the discovery of five new unbound hypervelocity stars (HVSs) in the outer Milky Way halo. Using a conservative estimate of Galactic escape velocity, our targeted spectroscopic survey has now identified 16 unbound HVSs as well as a comparable number of HVSs ejected on bound trajectories. A Galactic center origin for the HVSs is supported by their unbound velocities, the observed number of unbound stars, their stellar nature, their ejection time distribution, and their Galactic latitude and longitude distribution. Other proposed origins for the unbound HVSs, such as runaway ejections from the disk or dwarf galaxy tidal debris, cannot be reconciled with the observations. An intriguing result is the spatial anisotropy of HVSs on the sky, which possibly reflects an anisotropic potential in the central 10-100 pc region of the Galaxy. Further progress requires measurement of the spatial distribution of HVSs over the southern sky. Our survey also identifies seven B supergiants associated with known star-forming galaxies; the absence of B supergiants elsewhere in the survey implies there are no new star-forming galaxies in our survey footprint to a depth of 1-2 Mpc.

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.

We develop statistical methods for identifying star streams in the halo of the Milky Way galaxy that exploit observed spatial and radial velocity distributions. Within a great circle, departures of the observed spatial distribution from random provide a measure of the likelihood of a potential star stream. Comparisons between the radial velocity distribution within a great circle and the radial velocity distribution of the entire sample also measure the statistical significance of potential streams. The radial velocities enable construction of a more powerful joint statistical test for identifying star streams in the Milky Way halo. Applying our method to halo stars in the Hypervelocity Star (HVS) survey, we detect the Sagittarius stream at high significance. Great circle counts and comparisons with theoretical models suggest that the Sagittarius stream comprises 10% to 17% of the halo stars in the HVS sample. The population of blue stragglers and blue horizontal branch stars varies along the stream and is a potential probe of the distribution of stellar populations in the Sagittarius dwarf galaxy prior to disruption.

Extremely low mass (ELM) white dwarfs (WDs) with masses <0.25 Msun are rare objects that result from compact binary evolution. Here, we present a targeted spectroscopic survey of ELM WD candidates selected by color. The survey is 71% complete and has uncovered 18 new ELM WDs. Of the 7 ELM WDs with follow-up observations, 6 are short-period binaries and 4 have merger times less than 5 Gyr. The most intriguing object, J1741+6526, likely has either a pulsar companion or a massive WD companion making the system a possible supernova Type Ia or .Ia progenitor. The overall ELM Survey has now identified 19 double degenerate binaries with <10 Gyr merger times. The significant absence of short orbital period ELM WDs at cool temperatures suggests that common envelope evolution creates ELM WDs directly in short period systems. At least one-third of the merging systems are halo objects, thus ELM WD binaries continue to form and merge in both the disk and the halo.

We report the discovery of a new detached, double white dwarf system with an orbital period of 39.8 min. We targeted SDSS J163030.58+423305.8 (hereafter J1630) as part of our radial velocity program to search for companions around low-mass white dwarfs using the 6.5m MMT. We detect peak-to-peak radial velocity variations of 576 km/s. The mass function and optical photometry rule out main-sequence companions. In addition, no milli-second pulsar companions are detected in radio observations. Thus the invisible companion is most likely another white dwarf. Unlike the other 39 min binary SDSS J010657.39-100003.3, follow-up high speed photometric observations of J1630 obtained at the McDonald 2.1m telescope do not show significant ellipsoidal variations, indicating a higher primary mass and smaller radius. The absence of eclipses constrain the inclination angle to <82deg. J1630 contains a pair of white dwarfs, 0.3 Msun primary + >0.3 Msun invisible secondary, at a separation of >0.32 Rsun. The two white dwarfs will merge in less than 31 Myr. Depending on the core composition of the companion, the merger will form either a single core-He burning subdwarf star or a rapidly rotating massive white dwarf. The gravitational wave strain from J1630 is detectable by instruments like the Laser Interferometer Space Antenna (LISA) within the first year of operation.

We have discovered a detached pair of white dwarfs (WDs) with a 12.75 min orbital period and a 1,315 km/s radial velocity amplitude. We measure the full orbital parameters of the system using its light curve, which shows ellipsoidal variations, Doppler boosting, and primary and secondary eclipses. The primary is a 0.25 Msun tidally distorted helium WD, only the second tidally distorted WD known. The unseen secondary is a 0.55 Msun carbon-oxygen WD. The two WDs will come into contact in 0.9 Myr due to loss of energy and angular momentum via gravitational wave radiation. Upon contact the systems may merge yielding a rapidly spinning massive WD, form a stable interacting binary, or possibly explode as an underluminous supernova type Ia. The system currently has a gravitational wave strain of 10^-22, about 10,000 times larger than the Hulse-Taylor pulsar; this system would be detected by the proposed LISA gravitational wave mission in the first week of operation. This system's rapid change in orbital period will provide a fundamental test of general relativity.

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 identify SDSS J010657.39-100003.3 (hereafter J0106-1000) as the shortest period detached binary white dwarf (WD) system currently known. We targeted J0106-1000 as part of our radial velocity program to search for companions around known extremely low-mass (ELM, ~ 0.2 Msol) WDs using the 6.5m MMT. We detect peak-to-peak radial velocity variations of 740 km/s with an orbital period of 39.1 min. The mass function and optical photometry rule out a main-sequence star companion. Follow-up high-speed photometric observations obtained at the McDonald 2.1m telescope reveal ellipsoidal variations from the distorted primary but no eclipses. This is the first example of a tidally distorted WD. Modeling the lightcurve, we constrain the inclination angle of the system to be 67 +- 13 deg. J0106-1000 contains a pair of WDs (0.17 Msol primary + 0.43 Msol invisible secondary) at a separation of 0.32 Rsol. The two WDs will merge in 37 Myr and most likely form a core He-burning single subdwarf star. J0106-1000 is the shortest timescale merger system currently known. The gravitational wave strain from J0106-1000 is at the detection limit of the Laser Interferometer Space Antenna (LISA). However, accurate ephemeris and orbital period measurements may enable LISA to detect J0106-1000 above the Galactic background noise.

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 spectroscopic observations of 21 low-mass (<0.45 M_sun) white dwarfs (WDs) from the Palomar-Green Survey obtained over four years. We use both radial velocities and infrared photometry to identify binary systems, and find that the fraction of single, low-mass WDs is <30%. We discuss the potential formation channels for these single stars including binary mergers of lower-mass objects. However, binary mergers are not likely to explain the observed number of single low-mass WDs. Thus additional formation channels, such as enhanced mass loss due to winds or interactions with substellar companions, are likely.

We describe new radial velocity and X-ray observations of extremely low-mass white dwarfs (ELM WDs, ~0.2 Msol) in the Sloan Digital Sky Survey Data Release 4 and the MMT Hypervelocity Star survey. We identify four new short period binaries, including two merger systems. These observations bring the total number of short period binary systems identified in our survey to 20. No main-sequence or neutron star companions are visible in the available optical photometry, radio, and X-ray data. Thus, the companions are most likely WDs. Twelve of these systems will merge within a Hubble time due to gravitational wave radiation. We have now tripled the number of known merging WD systems. We discuss the characteristics of this merger sample and potential links to underluminous supernovae, extreme helium stars, AM CVn systems, and other merger products. We provide new observational tests of the WD mass-period distribution and cooling models for ELM WDs. We also find evidence for a new formation channel for single low-mass WDs through binary mergers of two lower mass objects.

We analyze radial velocity observations of the 12 extremely low-mass <0.25 Msol white dwarfs (WDs) in the MMT Hypervelocity Star Survey. Eleven of the 12 WDs are binaries with orbital periods shorter than 14 hours; the one non-variable WD is possibly a pole-on system among our non-kinematically selected targets. Our sample is unique: it is complete in a well-defined range of apparent magnitude and color. The orbital mass functions imply that the unseen companions are most likely other WDs, although neutron star companions cannot be excluded. Six of the 11 systems with orbital solutions will merge within a Hubble time due to the loss of angular momentum through gravitational wave radiation. The quickest merger is J0923+3028, a g=15.7 ELM WD binary with a 1.08 hr orbital period and a <130 Myr merger time. The chance of a supernova Ia event among our ELM WDs is only 1%-7%, however. Three binary systems (J0755+4906, J1233+1602, and J2119-0018) have extreme mass ratios and will most likely form stable mass-transfer AM CVn systems. Two of these objects, SDSS J1233+1602 and J2119-0018, are the lowest surface gravity WDs ever found; both show Ca II absorption likely from accretion of circumbinary material. We predict that at least one of our WDs is an eclipsing detached double WD system, important for constraining helium core WD models.

We study a complete, colour-selected sample of double-degenerate binary systems containing extremely low mass (ELM) <0.25 Msol white dwarfs (WDs). We show, for the first time, that Milky Way disk ELM WDs have a merger rate of approximately 4 x 10^(-5)/yr due to gravitational wave radiation. The merger end-product depends on the mass ratio of the binary. The ELM WD systems that undergo stable mass transfer can account for >3% of AM CVn stars. More importantly, the ELM WD systems that may detonate merge at a rate comparable to the estimated rate of underluminous SNe, rare explosions estimated to produce only ~0.2 Msol worth of ejecta. At least 25% of our ELM WD sample belong to the old thick disk and halo components of the Milky Way. Thus, if merging ELM WD systems are the progenitors of underluminous SNe, transient surveys must find them in both elliptical and spiral galaxies.

We measure the radial velocity curve of the eclipsing detached white dwarf binary NLTT 11748. The primary exhibits velocity variations with a semi-amplitude of 273 km/s and an orbital period of 5.641 hr. We do not detect any spectral features from the secondary star, or any spectral changes during the secondary eclipse. We use our composite spectrum to constrain the temperature and surface gravity of the primary to be Teff = 8690 (140) K and log g = 6.54 (0.05), which correspond to a mass of 0.18 Msol. For an inclination angle of 89.9 deg. derived from the eclipse modeling, the mass function requires a 0.76 Msol companion. The merger time for the system is 7.2 Gyr. However, due to the extreme mass ratio of 0.24, the binary will most likely create an AM CVn system instead of a merger.

The Hypervelocity Star survey presents the currently largest sample of radial velocity measurements of halo stars out to 80 kpc. We apply spherical Jeans modeling to these data in order to derive the mass profile of the Galaxy. We restrict the analysis to distances larger than 25 kpc from the Galactic center, where the density profile of halo stars is well approximated by a single power law with logarithmic slope between -3.5 and -4.5. With this restriction, we also avoid the complication of modeling a flattened Galactic disk. In the range 25 < r < 80 kpc, the radial velocity dispersion declines remarkably little; a robust measure of its logarithmic slope is between -0.05 and -0.1. The circular velocity profile also declines remarkably little with radius. The allowed range of V_c(80 kpc) lies between 175 and 231 km/s, with the most likely value 193 km/s. Compared with the value at the solar location, the Galactic circular velocity declines by less than 20% over an order of magnitude in radius. Such a flat profile requires a massive and extended dark matter halo. The mass enclosed within 80 kpc is 6.9(+3.0-1.2) 10^11 Msun. Our sample of radial velocities is large enough that the biggest uncertainty in the mass is not statistical but systematic, dominated by the density slope and anisotropy of the tracer population. Further progress requires modeling observed datasets within realistic simulations of galaxy formation.

We use Hubble Space Telescope imaging to measure the absolute proper motion of the hypervelocity star (HVS) HE 0437-5439, a short-lived B star located in the direction of the Large Magellanic Cloud (LMC). We observe (\mu_\alpha, \mu_\delta)=(+0.53+-0.25(stat)+-0.33(sys), +0.09+-0.21(stat)+-0.48(sys)) mas/yr. The velocity vector points directly away from the center of the Milky Way; an origin from the center of the LMC is ruled out at the 3-sigma level. The flight time of the HVS from the Milky Way exceeds its main-sequence lifetime, thus its stellar nature requires it to be a blue straggler. The large space velocity rules out a Galactic-disk ejection. Combining the HVS's observed trajectory, stellar nature, and required initial velocity, we conclude that HE 0437-5439 was most likely a compact binary ejected by the Milky Way's central black hole.

Turbulent, two-dimensional, hydrodynamic flows are characterized by the emergence of coherent, long-lived vortices without a need to invoke special initial conditions. Vortices have the ability to sequester particles, with typical radii from about 1 mm to 10 cm, that are slightly decoupled from the gas. A generic feature of discs with surface density and effective temperature profiles that are decreasing, power-law functions of radial distance is that four vortex zones exist for a fixed particle size. In particular, two of the zones form an annulus at intermediate radial distances within which small particles reside. Particle capture by vortices occurs on a dynamical time scale near and at the boundaries of this annulus. As the disc ages and the particles grow via coagulation, the size of the annulus shrinks. Older discs prefer to capture smaller particles because the gas surface density decreases with time, a phenomenon we term "vortex aging". More viscous, more dust-opaque and/or less massive discs can have vortices that age faster and trap a broader range of particle sizes throughout the lifetime of the disc. Thus, how efficiently a disc retains its mass in solids depends on the relative time scales between coagulation and vortex aging. If vortices form in protoplanetary discs, they are important in discs with typical masses and for particles that are likely to condense out of the protostellar nebula. Particle capture also occurs at distances relevant to planet formation. Future infrared, submillimetre and centimetre observations of grain opacity as a function of radial distance will test the hypothesis that vortices serve as nurseries for particle growth in protoplanetary discs.

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.

We present radial velocity observations of four extremely low-mass (0.2 Msol) white dwarfs. All four stars show peak-to-peak radial velocity variations of 540 - 710 km/s with 1.0 - 5.9 hr periods. The optical photometry rules out main-sequence companions. In addition, no milli-second pulsar companions are detected in radio observations. Thus the invisible companions are most likely white dwarfs. Two of the systems are the shortest period binary white dwarfs yet discovered. Due to the loss of angular momentum through gravitational radiation, three of the systems will merge within 500 Myr. The remaining system will merge within a Hubble time. The mass functions for three of the systems imply companions more massive than 0.46 Msol; thus those are carbon/oxygen core white dwarfs. The unknown inclination angles prohibit a definitive conclusion about the future of these systems. However, the chance of a supernova Ia event is only 1% to 5%. These systems are likely to form single R Coronae Borealis stars, providing evidence for a white dwarf + white dwarf merger mechanism for these unusual objects. One of the systems, SDSS J105353.89+520031.0 has a 70% chance of having a low-mass white dwarf companion. This system will probably form a single helium-enriched subdwarf O star. All four white dwarf systems have unusal mass ratios of < 0.2-0.8 that may also lead to the formation of AM CVn systems.

(Abridged) From photometric observations of $\sim$ 47,000 stars and spectroscopy of $\sim$ 11,000 stars, we describe the first extensive study of the stellar population of the famous Double Cluster, h and $\chi$ Persei, down to subsolar masses. Both clusters have E(B-V) $\sim$ 0.52--0.55 and dM = 11.8--11.85; the halo population, while more poorly constrained, likely has identical properties. As determined from the main sequence turnoff, the luminosity of M supergiants, and pre-main sequence isochrones, ages for h Persei, $\chi$ Persei and the halo population all converge on $\approx$ 14 Myr. From these data, we establish the first spectroscopic and photometric membership lists of cluster stars down to early/mid M dwarfs. At minimum, there are $\sim$ 5,000 members within 10' of the cluster centers, while the entire h and $\chi$ Persei region has at least $\sim$ 13,000 and as many as 20,000 members. The Double Cluster contains $\approx$ 8,400 M$_{\odot}$ of stars within 10' of the cluster centers. We estimate a total mass of at least 20,000 M$_{\odot}$. We conclude our study by outlining outstanding questions regarding the properties of h and $\chi$ Persei. From comparing recent work, we compile a list of intrinsic colors and derive a new effective temperature scale for O--M dwarfs, giants, and supergiants.

We report on a single star, B030D, observed as part of a large survey of objects in M31, which has the unusual radial velocity of -780 km/s. Based on details of its spectrum, we find that the star is an F supergiant, with a circumstellar shell. The evolutionary status of the star could be one of a post-mainsequence close binary, a symbiotic nova, or less likely, a post-AGB star, which additional observations could help sort out. Membership of the star in the Andromeda Giant Stream can explain its highly negative velocity.

We present a spectroscopic sample of 910 distant halo stars from the Hypervelocity Star survey from which we derive the velocity dispersion profile of the Milky Way halo. The sample is a mix of 74% evolved horizontal branch stars and 26% blue stragglers. We estimate distances to the stars using observed colors, metallicities, and stellar evolution tracks. Our sample contains twice as many objects with R>50 kpc as previous surveys. We compute the velocity dispersion profile in two ways: with a parametric method based on a Milky Way potential model, and with a non-parametric method based on the caustic technique originally developed to measure galaxy cluster mass profiles. The resulting velocity dispersion profiles are remarkably consistent with those found by two independent surveys based on other stellar populations: the Milky Way halo exhibits a mean decline in radial velocity dispersion of -0.38+-0.12 km/s/kpc over 15<R<75 kpc. This measurement is a useful basis for calculating the total mass and mass distribution of the Milky Way halo.

We present results of optical spectroscopic and BVR_CI_C photometric observations of 77 pre-main sequence (PMS) stars in the Cepheus flare region. A total of 64 of these are newly confirmed PMS stars, originally selected from various published candidate lists. We estimate effective temperatures and luminosities for the PMS stars, and comparing the results with pre-main sequence evolutionary models we estimate stellar masses of 0.2-2.4M_sun and stellar ages of 0.1-15 Myr. Among the PMS stars, we identify 15 visual binaries with separations of 2-10 arcsec. From archival IRAS, 2MASS, and Spitzer data, we construct their spectral energy distributions and classify 5% of the stars as Class I, 10% as Flat SED, 60% as Class II, and 3% as Class III young stellar objects (YSOs). We identify 12 CTTS and 2 WTTS as members of NGC 7023, with mean age of 1.6 Myr. The 13 PMS stars associated with L1228 belong to three small aggregates: RNO 129, L1228A, and L1228S. The age distribution of the 17 PMS stars associated with L1251 suggests that star formation has propagated with the expansion of the Cepheus flare shell. We detect sparse aggregates of 6-7 Myr old PMS stars around the dark clouds L1177 and L1219, at a distance of 400 pc. Three T Tauri stars appear to be associated with the Herbig Ae star SV Cep at a distance of 600 pc. Our results confirm that the molecular complex in the Cepheus flare region contains clouds of various distances and star forming histories.

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 new, deep MIPS photometry and new high signal-to-noise optical spectroscopy of the 2.5 Myr-old IC 348 Nebula. To probe the properties of the IC 348 disk population, we combine these data with previous optical/infrared photometry and spectroscopy to identify stars with gas accretion, to examine their mid-IR colors,and to model their spectral energy distributions. IC 348 contains many sources in different evolutionary states, including protostars and stars surrounded by primordial disks, two kinds of transitional disks, and debris disks. Most disks surrounding eary/intermediate spectral-type stars (> 1.4 Msun at 2.5 Myr) are debris disks; most disks surrounding solar and subsolar-mass stars are primordial disks. At the 1--2 sigma level, more massive stars also have a smaller frequency of gas accretion and smaller mid-IR luminosities than lower-mass stars. These trends are suggestive of a stellar mass-dependent evolution of disks, where most disks around high/intermediate-mass stars shed their primordial disks on rapid, 2.5 Myr timescales. The frequency of MIPS-detected transitional disks is ~ 15--35% for stars plausibly more massive than 0.5 Msun. The relative frequency of transitional disks in IC 348 compared to that for 1 Myr-old Taurus and 5 Myr-old NGC 2362 is consistent with a transition timescale that is a significant fraction of the total primordial disk lifetime.

(Abridged) We describe Spitzer IRAC and MIPS observations of the populous, 5 Myr-old open cluster NGC 2362. Early/intermediate-type confirmed/candidate cluster members either have photospheric mid-IR emission or weak, optically-thin infrared excess emission at < 24 microns consistent with debris disks. Few late-type, solar/subsolar-mass stars have primordial disks. The disk population around late-type stars is dominated by disks with inner holes (canonical 'transition disks') and 'homologously depleted' disks. Both types of disks represent an intermediate stage between primordial disks and debris disks. Thus, we find that multiple paths for the primordial-to-debris disk transition exist. Our results undermine standard arguments in favor of a ~ 0.01 Myr year timescale for the transition based on data from Taurus-Auriga and rule out standard UV photoevaporation scenarios as the primary mechanism to explain the transition. Combining our data with other Spitzer surveys, we investigate the evolution of debris disks around high/intermediate-mass stars and investigate timescales for giant planet formation. If the gas and dust in disks evolve on similar timescales, the formation timescale for gas giant planets surrounding early-type, high/intermediate-mass stars is likely 1--5 Myr. Most solar/subsolar-mass stars detected by Spitzer have SEDs that indicate their disks may be actively leaving the primordial disk phase. Thus, gas giant planet formation may also occur by 5 Myr around solar/subsolar-mass stars as well.

We report the detection of a radial velocity companion to the extremely low mass white dwarf LP400-22. The radial velocity of the white dwarf shows variations with a semi-amplitude of 119 km/s and a 0.98776 day period, which implies a companion mass of M > 0.37 Msun. The optical photometry rules out a main sequence companion. Thus the invisible companion is another white dwarf or a neutron star. Using proper motion measurements and the radial velocity of the binary system, we find that it has an unusual Galactic orbit. LP400-22 is moving away from the Galactic center with a velocity of 396 km/s, which is very difficult to explain by supernova runaway ejection mechanisms. Dynamical interactions with a massive black hole like that in the Galactic center can in principle explain its peculiar velocity, if the progenitor was a triple star system comprised of a close binary and a distant tertiary companion. Until better proper motions become available, we consider LP400-22 to be most likely a halo star with a very unusual orbit.

The detection of lower mass planets now being reported via radial velocity and microlensing surveys suggests that they may be ubiquitous. If missions such as Kepler are able to confirm this, the detection and study of rocky planets via direct imaging with ground-based telescopes of apertures > 20 m is feasible in the thermal infrared. We discuss two cases for detecting rocky planets, the first via detection of molten Earths formed though an Earth-Moon like impact event, and the second via detection of planets around very nearby stars. These observations have the potential to give us a first look at a rocky planet similar to the Earth.

We use published optical spectral and IR excess data from nine young clusters and associations to study the stellar mass dependent dispersal of circumstellar disks. All clusters older than ~3 Myr show a decrease in disk fraction with increasing stellar mass for Solar to higher mass stars. This result is significant at about the 1 $\sigma$ level in each cluster. For the complete set of clusters we reject the null hypothesis--that Solar and intermediate-mass stars lose their disks at the same rate--with 95-99.9% confidence. To interpret this behaviour, we investigate the impact of grain growth, binary companions, and photoevaporation on the evolution of disk signatures. Changes in grain growth timescales at fixed disk temperature may explain why early-type stars with IR excesses appear to evolve faster than their later-type counterparts. Little evidence that binary companions affect disk evolution suggests that photoevaporation is the more likely mechanism for disk dispersal. A simple photoevaporation model provides a good fit to the observed disk fractions for Solar and intermediate-mass stars. Although the current mass-dependent disk dispersal signal is not strong, larger and more complete samples of clusters with ages of 3-5 Myr can improve the significance and provide better tests of theoretical models. In addition, the orbits of extra-Solar planets can constrain models of disk dispersal and migration. We suggest that the signature of stellar mass dependent disk dispersal due to photoevaporation may be present in the orbits of observed extra Solar planets. Planets orbiting hosts more massive than ~1.6 Msun may have larger orbits because the disks in which they formed were dispersed before they could migrate.

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 describe Chandra/ACIS-I observations of the massive ~ 13--14 Myr-old cluster, h Persei, part of the famous Double Cluster (h and chi Persei) in Perseus. Combining the list of Chandra-detected sources with new optical/IR photometry and optical spectroscopy reveals ~ 165 X-ray bright stars with V < 23. Roughly 142 have optical magnitudes and colors consistent with cluster membership. The observed distribution of Lx peaks at Lx ~ 10^30.3 ergs/s and likely traces the bright edge of a far larger population of ~ 0.4--2 Msun X-ray active stars. From a short list of X-ray active stars with IRAC 8 micron excess from warm, terrestrial-zone dust, we derive a maximum X-ray flux incident on forming terrestrial planets. Although there is no correlation between X-ray activity and IRAC excess, the fractional X-ray luminosity correlates with optical colors and spectral type. By comparing the distribution of Lx/L* vs. spectral type and V-I in h Per with results for other 1--100 Myr-old clusters, we show that stars slightly more massive than the Sun (> 1.5 Msun) fall out of X-ray saturation by ~ 10--15 Myr. Changes in stellar structure for > 1.5 Msun stars likely play an important role in this decline of X-ray emission.

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 review the history and structure of star formation in the Taurus-Auriga dark clouds. Our discussion includes a summary of the macroscopic cloud properties, the population of single and binary pre-main sequence stars, the properties of jets and outflows, and detailed summaries of selected individual objects. We include comprehensive tables of dark clouds, young stars, and jets in the clouds.

We present results of optical spectroscopic and photometric observation of the pre-main sequence stars associated with the cometary shaped dark cloud Lynds 1622, and 12CO and 13CO observations of the cloud. We determined the effective temperatures and luminosities of 14 pre-main sequence stars associated with the cloud from their positions in the Hertzsprung--Russell diagram, as well as constructed their spectral energy distributions using optical, 2MASS and Spitzer IRAC and MIPS data. We derived physical parameters of L1622 from the molecular observations. Our results are not compatible with the assumption that L1622 lies on the near side of the Orion-Eridanus loop, but suggest that L1622 is as distant as Orion B. At a distance of 400 pc the mass of the cloud, derived from our CO data, is 1100 solar masses, its star formation efficiency is 1.8%, and the average age of its low-mass pre-main sequence star population is about 1 million years.

We describe a new survey for unbound hypervelocity stars (HVSs), stars traveling with such extreme velocities that dynamical ejection from a massive black hole (MBH) is their most likely origin. We investigate the possible contribution of unbound runaway stars, and show that the physical properties of binaries constrain low mass runaways to bound velocities. We measure radial velocities for HVS candidates with the colors of early A-type and late B-type stars. We report the discovery of 6 unbound HVSs with velocities and distances exceeding the conservative escape velocity estimate of Kenyon and collaborators. We additionally report 4 possibly unbound HVSs with velocities and distances exceeding the lower escape velocity estimate of Xue and collaborators. These discoveries increase the number of known HVSs by 60%-100%. Other survey objects include 19 newly identified z~2.4 quasars. One of the HVSs may be a horizontal branch star, consistent with the number of evolved HVSs predicted by Galactic center ejection models. Finding more evolved HVSs will one day allow a probe of the low-mass regime of HVSs and will constrain the mass function of stars in the Galactic center.

We report the discovery of four dusty cometary tails around low mass stars in two young clusters belonging to the W5 star forming region. Fits to the observed emission profiles from 24 micron observations with the Spitzer Space Telescope give tail lifetimes < 30 Myr, but more likely < 5 Myr. This result suggests that the cometary phase is a short lived phenomenon, occurring after photoevaporation by a nearby O star has removed gas from the outer disk of a young low mass star (see also Balog et al. 2006; Balog et al. 2008).

We describe Spitzer IRAC and MIPS observations of the nearby 25 Myr-old open cluster NGC 2232. Combining these data with ROSAT All-Sky Survey observations, proper motions, and optical photometry/spectroscopy, we construct a list of highly probable cluster members. We identify 1 A-type star, HD 45435, with definite excess emission at 4.5--24 micron indicative of debris from terrestrial planet formation. We also identify 2--4 late-type stars with possible 8 micron excesses, and 8 early-type stars with definite 24 micron excesses. Constraints on the dust luminosity and temperature suggest that the detected excesses are produced by debris disks. From our sample of B and A stars, stellar rotation appears correlated with 24 micron excess, a result expected if massive primordial disks evolve into massive debris disks. To explore the evolution of the frequency and magnitude of debris around A-type stars, we combine our results with data for other young clusters. The frequency of debris disks around A-type stars appears to increase from ~ 25% at 5 Myr to ~ 50--60% at 20--25 Myr. Older A-type stars have smaller debris disk frequencies: ~ 20% at 50--100 Myr. For these ages, the typical level of debris emission rises from 5--20 Myr and then declines. Because 24 micron dust emission probes icy planet formation around A-type stars, our results suggest that the frequency of icy planet formation is eta(i) > 0.5--0.6. Thus, most A-type stars (approx. 1.5--3 Msun) produce icy planets.

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.

We consider trends resulting from two formation mechanisms for short-period super-Earths: planet-planet scattering and migration. We model scenarios where these planets originate near the snow line in ``cold finger'' circumstellar disks. Low-mass planet-planet scattering excites planets to low periastron orbits only for lower mass stars. With long circularisation times, these planets reside on long-period eccentric orbits. Closer formation regions mean planets that reach short-period orbits by migration are most common around low-mass stars. Above ~1 Solar mass, planets massive enough to migrate to close-in orbits before the gas disk dissipates are above the critical mass for gas giant formation. Thus, there is an upper stellar mass limit for short-period super-Earths that form by migration. If disk masses are distributed as a power law, planet frequency increases with metallicity because most disks have low masses. For disk masses distributed around a relatively high mass, planet frequency decreases with increasing metallicity. As icy planets migrate, they shepherd interior objects toward the star, which grow to ~1 Earth mass. In contrast to icy migrators, surviving shepherded planets are rocky. Upon reaching short-period orbits, planets are subject to evaporation processes. The closest planets may be reduced to rocky or icy cores. Low-mass stars have lower EUV luminosities, so the level of evaporation decreases with decreasing stellar mass.

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 present a complete spectroscopic survey of 2414 2MASS-selected blue horizontal branch (BHB) candidates selected over 4300 deg^2 of the sky. We identify 655 BHB stars in this non-kinematically selected sample. We calculate the luminosity function of field BHB stars and find evidence for very few hot BHB stars in the field. The BHB stars located at a distance from the Galactic plane |Z|<4 kpc trace what is clearly a metal-weak thick disk population, with a mean metallicity of [Fe/H]= -1.7, a rotation velocity gradient of dv_{rot}/d|Z|= -28+-3.4 km/s in the region |Z|<6 kpc, and a density scale height of h_Z= 1.26+-0.1 kpc. The BHB stars located at 5<|Z|<9 kpc are a predominantly inner-halo population, with a mean metallicity of [Fe/H]= -2.0 and a mean Galactic rotation of -4+-31 km/s. We infer the density of halo and thick disk BHB stars is 104+-37 kpc^-3 near the Sun, and the relative normalization of halo to thick-disk BHB stars is 4+-1% near the Sun.

We use a semi-analytic circumstellar disk model that considers movement of the snow line through evolution of accretion and the central star to investigate how gas giant frequency changes with stellar mass. The snow line distance changes weakly with stellar mass; thus giant planets form over a wide range of spectral types. The probability that a given star has at least one gas giant increases linearly with stellar mass from 0.4 M_sun to 3 M_sun. Stars more massive than 3 M_sun evolve quickly to the main-sequence, which pushes the snow line to 10-15 AU before protoplanets form and limits the range of disk masses that form giant planet cores. If the frequency of gas giants around solar-mass stars is 6%, we predict occurrence rates of 1% for 0.4 M_sun stars and 10% for 1.5 M_sun stars. This result is largely insensitive to our assumed model parameters. Finally, the movement of the snow line as stars >2.5 M_sun move to the main-sequence may allow the ocean planets suggested by Leger et. al. to form without migration.

We report the discovery of accretion disks associated with ~ 13 Myr-old intermediate/low-mass stars in h and chi Persei. Optical spectroscopy of ~ 5000 stars in these clusters and a surrounding halo population reveal 32 A-K stars with H(alpha) emission. Matching these stars with 2MASS and optical photometry yields 25 stars with the highest probability of cluster membership and EW(H(alpha)) > 5 angstroms. Sixteen of these sources have EW(H(alpha)) > 10 angstroms. The population of accreting sources is strongly spectral type dependent: H(alpha) emission characteristic of accretion, especially strong accretion (EW(H(alpha)) > 10 angstroms), is much more prevalent around stars later than G0. Strong H(alpha) emission from accretion is typically associated with redder Ks-[8] colors. The existence of accreting pre-main sequence stars in h and chi Persei implies that circumstellar gas in some systems, especially those with primaries later than G5 spectral type, can last longer than 10-15 Myr.

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.

We measure the first stellar velocity dispersion of the Leo A dwarf galaxy, \sigma = 9.3 +- 1.3 km/s. We derive the velocity dispersion from the radial velocities of ten young B supergiants and two HII regions in the central region of Leo A. We estimate a projected mass of 8 +- 2.7 x10^7 solar masses within a radius of 2 arcmin, and a mass to light ratio of at least 20 +- 6 M_sun/L_sun. These results imply Leo A is at least ~80% dark matter by mass.

We report the detection of a radial velocity companion to SDSS J091709.55+463821.8, the lowest mass white dwarf currently known with M~0.17Msun. The radial velocity of the white dwarf shows variations with a semi-amplitude of 148.8 km/s and a period of 7.5936 hours, which implies a companion mass of M > 0.28Msun. The lack of evidence of a companion in the optical photometry forces any main-sequence companion to be smaller than 0.1Msun, hence a low mass main sequence star companion is ruled out for this system. The companion is most likely another white dwarf, and we present tentative evidence for an evolutionary scenario which could have produced it. However, a neutron star companion cannot be ruled out and follow-up radio observations are required to search for a pulsar companion.

We describe IRAC 3.6-8 mu_m observations and ground-based near-IR JHKs photometry from Mimir and 2MASS of the massive double cluster h & chi Persei complete to J=15.5 (M ~ 1.3 Msun. Within 25' of the cluster centers we detect \~11,000 sources with J < 15.5, ~7,000 sources with [4.5] < 15, and ~ 5000 sources with [8] < 14.5. In both clusters, the surface density profiles derived from the 2MASS data decline with distance from the cluster centers as expected for a bound cluster. Within 15' of the cluster centers, ~ 50% of the stars lie on a reddened ~ 13 Myr isochrone; at 15'-25' from the cluster centers, ~ 40% lie on this isochrone. Thus, the optical/2MASS color-magnitude diagrams indicate that h & chi Per are accompanied by a halo population with roughly the same age and distance as the two dense clusters. The double cluster lacks any clear IR excess sources for J < 13.5 (~ 2.7 Msun). Therefore, disks around high-mass stars disperse prior to ~ 10^{7} yr. At least 2-3% of the fainter cluster stars have strong IR excess at both [5.8] and [8]. About 4-8% of sources slightly more massive than the Sun (~ 1.4 Msun) have IR excesses at [8]. Combined with the lack of detectable excesses for brighter stars, this result suggests that disks around lower-mass stars have longer lifetimes. The IR excess population also appears to be larger at longer IRAC bands ([5.8], [8]) than at shorter IRAC/2MASS bands ($K_{s}$, [4.5]), a result consistent with an inside-out clearing of disks.

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 discuss the serendipitous discovery of an extremely metal poor galaxy, SDSS 0809+1729, classified as a star in the Sloan Digital Sky Survey (SDSS). The galaxy has a redshift z = 0.0441 and a B-band absolute magnitude M_B = -17.1. With a metallicity of log(O/H)+12~7.44 or ~1/20 solar, this galaxy is among the 10 most metal poor emission-line galaxies known. SDSS 0809+1729 is a blue compact galaxy (BCG) with a stellar age of ~4.5 Myr, a star formation rate of 0.18 Msun/yr, and a large gas-phase electron density ~367/cm^3. Similar values of these parameters are common among other extremely metal poor galaxies, including IZw18. SDSS 0809+1729 is, however, unusual among BCGs because it lies in the same region of the luminosity-metallicity diagram as the two lowest metallicity long duration gamma ray burst (GRB) hosts. For a given B-band luminosity, both nearby GRB hosts and SDSS 0809+1729 have systematically lower metallicities than dwarf irregulars, the majority of BCGs, and normal star-forming galaxies. Because the star formation properties of SDSS 0809+1729 are similar to nearby long duration GRB hosts, SDSS 0809+1729 may be a potential GRB host. Identification of larger samples of similar extremely metal poor objects may provide important insights into the conditions required to produce long duration GRBs.

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 discuss spectroscopy and infrared photometry for a complete sample of ~ 800 galaxies in close pairs objectively selected from the CfA2 redshift survey. We use 2MASS to compare near infrared color-color diagrams for our sample with the Nearby Field Galaxy Sample and with a set of IRAS flux-limited pairs from Surace et al. We construct a basic statistical model to explore the physical sources of the substantial differences among these samples. The model explains the spread of near infrared colors and is consistent with a picture where central star formation is triggered by the galaxy-galaxy interaction before a merger occurs. For 160 galaxies we report new, deep JHK photometry within our spectroscopic aperture and we use the combined spectroscopic and photometric data to explore the physical conditions in the central bursts. We find a set of objects with H-K >= 0.45 and with a large F(FIR)/F(H). We interpret the very red H-K colors as evidence for 600-1000 K dust within compact star-forming regions, perhaps similar to super-star clusters identified in individual well-studied interacting galaxies. The galaxies in our sample are candidate ``hidden'' bursts or, possibly, ``hidden'' AGN. Over the entire pair sample, both spectroscopic and photometric data show that the specific star formation rate decreases with the projected separation of the pair. The data suggest that the near infrared color-color diagram is also a function of the projected separation; all of the objects with central near infrared colors indicative of bursts of star formation lie at small projected separation.

Stellar explosions such as novae and supernovae produce most of the heavy elements in the Universe. Although the onset of novae from runaway thermonuclear fusion reactions on the surface of a white dwarf in a binary star system is understood[1], the structure, dynamics, and mass of the ejecta are not well known. In rare cases, the white dwarf is embedded in the wind nebula of a red-giant companion; the explosion products plow through the nebula and produce X-ray emission. Early this year, an eruption of the recurrent nova RS Ophiuchi[2,3] provided the first opportunity to perform comprehensive X-ray observations of such an event and diagnose conditions within the ejecta. Here we show that the hard X-ray emission from RS Ophiuchi early in the eruption emanates from behind a blast wave, or outward-moving shock wave, that expanded freely for less than 2 days and then decelerated due to interaction with the nebula. The X-rays faded rapidly, suggesting that the blast wave deviates from the standard spherical shell structure[4-6]. The early onset of deceleration indicates that the ejected shell had a low mass, the white dwarf has a high mass[7], and that RS Ophiuchi is a progenitor of the type of supernova integral to studies of the expansion of the universe.

We present low-resolution turbulence profiles of the atmosphere above Dome C, Antarctica, measured with the MASS instrument during 25 nights in March-May 2004. Except for the lowest layer, Dome C has significantly less turbulence than Cerro Tololo and Cerro Pachon. In particular, the integrated turbulence at 16 km is always less than the median values at the two Chilean sites. From these profiles we evaluate the photometric noise produced by scintillation, and the atmospheric contribution to the error budget in narrow-angle differential astrometry. In comparison with the two mid-latitude sites in Chile, Dome C offers a potential gain of about 3.6 in both photometric precision (for long integrations) and narrow-angle astrometry precision. These gain estimates are preliminary, being computed with average wind-speed profiles, but the validity of our approach is confirmed by independent data. Although the data from Dome C cover a fairly limited time frame, they lend strong support to expectations that Dome C will offer significant advantages for photometric and astrometric studies.

We discuss our targeted search for hypervelocity stars (HVSs), stars traveling with velocities so extreme that dynamical ejection from a massive black hole is their only suggested origin. Our survey, now half complete, has successfully identified a total of four probable HVSs plus a number of other unusual objects. Here we report the most recently discovered two HVSs: SDSS J110557.45+093439.5 and possibly SDSS J113312.12+010824, traveling with Galactic rest-frame velocities at least +508+-12 and +418+-10 km/s, respectively. The other late B-type objects in our survey are consistent with a population of post main-sequence stars or blue stragglers in the Galactic halo, with mean metallicity [Fe/H]=-1.3 and velocity dispersion 108+-5 km/s. Interestingly, the velocity distribution shows a tail of objects with large positive velocities that may be a mix of low-velocity HVSs and high-velocity runaway stars. Our survey also includes a number of DA white dwarfs with unusually red colors, possibly extremely low mass objects. Two of our objects are B supergiants in the Leo A dwarf, providing the first spectroscopic evidence for star formation in this dwarf galaxy within the last ~30 Myr.

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.

Hypervelocity stars (HVSs) travel with velocities so extreme that dynamical ejection from a massive black hole is their only suggested origin. Following our discovery of the first HVS, we have undertaken a dedicated survey for more HVSs in the Galactic halo and present here the resulting discovery of two new HVSs: SDSS J091301.0+305120 and SDSS J091759.5+672238, traveling with Galactic rest-frame velocities at least +558+-12 and +638+-12 km/s, respectively. Assuming the HVSs are B8 main sequence stars, they are at distances ~75 and ~55 kpc, respectively, and have travel times from the Galactic Center consistent with their lifetimes. The existence of two B8 HVSs in our 1900 deg^2 survey, combined with the Yu & Tremaine HVS rate estimates, is consistent with HVSs drawn from a standard initial mass function but inconsistent with HVS drawn from a truncated mass function like the one in the top-heavy Arches cluster. The travel times of the five currently known HVSs provide no evidence for a burst of HVSs from a major in-fall event at the Galactic Center in the last \~160 Myr.

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.

The recent discovery of exceptional seeing conditions at Dome C, Antarctica, raises the possibility of constructing an optical observatory there with unique capabilities. However, little is known from an astronomer's perspective about the optical sky brightness and extinction at Antarctic sites. We review the contributions to sky brightness at high-latitude sites, and calculate the amount of usable dark time at Dome C. We also explore the implications of the limited sky coverage of high-latitude sites and review optical extinction data from the South Pole. Finally, we examine the proposal of Baldry and Bland-Hawthorn (2001) to extend the amount of usable dark time through the use of polarising filters.

We describe observational evidence for a new kind of interacting-binary-star outburst that involves both an accretion instability and an increase in thermonuclear shell burning on the surface of an accreting white dwarf. We refer to this new type of eruption as a combination nova. In late 2000, the prototypical symbiotic star Z Andromedae brightened by roughly two magnitudes in the optical. We observed the outburst in the radio with the VLA and MERLIN, in the optical both photometrically and spectroscopically, in the far ultraviolet with FUSE, and in the X-rays with both Chandra and XMM. The two-year-long event had three distinct stages. During the first stage, the optical rise closely resembled an earlier, small outburst that was caused by an accretion-disk instability. In the second stage, the hot component ejected an optically thick shell of material. In the third stage, the shell cleared to reveal a white dwarf whose luminosity remained on the order of 10^4 Lsun for approximately one year. The eruption was thus too energetic to have been powered by accretion alone. We propose that the initial burst of accretion was large enough to trigger enhanced nuclear burning on the surface of the white dwarf and the ejection of an optically thick shell of material. This outburst therefore combined elements of both a dwarf nova and a classical nova. Our results have implications for the long-standing problem of producing shell flashes with short recurrence times on low-mass white dwarfs in symbiotic stars.

We discuss a 175 deg^2 spectroscopic survey for blue horizontal branch (BHB) stars in the Galactic halo. We use the Two Micron All Sky Survey (2MASS) and the Sloan Digital Sky Survey (SDSS) to select BHB candidates, and find that the 2MASS and SDSS color-selection is 38% and 50% efficient, respectively, for BHB stars. Our samples include one likely run-away B7 star 6 kpc below the Galactic plane. The global properties of the BHB samples are consistent with membership in the halo population: the median metallicity is [Fe/H]=-1.7, the velocity dispersion is 108 km/s, and the mean Galactic rotation of the BHB stars 3<|z|<15 kpc is -4 +- 30 km/s. We discuss the theoretical basis of the Preston, Shectman & Beers M_V-color relation for BHB stars, and conclude that intrinsic shape of the BHB M_V-color relation results from the physics of stars on the horizontal branch. We calculate the luminosity function for the field BHB star samples using the Efstathiou, Ellis, & Peterson maximum-likelihood method which is unbiased by density variations. The field BHB luminosity function exhibits a steep rise at bright luminosities, a peak between 0.8 < M_V < 1.0, and a tail at faint luminosities. We compare the field BHB luminosity functions with the luminosity functions derived from sixteen different globular cluster BHBs. Kolmogorov-Smirnov tests suggest that field BHB stars and BHB stars in globular clusters share a common distribution of luminosities, with the exception of globular clusters with extended BHBs.

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.

We have discovered a star, SDSS J090745.0+024507, leaving the Galaxy with a heliocentric radial velocity of +853+-12 km/s, the largest velocity ever observed in the Milky Way halo. The star is either a hot blue horizontal branch star or a B9 main sequence star with a heliocentric distance ~55 kpc. Corrected for the solar reflex motion and to the local standard of rest, the Galactic rest-frame velocity is +709 km/s. Because its radial velocity vector points 173.8 deg from the Galactic center, we suggest that this star is the first example of a hyper-velocity star ejected from the Galactic center as predicted by Hills and later discussed by Yu & Tremaine. The star has [Fe/H]~0, consistent with a Galactic center origin, and a travel time of <80 Myr from the Galactic center, consistent with its stellar lifetime. If the star is indeed traveling from the Galactic center, it should have a proper motion of 0.3 mas/yr observable with GAIA. Identifying additional hyper-velocity stars throughout the halo will constrain the production rate history of hyper-velocity stars at the Galactic center.

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 performed extensive, multi-wavelength observations of the prototypical symbiotic star Z Andromedae between 2000 and 2003, during a large eruption. The rise to optical maximum occurred in three distinct stages. During the first stage, the rise was very similar to an earlier, small outburst which we determined was due to an accretion-disk instability. In the second stage, an optically thick shell of material was ejected, and in the third stage, the shell cleared to reveal a white dwarf whose luminosity was roughly 10^4 Lsun. We suggest that the outburst was powered by an increase in the rate of nuclear burning on the white-dwarf surface, triggered by a sudden burst of accretion. This outburst thus combined elements of both dwarf novae and classical novae.

We describe the results of two near infrared (K-band) imaging surveys and a three color (JHK) survey of the vicinity of NGC 7538. The limiting magnitudes are K ~ 16.5 and K ~ 17.5 mag for the K-band surveys and K ~ 15 mag for the JHK survey. We identify more than 2000 and 9000 near-infrared (NIR) sources on the images of the two K-band surveys and 786 NIR sources in the JHK survey. From color-color diagrams, we derive a reddening law for background stars and identify 238 stars with NIR excesses. Contour maps indicate a high density peak coincident with a concentration of stars with NIR excesses. We identify this peak as a young, embedded cluster and confirm this result with the K-band luminosity function, color histograms, and color-magnitude diagrams. The center of the cluster is at RA = 23:13:39.34, DEC = 61:29:18.9. The cluster radius is $\sim$ 3' ~ 2.5 pc for an adopted distance, d ~ 2.8 kpc. For d = 2.8 kpc, and reddening, E_{J-K} = 0.55 mag, the slope of the logarithmic K-band luminosity function (KLF) of the cluster, s ~ 0.32 +- 0.03, agrees well with previous results for L1630 (s = 0.34) and M17 (s = 0.26).

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.

Spectra of the bright symbiotic star AG Draconis (BD+67D922) in the wavelength range 905-1187 A obtained with the Far Ultraviolet Spectroscopic Explorer (FUSE) are presented. The spectra show a number of narrow, nebular emission lines, together with a uniform continuum from the hot component of the system, and numerous interstellar absorption lines. We infer the existence of Ne VIII in the AG Dra nebula through the identification of the Ne VII 973.3 A recombination line. The emission line spectrum is dominated by intense lines of O VI, but also shows weaker lines from highly-ionized ions including Ne V, Ne VI, S IV, S VI. Members of the He II Balmer series can be identified up to n=20. Lines of Fe II and Fe III fluoresced by O VI 1032 are identified at wavelengths 1141.172 A and 1142.429 A, respectively. The emission lines are shown to be produced in a plasma with an electron temperature of 20-30,000 K, photoionized by the white dwarf. The Ne VI 997/999 ratio shows that this ion and all others except perhaps Ne VII are formed at least 300 white dwarf radii from the white dwarf. Revised wavelengths for the Ne V 2s^2 2p^2 ^3P_0,1 - 2s 2p^3 ^5S_2 and Ne VI 2s^2 2p ^2P_J - 2s 2p^2 ^4P_J' transitions are published.

Collimated outflows from accreting white dwarfs have an important role to play in the study of astrophysical jets. Observationally, collimated outflows are associated with systems in which material is accreted though a disk. Theoretically, accretion disks provide the foundation for many jet models. Perhaps the best-understood of all accretion disks are those in cataclysmic variable stars (CVs). Since the disks in other accreting white-dwarf (WD) binaries are probably similar to CV disks (at least to the extent that one does not expect complications such as, for example, advection-dominated flows), with WD accretors one has the advantage of a relatively good grasp of the region from which the outflows are likely to originate. We briefly compare the properties of the three main classes of WD accretors, two of which have members that produce jets, and review the cases of three specific jet-producing WD systems.

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 2MASS photometry to select blue horizontal branch (BHB) candidates covering the sky |b|>15 deg. A 12.5<J<15.5 sample of BHB stars traces the thick disk and inner halo to d<9 kpc, with a density comparable to that of M giant stars. We base our sample selection strategy on the Century Survey Galactic Halo Project, a survey that provides a complete, spectroscopically-identified sample of blue stars to a similar depth as the 2MASS catalog. We show that a -0.20<(J-H)_0<0.10, -0.10<(H-K)_0<0.10 color-selected sample of stars is 65% complete for BHB stars, and is composed of 47% BHB stars. We apply this photometric selection to the full 2MASS catalog, and see no spatial overdensities of BHB candidates at high Galactic latitude |b|>50 deg. We insert simulated star streams into the data and conclude that the high Galactic latitude BHB candidates are consistent with having no ~5 deg wide star stream with density greater than 0.33 objects deg^-2 at the 95% confidence level. The absence of structure suggests there have been no major accretion events in the inner halo in the last few Gyr. However, at low Galactic latitudes a two-point angular correlation analysis reveals structure on angular scales <1 deg. This structure is apparently associated with stars in the thick disk, and has a physical scale of 10-100 pc. Interestingly, such structures are expected by cosmological simulations that predict the majority of the thick disk may arise from accretion and disruption of satellite mergers.

We review the current status of QSO clustering measurements, particular with respect to their relevance in understanding AGN physics. Measurements based on the 2dF QSO Redshift Survey (2QZ) find a scale length for QSO clustering of s_0=5.76(+0.17-0.27) h-1 Mpc at a redshift <z>~1.5, very similar to low redshift galaxies. There is no evidence of evolution in the clustering of QSOs from z~0.5 to z~2.2. This lack of evolution and low clustering amplitude suggests a short life time for AGN activity of the order ~10^6-10^7 years. Large surveys such at the 2QZ and SDSS also allow the the study of QSO environments in 3D for the first time (at least at low redshift), early results from this work seem to show no difference between the environments of QSOs and normal galaxies. Future studies e.g. measuring clustering as a function of black hole mass, and deep QSO surveys should provide further insight into the formation and evolution of AGN.

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.

The Century Survey Galactic Halo Project is a photometric and spectroscopic survey from which we select relatively blue stars (V-R<0.30 mag) as probes of the Milky Way halo. The Survey strip spans the range of Galactic latitude 35<b<88 deg, allowing us to study the nature of populations of stars and their systematic motions as a function of Galactic latitude. One of our primary goals is to use blue horizontal-branch stars to trace potential star streams in the halo, and to test the hierarchical model for the formation of the Galaxy. In this paper we discuss spectroscopy and multi-passband photometry for a sample of 764 blue stars in the Century Survey region. Our sample consists predominantly of A- and F-type stars. We describe our techniques for determination of radial velocities, effective temperatures, metallicities, and surface gravities. Based on these measurements, we derive distance estimates by comparison with a set of calibrated isochrones. We devote special attention to the classification of blue horizontal-branch stars, and compare the results obtained from the application of the techniques of Kinman et al., Wilhelm et al., and Clewley et al. We identify 55 blue horizontal-branch stars. Our large sample of stars also uncovers a number of unusual objects, including three carbon-enhanced stars, a late B-type star located 0.8 kpc above the Galactic plane, and a DZ white dwarf.

CH Cygni began producing rapid, stochastic optical variations with the onset of symbiotic activity in 1963. We use changes in this flickering between 1997 and 2000 to diagnose the state of the accretion disk during this time. In the 1998 high state, the luminosity of the B-band flickering component was typically more than 20 times higher than in the 1997 and 2000 low states. Therefore, the physical process or region that produces the flickering was also primarily responsible for the large optical flux increase in the 1998 high state. Assuming that the rapid, stochastic optical variations in CH Cygni come from the accretion disk, as in cataclysmic variable stars, a change in the accretion rate through the disk led to the 1998 bright state. All flickering disappeared in 1999, when the accreting WD was eclipsed by the red giant orbiting with a period of approximately 14 yr, according to the ephemeris of Hinkle et al. and the interpretation of Eyres et al. We did not find any evidence for periodic or quasi-periodic oscillations in the optical emission from CH Cygni in either the high or low state, and we discuss the implications for magnetic propeller models of this system. As one alternative to propeller models, we propose that the activity in CH Cygni is driven by accretion through a disk with a thermal-viscous instability, similar to the instabilities believed to exist in dwarf novae and suggested for FU Ori pre-main-sequence stars and soft X-ray transients.

We investigate the role of accretion in the production of jets in the symbiotic star CH Cygni. Assuming that the rapid stochastic optical variations in CH Cygni come from the accretion disk, as in cataclysmic variables, we use changes in this flickering to diagnose the state of the disk in 1997. At that time, CH Cyg dropped to a very low optical state, and Karovska et al. report that a radio jet was produced. For approximately one year after the jet production, the amplitude of the fastest (time scales of minutes) variations was significantly reduced, although smooth, hour-time-scale variations were still present. This light curve evolution indicates that the inner disk may have been disrupted, or emission from this region suppressed, in association with the mass-ejection event. We describe optical spectra which support this interpretation of the flickering changes. The simultaneous state change, jet ejection, and disk disruption suggests a comparison between CH Cygni and some black-hole-candidate X-ray binaries that show changes in the inner disk radius in conjunction with discrete ejection events on a wide range of time scales (e.g., the microquasar GRS 1915+105 and XTE J1550-564).

We report the results of a spectroscopic survey of 20 close T Tauri binaries in the Taurus-Auriga dark cloud where the separations between primaries and their secondaries are less than the typical size of a circumstellar disk around a young star. Analysis of low-resolution and medium-resolution STIS spectra yields the stellar luminosities, reddenings, ages, masses, mass accretion rates, IR excesses, and emission line luminosities for each star in each pair. We examine the ability of IR color excesses, H-alpha equivalent widths, [O I] emission, and veiling to distinguish between weak emission and classical T Tauri stars. Four pairs have one cTTs and one wTTs; the cTTs is the primary in three of these systems. This frequency of mixed pairs among the close T Tauri binaries is similar to the frequency of mixed pairs in wider young binaries. Extinctions within pairs are usually similar; however, the secondary is more heavily reddened than the primary in some systems, where it may be viewed through the primary's disk. Mass accretion rates of primaries and secondaries are strongly correlated, and H-alpha luminosities, IR excesses, and ages also correlate within pairs. Primaries tend to have somewhat larger accretion rates than their secondaries do, and are typically slightly older than their secondaries according to three different sets of modern pre-main-sequence evolutionary tracks. Age differences for XZ Tau and FS Tau, systems embedded in reflection nebulae, are striking; the secondary in each pair is less massive but more luminous than the primary. The stellar masses of the UY Aur and GG Tau binaries measured from their rotating molecular disks are about 30% larger than the masses inferred from the spectra and evolutionary tracks.

Galaxy-galaxy interactions rearrange the baryons in galaxies and trigger substantial star formation; the aggregate effects of these interactions on the evolutionary histories of galaxies in the Universe are poorly understood. We combine B and R-band photometry and optical spectroscopy to estimate the strengths and timescales of bursts of triggered star formation in the centers of 190 galaxies in pairs and compact groups. Based on an analysis of the measured colors and EW(H-alpha), we characterize the pre-existing and triggered populations separately. The best-fitting burst scenarios assume stronger reddening corrections for line emission than for the continuum and continuous star formation lasting for \gtrsim a hundred Myr. The most realistic scenarios require an initial mass function that is deficient in the highest-mass stars. The color of the pre-existing stellar population is the most significant source of uncertainty. Triggered star formation contributes substantially (probably >= 50%) to the R-band flux in the central regions of several galaxies; tidal tails do not necessarily accompany this star formation. Many of the galaxies in our sample have bluer centers than outskirts, suggesting that pre- or non-merger interactions may lead to evolution along the Hubble sequence. These objects would appear blue and compact at higher redshifts; the older, redder outskirts of the disks would be difficult to detect. Our data indicate that galaxies with larger separations on the sky contain weaker, and probably older, bursts of star formation on average. However, confirmation of these trends requires further constraints on the colors of the older stellar populations and on the reddening for individual galaxies.

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 present evidence for approximately 400-d variations in the radial velocity of HD 154791 (V934 Her), the suggested optical counterpart of 4U 1700+24. The variations are correlated with the previously reported approximately 400 d variations in the X-ray flux of 4U 1700+24, which supports the association of these two objects, as well as the identification of this system as the second known X-ray binary in which a neutron star accretes from the wind of a red giant. The HD 154791 radial velocity variations can be fit with an eccentric orbit with period 404 +/- 3 d, amplitude K=0.75 +/- 0.12 km/s and eccentricity e=0.26 +/- 0.15. There are also indications of variations on longer time scales >~ 2000 d. We have re-examined all available ASM data following an unusually large X-ray outburst in 1997-98, and confirm that the 1-d averaged 2-10 keV X-ray flux from 4U 1700+24 is modulated with a period of 400 +/- 20 d. The mean profile of the persistent X-ray variations was approximately sinusoidal, with an amplitude of 0.108 +/- 0.012 ASM count/s (corresponding to 31% rms). The epoch of X-ray maximum was approximately 40 d after the time of periastron according to the eccentric orbital fit. If the 400 d oscillations from HD 154791/4U 1700+24 are due to orbital motion, then the system parameters are probably close to those of the only other neutron-star symbiotic-like binary, GX 1+4. We discuss the similarities and differences between these two systems.

We analyze spectra of 44 emission line stars detected in a low resolution optical spectroscopic survey of the galactic open cluster NGC 6871. The survey of 1217 stars is complete to V = 14.9 and includes stars with V < 16.5 between the zero age main sequence (ZAMS) and the 10^7 yr pre-main sequence (PMS) isochrone. Of the 44 emission line stars in this survey, 28 show obvious emission in H alpha and 16 have weak H alpha absorption (compared to H gamma). We use the reddening to separate foreground and background stars from the cluster members; the position in the HR diagram (HRD) or the presence of forbidden emission lines ([NII], [SII]) then yields the evolutionary status of the emission line stars. A comparison of the H alpha spectral index distribution in NGC 6871 with the distribution of young stars in the Taurus-Auriga molecular cloud indicates that the late type PMS stars in our sample may be weak line T Tauri stars (wTTs). Many of these stars show [SII] (36%) and [NII] (45%) emission.

We describe a near-infrared imaging survey of Globule 2 in the Coalsack. This Bok globule is the highest density region of this southern hemisphere molecular cloud and is the most likely location for young stars in this complex. The survey is complete for K < 14.0, H < 14.5, and J < 15.5, several magnitudes more sensitive than previous observations of this globule. From the large number of background stars, we derive an accurate near-infrared extinction law for the cloud. Our result, E_{J-H}/E_{H-K} = 2.08 \pm 0.03, is significantly steeper than results for other southern clouds. We use the J-H/H-K color-color diagram to identify two potential young stars with K < 14.0 in the region. We apply H-band star counts to derive the density profile of the Coalsack Globule 2 and use a polytropic model to describe the internal structure of this small cloud. For a gas temperature T \sim 15 K, this globule is moderately unstable.

We describe new optical photometric and spectroscopic observations of the semi-detached eclipsing binary BG Geminorum. A large change in the amount of Mg I absorption at secondary maximum indicates the presence of cool material in the outer edge of the disk surrounding the unseen primary star. Detection of weak He I emission implies a hot radiation source at the inner edge of the disk. If the velocity variations in the Hbeta emission line track the orbital motion of the primary star, the primary star has an orbital semiamplitude of K1 = K(Hbeta) = 16.0 +- 4.6 km/sec. This result yields a mass ratio, q = 0.22 +- 0.07, consistent with the q = 0.1 derived from the large ellipsoidal variation. Despite this progress, the nature of the primary star - B-type star or black hole - remains uncertain.

We describe near-IR imaging data for a sample of 23 class I sources in the Taurus-Auriga dark clouds. Combining our data with previous photometry, we detect brightness variations of 0.1-0.5 mag in many sources. The near-IR morphologies are consistent with mm continuum measurements. Most (60%) of the sample are true protostars; the rest may be objects in transition between class I and class II, T Tauri stars with edge-on disks, or heavily reddened T Tauri stars.

This paper reviews coagulation models for planet formation in the Kuiper Belt, emphasizing links to recent observations of our and other solar systems. At heliocentric distances of 35-50 AU, single annulus and multiannulus planetesimal accretion calculations produce several 1000 km or larger planets and many 50-500 km objects on timescales of 10-30 Myr in a Minimum Mass Solar Nebula. Planets form more rapidly in more massive nebulae. All models yield two power law cumulative size distributions, N_C propto r^{-q} with q = 3.0-3.5 for radii larger than 10 km and N_C propto r^{-2.5} for radii less than 1 km. These size distributions are consistent with observations of Kuiper Belt objects acquired during the past decade. Once large objects form at 35-50 AU, gravitational stirring leads to a collisional cascade where 0.1-10 km objects are ground to dust. The collisional cascade removes 80% to 90% of the initial mass in the nebula in roughly 1 Gyr. This dust production rate is comparable to rates inferred for alpha Lyr, beta Pic, and other extrasolar debris disk systems.

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.

We investigate the observational signatures of a straightforward evolutionary scenario for protoplanetary disks in which the disk mass of small (50 micron) particles decreases homologously with time, but the disk structure and stellar parameters do not change. Our goal is to identify optimal infrared spectral indicators of the existence of disks, their structure, and mass evolution that may be tested with the upcoming SIRTF mission. We present simulated spectral energy distributions and colors over a wide range of masses. The SED is most sensitive to disk mass in the far-IR and longer wavelengths, which is already known from millimeter and radio observations. As the disk mass decreases, the excess emission of the disk over the stellar photosphere diminishes more rapidly at the longest rather than at short wavelengths. At near-infrared wavelengths, the disk remains optically thick to stellar radiation over a wide range of disk mass, resulting in a slower decline in the SED in this spectral regime. Therefore, near-IR excesses (K-L) provide a robust means of detecting disks in star clusters down to 1E-7 solar masses, while the far-IR excess probes the disk mass. Reducing the disk mass results in a clear progression in color-color diagrams with low mass disks displaying the bluest colors. We interpret color-color diagrams for Taurus-Auriga sources in the context of decreasing disk mass.

We describe the first results from a program of multi-wavelength monitoring of symbiotic binary stars in outburst. The prototypical symbiotic star Z Andromedae was observed extensively at radio, optical, far-UV, and X-ray wavelengths during its current (ongoing) outburst. These observations provide evidence for an initial obscuration or the source of ionizing radiation by an optically thick shell of material, outflow with velocities of 100s of km/s, plus collisional excitation and shock heating of the red-giant wind. In addition, the optical light curve reveals variability on time scales as short as days on both the rise to and decline from optical maximum.

Ground based imaging, imaging polarimetry, and recent Hubble Space Telescope WFPC2 and NICMOS images of protostars have revealed very complex scattered light patterns that cannot be entirely explained by 2-D radiation transfer models. We present here for the first time radiation transfer models of T Tau and IRAS 04016+2610 that are fully 3-D, with the aim of investigating the effects on image morphology of multiple illuminating sources and infalling envelopes that have been shaped by multiple outflows. For T Tau we have constructed scattered light models where the illumination of the surrounding envelope is by a binary with each source surrounded by its own small circumstellar disk or envelope. We find that the asymmetries in the WFPC2 image of T Tau can be reproduced if the disks in the binary system are misaligned, consistent with a recently discovered bipolar outflow believed to originate from the secondary. For IRAS 04016+2610 we find that the observed scattered light pattern can be reproduced by scattering in an envelope with cavities carved by two sets of bipolar outflows, suggestive of an embedded binary system.

We discuss H-band (1.65 micrometer) near-infrared photometry of the central 3x3h^{-2} Mpc^2 of Abell 1644 to a limiting M_H~M^*_H+3 (throughout this paper H_0=100 h km s^{-1} Mpc^{-1}). There are 861 galaxies in the photometric survey region. We also measured radial velocities of 155 galaxies; 141 of these are cluster members within 2.44 h^{-1} Mpc of the cluster center. The cluster velocity dispersion of sigma~1000 km s^{-1} remains constant out to the limiting radius. We find no evidence for substructure in the cluster. The cluster mass within R=2.4 h^{-1} Mpc is 7.6+-1.3x10^{14} h^{-1} M_\odot. We compute the cluster luminosity function; the Schechter parameters \alpha=-1.14\pm0.08 and M^*_H=-24.3\pm0.2 (with h=0.5) agree well with other H-band luminosity functions. From the virial theorem and the caustic method we compute one of the first mass-to-light ratios at H; the result is M/L_H=82-127 h M_\odot/L_\odot within 1.5 h^{-1} Mpc. This ratio corresponds to 374-579 h M_\odot/L_\odot at R. The agreement of our IR measurement with previous $M/L$ determinations indicates that at low redshift dust and young stellar populations may produce only negligible systematic errors in optical mass-to-light ratios.

We analyze optical and ultraviolet observations of the symbiotic binary AG Pegasi acquired during 1992-97. The bolometric luminosity of the hot component declined by a factor of 2-3 from 1980-1985 to 1997. Since 1992, the effective temperature of the hot component may have declined by 10%-20%, but this decline is comparable to the measurement errors. Optical observations of H-beta and He I emission show a clear illumination effect, where high energy photons from the hot component ionize the outer atmosphere of the red giant. Simple illumination models generally account for the magnitude of the optical and ultraviolet emission line fluxes. High ionization emission lines - [Ne V], [Mg V], and [Fe VII] - suggest mechanical heating in the outer portions of the photoionized red giant wind. This emission probably originates in a low density region $\sim$ 30-300 AU from the central binary.

We describe an L-band photometric survey of 0.5 square deg of the Cha I dark cloud. The survey has a completeness limit of L < 11.0. Our survey detects 124 sources, including all known pre-main sequence stars with L < 11. The fraction of sources with near-IR excess emission is 58% +- 4% for K = 9-11. Cha I sources have bluer H-K and K-L colors than pre-main sequence stars in Taurus-Auriga. These sources also have a strong correlation between EW(H-alpha) and K-L. Stars with K-L < 0.6 have weak H-alpha emission; stars with K-L > 0.6 have strong H-alpha emission. Because many Cha I sources are heavily reddened, this division between weak emission T Tauri stars and classical T Tauri stars occurs at a redder K-L than in Taurus-Auriga.

We describe a near-infrared imaging survey covering approximately 1 square deg of the Chamaeleon I dark cloud. The survey is complete for K < 15.0, H < 16.0, and J < 16.5, roughly two magnitudes more sensitive than previous large scale surveys. We use the large number of background stars detected to derive an accurate near-infrared extinction law for the cloud and select new candidate members with near-infrared color excesses. We list about 100 candidates of the cloud with K > 12.0, based on their positions in the J-H, H-K color-color diagram. These new stars have low luminosities (K > 12 -- 16, H-K > 0.5 -- 1.5) and may have masses close to or even below the hydrogen burning limit.

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.)

We investigate the constraints that Olbers Paradox, applied to the Zodiacal Background as measured from space, sets on outer solar system objects. If extended to very faint limits, R = 40--50 mag, the steep optical luminosity function (LF) of Kuiper Belt objects (KBOs) at R less than 26 mag implies an infinitely bright night sky. Small KBOs with radii of 1 mm to 1 km must have a size distribution n(r) proportional to r^{-a}, with a = 3.5 or smaller to satisfy the known limits on the sky-surface brightness at optical and far-infrared wavelengths. Improved limits on the measured KBO surface brightness can yield direct estimates of the albedo, temperature, and size distribution for small KBOs in the outer solar system.

This paper reviews the theory of protostellar debris disks. After a brief introduction to accretion disk physics, I describe coagulation models of planet formation in the outer regions of planetesimal disks. Coagulation models for the Kuiper Belt produce Pluto-sized objects on timescales of 10-40 Myr. These models yield size distributions which agree with observations of Kuiper Belt objects with red magnitudes, R = 20-27. Velocity stirring models for other debris disk systems demonstrate that 500 km or larger objects can stir the velocities of small objects up to the shattering limit needed to begin a collisional cascade.

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.

We discuss V and R band photometry for 67% of the Sullivan et al. 2000 SA57 ultraviolet-selected galaxy sample. In a sample of 176 UV-selected galaxies, Sullivan et al. 2000 find that 24% have (UV-B) colors too blue for consistency with starburst spectral synthesis models. We propose that these extreme blue, UV excess galaxies are Wolf-Rayet (WR) galaxies, starburst galaxies with strong UV emission from WR stars. We measure a median (V-R)=0.38+-0.06 for the UV-selected sample, bluer than a sample optically selected at R but consistent with starburst and WR galaxy colors. We demonstrate that redshifted WR emission lines can double or triple the flux through the UV bandpass at high redshifts. Thus the (UV-B) color of a WR galaxy can be up to 1.3 mag bluer at high redshift, and the expected selection function is skewed to larger redshifts. The redshift distribution of the extreme blue, UV excess galaxies matches the selection function we predict from the properties of WR galaxies.

We analyze new optical spectra of a sample of 502 galaxies in close pairs and n-tuples, separated by <= 50/h kpc. We extracted the sample objectively from the CfA2 redshift survey, without regard to the surroundings of the tight systems. We probe the relationship between star formation and the dynamics of the systems of galaxies. The equivalent widths of H\alpha (EW(H\alpha) and other emission lines anti-correlate strongly with pair spatial separation (\Delta D) and velocity separation. We use the measured EW(H\alpha) and the starburst models of Leitherer et al. to estimate the time since the most recent burst of star for- mation began for each galaxy. In the absence of a large contribution from an old stellar population to the continuum around H\alpha, the observed \Delta D -- EW(H\alpha) correlation signifies that starbursts with larger separations on the sky are, on average, older. By matching the dynamical timescale to the burst timescale, we show that the data support a simple picture in which a close pass initiates a starburst; EW(H\alpha) decreases with time as the pair separation increases, accounting for the anti-correlation. This picture leads to a method for measuring the duration and the initial mass function of interaction-induced starbursts: our data are compatible with the starburst and orbit models in many respects, as long as the starburst lasts longer than \sim10^8 years and the delay between the close pass and the initiation of the starburst is less than a few \times 10^7 years. If there is no large contribution from an old stellar population to the continuum around H\alpha the Miller-Scalo and cutoff (M <= 30 M_\sun) Salpeter initial mass functions fit the data much better than a standard Salpeter IMF. (Abridged.)

We describe optical photometric and spectroscopic observations of the bright variable BG Geminorum. Optical photometry shows a pronounced ellipsoidal variation of the K0 I secondary, with amplitudes of ~0.5 mag at VRI and a period of 91.645 days. A deep primary eclipse is visible for wavelengths < 4400A; a shallower secondary eclipse is present at longer wavelengths. Eclipse timings and the radial velocity curve of the K0 secondary star indicate an interacting binary where a lobe-filling secondary, M_2 ~ 0.5 Msun, transfers material into a extended disk around a massive primary, M_1 ~ 4.5 Msun. The primary star is either an early B-type star or a black hole. If it did contain a black hole, BG Gem would be the longest period black hole binary known by a factor of 10, as well as the only eclipsing black hole binary system.

We analyze new and published optical photometric data of FU Orionis, an eruptive pre-main sequence star. The outburst consists of a 5.5 mag rise at B with an e-folding timescale of roughly 50 days. The rates of decline at B and V are identical, 0.015 +- 0.001 mag per yr. Random fluctuations superimposed on this decline have an amplitude of 0.035 +- 0.005 mag at V and occur on timescales of 1 day or less. Correlations between V and the color indices U-B, B-V, and V-R indicate that the variable source has the optical colors of a G0 supergiant. We associate this behavior with small amplitude flickering of the inner accretion disk.

We describe planetesimal accretion calculations for the dusty ring observed in the nearby A0 star HR 4796A. Models with initial masses of 10-20 times the minimum mass solar nebula produce a ring of width 7-15 AU and height 0.3-0.6 AU at 70 AU in roughly 10 Myr. The ring has a radial optical depth of 1. These results agree with limits derived from infrared images and from the excess infrared luminosity.

We describe calculations of the evolution of an ensemble of small planetesimals in the outer solar system. In a solar nebula with a mass of several times the Minimum Mass Solar Nebula, objects with radii of 100-1000 km can form on timescales of 10-100 Myr. Model luminosity functions derived from these calculations agree with current observations of bodies beyond the orbit of Neptune (Kuiper Belt objects). New surveys with current and planned instruments can place better constraints on the mass and dynamics of the solar nebula by measuring the luminosity function at red magnitudes of 28 or larger.

We describe new planetesimal accretion calculations in the Kuiper Belt that include fragmentation and velocity evolution. All models produce two power law cumulative size distributions, N_C propto r^{-q}, with q = 2.5 for radii less than 0.3-3 km and q = 3 for radii exceeding 1-3 km. The power law indices are nearly independent of the initial mass in the annulus, the initial eccentricity of the planetesimal swarm, and the initial size distribution of the planetesimal swarm. The transition between the two power laws moves to larger radii as the initial eccentricity increases. The maximum size of objects depends on their intrinsic tensile strength; Pluto formation requires a strength exceeding 300 erg per gram. Our models yield formation timescales for Pluto-sized objects of 30-40 Myr for a minimum mass solar nebula. The production of several `Plutos' and more than 10^5 50 km radius Kuiper Belt objects leaves most of the initial mass in 0.1-10 km radius objects that can be collisionally depleted over the age of the solar system. These results resolve the puzzle of large Kuiper Belt objects in a small mass Kuiper Belt.

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.

We report results of a near--infrared imaging survey of L977, a dark cloud in Cygnus seen in projection against the plane of the Milky Way. We use measurements of the near--infrared color excess and positions of the 1628 brightest stars in our survey to measure directly dust extinction through the cloud following the method described by Lada et al. (1994). We spatially convolve the individual extinction measurements with a square filter 90" in size to construct a large-scale map of extinction in the cloud. We derive a total mass of M(L977)= (660 +/- 30)(D/500 pc)^2 Msun and, via a comparison of source counts with predictions of a galactic model, estimate a distance to L977 of 500 +/- 100 pc. We find a correlation between the measured dispersion in our extinction determinations and the extinction. We interpret this as evidence for the presence of structure on scales smaller than the 90" resolution of our extinction map. To further investigate the structure of the cloud we construct the frequency distribution of the 1628 individual extinction measurements in the L977 cloud. The shape of the distribution is similar to that of the IC 5146 cloud. Monte Carlo modeling of this distribution suggests that between 2 < Av < 40 mag (or roughly 1 < r < 0.1 pc) the material inside L977 is characterized by a density profile n(r) \propto r^(-2). Direct measurement of the radial profile of a portion of the cloud confirms this result. (more...)

We describe planetesimal accretion calculations in the Kuiper Belt. Our evolution code simulates planetesimal growth in a single annulus and includes velocity evolution but not fragmentation. Test results match analytic solutions and duplicate previous simulations at 1 AU. In the Kuiper Belt, simulations without velocity evolution produce a single runaway body with a radius of 1000 km on a time scale inversely proportional to the initial mass in the annulus. Runaway growth occurs in 100 Myr for 10 earth masses and an initial eccentricity of 0.001 in a 6 AU annulus centered at 35 AU. This mass is close to the amount of dusty material expected in a minimum mass solar nebula extrapolated into the Kuiper Belt. Simulations with velocity evolution produce runaway growth on a wide range of time scales. Dynamical friction and viscous stirring increase particle velocities in models with large (8 km radius) initial bodies. This velocity increase delays runaway growth by a factor of two compared to models without velocity evolution. In contrast, collisional damping dominates over dynamical friction and viscous stirring in models with small (80--800 m radius) initial bodies. Collisional damping decreases the time scale to runaway growth by factors of 4--10 relative to constant velocity calculations. Simulations with minimum mass solar nebulae, 10 earth masses, reach runaway growth on time scales of 20-40 Myr with 80 m initial bodies, 50-100 Myr with 800 m bodies, and 75-250 Myr for 8 km initial bodies. These growth times vary linearly with the mass of the annulus but are less sensitive to the initial eccentricity than constant velocity models.

This paper describes the first optical spectroscopic survey of class I sources (also known as embedded sources and protostars) in the Taurus-Auriga dark cloud. We detect 10 of the 24 known class I sources in the cloud at 5500-9000 A. All detected class I sources have strong H-alpha emission; most also have strong [O I] and [S II] emission. These data - together with high quality optical spectra of T Tauri stars in the Taurus-Auriga cloud - demonstrate that forbidden emission lines are stronger and more common in class I sources than in T Tauri stars. Our results also provide a clear discriminant in the frequency of forbidden line emission between weak-emission and classical T Tauri stars. In addition to strong emission lines, three class I sources have prominent TiO absorption bands. The M-type central stars of these sources mingle with optically visible T Tauri stars in the HR diagram and lie somewhat below both the birthline for spherical accretion and the deuterium burning sequence for disc accretion.

We describe a non-LTE photoionization code to calculate the wind structure and emergent spectrum of a red giant wind illuminated by the hot component of a symbiotic binary system. We consider spherically symmetric winds with several different velocity and temperature laws and derive predicted line fluxes as a function of the red giant mass loss rate, \mdot. Our models generally match observations of the symbiotic stars EG And and AG Peg for \mdot about 10^{-8} \msunyr to 10^{-7} \msunyr. The optically thick cross- section of the red giant wind as viewed from the hot component is a crucial parameter in these models. Winds with cross-sections of 2--3 red giant radii reproduce the observed fluxes, because the wind density is then high, about 10^9 cm^{-3}. Our models favor winds with acceleration regions that either lie far from the red giant photosphere or extend for 2--3 red giant radii.

We analyze the ultraviolet spectrum of the cataclysmic variable GK Per at maximum light. The flat ultraviolet spectrum in this system requires a truncated inner accretion disk and an unusually flat radial temperature profile. This requirement is not satisfied by any non-magnetic steady or non-steady disk model. We consider a magnetized accretion disk model to explain the ultraviolet spectrum. The available data on the white dwarf spin and possible quasi-periodic oscillations constrain the magnetic field, $B_{*}$, and the disk accretion rate, ${\dot M}$, to lie along a well-defined spin-equilibrium condition $({\dot M}/10^{17} ~ g ~ s^{-1}) \sim 100(B_{*}/10^7G)^2$. Our self-consistent treatment of the magnetic torque on the disk flattens the disk temperature distribution outside the disk truncation radius. This modified temperature distribution is too steep to explain the UV spectrum for reasonable field strengths. X-ray heating is a plausible alternative to magnetic heating in GK Per. We estimate that the disk intercepts $\sim$ 5% of the accretion energy in outburst, which results in an extra disk luminosity of $\sim$ 5--10 $L_{\odot}$. Model spectra of optically thick disks are too blue to match observations. The UV spectrum of an optically thick disk with an optically thin, X-ray heated corona resembles the observed spectrum. The X-ray luminosity observed during the outburst indicates ${\dot M}<10^{18} ~ g ~ s^{-1}$, which is a factor of 10 lower than that required to explain the ultraviolet luminosity. Radiation drag on material flowing inward along the accretion column lowers the shock temperature and reduces the X-ray luminosity. Most of the accretion energy is then radiated at extreme ultraviolet wavelengths.

We present solutions for the accretion disks and boundary layers in pre-main-sequence stars undergoing FU Orionis outbursts. These solutions differ from earlier disk solutions in that they include a self-consistent treatment of the boundary layer region. In a previous paper (Popham 1996), we showed that these stars should stop accreting angular momentum once they spin up to modest rotation rates. Here we show that for reasonable values of $\alpha$, these low angular momentum accretion rate solutions fit the spectra and line profiles observed in FU Orionis objects better than solutions with high rates of angular momentum accretion. We find solutions which fit the observations of FU Orionis and V1057 Cygni. These solutions have mass accretion rates of 2 and $1 \times 10^{-4} \msyr$, stellar masses of 0.7 and $0.5 \msun$, and stellar radii of 5.75 and $5.03 \rsun$, respectively. They also have modest stellar rotation rates $8-9 \xss$, comparable to the observed rotation rates of T Tauri stars, and angular momentum accretion rates of zero. This supports our earlier suggestion that FU Orionis outbursts may regulate the rotation rates of T Tauri stars.