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.

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.

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.

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.

The migration of Neptune's resonances through the proto-Kuiper belt has been imprinted in the distribution of small bodies in the outer Solar System. Here we analyze five published Neptune migration models in detail, focusing on the high pericenter distance (high-q) trans-Neptunian Objects (TNOs) near Neptune's 5:2 and 3:1 mean-motion resonances, because they have large resonant populations, are outside the main classical belt, and are relatively isolated from other strong resonances. We compare the observationally biased output from these dynamical models with the detected TNOs from the Outer Solar System Origins Survey, via its Survey Simulator. All of the four new OSSOS detections of high-q non-resonant TNOs are on the Sunward side of the 5:2 and 3:1 resonances. We show that even after accounting for observation biases, this asymmetric distribution cannot be drawn from a uniform distribution of TNOs at 2sigma confidence. As shown by previous work, our analysis here tentatively confirms that the dynamical model that uses grainy slow Neptune migration provides the best match to the real high-q TNO orbital data. However, due to extreme observational biases, we have very few high-q TNO discoveries with which to statistically constrain the models. Thus, this analysis provides a framework for future comparison between the output from detailed, dynamically classified Neptune migration simulations and the TNO discoveries from future well-characterized surveys. We show that a deeper survey (to a limiting r-magnitude of 26.0) with a similar survey area to OSSOS could statistically distinguish between these five Neptune migration models.

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.

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.

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.

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.

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.

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.

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.

The scattering trans-Neptunian Objects (TNOs) can be measured to smaller sizes than any other distant small-body population. We use the largest sample yet obtained, 68 discoveries, primarily by the Outer Solar System Origins Survey (OSSOS), to constrain the slope of its luminosity distribution, with sensitivity to much fainter absolute $H$ magnitudes than previous work. Using the analysis technique in Shankman et al. (2016), we confirm that a single slope for the $H$-distribution is not an accurate representation of the scattering TNOs and Centaurs, and that a break in the distribution is required, in support of previous conclusions. A bright-end slope of $\alpha_b=0.9$ transitioning to a faint-end slope $\alpha_f$ of 0.4-0.5 with a differential number contrast $c$ from 1 (a knee) to 10 (a divot) provides an acceptable match to our data. We find that break magnitudes $H_b$ of 7.7 and 8.3, values both previously suggested for dynamically hot Kuiper belt populations, are equally non-rejectable for a range of $\alpha_f$ and $c$ in our statistical analysis. Our preferred divot $H$-distribution transitions to $\alpha_f=0.5$ with a divot of contrast $c=3$ at $H_b=8.3$, while our preferred knee $H$-distribution transitions to $\alpha_f=0.4$ at $H_b=7.7$. The intrinsic population of scattering TNOs required to match the OSSOS detections is $3\times10^6$ for $H_r<12$, and $9\times10^4$ for $H_r<8.66$ ($D\gtrsim100$~km), with Centaurs having an intrinsic population two orders of magnitude smaller.

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.

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.

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.

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.

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.

Comets are icy objects that orbitally evolve from the trans-Neptunian region (the Kuiper belt and beyond) into the inner Solar System, where they are heated by solar radiation and become active due to sublimation of water ice. Here we perform end-to-end simulations in which cometary reservoirs are formed in the early Solar System and evolved over 4.5 Gyr. The gravitational effects of Planet 9 (P9), hypothesized to circle the Sun on a wide orbit, are included in some of our simulations. Different models are considered for comets to be active, including a simple assumption that comets remain active for Np(q) perihelion passages with perihelion distance q<2.5 au. The orbital distribution and number of active comets produced in our model is compared to observations. The orbital distribution of ecliptic comets (ECs) is well reproduced in models with Np(2.5)=500 and without P9. With P9, the inclination distribution of model ECs is wider than the observed one. We find that the known Halley-type comets (HTCs) have a nearly isotropic inclination distribution (with only a slight preference for prograde orbits). In our model, the HTCs appear to be an extension of the population of returning Oort-cloud comets (OCCs) to shorter orbital periods. The inclination distribution of model HTCs becomes broader with increasing Np, but the existing observational data are not good enough to constrain Np from orbital fits. Np(2.5)>1000 is required to obtain a steady-state population of large active HTCs that is consistent with observations. To fit the ratio of the returning-to-new OCCs, by contrast, our model implies that Np(2.5)<10, possibly because the detected long-period comets are smaller and much easier to disrupt than observed HTCs.

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.

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.

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.

