More than a dozen young stars host spiral arms in their surrounding protoplanetary disks. The excitation mechanisms of such arms are under debate. The two leading hypotheses -- companion-disk interaction and gravitational instability (GI) -- predict distinct motion for spirals. By imaging the MWC 758 spiral arm system at two epochs spanning ${\sim}5$ yr using the SPHERE instrument on the Very Large Telescope (VLT), we test the two hypotheses for the first time. We find that the pattern speeds of the spirals are not consistent with the GI origin. Our measurements further evince the existence of a faint "missing planet" driving the disk arms. The average spiral pattern speed is $0.\!^\circ22\pm0.\!^\circ03$ yr$^{-1}$, pointing to a driver at $172_{-14}^{+18}$ au around a $1.9$ $M_\odot$ central star if it is on a circular orbit. In addition, we witness time varying shadowing effects on a global scale that are likely originated from an inner disk.

We report the results of a ${\sim}4$-year direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager Exoplanet Survey. We targeted nearby (${\lesssim}150$ pc), young (${\lesssim}500$ Myr) stars with high infrared excesses ($L_{\mathrm{IR}} / L_\star > 10^{-5}$), including 38 with previously resolved disks. Observations were made using the Gemini Planet Imager high-contrast integral field spectrograph in $H$-band (1.6 $\mu$m) coronagraphic polarimetry mode to measure both polarized and total intensities. We resolved 26 debris disks and three protoplanetary/transitional disks. Seven debris disks were resolved in scattered light for the first time, including newly presented HD 117214 and HD 156623, and we quantified basic morphologies of five of them using radiative transfer models. All of our detected debris disks but HD 156623 have dust-poor inner holes, and their scattered-light radii are generally larger than corresponding radii measured from resolved thermal emission and those inferred from spectral energy distributions. To assess sensitivity, we report contrasts and consider causes of non-detections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low inclination disks (${\lesssim} 70 \deg$). Based on post-survey statistics, we improved upon our pre-survey target prioritization metric predicting polarimetric disk detectability. We also examined scattered-light disks in the contexts of gas, far-IR, and millimeter detections. Comparing $H$-band and ALMA fluxes for two disks revealed tentative evidence for differing grain properties. Finally, we found no preference for debris disks to be detected in scattered light if wide-separation substellar companions were present.

Some Jupiter-mass exoplanets contain $\sim$$100\, M_\oplus$ of metals, well above the $\sim$$10\, M_\oplus$ typically needed in a solid core to trigger giant planet formation by runaway gas accretion. We demonstrate that such `heavy-metal Jupiters' can result from planetary mergers near $\sim$10 au. Multiple cores accreting gas at runaway rates gravitationally perturb one another onto crossing orbits such that the average merger rate equals the gas accretion rate. Concurrent mergers and gas accretion implies the core mass scales with the total planet mass as $M_{\rm core} \propto M^{1/5}$ - heavier planets harbour heavier cores, in agreement with the observed mass-metallicity relation. While the average gas giant merges about once to double its core, others may merge multiple times, as merger trees grow chaotically. We show that the dispersion of outcomes inherent in mergers can reproduce the large scatter in observed planet metallicities, assuming $3-30\, M_\oplus$ pre-runaway cores. Mergers potentially correlate metallicity, eccentricity, and spin.

Planets open gaps in discs. Gap opening is typically modeled by considering the planetary Lindblad torque which repels disc gas away from the planet's orbit. But gaps also clear because the planet consumes local material. We present a simple, easy-to-use, analytic framework for calculating how gaps deplete and how the disc's structure as a whole changes by the combined action of Lindblad repulsion and planetary consumption. The final mass to which a gap-embedded gas giant grows is derived in tandem. The analytics are tested against 1D numerical experiments and calibrated using published multi-dimensional simulations. In viscous alpha discs, the planet, while clearing a gap, initially accretes practically all of the gas that tries to diffuse past, rapidly achieving super-Jupiter if not brown dwarf status. By contrast, in inviscid discs---that may still accrete onto their central stars by, say, magnetized winds---planets open deep, repulsion-dominated gaps. Then only a small fraction of the disc accretion flow is diverted onto the planet, which grows to a fraction of a Jupiter mass. Transitional disc cavities might be cleared by families of such low-mass objects opening inviscid, repulsion-dominated, overlapping gaps which allow most of the outer disc gas to flow unimpeded onto host stars.

We present $Spitzer$ 4.5$\mu$m observations of the transit of TOI-700 d, a habitable zone Earth-sized planet in a multiplanet system transiting a nearby M-dwarf star (TIC 150428135, 2MASS J06282325-6534456). TOI-700 d has a radius of $1.144^{+0.062}_{-0.061}R_\oplus$ and orbits within its host star's conservative habitable zone with a period of 37.42 days ($T_\mathrm{eq} \sim 269$K). TOI-700 also hosts two small inner planets (R$_b$=$1.037^{+0.065}_{-0.064}R_\oplus$ & R$_c$=$2.65^{+0.16}_{-0.15}R_\oplus$) with periods of 9.98 and 16.05 days, respectively. Our $Spitzer$ observations confirm the TESS detection of TOI-700 d and remove any remaining doubt that it is a genuine planet. We analyze the $Spitzer$ light curve combined with the 11 sectors of TESS observations and a transit of TOI-700 c from the LCOGT network to determine the full system parameters. Although studying the atmosphere of TOI-700 d is not likely feasible with upcoming facilities, it may be possible to measure the mass of TOI-700 d using state-of-the-art radial velocity instruments (expected RV semi-amplitude of $\sim$70 cm/s).

We place the first constraints on the obliquity of a planetary-mass companion (PMC) outside of the Solar System. Our target is the directly imaged system 2MASS J01225093-2439505 (2M0122), which consists of a 120 Myr 0.4 M_sun star hosting a 12-27 M_J companion at 50 AU. We constrain all three of the system's angular momentum vectors: how the companion spin axis, the stellar spin axis, and the orbit normal are inclined relative to our line of sight. To accomplish this, we measure projected rotation rates (vsini) for both the star and the companion using new near-infrared high-resolution spectra with NIRSPEC at Keck Observatory. We combine these with a new stellar photometric rotation period from TESS and a published companion rotation period from HST to obtain spin axis inclinations for both objects. We also fitted multiple epochs of astrometry, including a new observation with NIRC2/Keck, to measure 2M0122b's orbital inclination. The three line-of-sight inclinations place limits on the true de-projected companion obliquity and stellar obliquity. We find that while the stellar obliquity marginally prefers alignment, the companion obliquity tentatively favors misalignment. We evaluate possible origin scenarios. While collisions, secular spin-orbit resonances, and Kozai-Lidov oscillations are unlikely, formation by gravitational instability in a gravito-turbulent disk - the scenario favored for brown dwarf companions to stars - appears promising.

The orbital period ratios of neighbouring sub-Neptunes are distributed asymmetrically near first-order resonances. There are deficits of systems---"troughs" in the period ratio histogram---just short of commensurability, and excesses---"peaks"---just wide of it. We reproduce quantitatively the strongest peak-trough asymmetries, near the 3:2 and 2:1 resonances, using dissipative interactions between planets and their natal discs. Disc eccentricity damping captures bodies into resonance and clears the trough, and when combined with disc-driven convergent migration, draws planets initially wide of commensurability into the peak. The migration implied by the magnitude of the peak is modest; reductions in orbital period are $\sim$10\%, supporting the view that sub-Neptunes complete their formation more-or-less in situ. Once captured into resonance, sub-Neptunes of typical mass $\sim$$5$--$15 M_{\oplus}$ stay captured (contrary to an earlier claim), as they are immune to the overstability that afflicts lower mass planets. Driving the limited, short-scale migration is a gas disc depleted in mass relative to a solar-composition disc by 3--5 orders of magnitude. Such gas-poor but not gas-empty environments are quantitatively consistent with sub-Neptune core formation by giant impacts (and not, e.g., pebble accretion). While disc-planet interactions at the close of the planet formation era adequately explain the 3:2 and 2:1 asymmetries at periods $\gtrsim$ $5$--$15$ days, subsequent modification by stellar tides appears necessary at shorter periods, particularly for the 2:1.

We present new $H$-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager (GPI). The disk is detected in total and polarized intensity down to a projected angular separation of 0.15", or 20au. On the other hand, the large scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We analyze the curvature of the disk spine and estimate a radius of $\approx$100au for the parent body belt, smaller than past scattered light studies but consistent with thermal emission maps of the system. We employ three different flux-preserving post-processing methods to suppress the residual starlight and evaluate the surface brightness and polarization profile along the disk spine. Unlike past studies of the system, our high fidelity images reveal the disk to be highly symmetric and devoid of morphological and surface brightness perturbations. We find the dust scattering properties of the system to be consistent with those observed in other debris disks, with the exception of HR 4796. Finally, we find no direct evidence for the presence of a planetary-mass object in the system.

Circumplanetary disks (CPDs) may be essential to the formation of planets, regulating their spin and accretion evolution. We perform a series of 3D hydrodynamics simulations in both the isothermal and adiabatic limits to systematically measure the rotation rates, sizes, and masses of CPDs as functions of $q_{\rm thermal}$, the ratio of the planet mass to the disk thermal mass. Our $q_{\rm thermal}$ ranges from 0.1 to 4; for our various disk temperatures, this corresponds to planet masses between 1 Earth mass and 4 Jupiter masses. Within this parameter space, we find that isothermal CPDs are disky and bound within $\sim$10\% of the planet's Bondi radius $r_{\rm B}$, with the innermost $\sim0.05\,r_{\rm B}$ in full rotational support. Adiabatic CPDs are spherical (therefore not actually "disks"), bound within $\sim0.2\,r_{\rm B}$, and mainly pressure-supported with rotation rates scaling linearly with $q_{\rm thermal}$; extrapolation suggests full rotational support of adiabatic envelopes at $\sim10\,q_{\rm thermal}$. Fast rotation and 3D supersonic flow render isothermal CPDs significantly different in structure from --- and orders of magnitude less massive than --- their 1D isothermal hydrostatic counterparts. Inside a minimum-mass solar nebula, even a maximally cooled, isothermal CPD around a 10 Earth-mass core may have less than 1 Earth mass, suggesting that gas giant formation may hinge on angular momentum transport processes in CPDs. Our CPD sizes and masses appear consistent with the regular satellites orbiting solar system giants.

During the runaway phase of their formation, gas giants fill their gravitational spheres of influence out to Bondi or Hill radii. When runaway ends, planets shrink several orders of magnitude in radius until they are comparable in size to present-day Jupiter; in 1D models, the contraction occurs on the Kelvin-Helmholtz time-scale $t_{\rm KH}$, which is initially a few thousand years. However, if angular momentum is conserved, contraction cannot complete, as planets are inevitably spun up to their breakup periods $P_{\rm break}$. We consider how a circumplanetary disc (CPD) can de-spin a primordially magnetized gas giant and remove the centrifugal barrier, provided the disc is hot enough to couple to the magnetic field, a condition that is easier to satisfy at later times. By inferring the planet's magnetic field from its convective cooling luminosity, we show that magnetic spin-down times are shorter than contraction times throughout post-runaway contraction: $t_{\rm mag}/t_{\rm KH}\sim(P_{\rm break}/t_{\rm KH})^{1/21}\lesssim 1$. Planets can spin down until they corotate with the CPD's magnetospheric truncation radius, at a period $P_{\rm max}/P_{\rm break} \sim (t_{\rm KH}/P_{\rm break})^{1/7}$. By the time the disc disperses, $P_{\rm max}/P_{\rm break}\sim 20-30$; further contraction at fixed angular momentum can spin planets back up to $\sim$$10 P_{\rm break}$, potentially explaining observed rotation periods of giant planets and brown dwarfs.

Giant planets are thought to form by runaway gas accretion onto solid cores. Growth must eventually stop running away, ostensibly because planets open gaps (annular cavities) in their surrounding discs. Typical models stop runaway by artificially capping the accretion rate and lowering it to zero over an arbitrarily short time-scale. In reality, post-runaway accretion persists as long as the disc remains. During this final and possibly longest phase of formation, when the planet is still emerging from the disc, its mass can more than double, and its radius contracts by orders of magnitude. By drawing from the theory of how gaps clear, we find that post-runaway accretion luminosities diverge depending on disc viscosity: luminosities fall in low-viscosity discs but continue to rise past runaway in high-viscosity discs. This divergence amounts to a factor of $10^2$ by the time the disc disperses. Irrespective of the specifics of how planets interact with discs, the observed luminosity and age of an accreting planet can be used to calculate its instantaneous mass, radius, and accretion rate. We perform this exercise for the planet candidates embedded within the discs orbiting PDS 70, HD 163296, and MWC 758, inferring masses of $1-10\, M_{\rm J}$, accretion rates of $0.1-10\, M_{\rm J}$/Myr, and radii of $1-10\, R_{\rm J}$. Our radii are computed self-consistently from the planet's concurrent contraction and accretion and do not necessarily equal the value of $2R_{\rm J}$ commonly assumed; in particular, the radius depends on the envelope opacity as $R \propto \kappa^{0.5}$.

Many debris disks seen in scattered light have shapes that imply their dust grains trace highly eccentric, apsidally aligned orbits. Apsidal alignment is surprising, especially for dust. Even when born from an apse-aligned ring of parent bodies, dust grains have their periastra dispersed in all directions by stellar radiation pressure. The periastra cannot be re-oriented by planets within the short dust lifetimes at the bottom of the collisional cascade. We propose that what re-aligns dust orbits is drag exerted by second-generation gas. Gas is largely immune to radiation pressure, and when released by photodesorption or collisions within an eccentric ring of parent bodies should occupy a similarly eccentric, apse-aligned ring. Dust grains launched onto misaligned orbits cross the eccentric gas ring supersonically and can become dragged into alignment within collisional lifetimes. The resultant dust configurations, viewed nearly but not exactly edge-on, with periastra pointing away from the observer, appear moth-like, with kinked wings and even doubled pairs of wings, explaining otherwise mysterious features in HD 61005 ("The Moth") and HD 32297, including their central bulbs when we account for strong forward scattering from irregularly shaped particles. Around these systems we predict gas at Kuiper-belt-like distances to move on highly elliptical streamlines that owe their elongation, ultimately, to highly eccentric planets. Unresolved issues and an alternative explanation for apsidal alignment are outlined.

All four circumbinary (CB) protoplanetary disks orbiting short-period ($P < 20$ day) double-lined spectroscopic binaries (SB2s)---a group that includes UZ Tau E, for which we present new ALMA data---exhibit sky-plane inclinations $i_{\rm disk}$ which match, to within a few degrees, the sky-plane inclinations $i_\star$ of their stellar hosts. Although for these systems the true mutual inclinations $\theta$ between disk and binary cannot be directly measured because relative nodal angles are unknown, the near-coincidence of $i_{\rm disk}$ and $i_\star$ suggests that $\theta$ is small for these most compact of systems. We confirm this hypothesis using a hierarchical Bayesian analysis, showing that 68% of CB disks around short-period SB2s have $\theta < 3.0^\circ$. Near co-planarity of CB disks implies near co-planarity of CB planets discovered by Kepler, which in turn implies that the occurrence rate of close-in CB planets is similar to that around single stars. By contrast, at longer periods ranging from $30-10^5$ days (where the nodal degeneracy can be broken via, e.g., binary astrometry), CB disks exhibit a wide range of mutual inclinations, from co-planar to polar. Many of these long-period binaries are eccentric, as their component stars are too far separated to be tidally circularized. We discuss how theories of binary formation and disk-binary gravitational interactions can accommodate all these observations.

Gas giants are thought to form by runaway accretion: an instability driven by the self-gravity of growing atmospheres that causes accretion rates to rise super-linearly with planet mass. Why runaway should stop at a Jupiter or any other mass is unknown. We consider the proposal that final masses are controlled by circumstellar disc gaps (cavities) opened by planetary gravitational torques. We develop a fully time-dependent theory of gap formation and couple it self-consistently to planetary growth rates. When gaps first open, planetary torques overwhelm viscous torques, and gas depletes as if it were inviscid. In low-viscosity discs, of the kind motivated by recent observations and theory, gaps stay predominantly in this inviscid phase and planet masses finalize at $M_{\rm final}/M_\star\sim(\Omega t_{\rm disc})^{0.07}(H/a)^{2.73}(G\rho_0/\Omega^2)^{1/3}$, with $M_\star$ the host stellar mass, $\Omega$ the planet's orbital angular velocity, $t_{\rm disc}$ the gas disc's lifetime, $H/a$ its aspect ratio, and $\rho_0$ its unperturbed density. This final mass is independent of the dimensionless viscosity $\alpha$ and applies to large orbital distances, typically beyond $\sim$10 AU, where disc scale heights exceed planet radii. It evaluates to a few Jupiter masses at 10-100 AU, increasing gradually with distance as gaps become harder to open.

We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^\alpha a^\beta$, finding $\alpha$ = -2.4 $\pm$ 0.8 and $\beta$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability.

The vertical distribution of dust in debris disks is sensitive to the number and size of large planetesimals dynamically stirring the disk, and is therefore well-suited for constraining the prevalence of otherwise unobservable Uranus and Neptune analogs. Information regarding stirring bodies has previously been inferred from infrared and optical observations of debris disk vertical structure, but theoretical works predict that the small particles traced by short-wavelength observations will be `puffed up' by radiation pressure, yielding only upper limits. The large grains that dominate the disk emission at millimeter wavelengths are much less sensitive to the effects of stellar radiation or stellar winds, and therefore trace the underlying mass distribution more directly. Here we present ALMA 1.3 mm dust continuum observations of the debris disk around the nearby M star AU Mic. The 3 au spatial resolution of the observations, combined with the favorable edge-on geometry of the system, allows us to measure the vertical thickness of the disk. We report a scale height-to-radius aspect ratio of $h = 0.031_{-0.004}^{+0.005}$ between radii of $\sim 23$ au and $\sim 41$ au. Comparing this aspect ratio to a theoretical model of size-dependent velocity distributions in the collisional cascade, we find that the perturbing bodies embedded in the local disk must be larger than about 400 km, and the largest perturbing body must be smaller than roughly 1.8 M$_\odot$. These measurements rule out the presence of a gas giant or Neptune analog near the $\sim 40$ au outer edge of the debris ring, but are suggestive of large planetesimals or an Earth-sized planet stirring the dust distribution.

