Search

Nearly one-third of objects occupying the most circular, coplanar Kuiper belt orbits (the cold classical belt) are binary, and several percent of them are "ultra-wide" binaries (UWBs): 100-km-sized companions spaced by tens of thousands of km. UWBs are dynamically fragile, and their existence is thought to constrain early Solar System processes and conditions. However, we demonstrate that UWBs can instead attain their wide architectures well after the Solar System's earliest epochs, when Neptune's orbital migration implants the modern non-cold, or "dynamic", Kuiper belt population. During this implantation, cold classical belt binaries are likely to have close encounters with many planetesimals scattered across the region, which can efficiently dissociate any existing UWBs and widen a small fraction of tighter binaries into UWB-like arrangements. Thus, today's UWBs may not be primordial and cannot be used to constrain the early Solar System as directly as previously surmised., Comment: 27 pages, 12 figures

We perform simulations here that include the gravitational effects of the primordial planetesimal belt consisting of ~10^5 massive bodies. In our simulations, Neptune unlocks from resonance with the other giant planets and begins to migrate outward due to interactions with planetesimals before a planetary orbital instability is triggered, and afterward, residual Neptunian migration completes the formation of the modern Kuiper belt. Our present work exhibits a number of notable differences from prior work. First, Neptune's planetary resonance unlocking requires the Neptunian 3:2 mean motion resonance to sweep much of the primordial disk interior to 30 au prior to the giant planet instability. The pre-instability population of planetesimals is consequently lower in semimajor axis, eccentricity, and inclination, and this effect persists after the instability. Second, direct scattering between Pluto-mass bodies and other small bodies removes material from Neptunian resonances more efficiently than resonant dropout resulting from small changes in Neptune's semimajor axis during scattering between Pluto-mass bodies and Neptune. Thus, the primordial population of Pluto-mass bodies may be as few as ~200 objects. Finally, our simulation end states display a wide variety of orbital distributions, and clear relationships between final bulk Kuiper belt properties and Neptune's migration or initial planetesimal properties largely elude us. In particular, we find that the rapid, stochastic planetary orbital evolution occurring during the giant planet instability can significantly alter final Kuiper belt properties such as its inclination dispersion and the prominence of resonant populations. This complicates using modern Kuiper belt properties to confidently constrain early solar system events and conditions, including planetary orbital migration and the primordial Kuiper belt's characteristics., Comment: Accepted to Icarus; 20 pages, 8 figures, 4 tables

Reconstructions of the paleoclimate indicate that ancient climatic fluctuations on Earth are often correlated with variations in its orbital elements. However, the chaos inherent in the solar system's orbital evolution prevents numerical simulations from confidently predicting Earth's past orbital evolution beyond 50-100 Myrs. Gravitational interactions among the Sun's planets and asteroids are believed to set this limiting time horizon, but most prior works approximate the solar system as an isolated system and neglect our surrounding Galaxy. Here we present simulations that include the Sun's nearby stellar population, and we find that close-passing field stars alter our entire planetary system's orbital evolution via their gravitational perturbations on the giant planets. This shortens the timespan over which Earth's orbital evolution can be definitively known by a further ~10%. In particular, in simulations that include an exceptionally close passage of the Sun-like star HD 7977 2.8 Myrs ago, new sequences of Earth's orbital evolution become possible in epochs before ~50 Myrs ago, which includes the Paleocene-Eocene Thermal Maximum. Thus, simulations predicting Earth's past orbital evolution before ~50 Myrs ago must consider the additional uncertainty from passing stars, which can open new regimes of past orbital evolution not seen in previous modeling efforts., Comment: Accepted to ApJL; 11 pages, 5 figures

We present a statistical characterization of circumstellar disk orientations toward 12 protostellar multiple systems in the Perseus molecular cloud using the Atacama Large Millimeter/submillimeter Array at Band 6 (1.3 mm) with a resolution of 25 mas (8 au). This exquisite resolution enabled us to resolve the compact inner disk structures surrounding the components of each multiple system and to determine the projected 3-D orientation of the disks (position angle and inclination) to high precision. We performed a statistical analysis on the relative alignment of disk pairs to determine whether the disks are preferentially aligned or randomly distributed. We considered three subsamples of the observations selected by the companion separations, a <100 au, a >500 au, and a < 10,000 au. We found for the compact (< 100 au) subsample, the distribution of orientation angles is best described by an underlying distribution of preferentially aligned sources (within 30deg) but does not rule out distributions with 40% misaligned sources. The wide companion (>500 au) subsample appears to be consistent with a distribution of 40%-80% preferentially aligned sources. Similarly, the full sample of systems with companions (a< 10, 000 au) is most consistent with a fractional ratio of at most 80% preferentially aligned source and rules out purely randomly aligned distributions. Thus our results imply the compact sources (<100 au) and the wide companions (>500 au) are statistically different.

Given the inexorable increase in the Sun's luminosity, Earth will exit the habitable zone in ~1 Gyr. There is a negligible chance that Earth's orbit will change during that time through internal Solar System dynamics. However, there is a ~1% chance per Gyr that a star will pass within 100 au of the Sun. Here, we use N-body simulations to evaluate the possible evolutionary pathways of the planets under the perturbation from a close stellar passage. We find a ~92% chance that all eight planets will survive on orbits similar to their current ones if a star passes within 100 au of the Sun. Yet a passing star may disrupt the Solar System, by directly perturbing the planets' orbits or by triggering a dynamical instability. Mercury is the most fragile, with a destruction rate (usually via collision with the Sun) higher than that of the four giant planets combined. The most probable destructive pathways for Earth are to undergo a giant impact (with the Moon or Venus) or to collide with the Sun. Each planet may find itself on a very different orbit than its present-day one, in some cases with high eccentricities or inclinations. There is a small chance that Earth could end up on a more distant (colder) orbit, through re-shuffling of the system's orbital architecture, ejection into interstellar space (or into the Oort cloud), or capture by the passing star. We quantify plausible outcomes for the post-flyby Solar System., Comment: 13 pages. Accepted for publication in MNRAS. Related blog post at https://planetplanet.net/2023/11/21/could-a-stellar-flyby-save-earth-from-impending-doom/

Dynamical instabilities among giant planets are thought to be nearly ubiquitous, and culminate in the ejection of one or more planets into interstellar space. Here we perform N-body simulations of dynamical instabilities while accounting for torques from the galactic tidal field. We find that a fraction of planets that would otherwise have been ejected are instead trapped on very wide orbits analogous to those of Oort cloud comets. The fraction of ejected planets that are trapped ranges from 1-10%, depending on the initial planetary mass distribution. The local galactic density has a modest effect on the trapping efficiency and the orbital radii of trapped planets. The majority of Oort cloud planets survive for Gyr timescales. Taking into account the demographics of exoplanets, we estimate that one in every 200-3000 stars could host an Oort cloud planet. This value is likely an overestimate, as we do not account for instabilities that take place at early enough times to be affected by their host stars' birth cluster, or planet stripping from passing stars. If the Solar System's dynamical instability happened after birth cluster dissolution, there is a ~7% chance that an ice giant was captured in the Sun's Oort cloud., Comment: MNRAS Letters, in press. Blog post about paper at https://planetplanet.net/2023/06/21/oort-cloud-exoplanets/