The orbital element distribution of trans-Neptunian objects (TNOs) with large pericenters has been suggested to be influenced by the presence of an undetected, large planet at >200 AU from the Sun. To find additional observables caused by this scenario, we here present the first detailed emplacement simulation in the presence of a massive ninth planet on the distant Kuiper Belt. We perform 4 Gyr N-body simulations with the currently known Solar System planetary architecture, plus a 10 Earth mass planet with similar orbital parameters to those suggested by Trujillo & Sheppard (2014) or Batygin & Brown (2016), and 10^5 test particles in an initial planetesimal disk. We find that including a distant superearth-mass planet produces a substantially different orbital distribution for the scattering and detached TNOs, raising the pericenters and inclinations of moderate semimajor axis (50<a<500 AU) objects. We test whether this signature is detectable via a simulator with the observational characteristics of four precisely characterized TNO surveys. We find that the qualitatively very distinct Solar System models that include a ninth planet are essentially observationally indistinguishable from an outer Solar System produced solely by the four giant planets. We also find that the mass of the Kuiper Belt's current scattering and detached populations is required to be 3-10 times larger in the presence of an additional planet. We do not find any evidence for clustering of orbital angles in our simulated TNO population. Wide-field, deep surveys targeting inclined high-pericenter objects will be required to distinguish between these different scenarios.

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.

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.

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.

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.

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.

We measure the absolute magnitude, $H$, distribution, $dN(H) \propto 10^{\alpha H}$ of the scattering Trans-Neptunian Objects (TNOs) as a proxy for their size-frequency distribution. We show that the H-distribution of the scattering TNOs is not consistent with a single-slope distribution, but must transition around $H_g \sim 9$ to either a knee with a shallow slope or to a divot, which is a differential drop followed by second exponential distribution. Our analysis is based on a sample of 22 scattering TNOs drawn from three different TNO surveys, the Canada-France Ecliptic Plane Survey (CFEPS, Petit et al. 2011), Alexandersen et al. (2014), and the Outer Solar System Origins Survey (OSSOS, Bannister et al. 2016), all of which provide well characterized detection thresholds, combined with a cosmogonic model for the formation of the scattering TNO population. Our measured absolute magnitude distribution result is independent of the choice of cosmogonic model. Based on our analysis, we estimate that number of scattering TNOs is (2.4-8.3)$\times 10^5$ for $H_r < 12$. A divot $H$-distribution is seen in a variety of formation scenarios and may explain several puzzles in Kuiper Belt science. We find that a divot $H$-distribution simultaneously explains the observed scattering TNO, Neptune Trojan, Plutino, and Centaur $H$-distributions while simultaneously predicting a large enough scattering TNO population to act as the sole supply of the Jupiter-Family Comets.

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.

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.

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.

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

We present Rossiter-McLaughlin observations of the transiting super-Earth 55 Cnc e collected during six transit events between January 2012 and November 2013 with HARPS and HARPS-N. We detect no radial-velocity signal above 35 cm/s (3-sigma) and confine the stellar v sin i to 0.2 +/- 0.5 km/s. The star appears to be a very slow rotator, producing a very low amplitude Rossiter-McLaughlin effect. Given such a low amplitude, the Rossiter-McLaughlin effect of 55 Cnc e is undetected in our data, and any spin-orbit angle of the system remains possible. We also performed Doppler tomography and reach a similar conclusion. Our results offer a glimpse of the capacity of future instrumentation to study low amplitude Rossiter-McLaughlin effects produced by super-Earths.