ALMA has found multiple dust gaps and rings in a number of protoplanetary disks in continuum emission at millimeter wavelengths. The origin of such structures is in debate. Recently, we documented how one super-Earth planet can open multiple (up to five) dust gaps in a disk with low viscosity ($\alpha\lesssim10^{-4}$). In this paper, we examine how the positions, depths, and total number of gaps opened by one planet depend on input parameters, and apply our results to real systems. Gap locations (equivalently, spacings) are the easiest metric to use when making comparisons between theory and observations, as positions can be robustly measured. We fit the locations of gaps empirically as functions of planet mass and disk aspect ratio. We find that the locations of the double gaps in HL Tau and TW Hya, and of all three gaps in HD 163296, are consistent with being opened by a sub-Saturn mass planet. This scenario predicts the locations of other gaps in HL Tau and TW Hya, some of which appear consistent with current observations. We also show how the Rossby wave instability may develop at the edges of several gaps and result in multiple dusty vortices, all caused by one planet. A planet as low in mass as Mars may produce multiple dust gaps in the terrestrial planet forming region.

The HR 8799 system uniquely harbors four young super-Jupiters whose orbits can provide insights into the system's dynamical history and constrain the masses of the planets themselves. Using the Gemini Planet Imager (GPI), we obtained down to one milliarcsecond precision on the astrometry of these planets. We assessed four-planet orbit models with different levels of constraints and found that assuming the planets are near 1:2:4:8 period commensurabilities, or are coplanar, does not worsen the fit. We added the prior that the planets must have been stable for the age of the system (40 Myr) by running orbit configurations from our posteriors through $N$-body simulations and varying the masses of the planets. We found that only assuming the planets are both coplanar and near 1:2:4:8 period commensurabilities produces dynamically stable orbits in large quantities. Our posterior of stable coplanar orbits tightly constrains the planets' orbits, and we discuss implications for the outermost planet b shaping the debris disk. A four-planet resonance lock is not necessary for stability up to now. However, planet pairs d and e, and c and d, are each likely locked in two-body resonances for stability if their component masses are above $6~M_{\rm{Jup}}$ and $7~M_{\rm{Jup}}$, respectively. Combining the dynamical and luminosity constraints on the masses using hot-start evolutionary models and a system age of $42 \pm 5$~Myr, we found the mass of planet b to be $5.8 \pm 0.5~M_{\rm{Jup}}$, and the masses of planets c, d, and e to be $7.2_{-0.7}^{+0.6}~M_{\rm{Jup}}$ each.

Pebble accretion refers to the assembly of rocky planet cores from particles whose velocity dispersions are damped by drag from circumstellar disc gas. Accretion cross-sections can approach maximal Hill-sphere scales for particles whose Stokes numbers approach unity. While fast, pebble accretion is also lossy. Gas drag brings pebbles to protocores but also sweeps them past; those particles with the largest accretion cross-sections also have the fastest radial drift speeds and are the most easily drained out of discs. We present a global model of planet formation by pebble accretion that keeps track of the disc's mass budget. Cores, each initialized with a lunar mass, grow from discs whose finite stores of mm-cm sized pebbles drift inward across all radii in viscously accreting gas. For every 1 $M_\oplus$ netted by a core, at least 10 $M_\oplus$ and possibly much more are lost to radial drift. Core growth rates are typically exponentially sensitive to particle Stokes number, turbulent Mach number, and solid surface density. This exponential sensitivity, when combined with disc migration, tends to generate binary outcomes from 0.1-30 AU: either sub-Earth cores remain sub-Earth, or explode into Jupiters, with the latter migrating inward to varying degrees. When Jupiter-breeding cores assemble from mm-cm sized pebbles, they do so in discs where such particles drain out in $\sim$10$^5$ yr or less; such fast-draining discs do not fit mm-wave observations.

Large-scale spiral arms have been revealed in scattered light images of a few protoplanetary disks. Theoretical models suggest that such arms may be driven by and co-rotate with giant planets, which has called for remarkable observational efforts to look for them. By examining the rotation of the spiral arms for the MWC 758 system over a 10-yr timescale, we are able to provide dynamical constraints on the locations of their perturbers. We present reprocessed Hubble Space Telescope (HST)/NICMOS F110W observations of the target in 2005, and the new Keck/NIRC2 $L'$-band observations in 2017. MWC~758's two well-known spiral arms are revealed in the NICMOS archive at the earliest observational epoch. With additional Very Large Telescope (VLT)/SPHERE data, our joint analysis leads to a pattern speed of ${0.6^\circ}^{+3.3^\circ}_{-0.6^\circ}\, \mathrm{yr}^{-1}$ at $3\sigma$ for the two major spiral arms. If the two arms are induced by a perturber on a near-circular orbit, its best fit orbit is at $89$ au ($0.59"$), with a $3\sigma$ lower limit of 30 au ($0.20"$). This finding is consistent with the simulation prediction of the location of an arm-driving planet for the two major arms in the system.

We present new high resolution imaging of a light-scattering dust ring and halo around the young star HD 35841. Using spectroscopic and polarimetric data from the Gemini Planet Imager in H-band (1.6 microns), we detect the highly inclined (i=85 deg) ring of debris down to a projected separation of ~12 au (~0.12") for the first time. Optical imaging from HST/STIS shows a smooth dust halo extending outward from the ring to >140 au (>1.4"). We measure the ring's scattering phase function and polarization fraction over scattering angles of 22-125 deg, showing a preference for forward scattering and a polarization fraction that peaks at ~30% near the ansae. Modeling of the scattered-light disk indicates that the ring spans radii of ~60-220 au, has a vertical thickness similar to that of other resolved dust rings, and contains grains as small as 1.5 microns in diameter. These models also suggest the grains have a low porosity, are more likely to consist of carbon than astrosilicates, and contain significant water ice. The halo has a surface brightness profile consistent with that expected from grains pushed by radiation pressure from the main ring onto highly eccentric but still bound orbits. We also briefly investigate arrangements of a possible inner disk component implied by our spectral energy distribution models, and speculate about the limitations of Mie theory for doing detailed analyses of debris disk dust populations.

Spatially resolved structures in protoplanetary disks hint at unseen planets. Previous imaging observations of the transitional disk around MWC 758 revealed an inner cavity, a ring-like outer disk, emission clumps, and spiral arms, all possibly generated by companions. We present ALMA dust continuum observations of MWC 758 at 0.87 millimeter (mm) wavelength with 43$\times$39 mas angular resolution (6.9$\times$6.2 AU) and 20 $\mu$Jy beam$^{-1}$ rms. The central sub-mm emission cavity is revealed to be eccentric; once deprojected, its outer edge can be well-fitted by an ellipse with an eccentricity of 0.1 and one focus on the star. The broad ring-like outer disk is resolved into three narrow rings with two gaps in between. The outer two rings tentatively show the same eccentricity and orientation as the innermost ring bounding the inner cavity. The two previously known dust emission clumps are resolved in both the radial and azimuthal directions, with radial widths equal to $\sim$4$\times$ the local scale height. Only one of the two spiral arms previously imaged in near-infrared (NIR) scattered light is revealed in ALMA dust emission, at a slightly larger stellocentric distance owing to projection effects. We also submit evidence of disk truncation at $\sim$100 AU based on comparing NIR imaging observations with models. The spirals, the north clump, and the truncated disk edge are all broadly consistent with the presence of one companion exterior to the spirals at roughly 100 AU.

Models of core accretion assume that in the radiative zones of accreting gas envelopes, radiation diffuses. But super-Earths/sub-Neptunes (1-4$R_\oplus$, 2-20$M_\oplus$) point to formation conditions that are optically thin: their modest gas masses are accreted from short-lived and gas-poor nebulae reminiscent of the transparent cavities of transitional disks. Planetary atmospheres born in such environments can be optically thin to both incident starlight and internally generated thermal radiation. We construct time-dependent models of such atmospheres, showing that super-Earths/sub-Neptunes can accrete their $\sim$1%-by-mass gas envelopes, and super-puffs/sub-Saturns their $\sim$20%-by-mass envelopes, over a wide range of nebular depletion histories requiring no fine tuning. Although nascent atmospheres can exhibit stratospheric temperature inversions effected by atomic Fe and various oxides that absorb strongly at visible wavelengths, the rate of gas accretion remains controlled by the radiative-convective boundary (rcb) at much greater pressures. For dusty envelopes, the temperature at the rcb $T_{\rm rcb} \simeq 2500$ K is still set by ${\rm H}_2$ dissociation; for dust-depleted envelopes, $T_{\rm rcb}$ tracks the temperature of the visible or thermal photosphere, whichever is deeper, out to at least $\sim$5 AU. The rate of envelope growth remains largely unchanged between the old radiative diffusion models and the new optically thin models, reinforcing how robustly super-Earths form as part of the endgame chapter in disk evolution.

The behavior of an interior test particle in the secular 3-body problem has been studied extensively. A well-known feature is the Lidov-Kozai resonance in which the test particle's argument of periapse librates about $\pm 90^\circ$ and large oscillations in eccentricity and inclination are possible. Less explored is the inverse problem: the dynamics of an exterior test particle and an interior perturber. We survey numerically the inverse secular problem, expanding the potential to hexadecapolar order and correcting an error in the published expansion. Four secular resonances are uncovered that persist in full $N$-body treatments (in what follows, $\varpi$ and $\Omega$ are the longitudes of periapse and of ascending node, $\omega$ is the argument of periapse, and subscripts 1 and 2 refer to the inner perturber and outer test particle): (i) an orbit-flipping quadrupole resonance requiring a non-zero perturber eccentricity $e_1$, in which $\Omega_2-\varpi_1$ librates about $\pm 90^\circ$; (ii) a hexadecapolar resonance (the "inverse Kozai" resonance) for perturbers that are circular or nearly so and inclined by $I \simeq 63^\circ/117^\circ$, in which $\omega_2$ librates about $\pm 90^\circ$ and which can vary the particle eccentricity by $\Delta e_2 \simeq 0.2$ and lead to orbit crossing; (iii) an octopole "apse-aligned" resonance at $I \simeq 46^\circ/107^\circ$ wherein $\varpi_2 - \varpi_1$ librates about $0^\circ$ and $\Delta e_2$ grows with $e_1$; and (iv) an octopole resonance at $I \simeq 73^\circ/134^\circ$ wherein $\varpi_2 + \varpi_1 - 2 \Omega_2$ librates about $0^\circ$ and $\Delta e_2$ can be as large as 0.3 for small $e_1 \neq 0$. The more eccentric the perturber, the more the particle's eccentricity and inclination vary; also, more polar orbits are more chaotic. Our inverse solutions may be applied to the Kuiper belt and debris disks, circumbinary planets, and stellar systems.

We explain the fast-moving, ripple-like features in the edge-on debris disk orbiting the young M dwarf AU Mic. The bright features are clouds of sub-micron dust repelled by the host star's wind. The clouds are produced by avalanches: radial outflows of dust that gain exponentially more mass as they shatter background disk particles in collisional chain reactions. The avalanches are triggered from a region a few AU across -- the "avalanche zone" -- located on AU Mic's primary "birth" ring, at a true distance of $\sim$35 AU from the star but at a projected distance more than a factor of 10 smaller: the avalanche zone sits directly along the line of sight to the star, on the side of the ring nearest Earth, launching clouds that disk rotation sends wholly to the southeast, as observed. The avalanche zone marks where the primary ring intersects a secondary ring of debris left by the catastrophic disruption of a progenitor up to Varuna in size, less than tens of thousands of years ago. Only where the rings intersect are particle collisions sufficiently violent to spawn the sub-micron dust needed to seed the avalanches. We show that this picture works quantitatively, reproducing the masses, sizes, and velocities of the observed escaping clouds. The Lorentz force exerted by the wind's magnetic field, whose polarity reverses periodically according to the stellar magnetic cycle, promises to explain the observed vertical undulations. The timescale between avalanches, about 10 yr, might be set by time variability of the wind mass-loss rate or, more speculatively, by some self-regulating limit cycle.

Gas kinematics are an important part of the planet formation process. Turbulence influences planetesimal growth and migration from the scale of sub-micron dust grains through gas-giant planets. Radio observations of resolved molecular line emission can directly measure this non-thermal motion and, taking advantage of the layered chemical structure of disks, different molecular lines can be combined to map the turbulence throughout the vertical extent of a protoplanetary disk. Here we present ALMA observations of three molecules (DCO$^+$(3-2), C$^{18}$O(2-1) and CO(2-1)) from the disk around HD 163296. We are able to place stringent upper limits ($v_{\rm turb}<$0.06c$_s$, $<$0.05c$_s$ and $<$0.04c$_s$ for CO(2-1), C$^{18}$O(2-1) and DCO$^+$(3-2) respectively), corresponding to $\alpha\lesssim$3$\times$10$^{-3}$, similar to our prior limit derived from CO(3-2). This indicates that there is little turbulence throughout the vertical extent of the disk, contrary to theoretical predictions based on the magneto-rotational instability and gravito-turbulence. In modeling the DCO$^+$ emission we also find that it is confined to three concentric rings at 65.7$\pm$0.9 au, 149.9$^{+0.5}_{-0.7}$ au and 259$\pm$1 au, indicative of a complex chemical environment.

We investigate the observational signatures of super-Earths (i.e., Earth-to-Neptune mass planets) in their natal disks of gas and dust. Combining two-fluid global hydrodynamics simulations with a radiative transfer code, we calculate the distributions of gas and of sub-mm-sized dust in a disk perturbed by a super-Earth, synthesizing images in near-infrared scattered light and the mm-wave thermal continuum for direct comparison with observations. In low viscosity gas ($\alpha\lesssim10^{-4}$), a super-Earth opens two annular gaps to either side of its orbit by the action of Lindblad torques. This double gap and its associated gas pressure gradients cause dust particles to be dragged by gas into three rings: one ring sandwiched between the two gaps, and two rings located at the gap edges farthest from the planet. Depending on system parameters, additional rings may manifest for a single planet. A double gap located at tens of AUs from a host star in Taurus can be detected in the dust continuum by the Atacama Large Millimeter Array (ALMA) at an angular resolution of ~0".03 after two hours of integration. Ring and gap features persist in a variety of background disk profiles, last for thousands of orbits, and change their relative positions and dimensions depending on the speed and direction of planet migration. Candidate double gaps have been observed by ALMA in systems like HL Tau (D5 and D6) and TW Hya (at 37 and 43 AU); we submit that each double gap is carved by one super-Earth in nearly inviscid gas.

