In two recent papers published in MNRAS, Namouni and Morais (2018, 2020) claimed evidence for the interstellar origin of some small Solar System bodies, including i) objects in retrograde co-orbital motion with the giant planets, and ii) the highly-inclined Centaurs. Here, we discuss the flaws of those papers that invalidate the authors' conclusions. Numerical simulations backwards in time are not representative of the past evolution of real bodies. Instead, these simulations are only useful as a means to quantify the short dynamical lifetime of the considered bodies and the fast decay of their population. In light of this fast decay, if the observed bodies were the survivors of populations of objects captured from interstellar space in the early Solar System, these populations should have been implausibly large (e.g. about 10 times the current main asteroid belt population for the retrograde coorbital of Jupiter). More likely, the observed objects are just transient members of a population that is maintained in quasi-steady state by a continuous flux of objects from some parent reservoir in the distant Solar System. We identify in the Halley type comets and the Oort cloud the most likely sources of retrograde coorbitals and highly-inclined Centaurs.

While the Solar System contains no planets between the sizes of Uranus and Saturn, our current exoplanet census includes several dozen such planets with well-measured masses and radii. These sub-Saturns exhibit a diversity of bulk densities, ranging from ~$0.1-3\ \rm{g\ cm}^{-3}$. When modeled simply as hydrogen/helium envelopes atop rocky cores, this diversity in densities translates to a diversity in planetary envelope fractions, $f_\rm{env}=M_\rm{env}/M_p$ ranging from ~$10\%$ to ~$50\%$. Planets with $f_\rm{env}\sim50\%$ pose a challenge to traditional models of giant planet formation by core-nucleated accretion, which predict the onset of runaway gas accretion when $M_\rm{env}\sim M_\rm{core}$. Here we show that many of these apparent $f_\rm{env}\sim50\%$ planets are less envelope rich than they seem, after accounting for tidal heating. We present a new framework for modeling sub-Saturn interiors that incorporates envelope inflation due to tides, which are driven by the observed non-zero eccentricities, as well as potential obliquities. Consequently, when we apply our models to known sub-Saturns, we infer lower $f_\rm{env}$ than tides-free estimates. We present a case study of K2-19 b, a moderately eccentric sub-Saturn. Neglecting tides, K2-19 b appears to have $f_\rm{env}\sim50\%$, poised precariously near the runaway threshold; by including tides, we find $f_\rm{env}\sim10\%$, resolving the tension. Through a systematic analysis of $4-8\ R_{\oplus}$ planets, we find that most (but not all) of the similarly envelope-rich planets have more modest envelopes of $f_\rm{env}\sim10\%-20\%$. Thus, many sub-Saturns may be understood as sub-Neptunes that have undergone significant radius inflation, rather than a separate class of objects. Tidal radius inflation likely plays an important role in other size classes of planets including ultra-low-density Jupiter-size planets like WASP-107 b.

Recent analyses have shown that the concluding stages of giant planet formation are accompanied by the development of large-scale meridional flow of gas inside the planetary Hill sphere. This circulation feeds a circumplanetary disk that viscously expels gaseous material back into the parent nebula, maintaining the system in a quasi-steady state. Here we investigate the formation of natural satellites of Jupiter and Saturn within the framework of this newly outlined picture. We begin by considering the long-term evolution of solid material, and demonstrate that the circumplanetary disk can act as a global dust trap, where $s_{\bullet}\sim0.1-10\,$mm grains achieve a hydrodynamical equilibrium, facilitated by a balance between radial updraft and aerodynamic drag. This process leads to a gradual increase in the system's metallicity, and eventually culminates in the gravitational fragmentation of the outer regions of the solid sub-disk into $\mathcal{R}\sim100\,$km satellitesimals. Subsequently, satellite conglomeration ensues via pairwise collisions, but is terminated when disk-driven orbital migration removes the growing objects from the satellitesimal feeding zone. The resulting satellite formation cycle can repeat multiple times, until it is brought to an end by photo-evaporation of the parent nebula. Numerical simulations of the envisioned formation scenario yield satisfactory agreement between our model and the known properties of the Jovian and Saturnian moons.

Rotationally-modulated variability of brown dwarfs and giant planets provides unique information about their surface brightness inhomogeneities, atmospheric circulation, cloud evolution, vertical atmospheric structure, and rotational angular momentum. We report results from Hubble Space Telescope/Wide Field Camera 3 near-infrared time-series spectroscopic observations of three companions with masses in or near the planetary regime: VHS J125601.92-125723.9 b, GSC 6214-210 B, and ROXs 42 B b. VHS J1256-1257 b exhibits strong total intensity and spectral variability with a brightness difference of 19.3% between 1.1-1.7 $\mu$m over 8.5 hours and even higher variability at the 24.7% level at 1.27 $\mu$m. The light curve of VHS J1256-1257 b continues to rise at the end of the observing sequence so these values represent lower limits on the full variability amplitude at this epoch. This observed variability rivals (and may surpass) the most variable brown dwarf currently known, 2MASS J21392676+0220226. The implied rotation period of VHS J1256-1257 b is $\approx$21-24 hr assuming sinusoidal modulations, which is unusually long for substellar objects. No significant variability is evident in the light curves of GSC 6214-210 B ($<$1.2%) and ROXs 42 B b ($<$15.6%). With a spectral type of L7, an especially red spectrum, and a young age, VHS J1256-1257 b reinforces emerging patterns between high variability amplitude, low surface gravity, and evolutionary phase near the L/T transition.

Observations indicate that members of close stellar binaries often have mass ratios close to unity, while longer-period systems exhibit a more uniform mass-ratio distribution. This paper provides a theoretical explanation for this finding by determining the tidal equilibrium states for binary star systems --- subject to the constraints of conservation of angular momentum and constant total mass. This work generalizes previous treatments by including the mass fraction as a variable in the optimization problem. The results show that the lowest energy state accessible to the system corresponds to equal mass stars on a circular orbit, where the stellar spin angular velocities are both synchronized and aligned with the orbit. These features are roughly consistent with observed properties of close binary systems. We also find the conditions required for this minimum energy state to exist: [1] The total angular momentum must exceed a critical value, [2] the orbital angular momentum must be three times greater than the total spin angular momentum, and [3] the semimajor axis is bounded from above. The last condition implies that sufficiently wide binaries are not optimized with equal mass stars, where the limiting binary separation occurs near $a_0\approx16R_\ast$.

Motivated by the trends found in the observed sample of extrasolar planets, this paper determines tidal equilibrium states for forming planetary systems --- subject to conservation of angular momentum, constant total mass, and fixed orbital spacing. In the low-mass limit, valid for superearth-class planets with masses of order $m_{\rm p}\sim10M_\oplus$, previous work showed that energy optimization leads to nearly equal mass planets, with circular orbits confined to a plane. The present treatment generalizes previous results by including the self-gravity of the planetary bodies. For systems with sufficiently large total mass $m_{\scriptstyle T}$ in planets, the optimized energy state switches over from the case of nearly equal mass planets to a configuration where one planet contains most of the material. This transition occurs for a critical mass threshold of approximately $m_{\scriptstyle T} \ge m_{\scriptstyle C}\sim40M_\oplus$ (where the value depends on the semimajor axes of the planetary orbits, the stellar mass, and other system properties). These considerations of energy optimization apply over a wide range of mass scales, from binary stars to planetary systems to the collection of moons orbiting the giant planets in our solar system.

Most planetary systems -- including our own -- are born within stellar clusters, where interactions with neighboring stars can help shape the system architecture. This paper develops an orbit-averaged formalism to characterize the cluster's mean-field effects as well as the physics of long-period stellar encounters. Our secular approach allows for an analytic description of the dynamical consequences of the cluster environment on its constituent planetary systems. We analyze special cases of the resulting Hamiltonian, corresponding to eccentricity evolution driven by planar encounters, as well as hyperbolic perturbations upon dissipative disks. We subsequently apply our results to the early evolution of our solar system, where the cluster's collective potential perturbs the solar system's plane, and stellar encounters act to increase the velocity dispersion of the Kuiper belt. Our results are two-fold: first, we find that cluster effects can alter the mean plane of the solar system by $\lesssim1\deg$, and are thus insufficient to explain the $\psi\approx6\deg$ obliquity of the sun. Second, we delineate the extent to which stellar flybys excite the orbital dispersion of the cold classical Kuiper belt, and show that while stellar flybys may grow the cold belt's inclination by the observed amount, the resulting distribution is incompatible with the data. Correspondingly, our calculations place an upper limit on the product of the stellar number density and residence time of the sun in its birth cluster, $\eta\,\tau\lesssim2\times10^4\,$Myr/pc$^3$.

K2-19b and c were among the first planets discovered by NASA's K2 mission and together stand in stark contrast with the physical and orbital properties of the solar system planets. The planets are between the size of Uranus and Saturn at 7.0$\pm$0.2 R_E and 4.1$\pm$0.2 R_E, respectively, and reside a mere 0.1% outside the nominal 3:2 mean-motion resonance. They represent a different outcome of the planet formation process than the solar system, as well as the vast majority of known exoplanets. We measured the physical and orbital properties of these planets using photometry from K2, Spitzer, and ground-based telescopes, along with radial velocities from Keck/HIRES. Through a joint photodynamical model, we found that the planets have moderate eccentricities of $e \approx0.20$ and well-aligned apsides $\Delta \varpi \approx 0$ deg. The planets occupy a strictly non-resonant configuration: the resonant angles circulate rather than librate. This defies the predictions of standard formation pathways that invoke convergent or divergent migration, both of which predict $\Delta \varpi \approx 180$ deg and eccentricities of a few percent or less. We measured masses of $M_{p,b}$ = 32.4$\pm$1.7 M_E and $M_{p,c}$ = 10.8$\pm$0.6 M_E. Our measurements, with 5% fractional uncertainties, are among the most precise of any sub-Jovian exoplanet. Mass and size reflect a planet's core/envelope structure. Despite having a relatively massive core of $M_{core} \approx15$ $M_E$, K2-19b is envelope-rich, with an envelope mass fraction of roughly 50%. This planet poses a challenge to standard models core-nucleated accretion, which predict that cores $\gtrsim 10$ $M_E$ will quickly accrete gas and trigger runaway accretion when the envelope mass exceeds that of the core.

The process of planet conglomeration, which primarily unfolds in a geometrically thin disk of gas and dust, is often accompanied by dynamical excitation of the forming planets and planetesimals. The ensuing orbital crossing can lead to large-scale collisional fragmentation, populating the system with icy and rocky debris. In a gaseous nebula, such leftover solid matter tends to spiral down towards the host star due to aerodynamic drag. Along the way, the inward drifting debris can encounter planets and gravitationally couple to them via mean-motion resonances, sapping them of their orbital energy and causing them to migrate. Here, we develop a simple theory for this migration mechanism, which we call "Aero-Resonant Migration" (ARM), in which small planetesimals ($10\,$m $\lesssim s \lesssim 10\,$km) undergo orbital decay due to aerodynamic drag and resonantly shepherd planets ahead of them. Using a combination of analytical calculations and numerical experiments, we show that ARM is a robust migration mechanism, able to significantly transport planets on timescales $\lesssim1$ Myr, and present simple formulae for the ARM rate.