We present numerical simulations modeling the orbital evolution of very wide binaries, pairs of stars separated by over ~1000 AU. Due to perturbations from other passing stars and the Milky Way's tide, the orbits of very wide binary stars occasionally become extremely eccentric, which forces close encounters between the companion stars (Kaib et al. 2013). We show that this process causes a stellar collision between very wide binary companion stars once every 1000-7500 years on average in the Milky Way. One of the main uncertainties in this collision rate is the amount of energy dissipated by dynamic tides during close (but not collisional) periastron passages. This dissipation presents a dynamical barrier to stellar collisions and can instead transform very wide binaries into close or contact binaries. However, for any plausible tidal dissipation model, very wide binary stars are an unrealized, and potentially the dominant, source of stellar collisions in our Galaxy. Such collisions should occur throughout the thin disk of the Milky Way. Stellar collisions within very wide binaries should yield a small population of single, Li-depleted, rapidly rotating massive stars.

We study the influence of outer Solar System architecture on the structural evolution of the Oort Cloud (OC) and the flux of Earth-crossing comets. In particular, we seek to quantify the role of the giant planets as "planetary protectors". To do so, we have run simulations in each of four different planetary mass configurations to understand the significance of each of the giant planets. Because the outer planets modify the structure of the OC throughout its formation, we integrate each simulation over the full age of the Solar System. Over this time, we follow the evolution of cometary orbits from their starting point in the protoplanetary disk to their injection into the OC to their possible re-entry into the inner planetary region. We find that the overall structure of the OC, including the location of boundaries and the relative number of comets in the inner and outer parts, does not change significantly between configurations; however, as planetary mass decreases, the trapping efficiency (TE) of comets into the OC and the flux of comets into the observable region increases. We determine that those comets that evolve onto Earth-crossing orbits come primarily from the inner OC but show no preference for initial protoplanetary disk location. We also find that systems that have at least a Saturn-mass object are effective at deflecting possible Earth-crossing comets but the difference in flux between systems with and without such a planet is less than an order of magnitude. We conclude by discussing the individual roles of the planets and the implications of incorporating more realistic planetary accretion and migration scenarios into simulations, particularly on existing discrepancies between low TE and the mass of the protoplanetary disk and on determining the structural boundaries of the OC.

We report the discovery and confirmation of a transiting circumbinary planet (PH1b) around KIC 4862625, an eclipsing binary in the Kepler field. The planet was discovered by volunteers searching the first six Quarters of publicly available Kepler data as part of the Planet Hunters citizen science project. Transits of the planet across the larger and brighter of the eclipsing stars are detectable by visual inspection every ~137 days, with seven transits identified in Quarters 1-11. The physical and orbital parameters of both the host stars and planet were obtained via a photometric-dynamical model, simultaneously fitting both the measured radial velocities and the Kepler light curve of KIC 4862625. The 6.18 +/- 0.17 Earth radii planet orbits outside the 20-day orbit of an eclipsing binary consisting of an F dwarf (1.734 +/- 0.044 Solar radii, 1.528 +/- 0.087 Solar masses) and M dwarf (0.378+/- 0.023 Solar radii, 0.408 +/- 0.024 Solar masses). For the planet, we find an upper mass limit of 169 Earth masses (0.531 Jupiter masses) at the 99.7% confidence level. With a radius and mass less than that of Jupiter, PH1b is well within the planetary regime. Outside the planet's orbit, at ~1000 AU,a previously unknown visual binary has been identified that is likely bound to the planetary system, making this the first known case of a quadruple star system with a transiting planet.