Sub-Neptunes around FGKM dwarfs are evenly distributed in log orbital period down to $\sim$10 days, but dwindle in number at shorter periods. Both the break at $\sim$10 days and the slope of the occurrence rate down to $\sim$1 day can be attributed to the truncation of protoplanetary disks by their host star magnetospheres at co-rotation. We demonstrate this by deriving planet occurrence rate profiles from empirical distributions of pre-main-sequence stellar rotation periods. Observed profiles are better reproduced when planets are distributed randomly in disks---as might be expected if planets formed in situ---rather than piled up near disk edges, as would be the case if they migrated in by disk torques. Planets can be brought from disk edges to ultra-short ($< 1$ day) periods by asynchronous equilibrium tides raised on their stars. Tidal migration can account for how ultra-short period planets (USPs) are more widely spaced than their longer period counterparts. Our picture provides a starting point for understanding why metal-rich stars tend to harbor more USPs, and why the sub-Neptune population drops at $\sim$10 days regardless of whether the host star is of type FGK or early M. We predict planet occurrence rates around A stars to also break at short periods, but at $\sim$1 day instead of $\sim$10 days because A stars rotate faster than lower mass stars (this prediction presumes that the planetesimal building blocks of planets can drift inside the dust sublimation radius).

Two longstanding problems in planet formation include (1) understanding how planets survive migration, and (2) articulating the process by which protoplanetary disks disperse---and in particular how they accrete onto their central stars. We can go a long way toward solving both problems if the disk gas surrounding planets has no intrinsic diffusivity ("viscosity"). In inviscid, laminar disks, a planet readily repels gas away from its orbit. On short timescales, zero viscosity gas accumulates inside a planet's orbit to slow Type I migration by orders of magnitude. On longer timescales, multiple super-Earths (distributed between, say, $\sim$0.1--10 AU) can torque inviscid gas out of interplanetary space, either inward to feed their stars, or outward to be blown away in a wind. We explore this picture with 2D hydrodynamics simulations of Earths and super-Earths embedded in inviscid disks, confirming their slow/stalled migration even under gas-rich conditions, and showing that disk transport rates range up to $\sim$$10^{-7} M_\odot~{\rm yr^{-1}}$ and scale as $\dot{M} \propto \Sigma M_{\rm p}^{3/2}$, where $\Sigma$ is the disk surface density and $M_{\rm p}$ is the planet mass. Gas initially sandwiched between two planets is torqued past both into the inner and outer disks. In sum, sufficiently compact systems of super-Earths can clear their natal disk gas, in a dispersal history that may be complicated and non-steady, but which conceivably leads over Myr timescales to large gas depletions similar to those characterizing transition disks.

We report ALMA Cycle 2 observations of 230 GHz (1.3 mm) dust continuum emission, and $^{12}$CO, $^{13}$CO, and C$^{18}$O J = 2-1 line emission, from the Upper Scorpius transitional disk [PZ99] J160421.7-213028, with an angular resolution of ~0".25 (35 AU). Armed with these data and existing H-band scattered light observations, we measure the size and depth of the disk's central cavity, and the sharpness of its outer edge, in three components: sub-$\mu$m-sized "small" dust traced by scattered light, millimeter-sized "big" dust traced by the millimeter continuum, and gas traced by line emission. Both dust populations feature a cavity of radius $\sim$70 AU that is depleted by factors of at least 1000 relative to the dust density just outside. The millimeter continuum data are well explained by a cavity with a sharp edge. Scattered light observations can be fitted with a cavity in small dust that has either a sharp edge at 60 AU, or an edge that transitions smoothly over an annular width of 10 AU near 60 AU. In gas, the data are consistent with a cavity that is smaller, about 15 AU in radius, and whose surface density at 15 AU is $10^{3\pm1}$ times smaller than the surface density at 70 AU; the gas density grades smoothly between these two radii. The CO isotopologue observations rule out a sharp drop in gas surface density at 30 AU or a double-drop model as found by previous modeling. Future observations are needed to assess the nature of these gas and dust cavities, e.g., whether they are opened by multiple as-yet-unseen planets or photoevaporation.

Giant planets can clear deep gaps when embedded in 2D (razor-thin) viscous circumstellar disks. We show by direct simulation that giant planets are just as capable of carving out gaps in 3D. Surface density maps are similar between 2D and 3D, even in detail. In particular, the scaling $\Sigma_{\rm gap} \propto q^{-2}$ of gap surface density with planet mass, derived from a global "zero-dimensional" balance of Lindblad and viscous torques, applies equally well to results obtained at higher dimensions. Our 3D simulations reveal extensive, near-sonic, meridional flows both inside and outside the gaps; these large-scale circulations might bear on disk compositional gradients, in dust or other chemical species. At high planet mass, gap edges are mildly Rayleigh unstable and intermittently shed streams of material into the gap - less so in 3D than in 2D.

The HD 61005 debris disk ("The Moth") stands out from the growing collection of spatially resolved circumstellar disks by virtue of its unusual swept-back morphology, brightness asymmetries, and dust ring offset. Despite several suggestions for the physical mechanisms creating these features, no definitive answer has been found. In this work, we demonstrate the plausibility of a scenario in which the disk material is shaped dynamically by an eccentric, inclined planet. We present new Keck NIRC2 scattered-light angular differential imaging of the disk at 1.2-2.3 microns that further constrains its outer morphology (projected separations of 27-135 AU). We also present complementary Gemini Planet Imager 1.6 micron total intensity and polarized light detections that probe down to projected separations less than 10 AU. To test our planet-sculpting hypothesis, we employed secular perturbation theory to construct parent body and dust distributions that informed scattered-light models. We found that this method produced models with morphological and photometric features similar to those seen in the data, supporting the premise of a planet-perturbed disk. Briefly, our results indicate a disk parent body population with a semimajor axis of 40-52 AU and an interior planet with an eccentricity of at least 0.2. Many permutations of planet mass and semimajor axis are allowed, ranging from an Earth mass at 35 AU to a Jupiter mass at 5 AU.

A "minimum model" for debris disks consists of a narrow ring of parent bodies, secularly forced by a single planet on a possibly eccentric orbit, colliding to produce dust grains that are perturbed by stellar radiation pressure. We demonstrate how this minimum model can reproduce a wide variety of disk morphologies imaged in scattered starlight. Five broad categories of disk shape can be captured: "rings," "needles," "ships-and-wakes," "bars," and "moths (a.k.a. fans)," depending on the viewing geometry. Moths can also sport "double wings." We explain the origin of morphological features from first principles, exploring the dependence on planet eccentricity, disk inclination dispersion, and the parent body orbital phases at which dust grains are born. A key determinant in disk appearance is the degree to which dust grain orbits are apsidally aligned. Our study of a simple steady-state (secularly relaxed) disk should serve as a reference for more detailed models tailored to individual systems. We use the intuition gained from our guidebook of disk morphologies to interpret, informally, the images of a number of real-world debris disks. These interpretations suggest that the farthest reaches of planetary systems are perturbed by eccentric planets, possibly just a few Earth masses each.

We report the discovery and confirmation of two sub-Saturn planets orbiting a bright (V = 11.3), metal-rich ([Fe/H] = 0.42 $\pm$ 0.04 dex) G3 dwarf in the K2 Campaign 2 field. The planets are 5.68 $\pm$ 0.56 Earth-radii and 7.82 $\pm$ 0.72 Earth-radii and have orbital periods of 20.8851 $\pm$ 0.0003 d and 42.3633$\pm$0.0006 d, near to the 2:1 mean-motion resonance. We obtained 32 radial velocities (RVs) with Keck/HIRES and detected the reflex motion due to EPIC-203771098b and c. These planets have masses of 21.0 $\pm$ 5.4 Earth-masses and 27.0 $\pm$ 6.9 Earth-masses, respectively. With low densities of 0.63 $\pm$ 0.25 g/cc and 0.31 $\pm$ 0.12 g/cc, respectively, the planets require thick envelopes of H/He to explain their large sizes and low masses. Interior structure models predict that the planets have fairly massive cores of 17.6 $\pm$ 4.3 Earth-masses and 16.1 $\pm$ 4.2 Earth-masses, respectively. They may have formed exterior to their present locations, accreted their H/He envelopes at large orbital distances, and migrated in as a resonant pair. The proximity to resonance, large transit depths, and host star brightness offer rich opportunities for TTV follow-up. Finally, the low surface gravities of the EPIC-203771098 planets make them favorable targets for transmission spectroscopy by HST, Spitzer, and JWST.

Direct imaging observations of protoplanetary disks at near-infrared (NIR) wavelengths have revealed structures of potentially planetary origin. Investigations of observational signatures from planet-induced features have so far focused on disks viewed face-on. Combining 3D hydrodynamics and radiative transfer simulations, we study how the appearance of the spiral arms and the gap produced in a disk by a companion varies with inclination and position angle in NIR scattered light. We compare the cases of a $3M_{\rm J}$ and a $0.1M_{\odot}$ companion, and make predictions suitable for testing with Gemini/GPI, VLT/NACO/SPHERE, and Subaru/HiCIAO/SCExAO. We find that the two trailing arms produced by an external perturber can have a variety of morphologies in inclined systems - they may appear as one trailing arm; two trailing arms on the same side of the disk; or two arms winding in opposite directions. The disk ring outside a planetary gap may also mimic spiral arms when viewed at high inclinations. We suggest potential explanations for the features observed in HH 30, HD 141569 A, AK Sco, HD 100546, and AB Aur. We emphasize that inclined views of companion-induced features cannot be converted into face-on views using simple and commonly practiced image deprojections.

Protostellar (class 0/I) disks, having masses comparable to those of their nascent host stars, and fed continuously from their natal infalling envelopes, are prone to gravitational instability (GI). Motivated by advances in near-infrared (NIR) adaptive optics imaging and mm-wave interferometry, we explore the observational signatures of GI in disks, using hydrodynamical and Monte Carlo radiative transfer simulations to synthesize NIR scattered light images and mm dust continuum maps. Spiral arms induced by GI, located at disk radii of hundreds of AUs, are local overdensities and have their photospheres displaced to higher altitudes above the disk midplane, arms therefore scatter more NIR light from their central stars than inter-arm regions, and are detectable at distances up to 1 kpc by Gemini/GPI, VLT/SPHERE, and Subaru/HiCIAO/SCExAO. By contrast, collapsed clumps formed by disk fragmentation have such strong local gravitational fields that their scattering photospheres are at lower altitudes, such fragments appear fainter than their surroundings in NIR scattered light. Spiral arms and streamers recently imaged in four FU Ori systems at NIR wavelengths resemble GI-induced structures and support the interpretation that FUors are gravitationally unstable protostellar disks. At mm wavelengths, both spirals and clumps appear brighter in thermal emission than the ambient disk and can be detected by ALMA at distances up to 0.4 kpc with one-hour integration times at ~0.1 arcsec resolution. Collapsed fragments having masses >=1 MJ can be detected by ALMA within ~10 minutes.

The giant impact phase of terrestrial planet formation establishes connections between super-Earths' orbital properties (semimajor axis spacings, eccentricities, mutual inclinations) and interior compositions (the presence or absence of gaseous envelopes). Using N-body simulations and analytic arguments, we show that spacings derive not only from eccentricities, but also from inclinations. Flatter systems attain tighter spacings, a consequence of an eccentricity equilibrium between gravitational scatterings, which increase eccentricities, and mergers, which damp them. Dynamical friction by residual disk gas plays a critical role in regulating mergers and in damping inclinations and eccentricities. Systems with moderate gas damping and high solid surface density spawn gas-enveloped super-Earths with tight spacings, small eccentricities, and small inclinations. Systems in which super-Earths coagulate without as much ambient gas, in disks with low solid surface density, produce rocky planets with wider spacings, larger eccentricities, and larger mutual inclinations. A combination of both populations can reproduce the observed distributions of spacings, period ratios, transiting planet multiplicities, and transit duration ratios exhibited by Kepler super-Earths. The two populations, both formed in situ, also help to explain observed trends of eccentricity vs. planet size, and bulk density vs. method of mass measurement (radial velocities vs. transit timing variations). Simplifications made in this study --- including the limited timespan of the simulations, and the approximate treatments of gas dynamical friction and gas depletion history --- should be improved upon in future work to enable a detailed quantitative comparison to the observations.

The dynamics of solid bodies in protoplanetary disks are subject to the properties of any underlying gas turbulence. Turbulence driven by disk self-gravity shows features distinct from those driven by the magnetorotational instability (MRI). We study the dynamics of solids in gravito-turbulent disks with two-dimensional (in the disk plane), hybrid (particle and gas) simulations. Gravito-turbulent disks can exhibit stronger gravitational stirring than MRI-active disks, resulting in greater radial diffusion and larger eccentricities and relative speeds for large particles (those with dimensionless stopping times $t_{stop} \Omega > 1$, where $\Omega$ is the orbital frequency). The agglomeration of large particles into planetesimals by pairwise collisions is therefore disfavored in gravito-turbulent disks. However, the relative speeds of intermediate-size particles $t_{stop} \Omega \sim 1$ are significantly reduced as such particles are collected by gas drag and gas gravity into coherent filament-like structures with densities high enough to trigger gravitational collapse. First-generation planetesimals may form via gravitational instability of dust in marginally gravitationally unstable gas disks.

Recent VLT/SPHERE near-infrared imaging observations revealed two spiral arms with a near m=2 rotational symmetry in the protoplanetary disk around the ~1.7 solar mass Herbig star HD 100453. A ~0.3 solar mass M dwarf companion, HD 100453 B, was also identified at a projected separation of 120 AU from the primary. In this Letter, we carry out hydrodynamic and radiative transfer simulations to examine the scattered light morphology of the HD 100453 disk as perturbed by the companion on a circular and coplanar orbit. We find that the companion truncates the disk at ~45 AU in scattered light images, and excites two spiral arms in the remaining (circumprimary) disk with a near m=2 rotational symmetry. Both the truncated disk size and the morphology of the spirals are in excellent agreement with the SPHERE observations at Y , J, H, and K1-bands, suggesting that the M dwarf companion is indeed responsible for the observed double-spiral-arm pattern. Our model suggests that the disk is close to face on (inclination angle ~5 degree), and that the entire disk-companion system rotates counterclockwise on the sky. The HD 100453 observations, along with our modeling work, demonstrate that double spiral arm patterns in near-infrared scattered light images can be generically produced by companions, and support future observations to identify the companions responsible for the arms observed in the MWC 758 and SAO 206462 systems.

The riddle posed by super-Earths (1-4$R_\oplus$, 2-20$M_\oplus$) is that they are not Jupiters: their core masses are large enough to trigger runaway gas accretion, yet somehow super-Earths accreted atmospheres that weigh only a few percent of their total mass. We show that this puzzle is solved if super-Earths formed late, as the last vestiges of their parent gas disks were about to clear. This scenario would seem to present fine-tuning problems, but we show that there are none. Ambient gas densities can span many (up to 9) orders of magnitude, and super-Earths can still robustly emerge after $\sim$0.1-1 Myr with percent-by-weight atmospheres. Super-Earth cores are naturally bred in gas-poor environments where gas dynamical friction has weakened sufficiently to allow constituent protocores to merge. So little gas is present at the time of core assembly that cores hardly migrate by disk torques: formation of super-Earths can be in situ. The picture --- that close-in super-Earths form in a gas-poor (but not gas-empty) inner disk, fed continuously by gas that bleeds inward from a more massive outer disk --- recalls the largely evacuated but still accreting inner cavities of transitional protoplanetary disks. We also address the inverse problem presented by super-puffs: an uncommon class of short-period planets seemingly too voluminous for their small masses (4-10$R_\oplus$, 2-6$M_\oplus$). Super-puffs easily acquire their thick atmospheres as dust-free, rapidly cooling worlds outside $\sim$1 AU where nebular gas is colder, less dense, and therefore less opaque. Unlike super-Earths which can form in situ, super-puffs migrated in to their current orbits; they are expected to form the outer links of mean-motion resonant chains, and to exhibit greater water content. We close by confronting observations and itemizing remaining questions.

The presence of debris disks around young main sequence stars hints at the existence and structure of planetary systems. Millimeter-wavelength observations probe large grains that trace the location of planetesimal belts. The FEPS (Formation and Evolution of Planetary Systems) $Spitzer$ Legacy survey of nearby young solar analogues yielded a sample of five debris disk-hosting stars with millimeter flux suitable for interferometric follow-up. We present observations with the Submillimeter Array (SMA) and the Combined Array for Research in Millimeter-wave Astronomy (CARMA) at ~2" resolution that spatially resolve the debris disks around these nearby ($d\sim$50 pc) stars. Two of the five disks (HD 377, HD 8907) are spatially resolved for the first time and one (HD 104860) is resolved at millimeter wavelengths for the first time. We combine our new observations with archival SMA and Atacama Large Millimeter/Submillimeter Array (ALMA) data to enable a uniform analysis of the full five-object sample. We simultaneously model the broad-band photometric data and resolved millimeter visibilities to constrain the dust temperatures and disk morphologies, and perform an MCMC analysis to fit for basic structural parameters. We find that the radii and widths of the cold outer belts exhibit properties consistent with scaled-up versions of the Solar System's Kuiper Belt. All the disks exhibit characteristic grain sizes comparable to the blowout size, and all the resolved observations of emission from large dust grains are consistent with an axisymmetric dust distribution to within the uncertainties. These results are consistent with comparable studies carried out at infrared wavelengths.

Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane and water vapor absorption. Modeling of the spectra and photometry yields a luminosity of L/LS=1.6-4.0 x 10-6 and an effective temperature of 600-750 K. For this age and luminosity, "hot-start" formation models indicate a mass twice that of Jupiter. This planet also has a sufficiently low luminosity to be consistent with the "cold- start" core accretion process that may have formed Jupiter.