The tilt of a planet's spin axis off its orbital axis ("obliquity") is a basic physical characteristic that plays a central role in determining the planet's global circulation and energy redistribution. Moreover, recent studies have also highlighted the importance of obliquities in sculpting not only the physical features of exoplanets but also their orbital architectures. It is therefore of key importance to identify and characterize the dominant processes of excitation of non-zero axial tilts. Here we highlight a simple mechanism that operates early on and is likely fundamental for many extrasolar planets and perhaps even Solar System planets. While planets are still forming in the protoplanetary disk, the gravitational potential of the disk induces nodal recession of the orbits. The frequency of this recession decreases as the disk dissipates, and when it crosses the frequency of a planet's spin axis precession, large planetary obliquities may be excited through capture into a secular spin-orbit resonance. We study the conditions for encountering this resonance and calculate the resulting obliquity excitation over a wide range of parameter space. Planets with semi-major axes in the range $0.3 \ \mathrm{AU} \lesssim a \lesssim 2 \ \mathrm{AU}$ are the most readily affected, but large-$a$ planets can also be impacted. We present a case study of Uranus and Neptune and show that this mechanism likely cannot help explain their high obliquities. While it could have played a role if finely tuned and envisioned to operate in isolation, large-scale obliquity excitation was likely inhibited by gravitational planet-planet perturbations.

Early dynamical evolution of close-in planetary systems is shaped by an intricate combination of planetary gravitational interactions, orbital migration, and dissipative effects. While the process of convergent orbital migration is expected to routinely yield resonant planetary systems, previous analyses have shown that the semi-major axes of initially resonant pairs of planets will gradually diverge under the influence of long-term energy damping, producing an overabundance of planetary period ratios in slight excess of exact commensurability. While this feature is clearly evident in the orbital distribution of close-in extrasolar planets, the existing theoretical picture is limited to the specific case of the planetary three-body problem. In this study, we generalise the framework of dissipative divergence of resonant orbits to multi-resonant chains, and apply our results to the current observational census of well-characterised three-planet systems. Focusing on the 2:1 and 3:2 commensurabilities, we identify three 3-planet systems, whose current orbital architecture is consistent with an evolutionary history wherein convergent migration first locks the planets into a multi-resonant configuration and subsequent dissipation repels the orbits away from exact commensurability. Nevertheless, we find that the architecture of the overall sample of multi planetary systems is incompatible with this simple scenario, suggesting that additional physical mechanisms must play a dominant role during the early stages of planetary systems' dynamical evolution.