The chemical compositions of Earth's core and mantle provide insight into the processes that led to their formation. N-body simulations, on the other hand, generally do not contain chemical information, and seek to only reproduce the masses and orbits of the terrestrial planets. These simulations can be grouped into four potentially viable scenarios of Solar System formation (Classical, Annulus, Grand Tack, and Early Instability) for which we compile a total of 433 N-body simulations. We relate the outputs of these simulations to the chemistry of Earth's core and mantle using a melt-scaling law combined with a multi-stage model of core formation. We find the compositions of Earth analogs to be largely governed by the fraction of equilibrating embryo cores and the initial embryo masses in N-body simulations. Simulation type may be important when considering magma ocean lifetimes, where Grand Tack simulations have the largest amounts of material accreted after the last giant impact. However, we cannot rule out any accretion scenarios or initial embryo masses due to the sensitivity of Earth's mantle composition to different parameters and the stochastic nature of N-body simulations. Comparing the last embryo impacts experienced by Earth analogs to specific Moon-forming scenarios, we find the characteristics of the Moon-forming impact are dependent on the initial conditions in N-body simulations where larger initial embryo masses promote larger and slower Moon-forming impactors. Mars-sized initial embryos are most consistent with the canonical hit-and-run scenario onto a solid mantle. Our results suggest that constraining the fraction of equilibrating impactor core and the initial embryo masses in N-body simulations could be significant for understanding both Earth's accretion history and characteristics of the Moon-forming impact.

The inner solar system's modern orbital architecture provides inferences into the epoch of terrestrial planet formation; a ~100 Myr time period of planet growth via collisions with planetesimals and other proto-planets. While classic numerical simulations of this scenario adequately reproduced the correct number of terrestrial worlds, their semi-major axes and approximate formation timescales, they struggled to replicate the Earth-Mars and Venus-Mercury mass ratios. In a series of past independent investigations, we demonstrated that Mars' mass is possibly the result of Jupiter and Saturn's early orbital evolution, while Mercury's diminutive size might be the consequence of a primordial mass deficit in the region. Here, we combine these ideas in a single modeled scenario designed to simultaneously reproduce the formation of all four terrestrial planets and the modern orbits of the giant planets in broad strokes. By evaluating our Mercury analogs' core mass fractions, masses, and orbital offsets from Venus, we favor a scenario where Mercury forms through a series of violent erosive collisions between a number of ~Mercury-mass embryos in the inner part of the terrestrial disk. We also compare cases where the gas giants begin the simulation locked in a compact 3:2 resonant configuration to a more relaxed 2:1 orientation and find the former to be more successful. In 2:1 cases, the entire Mercury-forming region is often depleted due to strong sweeping secular resonances that also tend to overly excite the orbits of Earth and Venus as they grow. While our model is quite successful at replicating Mercury's massive core and dynamically isolated orbit, the planets' low mass remains extremely challenging to match. Finally, we discuss the merits and drawbacks of alternative evolutionary scenarios and initial disk conditions., Comment: 23 pages, 9 figures, 3 tables, accepted for publication in Icarus

Within the dynamically cold low inclination portion of the Classical Kuiper Belt, there exists a population of weakly bound binary systems with a number of unusual properties; most notable of which is their extremely wide orbital separations; beyond 7% of their Hill radii. The stability and evolution of these Ultra-Wide Trans-Neptunian Binaries (TNBs) have, in the past, been studied extensively under the assumption that the primary evolving mechanisms are interactions between the binary components and impacting Trans-Neptunian Objects (TNOs). Here, we instead study their evolution as driven by the gravitational perturbations of close passing but non-impacting TNOs. By simulating these passages, we show that the aggregate effects of encounters over billions of years have a significant effect on Kuiper Belt binary evolution. Such processes can lead to tight binaries widening significantly over time, approaching and sometimes surpassing the separation of the widest known TNBs. We also find that the eccentricity and inclination distributions of observed Ultra-Wide TNBs can be sampled from such widened binaries. While we are unable to produce enough wide binaries to explain their abundance, the orbital properties of ones we do produce are consistent with known wide binaries.

The majority of binary star systems that host exoplanets will spend the first portion of their lives within a star-forming cluster that may drive dynamical evolution of the binary-planet system. We perform numerical simulations of S-type planets, with masses and orbital architecture analogous to the solar system's 4 gas giants, orbiting within the influence of a 0.5 solar-mass binary companion. The binary-planet system is integrated simultaneously with an embedded stellar cluster environment. ~10% of our planetary systems are destabilized when perturbations from our cluster environment drive the binary periastron toward the planets. This destabilization occurs despite all of our systems being initialized with binary orbits that would allow stable planets in the absence of the cluster. The planet-planet scattering triggered in our systems typically results in the loss of lower mass planets and the excitement of the eccentricities of surviving higher mass planets. Many of our planetary systems that go unstable also lose their binary companions prior to cluster dispersal and can therefore masquerade as hosts of eccentric exoplanets that have spent their entire histories as isolated stars. The cluster-driven binary orbital evolution in our simulations can also generate planetary systems with misaligned spin-orbit angles. This is typically done as the planetary system precesses as a rigid disk under the influence of an inclined binary, and those systems with the highest spin-orbit angles should often retain their binary companion and possess multiple surviving planets., Comment: 20 pages, 7 figures, 3 tables, accepted for publication in MNRAS

The Oort cloud and the scattered disk are the two primary reservoirs for long-period and short-period comets, respectively. In this review, we assess the known observational constraints on these reservoirs' properties and their formation. In addition, we discuss how the early orbital evolution of the giant planets generated the modern scattered disk from the early, massive planetesimal disk and how $\sim$5\% of this material was captured into the Oort cloud. We review how the Sun's birth environment and dynamical history within the Milky Way alters the formation and modern structure of the Oort cloud. Finally, we assess how the coming decade's anticipated observing campaigns may provide new insights into the formation and properties of the Oort cloud and scattered disk., Comment: 25 pages, 13 figures; Chapter in press for the book Comets III, edited by K. Meech and M. Combi, University of Arizona Press

The discovery probability of long-period comets (LPCs) passing near the Sun is highest during their first passage and then declines, or fades, during subsequent return passages. Comet fading is largely attributed to devolatilization and fragmentation via thermal processing within 2-3 au of the Sun (1 au being the Earth-Sun distance). Here our numerical simulations show that comet observing campaigns miss vast numbers of LPCs making returning passages through the Saturn region (near 10 au) because these comets fade during prior, even more distant passages exterior to Saturn and thus elude detection. Consequently, comet properties significantly evolve at solar distances much larger than previously considered, and this offers new insights into the physical and dynamical properties of LPCs, both near and far from Earth., Comment: 23 pages, 5 figures, 1 table; published in Science Advances

A distant, massive planet in the outer solar system has recently been proposed to explain some observed features of extreme trans-Neptunian objects (TNOs). Here we use N-body simulations of the formation of the Kuiper belt and Oort cloud as well as a survey simulator to compare models of the solar system with and without a 9th planet to one another as well as to observations. The main mechanism for TNOs to be deposited into the distant ($a$ > au), detached ($q$ > au) region of the Kuiper Belt in the 8-planet model is Kozai-Lidov oscillation of objects in mean motion resonances (MMR) with Neptune. This effect does not deposit low-inclination ($i \lesssim$ 20{\deg}) objects into this region. However, we find that the 9th planet generates a group of distant, detached TNOs at low inclinations that are not present in the 8-planet model. This disparity between the 8-planet and the 9-planet models could provide a strong constraint on a possible planet 9 with further detections of TNOs in the distant, detached region of the Kuiper Belt., Comment: 8 pages, 3 figures, accepted to ApJL