Turbulence can transport angular momentum in protoplanetary disks and influence the growth and evolution of planets. With spatially and spectrally resolved molecular emission line measurements provided by (sub)millimeter interferometric observations, it is possible to directly measure non-thermal motions in the disk gas that can be attributed to this turbulence. We report a new constraint on the turbulence in the disk around HD 163296, a nearby young A star, determined from ALMA Science Verification observations of four CO emission lines (the CO(3-2), CO(2-1), 13CO(2-1), and C18O(2-1) transitions). The different optical depths for these lines permit probes of non-thermal line-widths at a range of physical conditions (temperature and density) and depths into the disk interior. We derive stringent limits on the non-thermal motions in the upper layers of the outer disk such that any contribution to the line-widths from turbulence is <3% of the local sound speed. These limits are approximately an order of magnitude lower than theoretical predictions for full-blown MHD turbulence driven by the magneto-rotational instability, potentially suggesting that this mechanism is less efficient in the outer (R>30AU) disk than has been previously considered.

Combining 3D smoothed-particle hydrodynamics and Monte Carlo radiative transfer calculations, we examine the morphology of spiral density waves induced by gravitational instability (GI) in protoplanetary disks, as they would appear in direct images at near-infrared (NIR) wavelengths. We find that systems with disk-to-star-mass ratios q=M_disk/M_star that are ~0.25 or more may produce prominent spiral arms in NIR imaging, remarkably resembling features observed in the MWC 758 and SAO 206462 systems. The contrast of GI-induced arms at NIR wavelengths can reach a factor of ~3, and their pitch angles are about 10-15 degree. The dominant azimuthal wavenumber of GI-induced spiral arms roughly obeys m~1/q in the range 2<~1/q<~8. In particular, a massive disk with q~0.5 can exhibit grand-design m=2 spirals. GI-induced arms are in approximate corotation with the local disk, and may therefore trap dust particles by pressure drag. Although GI can produce NIR spiral arms with morphologies, contrasts, and pitch angles similar to those reported in recent observations, it also makes other demands that may or may not be satisfied in any given system. A GI origin requires that the spirals be relatively compact, on scales <~100 AU; that the disk be massive, q>~0.25; and that the accretion rate Mdot be high, on the order of 1e-6 solar mass per year.

Numerical hydrodynamics calculations are performed to determine conditions under which giant planet eccentricities can be excited by parent gas disks. Unlike in other studies, Jupiter-mass planets are found to have their eccentricities amplified --- provided their orbits start eccentric. We disentangle the web of co-rotation, co-orbital, and external resonances to show that this finite-amplitude instability is consistent with that predicted analytically. Ellipticities can grow until they reach of order the disk's aspect ratio, beyond which the external Lindblad resonances that excite eccentricity are weakened by the planet's increasingly supersonic epicyclic motion. Forcing the planet to still larger eccentricities causes catastrophic eccentricity damping as the planet collides into gap walls. For standard parameters, the range of eccentricities for instability is modest; the threshold eccentricity for growth ($\sim$$0.04$) is not much smaller than the final eccentricity to which orbits grow ($\sim$$0.07$). If this threshold eccentricity can be lowered (perhaps by non-barotropic effects), and if the eccentricity driving documented here survives in 3D, it may robustly explain the low-to-moderate eccentricities $\lesssim 0.1$ exhibited by many giant planets (including Jupiter and Saturn), especially those without planetary or stellar companions.

Planets acquire atmospheres from their parent circumstellar disks. We derive a general analytic expression for how the atmospheric mass grows with time $t$, as a function of the underlying core mass $M_{\rm core}$ and nebular conditions, including the gas metallicity $Z$. Planets accrete as much gas as can cool: an atmosphere's doubling time is given by its Kelvin-Helmholtz time. Dusty atmospheres behave differently from atmospheres made dust-free by grain growth and sedimentation. The gas-to-core mass ratio (GCR) of a dusty atmosphere scales as GCR $\propto t^{0.4} M_{\rm core}^{1.7} Z^{-0.4} \mu_{\rm rcb}^{3.4}$, where $\mu_{\rm rcb} \propto 1/(1-Z)$ (for $Z$ not too close to 1) is the mean molecular weight at the innermost radiative-convective boundary. This scaling applies across all orbital distances and nebular conditions for dusty atmospheres; their radiative-convective boundaries, which regulate cooling, are not set by the external environment, but rather by the internal microphysics of dust sublimation, H$_2$ dissociation, and the formation of H$^-$. By contrast, dust-free atmospheres have their radiative boundaries at temperatures $T_{\rm rcb}$ close to nebular temperatures $T_{\rm out}$, and grow faster at larger orbital distances where cooler temperatures, and by extension lower opacities, prevail. At 0.1 AU in a gas-poor nebula, GCR $\propto t^{0.4} T_{\rm rcb}^{-1.9} M_{\rm core}^{1.6} Z^{-0.4} \mu_{\rm rcb}^{3.3}$, while beyond 1 AU in a gas-rich nebula, GCR $\propto t^{0.4} T_{\rm rcb}^{-1.5} M_{\rm core}^1 Z^{-0.4}\mu_{\rm rcb}^{2.2}$. We confirm our analytic scalings against detailed numerical models for objects ranging in mass from Mars (0.1 $M_\oplus$) to the most extreme super-Earths (10-20 $M_\oplus$), and explain why heating from planetesimal accretion cannot prevent the latter from undergoing runaway gas accretion.

We present $H$-band observations of $\beta$ Pic with the Gemini Planet Imager's (GPI's) polarimetry mode that reveal the debris disk between ~0.3" (~6 AU) and ~1.7" (~33 AU), while simultaneously detecting $\beta$ Pic $b$. The polarized disk image was fit with a dust density model combined with a Henyey-Greenstein scattering phase function. The best fit model indicates a disk inclined to the line of sight ($\phi=85.27{\deg}^{+0.26}_{-0.19}$) with a position angle $\theta_{PA}=30.35{\deg}^{+0.29}_{-0.28}$ (slightly offset from the main outer disk, $\theta_{PA}\approx29{\deg}$), that extends from an inner disk radius of $23.6^{+0.9}_{-0.6}$ AU to well outside GPI's field of view. In addition, we present an updated orbit for $\beta$ Pic $b$ based on new astrometric measurements taken in GPI's spectroscopic mode spanning 14 months. The planet has a semi-major axis of $a=9.2^{+1.5}_{-0.4}$AU, with an eccentricity $e\leq 0.26$. The position angle of the ascending node is $\Omega=31.75{\deg}\pm0.15$, offset from both the outer main disk and the inner disk seen in the GPI image. The orbital fit constrains the stellar mass of $\beta$ Pic to $1.60\pm0.05 M_{\odot}$. Dynamical sculpting by $\beta$ Pic $b$ cannot easily account for the following three aspects of the inferred disk properties: 1) the modeled inner radius of the disk is farther out than expected if caused by $\beta$ Pic b; 2) the mutual inclination of the inner disk and $\beta$ Pic $b$ is $4{\deg}$, when it is expected to be closer to zero; and 3) the aspect ratio of the disk ($h_0 = 0.137^{+0.005}_{-0.006}$) is larger than expected from interactions with $\beta$ Pic $b$ or self-stirring by the disk's parent bodies.

We present Gemini Planet Imager (GPI) observations of AU Microscopii, a young M dwarf with an edge-on, dusty debris disk. Integral field spectroscopy and broadband imaging polarimetry were obtained during the commissioning of GPI. In our broadband imaging polarimetry observations, we detect the disk only in total intensity and find asymmetries in the morphology of the disk between the southeast and northwest sides. The southeast side of the disk exhibits a bump at 1$''$ (10 AU projected separation) that is three times more vertically extended and three times fainter in peak surface brightness than the northwest side at similar separations. This part of the disk is also vertically offset by 69$\pm$30 mas to the northeast at 1$''$ when compared to the established disk mid-plane and consistent with prior ALMA and Hubble Space Telescope/STIS observations. We see hints that the southeast bump might be a result of detecting a horizontal sliver feature above the main disk that could be the disk backside. Alternatively when including the morphology of the northwest side, where the disk mid-plane is offset in the opposite direction $\sim$50 mas between 0$.''$4 and 1$.''$2, the asymmetries suggest a warp-like feature. Using our integral field spectroscopy data to search for planets, we are 50% complete for $\sim$4 $M_\mathrm{Jup}$ planets at 4 AU. We detect a source, resolved only along the disk plane, that could either be a candidate planetary mass companion or a compact clump in the disk.

Planets with sizes between those of Earth and Neptune divide into two populations: purely rocky bodies whose atmospheres contribute negligibly to their sizes, and larger gas-enveloped planets possessing voluminous and optically thick atmospheres. We show that whether a planet forms rocky or gas-enveloped depends on the solid surface density of its parent disk. Assembly times for rocky cores are sensitive to disk solid surface density. Lower surface densities spawn smaller planetary embryos; to assemble a core of given mass, smaller embryos require more mergers between bodies farther apart and therefore exponentially longer formation times. Gas accretion simulations yield a rule of thumb that a rocky core must be at least 2$M_\oplus$ before it can acquire a volumetrically significant atmosphere from its parent nebula. In disks of low solid surface density, cores of such mass appear only after the gas disk has dissipated, and so remain purely rocky. Higher surface density disks breed massive cores more quickly, within the gas disk lifetime, and so produce gas-enveloped planets. We test model predictions against observations, using planet radius as an observational proxy for gas-to-rock content and host star metallicity as a proxy for disk solid surface density. Theory can explain the observation that metal-rich stars host predominantly gas-enveloped planets.

Close-in super-Earths having radii 1--4 $R_\oplus$ may possess hydrogen atmospheres comprising a few percent by mass of their rocky cores. We determine the conditions under which such atmospheres can be accreted by cores from their parent circumstellar disks. Accretion from the nebula is problematic because it is too efficient: we find that 10-$M_\oplus$ cores embedded in solar metallicity disks tend to undergo runaway gas accretion and explode into Jupiters, irrespective of orbital location. The threat of runaway is especially dire at $\sim$0.1 AU, where solids may coagulate on timescales orders of magnitude shorter than gas clearing times; thus nascent atmospheres on close-in orbits are unlikely to be supported against collapse by planetesimal accretion. The time to runaway accretion is well approximated by the cooling time of the atmosphere's innermost convective zone, whose extent is controlled by where H$_2$ dissociates. Insofar as the temperatures characterizing H$_2$ dissociation are universal, timescales for core instability tend not to vary with orbital distance --- and to be alarmingly short for 10-$M_\oplus$ cores. Nevertheless, in the thicket of parameter space, we identify two scenarios, not mutually exclusive, that can reproduce the preponderance of percent-by-mass atmospheres for super-Earths at $\sim$0.1 AU, while still ensuring the formation of Jupiters at $\gtrsim 1$ AU. Scenario (a): planets form in disks with dust-to-gas ratios that range from $\sim$20$\times$ solar at 0.1 AU to $\sim$2$\times$ solar at 5 AU. Scenario (b): the final assembly of super-Earth cores from mergers of proto-cores --- a process that completes quickly at $\sim$0.1 AU once begun --- is delayed by gas dynamical friction until just before disk gas dissipates completely.

The brightest transitional protoplanetary disks are often azimuthally asymmetric: their mm-wave thermal emission peaks strongly on one side. Dust overdensities can exceed $\sim$100:1, while gas densities vary by factors less than a few. We propose that these remarkable ALMA observations---which may bear on how planetesimals form---reflect a gravitational global mode in the gas disk. The mode is (1) fast---its pattern speed equals the disk's mean Keplerian frequency; (2) of azimuthal wavenumber $m=1$, displacing the host star from the barycenter; and (3) Toomre-stable. We solve for gas streamlines including the indirect stellar potential in the frame rotating with the pattern speed, under the drastic simplification that gas does not feel its own gravity. Near co-rotation, the gas disk takes the form of a horseshoe-shaped annulus. Dust particles with aerodynamic stopping times much shorter or much longer than the orbital period are dragged by gas toward the horseshoe center. For intermediate stopping times, dust converges toward a $\sim$45$^\circ$-wide arc on the co-rotation circle. Particles that do not reach their final accumulation points within disk lifetimes, either because of gas turbulence or long particle drift times, conform to horseshoe-shaped gas streamlines. Our mode is not self-consistent because we neglect gas self-gravity; still, we expect that trends between accumulation location and particle size, similar to those we have found, are generically predicted by fast modes and are potentially observable. Unlike vortices, global modes are not restricted in radial width to the pressure scale height; their large radial and azimuthal extents may better match observations.

The orbits of giant extrasolar planets often have surprisingly small semi-major axes, large eccentricities, or severe misalignments between their normals and their host stars' spin axes. In some formation scenarios invoking Kozai-Lidov oscillations, an external planetary companion drives a planet onto an orbit having these properties. The mutual inclinations for Kozai-Lidov oscillations can be large and have not been confirmed observationally. Here we deduce that observed eccentric warm Jupiters with eccentric giant companions have mutual inclinations that oscillate between 35-65 deg. Our inference is based on the pairs' observed apsidal separations, which cluster near 90 deg. The near-orthogonality of periapse directions is effected by the outer companion's quadrupolar and octupolar potentials. These systems may be undergoing a stalled version of tidal migration that produces warm Jupiters over hot Jupiters, and provide evidence for a population of multi-planet systems that are not flat and have been sculpted by Kozai-Lidov oscillations.

Characterizing turbulence in protoplanetary disks is crucial for understanding how they accrete and spawn planets. Recent measurements of spectral line broadening promise to diagnose turbulence, with different lines probing different depths. We use 3D local hydrodynamic simulations of cooling, self-gravitating disks to resolve how motions driven by "gravito-turbulence" vary with height. We find that gravito-turbulence is practically as vigorous at altitude as at depth: even though gas at altitude is much too rarefied to be itself self-gravitating, it is strongly forced by self-gravitating overdensities at the midplane. The long-range nature of gravity means that turbulent velocities are nearly uniform vertically, increasing by just a factor of 2 from midplane to surface, even as the density ranges over nearly three orders of magnitude. The insensitivity of gravito-turbulence to height contrasts with the behavior of disks afflicted by the magnetorotational instability (MRI); in the latter case, non-circular velocities increase by at least a factor of 15 from midplane to surface, with various non-ideal effects only magnifying this factor. The distinct vertical profiles of gravito-turbulence vs. MRI turbulence may be used in conjunction with measurements of non-thermal linewidths at various depths to identify the source of transport in protoplanetary disks.

We present multiwavelength photometry, high angular resolution imaging, and radial velocities, of the unique and confounding disintegrating low-mass planet candidate KIC 12557548b. Our high angular resolution imaging, which includes spacebased HST/WFC3 observations in the optical, and groundbased Keck/NIRC2 observations in K'-band, allow us to rule-out background and foreground candidates at angular separations greater than 0.2 arcsec that are bright enough to be responsible for the transits we associate with KIC 12557548. Our radial velocity limit from Keck/HIRES allows us to rule-out bound, low-mass stellar companions to KIC 12557548 on orbits less than 10 years, as well as placing an upper-limit on the mass of the candidate planet of 1.2 Jupiter masses; therefore, the combination of our radial velocities, high angular-resolution imaging, and photometry are able to rule-out most false positive interpretations of the transits. Our precise multiwavelength photometry includes two simultaneous detections of the transit of KIC 12557548b using CFHT/WIRCam at 2.15 microns and the Kepler space telescope at 0.6 microns, as well as simultaneous null-detections of the transit by Kepler and HST/WFC3 at 1.4 microns. Our simultaneous HST/WFC3 and Kepler null-detections, provide no evidence for radically different transit depths at these wavelengths. Our simultaneous CFHT/WIRCam detections in the near-infrared and with Kepler in the optical reveal very similar transit depths (the average ratio of the transit depths at ~2.15 microns compared to ~0.6 microns is: 1.02 +/- 0.20). This suggests that if the transits we observe are due to scattering from single-size particles streaming from the planet in a comet-like tail, then the particles must be ~0.5 microns in radius or larger, which would favor that KIC 12557548b is a sub-Mercury, rather than super-Mercury, mass planet.