We show that the $P\sim8\,{\rm h}$ photometric period and the astrometrically measured $A_{\rm ng}\sim2.5\times10^{-4}\,{\rm cm\,s^{-2}}$ non-gravitational acceleration (at $r\sim1.4\,{\rm AU}$) of the interstellar object 1I/2017 (`Oumuamua) can be explained by a nozzle-like venting of volatiles whose activity migrated to track the sub-solar location on the object's surface. Adopting the assumption that `Oumuamua was an elongated $a\times b \times c$ ellipsoid, this model produces a pendulum-like rotation of the body and implies a long semi-axis $a\sim 5A_{\rm ng}P^2/4\pi^2 \sim 260\,{\rm m}$. This scale agrees with the independent estimates of `Oumuamua's size that stem from its measured brightness, assuming an albedo of $p\sim0.1$, appropriate to ices that have undergone long-duration exposure to the interstellar cosmic ray flux. Using ray-tracing, we generate light curves for ellipsoidal bodies that are subject to both physically consistent sub-solar torques and to the time-varying geometry of the Sun-Earth-`Oumuamua configuration. Our synthetic light curves display variations from chaotic tumbling and changing cross-sectional illumination that are consistent with the observations, while avoiding significant secular changes in the photometric periodicity. If our model is correct, `Oumuamua experienced mass loss that wasted $\sim 10\%$ of its total mass during the $\sim 100\,{\rm d}$ span of its encounter with the inner Solar System and had an icy composition with a very low $[{\rm C}/{\rm O}]\lesssim 0.003$. Our interpretation of `Oumuamua's behavior is consistent with the hypothesis that it was ejected from either the outer regions of a planetesimal disk after an encounter with an embedded $M_{\rm p} \sim M_{\rm Nep}$ planet or from an exo-Oort cloud.

Over the course of the past two decades, observational surveys have unveiled the intricate orbital structure of the Kuiper Belt, a field of icy bodies orbiting the Sun beyond Neptune. In addition to a host of readily-predictable orbital behavior, the emerging census of trans-Neptunian objects displays dynamical phenomena that cannot be accounted for by interactions with the known eight-planet solar system alone. Specifically, explanations for the observed physical clustering of orbits with semi-major axes in excess of $\sim250\,$AU, the detachment of perihelia of select Kuiper belt objects from Neptune, as well as the dynamical origin of highly inclined/retrograde long-period orbits remain elusive within the context of the classical view of the solar system. This newly outlined dynamical architecture of the distant solar system points to the existence of a new planet with mass of $m_9\sim 5-10\,M_{\oplus}$, residing on a moderately inclined orbit ($i_9\sim15-25\deg$) with semi-major axis $a_9\sim 400 - 800\,$AU and eccentricity between $e_9 \sim 0.2 - 0.5$. This paper reviews the observational motivation, dynamical constraints, and prospects for detection of this proposed object known as Planet Nine.

We use radial velocity observations to search for long-period gas giant companions in systems hosting inner super-Earth (1-4 R_Earth, 1-10 M_Earth) planets to constrain formation and migration scenarios for this population. We consistently re-fit published RV datasets for 65 stars and find 9 systems with statistically significant trends indicating the presence of an outer companion. We combine these RV data with AO images to constrain the masses and semi-major axes of these companions. We quantify our sensitivity to the presence of long-period companions by fitting the sample with a power law distribution and find an occurrence rate of 39+/-7% for companions 0.5-20 M_Jup and 1-20 AU. Half of our systems were discovered by the transit method and half were discovered by the RV method. While differences in RV baselines and number of data points between the two samples lead to different sensitivities to distant companions, we find that occurrence rates of gas giant companions in each sample are consistent at the 0.5$\sigma$ level. We compare the frequency of Jupiter analogs in these systems to the equivalent rate from field star surveys and find that Jupiter analogs are more common around stars hosting super-Earths. We conclude that the presence of outer gas giants does not suppress the formation of inner super-Earths, and that these two populations of planets instead appear to be correlated. We also find that the stellar metallicities of systems with gas giant companions are higher than those without companions, in agreement with the well-established metallicity correlation from RV surveys of field stars.

The most distant Kuiper belt objects appear to be clustered in longitude of perihelion and in orbital pole position. To date, the only two suggestions for the cause of these apparent clusterings have been either the effects of observational bias or the existence of the distant giant planet in an eccentric inclined orbit known as Planet Nine. To determine if observational bias can be the cause of these apparent clusterings, we develop a rigorous method of quantifying the observational biases in the observations of longitude of perihelion and orbital pole position. From this now more complete understanding of the biases we calculate that the probability that these distant Kuiper belt objects would be clustered as strongly as observed in both longitude of perihelion and in orbital pole position is only 0.2%. While explanations other than Planet Nine may someday be found, the statistical significance of this clustering is now difficult to discount.

More than two decades after the widespread detection of Jovian-class planets on short-period orbits around other stars, their dynamical origins remain imperfectly understood. In the traditional narrative, these highly irradiated giant planets, like Jupiter and Saturn, are envisioned to have formed at large stello-centric distances and to have subsequently undergone large-scale orbital decay. Conversely, more recent models propose that a large fraction of hot Jupiters could have formed via rapid gas accretion in their current orbital neighborhood. In this study, we examine the period-mass distribution of close-in giant planets, and demonstrate that the inner boundary of this population conforms to the expectations of the in-situ formation scenario. Specifically, we show that if conglomeration unfolds close to the disk's inner edge, the semi-major axis - mass relation of the emergent planets should follow a power law $a \propto M^{-2/7}$ - a trend clearly reflected in the data. We further discuss corrections to this relationship due to tidal decay of planetary orbits. Although our findings do not discount orbital migration as an active physical process, they suggest that the characteristic range of orbital migration experienced by giant planets is limited.

It has been proposed that mean motion resonances (MMRs) between Planet Nine and distant objects of the scattered disk might inform the semimajor axis and instantaneous position of Planet Nine. Within the context of this hypothesis, the specific distribution of occupied MMRs largely determines the available constraints. Here we characterize the behavior of scattered Kuiper Belt objects arising in the presence of an eccentric Planet Nine ($e_9 \in 0.1$, $0.7$), focusing on relative sizes of populations occupying particular commensurabilities. Highlighting the challenge of predicting the exact MMR of a given object, we find that the majority of resonant test particles have period ratios with Planet Nine other than those of the form $P_9/P=N/1$, $N/2$ $(N \in \mathbb{Z}^+)$. Taking into account the updated prior distribution of MMRs outlined in this work, we find that the close spacing of high-order resonances, as well as chaotic transport, preclude resonance-based Planet Nine constraints from current observational data.

K2-55b is a Neptune-sized planet orbiting a K7 dwarf with a radius of $0.715^{+0.043}_{-0.040}R_\odot$, a mass of $0.688\pm0.069 M_\odot$, and an effective temperature of $4300^{+107}_{-100}$K. Having characterized the host star using near-infrared spectra obtained at IRTF/SpeX, we observed a transit of K2-55b with Spitzer/IRAC and confirmed the accuracy of the original K2 ephemeris for future follow-up transit observations. Performing a joint fit to the Spitzer/IRAC and K2 photometry, we found a planet radius of $4.41^{+0.32}_{-0.28} R_\oplus$, an orbital period of $2.84927265_{-6.42\times10^{-6}}^{+6.87\times10^{-6}}$ days, and an equilibrium temperature of roughly 900K. We then measured the planet mass by acquiring twelve radial velocity (RV) measurements of the system using HIRES on the 10m Keck I Telescope. Our RV data set precisely constrains the mass of K2-55b to $43.13^{+5.98}_{-5.80} M_\oplus$, indicating that K2-55b has a bulk density of $2.8_{-0.6}^{+0.8}$ g cm$^{-3}$ and can be modeled as a rocky planet capped by a modest H/He envelope ($M_{\rm envelope} = 12\pm3\% M_p$). K2-55b is denser than most similarly sized planets, raising the question of whether the high planetary bulk density of K2-55b could be attributed to the high metallicity of K2-55. The absence of a substantial volatile envelope despite the large mass of K2-55b poses a challenge to current theories of gas giant formation. We posit that K2-55b may have escaped runaway accretion by migration, late formation, or inefficient core accretion or that K2-55b was stripped of its envelope by a late giant impact.

While planets between the size of Uranus and Saturn are absent within the Solar System, the star K2-24 hosts two such planets, K2-24b and c, with radii equal to $5.4~R_E$ and $7.5~R_E$, respectively. The two planets have orbital periods of 20.9 days and 42.4 days, residing only 1% outside the nominal 2:1 mean-motion resonance. In this work, we present results from a coordinated observing campaign to measure planet masses and eccentricities that combines radial velocity (RV) measurements from Keck/HIRES and transit-timing measurements from K2 and Spitzer. K2-24b and c have low, but non-zero, eccentricities of $e_1 \sim e_2 \sim 0.08$. The low observed eccentricities provide clues regarding the formation and dynamical evolution of K2-24b and K2-24c, suggesting that they could be the result of stochastic gravitational interactions with a turbulent protoplanetary disk, among other mechanisms. K2-24b and c are $19\pm2~M_E$ and $15\pm2~M_E$, respectively; K2-24c is 20% less massive than K2-24b, despite being 40% larger. Their large sizes and low masses imply large envelope fractions, which we estimate at $26^{+3}_{-3}\%$ and $52^{+5}_{-3}\%$. In particular, K2-24c's large envelope presents an intriguing challenge to the standard model of core nucleated accretion that predicts the onset of runaway accretion when $f_{env} \approx 50\%$.

HAT-P-11 is a mid-K dwarf that hosts one of the first Neptune-sized planets found outside the solar system. The orbit of HAT-P-11b is misaligned with the star's spin --- one of the few known cases of a misaligned planet orbiting a star less massive than the Sun. We find an additional planet in the system based on a decade of precision radial velocity (RV) measurements from Keck/HIRES. HAT-P-11c is similar to Jupiter in its mass ($M_P \sin{i} = 1.6\pm0.1$ $M_J$) and orbital period ($P = 9.3^{+1.0}_{-0.5}$ year), but has a much more eccentric orbit ($e=0.60\pm0.03$). In our joint modeling of RV and stellar activity, we found an activity-induced RV signal of $\sim$7 m s$^{-1}$, consistent with other active K dwarfs, but significantly smaller than the 31 m s$^{-1}$ reflex motion due to HAT-P-11c. We investigated the dynamical coupling between HAT-P-11b and c as a possible explanation for HAT-P-11b's misaligned orbit, finding that planet-planet Kozai interactions cannot tilt planet b's orbit due to general relativistic precession; however, nodal precession operating on million year timescales is a viable mechanism to explain HAT-P-11b's high obliquity. This leaves open the question of why HAT-P-11c may have such a tilted orbit. At a distance of 38 pc, the HAT-P-11 system offers rich opportunities for further exoplanet characterization through astrometry and direct imaging.

The observation that the orbits of long-period Kuiper Belt objects are anomalously clustered in physical space has recently prompted the Planet Nine hypothesis - the proposed existence of a distant and eccentric planetary member of our solar system. Within the framework of this model, a Neptune-like perturber sculpts the orbital distribution of distant Kuiper Belt objects through a complex interplay of resonant and secular effects, such that in addition to perihelion-circulating objects, the surviving orbits get organized into apsidally aligned and anti-aligned configurations with respect to Planet Nine's orbit. In this work, we investigate the role of Kuiper Belt initial conditions on the evolution of the outer solar system using numerical simulations. Intriguingly, we find that the final perihelion distance distribution depends strongly on the primordial state of the system, and demonstrate that a bimodal structure corresponding to the existence of both aligned and anti-aligned clusters is only reproduced if the initial perihelion distribution is assumed to extend well beyond $\sim 36$ AU. The bimodality in the final perihelion distance distribution is due to the existence of permanently stable objects, with the lower perihelion peak corresponding to the anti-aligned orbits and the higher perihelion peak corresponding to the aligned orbits. We identify the mechanisms which enable the persistent stability of these objects and locate the regions of phase space in which they reside. The obtained results contextualize the Planet Nine hypothesis within the broader narrative of solar system formation, and offer further insight into the observational search for Planet Nine.

Within the general framework of core-nucleated accretion theory of giant planet formation, the conglomeration of massive gaseous envelopes is facilitated by a transient period of rapid accumulation of nebular material. While the concurrent buildup of angular momentum is expected to leave newly formed planets spinning at near-breakup velocities, Jupiter and Saturn, as well as super-Jovian long-period extrasolar planets, are observed to rotate well below criticality. In this work, we demonstrate that the large luminosity of a young giant planet simultaneously leads to the generation of a strong planetary magnetic field, as well as thermal ionization of the circumplanetary disk. The ensuing magnetic coupling between the planetary interior and the quasi-Keplerian motion of the disk results in efficient braking of planetary rotation, with hydrodynamic circulation of gas within the Hill sphere playing the key role of expelling spin angular momentum to the circumstellar nebula. Our results place early-stage giant planet and stellar rotation within the same evolutionary framework, and motivate further exploration of magnetohydrodynamic phenomena in the context of the final stages of giant planet formation.

An understanding of the long-term evolution of self-gravitating disks ranks among the classic outstanding problems of astrophysics. In this work, we show that the secular inclination dynamics of a geometrically thin quasi-Keplerian disk, with a surface density profile that scales as the inverse square-root of the orbital radius, are described by the time-dependent Schr\"odinger equation. Within the context of this formalism, nodal bending waves correspond to the eigenmodes of a quasiparticle's wavefunction, confined in an infinite square well with boundaries given by the radial extent of the disk. We further show that external secular perturbations upon self-gravitating disks exhibit a mathematical similarity to quantum scattering theory. Employing this framework, we derive an analytic criterion for the gravitational rigidity of a nearly-Keplerian disk under external perturbations. Applications of the theory to circumstellar disks and Galactic nuclei are discussed.

The Kepler mission and its successor K2 have brought forth a cascade of transiting planets. Many of these planetary systems exhibit multiple members, but a large fraction possess only a single transiting example. This overabundance of singles has lead to the suggestion that up to half of Kepler systems might possess significant mutual inclinations between orbits, reducing the transiting number (the so-called "Kepler Dichotomy"). In a recent paper, Spalding & Batygin (2016) demonstrated that the quadrupole moment arising from a young, oblate star is capable of misaligning the constituent orbits of a close-in planetary system enough to reduce their transit number, provided that the stellar spin axis is sufficiently misaligned with respect to the planetary orbital plane. Moreover, tightly packed planetary systems were shown to be susceptible to becoming destabilized during this process. Here, we investigate the ubiquity of the stellar obliquity-driven instability within systems with a range of multiplicities. We find that most planetary systems analysed, including those possessing only 2 planets, underwent instability for stellar spin periods below ~3 days and stellar tilts of order 30 degrees. Moreover, we are able to place upper limits on the stellar obliquity in systems such as K2-38 (obliquity <20 degrees), where other methods of measuring spin-orbit misalignment are not currently available. Given the known parameters of T-Tauri stars, we predict that up to 1/2 of super-Earth mass systems may encounter the instability, in general agreement with the fraction typically proposed to explain the observed abundance of single-transiting systems.

The M-dwarf GJ 876 is the closest known star to harbor a multi-planetary system. With three outer planets locked in a chaotic Laplace-type resonance and an appreciably eccentric short-period Super-Earth, this system represents a unique exposition of extrasolar planetary dynamics. A key question that concerns the long-term evolution of this system, and the fate of close-in planets in general, is how the significant eccentricity of the inner-most planet is maintained against tidal circularization on timescales comparable to the age of the universe. Here, we employ stochastic secular perturbation theory and N-body simulations to show that the orbit of the inner-most planet is shaped by a delicate balance between extrinsic chaotic forcing and tidal dissipation. As such, the planet's orbital eccentricity represents an indirect measure of its tidal quality factor. Based on the system's present-day architecture, we estimate that the extrasolar Super-Earth GJ 876d has a tidal Q~10,000-100,000 - a value characteristic of solar system gas giants.

Over the course of last decade, the Nice model has dramatically changed our view of the solar system's formation and early evolution. Within the context of this model, a transient period of planet-planet scattering is triggered by gravitational interactions between the giant planets and a massive primordial planetesimal disk, leading to a successful reproduction of the solar system's present-day architecture. In typical realizations of the Nice model, self-gravity of the planetesimal disk is routinely neglected, as it poses a computational bottleneck to the calculations. Recent analyses have shown, however, that a self-gravitating disk can exhibit behavior that is dynamically distinct, and this disparity may have significant implications for the solar system's evolutionary path. In this work, we explore this discrepancy utilizing a large suite of Nice odel simulations with and without a self-gravitating planetesimal disk, taking advantage of the inherently parallel nature of graphic processing units. Our simulations demonstrate that self-consistent modeling of particle interactions does not lead to significantly different final planetary orbits from those obtained within conventional simulations. Moreover, self-gravitating calculations show similar planetesimal evolution to non-self-gravitating numerical experiments after dynamical instability is triggered, suggesting that the orbital clustering observed in the distant Kuiper belt is unlikely to have a self-gravitational origin.

Surveys of young star-forming regions have discovered a growing population of planetary-mass (<13 M_Jup) companions around young stars. There is an ongoing debate as to whether these companions formed like planets (that is, from the circumstellar disk), or if they represent the low-mass tail of the star formation process. In this study we utilize high-resolution spectroscopy to measure rotation rates of three young (2-300 Myr) planetary-mass companions and combine these measurements with published rotation rates for two additional companions to provide a look at the spin distribution of these objects. We compare this distribution to complementary rotation rate measurements for six brown dwarfs with masses <20 M_Jup, and show that these distributions are indistinguishable. This suggests that either that these two populations formed via the same mechanism, or that processes regulating rotation rates are independent of formation mechanism. We find that rotation rates for both populations are well below their break-up velocities and do not evolve significantly during the first few hundred million years after the end of accretion. This suggests that rotation rates are set during late stages of accretion, possibly by interactions with a circumplanetary disk. This result has important implications for our understanding of the processes regulating the angular momentum evolution of young planetary-mass objects, and of the physics of gas accretion and disk coupling in the planetary-mass regime.

The arrival of the robustly hyperbolic asteroid A/2017 U1 has potentially interesting ramifications for the planet-formation process. Although extrapolations from a sample size of one are necessarily uncertain, order-of-magnitude estimates suggest that the Galaxy contains a substantial mass in similar bodies. We argue that despite its lack of Coma, A/2017 U1 likely contained a significant mass fraction of volatile components, and we argue that its presence can be used to infer a potentially large population of as-yet undetected Neptune-like extrasolar planets.

We describe a joint high contrast imaging survey for planets at Keck and VLT of the last large sample of debris disks identified by the Spitzer Space Telescope. No new substellar companions were discovered in our survey of 30 Spitzer-selected targets. We combine our observations with data from four published surveys to place constraints on the frequency of planets around 130 debris disk single stars, the largest sample to date. For a control sample, we assembled contrast curves from several published surveys targeting 277 stars which do not show infrared excesses. We assumed a double power law distribution in mass and semi-major axis of the form f(m,a) = $Cm^{\alpha}a^{\beta}$, where we adopted power law values and logarithmically flat values for the mass and semi-major axis of planets. We find that the frequency of giant planets with masses 5-20 $M_{\rm Jup}$ and separations 10-1000 AU around stars with debris disks is 6.27% (68% confidence interval 3.68 - 9.76%), compared to 0.73% (68% confidence interval 0.20 - 1.80%) for the control sample of stars without disks. These distributions differ at the 88% confidence level, tentatively suggesting distinctness of these samples.

The observational census of trans-Neptunian objects with semi-major axes greater than ~250 AU exhibits unexpected orbital structure that is most readily attributed to gravitational perturbations induced by a yet-undetected, massive planet. Although the capacity of this planet to (i) reproduce the observed clustering of distant orbits in physical space, (ii) facilitate dynamical detachment of their perihelia from Neptune, and (iii) excite a population of long-period centaurs to extreme inclinations is well established through numerical experiments, a coherent theoretical description of the dynamical mechanisms responsible for these effects remains elusive. In this work, we characterize the dynamical processes at play, from semi-analytic grounds. We begin by considering a purely secular model of orbital evolution induced by Planet Nine, and show that it is at odds with the ensuing stability of distant objects. Instead, the long-term survival of the clustered population of long-period KBOs is enabled by a web of mean-motion resonances driven by Planet Nine. Then, by taking a compact-form approach to perturbation theory, we show that it is the secular dynamics embedded within these resonances that regulates the orbital confinement and perihelion detachment of distant Kuiper belt objects. Finally, we demonstrate that the onset of large-amplitude oscillations of orbital inclinations is accomplished through capture of low-inclination objects into a high-order secular resonance and identify the specific harmonic that drives the evolution. In light of the developed qualitative understanding of the governing dynamics, we offer an updated interpretation of the current observational dataset within the broader theoretical framework of the Planet Nine hypothesis.

Despite decades of inquiry, the origin of giant planets residing within a few tenths of an astronomical unit from their host stars remains unclear. Traditionally, these objects are thought to have formed further out before subsequently migrating inwards. However, the necessity of migration has been recently called into question with the emergence of in-situ formation models of close-in giant planets. Observational characterization of the transiting sub-sample of close-in giants has revealed that "warm" Jupiters, possessing orbital periods longer than roughly 10 days more often possess close-in, co-transiting planetary companions than shorter period "hot" Jupiters, that are usually lonely. This finding has previously been interpreted as evidence that smooth, early migration or in situ formation gave rise to warm Jupiter-hosting systems, whereas more violent, post-disk migration pathways sculpted hot Jupiter-hosting systems. In this work, we demonstrate that both classes of planet may arise via early migration or in-situ conglomeration, but that the enhanced loneliness of hot Jupiters arises due to a secular resonant interaction with the stellar quadrupole moment. Such an interaction tilts the orbits of exterior, lower mass planets, removing them from transit surveys where the hot Jupiter is detected. Warm Jupiter-hosting systems, in contrast, retain their coplanarity due to the weaker influence of the host star's quadrupolar potential relative to planet-disk interactions. In this way, hot Jupiters and warm Jupiters are placed within a unified theoretical framework that may be readily validated or falsified using data from upcoming missions such as TESS.

HD 3167 is a bright (V = 8.9), nearby K0 star observed by the NASA K2 mission (EPIC 220383386), hosting two small, short-period transiting planets. Here we present the results of a multi-site, multi-instrument radial velocity campaign to characterize the HD 3167 system. The masses of the transiting planets are 5.02+/-0.38 MEarth for HD 3167 b, a hot super-Earth with a likely rocky composition (rho_b = 5.60+2.15-1.43 g/cm^3), and 9.80+1.30-1.24 MEarth for HD 3167 c, a warm sub-Neptune with a likely substantial volatile complement (rho_c = 1.97+0.94-0.59 g/cm^3). We explore the possibility of atmospheric composition analysis and determine that planet c is amenable to transmission spectroscopy measurements, and planet b is a potential thermal emission target. We detect a third, non-transiting planet, HD 3167 d, with a period of 8.509+/-0.045 d (between planets b and c) and a minimum mass of 6.90+/-0.71 MEarth. We are able to constrain the mutual inclination of planet d with planets b and c: we rule out mutual inclinations below 1.3 degrees as we do not observe transits of planet d. From 1.3-40 degrees, there are viewing geometries invoking special nodal configurations which result in planet d not transiting some fraction of the time. From 40-60 degrees, Kozai-Lidov oscillations increase the system's instability, but it can remain stable for up to 100Myr. Above 60 degrees, the system is unstable. HD 3167 promises to be a fruitful system for further study and a preview of the many exciting systems expected from the upcoming NASA TESS mission.

Our Keck/NIRC2 imaging survey searches for stellar companions around 144 systems with radial velocity (RV) detected giant planets to determine whether stellar binaries influence the planets' orbital parameters. This survey, the largest of its kind to date, finds eight confirmed binary systems and three confirmed triple systems. These include three new multi-stellar systems (HD 30856, HD 86081, and HD 207832) and three multi-stellar systems with newly confirmed common proper motion (HD 43691, HD 116029, and HD 164509). We combine these systems with seven RV planet-hosting multi-stellar systems from the literature in order to test for differences in the properties of planets with semimajor axes ranging between 0.1-5 au in single vs multi-stellar systems. We find no evidence that the presence or absence of stellar companions alters the distribution of planet properties in these systems. Although the observed stellar companions might influence the orbits of more distant planetary companions in these systems, our RV observations currently provide only weak constraints on the masses and orbital properties of planets beyond 5 au. In order to aid future efforts to characterize long period RV companions in these systems, we publish our contrast curves for all 144 targets. Using four years of astrometry for six hierarchical triple star systems hosting giant planets, we fit the orbits of the stellar companions in order to characterize the orbital architecture in these systems. We find that the orbital plane of the secondary and tertiary companions are inconsistent with an edge-on orbit in four out of six cases.

Extrasolar planets on eccentric short-period orbits provide a laboratory in which to study radiative and tidal interactions between a planet and its host star under extreme forcing conditions. Studying such systems probes how the planet's atmosphere redistributes the time-varying heat flux from its host and how the host star responds to transient tidal distortion. Here, we report the insights into the planet-star interactions in HAT-P-2's eccentric planetary system gained from the analysis of 350 hr of 4.5 micron observations with the Spitzer Space Telescope. The observations show no sign of orbit-to-orbit variability nor of orbital evolution of the eccentric planetary companion, HAT-P-2 b. The extensive coverage allows us to better differentiate instrumental systematics from the transient heating of HAT-P-2 b's 4.5 micron photosphere and yields the detection of stellar pulsations with an amplitude of approximately 40 ppm. These pulsation modes correspond to exact harmonics of the planet's orbital frequency, indicative of a tidal origin. Transient tidal effects can excite pulsation modes in the envelope of a star, but, to date, such pulsations had only been detected in highly eccentric stellar binaries. Current stellar models are unable to reproduce HAT-P-2's pulsations, suggesting that our understanding of the interactions at play in this system is incomplete.

Mean motion commensurabilities in multi-planet systems are an expected outcome of protoplanetary disk-driven migration, and their relative dearth in the observational data presents an important challenge to current models of planet formation and dynamical evolution. One natural mechanism that can lead to the dissolution of commensurabilities is stochastic orbital forcing, induced by turbulent density fluctuations within the nebula. While this process is qualitatively promising, the conditions under which mean motion resonances can be broken are not well understood. In this work, we derive a simple analytic criterion that elucidates the relationship among the physical parameters of the system, and find the conditions necessary to drive planets out of resonance. Subsequently, we confirm our findings with numerical integrations carried out in the perturbative regime, as well as direct N-body simulations. Our calculations suggest that turbulent resonance disruption depends most sensitively on the planet-star mass ratio. Specifically, for a disk with properties comparable to the early solar nebula with $\alpha=0.01$, only planet pairs with cumulative mass ratios smaller than $(m_1+m_2)/M\lesssim10^{-5}\sim3M_{\oplus}/M_{\odot}$ are susceptible to breaking resonance at semi-major axis of order $a\sim0.1\,$AU. Although turbulence can sometimes compromise resonant pairs, an additional mechanism (such as suppression of resonance capture probability through disk eccentricity) is required to adequately explain the largely non-resonant orbital architectures of extrasolar planetary systems.

The trans-Neptunian region of the solar system exhibits an intricate dynamical structure, much of which can be explained by an instability-driven orbital history of the giant planets. However, the origins of a highly inclined, and in certain cases retrograde, population of trans-Neptunian objects remain elusive within the framework of this evolutionary picture. In this work, we show that the existence of a distant, Neptune-like planet that resides on an eccentric and mildly inclined orbit fully accounts for the anomalous component the trans-Neptunian orbital distribution. Adopting the same parameters for Planet Nine as those previously invoked to explain the clustering of distant Kuiper belt orbits in physical space, we carry out a series of numerical experiments which elucidate the physical process though which highly inclined Kuiper belt objects with semi-major axes smaller than 100 AU are generated. The identified dynamical pathway demonstrates that enigmatic members of the Kuiper belt such as Drac and Niku are derived from the extended scattered disk of the solar system.

Hot Jupiters, giant extrasolar planets with orbital periods shorter than ~10 days, have long been thought to form at large radial distances, only to subsequently experience long-range inward migration. Here, we propose that in contrast with this picture, a substantial fraction of the hot Jupiter population formed in situ via the core accretion process. We show that under conditions appropriate to the inner regions of protoplanetary disks, rapid gas accretion can be initiated by Super-Earth type planets, comprising 10-20 Earth masses of refractory composition material. An in situ formation scenario leads to testable consequences, including the expectation that hot Jupiters should frequently be accompanied by additional low-mass planets with periods shorter than ~100 days. Our calculations further demonstrate that dynamical interactions during the early stages of planetary systems' lifetimes should increase the inclinations of such companions, rendering transits rare. High-precision radial velocity monitoring provides the best prospect for their detection.

The six-degree obliquity of the sun suggests that either an asymmetry was present in the solar system's formation environment, or an external torque has misaligned the angular momentum vectors of the sun and the planets. However, the exact origin of this obliquity remains an open question. Batygin & Brown (2016) have recently shown that the physical alignment of distant Kuiper Belt orbits can be explained by a 5-20 Earth-mass planet on a distant, eccentric, and inclined orbit, with an approximate perihelion distance of ~250 AU. Using an analytic model for secular interactions between Planet Nine and the remaining giant planets, here we show that a planet with similar parameters can naturally generate the observed obliquity as well as the specific pole position of the sun's spin axis, from a nearly aligned initial state. Thus, Planet Nine offers a testable explanation for the otherwise mysterious spin-orbit misalignment of the solar system.

During its 5 year mission, the Kepler spacecraft has uncovered a diverse population of planetary systems with orbital configurations ranging from single-transiting planets to systems of multiple planets co-transiting the parent star. By comparing the relative occurrences of multiple to single-transiting systems, recent analyses have revealed a significant over-abundance of singles. Dubbed the "Kepler Dichotomy," this feature has been interpreted as evidence for two separate populations of planetary systems: one where all orbits are confined to a single plane, and a second where the constituent planetary orbits possess significant mutual inclinations, allowing only a single member to be observed in transit at a given epoch. In this work, we demonstrate that stellar obliquity, excited within the disk-hosting stage, can explain this dichotomy. Young stars rotate rapidly, generating a significant quadrupole moment which torques the planetary orbits, with inner planets influenced more strongly. Given nominal parameters, this torque is sufficiently strong to excite significant mutual inclinations between planets, enhancing the number of single-transiting planets, sometimes through a dynamical instability. Furthermore, as hot stars appear to possess systematically higher obliquities, we predict that single-transiting systems should be relatively more prevalent around more massive stars. We analyze the Kepler data and confirm this signal to be present.

Stellar companions can influence the formation and evolution of planetary systems, but there are currently few observational constraints on the properties of planet-hosting binary star systems. We search for stellar companions around 77 transiting hot Jupiter systems to explore the statistical properties of this population of companions as compared to field stars of similar spectral type. After correcting for survey incompleteness, we find that $47\%\pm7\%$ of hot Jupiter systems have stellar companions with semi-major axes between 50-2000 AU. This is 2.9 times larger than the field star companion fraction in this separation range, with a significance of $4.4\sigma$. In the 1-50AU range, only $3.9^{+4.5}_{-2.0}\%$ of hot Jupiters host stellar companions compared to the field star value of $16.4\%\pm0.7\%$, which is a $2.7\sigma$ difference. We find that the distribution of mass ratios for stellar companions to hot Jupiter systems peaks at small values and therefore differs from that of field star binaries which tend to be uniformly distributed across all mass ratios. We conclude that either wide separation stellar binaries are more favorable sites for gas giant planet formation at all separations, or that the presence of stellar companions preferentially causes the inward migration of gas giant planets that formed farther out in the disk via dynamical processes such as Kozai-Lidov oscillations. We determine that less than 20% of hot Jupiters have stellar companions capable of inducing Kozai-Lidov oscillations assuming initial semi-major axes between 1-5 AU, implying that the enhanced companion occurrence is likely correlated with environments where gas giants can form efficiently.

We use an extensive suite of numerical simulations to constrain the mass and orbit of Planet Nine, the recently proposed perturber in a distant eccentric orbit in the outer solar system. We compare our simulations to the observed population of aligned eccentric high semimajor axis Kuiper belt objects and determine which simulation parameters are statistically compatible with the observations. We find that only a narrow range of orbital elements can reproduce the observations. In particular, the combination of semimajor axis, eccentricity, and mass of Planet Nine strongly dictates the semimajor axis range of the orbital confinement of the distant eccentric Kuiper belt objects. Allowed orbits, which confine Kuiper belt objects with semimajor axis beyond 380~AU, have perihelia roughly between 150 and 350~AU, semimajor axes between 380 and 980~AU, and masses between 5 and 20 Earth masses. Orbitally confined objects also generally have orbital planes similar to that of the planet, suggesting that the planet is inclined approximately 30~degrees to the ecliptic. We compare the allowed orbital positions and estimated brightness of Planet Nine to previous and ongoing surveys which would be sensitive to the planet's detection and use these surveys to rule out approximately two-thirds of the planet's orbit. Planet Nine is likely near aphelion with an approximate brightness of $22<V<25$. At opposition, its motion, mainly due to parallax, can easily be detected within 24 hours.

Two decades after the discovery of 51 Peg b, the formation processes and atmospheres of short-period gas giants remain poorly understood. Observations of eccentric systems provide key insights on those topics as they can illuminate how a planet's atmosphere responds to changes in incident flux. We report here the analysis of multi-day multi-channel photometry of the eccentric (e~ 0.93) hot Jupiter HD 80606 b obtained with the Spitzer Space Telescope. The planet's extreme eccentricity combined with the long coverage and exquisite precision of new periastron-passage observations allow us to break the degeneracy between the radiative and dynamical timescales of HD 80606 b's atmosphere and constrain its global thermal response. Our analysis reveals that the atmospheric layers probed heat rapidly (~4 hr radiative timescale) from $\lt$500 to 1400 K as they absorb ~ 20% of the incoming stellar flux during the periastron passage, while the planet's rotation period is 93$\pm_{35}^{85}$ hr, which exceeds the predicted pseudo-synchronous period (40 hr).

HAT-P-13b is a Jupiter-mass transiting exoplanet that has settled onto a stable, short-period, and mildly eccentric orbit as a consequence of the action of tidal dissipation and perturbations from a second, highly eccentric, outer companion. Due to the special orbital configuration of the HAT-P-13 system, the magnitude of HAT-P-13b's eccentricity ($e_b$) is in part dictated by its Love number ($k_{2_b}$), which is in turn a proxy for the degree of central mass concentration in its interior. Thus, the measurement of $e_b$ constrains $k_{2_b}$ and allows us to place otherwise elusive constraints on the mass of HAT-P-13b's core ($M_{\rm{core,b}}$). In this study we derive new constraints on the value of $e_b$ by observing two secondary eclipses of HAT-P-13b with the Infrared Array Camera on board the $\textit{Spitzer Space Telescope}$. We fit the measured secondary eclipse times simultaneously with radial velocity measurements and find that $e_b = 0.00700 \pm 0.00100$. We then use octupole-order secular perturbation theory to find the corresponding $k_{2_b} = 0.31^{+0.08}_{-0.05}$. Applying structural evolution models, we then find, with 68\% confidence, that $M_{\rm{core,b}}$ is less than 25 Earth masses ($M_{\oplus}$). The most likely value of $M_{\rm{core,b}} = 11 M_{\oplus}$, which is similar to the core mass theoretically required for runaway gas accretion. This is the tightest constraint to date on the core mass of a hot Jupiter. Additionally, we find that the measured secondary eclipse depths, which are in the 3.6 $\mu$m and 4.5 $\mu$m bands, best match atmospheric model predictions with a dayside temperature inversion and relatively efficient day-night circulation.

We conducted a Doppler survey at Keck combined with NIRC2 K-band AO imaging to search for massive, long-period companions to 123 known exoplanet systems with one or two planets detected using the radial velocity (RV) method. Our survey is sensitive to Jupiter mass planets out to 20 AU for a majority of stars in our sample, and we report the discovery of eight new long-period planets, in addition to 20 systems with statistically significant RV trends indicating the presence of an outer companion beyond 5 AU. We combine our RV observations with AO imaging to determine the range of allowed masses and orbital separations for these companions, and account for variations in our sensitivity to companions among stars in our sample. We estimate the total occurrence rate of companions in our sample to be 52 +/- 5% over the range 1 - 20 M_Jup and 5 - 20 AU. Our data also suggest a declining frequency for gas giant planets in these systems beyond 3-10 AU, in contrast to earlier studies that found a rising frequency for giant planets in the range 0.01-3 AU. This suggests either that the frequency of gas giant planets peaks between 3-10 AU, or that outer companions in these systems have a different semi-major axis distribution than the overall gas giant planet population. Our results also suggest that hot gas giants may be more likely to have an outer companion than cold gas giants. We find that planets with an outer companion have higher average eccentricities than their single counterparts, suggesting that dynamical interactions between planets may play an important role in these systems.

Recent analyses have shown that distant orbits within the scattered disk population of the Kuiper belt exhibit an unexpected clustering in their respective arguments of perihelion. While several hypotheses have been put forward to explain this alignment, to date, a theoretical model that can successfully account for the observations remains elusive. In this work we show that the orbits of distant Kuiper belt objects cluster not only in argument of perihelion, but also in physical space. We demonstrate that the perihelion positions and orbital planes of the objects are tightly confined and that such a clustering has only a probability of 0.007% to be due to chance, thus requiring a dynamical origin. We find that the observed orbital alignment can be maintained by a distant eccentric planet with mass greater than ~10 Earth masses, whose orbit lies in approximately the same plane as those of the distant Kuiper belt objects, but whose perihelion is 180 degrees away from the perihelia of the minor bodies. In addition to accounting for the observed orbital alignment, the existence of such a planet naturally explains the presence of high perihelion Sedna-like objects, as well as the known collection of high semimajor axis objects with inclinations between 60 and 150 degrees whose origin was previously unclear. Continued analysis of both distant and highly inclined outer solar system objects provides the opportunity for testing our hypothesis as well as further constraining the orbital elements and mass of the distant planet.

Jupiter and Saturn play host to an impressive array of satellites, making it reasonable to suspect that similar systems of moons might exist around giant extrasolar planets. Furthermore, a significant population of such planets is known to reside at distances of several Astronomical Units (AU), leading to speculation that some moons thereof might support liquid water on their surfaces. However, giant planets are thought to undergo inward migration within their natal protoplanetary disks, suggesting that gas giants currently occupying their host star's habitable zone formed further out. Here we show that when a moon-hosting planet undergoes inward migration, dynamical interactions may naturally destroy the moon through capture into a so-called "evection resonance." Within this resonance, the lunar orbit's eccentricity grows until the moon eventually collides with the planet. Our work suggests that moons orbiting within about 10 planetary radii are susceptible to this mechanism, with the exact number dependent upon the planetary mass, oblateness and physical size. Whether moons survive or not is critically related to where the planet began its inward migration as well as the character of inter-lunar perturbations. For example, a Jupiter-like planet currently residing at 1AU could lose moons if it formed beyond 5AU. Cumulatively, we suggest that an observational census of exomoons could potentially inform us on the extent of inward planetary migration, for which no reliable observational proxy currently exists.

Surveys of nearby field stars indicate that stellar binaries are common, yet little is known about the effects that these companions may have on planet formation and evolution. The Friends of Hot Jupiters project uses three complementary techniques to search for stellar companions to known planet-hosting stars: radial velocity monitoring, adaptive optics imaging, and near-infrared spectroscopy. In this paper, we examine high-resolution K band infrared spectra of fifty stars hosting gas giant planets on short-period orbits. We use spectral fitting to search for blended lines due to the presence of cool stellar companions in the spectra of our target stars, where we are sensitive to companions with temperatures between 3500-5000 K and projected separations less than 100 AU in most systems. We identify eight systems with candidate low-mass companions, including one companion that was independently detected in our AO imaging survey. For systems with radial velocity accelerations, a spectroscopic non-detection rules out scenarios involving a stellar companion in a high inclination orbit. We use these data to place an upper limit on the stellar binary fraction at small projected separations, and show that the observed population of candidate companions is consistent with that of field stars and also with the population of wide-separation companions detected in our previous AO survey. We find no evidence that spectroscopic stellar companions are preferentially located in systems with short-period gas giant planets on eccentric and/or misaligned orbits.

Detailed observational characterization of transiting exoplanet systems has revealed that the spin-axes of massive (M > ~1.2 solar masses) stars often exhibit substantial misalignments with respect to the orbits of the planets they host. Conversely, lower-mass stars tend to only have limited obliquities. A similar trend has recently emerged within the observational dataset of young stars' magnetic field strengths: massive T-Tauri stars tend to have dipole fields that are ~10 times weaker than their less-massive counterparts. Here we show that the associated dependence of magnetic star-disk torques upon stellar mass naturally explains the observed spin-orbit misalignment trend, provided that misalignments are obtained within the disk-hosting phase. Magnetic torques act to realign the stellar spin-axes of lower-mass stars with the disk plane on a timescale significantly shorter than the typical disk lifetime, whereas the same effect operates on a much longer timescale for massive stars. Cumulatively, our results point to a primordial excitation of extrasolar spin-orbit misalignment, signalling consistency with disk-driven migration as the dominant transport mechanism for short-period planets. Furthermore, we predict that spin-orbit misalignments in systems where close-in planets show signatures of dynamical, post-nebular emplacement will not follow the observed correlation with stellar mass.

The richness of dynamical behavior exhibited by the rotational states of various solar system objects has driven significant advances in the theoretical understanding of their evolutionary histories. An important factor that determines whether a given object is prone to exhibiting non-trivial rotational evolution is the extent to which such an object can maintain a permanent aspheroidal shape, meaning that exotic behavior is far more common among the small body populations of the solar system. Gravitationally bound binary objects constitute a substantial fraction of asteroidal and TNO populations, comprising systems of triaxial satellites that orbit permanently deformed central bodies. In this work, we explore the rotational evolution of such systems with specific emphasis on quadrupole-quadrupole interactions, and show that for closely orbiting, highly deformed objects, both prograde and retrograde spin-spin resonances naturally arise. Subsequently, we derive capture probabilities for leading order commensurabilities and apply our results to the illustrative examples of (87) Sylvia and (216) Kleopatra asteroid systems. Cumulatively, our results suggest that spin-spin coupling may be consequential for highly elongated, tightly orbiting binary objects.

We consider the dynamical evolution of two planets orbiting in the vicinity of a first order mean motion reso- nance while simultaneously undergoing eccentricity damping and convergent migration. Following Goldreich & Schlichting (2014), we include a coupling between the dissipative semimajor axis evolution and the damping of the eccentricities. In agreement with past studies, we find that this coupling can lead to overstability of the resonance and that for a certain range of parameters capture into resonance is only temporary. Using a more general model, we show that whether overstable motion can occur depends in a characteristic way on the mass ratio between the two planets as well as their relative eccentricity damping timescales. Moreover, we show that even when escape from resonance does occur, the timescale for escape is long enough such at any given time a pair of planets is more likely to be found in a resonance rather than migrating between them. Thus, we argue that overstability of resonances cannot singlehandedly reconcile convergent migration with the observed lack of Kepler planet pairs found near resonances. However, it is possible that overstable motion in combination with other effects such as large scale orbital instability could produce the observed period ratio distribution.

The early stages of dynamical evolution of planetary systems are often shaped by dissipative processes that drive orbital migration. In multi-planet systems, convergent amassing of orbits inevitably leads to encounters with rational period ratios, which may result in establishment of mean motion resonances. The success or failure of resonant capture yields exceedingly different subsequent evolutions, and thus plays a central role in determining the ensuing orbital architecture of planetary systems. In this work, we employ an integrable Hamiltonian formalism for first order planetary resonances that allows both secondary bodies to have finite masses and eccentricities, and construct a comprehensive theory for resonant capture. Particularly, we derive conditions under which orbital evolution lies within the adiabatic regime, and provide a generalized criterion for guaranteed resonant locking as well as a procedure for calculating capture probabilities when capture is not certain. Subsequently, we utilize the developed analytical model to examine the evolution of Jupiter and Saturn within the protosolar nebula, and investigate the origins of the dominantly non-resonant orbital distribution of sub-Jovian extrasolar planets. Our calculations show that the commonly observed extrasolar orbital structure can be understood if planet pairs encounter mean motion commensurabilities on slightly eccentric (e~0.02) orbits. Accordingly, we speculate that resonant capture among low-mass planets is typically rendered unsuccessful due to subtle axial asymmetries inherent to the global structure of protoplanetary disks.

The statistics of extrasolar planetary systems indicate that the default mode of planet formation generates planets with orbital periods shorter than 100 days, and masses substantially exceeding that of the Earth. When viewed in this context, the Solar System is unusual. Here, we present simulations which show that a popular formation scenario for Jupiter and Saturn, in which Jupiter migrates inward from a > 5 AU to a ~ 1.5 AU before reversing direction, can explain the low overall mass of the Solar System's terrestrial planets, as well as the absence of planets with a < 0.4 AU. Jupiter's inward migration entrained s ~ 10-100 km planetesimals into low-order mean-motion resonances, shepherding and exciting their orbits. The resulting collisional cascade generated a planetesimal disk that, evolving under gas drag, would have driven any pre-existing short-period planets into the Sun. In this scenario, the Solar System's terrestrial planets formed from gas-starved mass-depleted debris that remained after the primary period of dynamical evolution.

The GJ876 system was among the earliest multi-planetary detections outside of the Solar System, and has long been known to harbor a resonant pair of giant planets. Subsequent characterization of the system revealed the presence of an additional Neptune mass object on an external orbit, locked in a three body Laplace mean motion resonance with the previously known planets. While this system is currently the only known extrasolar example of a Laplace resonance, it differs from the Galilean satellites in that the orbital motion of the planets is known to be chaotic. In this work, we present a simple perturbative model that illuminates the origins of stochasticity inherent to this system and derive analytic estimates of the Lyapunov time as well as the chaotic diffusion coefficient. We then address the formation of the multi-resonant structure within a protoplanetary disk and show that modest turbulent forcing in addition to dissipative effects is required to reproduce the observed chaotic configuration. Accordingly, this work places important constraints on the typical formation environments of planetary systems and informs the attributes of representative orbital architectures that arise from extended disk-driven evolution.

Multi-star systems are common, yet little is known about a stellar companion's influence on the formation and evolution of planetary systems. For instance, stellar companions may have facilitated the inward migration of hot Jupiters towards to their present day positions. Many observed short period gas giant planets also have orbits that are misaligned with respect to their star's spin axis, which has also been attributed to the presence of a massive outer companion on a non-coplanar orbit. We present the results of a multi-band direct imaging survey using Keck NIRC2 to measure the fraction of short period gas giant planets found in multi-star systems. Over three years, we completed a survey of 50 targets ("Friends of Hot Jupiters") with 27 targets showing some signature of multi-body interaction (misaligned or eccentric orbits) and 23 targets in a control sample (well-aligned and circular orbits). We report the masses, projected separations, and confirmed common proper motion for the 19 stellar companions found around 17 stars. Correcting for survey incompleteness, we report companion fractions of $48\%\pm9\%$, $47\%\pm12\%$, and $51\%\pm13\%$ in our total, misaligned/eccentric, and control samples, respectively. This total stellar companion fraction is $2.8\,\sigma$ larger than the fraction of field stars with companions approximately $50-2000\,$AU. We observe no correlation between misaligned/eccentric hot Jupiter systems and the incidence of stellar companions. Combining this result with our previous radial velocity survey, we determine that $72\% \pm 16\%$ of hot Jupiters are part of multi-planet and/or multi-star systems.

Star formation proceeds via the collapse of a molecular cloud core over multiple dynamical timescales. Turbulence within cores results in a spatially non-uniform angular momentum of the cloud, causing a stochastic variation in orientation of the disk forming from the collapsing material. In the absence of star-disk angular momentum coupling, such disk-tilting would provide a natural mechanism for production of primordial spin-orbit misalignments in the resulting planetary systems. However, owing to high accretion rates in the embedded phase of star formation, the inner edge of the circumstellar disk extends down to the stellar surface, resulting in efficient gravitational and accretional angular momentum transfer between the star and the disk. Here, we demonstrate that the resulting gravitational coupling is sufficient to suppress any significant star-disk misalignment, with accretion playing a secondary role. The joint tilting of the star-disk system leads to a stochastic wandering of star-aligned bipolar outflows. Such wandering widens the effective opening angle of stellar outflows, allowing for more efficient clearing of the remainder of the protostar's gaseous envelope. Accordingly, the processes described in this work provide an additional mechanism responsible for sculpting the stellar Initial Mass Function (IMF).

On timescales that greatly exceed an orbital period, typical planetary orbits evolve in a stochastic yet stable fashion. On even longer timescales, however, planetary orbits can spontaneously transition from bounded to unbound chaotic states. Large-scale instabilities associated with such behavior appear to play a dominant role in shaping the architectures of planetary systems, including our own. Here we show how such transitions are possible, focusing on the specific case of the long-term evolution of Mercury. We develop a simple analytical model for Mercury's dynamics and elucidate the origins of its short term stochastic behavior as well as of its sudden progression to unbounded chaos. Our model allows us to estimate the timescale on which this transition is likely to be triggered, i.e. the dynamical lifetime of the Solar System as we know it. The formulated theory is consistent with the results of numerical simulations and is broadly applicable to extrasolar planetary systems dominated by secular interactions. These results constitute a significant advancement in our understanding of the processes responsible for sculpting of the dynamical structures of generic planetary systems.

The Earth's obliquity is stabilized by the Moon, which facilitates a rapid precession of the Earth's spin-axis, de-tuning the system away from resonance with orbital modulation. It is however, likely that the architecture of the Solar System underwent a dynamical instability-driven transformation, where the primordial configuration was more compact. Hence, the characteristic frequencies associated with orbital perturbations were likely faster in the past, potentially allowing for secular resonant encounters. In this work we examine if at any point in the Earth's evolutionary history, the obliquity varied significantly. Our calculations suggest that even though the orbital perturbations were different, the system nevertheless avoided resonant encounters throughout its evolution. This indicates that the Earth obtained its current obliquity during the formation of the Moon.

Giant planets that reside in close proximity to their host stars are subject to extreme irradiation, which gives rise to thermal ionization of trace Alkali metals in their atmospheres. On objects where the atmospheric electrical conductivity is substantial, the global circulation couples to the background magnetic field, inducing supplementary fields and altering the nature of the flow. To date, a number of authors have considered the influence of a spin-pole aligned dipole magnetic field on the dynamical state of a weakly-ionized atmosphere and found that magnetic breaking may lead to significantly slower winds than predicted within a purely hydrodynamical framework. Here, we consider the effect of a tilted dipole magnetic field on the circulation and demonstrate that in addition to regulating wind velocities, an oblique field generates stationary non-axisymmetric structures that adhere to the geometry of the magnetic pole. Using a kinematic perturbative approach, we derive a closed-form solution for the perturbed circulation and show that the fractional distortion of zonal jets scales as the product of the field obliquity and the Elsasser number. The results obtained herein suggest that on planets with oblique magnetic fields, advective shifts of dayside hotspots may have substantial latitudinal components. This prediction may be tested observationally using the eclipse mapping technique.

Continued observational characterization of transiting planets that reside in close proximity to their host stars has shown that a substantial fraction of such objects posses orbits that are inclined with respect to the spin axes of their stars. Mounting evidence for the wide-spread nature of this phenomenon has challenged the conventional notion that large-scale orbital transport occurs during the early epochs of planet formation and is accomplished via planet-disk interactions. However, recent work has shown that the excitation of spin-orbit misalignment between protoplanetary nebulae and their host stars can naturally arise from gravitational perturbations in multi-stellar systems as well as magnetic disk-star coupling. In this work, we examine these processes in tandem. We begin with a thorough exploration of the gravitationally-facilitated acquisition of spin-orbit misalignment and analytically show that the entire possible range of misalignments can be trivially reproduced. Moreover, we demonstrate that the observable spin-orbit misalignment only depends on the primordial disk-binary orbit inclination. Subsequently, we augment our treatment by accounting for magnetic torques and show that more exotic dynamical evolution is possible, provided favorable conditions for magnetic tilting. Cumulatively, our results suggest that observed spin-orbit misalignments are fully consistent with disk-driven migration as a dominant mechanism for the origin of close-in planets.

For 61 transiting hot Jupiters, the projection of the angle between the orbital plane and the stellar equator (called the spin-orbit angle) has been measured. For about half of them, a significant misalignment is detected, and retrograde planets have been observed. This challenges scenarios of the formation of hot Jupiters. In order to better constrain formation models, we relate the distribution of the real spin-orbit angle $\Psi$ to the projected one $\beta$. Then, a comparison with the observations is relevant. We analyse the geometry of the problem to link analytically the projected angle $\beta$ to the real spin-orbit angle $\Psi$. The distribution of $\Psi$ expected in various models is taken from the literature, or derived with a simplified model and Monte-Carlo simulations in the case of the disk-torquing mechanism. An easy formula to compute the probability density function (PDF) of $\beta$ knowing the PDF of $\Psi$ is provided. All models tested here look compatible with the observed distribution beyond 40 degrees, which is so far poorly constrained by only 18 observations. But only the disk-torquing mechanism can account for the excess of aligned hot Jupiters, provided that the torquing is not always efficient. This is the case if the exciting binaries have semi-major axes as large as 10000 AU. Based on comparison with the set of observations available today, scattering models and the Kozai cycle with tidal friction models can not be solely responsible for the production of all hot Jupiters. Conversely, the presently observed distribution of the spin-orbit angles is compatible with most hot Jupiters having been transported by smooth migration inside a proto-planetary disk, itself possibly torqued by a companion.

The variation of a planet's obliquity is influenced by the existence of satellites with a high mass ratio. For instance, the Earth's obliquity is stabilized by the Moon, and would undergo chaotic variations in the Moon's absence. In turn, such variations can lead to large-scale changes in the atmospheric circulation, rendering spin-axis dynamics a central issue for understanding climate. The relevant quantity for dynamically-forced climate change is the rate of chaotic diffusion. Accordingly, here we reexamine the spin-axis evolution of a Moonless Earth within the context of a simplified perturbative framework. We present analytical estimates of the characteristic Lyapunov coefficient as well as the chaotic diffusion rate and demonstrate that even in absence of the Moon, the stochastic change in the Earth's obliquity is sufficiently slow to not preclude long-term habitability. Our calculations are consistent with published numerical experiments and illustrate the putative system's underlying dynamical structure in a simple and intuitive manner.

We present high spatial resolution images that demonstrate the hot Jupiters WASP-12b and HAT-P-8b orbit the primary star of hierarchical triple star systems. In each case, two distant companions with colors and brightness consistent with M dwarfs co-orbit the planet host as well as one another. Our adaptive optics images spatially resolve the secondary around WASP-12, previously identified by Bergfors et al. 2011 and Crossfield et al. 2012, into two distinct sources separated by 84.3+/-0.6 mas (21 +/- 3 AU). We find that the secondary to HAT-P-8, also identified by Bergfors et al. 2011, is in fact composed of two stars separated by 65.3+/-0.5 mas (15+/-1 AU). Our follow-up observations demonstrate physical association through common proper-motion. HAT-P-8 C has a particularly low mass, which we estimate to be 0.18+/-0.02Msun using photometry. Due to their hierarchy, WASP-12 BC and HAT-P-8 BC will enable the first dynamical mass determination for hot Jupiter stellar companions. These previously well-studied planet hosts now represent higher-order multi-star systems with potentially complex dynamics, underscoring the importance of diffraction-limited imaging and providing additional context for understanding the migrant population of transiting hot Jupiters.

In this paper we search for distant massive companions to known transiting gas giant planets that may have influenced the dynamical evolution of these systems. We present new radial velocity observations for a sample of 51 planets obtained using the Keck HIRES instrument, and find statistically significant accelerations in fifteen systems. Six of these systems have no previously reported accelerations in the published literature: HAT-P-10, HAT-P-22, HAT-P-29, HAT-P-32, WASP-10, and XO-2. We combine our radial velocity fits with Keck NIRC2 adaptive optics (AO) imaging data to place constraints on the allowed masses and orbital periods of the companions responsible for the detected accelerations. The estimated masses of the companions range between 1-500 M_Jup, with orbital semi-major axes typically between 1-75 AU. A significant majority of the companions detected by our survey are constrained to have minimum masses comparable to or larger than those of the transiting planets in these systems, making them candidates for influencing the orbital evolution of the inner gas giant. We estimate a total occurrence rate of 51 +/- 10% for companions with masses between 1-13 M_Jup and orbital semi-major axes between 1-20 AU in our sample. We find no statistically significant difference between the frequency of companions to transiting planets with misaligned or eccentric orbits and those with well-aligned, circular orbits. We combine our expanded sample of radial velocity measurements with constraints from transit and secondary eclipse observations to provide improved measurements of the physical and orbital characteristics of all of the planets included in our survey.

Here we measure the absolute magnitude distributions (H-distribution) of the dynamically excited and quiescent (hot and cold) Kuiper Belt objects (KBOs), and test if they share the same H-distribution as the Jupiter Trojans. From a compilation of all useable ecliptic surveys, we find that the KBO H-distributions are well described by broken power-laws. The cold population has a bright-end slope, $\alpha_{\textrm{1}}=1.5_{-0.2}^{+0.4}$, and break magnitude, $H_{\textrm{B}}=6.9_{-0.2}^{+0.1}$ (r'-band). The hot population has a shallower bright-end slope of, $\alpha_{\textrm{1}}=0.87_{-0.2}^{+0.07}$, and break magnitude $H_{\textrm{B}}=7.7_{-0.5}^{+1.0}$. Both populations share similar faint end slopes of $\alpha_2\sim0.2$. We estimate the masses of the hot and cold populations are $\sim0.01$ and $\sim3\times10^{-4} \mbox{ M$_{\bigoplus}$}$. The broken power-law fit to the Trojan H-distribution has $\alpha_\textrm{1}=1.0\pm0.2$, $\alpha_\textrm{2}=0.36\pm0.01$, and $H_{\textrm{B}}=8.3$. The KS test reveals that the probability that the Trojans and cold KBOs share the same parent H-distribution is less than 1 in 1000. When the bimodal albedo distribution of the hot objects is accounted for, there is no evidence that the H-distributions of the Trojans and hot KBOs differ. Our findings are in agreement with the predictions of the Nice model in terms of both mass and H-distribution of the hot and Trojan populations. Wide field survey data suggest that the brightest few hot objects, with $H_{\textrm{r'}}\lesssim3$, do not fall on the steep power-law slope of fainter hot objects. Under the standard hierarchical model of planetesimal formation, it is difficult to account for the similar break diameters of the hot and cold populations given the low mass of the cold belt.

