In this work, we investigate configuration formation of two inner terrestrial planets near mean motion resonance (MMRs) induced by the perturbation of a distant gas-giant for the Kepler-68 system, by conducting thousands of numerical simulations. The results show that the formation of terrestrial planets is relevant to the speed of Type I migration, the mass of planets, and the existence of giant planet. The mass and eccentricity of the giant planet may play a crucial role in shaping the final configuration of the system. The inner planet pair can be trapped in 5:3 or 7:4 MMRs if the giant planet revolves the central star with an eccentric orbit, which is similar to the observed configuration of Kepler-68. Moreover, we find that the eccentricity of the middle planet can be excited to roughly 0.2 if the giant planet is more massive than 5 $M_J$,otherwise the terrestrial planets are inclined to remain near-circular orbits. Our study may provide a likely formation scenario for the planetary systems that harbor several terrestrial planets near MMRs inside and one gas-giant exterior to them.

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.

We present ALMA mosaic observations at 1.3 mm (223 GHz) of the Fomalhaut system with a sensitivity of 14 $\mu$Jy/beam. These observations provide the first millimeter map of the continuum dust emission from the complete outer debris disk with uniform sensitivity, enabling the first conclusive detection of apocenter glow. We adopt a MCMC modeling approach that accounts for the eccentric orbital parameters of a collection of particles within the disk. The outer belt is radially confined with an inner edge of $136.3\pm0.9$ AU and width of $13.5\pm1.8$ AU. We determine a best-fit eccentricity of $0.12\pm0.01$. Assuming a size distribution power law index of $q=3.46\pm 0.09$, we constrain the dust absorptivity power law index $\beta$ to be $0.9<\beta<1.5$. The geometry of the disk is robustly constrained with inclination $65.\!\!^\circ6\pm0.\!\!^\circ3$, position angle $337.\!\!^\circ9\pm0.\!\!^\circ3$, and argument of periastron $22.\!\!^\circ5\pm4.\!\!^\circ3$. Our observations do not confirm any of the azimuthal features found in previous imaging studies of the disk with HST, SCUBA, and ALMA. However, we cannot rule out structures $\leq10$ AU in size or which only affect smaller grains. The central star is clearly detected with a flux density of $0.75\pm0.02$ mJy, significantly lower than predicted by current photospheric models. We discuss the implications of these observations for the directly imaged Fomalhaut b and the inner dust belt detected at infrared wavelengths.

Recent ALMA observations present mounting evidence for the presence of exocometary gas released within Kuiper belt analogues around nearby main sequence stars. This represents a unique opportunity to study their ice reservoir at the younger ages when volatile delivery to planets is most likely to occur. We here present the detection of CO J=2-1 emission co-located with dust emission from the cometary belt in the 440 Myr-old Fomalhaut system. Through spectro-spatial filtering, we achieve a 5.4$\sigma$ detection and determine that the ring's sky-projected rotation axis matches that of the star. The CO mass derived ($0.65-42 \times10^{-7}$ M$_{\oplus}$) is the lowest of any circumstellar disk detected to date, and must be of exocometary origin. Using a steady state model, we estimate the CO+CO$_2$ mass fraction of exocomets around Fomalhaut to be between 4.6-76%, consistent with Solar System comets and the two other belts known to host exocometary gas. This is the first indication of a similarity in cometary compositions across planetary systems that may be linked to their formation scenario and is consistent with direct ISM inheritance. In addition, we find tentative evidence that $(49\pm 27)$% of the detected flux originates from a region near the eccentric belt's pericentre. If confirmed, the latter may be explained through a recent impact event or CO pericentre glow due to exocometary release within a steady state collisional cascade. In the latter scenario, we show how the azimuthal dependence of the CO release rate leads to asymmetries in gas observations of eccentric exocometary belts.

A commonly noted feature of the population of multi-planet extrasolar systems is the rarity of planet pairs in low-order mean-motion resonances. We revisit the physics of resonance capture via convergent disk-driven migration. We point out that for planet spacings typical of stable configurations for Kepler systems, the planets can routinely maintain a small but nonzero eccentricity due to gravitational perturbations from their neighbors. Together with the upper limit on the migration rate needed for capture, the finite eccentricity can make resonance capture difficult or impossible in Sun-like systems for planets smaller than ~Neptune-sized. This mass limit on efficient capture is broadly consistent with observed exoplanet pairs that have mass determinations: of pairs with the heavier planet exterior to the lighter planet -- which would have been undergoing convergent migration in their disks -- those in or nearly in resonance are much more likely to have total mass greater than two Neptune masses than to have smaller masses. The agreement suggests that the observed paucity of resonant pairs around sun-like stars may simply arise from a small resonance capture probability for lower-mass planets. Planet pairs that thereby avoid resonance capture are much less likely to collide in an eventual close approach than to simply migrate past one another to become a divergently migrating pair with the lighter planet exterior. For systems around M stars we expect resonant pairs to be much more common, since there the minimum mass threshhold for efficient capture is about an Earth mass.