Gap clearing by giant planets has been proposed to explain the optically thin cavities observed in many protoplanetary disks. How much material remains in the gap determines not only how detectable young planets are in their birth environments, but also how strong corotation torques are, which impacts how planets can survive fast orbital migration. We determine numerically how the average surface density inside the gap, sigma_gap, depends on planet-to-star mass ratio q, Shakura-Sunyaev viscosity parameter alpha, and disk height-to-radius aspect ratio h/r. Our results are derived from our new GPU-accelerated Lagrangian hydrodynamical code PEnGUIn, and are verified by independent simulations with ZEUS90. For Jupiter-like planets, we find sigma_gap \propto q^-2.2 alpha^1.4 (h/r)^6.6, and for near brown dwarf masses, sigma_gap \propto q^-1 alpha^1.3 (h/r)^6.1. Surface density contrasts inside and outside gaps can be as large as 10^4, even when the planet does not accrete. We derive a simple analytic scaling, sigma_gap \propto q^-2 alpha^1 (h/r)^5, that compares reasonably well to empirical results, especially at low Neptune-like masses, and use discrepancies to highlight areas for progress.

Protoplanetary "transition" disks have large, mass-depleted central cavities, yet also deliver gas onto their host stars at rates comparable to disks without holes. The paradox of simultaneous transparency and accretion can be explained if gas flows inward at much higher radial speeds inside the cavity than outside the cavity, since surface density (and by extension optical depth) varies inversely with inflow velocity at fixed accretion rate. Radial speeds within the cavity might even have to approach free-fall values to explain the huge surface density contrasts inferred for transition disks. We identify observational diagnostics of fast radial inflow in channel maps made in optically thick spectral lines. Signatures include (1) twisted isophotes in maps made at low systemic velocities and (2) rotation of structures observed between maps made in high-velocity line wings. As a test case, we apply our new diagnostic tools to archival ALMA data on the transition disk HD 142527, and uncover evidence for free-fall radial velocities inside its cavity. Although the observed kinematics are also consistent with a disk warp, the radial inflow scenario is preferred because it predicts low surface densities that appear consistent with recent observations of optically thin CO isotopologues in this disk. How material in the disk cavity sheds its angular momentum wholesale to fall freely onto the star is an unsolved problem; gravitational torques exerted by giant planets or brown dwarfs are briefly discussed as a candidate mechanism.

Short-period exoplanets can have dayside surface temperatures surpassing 2000 K, hot enough to vaporize rock and drive a thermal wind. Small enough planets evaporate completely. We construct a radiative-hydrodynamic model of atmospheric escape from strongly irradiated, low-mass rocky planets, accounting for dust-gas energy exchange in the wind. Rocky planets with masses < 0.1 M_Earth (less than twice the mass of Mercury) and surface temperatures > 2000 K are found to disintegrate entirely in < 10 Gyr. When our model is applied to Kepler planet candidate KIC 12557548b --- which is believed to be a rocky body evaporating at a rate of dM/dt > 0.1 M_Earth/Gyr --- our model yields a present-day planet mass of < 0.02 M_Earth or less than about twice the mass of the Moon. Mass loss rates depend so strongly on planet mass that bodies can reside on close-in orbits for Gyrs with initial masses comparable to or less than that of Mercury, before entering a final short-lived phase of catastrophic mass loss (which KIC 12557548b has entered). Because this catastrophic stage lasts only up to a few percent of the planet's life, we estimate that for every object like KIC 12557548b, there should be 10--100 close-in quiescent progenitors with sub-day periods whose hard-surface transits may be detectable by Kepler --- if the progenitors are as large as their maximal, Mercury-like sizes (alternatively, the progenitors could be smaller and more numerous). According to our calculations, KIC 12557548b may have lost ~70% of its formation mass; today we may be observing its naked iron core.

Close-in super-Earths, with radii R = 2-5 R_Earth and orbital periods P < 100 days, orbit more than half, and perhaps nearly all Sun-like stars in the universe. We use this omnipresent population to construct the minimum-mass extrasolar nebula (MMEN), the circumstellar disk of solar-composition solids and gas from which such planets formed, if they formed near their current locations and did not migrate. In a series of back-of-the-envelope calculations, we demonstrate how in-situ formation in the MMEN is fast, efficient, and can reproduce many of the observed properties of close-in super-Earths, including their gas-to-rock fractions. Testable predictions are discussed.

Dust particles sediment toward the midplanes of protoplanetary disks, forming dust-rich sublayers encased in gas. What densities must the particle sublayer attain before it can fragment by self-gravity? We describe various candidate threshold densities. One of these is the Roche density, which is that required for a strengthless satellite to resist tidal disruption by its primary. Another is the Toomre density, which is that required for de-stabilizing self-gravity to defeat the stabilizing influences of pressure and rotation. We show that for sublayers containing aerodynamically well-coupled dust, the Toomre density exceeds the Roche density by many (up to about 4) orders of magnitude. We present 3D shearing box simulations of self-gravitating, stratified, dust-gas mixtures to test which of the candidate thresholds is relevant for collapse. All our simulations indicate that the larger Toomre density is required for collapse. This result is sensible because sublayers are readily stabilized by pressure. Sound-crossing times for thin layers are easily shorter than free-fall times, and the effective sound speed in dust-gas suspensions decreases only weakly with the dust-to-gas ratio (as the inverse square root). Our findings assume that particles are small enough that their stopping times in gas are shorter than all other timescales. Relaxing this assumption may lower the threshold for gravitational collapse back down to the Roche criterion. In particular, if the particle stopping time becomes longer than the sound-crossing time, sublayers may lose pressure support and become gravitationally unstable.

We present 1.3 millimeter ALMA Cycle 0 observations of the edge-on debris disk around the nearby, ~10 Myr-old, M-type star AU Mic. These observations obtain 0.6 arcsec (6 AU) resolution and reveal two distinct emission components: (1) the previously known dust belt that extends to a radius of 40 AU, and (2) a newly recognized central peak that remains unresolved. The cold dust belt of mass about 1 lunar mass is resolved in the radial direction with a rising emission profile that peaks sharply at the location of the outer edge of the "birth ring" of planetesimals hypothesized to explain the midplane scattered light gradients. No significant asymmetries are discerned in the structure or position of this dust belt. The central peak identified in the ALMA image is ~6 times brighter than the stellar photosphere, which indicates an additional emission process in the inner regions of the system. Emission from a stellar corona or activity may contribute, but the observations show no signs of temporal variations characteristic of radio-wave flares. We suggest that this central component may be dominated by dust emission from an inner planetesimal belt of mass about 0.01 lunar mass, consistent with a lack of emission shortward of 25 microns and a location <3 AU from the star. Future millimeter observations can test this assertion, as an inner dust belt should be readily separated from the central star at higher angular resolution.

When transiting their host stars, hot Jupiters absorb about 10% of the light in the wings of the stellar Lyman-alpha emission line. The absorption occurs at wavelengths Doppler-shifted from line center by +/- 100 km/s - larger than the thermal speeds with which partially neutral, 10E4 K hydrogen escapes from hot Jupiter atmospheres. It has been proposed that the absorption arises from 10E6 K hydrogen from the host stellar wind, made momentarily neutral by charge exchange with planetary H I. The +/-100 km/s velocities would then be attributed to the typical velocity dispersions of protons in the stellar wind - as inferred from spacecraft measurements of the Solar wind. To test this proposal, we perform 2D hydrodynamic simulations of colliding hot Jupiter and stellar winds, augmented by a chemistry module to compute the amount of hot neutral hydrogen produced by charge exchange. The Kelvin-Helmholtz instability mixes the two winds; in the mixing layer, charge exchange reactions establish, within tens of seconds. Enough hot neutral hydrogen is generated to reproduce the transit observations, and the amount of absorption converges with both spatial resolution and time. Our calculations support the idea that charge transfer between colliding winds correctly explains Lyman-alpha transit observations - modulo the effects of magnetic fields, which we do not model but which may suppress mixing. Other neglected effects include, in order of decreasing importance, rotational forces related to orbital motion, gravity, and stellar radiation pressure; we discuss quantitatively the errors introduced by our approximations. How hot stellar hydrogen cools when it collides with cold planetary hydrogen is also considered; a more careful treatment of how the mixing layer thermally equilibrates might explain the recent detection of Balmer H-alpha absorption in transiting hot Jupiters.

Kilometer-sized moonlets in Saturn's A ring create S-shaped wakes called "propellers" in surrounding material. The Cassini spacecraft has tracked the motions of propellers for several years and finds that they deviate from Keplerian orbits having constant semimajor axes. The inferred orbital migration is known to switch sign. We show using a statistical test that the time series of orbital longitudes of the propeller Bl\'eriot is consistent with that of a time-integrated Gaussian random walk. That is, Bl\'eriot's observed migration pattern is consistent with being stochastic. We further show, using a combination of analytic estimates and collisional N-body simulations, that stochastic migration of the right magnitude to explain the Cassini observations can be driven by encounters with ring particles 10-20 m in radius. That the local ring mass is concentrated in decameter-sized particles is supported on independent grounds by occultation analyses.

Few scientific discoveries have captured the public imagination like the explosion of exoplanetary science during the past two decades. This work has fundamentally changed our picture of Earth's place in the Universe and led NASA to make significant investments towards understanding the demographics of exoplanetary systems and the conditions that lead to their formation. The story of the formation and evolution of exoplanetary systems is essentially the story of the circumstellar gas and dust that are initially present in the protostellar environment; in order to understand the variety of planetary systems observed, we need to understand the life cycle of circumstellar gas from its initial conditions in protoplanetary disks to its endpoint as planets and their atmospheres. In this white paper response to NASA's Request for Information "Science Objectives and Requirements for the Next NASA UV/Visible Astrophysics Mission Concepts (NNH12ZDA008L)", we describe scientific programs that would use the unique capabilities of a future NASA ultraviolet (UV)/visible space observatory to make order-of-magnitude advances in our understanding of the life cycle of circumstellar gas.

"Propellers" are features in Saturn's A ring associated with moonlets that open partial gaps. They exhibit non-Keplerian motion (Tiscareno 2010); the longitude residuals of the best-observed propeller, "Bl\'eriot," appear consistent with a sinusoid of period ~4 years. Pan and Chiang (2010) proposed that propeller moonlets librate in "frog resonances" with co-orbiting ring material. By analogy with the restricted three-body problem, they treated the co-orbital material as stationary in the rotating frame and neglected non-co-orbital material. Here we use simple numerical experiments to extend the frog model, including feedback due to the gap's motion, and drag associated with the Lindblad disk torques that cause Type I migration. Because the moonlet creates the gap, we expect the gap centroid to track the moonlet, but only after a time delay t_diff, the time for a ring particle to travel from conjunction with the moonlet to the end of the gap. We find that frog librations can persist only if t_diff exceeds the frog libration period P_lib, and if damping from Lindblad torques balances driving from co-orbital torques. If t_diff << P_lib, then the libration amplitude damps to zero. In the case of Bl\'eriot, the frog resonance model can reproduce the observed libration period P_lib ~ 4 yr. However, our simple feedback prescription suggests that Bl\'eriot's t_diff ~ 0.01P_lib, which is inconsistent with the observed libration amplitude of 260 km. We urge more accurate treatments of feedback to test the assumptions of our toy models.

We report here on the discovery of stellar occultations, observed with Kepler, that recur periodically at 15.685 hour intervals, but which vary in depth from a maximum of 1.3% to a minimum that can be less than 0.2%. The star that is apparently being occulted is KIC 12557548, a K dwarf with T_eff = 4400 K and V = 16. Because the eclipse depths are highly variable, they cannot be due solely to transits of a single planet with a fixed size. We discuss but dismiss a scenario involving a binary giant planet whose mutual orbit plane precesses, bringing one of the planets into and out of a grazing transit. We also briefly consider an eclipsing binary, that either orbits KIC 12557548 in a hierarchical triple configuration or is nearby on the sky, but we find such a scenario inadequate to reproduce the observations. We come down in favor of an explanation that involves macroscopic particles escaping the atmosphere of a slowly disintegrating planet not much larger than Mercury. The particles could take the form of micron-sized pyroxene or aluminum oxide dust grains. The planetary surface is hot enough to sublimate and create a high-Z atmosphere; this atmosphere may be loaded with dust via cloud condensation or explosive volcanism. Atmospheric gas escapes the planet via a Parker-type thermal wind, dragging dust grains with it. We infer a mass loss rate from the observations of order 1 M_E/Gyr, with a dust-to-gas ratio possibly of order unity. For our fiducial 0.1 M_E planet, the evaporation timescale may be ~0.2 Gyr. Smaller mass planets are disfavored because they evaporate still more quickly, as are larger mass planets because they have surface gravities too strong to sustain outflows with the requisite mass-loss rates. The occultation profile evinces an ingress-egress asymmetry that could reflect a comet-like dust tail trailing the planet; we present simulations of such a tail.

Clumpy structure in the debris disk around Vega has been previously reported at millimeter wavelengths and attributed to concentrations of dust grains trapped in resonances with an unseen planet. However, recent imaging at similar wavelengths with higher sensitivity has disputed the observed structure. We present three new millimeter-wavelength observations that help to resolve the puzzling and contradictory observations. We have observed the Vega system with the Submillimeter Array (SMA) at a wavelength of 880 um and angular resolution of 5"; with the Combined Array for Research in Millimeter-wave Astronomy (CARMA) at a wavelength of 1.3 mm and angular resolution of 5"; and with the Green Bank Telescope (GBT) at a wavelength of 3.3 mm and angular resolution of 10". Despite high sensitivity and short baselines, we do not detect the Vega debris disk in either of the interferometric data sets (SMA and CARMA), which should be sensitive at high significance to clumpy structure based on previously reported observations. We obtain a marginal (3-sigma) detection of disk emission in the GBT data; the spatial distribution of the emission is not well constrained. We analyze the observations in the context of several different models, demonstrating that the observations are consistent with a smooth, broad, axisymmetric disk with inner radius 20-100 AU and width >50 AU. The interferometric data require that at least half of the 860 um emission detected by previous single-dish observations with the James Clerk Maxwell Telescope be distributed axisymmetrically, ruling out strong contributions from flux concentrations on spatial scales of <100 AU. These observations support recent results from the Plateau de Bure Interferometer indicating that previous detections of clumpy structure in the Vega debris disk were spurious.

Whether protoplanetary disks accrete at observationally significant rates by the magnetorotational instability (MRI) depends on how well ionized they are. Disk surface layers ionized by stellar X-rays are susceptible to charge neutralization by small condensates, ranging from ~0.01-micron-sized grains to angstrom-sized polycyclic aromatic hydrocarbons (PAHs). Ion densities in X-ray-irradiated surfaces are so low that ambipolar diffusion weakens the MRI. Here we show that ionization by stellar far-ultraviolet (FUV) radiation enables full-blown MRI turbulence in disk surface layers. Far-UV ionization of atomic carbon and sulfur produces a plasma so dense that it is immune to ion recombination on grains and PAHs. The FUV-ionized layer, of thickness 0.01--0.1 g/cm^2, behaves in the ideal magnetohydrodynamic limit and can accrete at observationally significant rates at radii > 1--10 AU. Surface layer accretion driven by FUV ionization can reproduce the trend of increasing accretion rate with increasing hole size seen in transitional disks. At radii < 1--10 AU, FUV-ionized surface layers cannot sustain the accretion rates generated at larger distance, and unless turbulent mixing of plasma can thicken the MRI-active layer, an additional means of transport is needed. In the case of transitional disks, it could be provided by planets.

Transitional T Tauri disks have optically thin holes with radii > 10 AU, yet accrete up to the median T Tauri rate. Multiple planets inside the hole can torque the gas to high radial speeds over large distances, reducing the local surface density while maintaining accretion. Thus multi-planet systems, together with reductions in disk opacity due to grain growth, can explain how holes can be simultaneously transparent and accreting. There remains the problem of how outer disk gas diffuses into the hole. Here it has been proposed that the magnetorotational instability (MRI) erodes disk surface layers ionized by stellar X-rays. In contrast to previous work, we find that the extent to which surface layers are MRI-active is limited not by ohmic dissipation but by ambipolar diffusion, the latter measured by Am: the number of times a neutral hydrogen molecule collides with ions in a dynamical time. Simulations by Hawley & Stone showed that Am ~ 100 is necessary for ions to drive MRI turbulence in neutral gas. We calculate that in X-ray-irradiated surface layers, Am typically varies from ~0.001 to 1, depending on the abundance of charge-adsorbing polycyclic aromatic hydrocarbons, whose properties we infer from Spitzer observations. We conclude that ionization of H_2 by X-rays and cosmic rays can sustain, at most, only weak MRI turbulence in surface layers 1 to 10 g/cm^2 thick, and that accretion rates in such layers are too small compared to observed accretion rates for the majority of disks.