The presence of giant gaseous planets that reside in close proximity to their host stars may be a consequence of large-scale radial migration through the proto-planetary nebulae. Within the context of this picture, significant orbital obliquities characteristic of a substantial fraction of such planets can be attributed to external torques that perturb the disks out of alignment with the spin axes of their host stars. Therefore, the acquisition of orbital obliquity exhibits sensitive dependence on the physics of disk-star interactions. Here, we analyze the primordial excitation of spin-orbit misalignment of Sun-like stars, in light of disk-star angular momentum transfer. We begin by calculating the stellar pre-main sequence rotational evolution, accounting for spin-up due to gravitational contraction and accretion as well as spin-down due to magnetic star-disk coupling. We devote particular attention to angular momentum transfer by accretion, and show that while generally subdominant to gravitational contraction, this process is largely controlled by the morphology of the stellar magnetic field (i.e. specific angular momentum accreted by stars with octupole-dominated surface fields is smaller than that accreted by dipole-dominated stars by an order of magnitude). Subsequently, we examine the secular spin-axis dynamics of disk-bearing stars, accounting for the time-evolution of stellar and disk properties and demonstrate that misalignments are preferentially excited in systems where stellar rotation is not overwhelmingly rapid. Moreover, we show that the excitation of spin-orbit misalignment occurs impulsively, through an encounter with a resonance between the stellar precession frequency and the disk-torquing frequency. Cumulatively, the model developed herein opens up a previously unexplored avenue towards understanding star-disk evolution and its consequences in a unified manner.