Debris disks often take the form of eccentric rings with azimuthal asymmetries in surface brightness. Such disks are often described as showing "pericenter glow", an enhancement of the disk brightness in regions nearest the central star. At long wavelengths, however, the disk apocenters should appear brighter than their pericenters: in the long wavelength limit, we find the apocenter/pericenter flux ratio scales as 1+e for disk eccentricity e. We produce new models of this "apocenter glow" to explore its causes and wavelength dependence and study its potential as a probe of dust grain properties. Based on our models, we argue that several far-infrared and (sub)millimeter images of the Fomalhaut and epsilon Eridani debris rings obtained with Herschel, JCMT, SHARC II, ALMA, and ATCA should be reinterpreted as suggestions or examples of apocenter glow. This reinterpretation yields new constraints on the disks' dust grain properties and size distributions.

Observations in 2013 and 2014 of the Centaur 10199 Chariklo and its ring system consistently indicated that the radial width of the inner, more massive ring varies with longitude. That strongly suggests that this ring has a finite eccentricity despite the fast differential precession that Chariklo's large quadrupole moment should induce. If the inferred apse alignment is maintained by the ring's self-gravity, as it is for the Uranian rings, we estimate a ring mass of a few times 10^16 g and a typical particle size of a few meters. These imply a short collisional spreading time of ~10^5 years, somewhat shorter than the typical Centaur dynamical lifetime of a few Myrs and much shorter than the age of the solar system. In light of this time constraint, we evaluate previously suggested ring formation pathways including collisional ejection and satellite disruption. We also investigate in detail a contrasting formation mechanism, the lofting of dust particles off Chariklo's surface into orbit via outflows of sublimating CO and/or N_2 triggered after Chariklo was scattered inward by giant planets. This latter scenario predicts that rings should be common among 100-km class Centaurs but rare among Kuiper belt objects and smaller Centaurs. It also predicts that Centaurs should show seasonal variations in cometary activity with activity maxima occurring shortly after equinox.

We present the Superparticle Model/Algorithm for Collisions in Kuiper belts and debris disks (SMACK), a new method for simultaneously modeling, in 3-D, the collisional and dynamical evolution of planetesimals in a debris disk with planets. SMACK can simulate azimuthal asymmetries and how these asymmetries evolve over time. We show that SMACK is stable to numerical viscosity and numerical heating over 10^7 yr, and that it can reproduce analytic models of disk evolution. We use SMACK to model the evolution of a debris ring containing a planet on an eccentric orbit. Differential precession creates a spiral structure as the ring evolves, but collisions subsequently break up the spiral, leaving a narrower eccentric ring.

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.

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

The standard theoretical treatment of collisional cascades derives a steady-state size distribution assuming a single constant velocity dispersion for all bodies regardless of size. Here we relax this assumption and solve self-consistently for the bodies' steady-state size and size-dependent velocity distributions. Specifically, we account for viscous stirring, dynamical friction, and collisional damping of the bodies' random velocities in addition to the mass conservation requirement typically applied to find the size distribution in a steady-state cascade. The resulting size distributions are significantly steeper than those derived without velocity evolution. For example, accounting self-consistently for the velocities can change the standard q=3.5 power-law index of the Dohnanyi (1969) differential size spectrum to an index as large as q=4. Similarly, for bodies held together by their own gravity, the corresponding power-law index range 2.88<q<3.14 of Pan & Sari (2005) can steepen to values as large as q=3.26. Our velocity results allow quantitative predictions of the bodies' scale heights as a function of size. Together with our predictions, observations of the scale heights for different sized bodies for the Kuiper belt, the asteroid belt, and extrasolar debris disks may constrain the mass and number of large bodies stirring the cascade as well as the colliding bodies' internal strengths.

"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 investigate potential biases in the measurements of exoplanet orbital parameters obtained from radial velocity observations for single-planet systems. We create a mock catalog of radial velocity data, choosing input planet masses, periods, and observing patterns from actual radial velocity surveys and varying input eccentricities. We apply Markov Chain Monte Carlo (MCMC) simulations and compare the resulting orbital parameters to the input values. We find that a combination of the effective signal-to-noise ratio of the data, the maximal gap in phase coverage, and the total number of periods covered by observations is a good predictor of the quality of derived orbit parameters. As eccentricity is positive definite, we find that eccentricities of planets on nearly circular orbits are preferentially overestimated, with typical bias of 1-2 times the median eccentricity uncertainty in a survey (e.g., 0.04 in the Butler et al. 2006 catalog). When performing population analysis, we recommend using the mode of the marginalized posterior eccentricity distribution to minimize potential biases. While the Butler et al. (2006) catalog reports eccentricities below 0.05 for just 17% of single-planet systems, we estimate that the true fraction of e < 0.05 orbits is about f(0.05)=38\pm 9%. For planets with P > 10 days, we find f(0.05)=28\pm 8% versus 10% from Butler et al. (2006). These planets either never acquired a large eccentricity or were circularized following any significant eccentricity excitation.