Dust at the midplane of a circumstellar disk can become gravitationally unstable and fragment into planetesimals if the local dust-to-gas density ratio mu is sufficiently high. We simulate how dust settles in passive disks and ask how high mu can become. We settle the dust using a 1D code and test for dynamical stability using a 3D shearing box code. This scheme allows us to explore the behavior of small particles having short but non-zero stopping times in gas: 0 < t_stop << the orbital period. The streaming instability is thereby filtered out. Dust settles until shearing instabilities in the edges of the dust layer threaten to overturn the entire layer. In this state of marginal stability, mu=2.9 for a disk whose bulk (height-integrated) metallicity is solar. For a disk whose bulk metallicity is 4x solar, mu reaches 26.4. These maximum values of mu, which depend on the background radial pressure gradient, are so large that gravitational instability of small particles is viable in disks whose bulk metallicities are just a few (<4) times solar. Earlier studies assumed that dust settles until the Richardson number Ri is spatially constant. Our simulations are free of this assumption but provide support for it within the dust layer's edges, with the proviso that Ri increases with bulk metallicity in the same way that we found in Paper I. Only modest enhancements in bulk metallicity are needed to spawn planetesimals directly from small particles.

Gravitational instability (GI) of a dust-rich layer at the midplane of a gaseous circumstellar disk is one proposed mechanism to form planetesimals, the building blocks of rocky planets and gas giant cores. Self-gravity competes against the Kelvin-Helmholtz instability (KHI): gradients in dust content drive a vertical shear which risks overturning the dusty subdisk and forestalling GI. To understand the conditions under which the disk can resist the KHI, we perform 3D simulations of stratified subdisks in the limit that dust particles are small and aerodynamically well coupled to gas. This limit screens out the streaming instability and isolates the KHI. Each subdisk is assumed to have a vertical density profile given by a spatially constant Richardson number Ri. We vary Ri and the midplane dust-to-gas ratio mu and find that the critical Richardson number dividing KH-unstable from KH-stable flows is not unique; rather Ri_crit grows nearly linearly with mu for mu=0.3-10. Only for disks of bulk solar metallicity is Ri_crit ~ 0.2, close to the classical value. Our results suggest that a dusty sublayer can gravitationally fragment and presumably spawn planetesimals if embedded within a solar metallicity gas disk ~4x more massive than the minimum-mass solar nebula; or a minimum-mass disk having ~3x solar metallicity; or some intermediate combination of these two possibilities. Gravitational instability seems possible without resorting to the streaming instability or to turbulent concentration of particles.

"Propellers" in planetary rings are disturbances in ring material excited by moonlets that open only partial gaps. We describe a new type of co-orbital resonance that can explain the observed non-Keplerian motions of propellers. The resonance is between the moonlet underlying the propeller, and co-orbiting ring particles downstream of the moonlet where the gap closes. The moonlet librates within the gap about an equilibrium point established by co-orbiting material and stabilized by the Coriolis force. In the limit of small libration amplitude, the libration period scales linearly with the gap azimuthal width and inversely as the square root of the co-orbital mass. The new resonance recalls but is distinct from conventional horseshoe and tadpole orbits; we call it the "frog" resonance, after the relevant term in equine hoof anatomy. For a ring surface density and gap geometry appropriate for the propeller Bl\'eriot in Saturn's A ring, our theory predicts a libration period of ~4 years, similar to the ~3.7 year period over which Bl\'eriot's orbital longitude is observed to vary. These librations should be subtracted from the longitude data before any inferences about moonlet migration are made.

We present Hubble Space Telescope optical coronagraphic polarization imaging observations of the dusty debris disk HD 61005. The scattered light intensity image and polarization structure reveal a highly inclined disk with a clear asymmetric, swept back component, suggestive of significant interaction with the ambient interstellar medium. The combination of our new data with the published 1.1 micron discovery image shows that the grains are blue scattering with no strong color gradient as a function of radius, implying predominantly sub-micron sized grains. We investigate possible explanations that could account for the observed swept back, asymmetric morphology. Previous work has suggested that HD 61005 may be interacting with a cold, unusually dense interstellar cloud. However, limits on the intervening interstellar gas column density from an optical spectrum of HD 61005 in the Na I D lines render this possibility unlikely. Instead, HD 61005 may be embedded in a more typical warm, low-density cloud that introduces secular perturbations to dust grain orbits. This mechanism can significantly distort the ensemble disk structure within a typical cloud crossing time. For a counterintuitive relative flow direction--parallel to the disk midplane--we find that the structures generated by these distortions can very roughly approximate the HD 61005 morphology. Future observational studies constraining the direction of the relative interstellar medium flow will thus provide an important constraint for future modeling. Independent of the interpretation for HD 61005, we expect that interstellar gas drag likely plays a role in producing asymmetries observed in other debris disk systems, such as HD 15115 and Delta-Velorum.

Planets are built from planetesimals: solids larger than a kilometer which grow by colliding pairwise. Planetesimals themselves are unlikely to form by two-body collisions; sub-km objects have gravitational fields individually too weak, and electrostatic attraction is too feeble for growth beyond a few cm. We review the possibility that planetesimals form when self-gravity brings together vast ensembles of small particles. Even when self-gravity is weak, aerodynamic processes can accumulate solids relative to gas, paving the way for gravitational collapse. Particles pile up as they drift radially inward. Gas turbulence stirs particles, but can also seed collapse by clumping them. While the feedback of solids on gas triggers vertical shear instabilities that obstruct self-gravity, this same feedback triggers streaming instabilities that strongly concentrate particles. Numerical simulations find that solids 10-100 cm in size gravitationally collapse in turbulent disks. We outline areas for progress, including the possibility that still smaller objects self-gravitate.

Following the optical imaging of the exoplanet candidate Fomalhaut b (Fom b), we present a numerical model of how Fomalhaut's debris disk is gravitationally shaped by an interior planet. The model is simple, adaptable to other debris disks, and can be extended to accommodate multiple planets. If Fom b is the dominant perturber of the belt, then to produce the observed disk morphology it must have a mass < 3 Jupiter masses. If the belt and planet orbits are apsidally aligned, our model predicts a planet mass of 0.5 Jupiter masses. The inner edge of the debris disk at 133 AU lies at the periphery of Fom b's chaotic zone, and the mean disk eccentricity of 0.11 is secularly forced by the planet, supporting predictions made prior to the discovery of Fom b. However, previous mass constraints based on disk morphology rely on several oversimplifications. We explain why our constraint is more reliable. It is based on a global model of the disk that is not restricted to the planet's chaotic zone boundary. Moreover, we screen disk parent bodies for dynamical stability over the system age of 100 Myr, and model them separately from their dust grain progeny; the latter's orbits are strongly affected by radiation pressure and their lifetimes are limited to 0.1 Myr by destructive grain-grain collisions. The single planet model predicts that planet and disk orbits be apsidally aligned. Fom b's nominal space velocity does not bear this out, but the astrometric uncertainties may be large. If the apsidal misalignment is real, our upper mass limit of 3 Jupiter masses still holds. The belt contains at least 3 Earth masses of solids that are grinding down to dust. Such a large mass in solids is consistent with Fom b having formed in situ.

We present observations of thermal emission from fifteen transneptunian objects (TNOs) made using the Spitzer Space Telescope. Thirteen of the targets are members of the Classical population: six dynamically hot Classicals, five dynamically cold Classicals, and two dynamically cold inner Classical Kuiper Belt Objects (KBOs). We fit our observations using thermal models to determine the sizes and albedos of our targets finding that the cold Classical TNOs have distinctly higher visual albedos than the hot Classicals and other TNO dynamical classes. The cold Classicals are known to be distinct from other TNOs in terms of their color distribution, size distribution, and binarity fraction. The Classical objects in our sample all have red colors yet they show a diversity of albedos which suggests that there is not a simple relationship between albedo and color. As a consequence of high albedos, the mass estimate of the cold Classical Kuiper Belt is reduced from approximately 0.01 Earth masses to approximately 0.001 Earth masses. Our results also increase significantly the sample of small Classical KBOs with known albedos and sizes from 21 to 32 such objects.

Fomalhaut is a bright star 7.7 parsecs (25 light years) from Earth that harbors a belt of cold dust with a structure consistent with gravitational sculpting by an orbiting planet. Here, we present optical observations of an exoplanet candidate, Fomalhaut b. In the plane of the belt, Fomalhaut b lies approximately 119 astronomical units (AU) from the star and 18 AU from the belt, matching predictions. We detect counterclockwise orbital motion using Hubble Space Telescope observations separated by 1.73 years. Dynamical models of the interaction between the planet and the belt indicate that the planet's mass is at most three times that of Jupiter for the belt to avoid gravitational disruption. The flux detected at 800 nm is also consistent with that of a planet with mass no greater than a few times that of Jupiter. The brightness at 600 nm and the lack of detection at longer wavelengths suggest that the detected flux may include starlight reflected off a circumplanetary disk, with dimension comparable to the orbits of the Galilean satellites. We also observed variability of unknown origin at 600 nm.

Photoionization heating from UV radiation incident on the atmospheres of hot Jupiters may drive planetary mass loss. We construct a model of escape that includes realistic heating and cooling, ionization balance, tidal gravity, and pressure confinement by the host star wind. We show that mass loss takes the form of a hydrodynamic ("Parker") wind, emitted from the planet's dayside during lulls in the stellar wind. When dayside winds are suppressed by the confining action of the stellar wind, nightside winds might pick up if there is sufficient horizontal transport of heat. A hot Jupiter loses mass at maximum rates of ~2 x 10^12 g/s during its host star's pre-main-sequence phase and ~2 x10^10 g/s during the star's main sequence lifetime, for total maximum losses of ~0.06% and ~0.6% of the planet's mass, respectively. For UV fluxes F_UV < 10^4 erg/cm^2/s, the mass loss rate is approximately energy-limited and is proportional to F_UV^0.9. For larger UV fluxes, such as those typical of T Tauri stars, radiative losses and plasma recombination force the mass loss rate to increase more slowly as F_UV^0.6. Dayside winds are quenched during the T Tauri phase because of confinement by overwhelming stellar wind pressure. We conclude that while UV radiation can indeed drive winds from hot Jupiters, such winds cannot significantly alter planetary masses during any evolutionary stage. They can, however, produce observable signatures. Candidates for explaining why the Lyman-alpha photons of HD 209458 are absorbed at Doppler-shifted velocities of +/- 100 km/s include charge-exchange in the shock between the planetary and stellar winds.

By numerically integrating the orbits of the giant planets and of test particles over a period of four billion years, we follow the evolution of the location of the midplane of the Kuiper belt. The Classical Kuiper belt conforms to a warped sheet that precesses with a 1.9 Myr period. The present-day location of the Kuiper belt plane can be computed using linear secular perturbation theory: the local normal to the plane is given by the theory's forced inclination vector, which is specific to every semimajor axis. The Kuiper belt plane does not coincide with the invariable plane, but deviates from it by up to a few degrees in stable zones. For example, at a semimajor axis of 38 AU, the local Kuiper belt plane has an inclination of 1.9 deg and a longitude of ascending node of 149.9 deg when referred to the mean ecliptic and equinox of J2000. At a semimajor axis of 43 AU, the local plane has an inclination of 1.9 deg and a nodal longitude of 78.3 deg. Only at infinite semimajor axis does the Kuiper belt plane merge with the invariable plane, whose inclination is 1.6 deg and nodal longitude is 107.7 deg. A Kuiper belt object keeps its inclination relative to the Kuiper belt plane nearly constant, even while the latter plane departs from the trajectory predicted by linear theory. The constancy of relative inclination reflects the undamped amplitude of free oscillation. Current observations of Classical Kuiper belt objects are consistent with the plane being warped by the giant planets alone, but the sample size will need to increase by a few times before confirmation exceeds 3-sigma in confidence. In principle, differences between the theoretically expected plane and the observed plane could be used to infer as yet unseen masses orbiting the Sun, but carrying out such a program would be challenging.

Gravitational instability of a vertically thin, dusty sheet near the midplane of a protoplanetary disk has long been proposed as a way of forming planetesimals. Before Roche densities can be achieved, however, the dust-rich layer, sandwiched from above and below by more slowly rotating dust-poor gas, threatens to overturn and mix by the Kelvin-Helmholtz instability (KHI). Whether such a threat is real has never been demonstrated: the Richardson criterion for the KHI is derived for 2-D Cartesian shear flow and does not account for rotational forces. Here we present 3-D numerical simulations of gas-dust mixtures in a shearing box, accounting for the full suite of disk-related forces: the Coriolis and centrifugal forces, and radial tidal gravity. Dust particles are assumed small enough to be perfectly entrained in gas; the two fluids share the same velocity field but obey separate continuity equations. We find that the Richardson number Ri does not alone determine stability. The critical value of Ri below which the dust layer overturns and mixes depends on the height-integrated metallicity Z (surface density ratio of dust to gas). Nevertheless, for Z between one and five times solar, the critical Ri is nearly constant at 0.1. Keplerian radial shear stabilizes those modes that would otherwise disrupt the layer at large Ri. If Z is at least 5 times greater than the solar value of 0.01, then midplane dust densities can approach Roche densities. Such an environment might be expected to produce gas giant planets having similarly super-solar metallicities.

The triple nucleus of M31 consists of a population of old red stars in an eccentric disk (P1 and P2) and another population of younger A stars in a circular disk (P3) around M31's central supermassive black hole (SMBH). We argue that P1 and P2 determine the maximal radial extent of the younger A star population and provide the gas that fueled the starburst that generated P3. The eccentric stellar disk creates an $m=1$ non-axisymmetric perturbation to the potential. This perturbed potential drives gas into the inner parsec around the SMBH, if the pattern speed of the eccentric stellar disk is $\Omega_p \lesssim 3-10 {\rm km s^{-1} pc^{-1}}$. We show that stellar mass loss from P1 and P2 is sufficient to create a gravitationally unstable gaseous disk of $\sim 10^5\Msun$ every $0.1-1$ Gyrs, consistent with the 200 Myr age of P3. Similar processes may act in other systems to produce very compact nuclear starbursts.

How do T Tauri disks accrete? The magneto-rotational instability (MRI) supplies one means, but protoplanetary disk gas is typically too poorly ionized to be magnetically active. Here we show that the MRI can, in fact, explain observed accretion rates for the sub-class of T Tauri disks known as transitional systems. Transitional disks are swept clean of dust inside rim radii of ~10 AU. Stellar coronal X-rays ionize material in the disk rim, activating the MRI there. Gas flows from the rim to the star, at a rate limited by the depth to which X-rays ionize the rim wall. The wider the rim, the larger the surface area that the rim wall exposes to X-rays, and the greater the accretion rate. Interior to the rim, the MRI continues to transport gas; the MRI is sustained even at the disk midplane by super-keV X-rays that Compton scatter down from the disk surface. Accretion is therefore steady inside the rim. Blown out by radiation pressure, dust largely fails to accrete with gas. Contrary to what is usually assumed, ambipolar diffusion, not Ohmic dissipation, limits how much gas is MRI-active. We infer values for the transport parameter alpha on the order of 0.01 for GM Aur, TW Hyd, and DM Tau. Because the MRI can only afflict a finite radial column of gas at the rim, disk properties inside the rim are insensitive to those outside. Thus our picture provides one robust setting for planet-disk interaction: a protoplanet interior to the rim will interact with gas whose density, temperature, and transport properties are definite and decoupled from uncertain initial conditions. Our study also supplies half the answer to how disks dissipate: the inner disk drains from the inside out by the MRI, while the outer disk photoevaporates by stellar ultraviolet radiation.

Binaries have played a crucial role many times in the history of modern astronomy and are doing so again in the rapidly evolving exploration of the Kuiper Belt. The large fraction of transneptunian objects that are binary or multiple, 48 such systems are now known, has been an unanticipated windfall. Separations and relative magnitudes measured in discovery images give important information on the statistical properties of the binary population that can be related to competing models of binary formation. Orbits, derived for 13 systems, provide a determination of the system mass. Masses can be used to derive densities and albedos when an independent size measurement is available. Angular momenta and relative sizes of the majority of binaries are consistent with formation by dynamical capture. The small satellites of the largest transneptunian objects, in contrast, are more likely formed from collisions. Correlations of the fraction of binaries with different dynamical populations or with other physical variables have the potential to constrain models of the origin and evolution of the transneptunian population as a whole. Other means of studying binaries have only begun to be exploited, including lightcurve, color, and spectral data. Because of the several channels for obtaining unique physical information, it is already clear that binaries will emerge as one of the most useful tools for unraveling the many complexities of transneptunian space.