Giant gaseous planets often reside on orbits in sufficient proximity to their host stars for the planetary quadrupole gravitational field to become non-negligible. In presence of an additional planetary companion, a precise characterization of the system's orbital state can yield meaningful constraints on the transiting planet's interior structure. However, such methods can require a very specific type of system. This paper explores the dynamic range of applicability of these methods and shows that interior structure calculations are possible for a wide array of orbital architectures. The HAT-P-13 system is used as a case study, and the implications of perturbations arising from a third distant companion on the feasibility of an interior calculation are discussed. We find that the method discussed here is likely to be useful in studying other planetary systems, allowing the possibility of an expanded survey of the interiors of exoplanets.

Through the process of thermal ionization, intense stellar irradiation renders Hot Jupiter atmospheres electrically conductive. Simultaneously, lateral variability in the irradiation drives the global circulation with peak wind speeds of order ~ km/s. In turn, the interactions between the atmospheric flows and the background magnetic field give rise to Lorentz forces that can act to perturb the flow away from its purely hydrodynamical counterpart. Using analytical theory and numerical simulations, here we show that significant deviations away from axisymmetric circulation are unstable in presence of a non-negligible axisymmetric magnetic field. Specifically, our results suggest that dayside-to-nightside flows, often obtained within the context of three-dimensional circulation models, only exist on objects with anomalously low magnetic fields, while the majority of highly irradiated exoplanetary atmospheres are entirely dominated by zonal jets.