Via joint analysis of a calibrated telescopic survey, which found scattering Kuiper Belt objects, and models of their expected orbital distribution, we measure the form of the scattering object's size distribution. Ruling out a single power-law at greater than 99% confidence, we constrain the form of the size distribution and find that, surprisingly, our analysis favours a very sudden decrease (a divot) in the number distribution as diameters decrease below 100 km, with the number of smaller objects then rising again as expected via collisional equilibrium. Extrapolating at this collisional equilibrium slope produced enough kilometer-scale scattering objects to supply the nearby Jupiter-Family comets. Our interpretation is that this divot feature is a preserved relic of the size distribution made by planetesimal formation, now "frozen in" to portions of the Kuiper Belt sharing a "hot" orbital inclination distribution, explaining several puzzles in Kuiper Belt science. Additionally, we show that to match today's scattering-object inclination distribution, the supply source that was scattered outward must have already been vertically heated to of order 10 degrees.

Nearly half of the exoplanets found within binary star systems reside in very wide binaries with average stellar separations beyond 1,000 AU (1 AU being the Earth-Sun distance), yet the influence of such distant binary companions on planetary evolution remains largely unstudied. Unlike their tighter counterparts, the stellar orbits of wide binaries continually change under the influence of the Galactic tide and impulses from other passing stars. Here we report numerical simulations demonstrating that the variable nature of wide binary star orbits dramatically reshapes the planetary systems they host, typically Gyrs after formation. Contrary to previous understanding, wide binary companions may often strongly perturb planetary systems, triggering planetary ejections and exciting orbital eccentricities of surviving planets. Indeed, observed exoplanet eccentricities offer evidence of this; giant exoplanet orbits within wide binaries are statistically more eccentric than those around isolated stars. Both eccentricity distributions are well-reproduced when we assume isolated stars and wide binaries host similar planetary systems whose outermost giant planets are scattered beyond ~10 AU from their parent stars via early internal instabilities. Consequently, our results suggest that although wide binaries eventually truncate their planetary systems, most isolated giant exoplanet systems harbor additional distant, still undetected planets.

Although the 55 Cnc system contains multiple, closely packed planets that are presumably in a coplanar configuration, we use numerical simulations to demonstrate that they are likely to be highly inclined to their parent star's spin axis. Due to perturbations from its distant binary companion, this planetary system precesses like a rigid body about its parent star. Consequently, the parent star's spin axis and the planetary orbit normal likely diverged long ago. Because only the projected separation of the binary is known, we study this effect statistically, assuming an isotropic distribution for wide binary orbits. We find that the most likely projected spin-orbit angle is ~50 degrees, with a ~30% chance of a retrograde configuration. Transit observations of the innermost planet - 55 Cnc e - may be used to verify these findings via the Rossiter-McLaughlin effect. 55 Cancri may thus represent a new class of planetary systems with well-ordered, coplanar orbits that are inclined with respect to the stellar equator.

Recent numerical simulations have demonstrated that the Sun's dynamical history within the Milky Way may be much more complex than that suggested by its current low peculiar velocity. In particular, the Sun may have radially migrated through the galactic disk by up to 5-6 kpc. This has important ramifications for the structure of the Oort Cloud, as it means that the solar system may have experienced tidal and stellar perturbations that were significantly different from its current local galactic environment. To characterize the effects of solar migration within the Milky Way, we use direct numerical simulations to model the formation of an Oort Cloud around stars that end up on solar-type orbits in a galactic-scale simulation of a Milky Way-like disk formation. Surprisingly, our simulations indicate that Sedna's orbit may belong to the classical Oort Cloud. Contrary to previous understanding, we show that field star encounters play a pivotal role in setting the Oort Cloud's extreme inner edge, and due to their stochastic nature this inner edge sometimes extends to Sedna's orbit. The Sun's galactic migration heightens the chance of powerful stellar passages, and Sedna production occurs around ~20-30% of the solar-like stars we study. Considering the entire Oort Cloud, we find its median distance depends on the minimum galactocentric distance attained during the Sun's orbital history. The inner edge also shows a similar dependence but with increased scatter due to the effects of powerful stellar encounters. Both of these Oort Cloud parameters can vary by an order of magnitude and are usually overestimated by an Oort Cloud formation model that assumes a fixed galactic environment. In addition, the amount of material trapped in outer Oort Cloud orbits can be extremely low and may present difficulties for models of Oort Cloud formation and long-period comet production.