As many as 5 ice giants--Neptune-mass planets composed of 90% ice and rock and 10% hydrogen--are thought to form at heliocentric distances of 10-25 AU on closely packed orbits spaced ~5 Hill radii apart. Such oligarchies are ultimately unstable. Once the parent disk of planetesimals is sufficiently depleted, oligarchs perturb one another onto crossing orbits. We explore both the onset and the outcome of the instability through numerical integrations, including dynamical friction cooling of planets by a planetesimal disk whose properties are held fixed. To trigger instability and the ejection of the first ice giant in systems having an original surface density in oligarchs of Sigma ~ 1 g/cm^2, the disk surface density s must fall below 0.1 g/cm^2. Ejections are predominantly by Jupiter and occur within 10 Myr. To eject more than 1 oligarch requires s < 0.03 g/cm^2. Systems starting with up to 4 oligarchs in addition to Jupiter and Saturn can readily yield solar-system-like outcomes in which 2 surviving ice giants lie inside 30 AU and have their orbits circularized by dynamical friction. Our numerical simulations support the idea that planetary systems begin in more crowded and compact configurations, like those of shear-dominated oligarchies. In contrast to previous studies, we identify s < 0.1 Sigma as the regime relevant for understanding the evolution of the outer solar system, and we encourage future studies to concentrate on this regime while relaxing our assumption of a fixed planetesimal disk.

We use spectrally dispersed near-IR interferometry data to constrain the temperature profiles of sub-AU-sized regions of 11 Herbig Ae/Be sources. We find that a single-temperature ring does not reproduce the data well. Rather, models incorporating radial temperature gradients are preferred. These gradients may arise in a dusty disk, or may reflect separate gas and dust components with different temperatures and spatial distributions. Comparison of our models with broadband spectral energy distributions suggests the latter explanation. The data support the view that the near-IR emission of Herbig Ae/Be sources arises from hot circumstellar dust and gas in sub-AU-sized disk regions. Intriguingly, our derived temperature gradients appear systematically steeper for disks around higher mass stars. It is not clear, however, whether this reflects trends in relative dust/gas contributions or gradients within individual components.

AU Microscopii is a 12 Myr old M dwarf that harbors an optically thin, edge-on disk of dust. The scattered light surface brightness falls with projected distance b from the star as b^-a; within b = 43 AU, a = 1-2, while outside 43 AU, a = 4-5. We devise a theory to explain this profile. At a stellocentric distance r = r_BR = 43 AU, we posit a ring of parent bodies on circular orbits: the "birth ring," wherein micron-sized grains are born from the collisional attrition of parent bodies. The "inner disk" at r < r_BR contains grains that migrate inward by corpuscular and Poynting-Robertson (CPR) drag. The "outer disk" at r > r_BR comprises grains just large enough to remain bound to the star, on orbits rendered highly eccentric by stellar wind and radiation pressure. How the vertical optical depth tau scales with r depends on the fraction of grains that migrate inward by CPR drag without suffering a collision. If this fraction is large, the inner disk and birth ring share the same optical depth, and tau scales as r^-5/2 in the outer disk. By contrast, under collision-dominated conditions, the inner disk is empty, and tau scales as r^-3/2 outside. These scaling relations, which we derive analytically and confirm numerically, are robust against uncertainties in the grain size distribution. By simultaneously modeling the surface brightness and thermal spectrum, we break model degeneracies to establish that the AU Mic system is collision-dominated, and that its narrow birth ring contains a lunar mass of decimeter-sized bodies. The inner disk is devoid of micron-sized grains; the surface brightness at b < 43 AU arises from light forward scattered by the birth ring. Inside b = 43 AU, the disk's V-H color should not vary with b; outside, the disk must become bluer as ever smaller grains are probed.

The Edgeworth-Kuiper belt encodes the dynamical history of the outer solar system. Kuiper belt objects (KBOs) bear witness to coagulation physics, the evolution of planetary orbits, and external perturbations from the solar neighborhood. We critically review the present-day belt's observed properties and the theories designed to explain them. Theories are organized according to a possible time-line of events. In chronological order, epochs described include (1) coagulation of KBOs in a dynamically cold disk, (2) formation of binary KBOs by fragmentary collisions and gravitational captures, (3) stirring of KBOs by Neptune-mass planets (``oligarchs''), (4) eviction of excess oligarchs, (5) continued stirring of KBOs by remaining planets whose orbits circularize by dynamical friction, (6) planetary migration and capture of Resonant KBOs, (7) creation of the inner Oort cloud by passing stars in an open stellar cluster, and (8) collisional comminution of the smallest KBOs. Recent work underscores how small, collisional, primordial planetesimals having low velocity dispersion permit the rapid assembly of ~5 Neptune-mass oligarchs at distances of 15-25 AU. We explore the consequences of such a picture. We propose that Neptune-mass planets whose orbits cross into the Kuiper belt for up to ~20 Myr help generate the high-perihelion members of the hot Classical disk and Scattered belt. By contrast, raising perihelia by sweeping secular resonances during Neptune's migration might fill these reservoirs too inefficiently when account is made of how little primordial mass might reside in bodies having sizes of order 100 km. These and other frontier issues in trans-Neptunian space are discussed quantitatively.

We present a new, simple, fast algorithm to numerically evolve disks of inelastically colliding particles surrounding a central star. Our algorithm adds negligible computational cost to the fastest existing collisionless N-body codes, and can be used to simulate, for the first time, the interaction of planets with disks over many viscous times. Though the algorithm is implemented in two dimensions-i.e., the motions of bodies need only be tracked in a plane-it captures the behavior of fully three-dimensional disks in which collisions maintain inclinations that are comparable to random eccentricities. We subject the algorithm to a battery of tests for the case of an isolated, narrow, circular ring. Numerical simulations agree with analytic theory with regards to how particles' random velocities equilibrate; how the ring viscously spreads; and how energy dissipation, angular momentum transport, and material transport are connected. We derive and measure the critical value of the coefficient of restitution above which viscous stirring dominates inelastic damping and the particles' velocity dispersion runs away.

A residual planetesimal disk of mass 10-100 Earth masses remained in the outer solar system following the birth of the giant planets, as implied by the existence of the Oort cloud, coagulation requirements for Pluto, and inefficiencies in planet formation. Upon gravitationally scattering planetesimal debris, planets migrate. Orbital migration can lead to resonance capture, as evidenced here in the Kuiper and asteroid belts, and abroad in extra-solar systems. Finite sizes of planetesimals render migration stochastic ("noisy"). At fixed disk mass, larger (fewer) planetesimals generate more noise. Extreme noise defeats resonance capture. We employ order-of-magnitude physics to construct an analytic theory for how a planet's orbital semi-major axis fluctuates in response to random planetesimal scatterings. To retain a body in resonance, the planet's semi-major axis must not random walk a distance greater than the resonant libration width. We translate this criterion into an analytic formula for the retention efficiency of the resonance as a function of system parameters, including planetesimal size. We verify our results with tailored numerical simulations. Application of our theory reveals that capture of Resonant Kuiper belt objects by a migrating Neptune remains effective if the bulk of the primordial disk was locked in bodies having sizes < O(100) km and if the fraction of disk mass in objects with sizes > 1000 km was less than a few percent. Coagulation simulations produce a size distribution of primordial planetesimals that easily satisfies these constraints. We conclude that stochasticity did not interfere with, nor modify in any substantive way, Neptune's ability to capture and retain Resonant Kuiper belt objects during its migration.

In a protoplanetary disk, the inner edge of the region where the temperature falls below the condensation temperature of water is referred to as the 'snow line'. Outside the snow line, water ice increases the surface density of solids by a factor of 4. The mass of the fastest growing planetesimal (the 'isolation mass') scales as the surface density to the 3/2 power. It is thought that ice-enhanced surface densities are required to make the cores of the gas giants (Jupiter and Saturn) before the disk gas dissipates. Observations of the Solar System's asteroid belt suggest that the snow line occurred near 2.7 AU. In this paper we revisit the theoretical determination of the snow line. In a minimum-mass disk characterized by conventional opacities and a mass accretion rate of 10^-8 solar masses per year, the snow line lies at 1.6-1.8 AU, just past the orbit of Mars. The minimum-mass disk, with a mass of 0.02 solar, has a life time of 2 million years with the assumed accretion rate. Moving the snow line past 2.7 AU requires that we increase the disk opacity, accretion rate, and/or disk mass by factors ranging up to an order of magnitude above our assumed baseline values.

We present Keck Interferometer observations of TW Hya that spatially resolve its emission at 2 micron wavelength. Analyzing these data together with existing K-band veiling and near-infrared photometric measurements, we conclude that the inner disk consists of optically thin, sub-micron-sized dust extending from ~4 AU to within 0.06 AU of the central star. The inner disk edge may be magnetospherically truncated. Even if we account for the presence of gas in the inner disk, these small dust grains have survival times against radiation blow-out that are orders of magnitude shorter than the age of the system, suggesting continual replenishment through collisions of larger bodies.

The problem of accretion in the Trojan 1:1 resonance is akin to the standard problem of planet formation, transplanted from a star-centered disk to a disk centered on the Lagrange point. The newly discovered class of Neptune Trojans promises to test theories of planet formation by coagulation. Neptune Trojans resembling the prototype 2001 QR322 (``QR'')--whose radius of ~100 km is comparable to that of the largest Jupiter Trojan--may outnumber their Jovian counterparts by a factor of ~10. We discover that seeding the 1:1 resonance with debris from planetesimal collisions and having the seed particles accrete in situ naturally reproduces the inferred number of QR-sized Trojans. We analyze accretion in the Trojan sub-disk by applying the two-groups method, accounting for kinematics specific to the resonance. We find that a Trojan sub-disk comprising decimeter-sized seed particles and having a surface density 1e-3 that of the local minimum-mass disk produces ~10 QR-sized objects in ~1 Gyr, in accord with observation. Further growth is halted by collisional diffusion of seed particles out of resonance. In our picture, the number and sizes of the largest Neptune Trojans represent the unadulterated outcome of dispersion-dominated, oligarchic accretion. Large Neptune Trojans, perhaps the most newly accreted objects in our Solar System, may today have a dispersion in orbital inclination of less than ~10 degrees, despite the existence of niches of stability at higher inclinations. Such a vertically thin disk, born of a dynamically cold environment necessary for accretion, and raised in minimal contact with external perturbations, contrasts with the thick disks of other minor body belts.

The spatial distribution of Kuiper belt objects (KBOs) in 2:1 exterior resonance with Neptune constrains that planet's migration history. Numerical simulations demonstrate that fast planetary migration generates a larger population of KBOs trailing rather than leading Neptune in orbital longitude. This asymmetry corresponds to a greater proportion of objects caught into asymmetric resonance such that their resonance angles librate about values > pi (trailing) as opposed to < pi (leading). We provide, for the first time, an explanation of this phenomenon, using physical, analytic, and semi-analytic arguments. Central to our understanding is how planetary migration shifts the equilibrium points of the superposed direct and indirect potentials. Large capture asymmetries appear for exponential timescales of migration shorter than 10 Myr. The observed distribution of 2:1 KBOs (2 trailing and 7 leading) excludes migration timescales < 1 Myr with 99.65% confidence.

The light curves of the pre-main-sequence star KH 15D from the years 1913--2003 can be understood if the star is a member of an eccentric binary that is encircled by a vertically thin, inclined ring of dusty gas. Eclipses occur whenever the reflex motion of a star carries it behind the circumbinary ring; the eclipses occur with period equal to the binary orbital period of 48.4 days. Features of the light curve--including the amplitude of central reversals during mid-eclipse, the phase of eclipse with respect to the binary orbit phase, the level of brightness out-of-eclipse, the depth of eclipse, and the eclipse duty cycle--are all modulated on the timescale of nodal regression of the obscuring ring, in accord with the historical data. The ring has a mean radius near 3 AU and a radial width that is likely less than this value. While the inner boundary could be shepherded by the central binary, the outer boundary may require an exterior planet to confine it against viscous spreading. The ring must be vertically warped to maintain a non-zero inclination. Thermal pressure gradients and/or ring self-gravity can readily enforce rigid precession. In coming years, as the node of the ring regresses out of our line-of-sight toward the binary, the light curve from the system should cycle approximately back through its previous behavior. Near-term observations should seek to detect a mid-infrared excess from this system; we estimate the flux densities from the ring to be 3 mJy at wavelengths of 10--100 microns.

Solid particles in protoplanetary disks that are sufficiently super-solar in metallicity overcome turbulence generated by vertical shear to gravitationally condense into planetesimals. Super-solar metallicities result if solid particles pile up as they migrate starward as a result of aerodynamic drag. Previous analyses of aerodynamic drift rates that account for mean flow differences between gas and particles yield particle pile-ups. We improve on these studies not only by accounting for the collective inertia of solids relative to that of gas, but also by including the transport of angular momentum by turbulent stresses within the particle layer. These turbulent stresses are derived in a physically self-consistent manner from the structure of marginally Kelvin-Helmholtz turbulent flows. They are not calculated using the usual plate drag formulae, whose use we explain is inappropriate. Accounting for the relative inertia of solids to gas retards, but does not prevent, particle pile-ups, and generates more spatially extended regions of metal enrichment. Turbulent transport hastens pile-ups. We conclude that particle pile-up is a robust outcome in sufficiently passive protoplanetary disks. Connections to observations of circumstellar disks, including the Kuiper Belt, and the architectures of planetary systems are made.

We review ongoing efforts to identify occupants of mean-motion resonances (MMRs) and collisional families in the Edgeworth-Kuiper belt. Direct integrations of trajectories of Kuiper belt objects (KBOs) reveal the 1:1 (Trojan), 5:4, 4:3, 3:2 (Plutino), 5:3, 7:4, 9:5, 2:1 (Twotino), and 5:2 MMRs to be inhabited. Apart from the Trojan, resonant KBOs typically have large orbital eccentricities and inclinations. The observed pattern of resonance occupation is consistent with resonant capture and adiabatic excitation by a migratory Neptune; however, the dynamically cold initial conditions prior to resonance sweeping that are typically assumed by migration simulations are probably inadequate. Given the dynamically hot residents of the 5:2 MMR and the substantial inclinations observed in all exterior MMRs, a fraction of the primordial belt was likely dynamically pre-heated prior to resonance sweeping. A pre-heated population may have arisen as Neptune gravitationally scattered objects into trans-Neptunian space. The spatial distribution of Twotinos offers a unique diagnostic of Neptune's migration history. The Neptunian Trojan population may rival the Jovian Trojan population, and the former's existence is argued to rule out violent orbital histories for Neptune. Finally, lowest-order secular theory is applied to several hundred non-resonant KBOs with well-measured orbits to update proposals of collisional families. No convincing family is detected.

Narrow planetary rings are eccentric and inclined. Particles within a given ring must therefore share the same pericenter and node. We solve for the three-dimensional geometries and mass distributions that enable the Uranian Alpha and Beta rings, and the Saturnian Maxwell and Colombo (Titan) rings, to maintain simultaneous apsidal and nodal lock. Ring self-gravity, interparticle collisions, and the quadrupole field of the host planet balance each other to achieve this equilibrium. We prove that such an equilibrium is linearly stable. Predictions for the Saturnian ringlets to be tested by the Cassini spacecraft include: (1) ringlet masses are of order 1e19 g, (2) surface mass densities should increase from ring midline to ring edges, and (3) rings are vertically warped such that the fractional variation of inclination across the ring is of order 10%. Analogous predictions are made for the Uranian rings. Simultaneous apsidal and nodal locking forces the narrowest portion of the ring--its ``pinch,'' where self-gravitational and collisional forces are strongest--to circulate relative to the node, and introduces previously unrecognized time-varying forces perpendicular to the planet's equator plane. We speculate that such periodic stressing might drive kilometer-scale bending waves at a frequency twice that of apsidal precession. Such flexing might be observed over a few weeks by Cassini.

We critically examine the best lines of evidence for grain growth in protoplanetary disks, based on modelling of observed spectral energy distributions and images of T Tauri and Herbig Ae stars. The data are consistent with millimeter-sized grains near disk midplanes, and micron-sized grains near disk surfaces. We review three channels by which grains can grow, including direct condensation from the vapor phase, grain-grain collisional sticking, and gravitational instability. The utility of dust in identifying as yet unseen extrasolar planets is highlighted.