The formation of the solar system's giant planets predated the ultimate epoch of massive impacts that concluded the process of terrestrial planet formation. Following their formation, the giant planets' orbits evolved through an episode of dynamical instability. Several qualities of the solar system have recently been interpreted as evidence of this event transpiring within the first ~100 Myr after the Sun's birth; around the same time as the final assembly of the inner planets. In a series of recent papers we argued that such an early instability could resolve several problems revealed in classic numerical studies of terrestrial planet formation; namely the small masses of Mars and the asteroid belt. In this paper, we revisit the early instability scenario with a large suite of simulations specifically designed to understand the degree to which Earth and Mars' formation are sensitive to the specific evolution of Jupiter and Saturn's orbits. By deriving our initial terrestrial disks directly from recent high-resolution simulations of planetesimal accretion, our results largely confirm our previous findings regarding the instability's efficiency of truncating the terrestrial disk outside of the Earth-forming region in simulations that best replicate the outer solar system. Moreover, our work validates the primordial 2:1 Jupiter-Saturn resonance within the early instability framework as a viable evolutionary path for the solar system. While our simulations elucidate the fragility of the terrestrial system during the epoch of giant planet migration, many realizations yield outstanding solar system analogs when scrutinized against a number of observational constraints. Finally, we highlight the inability of models to form adequate Mercury-analogs and the low eccentricities of Earth and Venus as the most significant outstanding problems for future numerical studies to resolve., Comment: 24 pages, 10 figures, 6 tables, accepted for publication in Icarus

In a recent paper we proposed that the giant planets' primordial orbits may have been eccentric (~0.05), and used a suite of dynamical simulations to show outcomes of the giant planet instability that are consistent with their present-day orbits. In this follow-up investigation, we present more comprehensive simulations incorporating superior particle resolution, longer integration times, and eliminating our prior means of artificially forcing instabilities to occur at specified times by shifting a planets' position in its orbit. While we find that the residual phase of planetary migration only minimally alters the the planets' ultimate eccentricities, our work uncovers several intriguing outcomes in realizations where Jupiter and Saturn are born with extremely large eccentricities (~0.10 and ~0.25, respectively). In successful simulations, the planets' orbits damp through interactions with the planetesimal disk prior to the instability, thus loosely replicating the initial conditions considered in our previous work. Our results therefore suggest an even wider range of plausible evolutionary pathways are capable of replicating Jupiter and Saturn's modern orbital architecture., Comment: 12 pages, 3 figures, 2 tables, accepted for publication in Icarus

The non-resonant Kuiper belt objects (KBOs) between the 3:2 and 2:1 Neptunian mean motion resonances can be largely divided between a cold classical belt (CCB) and a hot classical belt (HCB). A notable difference between these two subpopulations is the prevalence of widely spaced, equal-mass binaries in the CCB and a much smaller but non-zero number in the HCB. The primary reason for this difference in binary rate remains unclear. Here using N-body simulations we examine whether close encounters with the giant planets during an early outer solar system instability may have disrupted primordial Kuiper Belt binaries that existed within the primordial Kuiper belt before they attained HCB orbits. We find that such encounters are very effective at disrupting binaries down to separations of ~1% of their Hill radius (as measured in the modern Kuiper belt), potentially explaining the paucity of widely spaced, equal mass binaries in the modern HCB. Moreover, we find that the widest binaries observed in the modern HCB are quite unlikely to survive planetary encounters, but these same planetary encounters can widen a small subset of tighter binaries to give rise to the small population of very wide binaries seen in today's HCB., Comment: Accepted 2021 April 28. Received 2021 April 26; in original form 2021 February 22. 5 pages, 3 figures, 2 tables

The orbital architecture of the Solar System is thought to have been sculpted by a dynamical instability among the giant planets. During the instability a primordial outer disk of planetesimals was destabilized and ended up on planet-crossing orbits. Most planetesimals were ejected into interstellar space but a fraction were trapped on stable orbits in the Kuiper belt and Oort cloud. We use a suite of N-body simulations to map out the diversity of planetesimals' dynamical pathways. We focus on two processes: tidal disruption from very close encounters with a giant planet, and loss of surface volatiles from repeated passages close to the Sun. We show that the rate of tidal disruption is more than a factor of two higher for ejected planetesimals than for surviving objects in the Kuiper belt or Oort cloud. Ejected planetesimals are preferentially disrupted by Jupiter and surviving ones by Neptune. Given that the gas giants contracted significantly as they cooled but the ice giants did not, taking into account the thermal evolution of the giant planets decreases the disruption rate of ejected planetesimals. The frequency of volatile loss and extinction is far higher for ejected planetesimals than for surviving ones and is not affected by the giant planets' contraction. Even if all interstellar objects were ejected from Solar System-like systems, our analysis suggests that their physical properties should be more diverse than those of Solar System small bodies as a result of their divergent dynamical histories. This is consistent with the characteristics of the two currently-known interstellar objects., Comment: ApJ Letters, in press. 6 pages, 4 figures

An episode of dynamical instability is thought to have sculpted the orbital structure of the outer solar system. When modeling this instability, a key constraint comes from Jupiter's fifth eccentric mode (quantified by its amplitude M55), which is an important driver of the solar system's secular evolution. Starting from commonly-assumed near-circular orbits, the present-day giant planets' architecture lies at the limit of numerically generated systems, and M55 is rarely excited to its true value. Here we perform a dynamical analysis of a large batch of artificially triggered instabilities, and test a variety of configurations for the giant planets' primordial orbits. In addition to more standard setups, and motivated by the results of modern hydrodynamical simulations of the giant planets' evolution within the primordial gaseous disk, we consider the possibility that Jupiter and Saturn emerged from the nebular gas locked in 2:1 resonance with non-zero eccentricities. We show that, in such a scenario, the modern Jupiter-Saturn system represents a typical simulation outcome, and M55 is commonly matched. Furthermore, we show that Uranus and Neptune's final orbits are determined by a combination of the mass in the primordial Kuiper belt and that of an ejected ice giant., Comment: 35 pages, 21 figures, 8 tables, accepted for publication in AJ

The longitudes of perihelia and orbital poles of the solar system's dozen or so most remote detected objects are clustered in a manner inconsistent with that of a random sample of uniformly distributed orbits. While small number statistics and observational biases may explain these features, the statistical significance of the clustering has led to the recent development of the "Planet Nine hypothesis." In the proposed scenario, orbits in the distant solar system are shepherded via secular perturbations from an undetected massive planet on an eccentric orbit. However, the precession of perihelia and nodes in the outer Kuiper Belt and inner Oort Cloud are also affected by the the giant planets, passing stars, and the galactic tide. We perform a large suite of numerical simulations designed to study the orbital alignment of Extreme Trans-Neptunian Objects (ETNOs) and Inner Oort Cloud Objects (IOCOs). In our various integrations that include Planet Nine, we consistently find that >60% of ETNOs and IOCOs that are detectable after 4 Gyr are also anti-aligned in perihelia with the distant massive perturber. However, when we randomly select 17 objects from this sample of remaining orbits, there is significant scatter in the degree of longitude of perihelion and orbital pole clustering that might be observed. Furthermore, we argue that, in the absence of Planet Nine, 17 randomly drawn orbits should still exhibit some clustering even if the underlying distribution is uniform. Thus, we find that still more ETNO and IOCO detections are required to confidently infer the presence of Planet Nine., Comment: 11 pages, 7 figures, 4 tables, accepted for publication in AJ