An ever-growing observational aggregate of extrasolar planets has revealed that systems of planets that reside in or near mean-motion resonances are relatively common. While the origin of such systems is attributed to protoplanetary disk-driven migration, a qualitative description of the dynamical evolution of resonant planets remains largely elusive. Aided by the pioneering works of the last century, we formulate an approximate, integrable theory for first-order resonant motion. We utilize the developed theory to construct an intuitive, geometrical representation of resonances within the context of the unrestricted three-body problem. Moreover, we derive a simple analytical criterion for the appearance of secondary resonances between resonant and secular motion. Subsequently, we demonstrate the onset of rapid chaotic motion as a result of overlap among neighboring first-order mean-motion resonances, as well as the appearance of slow chaos as a result of secular modulation of the planetary orbits. Finally, we take advantage of the integrable theory to analytically show that, in the adiabatic regime, divergent encounters with first-order mean-motion resonances always lead to persistent apsidal anti-alignment.

Recently, the Kepler spacecraft has detected a sizable aggregate of objects, characterized by giant-planet-like radii and modest levels of stellar irradiation. With the exception of a handful of objects, the physical nature, and specifically the average densities, of these bodies remain unknown. Here, we propose that the detected giant planet radii may partially belong to planets somewhat less massive than Uranus and Neptune. Accordingly, in this work, we seek to identify a physically sound upper limit to planetary radii at low masses and moderate equilibrium temperatures. As a guiding example, we analyze the interior structure of the Neptune-mass planet Kepler-30d and show that it is acutely deficient in heavy elements, especially compared with its solar system counterparts. Subsequently, we perform numerical simulations of planetary thermal evolution and in agreement with previous studies, show that generally, 10 - 20 Earth mass, multi-billion year old planets, composed of high density cores and extended H/He envelopes can have radii that firmly reside in the giant planet range. We subject our results to stability criteria based on extreme ultraviolet radiation, as well as Roche-lobe overflow driven mass-loss and construct mass-radius relationships for the considered objects. We conclude by discussing observational avenues that may be used to confirm or repudiate the existence of putative low mass, gas-dominated planets.