In order to attract and retain excellent researchers and diverse individuals in astrophysics, we recommend action be taken in several key areas impacting young scientists: (1) Maintain balance between large collaborations and individual projects through distribution of funding; encourage public releases of observational and simulation data for use by a broader community. (2) Improve the involvement of women, particularly at leading institutions. (3) Address the critical shortage of child care options and design reasonable profession-wide parental leave policies. (4) Streamline the job application and hiring process. We summarize our reasons for bringing these areas to the attention of the committee, and we suggest several practical steps that can be taken to address them.

The flow resulting from a strong ultrarelativistic shock moving through a stellar envelope with a polytrope-like density profile has been studied analytically and numerically at early times while the fluid temperature is relativistic--that is, just before and just after the shock breaks out of the star. Such a flow should expand and accelerate as its internal energy is converted to bulk kinetic energy; at late enough times, the assumption of relativistic temperatures becomes invalid. Here we present a new self-similar solution for the post-breakout flow when the accelerating fluid has bulk kinetic Lorentz factors much larger than unity but is cooling through $p/n$ of order unity to subrelativistic temperatures. This solution gives a relation between a fluid element's terminal Lorentz factor and that element's Lorentz factor just after it is shocked. Our numerical integrations agree well with the solution. While our solution assumes a planar flow, we show that corrections due to spherical geometry are important only for extremely fast ejecta originating in a region very close to the stellar surface. This region grows if the shock becomes relativistic deeper in the star.

We consider a strong ultrarelativistic shock moving through a star whose envelope has a polytrope-like density profile. When the shock is close to the star's outer boundary, its behavior follows the self-similar solution given by Sari (2005) for implosions in planar geometry. Here we outline this solution and find the asymptotic solution as the shock reaches the star's edge. We then show that the motion after the shock breaks out of the star is described by a self-similar solution remarkably like the solution for the motion inside the star. In particular, the characteristic Lorentz factor, pressure, and density vary with time according to the same power laws both before and after the shock breaks out of the star. After emergence from the star, however, the self-similar solution's characteristic position corresponds to a point behind the leading edge of the flow rather than at the shock front, and the relevant range of values for the similarity variable changes. Our numerical integrations agree well with the analytic results both before and after the shock reaches the star's edge.

The observed size distribution of Kuiper belt objects (KBOs)--small icy and rocky solar system bodies orbiting beyond Neptune--is well described by a power law at large KBO sizes. However, recent work by Bernstein et al. (2003) indicates that the size spectrum breaks and becomes shallower for KBOs smaller than about 70 km in size. Here we show that we expect such a break at KBO radius ~40 km since destructive collisions are frequent for smaller KBOs. Specifically, we assume that KBOs are rubble piles with low material strength rather than solid monoliths. This gives a power-law slope q~3 where the number N(r) of KBOs larger than a size r is given by $N(r) \propto r^{1-q}$; the break location follows from this slope through a self-consistent calculation. The existence of this break, the break's location, and the power-law slope we expect below the break are consistent with the findings of Bernstein et al. (2003). The agreement with observations indicates that KBOs are effectively strengthless rubble piles.

We develop a framework based on energy kicks for the evolution of high-eccentricity long-period orbits with Jacobi constant close to 3 in the restricted circular planar three-body problem where the secondary and primary masses have mass ratio mu << 1. We use this framework to explore mean-motion resonances between the test particle and the secondary mass. This approach leads to (i) a redefinition of resonance orders to reflect the importance of interactions at periapse; (ii) a pendulum-like equation describing the librations of resonance orbits; (iii) an analogy between these new fixed points and the Lagrangian points as well as between librations around the fixed points and the well known tadpole and horseshoe orbits; (iv) a condition a ~ mu^(-2/5) for the onset of chaos at large semimajor axis a; (v) the existence of a range mu < ~5 x 10^(-6) in secondary mass for which a test particle initially close to the secondary cannot escape from the system, at least in the planar problem; (vi) a simple explanation for the presence of asymmetric librations in exterior 1:N resonances.

We report spectroscopic and interferometric observations of the moderately metal-poor double-lined binary system HD 195987, with an orbital period of 57.3 days. By combining our radial-velocity and visibility measurements we determine the orbital elements and derive absolute masses for the components of M(A) = 0.844 +/- 0.018 Msun and M(B) = 0.6650 +/- 0.0079 Msun, with relative errors of 2% and 1%, respectively. We also determine the orbital parallax, pi(orb) = 46.08 +/- 0.27 mas, corresponding to a distance of 21.70 +/- 0.13 pc. The parallax and the measured brightness difference between the stars in V, H, and K yield the component absolute magnitudes in those bands. We also estimate the effective temperatures of the stars as Teff(A) = 5200 +/- 100 K and Teff(B) = 4200 +/- 200 K. Together with detailed chemical abundance analyses from the literature giving [Fe/H] approximately -0.5 (corrected for binarity) and [alpha/Fe] = +0.36, we use these physical properties to test current models of stellar evolution for metal-poor stars. Among the four that we considered, we find that no single model fits all observed properties at the measured composition, although we identify the assumptions in each one that account for the discrepancy and we conclude that a model with the proper combination of assumptions should be able to reproduce all the radiative properties. The indications from the isochrone fits and the pattern of enhancement of the metals in HD 195987 are consistent with this being a thick disk object, with an age of 10-12 Gyr.