A number of giant planet pairs discovered by the radial velocity method with period ratios $\lesssim 2 $ may reside in mean motion resonances. Convergent orbital migration and resonant capture at the time of formation would naturally explain the present-day resonant orbital configurations of these systems. Planets that experience smooth migration and eccentricity damping forces due to a proto-planetary disk should not only capture into mean motion resonances but also end up in a specific dynamical configuration within the resonance, sometimes referred to as apsidal corotation resonance (ACR). Here we develop a method for testing the hypothesis that a planet pair resides in an ACR by directly fitting radial velocity data. The ACR hypothesis strongly restricts the number of free parameters describing the radial velocity signal and we compare fits using this highly restricted model to fits using a more conventional two-planet RV model by using nested sampling simulations. We apply our method to HD 45364 and HD 33844, two systems hosting giant planet pairs in 3:2 and 5:3 resonances, respectively. We demonstrate that the observations of both systems support an ACR configuration and we use the results of our ACR model fits to constrain possible migration histories of these systems.

As the NASA Transiting Exoplanet Survey Satellite (TESS) fulfills its primary mission, it is executing an unprecedented survey of almost the entire sky: TESS's approved extended mission will likely extend sky coverage to ~94%, including ~60% of the ecliptic. In an accompanying note we demonstrated that `digital tracking' techniques can be used to efficiently `shift-and-stack' TESS full frame images (FFIs) and showed that combining ~1,300 exposures from a TESS sector gives a 50% detection threshold of $I_C\sim 22.0\pm0.5$, raising the possibility that TESS could discover the hypothesized Planet Nine. In this note, we estimate the yield for such a survey and demonstrate that this technique has the potential to discover hundreds of Kuiper Belt Objects, Scattered Disk Objects and Centaurs in TESS FFI data.

As the NASA Transiting Exoplanet Survey Satellite (TESS) fulfills its primary mission it is executing an unprecedented all-sky survey with the potential to discover distant planets in our own solar system, as well as hundreds of Transneptunian Objects (TNOs) and Centaurs. We demonstrate that shift-and-stack techniques can be used to efficiently search the Full-Frame Image (FFI) data from the TESS mission and survey the entire sky for outer Solar System objects down to $\sim22^{nd}$ magnitude.

The third flight of the High-Resolution Coronal Imager (Hi-C 2.1) occurred on May 29, 2018, the Sounding Rocket was launched from White Sands Missile Range in New Mexico. The instrument has been modified from its original configuration (Hi-C 1) to observe the solar corona in a passband that peaks near 172 Angstrom and uses a new, custom-built low-noise camera. The instrument targeted Active Region 12712, and captured 78 images at a cadence of 4.4 sec (18:56:22 - 19:01:57 UT; 5 min and 35 sec observing time). The image spatial resolution varies due to quasi-periodic motion blur from the rocket; sharp images contain resolved features of at least 0.47 arcsec. There are coordinated observations from multiple ground- and space-based telescopes providing an unprecedented opportunity to observe the mass and energy coupling between the chromosphere and the corona. Details of the instrument and the data set are presented in this paper.

We consider the potential for the Transiting Exoplanet Survey Satellite (TESS) to detect transit timing variations (TTVs) during both its nominal and extended mission phases. Building on previous estimates of the overall yield of planetary systems from the TESS mission, we predict that during its nominal two-year mission, TESS will observe measurable TTVs in $\sim30$ systems, from which $\mathcal{O}(10)$ planet will get precise mass measurements from TTVs alone, $\sim 5$ planets will have significant constraints placed on their masses from TTVs, and over a dozen systems will be singly-transiting TTV systems. We consider a number of different extended mission scenarios, and predict that in a typical scenario, an extended mission will allow TESS to increase the number of systems with measurable TTVs to a total of $\sim90$, from which $\sim 15$ planets will have precise mass measurements, another $\sim 15$ will have significant constraints placed on their masses, and $\sim60$ will be singly-transiting TTV systems. We also describe how follow-up transit observations of multi-planet systems discovered by the TESS mission can be optimally planned to maximize TTV mass and eccentricity constraints.

We describe the digest2 software package, a fast, short-arc orbit classifier for small Solar System bodies. The digest2 algorithm has been serving the community for more than 13 years. The code provides a score, D2, which represents a pseudo-probability that a tracklet belongs to a given Solar System orbit type. digest2 is primarily used as a classifier for Near-Earth Object (NEO) candidates, to identify those to be prioritized for follow-up observation. We describe the historical development of digest2 and demonstrate its use on real and synthetic data. We find that digest2 can accurately and precisely distinguish NEOs from non- NEOs. At the time of detection, 14% of NEO tracklets and 98.5% of non-NEOs tracklets have D2 below the critical value of D2 = 65. 94% of our simulated NEOs achieved the maximum D2 = 100 and 99.6% of NEOs achieved $D2 \ge 65$ at least once during the simulated 10-year timeframe. We demonstrate that D2 varies as a function of time, rate of motion, magnitude and sky-plane location, and show that NEOs tend to have lower D2 at low Solar elongations close to the ecliptic. We use our findings to recommend future development directions for the digest2 code.

We present results of the largest, most comprehensive study ever done of the stellar multiplicity of the most common stars in the Galaxy, the red dwarfs. We have conducted an all-sky, volume-limited survey for stellar companions to 1120 M dwarf primaries known to lie within 25 pc of the Sun via trigonometric parallaxes. In addition to a comprehensive literature search, stars were explored in new surveys for companions at separations of 2" to 300". A reconnaissance of wide companions to separations of 300" was done via blinking archival images. I-band images were used to search our sample for companions at separations of 2" to 180". Various astrometric and photometric methods were used to probe the inner 2" to reveal close companions. We report the discovery of 20 new companions and identify 56 candidate multiple systems. We find a stellar multiplicity rate of 26.8 +/- 1.4% and a stellar companion rate of 32.4 +/- 1.4% for M dwarfs. There is a broad peak in the separation distribution of the companions at 4 -- 20 AU, with a weak trend of smaller projected linear separations for lower mass primaries. A hint that M dwarf multiplicity may be a function of tangential velocity is found, with faster moving, presumably older, stars found to be multiple somewhat less often. We calculate that stellar companions make up at least 17% of mass attributed to M dwarfs in the solar neighborhood, with roughly 11% of M dwarf mass hidden as unresolved companions. Finally, when considering all M dwarf primaries and companions, we find that the mass distribution for M dwarfs increases to the end of the stellar main sequence.

The existence of Planet Nine has been suggested to explain the pericenter clustering of extreme trans-Neptunian objects (TNOs). However, the underlying dynamics involving Planet Nine, test particles and Neptune is rich, and it remains unclear which dynamical processes lead to the alignment and how they depends on the properties of Planet Nine. Here, we investigate the secular interactions between an eccentric outer perturber and TNOs starting in a near-coplanar configuration. We find that a large number of TNOs could survive outside of mean motion resonances at 4Gyr, which differs from previous results obtained in the exact coplanar case with Neptune being treated as a quadrupole potential. In addition, secular dynamics leads to the orbital clustering seen in N-body simulations. We find that a near coplanar Planet Nine can flip TNO orbital planes, and when this happens, the geometrical longitudes of pericenter of the TNOs librate around $180^\circ$ during the flip. Orbital precession caused by the inner giant planets can suppress the flips while keeping the longitude of pericenter librating when $30\lesssim r_p \lesssim 80$ AU $\&$ $a\gtrsim 250$ AU. This results in the alignment of the pericenter of the low inclination TNOs ($i \lesssim 40^\circ$). We find the anti-aligned population and the flipped orbits could be produced by an eccentric ($e_9 \gtrsim 0.4$) outer planet of $\sim 10M_{\oplus}$ in a wide $a_9 \gtrsim 400 \sim 800$ AU orbit. Future surveys on the high inclination TNOs will help further constrain the properties of possible outer planets.

In this paper we investigate systems previously identified to exhibit transit timing variations (TTVs) in Kepler data, with the goal of predicting the expected improvements to the mass and eccentricity constraints that will arise from combining Kepler data with future data from the TESS mission. We advocate for the use of the Kullback-Leibler (KL) divergence as a means to quantify improvements in the measured constraints. Compared to the original Kepler data, the TESS data will have a lower signal-to-noise ratio, rendering some of the planetary transits undetectable, and lowering the accuracy with which the transit mid-time can be estimated. Despite these difficulties, out of the 55 systems (containing 143 planets) investigated, we predict that the collection of short-cadence data by TESS will be of significant value (i.e. it will improve the mass uncertainty such that the KL divergence is > 0.1) for approximately 6 - 14 planets during the nominal mission, with the range primarily driven by the uncertain precision with which transit mid-times will be recovered from TESS data. In an extended mission this would increase to a total of approximately 12 - 25 planets.

We studied the Near-Earth Asteroid (NEA) candidates posted on the Minor Planet Center's Near-Earth Object Confirmation Page (NEOCP) between years 2013 and 2016. Out of more than 17,000 NEA candidates, while the majority became either new discoveries or were associated with previously known objects, about 11% were unable to be followed-up or confirmed. We further demonstrate that of the unconfirmed candidates, 926+/-50 are likely to be NEAs, representing 18% of discovered NEAs in that period. Only 11% (~93) of the unconfirmed NEA candidates were large (having absolute magnitude H<22). To identify the reasons why these NEAs were not recovered, we analyzed those from the most prolific asteroid surveys: Pan-STARRS, the Catalina Sky Survey, the Dark Energy Survey, and the Space Surveillance Telescope. We examined the influence of plane-of-sky positions and rates of motion, brightnesses, submission delays, and computed absolute magnitudes, as well as correlations with the phase of the moon and seasonal effects. We find that delayed submission of newly discovered NEA candidate to the NEOCP drove a large fraction of the unconfirmed NEA candidates. A high rate of motion was another significant contributing factor. We suggest that prompt submission of suspected NEA discoveries and rapid response to fast moving targets and targets with fast growing ephemeris uncertainty would allow better coordination among dedicated follow-up observers, decrease the number of unconfirmed NEA candidates, and increase the discovery rate of NEAs.

We present a novel approach to the minor planet linking problem. Our heliocentric transformation-and-propagation algorithm clusters tracklets at common epochs, allowing for the efficient identification of tracklets that represent the same minor planet. This algorithm scales as $O(N log N)$, with the number of tracklets $N$, a significant advance over standard methods, which scale as $O(N^3)$. This overcomes one of the primary computational bottlenecks faced by current and future asteroid surveys. We apply our algorithm to the Minor Planet Center's Isolated Tracklet File, establishing orbits for $\sim41,000$ new minor planets.

Observations of clustering among the orbits of the most distant trans-Neptunian objects (TNOs) has inspired interest in the possibility of an undiscovered ninth planet lurking in the outskirts of the solar system. Numerical simulations by a number of authors have demonstrated that, with appropriate choices of planet mass and orbit, such a planet can maintain clustering in the orbital elements of the population of distant TNOs, similar to the observed sample. However, many aspects of the rich underlying dynamical processes induced by such a distant eccentric perturber have not been fully explored. We report the results of our investigation of the dynamics of coplanar test-particles that interact with a massive body on an circular orbit (Neptune) and a massive body on a more distant, highly eccentric orbit (the putative Planet Nine). We find that a detailed examination of our idealized simulations affords tremendous insight into the rich test-particle dynamics that are possible. In particular, we find that chaos and resonance overlap plays an important role in particles' dynamical evolution. We develop a simple mapping model that allows us to understand, in detail, the web of overlapped mean-motion resonances explored by chaotically evolving particles. We also demonstrate that gravitational interactions with Neptune can have profound effects on the orbital evolution of particles. Our results serve as a starting point for a better understanding of the dynamical behavior observed in more complicated simulations that can be used to constrain the mass and orbit of Planet 9.

We report the discovery of a $H_r = 3.4\pm0.1$ dwarf planet candidate by the Pan-STARRS Outer Solar System Survey. 2010 JO$_{179}$ is red with $(g-r)=0.88 \pm 0.21$, roughly round, and slowly rotating, with a period of $30.6$ hr. Estimates of its albedo imply a diameter of 600--900~km. Observations sampling the span between 2005--2016 provide an exceptionally well-determined orbit for 2010 JO$_{179}$, with a semi-major axis of $78.307\pm0.009$ au, distant orbits known to this precision are rare. We find that 2010 JO$_{179}$ librates securely within the 21:5 mean-motion resonance with Neptune on hundred-megayear time scales, joining the small but growing set of known distant dwarf planets on metastable resonant orbits. These imply a substantial trans-Neptunian population that shifts between stability in high-order resonances, the detached population, and the eroding population of the scattering disk.

Roughly 1000 white dwarfs are known to be polluted with planetary material, and the progenitors of this material are typically assumed to be asteroids. The dynamical architectures which perturb asteroids into white dwarfs are still unknown, but may be crucially dependent on moons liberated from parent planets during post-main-sequence gravitational scattering. Here, we trace the fate of these exomoons, and show that they more easily achieve deep radial incursions towards the white dwarf than do scattered planets. Consequently, moons are likely to play a significant role in white dwarf pollution, and in some cases may be the progenitors of the pollution itself.

In this work we report the detection of seven Neptune Trojans (NTs) in the Pan-STARRS 1 (PS1) survey. Five of these are new discoveries, consisting of four L4 Trojans and one L5 Trojan. Our orbital simulations show that the L5 Trojan stably librates for only several million years. This suggests that the L5 Trojan must be of recent capture origin. On the other hand, all four new L4 Trojans stably occupy the 1:1 resonance with Neptune for more than 1 Gyr. They can, therefore, be of primordial origin. Our survey simulation results show that the inclination width of the Neptune Trojan population should be between $7^{\circ}$ and $27^{\circ}$ at $>$ 95% confidence, and most likely $\sim 11^{\circ}$. In this paper, we describe the PS1 survey, the Outer Solar System pipeline, the confirming observations, and the orbital/physical properties of the new Neptune Trojans.

Although the majority of Centaurs are thought to have originated in the scattered disk, with the high-inclination members coming from the Oort cloud, the origin of the high inclination component of trans-Neptunian objects (TNOs) remains uncertain. We report the discovery of a retrograde TNO, which we nickname "Niku", detected by the Pan-STARRS 1 Outer Solar System Survey. Our numerical integrations show that the orbital dynamics of Niku are very similar to that of 2008 KV$_{42}$ (Drac), with a half-life of $\sim 500$ Myr. Comparing similar high inclination TNOs and Centaurs ($q > 10$ AU, $a < 100$ AU and $i > 60^\circ$), we find that these objects exhibit a surprising clustering of ascending node, and occupy a common orbital plane. This orbital configuration has high statistical significance: 3.8-$\sigma$. An unknown mechanism is required to explain the observed clustering. This discovery may provide a pathway to investigate a possible reservoir of high-inclination objects.

We present ALMA band 7 (345 GHz) continuum and $^{12}$CO(J = 3-2) observations of the circumstellar disk surrounding HD141569. At an age of about 5 Myr, the disk has a complex morphology that may be best interpreted as a nascent debris system with gas. Our $870\rm~\mu m$ ALMA continuum observations resolve a dust disk out to approximately $ 56 ~\rm au$ from the star (assuming a distance of 116 pc) with $0."38$ resolution and $0.07 ~ \rm mJy~beam^{-1}$ sensitivity. We measure a continuum flux density for this inner material of $3.8 \pm 0.4 ~ \rm mJy$ (including calibration uncertainties). The $^{12}$CO(3-2) gas is resolved kinematically and spatially from about 30 to 210 au. The integrated $^{12}$CO(3-2) line flux density is $15.7 \pm 1.6~\rm Jy~km~s^{-1}$. We estimate the mass of the millimeter debris and $^{12}$CO(3-2) gas to be $\gtrsim0.04~\rm M_{\oplus}$ and $\sim2\times 10^{-3}~\rm M_{\oplus}$, respectively. If the millimeter grains are part of a collisional cascade, then we infer that the inner disk ($<50$ au) has $\sim 160~\rm M_{\oplus}$ contained within objects less than 50 km in radius, depending on the planetesimal size distribution and density assumptions. MCMC modeling of the system reveals a disk morphology with an inclination of $53.4^{\circ}$ centered around a $\rm M=2.39~ M_{\odot}$ host star ($\rm Msin(i)=1.92~ M_{\odot}$). We discuss whether the gas in HD141569's disk may be second generation. If it is, the system can be used to study the clearing stages of planet formation.

We significantly constrain the sky position, distance, and mass of a possible additional, distant planet in the solar system by examining its influence on the distance between Earth and the Cassini Spacecraft. Our preferred region is approximately centered on (RA, Dec) = ($40\arcdeg$, $-15\arcdeg$), extending approximately 20 degrees in all directions.

Previous studies indicate that more than a quarter of all white dwarf (WD) atmospheres are polluted by remnant planetary material, with some WDs being observed to accrete the mass of Pluto in 10^6 years. The short sinking timescale for the pollutants indicate that the material must be frequently replenished. Moons may contribute decisively to this pollution process if they are liberated from their parent planets during the post-main-sequence evolution of the planetary systems. Here, we demonstrate that gravitational scattering events among planets in WD systems easily triggers moon ejection. Repeated close encounters within tenths of a planetary Hill radii are highly destructive to even the most massive, close-in moons. Consequently, scattering increases both the frequency of perturbing agents in WD systems, as well as the available mass of polluting material in those systems, thereby enhancing opportunities for collision and fragmentation and providing more dynamical pathways for smaller bodies to reach the WD. Moreover, during intense scattering, planets themselves have pericenters with respect to the WD of only a fraction of an AU, causing extreme Hill-sphere contraction, and the liberation of moons into WD-grazing orbits. Many of our results are directly applicable to exomoons orbiting planets around main sequence stars.

We use astrometry of Pluto and other TNOs to constrain the sky location, distance, and mass of the possible additional planet (Planet Nine) hypothesized by Batygin and Brown (2016). We find that over broad regions of the sky, the inclusion of a massive, distant planet degrades the fits to the observations. However, in other regions, the fits are significantly improved by the addition of such a planet. Our best fits suggest a planet that is either more massive or closer than argued for by Batygin and Brown (2016) based on the orbital distribution of distant trans-neptunian objects (or by Fienga et al. (2016) based on range measured to the Cassini spacecraft). The trend to favor larger and closer perturbing planets is driven by the residuals to the astrometry of Pluto, remeasured from photographic plates using modern stellar catalogs (Buie and Folkner 2015), which show a clear trend in declination, over the course of two decades, that drive a preference for large perturbations. Although this trend may be the result of systematic errors of unknown origin in the observations, a possible resolution is that the declination trend may be due to perturbations from a body, additional to Planet Nine , that is closer to Pluto, but less massive than, Planet Nine .

We report constraints on the three-dimensional orbital architecture for all four planets known to orbit the nearby M dwarf Gliese 876 based solely on Doppler measurements and demanding long-term orbital stability. Our dataset incorporates publicly available radial velocities taken with the ELODIE and CORALIE spectrographs, HARPS, and Keck HIRES as well as previously unpublished HIRES velocities. We first quantitatively assess the validity of the planets thought to orbit GJ 876 by computing the Bayes factors for a variety of different coplanar models using an importance sampling algorithm. We find that a four-planet model is preferred over a three-planet model. Next, we apply a Newtonian MCMC algorithm to perform a Bayesian analysis of the planet masses and orbits using an n-body model in three-dimensional space. Based on the radial velocities alone, we find that a 99% credible interval provides upper limits on the mutual inclinations for the three resonant planets ($\Phi_{cb}<6.20^\circ$ for the "c" and "b" pair and $\Phi_{be}<28.5^\circ$ for the "b" and "e" pair). Subsequent dynamical integrations of our posterior sample find that the GJ 876 planets must be roughly coplanar ($\Phi_{cb}<2.60^\circ$ and $\Phi_{be}<7.87^\circ$), suggesting the amount of planet-planet scattering in the system has been low. We investigate the distribution of the respective resonant arguments of each planet pair and find that at least one argument for each planet pair and the Laplace argument librate. The libration amplitudes in our three-dimensional orbital model supports the idea of the outer-three planets having undergone significant past disk migration.

We report the discovery of Kepler-432b, a giant planet ($M_b = 5.41^{+0.32}_{-0.18} M_{\rm Jup}, R_b = 1.145^{+0.036}_{-0.039} R_{\rm Jup}$) transiting an evolved star $(M_\star = 1.32^{+0.10}_{-0.07} M_\odot, R_\star = 4.06^{+0.12}_{-0.08} R_\odot)$ with an orbital period of $P_b = 52.501129^{+0.000067}_{-0.000053}$ days. Radial velocities (RVs) reveal that Kepler-432b orbits its parent star with an eccentricity of $e = 0.5134^{+0.0098}_{-0.0089}$, which we also measure independently with asterodensity profiling (AP; $e=0.507^{+0.039}_{-0.114}$), thereby confirming the validity of AP on this particular evolved star. The well-determined planetary properties and unusually large mass also make this planet an important benchmark for theoretical models of super-Jupiter formation. Long-term RV monitoring detected the presence of a non-transiting outer planet (Kepler-432c; $M_c \sin{i_c} = 2.43^{+0.22}_{-0.24} M_{\rm Jup}, P_c = 406.2^{+3.9}_{-2.5}$ days), and adaptive optics imaging revealed a nearby (0\farcs87), faint companion (Kepler-432B) that is a physically bound M dwarf. The host star exhibits high signal-to-noise asteroseismic oscillations, which enable precise measurements of the stellar mass, radius and age. Analysis of the rotational splitting of the oscillation modes additionally reveals the stellar spin axis to be nearly edge-on, which suggests that the stellar spin is likely well-aligned with the orbit of the transiting planet. Despite its long period, the obliquity of the 52.5-day orbit may have been shaped by star-planet interaction in a manner similar to hot Jupiter systems, and we present observational and theoretical evidence to support this scenario. Finally, as a short-period outlier among giant planets orbiting giant stars, study of Kepler-432b may help explain the distribution of massive planets orbiting giant stars interior to 1 AU.

Initially designed to discover short-period planets, the N2K campaign has since evolved to discover new worlds at large separations from their host stars. Detecting such worlds will help determine the giant planet occurrence at semi-major axes beyond the ice line, where gas giants are thought to mostly form. Here we report four newly-discovered gas giant planets (with minimum masses ranging from 0.4 to 2.1 MJup) orbiting stars monitored as part of the N2K program. Two of these planets orbit stars already known to host planets: HD 5319 and HD 11506. The remaining discoveries reside in previously-unknown planetary systems: HD 10442 and HD 75784. The refined orbital period of the inner planet orbiting HD 5319 is 641 days. The newly-discovered outer planet orbits in 886 days. The large masses combined with the proximity to a 4:3 mean motion resonance make this system a challenge to explain with current formation and migration theories. HD 11506 has one confirmed planet, and here we confirm a second. The outer planet has an orbital period of 1627.5 days, and the newly-discovered inner planet orbits in 223.6 days. A planet has also been discovered orbiting HD 75784 with an orbital period of 341.7 days. There is evidence for a longer period signal; however, several more years of observations are needed to put tight constraints on the Keplerian parameters for the outer planet. Lastly, an additional planet has been detected orbiting HD 10442 with a period of 1043 days.

Gas giant planets in the Solar system host large satellite systems with multiple regular and irregular moons. Regular moons revolve around their host planet in circular, low inclination short period orbits, and are thought to form in-situ through coagulation processes. In contrast, irregular moons have highly inclined (and even retrograde), typically eccentric and long period orbits around their host planet. Irregular moons are therefore often thought to be unbound objects in helio-centric orbits that were later captured to their current orbits around the planet. Here we study the a different possibility in which regular moons form in-situ, and dynamically evolve through mutual moon-moon scattering. We find that such evolution can excite the satellites into high eccentricities and inclinations. We find that moons are either ejected from the host planet to become runaway moons, or stay bound and become prograde orbiting irregular moons with inclined and eccentric orbits around their host planet. Ejected moons, unbound to the planet, can later be temporarily re-captured by the host planet even at retrograde orbits. Such moons are eventually re-ejected from the system or collide with the planet, at least in the absence of dissipative processes (e.g. collisions with existing bound moons, a debris disk or through tidal interactions with the host planet), not currently modeled. Uncaptured runaway moons may eventually be ejected from the Solar system, or be captured into stable helio-centric orbits and contribute to the populations of asteroidal or trans-Neptunian objects. Such scenarios are potentially relevant both for the gas-giant satellites in the Solar system and for the dynamical evolution of exomoons.

We present an updated study of the planets known to orbit 55 Cancri A using 1,418 high-precision radial velocity observations from four observatories (Lick, Keck, Hobby-Eberly Telescope, Harlan J. Smith Telescope) and transit time/durations for the inner-most planet, 55 Cancri "e" (Winn et al. 2011). We provide the first posterior sample for the masses and orbital parameters based on self-consistent n-body orbital solutions for the 55 Cancri planets, all of which are dynamically stable (for at least $10^8$ years). We apply a GPU version of Radial velocity Using N-body Differential evolution Markov Chain Monte Carlo (RUN DMC; B. Nelson et al. 2014) to perform a Bayesian analysis of the radial velocity and transit observations. Each of the planets in this remarkable system has unique characteristics. Our investigation of high-cadence radial velocities and priors based on space-based photometry yields an updated mass estimate for planet "e" ($8.09\pm0.26$ M$_\oplus$), which affects its density ($5.51\pm^{1.32}_{1.00}$ g cm$^{-3}$) and inferred bulk composition. Dynamical stability dictates that the orbital plane of planet "e" must be aligned to within $60^o$ of the orbital plane of the outer planets (which we assume to be coplanar). The mutual interactions between the planets "b" and "c" may develop an apsidal lock about $180^o$. We find 36-45% of all our model systems librate about the anti-aligned configuration with an amplitude of $51^o\pm^{6^o}_{10^o}$. Other cases showed short-term perturbations in the libration of $\varpi_b-\varpi_c$, circulation, and nodding, but we find the planets are not in a 3:1 mean-motion resonance. A revised orbital period and eccentricity for planet "d" pushes it further toward the closest known Jupiter analog in the exoplanet population.

In the 20+ years of Doppler observations of stars, scientists have uncovered a diverse population of extrasolar multi-planet systems. A common technique for characterizing the orbital elements of these planets is Markov chain Monte Carlo (MCMC), using a Keplerian model with random walk proposals and paired with the Metropolis-Hastings algorithm. For approximately a couple of dozen planetary systems with Doppler observations, there are strong planet-planet interactions due to the system being in or near a mean-motion resonance (MMR). An N-body model is often required to accurately describe these systems. Further computational difficulties arise from exploring a high-dimensional parameter space ($\sim$7 x number of planets) that can have complex parameter correlations. To surmount these challenges, we introduce a differential evolution MCMC (DEMCMC) applied to radial velocity data while incorporating self-consistent N-body integrations. Our Radial velocity Using N-body DEMCMC (RUN DMC) algorithm improves upon the random walk proposal distribution of the traditional MCMC by using an ensemble of Markov chains to adaptively improve the proposal distribution. We describe the methodology behind the algorithm, along with results of tests for accuracy and performance. We find that most algorithm parameters have a modest effect on the rate of convergence. However, the size of the ensemble can have a strong effect on performance. We show that the optimal choice depends on the number of planets in a system, as well as the computer architecture used and the resulting extent of parallelization. While the exact choices of optimal algorithm parameters will inevitably vary due to the details of individual planetary systems, we offer recommendations for choosing the DEMCMC algorithm's algorithmic parameters that result in excellent performance for a wide variety of planetary systems.

We have analyzed new and previously published radial velocity observations of MARVELS-1, known to have an ostensibly substellar companion in a ~6- day orbit. We find significant (~100 m/s) residuals to the best-fit model for the companion, and these residuals are naively consistent with an interior giant planet with a P = 1.965d in a nearly perfect 3:1 period commensuribility (|Pb/Pc - 3| < 10^{-4}). We have performed several tests for the reality of such a companion, including a dynamical analysis, a search for photometric variability, and a hunt for contaminating stellar spectra. We find many reasons to be critical of a planetary interpretation, including the fact that most of the three-body dynamical solutions are unstable. We find no evidence for transits, and no evidence of stellar photometric variability. We have discovered two apparent companions to MARVELS-1 with adaptive optics imaging at Keck; both are M dwarfs, one is likely bound, and the other is likely a foreground object. We explore false-alarm scenarios inspired by various curiosities in the data. Ultimately, a line profile and bisector analysis lead us to conclude that the ~100 m/s residuals are an artifact of spectral contamination from a stellar companion contributing ~15-30% of the optical light in the system. We conclude that origin of this contamination is the previously detected radial velocity companion to MARVELS-1, which is not, as previously reported, a brown dwarf, but in fact a G dwarf in a face-on orbit.

We present an updated analysis of radial velocity data of the HD 82943 planetary system based on 10 years of measurements obtained with the Keck telescope. Previous studies have shown that the HD 82943 system has two planets that are likely in 2:1 mean-motion resonance (MMR), with the orbital periods about 220 and 440 days (Lee et al. 2006). However, alternative fits that are qualitatively different have also been suggested, with two planets in a 1:1 resonance (Gozdziewski & Konacki 2006) or three planets in a Laplace 4:2:1 resonance (Beauge et al. 2008). Here we use \c{hi}2 minimization combined with parameter grid search to investigate the orbital parameters and dynamical states of the qualitatively different types of fits, and we compare the results to those obtained with the differential evolution Markov chain Monte Carlo method. Our results support the coplanar 2:1 MMR configuration for the HD 82943 system, and show no evidence for either the 1:1 or 3-planet Laplace resonance fits. The inclination of the system with respect to the sky plane is well constrained at about 20(+4.9 -5.5) degree, and the system contains two planets with masses of about 4.78 MJ and 4.80 MJ (where MJ is the mass of Jupiter) and orbital periods of about 219 and 442 days for the inner and outer planet, respectively. The best fit is dynamically stable with both eccentricity-type resonant angles {\theta}1 and {\theta}2 librating around 0 degree.

We propose a scientific program to complete a census of planets, characterizing their masses, orbital properties, and dynamical histories using continued observations of the Kepler field of view with the Kepler spacecraft in a two reaction wheel mission. Even with a significantly reduced photometric precision, extending time-domain observations of this field is uniquely capable of pursuing several critical science goals: 1) measuring the architectures of planetary systems by identifying non-transiting planets interleaved among known transiting planets, 2) establishing the mass-radius relationship for planets in the important transition region between small, gas-rich sub-Neptune planets and large, rocky super-Earths, and 3) uncovering dynamical evidence of the formation and evolution of the inner regions of planetary systems. To meet these objectives, the unique multi-object observing capabilities of Kepler will be used in a set of concurrent campaigns with specific motivations. These campaigns focus largely on the ability to interpret Transit Timing Variations (TTVs) that result from dynamical interactions among planets in a system and include: 1) observations of systems that exhibit large TTVs and are particularly rich in dynamical information, 2) observations of systems where additional transit times will yield mass measurements of the constituent planets, 3) observations of systems where the TTV signal evolves over very long timescales, and 4) observations of systems with long-period planet candidates where additional transits will remove orbital period ambiguities caused by gaps in the original Kepler data.

We propose to use Kepler in 2-wheel mode to conduct a detailed search for Earth-sized planets orbiting ultra-cool stars and brown dwarfs (spectral types from M7 to L3). This population of objects presents several advantages for exoplanet surveys. First, ultra-cool stars and brown dwarfs are small and thus result in favorable planet-to-star area ratios. Second, because of their low effective temperature, the inner edge of their habitable zone is extremely close (2 to 3 days only). Third, our targets are bright at infrared wavelengths, which will enable detailed follow-up studies. Our program therefore represents a unique opportunity to find a transiting Earth-size exoplanet for which atmospheric features (including biosignatures) could be detected with near-to-come facilities such as JWST. Such exoplanet has not been discovered yet. Kepler in 2-wheel mode provides the required stability and photometric precision to make this survey successful. Our initial target sample includes 60 ultra-cool stars and brown dwarfs from which we expect to detect at least one transiting planet. We propose to monitor each source for 4 days, resulting in a total program duration of ~240 days.

We perform numerical integrations of four-body (star, planet, planet, satellite) systems to investigate the stability of satellites in planetary Systems with Tightly-packed Inner Planets (STIPs). We find that the majority of closely-spaced stable two-planet systems can stably support satellites across a range of parameter-space which is only slightly decreased compared to that seen for the single-planet case. In particular, circular prograde satellites remain stable out to $\sim 0.4 R_H$ (where $R_H$ is the Hill Radius) as opposed to $\sim 0.5 R_H$ in the single-planet case. A similarly small restriction in the stable parameter-space for retrograde satellites is observed, where planetary close approaches in the range 2.5 to 4.5 mutual Hill radii destabilize most satellites orbits only if $a\sim 0.65 R_H$. In very close planetary pairs (e.g. the 12:11 resonance) the addition of a satellite frequently destabilizes the entire system, causing extreme close-approaches and the loss of satellites over a range of circumplanetary semi-major axes. The majority of systems investigated stably harbored satellites over a wide parameter-space, suggesting that STIPs can generally offer a dynamically stable home for satellites, albeit with a slightly smaller stable parameter-space than the single-planet case. As we demonstrate that multi-planet systems are not a priori poor candidates for hosting satellites, future measurements of satellite occurrence rates in multi-planet systems versus single-planet systems could be used to constrain either satellite formation or past periods of strong dynamical interaction between planets.

Motivated by the population of multi-planet systems with orbital period ratios 1<P2/P1<2, we study the long-term stability of packed two planet systems. The Hamiltonian for two massive planets on nearly circular and nearly coplanar orbits near a first order mean motion resonance can be reduced to a one degree of freedom problem (Sessin & Ferraz Mello (1984), Wisdom (1986), Henrard et al. (1986)). Using this analytically tractable Hamiltonian, we apply the resonance overlap criterion to predict the onset of large scale chaotic motion in close two planet systems. The reduced Hamiltonian has only a weak dependence on the planetary mass ratio, and hence the overlap criterion is independent of the planetary mass ratio at lowest order. Numerical integrations confirm that the planetary mass ratio has little effect on the structure of the chaotic phase space for close orbits in the low eccentricity (e <~0.1) regime. We show numerically that orbits in the chaotic web produced primarily by first order resonance overlap eventually experience large scale erratic variation in semimajor axes and are Lagrange unstable. This is also true of the orbits in this overlap region which are Hill stable. As a result, we can use the first order resonance overlap criterion as an effective stability criterion for pairs of observed planets. We show that for low mass (<~10 M_Earth) planetary systems with initially circular orbits the period ratio at which complete overlap occurs and widespread chaos results lies in a region of parameter space which is Hill stable. Our work indicates that a resonance overlap criterion which would apply for initially eccentric orbits needs to take into account second order resonances. Finally, we address the connection found in previous work between the Hill stability criterion and numerically determined Lagrange instability boundaries in the context of resonance overlap.

We report the radial-velocity discovery of a second planetary mass companion to the K0 V star HD 37605, which was already known to host an eccentric, P~55 days Jovian planet, HD 37605b. This second planet, HD 37605c, has a period of ~7.5 years with a low eccentricity and an Msini of ~3.4 MJup. Our discovery was made with the nearly 8 years of radial velocity follow-up at the Hobby-Eberly Telescope and Keck Observatory, including observations made as part of the Transit Ephemeris Refinement and Monitoring Survey (TERMS) effort to provide precise ephemerides to long-period planets for transit follow-up. With a total of 137 radial velocity observations covering almost eight years, we provide a good orbital solution of the HD 37605 system, and a precise transit ephemeris for HD 37605b. Our dynamic analysis reveals very minimal planet-planet interaction and an insignificant transit time variation. Using the predicted ephemeris, we performed a transit search for HD 37605b with the photometric data taken by the T12 0.8-m Automatic Photoelectric Telescope (APT) and the Microvariability and Oscillations of Stars (MOST) satellite. Though the APT photometry did not capture the transit window, it characterized the stellar activity of HD 37605, which is consistent of it being an old, inactive star, with a tentative rotation period of 57.67 days. The MOST photometry enabled us to report a dispositive null detection of a non-grazing transit for this planet. Within the predicted transit window, we exclude an edge-on predicted depth of 1.9% at >>10sigma, and exclude any transit with an impact parameter b>0.951 at greater than 5sigma. We present the BOOTTRAN package for calculating Keplerian orbital parameter uncertainties via bootstrapping. We found consistency between our orbital parameters calculated by the RVLIN package and error bars by BOOTTRAN with those produced by a Bayesian analysis using MCMC.

At least two multi-planetary systems in a 4:3 mean motion resonance have been found by radial velocity surveys. These planets are gas giants and the systems are only stable when protected by a resonance. Additionally the Kepler mission has detected at least 4 strong candidate planetary systems with a period ratio close to 4:3. This paper investigates traditional dynamical scenarios for the formation of these systems. We systematically study migration scenarios with both N-body and hydro-dynamic simulations. We investigate scenarios involving the in-situ formation of two planets in resonance. We look at the results from finely tuned planet-planet scattering simulations with gas disk damping. Finally, we investigate a formation scenario involving isolation-mass embryos. Although the combined planet-planet scattering and damping scenario seems promising, none of the above scenarios is successful in forming enough systems in 4:3 resonance with planetary masses similar to the observed ones. This is a negative result but it has important implications for planet formation. Previous studies were successful in forming 2:1 and 3:2 resonances. This is generally believed to be evidence of planet migration. We highlight the main differences between those studies and our failure in forming a 4:3 resonance. We also speculate on more exotic and complicated ideas. These results will guide future investigators toward exploring the above scenarios and alternative mechanisms in a more general framework.

The chance that a planetary system will interact with another member of its host star's nascent cluster would be greatly increased if gas giant planets form in situ on wide orbits. In this paper, we explore the outcomes of planet-planet scattering for a distribution of multiplanet systems that all have one of the planets on an initial orbit of 100 AU. The scattering experiments are run with and without stellar flybys. We convolve the outcomes with distributions for protoplanetary disk and stellar cluster sizes to generalize the results where possible. We find that the frequencies of large mutual inclinations and high eccentricities are sensitive to the number of planets in a system, but not strongly to stellar flybys. However, flybys do play a role in changing the low and moderate portions of the mutual inclination distributions, and erase dynamically cold initial conditions on average. Wide-orbit planets can be mixed throughout the planetary system, and in some cases, can potentially become hot Jupiters, which we demonstrate using scattering experiments that include a tidal damping model. If planets form on wide orbits in situ, then there will be discernible differences in the proper motion distributions of a sample of wide-orbit planets compared with a pure scattering formation mechanism. Stellar flybys can enhance the frequency of ejections in planetary systems, but auto-ionization is likely to remain the dominant source of free-floating planets.

The dynamical evolution of planetary systems leaves observable signatures in debris disks. Optical images trace micron-sized grains, which are strongly affected by stellar radiation and need not coincide with their parent body population. Observations of mm-size grains accurately trace parent bodies, but previous images lack the resolution and sensitivity needed to characterize the ring's morphology. Here we present ALMA 350 GHz observations of the Fomalhaut debris ring. These observations demonstrate that the parent body population is 13-19 AU wide with a sharp inner and outer boundary. We discuss three possible origins for the ring, and suggest that debris confined by shepherd planets is the most consistent with the ring's morphology.

Most detected planet-bearing binaries are in wide orbits, for which a high inclination, $i_B$, between the binary orbital plane and the plane of the planetary disk around the primary is likely to be common. In this paper, we investigate the intermediate stages - from planetesimals to planetary embryos/cores - of planet formation in such highly inclined cases. Our focus is on the effects of gas drag on the planetesimals' orbital evolution, in particular on the evolution of the planetesimals' semimajor axis distribution and their mutual relative velocities. We first demonstrate that a non-evolving axisymmetric disk model is a good approximation for studying the effects of gas drag on a planetesimal in the highly inclined case ($30^\circ <i_B<150^\circ$). We then find that gas drag plays a crucial role, and the results can be generally divided into two categories, i.e., the Kozai-on regime and the Kozai-off regime, depending on the specific value of $i_B$. For both regimes, a robust outcome over a wide range of parameters is that, planetesimals migrate/jump inwards and pile up, leading to a severely truncated and dense planetesimal disk around the primary. In this compact and dense disk, collision rates are high but relative velocities are low, providing conditions which are favorable for planetesimal growth, and potentially allow for the subsequent formation of planets.

Transit surveys combined with Doppler data have revealed a class of gas giant planets that are massive and highly enriched in heavy elements (e.g., HD149026b, GJ436b, and HAT-P-20b). It is tempting to consider these planets as validation of core accretion plus gas capture because it is often assumed that disk instability planets should be of nebular composition. We show in this paper, to the contrary, that gas giants that form by disk instability can have a variety of heavy element compositions, ranging from sub- to super-nebular values. High levels of enrichment can be achieved through one or multiple mechanisms, including enrichment at birth, planetesimal capture, and differentiation plus tidal stripping. As a result, the metallicity of an individual gas giant cannot be used to discriminate between gas giant formation modes.

If the two planets in the HAT-P-13 system are coplanar, the orbital states provide a probe of the internal planetary structure. Previous analyses of radial velocity and transit timing data of the system suggested that the observational constraints on the orbital states were rather small. We reanalyze the available data, treating the jitter as an unknown MCMC parameter, and find that a wide range of jitter values are plausible, hence the system parameters are less well constrained than previously suggested. For slightly increased levels of jitter ($\sim 4.5\,m\,s^{-1}$) the eccentricity of the inner planet can be in the range $0<e_{inner}<0.07$, the period and eccentricity of the outer planet can be $440<P_{outer}<470$ days and $0.55<e_{outer}<0.85$ respectively, while the relative pericenter alignment, $\eta$, of the planets can take essentially any value $-180^{\circ}<\eta<+180^{\circ}$. It is therefore difficult to determine whether $e_{inner}$ and $\eta$ have evolved to a fixed-point state or a limit cycle, or to use $e_{inner}$ to probe the internal planetary structure. We perform various transit timing variation (TTV) analyses, demonstrating that current constraints merely restrict $e_{outer}<0.85$, and rule out relative planetary inclinations within $\sim 2^{\circ}$ of $i_{rel}=90^{\circ}$, but that future observations could significantly tighten the restriction on both these parameters. We demonstrate that TTV profiles can readily distinguish the theoretically favored inclinations of $i_{rel}=0^{\circ}\,&\,45^{\circ}$, provided that sufficiently precise and frequent transit timing observations of HAT-P-13b can be made close to the pericenter passage of HAT-P-13c. We note the relatively high probability that HAT-P-13c transits and suggest observational dates and strategies.

Both ground and space-based transit observatories are poised to significantly increase the number of known transiting planets and the number of precisely measured transit times. The variation in a planet's transit times may be used to infer the presence of additional planets. Deducing the masses and orbital parameters of such planets from transit time variations (TTVs) alone is a rich and increasingly relevant dynamical problem. In this work, we evaluate the extent of the degeneracies in this process, systematically explore the dependence of TTV signals on several parameters and provide phase space plots that could aid observers in planning future observations. Our explorations are focused on a likely-to-be prevalent situation: a known transiting short-period Neptune or Jupiter-sized planet and a suspected external low-mass perturber on a nearly-coplanar orbit. Through approximately 10^7 N-body simulations, we demonstrate how TTV signal amplitudes may vary by orders of magnitude due to slight variations in any one orbital parameter (0.001 AU in semimajor axis, 0.005 in eccentricity, or a few degrees in orbital angles), and quantify the number of consecutive transit observations necessary in order to obtain a reasonable opportunity to characterize the unseen planet (approximately greater or equal to 50 observations). Planets in or near period commensurabilities of the form p:q, where p < 21 and q < 4, produce distinct TTV signatures, regardless of whether the planets are actually locked in a mean motion resonance. We distinguish these systems from the secular systems in our explorations. Additionally, we find that computing the autocorrelation function of a TTV signal can provide a useful diagnostic for identifying possible orbits for additional planets and suggest that this method could aid integration of TTV signals in future studies of particular exosystems.

The standard model of planet formation considers an initial phase in which planetesimals form from a dust disk, followed by a phase of mutual planetesimal-planetesimal collisions, leading eventually to the formation of planetary embryos. However, there is a potential transition phase (which we call the "snowball phase"), between the formation of the first planetesimals and the onset of mutual collisions amongst them, which has often been either ignored or underestimated in previous studies. In this snowball phase, isolated planetesimals move on Keplerian orbits and grow solely via the direct accretion of sub-cm sized dust entrained with the gas in the protoplanetary disk. Using a simplified model in which planetesimals are progressively produced from the dust, we consider the expected sizes to which the planetesimals can grow before mutual collisions commence and derive the dependence of this size on a number of critical parameters, including the degree of disk turbulence, the planetesimal size at birth and the rate of planetesimal creation. For systems in which turbulence is weak and the planetesimals are created at a low rate and with relatively small birth size, we show that the snowball growth phase can be very important, allowing planetesimals to grow by a factor of 10^6 in mass before mutual collisions take over. In such cases, the snowball growth phase can be the dominant mode to transfer mass from the dust to planetesimals. Moreover, such growth can take place within the typical lifetime of a protoplanetary gas disk. A noteworthy result is that ... ...(see the paper). For the specific case of close binaries such as Alpha Centauri ... ... (see the paper). From a more general perspective, these preliminary results suggest that an efficient snowball growth phase provides a large amount of "room at the bottom" for theories of planet formation.

We report radial velocity measurements of the G-type subgiants 24 Sextanis (=HD90043) and HD200964. Both are massive, evolved stars that exhibit periodic variations due to the presence of a pair of Jovian planets. Photometric monitoring with the T12 0.80m APT at Fairborn Observatory demonstrates both stars to be constant in brightness to <= 0.002 mag, thus strengthening the planetary interpretation of the radial velocity variations. 24 Sex b,c have orbital periods of 453.8 days and 883~days, corresponding to semimajor axes 1.333 AU and 2.08 AU, and minimum masses (Msini) 1.99 Mjup and 0.86 Mjup, assuming a stellar mass 1.54 Msun. HD200964 b,c have orbital periods of 613.8 days and 825 days, corresponding to semimajor axes 1.601 AU and 1.95 AU, and minimum masses 1.85 Mjup and 0.90 Mjup, assuming M* = 1.44 Msun. We also carry out dynamical simulations to properly account for gravitational interactions between the planets. Most, if not all, of the dynamically stable solutions include crossing orbits, suggesting that each system is locked in a mean motion resonance that prevents close encounters and provides long-term stability. The planets in the 24 Sex system likely have a period ratio near 2:1, while the HD200964 system is even more tightly packed with a period ratio close to 4:3. However, we caution that further radial velocity observations and more detailed dynamical modelling will be required to provide definitive and unique orbital solutions for both cases, and to determine whether the two systems are truly resonant.

We perform numerical calculations of the expected transit timing variations (TTVs) induced on a Hot-Jupiter by an Earth-mass perturber. Motivated by the recent discoveries of retrograde transiting planets, we concentrate on an investigation of the effect of varying relative planetary inclinations, up to and including completely retrograde systems. We find that planets in low order (E.g. 2:1) mean-motion resonances (MMRs) retain approximately constant TTV amplitudes for $0<\,^{\circ}i<170\,^{\circ}$, only reducing in amplitude for $i > 170\,^{\circ}$. Systems in higher order MMRs (E.g. 5:1) increase in TTV amplitude as inclinations increase towards $45\,^{\circ}$, becoming approximately constant for $45 < i < 135$, and then declining for $i > 135\,^{\circ}$. Planets away from resonance slowly decrease in TTV amplitude as inclinations increase from 0 to 180, where-as planets adjacent to resonances can exhibit a huge range of variability in TTV amplitude as a function of both eccentricity and inclination. For highly retrograde systems ($135\,^{\circ} < i \leq 180\,^{\circ}$), TTV signals will be undetectable across almost the entirety of parameter space, with the exceptions occurring when the perturber has high eccentricity or is very close to a MMR. This high inclination decrease in TTV amplitude (on and away from resonance) is important for the analysis of the known retrograde and multi-planet transiting systems, as inclination effects need to be considered if TTVs are to be used to exclude the presence of any putative planetary companions: absence of evidence is not evidence of absence.

Models for the steady state collisional evolution of low eccentricity planetesimal belts identify debris disks with hot dust at 1AU, like eta Corvi and HD69830, as anomalous since collisional processing should have removed most of the planetesimal mass over their >1 Gyr lifetimes. This paper looks at the effect of large planetesimal eccentricities (e>>0.3) on their collisional lifetime and the amount of mass that can remain at late times M_{late}. For an axisymmetric planetesimal disk with common pericentres and eccentricities e, we find that M_{late} \propto e^{-5/3}(1+e)^{4/3}(1-e)^{-3}. For a scattered disk-like population (i.e., common pericentres), in the absence of dynamical evolution, the mass evolution at late times would be as if only planetesimals with the largest eccentricity were present. Despite the increased remaining mass, higher eccentricities do not increase the hot emission from the collisional cascade until e>0.99, partly because most collisions occur near pericentre thus increasing the dust blow-out diameter. However, at high eccentricities (e>0.97) the blow-out population extending out from pericentre may be detectable above the collisional cascade; higher eccentricities also increase the probability of witnessing a recent collision. All of the imaging and spectroscopic constraints for eta Corvi can be explained with a single planetesimal population with pericentre at 0.75AU, apocentre at 150AU, and mass 5M_\oplus; however, the origin of such a high eccentricity population remains challenging. The mid-IR excess to HD69830 can be explained by the ongoing destruction of a debris belt produced in a recent collision in an eccentric planetesimal belt, but the lack of far-IR emission requires small bound grains to be absent from the parent planetesimal belt, possibly due to sublimation.

We consider the formation and migration of protoplanetary embryos in disks around the stars in tight binary systems (separations ~ 20 AU. In such systems, the initial stages of runaway embryo formation are expected to only take place within some critical disk radius a_{crit}, due to the perturbing effect of the binary companions (Thebault et al. 2009). We perform n-body simulations of the evolution of such a population of inner-disk embryos surrounded by an outer-disk of smaller planetesimals. Taking Alpha Centauri-B as our fiducial reference example in which a_{crit} ~ 0.7 AU, and using a Minimum Mass Nebular Model with $\Sigma \propto a^{-3/2}$, we find that within 10^6 yrs (10^7 yrs), systems will on average contain embryos which have migrated out to 0.9 AU (1.2 AU), with the average outer-most body having a mass of 0.2 M_{earth} 0.4 M_{earth}. Changes to increase the surface density of solids or to use a flatter profile both produce increased embryo migration and growth. At a given time, the relative change in semi-major axis of the outer-most embryo in these simulations is found to be essentially independent of a_{crit}, and we note that little further embryo migration takes place beyond 10^7 years. We conclude that the suppression of runaway growth outside a_{crit} does not mean that the habitable zones in such tight binary systems will be devoid of detectable, terrestrial mass planets, even if a_{crit} lies significantly interior to the inner edge of the habitable zone.

HD 69830 exhibits radial velocity variations attributed to three planets as well as infrared emission attributed to a warm debris disk. Previous studies have developed models for the planet migration and mass growth (Alibert et al. 2006) and the replenishment of warm grains (Wyatt et al. 2007). We perform n-body integrations in order to explore the implications of these models for: 1) the excitation of planetary eccentricity, 2) the accretion and clearing of a putative planetesimal disk, 3) the distribution of planetesimal orbits following migration, and 4) the implications for the origin of the IR emission. We find that: i) It is not possible to explain the observed planetary eccentricities (e~0.1) purely as the result of planetary perturbations during migration unless the planetary system is nearly face-on. ii) The rate of accretion of planetesimals onto planets in our n-body simulations is significantly different to that assumed in the semi-analytic models, suggesting that one cannot successfully treat planetesimal accretion in the simplified manner of Alibert et al. (2006). iii) Eccentricity damping of planetesimals does not act as an insurmountable obstacle to the existence of an excited eccentric disk: All simulations result in at least 25 Earth-masses of material remaining bound in the region ~1-9 AU, even after all three planets have migrated through the region. iv) Gas drag works to size-sort the planetesimals, with the largest bodies preferentially occupying the highest eccentricity and longest-lived orbits. Further work will be required to understand whether these eccentricity distributions are high enough to explain the level of dust emission observed despite mass loss via steady state collisional evolution. [abridged]

Recent observations point to the presence of structured dust grains in the discs surrounding young brown dwarfs, thus implying that the first stages of planet formation take place also in the sub-stellar regime. Here, we investigate the potential for planet formation around brown dwarfs and very low mass stars according to the sequential core accretion model of planet formation. We find that, for a brown dwarfs of mass 0.05M_{\odot}, our models predict a maximum planetary mass of ~5M_{\oplus}, orbiting with semi-major axis ~1AU. However, we note that the predictions for the mass - semi-major axis distribution are strongly dependent upon the models chosen for the disc surface density profiles and the assumed distribution of disc masses. In particular, if brown dwarf disc masses are of the order of a few Jupiter masses, Earth-mass planets might be relatively frequent, while if typical disc masses are only a fraction of Jupiter mass, we predict that planet formation would be extremely rare in the sub-stellar regime. As the observational constraints on disc profiles, mass dependencies and their distributions are poor in the brown dwarf regime, we advise caution in validating theoretical models only on stars similar to the Sun and emphasise the need for observational data on planetary systems around a wide range of stellar masses. We also find that, unlike the situation around solar-like stars, Type-II migration is totally absent from the planet formation process around brown dwarfs, suggesting that any future observations of planets around brown dwarfs would provide a direct measure of the role of other types of migration.

We present the results of photometric measurements from images of the LMC cluster ESO 121-SC03 taken with HST/ACS. Our resulting CMD reaches 3 magnitudes below the main-sequence turn-off, representing the deepest observation of this cluster to date. We also present similar photometry from ACS imaging of the accreted Sagittarius dSph cluster Pal. 12, used in this work as a comparison cluster. From analysis of its CMD, we obtain estimates for the metallicity and reddening of ESO 121-SC03: [Fe/H] = -0.97 +/- 0.10 and E(V-I) = 0.04 +/- 0.02, in excellent agreement with previous studies. The observed horizontal branch level in ESO 121-SC03 suggests this cluster may lie 20 per cent closer to us than does the centre of the LMC. ESO 121-SC03 also possesses a significant population of blue stragglers, which we briefly discuss. Our new photometry allows us to undertake a detailed study of the age of ESO 121-SC03 relative to Pal. 12 and the Galactic globular cluster 47 Tuc. Our main result is that ESO 121-SC03 is 73 +/- 4 per cent the age of 47 Tuc, and 91 +/- 5 per cent the age of Pal. 12. Pal. 12 is 79 +/- 6 per cent as old as 47 Tuc, consistent with previous work. Our result corresponds to an absolute age for ESO 121-SC03 in the range 8.3-9.8 Gyr depending on the age assumed for 47 Tuc, therefore confirming ESO 121-SC03 as the only known cluster to lie squarely within the LMC age gap. We briefly discuss a suggestion from earlier work that ESO 121-SC03 may have been accreted into the LMC system.