The Deep Ecliptic Survey is a project whose goal is to survey a large area of the near-ecliptic region to a faint limiting magnitude (R ~ 24) in search of objects in the outer solar system. We are collecting a large homogeneous data sample from the Kitt Peak Mayall 4-m and Cerro Tololo Blanco 4-m telescopes with the Mosaic prime-focus CCD cameras. Our goal is to collect a sample of 500 objects with good orbits to further our understanding of the dynamical structure of the outer solar system. This survey has been in progress since 1998 and is responsible for 272 designated discoveries as of March 2003. We summarize our techniques, highlight recent results, and describe publically available resources.

As part of our ongoing Deep Ecliptic Survey (DES) of the Kuiper belt, we report on the occupation of the 1:1 (Trojan), 4:3, 3:2, 7:4, 2:1, and 5:2 Neptunian mean-motion resonances (MMRs). The occupation of the 1:1 and 5:2 MMRs is not easily understood within the standard model of resonance sweeping by a migratory Neptune over an initially dynamically cold belt. Our dynamically hot, 5:2 resonant objects can librate with modest amplitudes of 90 deg within the resonance for at least 1 Gyr. Their trajectories cannot be explained by close encounters with Neptune alone, given the latter's current orbit. The hot orbits of such 5:2 resonant KBOs, unlike hot orbits of previously known resonant KBOs, may imply that these objects were pre-heated to large inclination and large eccentricity prior to resonance capture by a migratory Neptune. Our first discovered Neptunian Trojan, 2001QR322, may not owe its existence to Neptune's migration at all. The trajectory of 2001QR322 is remarkably stable; the object can undergo tadpole-type libration about Neptune's leading Lagrange (L4) point for at least 1 Gyr with a libration amplitude of 24 deg. Trojan capture probably occurred while Neptune accreted the bulk of its mass. For an assumed albedo of 12--4%, our Trojan is 130--230 km in diameter. Model-dependent estimates place the total number of Neptune Trojans resembling 2001QR322 at 20--60. Their existence might rule out violent orbital histories for Neptune.

We review the orbital dynamics exhibited by the first extra-solar planetary system discovered, Upsilon Andromedae. This system is unique in combining all of the surprising architectural features displayed individually by extrasolar planetary systems found today: (1) a hot Jupiter, (2) two planets on highly eccentric orbits, and (3) a stellar companion. We discuss the system's stability properties and its possible origin. Planet-disk interactions seem critical to the emerging story.

We illuminate dynamical properties of Kuiper Belt Objects (KBOs) in the 3:2 (``Plutino'') and 2:1 (``Twotino'') Neptunian resonances within the model of resonant capture and migration. We analyze a series of numerical integrations, each involving the 4 migratory giant planets and 400 test particles distributed throughout trans-Neptunian space, to measure efficiencies of capture as functions of migration speed. Snapshots of the spatial distribution of resonant KBOs reveal that Twotinos cluster +/- 75 degrees away from Neptune's longitude, while Plutinos cluster +/- 90 degrees away. Longitudinal clustering persists even for surveys that are not volume-limited in their ability to detect resonant KBOs. Remarkably, between -90 degrees and -60 degrees of Neptune's longitude, we find the sky density of Twotinos to nearly equal that of Plutinos, despite the greater average distance of Twotinos. We couple our findings to observations to crudely estimate that the intrinsic Twotino population is within a factor of 3 of the Plutino population. Most strikingly, the migration model predicts that more Twotinos may lie at longitudes behind that of Neptune than ahead of it. The magnitude of the asymmetry amplifies dramatically with faster rates of migration and can be as large as 300%. A differential measurement of the sky density of 2:1 resonant objects behind of and in front of Neptune's longitude would powerfully constrain the migration history of that planet.

Divergent migration of planets within a viscous circumstellar disk can engender resonance crossings and dramatic excitation of orbital eccentricities. We provide quantitative criteria for the viability of this mechanism. For the orbits of two bodies to diverge, a ring of viscous material must be shepherded between them. As the ring diffuses in radius by virtue of its intrinsic viscosity, the two planets are wedged further apart. The ring mass must be smaller than the planetary masses so that the crossing of an individual resonance lasts longer than the resonant libration period. At the same time, the crossing cannot be of such long duration that the disk's direct influence on the bodies' eccentricities interferes with the resonant interaction between the two planets. This last criterion is robustly satisfied because resonant widths are typically tiny fractions of the orbital radius. We evaluate our criteria not only for giant planets within gaseous protoplanetary disks, but also for shepherd moons that bracket narrow planetary rings in the solar system. A shepherded ring of gas orbiting at a distance of 1 AU from a solar-type star and having a surface density of less than 500 g/cm^2, a dimensionless alpha viscosity of 0.1, and a height-to-radius aspect ratio of 0.05 can drive two Jovian-mass planets through the 2:1 and higher-order resonances so that their eccentricities amplify to values of several tenths. Because of the requirement that the disk mass in the vicinity of the planets be smaller than the planet masses, divergent resonance crossings may figure significantly into the orbital evolution of planets during the later stages of disk evolution, including the debris disk phase.

The dynamical evolution of Classical Kuiper Belt Objects (CKBOs) divides into two parts, according to the secular theory of test particle orbits. The first part is a forced oscillation driven by the planets, while the second part is a free oscillation whose amplitude is determined by the initial orbit of the test particle. We extract the free orbital inclinations and free orbital eccentricities from the osculating elements of 125 known CKBOs. The free inclinations of 32 CKBOs strongly cluster about 2 degrees at orbital semi-major axes between 44 and 45 AU. We propose that these objects comprise a collisional family, the first so identified in the Kuiper Belt. Members of this family are plausibly the fragments of an ancient parent body having a minimum diameter of \~800 km. This body was disrupted upon colliding with a comparably sized object, and generated ejecta having similar free inclinations. Our candidate family is dynamically akin to a sub-family of Koronis asteroids located at semi-major axes less than 2.91 AU; both families exhibit a wider range in free eccentricity than in free inclination, implying that the relative velocity between parent and projectile prior to impact lay mostly in the invariable plane of the solar system. We urge more discoveries of new CKBOs to test the reality of our candidate family and physical studies of candidate family members to probe the heretofore unseen interior of a massive, primitive planetesimal.

Beneath the fusion-encrusted surfaces of the most primitive stony meteorites lies not homogeneous rock, but a profusion of millimeter-sized igneous spheres. These chondrules, and their centimeter-sized counterparts, the calcium-aluminum-rich inclusions, comprise more than half of the volume fraction of chondritic meteorites. They are the oldest creations of the solar system. Their chemical composition matches that of the solar photosphere in all but the most volatile of elements, reflecting their condensation from the same pristine gas that formed the sun. In this invited editorial, we review the nebular shock wave model of Desch and Connolly (Meteoritics and Planetary Science 2002, 37, 183) that seeks to explain their origin. While the model succeeds in reproducing the unique petrological signatures of chondrules, the origin of the required shock waves in protoplanetary disks remains a mystery. Outstanding questions are summarized, with attention paid briefly to competing models.

The orbits of the outer two known planets orbiting Upsilon Andromedae are remarkably eccentric. Planet C possesses an orbital eccentricity of e1 = 0.253. For the more distant planet D, e2 = 0.308. Previous dynamical analyses strongly suggest that the two orbits are nearly co-planar and are trapped in an apsidal resonance in which the difference between their longitudes of periastron undergoes a bounded oscillation about 0 degrees. Here we elucidate the origin of these large eccentricities and of the apsidal alignment. Resonant interactions between a remnant circumstellar disk of gas lying exterior to the orbits of both planets can smoothly grow e2. Secular interactions between planets D and C can siphon off the eccentricity of the former to grow that of the latter. Externally amplifying e2 during the phase of the apsidal oscillation when e2/e1 is smallest drives the oscillation amplitude towards zero. Thus, the substantial eccentricity of planet C and the locking of orbital apsides are both consequences of externally pumping the eccentricity of planet D over timescales exceeding apsidal precession periods of order 1e4 yr. We explain why the recently detected stellar companion to Upsilon Andromedae is largely dynamically decoupled from the planetary system.

Orbits of known extrasolar planets that are located outside the tidal circularization regions of their parent stars are often substantially eccentric. By contrast, planetary orbits in our Solar System are approximately circular, reflecting planet formation within a nearly axisymmetric, circumsolar disk. We propose that orbital eccentricities may be generated by divergent orbital migration of two planets in a viscously accreting circumstellar disk. The migration is divergent in the sense that the ratio of the orbital period of the outer planet to that of the inner planet grows. As the period ratio diverges, the planets traverse, but are not captured into, a series of mean-motion resonances that amplify their orbital eccentricities in rough inverse proportion to their masses. Strong viscosity gradients in protoplanetary disks offer a way to reconcile the circular orbits of Solar System gas giants with the eccentric orbits of currently known extrasolar planets.

The face-on disk around TW Hya is imaged in scattered light at wavelengths of 1.1 and 1.6 micron using the coronagraph in the Near Infrared Camera and Multi Object Spectrometer aboard the Hubble Space Telescope. Stellar light scattered from the optically thick dust disk is seen from 20-230 AU. The surface brightness declines as a power law of r^(-2.6+/-0.1) between 45 and 150 AU. The scattering profile indicates that the disk is flared, not geometrically flat. The disk, while spatially unresolved in thermal radiation at wavelengths of 12 and 18 micron in observations from the W. M. Keck Observatory, shows amorphous and crystalline silicate emission in its spectrum. A disk with silicate grains of a ~1 micron in size in its surface layers can explain the shape of the mid-infrared spectrum. Much larger grains in the disk interior are necessary to fit the millimeter-wave spectral energy distribution, and hence grain growth from an original interstellar size population may have occurred.

One of the parameters fitted by Doppler radial velocity measurements of extrasolar planetary systems is omega, the argument of pericenter of a given planet's orbit referenced to the plane of the sky. Curiously, the omega's of the outer two planets orbiting Upsilon Andromedae are presently nearly identical: Delta-omega = omega_D - omega_C = 4.8 deg +/- 4.8 deg (1 sigma). This observation is least surprising if planets C and D occupy orbits that are seen close to edge-on (sin i_C, sin i_D > 0.5) and whose mutual inclination Theta does not exceed 20 deg. In this case, planets C and D inhabit a secular resonance in which Delta-omega librates about 0 deg with an amplitude of 30 deg and a period of 4000 yr. The resonant configuration spends about one-third of its time with |Delta-omega| < 10 deg. If Theta > 40 deg, either Delta-omega circulates or the system is unstable. This instability is driven by the Kozai mechanism which couples the eccentricity of planet C to Theta to drive the former quantity to values approaching unity. Our expectation that Theta < 20 deg suggests that planets C and D formed in a flattened, circumstellar disk, and may be tested by upcoming astrometric measurements with the FAME satellite.

We improve upon the radiative, hydrostatic equilibrium models of passive circumstellar disks constructed by Chiang & Goldreich (1997). New features include (1) account for a range of particle sizes, (2) employment of laboratory-based optical constants of representative grain materials, and (3) numerical solution of the equations of radiative and hydrostatic equilibrium within the original 2-layer (disk surface + disk interior) approximation. We explore how the spectral energy distribution (SED) of a face-on disk depends on grain size distributions, disk geometries and surface densities, and stellar photospheric temperatures. Observed SEDs of 3 Herbig Ae and 2 T Tauri stars, including spectra from the Long Wavelength Spectrometer (LWS) aboard the Infrared Space Observatory (ISO), are fitted with our models. Silicate emission bands from optically thin, superheated disk surface layers appear in nearly all systems. Water ice emission bands appear in LWS spectra of 2 of the coolest stars. Infrared excesses in several sources are consistent with vertical settling of photospheric grains. While this work furnishes further evidence that passive reprocessing of starlight by flared disks adequately explains the origin of infrared-to-millimeter wavelength excesses of young stars, we emphasize how the SED alone does not provide sufficient information to constrain particle sizes and disk masses uniquely.

We examine the angular momentum transport properties of disks composed of macroscopic particles whose velocity dispersions are externally enhanced (``stirred''). Our simple Boltzmann equation model serves as an analogy for unmagnetized fluid disks in which turbulence may be driven by thermal convection. We show that interparticle collisions in particle disks play the same role as fluctuating pressure forces and viscous dissipation in turbulent disks: both transfer energy in random motions associated with one direction to those associated with another, and convert kinetic energy into heat. The direction of angular momentum transport in stirred particle and fluid disks is determined by the direction of external stirring and by the properties of the collision term in the Boltzmann equation (or its analogue in the fluid problem). In particular, our model problem yields inward transport for vertically or radially stirred disks, provided collisions are suitably inelastic; the transport is outwards in the elastic limit. Numerical simulations of hydrodynamic turbulence driven by thermal convection find inward transport; this requires that fluctuating pressure forces do little to no work, and is analogous to an externally stirred particle disk in which collisions are highly inelastic.

The boundaries of the Uranian epsilon, alpha, and beta rings can be fitted by Keplerian ellipses. The pair of ellipses that outline a given ring share a common line of apsides. Apse alignment is surprising because the quadrupole moment of Uranus induces differential precession. We propose that rigid precession is maintained by a balance of forces due to ring self-gravity, planetary oblateness, and interparticle collisions. Collisional impulses play an especially dramatic role near ring edges. Pressure-induced accelerations are maximal near edges because there (1) velocity dispersions are enhanced by resonant satellite perturbations, and (2) the surface density declines steeply. Remarkably, collisional forces felt by material in the last 100 m of a 10 km wide ring can increase equilibrium masses up to a factor of 100. New ring surface densities are derived which accord with Voyager radio measurements. In contrast to previous models, collisionally modified self-gravity appears to allow for both negative and positive eccentricity gradients; why all narrow planetary rings exhibit positive eccentricity gradients remains an open question.

We present the results of a pencil-beam survey of the Kuiper Belt using the Keck 10-m telescope. A single 0.01 square degree field is imaged 29 times for a total integration time of 4.8 hr. Combining exposures in software allows the detection of Kuiper Belt Objects (KBOs) having visual magnitude V < 27.9. Two new KBOs are discovered. One object having V = 25.5 lies at a probable heliocentric distance d = 33 AU. The second object at V = 27.2 is located at d = 44 AU. Both KBOs have diameters of about 50 km, assuming comet-like albedos of 4%. Data from all surveys are pooled to construct the luminosity function from red magnitude R = 20 to 27. The cumulative number of objects per square degree, N (< R), is fitted to a power law of the form log_(10) N = 0.52 (R - 23.5). Differences between power laws reported in the literature are due mainly to which survey data are incorporated, and not to the method of fitting. The luminosity function is consistent with a power-law size distribution for objects having diameters s = 50 to 500 km; dn ~ s^(-q) ds, where the differential size index q = 3.6 +/- 0.1. The distribution is such that the smallest objects possess most of the surface area, but the largest bodies contain the bulk of the mass. Though our inferred size index nearly matches that derived by Dohnanyi (1969), it is unknown whether catastrophic collisions are responsible for shaping the size distribution. Implications of the absence of detections of classical KBOs beyond 50 AU are discussed.

We compute spectral energy distributions (SEDs) for passive T Tauri disks viewed at arbitrary inclinations. Semi-analytic models of disks in radiative and hydrostatic equilibrium are employed. Over viewing angles for which the flared disk does not occult the central star, the SED varies negligibly with inclination. For such aspects, the SED shortward of ~80 microns is particularly insensitive to orientation, since short wavelength disk emission is dominated by superheated surface layers which are optically thin. The SED of a nearly edge-on disk is that of a class I source. The outer disk occults inner disk regions, and emission shortward of ~30 microns is dramatically extinguished. Spectral features from dust grains may appear in absorption. However, millimeter wavelength fluxes decrease by at most a factor of 2 from face-on to edge-on orientations. We present illustrative applications of our SED models. The class I source 04108+2803B is considered a T Tauri star hidden from view by an inclined circumstellar disk. Fits to its observed SED yield model-dependent values for the disk mass of ~0.015 solar masses and a disk inclination of ~65 degrees relative to face-on. The class II source GM Aur represents a T Tauri star unobscured by its circumstellar disk. Fitted parameters include a disk mass of \~0.050 solar masses and an inclination of ~60 degrees.

We derive hydrostatic, radiative equilibrium models for passive disks surrounding T Tauri stars. Each disk is encased by an optically thin layer of superheated dust grains. This layer re-emits directly to space about half the stellar energy it absorbs. The other half is emitted inward and regulates the interior temperature of the disk. The heated disk flares. As a consequence, it absorbs more stellar radiation, especially at large radii, than a flat disk would. The portion of the spectral energy distribution contributed by the disk is fairly flat throughout the thermal infrared. At fixed frequency, the contribution from the surface layer exceeds that from the interior by about a factor 3 and is emitted at more than an order of magnitude greater radius. Spectral features from dust grains in the superheated layer appear in emission if the disk is viewed nearly face-on.