The solar system's terrestrial planets are thought to have accreted over millions of years out of a sea of smaller embryos and planetesimals. Because it is impossible to know the surface density profile for solids and size frequency distribution in the primordial solar nebula, distinguishing between the various proposed evolutionary schemes has been historically difficult. Nearly all previous simulations of terrestrial planet formation assume that Moon to Mars massed embryos formed throughout the inner solar system during the primordial gas-disk phase. However, validating this assumption through models of embryo accretion is computationally challenging because of the large number of bodies required. Here, we reevaluate this problem with GPU-accelerated, direct N-body simulations of embryo growth starting from r~100 km planetesimals. We find that embryos emerging from the primordial gas phase at a given radial distance already have masses similar to the largest objects at the same semi-major axis in the modern solar system. Thus, Earth and Venus attain ~50% of their modern mass, Mars-massed embryos form in the Mars region, and Ceres-massed objects are prevalent throughout asteroid belt. Consistent with other recent work, our new initial conditions for terrestrial accretion models produce markedly improved solar system analogs when evolved through the giant impact phase of planet formation. However, we still conclude that an additional dynamical mechanism such as giant planet migration is required to prevent Earth-massed Mars analogs from growing., Comment: 16 pages, 13 figures, 4 tables, accepted for publication in The Planetary Science Journal

The asteroid belt is characterized by an extreme low total mass of material on dynamically excited orbits. The Nice Model explains many peculiar qualities of the solar system, including the belt's excited state, by invoking an orbital instability between the outer planets. However, previous studies of the Nice Model's effect on the belt's structure struggle to reproduce the innermost asteroids' orbital inclination distribution. Here, we show how the final phase of giant planet migration sculpts the asteroid belt, in particular its inclination distribution. As interactions with leftover planetesimals cause Saturn to move away from Jupiter, its rate of orbital precession slows as the two planets' mutual interactions weaken. When the planets approach their modern separation, where Jupiter completes just short of five orbits for every two of Saturn's, Jupiter's eccentric forcing on Saturn strengthens. We use numerical simulations to show that the absence of asteroids with orbits that precess between 24-28 arcsec/yr is related to the inclination problem. As Saturn's precession speeds back up, high inclination asteroids are excited on to planet crossing orbits and removed from the inner main belt. Through this process, the asteroid belt's orbital structure is reshaped, leading to markedly improved simulation outcomes., Comment: 5 pages, 4 figures, 2 tables, accepted for publication in MRNAS Letters

Most known trans-Neptunian objects (TNOs) gravitationally scattering off the giant planets have orbital inclinations consistent with an origin from the classical Kuiper belt, but a small fraction of these "scattering TNOs" have inclinations that are far too large (i > 45 deg) for this origin. These scattering outliers have previously been proposed to be interlopers from the Oort cloud or evidence of an undiscovered planet. Here we test these hypotheses using N-body simulations and the 69 centaurs and scattering TNOs detected in the Outer Solar Systems Origins Survey and its predecessors. We confirm that observed scattering objects cannot solely originate from the classical Kuiper belt, and we show that both the Oort cloud and a distant planet generate observable highly inclined scatterers. Although the number of highly inclined scatterers from the Oort Cloud is ~3 times less than observed, Oort cloud enrichment from the Sun's galactic migration or birth cluster could resolve this. Meanwhile, a distant, low-eccentricity 5 Earth-mass planet replicates the observed fraction of highly inclined scatterers, but the overall inclination distribution is more excited than observed. Furthermore, the distant planet generates a longitudinal asymmetry among detached TNOs that is less extreme than often presumed, and its direction reverses across the perihelion range spanned by known TNOs. More complete models that explore the dynamical origins of the planet are necessary to further study these features. With observational biases well-characterized, our work shows that the orbital distribution of detected scattering bodies is a powerful constraint on the unobserved distant solar system., Comment: 17 pages, 11 figures, accepted to AJ

Of the solar system's four terrestrial planets, the origin of Mercury is perhaps the most mysterious. Modern numerical simulations designed to model the dynamics of terrestrial planet formation systematically fail to replicate Mercury; which possesses just 5% the mass of Earth and the highest orbital eccentricity and inclination among the planets. However, Mercury's large iron-rich core and low volatile inventory stand out among the inner planets, and seem to imply a violent collisional origin. Because most algorithms used for simulating terrestrial accretion do not consider the effects of collisional fragmentation, it has been difficult to test these collisional hypotheses within the larger context of planet formation. Here, we analyze a large suite of terrestrial accretion models that account for the fragmentation of colliding bodies. We find that planets with core mass fractions boosted as a result of repeated hit-and-run collisions are produced in 90% of our simulations. While many of these planets are similar to Mercury in mass, they rarely lie on Mercury-like orbits. Furthermore, we perform an additional batch of simulations designed to specifically test the single giant impact origin scenario. We find less than a 1% probability of simultaneously replicating the Mercury-Venus dynamical spacing and the terrestrial system's degree of orbital excitation after such an event. While dynamical models have made great strides in understanding Mars' low mass, their inability to form accurate Mercury analogs remains a glaring problem., Comment: 11 pages, 7 figures, 2 tables, accepted for publication in AJ

Modern studies of the early solar system routinely invoke the possibility of an orbital instability among the giant planets triggered by gravitational interactions between the planets and a massive exterior disk of planetesimals. Previous works have suggested that this instability can be substantially delayed ($\sim$100s Myr) after the formation of the giant planets. Bodies in the disk are typically treated in a semi-active manner, wherein their gravitational force on the planets is included, but interactions between the planetesimals are ignored. We perform $N$-body numerical simulations using GENGA, which makes use of GPUs to allow for the inclusion of all gravitational interactions between bodies. Although our simulated Kuiper belt particles are more massive than the probable masses of real primordial Kuiper belt objects, our simulations indicate that the self-stirring of the primordial Kuiper belt is very important to the dynamics of the giant planet instability. We find that interactions between planetesimals dynamically heat the disk and typically prevent the outer solar system instability from being delayed by more than a few tens of million years after giant planet formation. Longer delays occur in a small fraction of systems that have at least 3.5 AU gaps between the planets and planetesimal disk. Our final planetary configurations match the solar system at a rate consistent with other previous works in most regards. Pre-instability heating of the disk typically yields final Jovian eccentricities comparable to the modern solar system value, which has been a difficult constraint to match in past works., Comment: 16 pages, 17 figures, 3 tables; accepted for publication in AJ

The solar system's dynamical state can be explained by an orbital instability among the giant planets. A recent model has proposed that the giant planet instability happened during terrestrial planet formation. This scenario has been shown to match the inner solar system by stunting Mars' growth and preventing planet formation in the asteroid belt. Here we present a large sample of new simulations of the "Early Instability" scenario. We use an N-body integration scheme that accounts for collisional fragmentation, and also perform a large set of control simulations that do not include an early giant planet instability. Since the total particle number decreases slower when collisional fragmentation is accounted for, the growing planets' orbits are damped more strongly via dynamical friction and encounters with small bodies that dissipate angular momentum (eg: hit-and-run impacts). Compared with simulations without collisional fragmentation, our fully evolved systems provide better matches to the solar system's terrestrial planets in terms of their compact mass distribution and dynamically cold orbits. Collisional processes also tend to lengthen the dynamical accretion timescales of Earth analogs, and shorten those of Mars analogs. This yields systems with relative growth timescales more consistent with those inferred from isotopic dating. Accounting for fragmentation is thus supremely important for any successful evolutionary model of the inner solar system., Comment: 16 pages, 7 figures, 4 tables, accepted for publication in Icarus