The presence of gaseous giant planets whose orbits lie in extreme proximity to their host stars ("hot Jupiters"), can largely be accounted for by planetary migration, associated with viscous evolution of proto-planetary nebulae. Recently, observations of the Rossiter-McLaughlin effect during planetary transits have revealed that a considerable fraction of detected hot Jupiters reside on orbits that are misaligned with respect to the spin-axes of their host stars. This observational fact has cast significant doubts on the importance of disk-driven migration as a mechanism for production of hot Jupiters, thereby reestablishing the origins of close-in planetary orbits as an open question. Here we show that misaligned orbits can be a natural consequence of disk migration. Our argument rests on an enhanced abundance of binary stellar companions in star formation environments, whose orbital plane is uncorrelated with the spin axes of the individual stars. We analyze the dynamical evolution of idealized proto-planetary disks under perturbations from massive distant bodies and demonstrate that the resulting gravitational torques act to misalign the orbital planes of the disks relative to the spin poles of their host stars. As a result, we predict that in the absence of strong disk-host star angular momentum coupling or sufficient dissipation that acts to realign the stellar spin axis and the planetary orbits, the fraction of planetary systems (including systems of hot Neptunes and Super-Earths), whose angular momentum vectors are misaligned with respect to their host-stars should be commensurate with the rate of primordial stellar multiplicity.

Here we present new adaptive optics observations of the Quaoar-Weywot system. With these new observations we determine an improved system orbit. Due to a 0.39 day alias that exists in available observations, four possible orbital solutions are available with periods of $\sim11.6$, $\sim12.0$, $\sim12.4$, and $\sim12.8$ days. From the possible orbital solutions, system masses of $1.3-1.5\pm0.1\times10^{21}$ kg are found. These observations provide an updated density for Quaoar of $2.7-5.0{g cm$^{-3}$}$. In all cases, Weywot's orbit is eccentric, with possible values $\sim0.13-0.16$. We present a reanalysis of the tidal orbital evolution of the Quoaor-Weywot system. We have found that Weywot has probably evolved to a state of synchronous rotation, and have likely preserved their initial inclinations over the age of the Solar system. We find that for plausible values of the effective tidal dissipation factor tides produce a very slow evolution of Weywot's eccentricity and semi-major axis. Accordingly, it appears that Weywot's eccentricity likely did not tidally evolve to its current value from an initially circular orbit. Rather, it seems that some other mechanism has raised its eccentricity post-formation, or Weywot formed with a non-negligible eccentricity.

We show that the existence of prograde equatorial satellites is consistent with a collisional tilting scenario for Uranus. In fact, if the planet was surrounded by a proto-satellite disk at the time of the tilting and a massive ring of material was temporarily placed inside the Roche radius of the planet by the collision, the proto-satellite disk would have started to precess incoherently around the equator of the planet, up to a distance greater than that of Oberon. Collisional damping would then have collapsed it into a thin equatorial disk, from which the satellites eventually formed. The fact that the orbits of the satellites are prograde requires Uranus to have had a non-negligible initial obliquity (comparable to that of Neptune) before it was finally tilted to 98 degrees.

A considerable fraction of multi-planet systems discovered by the observational surveys of extrasolar planets reside in mild proximity to first-order mean motion resonances. However, the relative remoteness of such systems from nominal resonant period ratios (e.g. 2:1, 3:2, 4:3) has been interpreted as evidence for lack of resonant interactions. Here we show that a slow divergence away from exact commensurability is a natural outcome of dissipative evolution and demonstrate that libration of critical angles can be maintained tens of percent away from nominal resonance. We construct an analytical theory for the long-term dynamical evolution of dissipated resonant planetary pairs and confirm our calculations numerically. Collectively, our results suggest that a significant fraction of the near-commensurate extrasolar planets are in fact resonant and have undergone significant dissipative evolution.

We present calculations of thermal evolution of Hot Jupiters with various masses and effective temperatures under Ohmic dissipation. The resulting evolutionary sequences show a clear tendency towards inflated radii for effective temperatures that give rise to significant ionization of alkali metals in the atmosphere, compatible with the trend of the data. The degree of inflation shows that Ohmic dissipation, along with the likely variability in heavy element content can account for all of the currently detected radius anomalies. Furthermore, we find that in absence of a massive core, low-mass hot Jupiters can over-flow their Roche-lobes and evaporate on Gyr time-scales, possibly leaving behind small rocky cores.