It has previously been shown that varying the numerical timestep during a symplectic orbital integration leads to a random walk in energy and angular momentum, destroying the phase space-conserving property of symplectic integrators. Here we show that when altering the timestep symplectic correctors can be used to reduce this error to a negligible level. Furthermore, these correctors can also be employed to avoid a large error introduction when changing the Hamiltonian's partitioning. We have constructed a numerical integrator using this technique that is nearly as accurate as widely used fixed-step routines. In addition, our algorithm is drastically faster for integrations of highly eccentricitic, large semimajor axis orbits, such as those found in the Oort Cloud.

We present numerical simulations to model the production of observable long-period comets (LPCs) from the Oort Cloud, a vast reservoir of icy bodies surrounding the Sun. We show that inner Oort Cloud objects can penetrate Jupiter's orbit via a largely unexplored dynamical pathway, and they are an important, if not the dominant, source of known LPCs. We use this LPC production to place observationally motivated constraints on the population and mass of the inner Oort Cloud, which are consistent with giant planet formation theory. These constraints indicate that only one comet shower producing late Eocene bombardment levels has likely occurred since the Cambrian Explosion, making these phenomena an improbable cause of additional extinction events.

CoRoT-7 b, a planet about 70% larger than the Earth orbiting a Sun-like star, is the first-discovered rocky exoplanet, and hence has been dubbed a "super-Earth". Some initial studies suggested that since the planet is so close to its host star, it receives enough insolation to partially melt its surface. However, these past studies failed to take into consideration the role that tides may play in this system. Even if the planet's eccentricity has always been zero, we show that tidal decay of semi-major axis could have been large enough that the planet formed on a wider orbit which received less insolation. Moreover, CoRoT-7 b could be tidally heated at a rate that dominates its geophysics and drives extreme volcanism. In this case, CoRoT-7 b is a "super-Io" that, like Jupiter's volcanic moon, is dominated by volcanism and rapid resurfacing. Such heating could occur with an eccentricity of just 10^-5. This small value could be driven by CoRoT-7 c if its own eccentricity is larger than ~10^-4. CoRoT-7 b may be the first of a class of planetary super-Ios likely to be revealed by the CoRoT and Kepler spacecraft.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

To date, no accretion model has succeeded in reproducing all observed constraints in the inner Solar System. These constraints include 1) the orbits, in particular the small eccentricities, and 2) the masses of the terrestrial planets -- Mars' relatively small mass in particular has not been adequately reproduced in previous simulations; 3) the formation timescales of Earth and Mars, as interpreted from Hf/W isotopes; 4) the bulk structure of the asteroid belt, in particular the lack of an imprint of planetary embryo-sized objects; and 5) Earth's relatively large water content, assuming that it was delivered in the form of water-rich primitive asteroidal material. Here we present results of 40 high-resolution (N=1000-2000) dynamical simulations of late-stage planetary accretion with the goal of reproducing these constraints, although neglecting the planet Mercury. We assume that Jupiter and Saturn are fully-formed at the start of each simulation, and test orbital configurations that are both consistent with and contrary to the "Nice model." We find that a configuration with Jupiter and Saturn on circular orbits forms low-eccentricity terrestrial planets and a water-rich Earth on the correct timescale, but Mars' mass is too large by a factor of 5-10 and embryos are often stranded in the asteroid belt. A configuration with Jupiter and Saturn in their current locations but with slightly higher initial eccentricities (e = 0.07-0.1) produces a small Mars, an embryo-free asteroid belt, and a reasonable Earth analog but rarely allows water delivery to Earth. None of the configurations we tested reproduced all the observed constraints. (abridged)