Containing only a few percent the mass of the moon, the current asteroid belt is around three to four orders of magnitude smaller that its primordial mass inferred from disk models. Yet dynamical studies have shown that the asteroid belt could not have been depleted by more than about an order of magnitude over the past ~4 Gyr. The remainder of the mass loss must have taken place during an earlier phase of the solar system's evolution. An orbital instability in the outer solar system occurring during the process of terrestrial planet formation can reproduce the broad characteristics of the inner solar system. Here, we test the viability of this model within the constraints of the main belt's low present-day mass and orbital structure. While previous studies modeled asteroids as massless test particles because of limited computing power, our work uses GPU (Graphics Processing Unit) acceleration to model a fully self-gravitating asteroid belt. We find that depletion in the main belt is related to the giant planets' exact evolution within the orbital instability. Simulations that produce the closest matches to the giant planets' current orbits deplete the main belt by two to three orders of magnitude. These simulated asteroid belts are also good matches to the actual asteroid belt in terms of their radial mixing and broad orbital structure., Comment: 12 pages, 7 figures, 2 tables, accepted for publication in AJ corrected minor typographical errors

Inner Oort Cloud objects (IOCs) are Trans-Plutonian for their entire orbits. They are beyond the strong gravitational influences of the known planets yet close enough to the Sun that outside forces are minimal. Here we report the discovery of the third known IOC after Sedna and 2012 VP113, called 2015 TG387. 2015 TG387 has a perihelion of $65 \pm 1$ au and semi-major axis of $1170 \pm 70$ au. The longitude of perihelion angle, $\bar{\omega}$, for 2015 TG387 is between that of Sedna and 2012 VP113, and thus similar to the main group of clustered extreme trans-Neptunian objects (ETNOs), which may be shepherded into similar orbital angles by an unknown massive distant planet, called Planet X or Planet Nine. 2015 TG387's orbit is stable over the age of the solar system from the known planets and Galactic tide. When including outside stellar encounters over 4 Gyrs, 2015 TG387's orbit is usually stable, but its dynamical evolution depends on the stellar encounter scenarios used. Surprisingly, when including a massive Planet X beyond a few hundred au on an eccentric orbit that is anti-aligned in longitude of perihelion with most of the known ETNOs, we find 2015 TG387 is typically stable for Planet X orbits that render the other ETNOs stable as well. Notably, 2015 TG387's argument of perihelion is constrained and its longitude of perihelion librates about 180 degs from Planet X's longitude of perihelion, keeping 2015 TG387 anti-aligned with Planet X over the age of the solar system. We find a power law slope near 3 for the semi-major axis distribution of IOCs, meaning there are many more high than low semi-major axis IOCs. There are about 2 million IOCs larger than 40 km, giving a mass of $10^{22}$ kg. The IOCs inclination distribution is similar to the scattered disk, with an average inclination of 19 degs., Comment: Accepted to The Astronomical Journal

The Outer Solar System Origins Survey (OSSOS), a wide-field imaging program in 2013-2017 with the Canada-France-Hawaii Telescope, surveyed 155 deg$^{2}$ of sky to depths of $m_r = 24.1$-25.2. We present 838 outer Solar System discoveries that are entirely free of ephemeris bias. This increases the inventory of trans-Neptunian objects (TNOs) with accurately known orbits by nearly 50%. Each minor planet has 20-60 Gaia/Pan-STARRS-calibrated astrometric measurements made over 2-5 oppositions, which allows accurate classification of their orbits within the trans-Neptunian dynamical populations. The populations orbiting in mean-motion resonance with Neptune are key to understanding Neptune's early migration. Our 313 resonant TNOs, including 132 plutinos, triple the available characterized sample and include new occupancy of distant resonances out to semi-major axis $a \sim 130$ au. OSSOS doubles the known population of the non-resonant Kuiper belt, providing 436 TNOs in this region, all with exceptionally high-quality orbits of $a$ uncertainty $\sigma_{a} \leq 0.1\%$; they show the belt exists from $a \gtrsim 37$ au, with a lower perihelion bound of $35$ au. We confirm the presence of a concentrated low-inclination $a\simeq 44$ au "kernel" population and a dynamically cold population extending beyond the 2:1 resonance. We finely quantify the survey's observational biases. Our survey simulator provides a straightforward way to impose these biases on models of the trans-Neptunian orbit distributions, allowing statistical comparison to the discoveries. The OSSOS TNOs, unprecedented in their orbital precision for the size of the sample, are ideal for testing concepts of the history of giant planet migration in the Solar System., Comment: Invited paper, special issue Data: Insights and Challenges in a Time of Abundance. Data tables and example survey simulator are in the supplementary materials (see arXiv source under Downloads > Other formats)

Many dynamical aspects of the solar system can be explained by the outer planets experiencing a period of orbital instability sometimes called the Nice Model. Though often correlated with a perceived delayed spike in the lunar cratering record known as the Late Heavy Bombardment (LHB), recent work suggests that this event may have occurred much earlier; perhaps during the epoch of terrestrial planet formation. While current simulations of terrestrial accretion can reproduce many observed qualities of the solar system, replicating the small mass of Mars requires modification to standard planet formation models. Here we use 800 dynamical simulations to show that an early instability in the outer solar system strongly influences terrestrial planet formation and regularly yields properly sized Mars analogs. Our most successful outcomes occur when the terrestrial planets evolve an additional 1-10 million years (Myr) following the dispersal of the gas disk, before the onset of the giant planet instability. In these simulations, accretion has begun in the Mars region before the instability, but the dynamical perturbation induced by the giant planets' scattering removes large embryos from Mars' vicinity. Large embryos are either ejected or scattered inward toward Earth and Venus (in some cases to deliver water), and Mars is left behind as a stranded embryo. An early giant planet instability can thus replicate both the inner and outer solar system in a single model., Comment: 18 pages, 9 figures, accepted for publication in Icarus

The stability limit for circumbinary planets (CBPs) is not well defined and can depend on initial parameters defining either the planetary orbit or the inner binary orbit. We expand on the work of Holman & Wiegert (1999, AJ 117, 621) to develop numerical tools for quick, easy, and accurate determination of the stability limit. The results of our simulations, as well as our numerical tools, are available to the community through $\texttt{Zenodo}$ and $\texttt{GitHub}$, respectively. We employ a grid interpolation method based on $\sim$150 million full N-body simulations of initially circular, coplanar systems and compare to the 9 known Kepler CBP systems. Using a formalism from planet packing studies, we find that 55% of the Kepler CBP systems allow for an additional equal-mass planet to potentially exist on an interior orbit relative to the observed planet. Therefore, we do $\textit{not}$ find strong evidence for a pile-up in the Kepler CBP systems and more detections are needed to adequately characterize the formation mechanisms for the CBP population. Observations from the Transiting Exoplanet Survey Satellite are expected to substantially increase the number of detections using the unique geometry of CBP systems, where multiple transits can occur during a single conjunction., Comment: 18 pages, 8 figures, 4 tables; Accepted in Astrophysical Journal