Over the last decade, evidence has mounted that the solar system's observed state can be favorably reproduced in the context of an instability-driven dynamical evolution model, such as the "Nice" model. To date, all successful realizations of instability models have concentrated on evolving the four giant planets onto their current orbits from a more compact configuration. Simultaneously, the possibility of forming and ejecting additional planets has been discussed, but never successfully implemented. Here we show that a large array of 5-planet (2 gas giants + 3 ice giants) multi-resonant initial states can lead to an adequate formation of the outer solar system, featuring an ejection of an ice giant during a phase of instability. Particularly, our simulations demonstrate that the eigenmodes which characterize the outer solar system's secular dynamics can be closely matched with a 5-planet model. Furthermore, provided that the ejection timescale of the extra planet is short, orbital excitation of a primordial cold classical Kuiper belt can also be avoided in this scenario. Thus the solar system is one of many possible outcomes of dynamical relaxation and can originate from a wide variety of initial states. This deems the construction of a unique model of solar system's early dynamical evolution impossible.

The cold classical population of the Kuiper belt exhibits a wide variety of unique physical characteristics, which collectively suggest that its dynamical coherence has been maintained through out the solar system's lifetime. Simultaneously, the retention of the cold population's relatively unexcited orbital state has remained a mystery, especially in the context of a solar system formation model, that is driven by a transient period of instability, where Neptune is temporarily eccentric. Here, we show that the cold belt can survive the instability, and its dynamical structure can be reproduced. We develop a simple analytical model for secular excitation of cold KBOs and show that comparatively fast apsidal precession and nodal recession of Neptune, during the eccentric phase, are essential for preservation of an unexcited state in the cold classical region. Subsequently, we confirm our results with self-consistent N-body simulations. We further show that contamination of the hot classical and scattered populations by objects of similar nature to that of cold classicals has been instrumental in shaping the vast physical diversity inherent to the Kuiper belt.

The presence of highly eccentric extrasolar planets in binary stellar systems suggests that the Kozai effect has played an important role in shaping their dynamical architectures. However, the formation of planets in inclined binary systems poses a considerable theoretical challenge, as orbital excitation due to the Kozai resonance implies destructive, high-velocity collisions among planetesimals. To resolve the apparent difficulties posed by Kozai resonance, we seek to identify the primary physical processes responsible for inhibiting the action of Kozai cycles in protoplanetary disks. Subsequently, we seek to understand how newly-formed planetary systems transition to their observed, Kozai-dominated dynamical states. We find that theoretical difficulties in planet formation arising from the presence of a distant companion star, posed by the Kozai effect and other secular perturbations, can be overcome by a proper account of gravitational interactions within the protoplanetary disk. In particular, fast apsidal recession induced by disk self-gravity tends to erase the Kozai effect, and ensure that the disk's unwarped, rigid structure is maintained. Subsequently, once a planetary system has formed, the Kozai effect can continue to be wiped out as a result of apsidal precession, arising from planet-planet interactions. However, if such a system undergoes a dynamical instability, its architecture may change in such a way that the Kozai effect becomes operative. The results presented here suggest that planetary formation in highly inclined binary systems is not stalled by perturbations, arising from the stellar companion. Consequently, planet formation in binary stars is probably no different from that around single stars on a qualitative level. Furthermore, it is likely that systems where the Kozai effect operates, underwent a transient phase of dynamical instability in the past.

As a result of resonance overlap, planetary systems can exhibit chaotic motion. Planetary chaos has been studied extensively in the Hamiltonian framework, however, the presence of chaotic motion in systems where dissipative effects are important, has not been thoroughly investigated. Here, we study the onset of stochastic motion in presence of dissipation, in the context of classical perturbation theory, and show that planetary systems approach chaos via a period-doubling route as dissipation is gradually reduced. Furthermore, we demonstrate that chaotic strange attractors can exist in mildly damped systems. The results presented here are of interest for understanding the early dynamical evolution of chaotic planetary systems, as they may have transitioned to chaos from a quasi-periodic state, dominated by dissipative interactions with the birth nebula.

Long-term orbital evolution of multi-planet systems under tidal dissipation often converges to a stationary state, known as the tidal fixed point. The fixed point is characterized by a lack of oscillations in the eccentricities and apsidal alignment among the orbits. Quantitatively, the nature of the fixed point is dictated by mutual interactions among the planets as well as non-Keplerian effects. We show that if a roughly coplanar system hosts a hot, sub-Saturn mass planet, and is tidally relaxed, separation of planet-planet interactions and non-Keplerian effects in the equations of motion leads to a direct determination of the true masses of the planets. Consequently, a "snap-shot" observational determination of the orbital state resolves the sin(I) degeneracy, and opens up a direct avenue towards identification of the true lowest-mass exo-planets detected. We present an approximate, as well as a general, mathematical framework for computation of the line of sight inclination of secular systems, and apply our models illustratively to the 61 Vir system. We conclude by discussing the observability of planetary systems to which our method is applicable and we set our analysis into a broader context by presenting a current summary of the various possibilities for determining the physical properties of planets from observations of their orbital states.

In the recent years, the "Nice" model of solar system formation has attained an unprecedented level of success in reproducing much of the observed orbital architecture of the solar system by evolving the planets to their current locations from a more compact configuration. Within the context of this model, the formation of the classical Kuiper belt requires a phase during which the ice giants have a high eccentricity. An outstanding question of this model is the initial configuration from which the Solar System started out. Recent work has shown that multi-resonant initial conditions can serve as good candidates, as they naturally prevent vigorous type-II migration. In this paper, we use analytical arguments, as well as self-consistent numerical N-body simulations to identify fully-resonant initial conditions, whose dynamical evolution is characterized by an eccentric phase of the ice-giants, as well as planetary scattering. We find a total of eight such initial conditions. Four of these primordial states are compatible with the canonical "Nice" model, while the others imply slightly different evolutions. The results presented here should prove useful in further development of a comprehensive model for solar system formation.

We present a new, magnetohydrodynamic mechanism for inflation of close-in giant extrasolar planets. The idea behind the mechanism is that current, which is induced through interaction of atmospheric winds and the planetary magnetic field, results in significant Ohmic dissipation of energy in the interior. We develop an analytical model for computation of interior Ohmic dissipation, with a simplified treatment of the atmosphere. We apply our model to HD209458b, Tres-4b and HD189733b. With conservative assumptions for wind speed and field strength, our model predicts a generated power that appears to be large enough to maintain the transit radii, opening an unexplored avenue towards solving a decade-old puzzle of extrasolar gas giant radius anomalies.

We present precision radial velocity (RV) data that reveal a multiple exoplanet system orbiting the bright nearby G5V star 61 Virginis. Our 4.6 years of combined Keck/HIRES and Anglo-Australian Telescope precision RVs indicate the hitherto unknown presence of at least three planets orbiting this well-studied star. These planets are all on low-eccentricity orbits with periods of 4.2, 38.0, and 124.0 days, and projected masses (M sini) of 5.1, 18.2, and 24.0 M_Earth, respectively. Test integrations of systems consistent with the RV data suggest that the configuration is dynamically stable. Depending on the effectiveness of tidal dissipation within the inner planet, the inner two planets may have evolved into an eccentricity fixed-point configuration in which the apsidal lines of all three planets corotate. This conjecture can be tested with additional observations. We present a 16-year time series of photometric observations of 61 Virginis, which comprise 1194 individual measurements, and indicate that it has excellent photometric stability. No significant photometric variations at the periods of the proposed planets have been detected. This new system is the first known example of a G-type Sun-like star hosting a Super-Earth mass planet. It joins HD 75732 (55 Cnc), HD 69830, GJ 581, HD 40307, and GJ 876 as a growing group of exoplanet systems that have multiple planets orbiting with periods less than an Earth-year. The ubiquity of such systems portends that space-based transit-search missions such as KEPLER and COROT will find many multi-transiting systems.

Tidal dissipation within a short-period transiting extrasolar planet perturbed by a companion object can drive orbital evolution of the system to a so-called tidal fixed point, in which the apsidal lines of the transiting planet and its perturber are aligned, and for which variations in the orbital eccentricities of both planet and perturber are damped out. Significant contributions to the apsidal precession rate are made by the secular planet-planet interaction, by general relativity, and by the gravitational quadropole fields created by the transiting planet's tidal and rotational distortions. The fixed-point orbital eccentricity of the inner planet is therefore a strong function of the planet's interior structure. We illustrate these ideas in the specific context of the recently discovered HAT-P-13 exo-planetary system, and show that one can already glean important insights into the physical properties of the inner transiting planet. We present structural models of the planet, which indicate that its observed radius can be maintained for a one-parameter sequence of models that properly vary core mass and tidal energy dissipation in the interior. We use an octopole-order secular theory of the orbital dynamics to derive the dependence of the inner planet's eccentricity, on its tidal Love number. We find that the currently measured eccentricity, implies 0.116 < k2_{b} < 0.425, 0 M_{Earth}<M_{core}<120 M_{Earth}$, and Q_{b} < 300,000. Improved measurement of the eccentricity will soon allow for far tighter limits to be placed on all three of these quantities, and will provide an unprecedented probe into the interior structure of an extrasolar planet.

We investigate the possibility that the large orbital eccentricity of the transiting Neptune-mass planet Gliese 436b is maintained in the face of tidal dissipation by a second planet in the system. We find that the currently observed configuration can be understood if Gliese 436b and a putative companion have evolved to a quasi-stationary fixed point in which the planets' orbital apses are co-linear and in which secular variations in the orbital eccentricities of the two planets have been almost entirely damped out. We adopt an octopole-order secular theory based on a Legendre expansion in the semi-major axis ratio to delineate well-defined regions of (P_c, M_c, e_c) space that can be occupied by a perturbing companion. We incorporate the evolutionary effect of tidal dissipation into our secular model of the system, and solve the resulting initial value problems for a large sample of the allowed configurations. We then polish the stationary configurations derived from secular theory with full numerical integrations. We present our results in the form of candidate companion planets to Gliese 436b. For these candidates, radial velocity half-amplitudes, K_c, are of order 3 m/s, and the maximum amplitude of orbit-to-orbit transit timing variations are of order Delta t=1 s to Delta t=5s. For the particular example case of a perturber with orbital period, P_c=40 d, mass, M_c=8.5 M_Earth, and eccentricity, e_c=0.58, we confirm our semi-analytic calculations with a full numerical 3-body integration of the orbital decay that includes tidal damping and spin evolution. Additionally, we discuss the possibility of many-perturber stationary configurations, utilizing modified Laplace-Lagrange secular theory.

A long-term numerical integration of the classical Newtonian approximation to the planetary orbital motions of the full Solar System (sun + 8 planets), spanning 20 Gyr, was performed. The results showed no severe instability arising over this time interval. Subsequently, utilizing a bifurcation method described by Jacques Laskar, two numerical experiments were performed with the goal of determining dynamically allowed evolutions for the Solar System in which the planetary orbits become unstable. The experiments yielded one evolution in which Mercury falls onto the Sun at ~1.261Gyr from now, and another in which Mercury and Venus collide in ~862Myr. In the latter solution, as a result of Mercury's unstable behavior, Mars was ejected from the Solar System at ~822Myr. We have performed a number of numerical tests that confirm these results, and indicate that they are not numerical artifacts. Using synthetic secular perturbation theory, we find that Mercury is destabilized via an entrance into a linear secular resonance with Jupiter in which their corresponding eigenfrequencies experience extended periods of commensurability. The effects of general relativity on the dynamical stability are discussed. An application of the bifurcation method to the outer Solar System (Jupiter, Saturn, Uranus, and Neptune) showed no sign of instability during the course of 24Gyr of integrations, in keeping with an expected Uranian dynamical lifetime of 10^(18) years.