We report the discovery of a minor planet (2006 SQ372) on an orbit with a perihelion of 24 AU and a semimajor axis of 796 AU. Dynamical simulations show that this is a transient orbit and is unstable on a timescale of 200 Myrs. Falling near the upper semimajor axis range of the scattered disk and the lower semimajor axis range of the Oort Cloud, previous membership in either class is possible. By modeling the production of similar orbits from the Oort Cloud as well as from the scattered disk, we find that the Oort Cloud produces 16 times as many objects on SQ372-like orbits as the scattered disk. Given this result, we believe this to be the most distant long-period comet ever discovered. Furthermore, our simulation results also indicate that 2000 OO67 has had a similar dynamical history. Unaffected by the "Jupiter-Saturn Barrier," these two objects are most likely long-period comets from the inner Oort Cloud.

We report the discovery and orbit determination of 14 trans-Neptunian objects (TNOs) from the ESSENCE Supernova Survey difference imaging dataset. Two additional objects discovered in a similar search of the SDSS-II Supernova Survey database were recovered in this effort. ESSENCE repeatedly observed fields far from the Solar System ecliptic (-21 deg < beta < -5 deg), reaching limiting magnitudes per observation of I~23.1 and R~23.7. We examine several of the newly detected objects in detail, including 2003 UC_414 which orbits entirely between Uranus and Neptune and lies very close to a dynamical region that would make it stable for the lifetime of the Solar System. 2003 SS_422 and 2007 TA_418 have high eccentricities and large perihelia, making them candidate members of an outer class of trans-Neptunian objects. We also report a new member of the ''extended'' or ''detached'' scattered disk, 2004 VN_112, and verify the stability of its orbit using numerical simulations. This object would have been visible to ESSENCE for only ~2% of its orbit, suggesting a vast number of similar objects across the sky. We emphasize that off-ecliptic surveys are optimal for uncovering the diversity of such objects, which in turn will constrain the history of gravitational influences that shaped our early Solar System.

We study the influence of an open cluster environment on the formation and current structure of the Oort cloud. To do this, we have run 19 different simulations of the formation of the Oort Cloud for 4.5 Gyrs. In each simulation, the solar system spends its first 100 Myrs in a different open cluster environment before transitioning to its current field environment. We find that, compared to forming in the field environment, the inner Oort Cloud is preferentially loaded with comets while the Sun resides in the open cluster and that most of this material remains locked in the interior of the cloud for the next 4.4 Gyrs. In addition, the outer Oort Cloud trapping efficiencies we observe in our simulations are lower than previous formation models by about a factor of 2, possibly implying an even more massive early planetesimal disk. Furthermore, some of our simulations reproduce the orbits of observed extended scattered disk objects, which may serve as an observational constraint on the Sun's early environment. Depending on the particular open cluster environment, the properties of the inner Oort Cloud and extended scattered disk can vary widely. On the other hand, the outer portions of the Oort Cloud in each of our simulations are all similar.