We extend our study of the extent of the regions within the $\alpha$ Centauri AB star system where small planets are able to orbit for billion-year timescales (Quarles & Lissauer 2016, AJ 151, 111) to investigate the effects of minimizing the forced eccentricity of initial trajectories. We find that initially prograde, circumstellar orbits require a piecewise quadratic function to accurately approximate forced eccentricity as a function of semimajor axis, but retrograde orbits can be modeled using a linear function. Circumbinary orbits in the $\alpha$ Centauri AB system are less affected by the forced eccentricity. {Planets on circumstellar orbits that begin with eccentricity vectors near their forced values are generally stable, up to $\sim$$10^9$ yr, out to a larger semimajor axis than are planets beginning on circular orbits. The amount by which the region of stability expands is much larger for retrograde orbits than it is for prograde orbits. The location of the stability boundary for two planet systems on prograde, circular orbits is much more sensitive to the initial eccentricity state than it is for analogous single planet systems., Comment: 11 pages, 10 figures, 3 tables; Accepted in Astronomical Journal

The galactic environment has been suspected to influence planetary habitability in many ways. Very metal-poor regions of the Galaxy, or those largely devoid of atoms more massive than H and He, are thought to be unable to form habitable planets. Moreover, if such planets do form, the young system is subjected to close stellar passages while it resides in its stellar birth cluster. Various potential hazards remain after clusters disperse. For instance, central galactic regions may present risks to habitability via nearby supernovae, gamma ray bursts (GRBs), and frequent comet showers. In addition, planets residing within very wide binary star systems are affected by the Galaxy, as local gravitational perturbations from the Galaxy can increase the binary's eccentricity until it destabilizes the planets it hosts. Here we review the most recent work on the main galactic influences over planetary habitability. Although there must be some metallicity limit below which rocky planets cannot form, recent exoplanet surveys show that they form around stars with a very large range of metallicities. Once formed, the probability of star clusters destabilizing planetary systems only becomes high for rare, extremely long-lived clusters. Regarding threats to habitability from supernovae, GRBs, and comet showers, many recent studies suggest that their hazards are more limited than originally thought. Finally, denser regions of the Galaxy enhance the threat that very wide binary companions pose to planetary habitability, but the probability that a very wide binary star disrupts habitability will always be substantially below 100% for any environment. While some Milky Way regions must be more hospitable to habitable planets than others, it is difficult to state that habitable planets are confined to any well-defined region of the Galaxy or that any other particular region of the Galaxy is uninhabitable., Comment: Invited review chapter, accepted for publication in the "Handbook of Exoplanets"; 19 pages; 2 figures

Fomalhaut A is among the most well-studied nearby stars and has been discovered to possess a putative planetary object as well as a remarkable eccentric dust belt. This eccentric dust belt has often been interpreted as the dynamical signature of one or more planets that elude direct detection. However, the system also contains two other stellar companions residing ~100,000 AU from Fomalhaut A. We have designed a new symplectic integration algorithm to model the evolution of Fomalhaut A's planetary dust belt in concert with the dynamical evolution of its stellar companions to determine if these companions are likely to have generated the dust belt's morphology. Using our numerical simulations, we find that close encounters between Fomalhaut A and B are expected, with a ~25% probability that the two stars have passed within at least 400 AU of each other at some point. Although the outcomes of such encounter histories are extremely varied, these close encounters nearly always excite the eccentricity of Fomalhaut A's dust belt and occasionally yield morphologies very similar to the observed belt. With these results, we argue that close encounters with Fomalhaut A's stellar companions should be considered a plausible mechanism to explain its eccentric belt, especially in the absence of detected planets capable of sculpting the belt's morphology. More broadly, we can also conclude from this work that very wide binary stars may often generate asymmetries in the stellar debris disks they host., Comment: Accepted to MNRAS, 22 pages, 15 figures, 2 appendices

The accumulating, but small, set of large semi-major axis trans-Neptunian objects (TNOs) shows an apparent clustering in the orientations of their orbits. This clustering must either be representative of the intrinsic distribution of these TNOs, or else arise as a result of observation biases and/or statistically expected variations for such a small set of detected objects. The clustered TNOs were detected across different and independent surveys, which has led to claims that the detections are therefore free of observational bias. This apparent clustering has led to the so-called "Planet 9" hypothesis that a super-Earth currently resides in the distant solar system and causes this clustering. The Outer Solar System Origins Survey (OSSOS) is a large program that ran on the Canada-France-Hawaii Telescope from 2013--2017, discovering more than 800 new TNOs. One of the primary design goals of OSSOS was the careful determination of observational biases that would manifest within the detected sample. We demonstrate the striking and non-intuitive biases that exist for the detection of TNOs with large semi-major axes. The eight large semi-major axis OSSOS detections are an independent dataset, of comparable size to the conglomerate samples used in previous studies. We conclude that the orbital distribution of the OSSOS sample is consistent with being detected from a uniform underlying angular distribution., Comment: Accepted for publication

We report the discovery of the minor planet 2013 SY$_{99}$, on an exceptionally distant, highly eccentric orbit. With a perihelion of 50.0 au, 2013 SY$_{99}$'s orbit has a semi-major axis of $730 \pm 40$ au, the largest known for a high-perihelion trans-Neptunian object (TNO), well beyond those of (90377) Sedna and 2012 VP$_{113}$. Yet, with an aphelion of $1420 \pm 90$ au, 2013 SY$_{99}$'s orbit is interior to the region influenced by Galactic tides. Such TNOs are not thought to be produced in the current known planetary architecture of the Solar System, and they have informed the recent debate on the existence of a distant giant planet. Photometry from the Canada-France-Hawaii Telescope, Gemini North and Subaru indicate 2013 SY$_{99}$ is $\sim 250$ km in diameter and moderately red in colour, similar to other dynamically excited TNOs. Our dynamical simulations show that Neptune's weak influence during 2013 SY$_{99}$'s perihelia encounters drives diffusion in its semi-major axis of hundreds of astronomical units over 4 Gyr. The overall symmetry of random walks in semi-major axis allow diffusion to populate 2013 SY$_{99}$'s orbital parameter space from the 1000-2000 au inner fringe of the Oort cloud. Diffusion affects other known TNOs on orbits with perihelia of 45 to 49 au and semi-major axes beyond 250 au, providing a formation mechanism that implies an extended population, gently cycling into and returning from the inner fringe of the Oort cloud., Comment: First reviewer report comments incorporated. Comments welcome

The formation of the solar system's terrestrial planets has been numerically modeled in various works, and many other studies have been devoted to characterizing our modern planets' chaotic dynamical state. However, it is still not known whether our planet's fragile chaotic state is an expected outcome of terrestrial planet accretion. We use a suite of numerical simulations to present a detailed analysis and characterization of the dynamical chaos in 145 different systems produced via terrestrial planet formation in Kaib and Cowan (2015). These systems were created in the presence of a fully formed Jupiter and Saturn, using a variety of different initial conditions. They are not meant to provide a detailed replication of the actual present solar system, but rather serve as a sample of similar systems for comparison and analysis. We find that dynamical chaos is prevalent in roughly half of the systems we form. We show that this chaos disappears in the majority of such systems when Jupiter is removed, implying that the largest source of chaos is perturbations from Jupiter. Chaos is most prevalent in systems that form 4 or 5 terrestrial planets. Additionally, an eccentric Jupiter and Saturn is shown to enhance the prevalence of chaos in systems. Furthermore, systems in our sample with a center of mass highly concentrated between 0.8 -1.2 AU generally prove to be less chaotic than systems with more exotic mass distributions. Through the process of evolving systems to the current epoch, we show that late instabilities are quite common in our systems. Of greatest interest, many of the sources of chaos observed in our own solar system (such as the secularly driven chaos between Mercury and Jupiter) are shown to be common outcomes of terrestrial planetary formation.

We analyze the evolution of the potentially habitable planet Proxima Centauri b to identify environmental factors that affect its long-term habitability. We consider physical processes acting on size scales ranging from the galactic to the stellar system to the planet's core. We find that there is a significant probability that Proxima Centauri has had encounters with its companion stars, Alpha Centauri A and B, that are close enough to destabilize an extended planetary system. If the system has an additional planet, as suggested by the discovery data, then it may perturb planet b's eccentricity and inclination, possibly driving those parameters to non-zero values, even in the presence of strong tidal damping. We also model the internal evolution of the planet, evaluating the roles of different radiogenic abundances and tidal heating and find that magnetic field generation is likely for billions of years. We find that if planet b formed in situ, then it experienced 169 +/- 13 million years in a runaway greenhouse as the star contracted during its formation. This early phase could remove up to 5 times as much water as in the modern Earth's oceans, possibly producing a large abiotic oxygen atmosphere. On the other hand, if Proxima Centauri b formed with a substantial hydrogen atmosphere (0.01 - 1% of the planet's mass), then this envelope could have shielded the water long enough for it to be retained before being blown off itself. After modeling this wide range of processes we conclude that water retention during the host star's pre-main sequence phase is the biggest obstacle for Proxima b's habitability. These results are all obtained with a new software package called VPLANET., Comment: 46 pages, 25 figures. Final version, includes corrected geophysical evolution and expanded analysis orbital/rotational evolution and atmosphere loss

