We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system.

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.

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.

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.

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 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 have investigated the effects of uniform rotation and a specific model for differential rotation on the pulsation frequencies of 10 \Msun\ stellar models. Uniform rotation decreases the frequencies for all modes. Differential rotation does not appear to have a significant effect on the frequencies, except for the most extreme differentially rotating models. In all cases, the large and small separations show the effects of rotation at lower velocities than do the individual frequencies. Unfortunately, to a certain extent, differential rotation mimics the effects o f more rapid rotation, and only the presence of some specific observed frequencies with well identified modes will be able to uniquely constrain the internal rotation of pulsating stars.