The Sloan Digital Sky Survey (SDSS) and Two Micron All Sky Survey (2MASS) are rich resources for studying stellar astrophysics and the structure and formation history of the Galaxy. As new surveys and instruments adopt similar filter sets, it is increasingly important to understand the properties of the ugrizJHKs stellar locus, both to inform studies of `normal' main sequence stars as well as for robust searches for point sources with unusual colors. Using a sample of ~600,000 point sources detected by SDSS and 2MASS, we tabulate the position and width of the ugrizJHKs stellar locus as a function of g-i color, and provide accurate polynomial fits. We map the Morgan-Keenan spectral type sequence to the median stellar locus by using synthetic photometry of spectral standards and by analyzing 3000 SDSS stellar spectra with a custom spectral typing pipeline. We develop an algorithm to calculate a point source's minimum separation from the stellar locus in a seven-dimensional color space, and use it to robustly identify objects with unusual colors, as well as spurious SDSS/2MASS matches. Analysis of a final catalog of 2117 color outliers identifies 370 white-dwarf/M dwarf (WDMD) pairs, 93 QSOs, and 90 M giant/carbon star candidates, and demonstrates that WDMD pairs and QSOs can be distinguished on the basis of their J-Ks and r-z colors. We also identify a group of objects with correlated offsets in the u-g vs. g-r and g-r vs. r-i color-color spaces, but subsequent follow-up is required to reveal the nature of these objects. Future applications of this algorithm to a matched SDSS-UKIDSS catalog may well identify additional classes of objects with unusual colors by probing new areas of color-magnitude space.

Recent results have shown that many of the known extrasolar planetary systems contain regions which are stable for both Earth-mass and Saturn-mass planets. Here we simulate the formation of terrestrial planets in four planetary systems -- 55 Cancri, HD 38529, HD 37124, and HD 74156 -- under the assumption that these systems of giant planets are complete and that their orbits are well-determined. Assuming the giant planets formed and migrated quickly, then terrestrial planets may form from a second generation of planetesimals. In each case, Moon- to Mars-sized planetary embryos are placed in between the giant planets and evolved for 100 Myr. We find that planets form relatively easily in 55 Cnc, with masses up to 0.6 Earth masses and in some cases substantial water contents and orbits in the habitable zone. HD 38529 is likely to support an asteroid belt but no terrestrial planets of significant mass. No terrestrial planets form in HD 37124 and HD 74156, although in some cases 1-2 lone embryos survive for 100 Myr. If migration occurred later, depleting the planetesimal disk, then massive terrestrial planets are unlikely to form in any of these systems.

We discuss the UV, optical, and IR properties of the SDSS sources detected by GALEX as part of its All-sky Imaging Survey Early Release Observations. Virtually all of the GALEX sources in the overlap region are detected by SDSS. GALEX sources represent ~2.5% of all SDSS sources within these fields and about half are optically unresolved. Most unresolved GALEX/SDSS sources are bright blue turn-off thick disk stars and are typically detected only in the GALEX near-UV band. The remaining unresolved sources include low-redshift quasars, white dwarfs, and white dwarf/M dwarf pairs, and these dominate the optically unresolved sources detected in both GALEX bands. Almost all the resolved SDSS sources detected by GALEX are fainter than the SDSS 'main' spectroscopic limit. These sources have colors consistent with those of blue (spiral) galaxies (u-r<2.2), and most are detected in both GALEX bands. Measurements of their UV colors allow much more accurate and robust estimates of star-formation history than are possible using only SDSS data. Indeed, galaxies with the most recent (<20 Myr) star formation can be robustly selected from the GALEX data by requiring that they be brighter in the far-UV than in the near-UV band. However, older starburst galaxies have UV colors similar to AGN, and thus cannot be selected unambiguously on the basis of GALEX fluxes alone. With the aid of 2MASS data, we construct and discuss median 10 band UV-optical-IR spectral energy distributions for turn-off stars, hot white dwarfs, low-redshift quasars, and spiral and elliptical galaxies. We point out the high degree of correlation between the UV color and the contribution of the UV flux to the UV-optical-IR flux of galaxies detected by GALEX.

We report the search for intracluster light in four Abell Type II/III (non-cD) galaxy clusters: Abell 801, 1234, 1553, & 1914. We find on average that these clusters contain $\sim$ 10% of their detected stellar luminosity in a diffuse component. We show that for two of the clusters the intracluster light closely follows the galaxy distribution, but in the other two cases, there are noticeable differences between the spatial distribution of the galaxies and the intracluster light. We report the results of a search for intracluster tidal debris in each cluster, and note that Abell 1914 in particular has a number of strong tidal features likely due to its status as a recent cluster merger. One of the Abell 1914 features appears to be spatially coincident with an extension seen in weak lensing maps, implying the feature traces a large amount of mass. We compare these results to numerical simulations of hierarchically-formed galaxy clusters, and find good general agreement between the observed and simulated images, although we also find that our observations sample only the brightest features of the intracluster light. Together, these results suggest that intracluster light can be a valuable tool in determining the evolutionary state of galaxy clusters.

We present some early results from our deep imaging survey of galaxy clusters intended to detect and study intracluster light (ICL). From our observations to date, we find that ICL is common in galaxy clusters, and that substructure in the ICL also appears to be common as well. We also discuss some initial comparisons of our imaging results to high-resolution numerical simulations of galaxy clusters, and give avenues for future research.