We report the discovery and orbit of a new dwarf planet candidate, 2015 RR$_{245}$, by the Outer Solar System Origins Survey (OSSOS). 2015 RR$_{245}$'s orbit is eccentric ($e=0.586$), with a semi-major axis near 82 au, yielding a perihelion distance of 34 au. 2015 RR$_{245}$ has $g-r = 0.59 \pm 0.11$ and absolute magnitude $H_{r} = 3.6 \pm 0.1$; for an assumed albedo of $p_V = 12$% the object has a diameter of $\sim670$ km. Based on astrometric measurements from OSSOS and Pan-STARRS1, we find that 2015 RR$_{245}$ is securely trapped on ten-Myr timescales in the 9:2 mean-motion resonance with Neptune. It is the first TNO identified in this resonance. On hundred-Myr timescales, particles in 2015 RR$_{245}$-like orbits depart and sometimes return to the resonance, indicating that 2015 RR$_{245}$ likely forms part of the long-lived metastable population of distant TNOs that drift between resonance sticking and actively scattering via gravitational encounters with Neptune. The discovery of a 9:2 TNO stresses the role of resonances in the long-term evolution of objects in the scattering disk, and reinforces the view that distant resonances are heavily populated in the current Solar System. This object further motivates detailed modelling of the transient sticking population., Comment: Accepted to AJ

Recently, Sheppard et al. (2016) presented the discovery of seven new trans-Neptunian objects with moderate eccentricities, perihelia beyond 40 AU, and semimajor axes beyond 50 AU. Like the few previously known objects on similar orbits, these objects' semimajor axes are just beyond the Kuiper belt edge and clustered around Neptunian mean motion resonances (MMRs). These objects likely obtained their observed orbits while trapped within MMRs, when the Kozai-Lidov mechanism raised their perihelia and weakened Neptune's dynamical influence. Using numerical simulations that model the production of this population, we find that high-perihelion objects near Neptunian MMRs can constrain the nature and timescale of Neptune's past orbital migration. In particular, the population near the 3:1 MMR (near 62 AU) is especially useful due to its large population and short dynamical evolution timescale. If Neptune finishes migrating within ~100 Myrs or less, we predict over 90% of high-perihelion objects near the 3:1 MMR will have semimajor axes within 1 AU of each other, very near the modern resonance's center. On the other hand, if Neptune's migration takes ~300 Myrs, we expect ~50% of this population to reside in dynamically fossilized orbits over ~1 AU closer to the Sun than the modern resonance. We highlight 2015 KH162 as a likely member of this fossilized 3:1 population. Under any plausible migration scenario, nearly all high-perihelion objects in resonances beyond the 4:1 MMR (near 76 AU) reach their orbits well after Neptune stops migrating and comprise a recently generated, dynamically active population., Comment: Accepted to ApJ; 15 pages, 13 figures, 1 table

Cosmology experiments at mm-wavelengths can detect Planet Nine if it is the size of Neptune, has an effective temperature of 40 K, and is 700 AU from the Sun. It would appear as a ~30 mJy source at 1 mm with an annual parallax of ~5 arcmin. The challenge is to distinguish it from the approximately 4000 foreground asteroids brighter than 30 mJy. Fortunately, these asteroids are known to the Minor Planet Center and can be identified because they move across a resolution element in a matter of hours, orders of magnitude faster than Planet Nine. If Planet Nine is smaller, colder, and/or more distant than expected, then it could be as faint as 1 mJy at 1 mm. There are roughly $10^6$ asteroids this bright and many are unknown, making current cosmology experiments confusion limited for moving sources. Nonetheless, it may still be possible to find the proverbial needle in the haystack using a matched filter. This would require mm telescopes with high angular resolution and high sensitivity in order to alleviate confusion and to enable the identification of moving sources with relatively short time baselines. Regardless of its mm flux density, searching for Planet Nine would require frequent radio measurements for large swaths of the sky, including the ecliptic and Galactic plane. Even if Planet Nine had already been detected by other means, measuring its mm-flux would constrain its internal energy budget, and therefore help resolve the mystery of Uranus and Neptune, which have vastly different internal heat., Comment: 5 pages, 1 figure, ApJL accepted

We present high-contrast Magellan adaptive optics (MagAO) images of HD 7449, a Sun-like star with one planet and a long-term radial velocity (RV) trend. We unambiguously detect the source of the long-term trend from 0.6-2.15 \microns ~at a separation of \about 0\fasec 54. We use the object's colors and spectral energy distribution to show that it is most likely an M4-M5 dwarf (mass \about 0.1-0.2 \msun) at the same distance as the primary and is therefore likely bound. We also present new RVs measured with the Magellan/MIKE and PFS spectrometers and compile these with archival data from CORALIE and HARPS. We use a new Markov chain Monte Carlo procedure to constrain both the mass ($> 0.17$ \msun ~at 99$\%$ confidence) and semimajor axis (\about 18 AU) of the M dwarf companion (HD 7449B). We also refine the parameters of the known massive planet (HD 7449Ab), finding that its minimum mass is $1.09^{+0.52}_{-0.19}$ \mj, its semimajor axis is $2.33^{+0.01}_{-0.02}$ AU, and its eccentricity is $0.8^{+0.08}_{-0.06}$. We use N-body simulations to constrain the eccentricity of HD 7449B to $\lesssim$ 0.5. The M dwarf may be inducing Kozai oscillations on the planet, explaining its high eccentricity. If this is the case and its orbit was initially circular, the mass of the planet would need to be $\lesssim$ 1.5 \mj. This demonstrates that strong constraints on known planets can be made using direct observations of otherwise undetectable long-period companions., Comment: Corrected planet mass error (7.8 Mj --> 1.09 Mj, in agreement with previous studies)

We report the discovery, tracking and detection circumstances for 85 trans-Neptunian objects (TNOs) from the first 42 deg$^{2}$ of the Outer Solar System Origins Survey (OSSOS). This ongoing $r$-band Solar System survey uses the 0.9 deg$^{2}$ field-of-view MegaPrime camera on the 3.6 m Canada-France-Hawaii Telescope. Our orbital elements for these TNOs are precise to a fractional semi-major axis uncertainty $<0.1\%$. We achieve this precision in just two oppositions, as compared to the normal 3-5 oppositions, via a dense observing cadence and innovative astrometric technique. These discoveries are free of ephemeris bias, a first for large trans-Neptunian surveys. We also provide the necessary information to enable models of TNO orbital distributions to be tested against our TNO sample. We confirm the existence of a cold "kernel" of objects within the main cold classical Kuiper belt, and infer the existence of an extension of the "stirred" cold classical Kuiper belt to at least several AU beyond the 2:1 mean motion resonance with Neptune. We find that the population model of Petit et al. (2011) remains a plausible representation of the Kuiper belt. The full survey, to be completed in 2017, will provide an exquisitely characterized sample of important resonant TNO populations, ideal for testing models of giant planet migration during the early history of the Solar System., Comment: Accepted to AJ, 27 April 2016. 59 pp

Many features of the outer solar system are replicated in numerical simulations if the giant planets undergo an orbital instability that ejects one or more ice giants. During this instability, Jupiter and Saturn's orbits diverge, crossing their 2:1 mean motion resonance (MMR), and this resonance-crossing can excite the terrestrial planet orbits. Using a large ensemble of simulations of this giant planet instability, we directly model the evolution of the terrestrial planet orbits during this process, paying special attention to systems that reproduce the basic features of the outer planets. In systems that retain four giant planets and finish with Jupiter and Saturn beyond their 2:1 MMR, we find at least an 85% probability that at least one terrestrial planet is lost. Moreover, systems that manage to retain all four terrestrial planets often finish with terrestrial planet eccentricities and inclinations larger than the observed ones. There is less than a ~5% chance that the terrestrial planet orbits will have a level of excitation comparable to the observed orbits. If we factor in the probability that the outer planetary orbits are well-replicated, we find a probability of 1% or less that the orbital architectures of the inner and outer planets are simultaneously reproduced in the same system. These small probabilities raise the prospect that the giant planet instability occurred before the terrestrial planets had formed. This scenario implies that the giant planet instability is not the source of the Late Heavy Bombardment and that terrestrial planet formation finished with the giant planets in their modern configuration., Comment: Accepted to MNRAS, 10 pages, 5 figures, 2 tables

Binaries are typically excluded from direct imaging exoplanet surveys. However, the recent findings of Kepler and radial velocity programs show that planets can and do form in binary systems. Here, we suggest that visual binaries offer unique advantages for direct imaging. We show that Binary Differential Imaging (BDI), whereby two stars are imaged simultaneously at the same wavelength within the isoplanatic patch at high Strehl ratio, offers improved point spread function (PSF) subtraction that can result in increased sensitivity to planets close to each star. We demonstrate this by observing a young visual binary separated by 4\asec ~with MagAO/Clio-2 at 3.9 \microns, where the Strehl ratio is high, the isoplanatic patch is large, and giant planets are bright. Comparing BDI to angular differential imaging (ADI), we find that BDI's 5$\sigma$ contrast is \about 0.5 mags better than ADI's within \about 1\asec ~for the particular binary we observed. Because planets typically reside close to their host stars, BDI is a promising technique for discovering exoplanets in stellar systems that are often ignored. BDI is also 2-4$\times$ more efficient than ADI and classical reference PSF subtraction, since planets can be detected around both the target and PSF reference simultaneously. We are currently exploiting this technique in a new MagAO survey for giant planets in 140 young nearby visual binaries. BDI on a space-based telescope would not be limited by isoplanatism effects and would therefore be an even more powerful tool for imaging and discovering planets., Comment: Accepted to ApJ on Aug. 30, 2015. 9 pages (emulateapj), 4 figures. Full-resolution version available upon request

Kaib & Cowan (2015) recently used terrestrial planet formation simulations to conclude that the moon-forming impactor (Theia) had only a ~5% or less chance of having the same oxygen isotope composition as Earth, while Mastrobuono-Battisti et al. (2015) used seemingly similar simulations and methods to arrive at a higher value of ~20% or more. Here we derive the results of both papers from a single set of simulations. Compared to Kaib & Cowan (2015), the analysis of Mastrobuono-Battisti et al. (2015) systematically yields more massive Theia analogs and imposes flatter isotopic gradients across the original protoplanetary disk. Both of these effects diminish isotopic differences between Earth and Theia analogs. While it is notoriously difficult to produce systems resembling our actual terrestrial planets, the analysis of Kaib & Cowan (2015) more often selects and analyzes Earth and Mars analogs at orbital locations near the real planets. Given this, we conclude that the greater isotopic differences between Earth and Theia found in Kaib & Cowan (2015) better reflect the predictions of terrestrial planet formation models. Finally, although simulation uncertainties and a terrestrial contribution to Moon formation enhance the fraction of Theia analogs consistent with the canonical giant impact hypothesis, this fraction still remains in the 5-8% range., Comment: 5 pages, 4 figures, accepted to Icarus

[Abridged] We present an extensive suite of terrestrial planet formation simulations that allows quantitative analysis of the stochastic late stages of planet formation. We quantify the feeding zone width, Delta a, as the mass-weighted standard deviation of the initial semi-major axes of the planetary embryos and planetesimals that make up the final planet. The size of a planet's feeding zone in our simulations does not correlate with its final mass or semi-major axis, suggesting there is no systematic trend between a planet's mass and its volatile inventory. Instead, we find that the feeding zone of any planet more massive than 0.1M_Earth is roughly proportional to the radial extent of the initial disk from which it formed: Delta a~0.25(a_max-a_min), where a_min and a_max are the inner and outer edge of the initial planetesimal disk. These wide stochastic feeding zones have significant consequences for the origin of the Moon, since the canonical scenario predicts the Moon should be primarily composed of material from Earth's last major impactor (Theia), yet its isotopic composition is indistinguishable from Earth. In particular, we find that the feeding zones of Theia analogs are significantly more stochastic than the planetary analogs. Depending on our assumed initial distribution of oxygen isotopes within the planetesimal disk, we find a ~5% or less probability that the Earth and Theia will form with an isotopic difference equal to or smaller than the Earth and Moon's. In fact we predict that every planetary mass body should be expected to have a unique isotopic signature. In addition, we find paucities of massive Theia analogs and high velocity moon-forming collisions, two recently proposed explanations for the Moon's isotopic composition. Our work suggests that there is still no scenario for the Moon's origin that explains its isotopic composition with a high probability event., Comment: 16 pages, 22 figures, accepted for publication in Icarus; fixed typos