Source author record

Konstantin Batygin

Konstantin Batygin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.EP astro-ph.SR math.DS nlin.CD astro-ph.GA

Catalog footprint

What is connected

72works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

A young progenitor for the most common planetary systems in the Galaxy

The Galaxy's most common known planetary systems have several Earth-to-Neptune-size planets in compact orbits. At small orbital separations, larger planets are less common than their smaller counterparts by an order of magnitude. The young star V1298 Tau hosts one such compact planetary system, albeit with four planets that are uncommonly large (5 to 10 Earth radii). The planets form a chain of near-resonances that result in transit-timing variations of several hours. Here we present a multi-year campaign to characterize this system with transit-timing variations, a method insensitive to the intense magnetic activity of the star. Through targeted observations, we first resolved the previously unknown orbital period of the outermost planet. The full 9-year baseline from these and archival data then enabled robust determination of the masses and orbital parameters for all four planets. We find the planets have low, sub-Neptune masses and nearly circular orbits, implying a dynamically tranquil history. Their low masses and large radii indicate that the inner planets underwent a period of rapid cooling immediately after dispersal of the protoplanetary disk. Still, they are much less dense than mature planets of comparable size. We predict the planets will contract to 1.5-4.0 Earth radii and join the population of super-Earths and sub-Neptunes that nature produces in abundance.

preprint2026arXiv

Callisto's Nonresonant Orbit as an Outcome of Circum-Jovian Disk Substructure

The Galilean moons of Io, Europa, and Ganymede exhibit a 4:2:1 commensurability in their mean motions, a configuration known as the Laplace resonance. The prevailing view for the origin of this three-body resonance involves the convergent migration of the moons, resulting from gas-driven torques in the circum-Jovian disk wherein they accreted. To account for Callisto's exclusion from the resonant chain, a late and/or slow accretion of the fourth and outermost Galilean moon is typically invoked, stalling its migration. Here, we consider an alternative scenario in which Callisto's nonresonant orbit is a consequence of disk substructure. Using a suite of N-body simulations that self-consistently account for satellite-disk interactions, we show that a pressure bump can function as a migration trap, isolating Callisto and alleviating constraints on its timing of accretion. Our simulations position the bump interior to the birthplaces of all four moons. In exploring the impact of bump structure on simulation outcomes, we find that it cannot be too sharp nor flat to yield the observed orbital architecture. In particular, a "Goldilocks" zone is mapped in parameter space, corresponding to a well-defined range in bump aspect ratio. Within this range, Io, Europa, and Ganymede are sequentially trapped at the bump, and ushered across it through resonant lockstep migration with their neighboring, exterior moon. The implications of our work are discussed in the context of uncertainties regarding Callisto's interior structure, arising from the possibility of non-hydrostatic contributions to its shape and gravity field, unresolved by the Galileo spacecraft.

preprint2026arXiv

Secular Excitation of Polar Neptune Orbits in Pure Disk-Planet Systems

The stellar spin-orbit angles of Neptune-sized planets present a primordial yet puzzling view of the planetary formation epoch. The striking dichotomy of aligned and perpendicular orbital configurations are suggestive of obliquity excitation through secular resonance -- a process where the precession of a hot Neptune becomes locked onto a forcing frequency, and is slowly guided into a perpendicular state. Previous models of resonant capture have involved the presence of companion perturbers to the star-planet-disk system, but in most cases, such companions are not confirmed to be present. In this work, we present a mechanism for exciting Neptunes to polar orbits in systems without giant perturbers, where photo-evaporation is the self-contained mechanism. Photo-evaporation opens a gap in the protoplanetary disk at ~1 au, and the inner disk continues to viscously accrete onto the host star, precessing quickly due to the perturbation of the outer disk. As the inner disk shrinks, it precesses more slowly, and encounters a resonance with the J2 precession of the Neptune, quickly exciting it to a polar configuration. While likely not applicable to more massive planets which trigger back-reactions onto the disk, this mechanism reproduces the obliquities of small planets in multiple respects.

preprint2023arXiv

Formation of Rocky Super-Earths From A Narrow Ring of Planetesimals

The formation of super-Earths, the most abundant planets in the Galaxy, remains elusive. These planets have masses that typically exceed that of the Earth by a factor of a few; appear to be predominantly rocky, although often surrounded by H/He atmospheres; and frequently occur in multiples. Moreover, planets that encircle the same star tend to have similar masses and radii, whereas those belonging to different systems exhibit remarkable overall diversity. Here, we advance a theoretical picture for rocky planet formation that satisfies the aforementioned constraints: building upon recent work - which demonstrates that planetesimals can form rapidly at discrete locations in the disk - we propose that super-Earths originate inside rings of silicate-rich planetesimals at approximately ~1 AU. Within the context of this picture, we show that planets grow primarily through pairwise collisions among rocky planetesimals, until they achieve terminal masses that are regulated by isolation and orbital migration. We quantify our model with numerical simulations and demonstrate that our synthetic planetary systems bear a close resemblance to compact, multi-resonant progenitors of the observed population of short-period extrasolar planets. Our results thus indicate that the absence of short-period super-Earths within the solar system can simply be attributed to the comparatively low mass of the primordial planetesimal ring within the protosolar nebula.

preprint2022arXiv

A Criterion for the Stability of Planets in Chains of Resonances

Uncovering the formation process that reproduces the distinct properties of compact super-Earth exoplanet systems is a major goal of planet formation theory. The most successful model argues that non-resonant systems begin as resonant chains of planets that later experience a dynamical instability. However, both the boundary of stability in resonant chains and the mechanism of the instability itself are poorly understood. Previous work postulated that a secondary resonance between the fastest libration frequency and a difference in synodic frequencies destabilizes the system. Here, we use that hypothesis to produce a simple and general criterion for resonant chain stability that depends only on planet orbital periods and masses. We show that the criterion accurately predicts the maximum mass of planets in synthetic resonant chains up to six planets. More complicated resonant chains produced in population synthesis simulations are found to be less stable than expected, although our criterion remains useful and superior to machine learning models.

preprint2022arXiv

A Multi-Planet System's Sole Super-Puff: Exploring Allowable Physical Parameters for the Cold Super-Puff HIP 41378 f

The census of known exoplanets exhibits a variety of physical parameters, including densities that are measured to span the range from less dense than styrofoam to more dense than iron. These densities represent a large diversity of interior structures. Despite this staggering diversity, recent analyses have shown that the densities of planets that orbit a common star exhibit remarkable uniformity. A fascinating exception to this is the system HIP 41378 (also known as K2-93), which contains a super-puff planet, HIP 41378 f, as well as several planets with more typical bulk densities. The range of densities in this system begs the question of what physical processes are responsible for the disparate planetary structures in this system. In this paper, we consider how the densities of the planets in the HIP 41378 system would have changed over time as the host star evolved and the planets' atmospheres were subsequently affected by the evolving insolation level. We also present a range of allowable core masses for HIP 41378 f based on the measured planet parameters, and comment on the feasibility of the proposed existence of planetary rings around HIP 41378 f as an explanation for its current low density.

preprint2022arXiv

Analytic Approach to the Late Stages of Giant Planet Formation

This paper constructs an analytic description for the late stages of giant planet formation. During this phase of evolution, the planet gains the majority of its final mass through gas accretion at a rapid rate. This work determines the density and velocity fields for material falling onto the central planet and its circumplanetary disk, and finds the corresponding column density of this infalling envelope. We derive a steady-state solution for the surface density of the disk as a function of its viscosity (including the limiting case where no disk accretion occurs). Planetary magnetic fields truncate the inner edge of the disk and determine the boundary conditions for mass accretion onto the planet from both direct infall and from the disk. The properties of the forming planet and its circumplanetary disk are determined, including the luminosity contributions from infall onto the planet and disk surfaces, and from disk viscosity. The radiative signature of the planet formation process is explored using a quasi-spherical treatment of the emergent spectral energy distributions. The analytic solutions developed herein show how the protoplanet properties (envelope density distribution, velocity field, column density, disk surface density, luminosity, and radiative signatures) vary with input parameters (instantaneous mass, orbital location, accretion rate, and planetary magnetic field strength).

preprint2022arXiv

Architectures of Compact Super-Earth Systems Shaped by Instabilities

Compact non-resonant systems of sub-Jovian planets are the most common outcome of the planet formation process. Despite exhibiting broad overall diversity, these planets also display dramatic signatures of intra-system uniformity in their masses, radii, and orbital spacings. Although the details of their formation and early evolution are poorly known, sub-Jovian planets are expected to emerge from their natal nebulae as multi-resonant chains, owing to planet-disk interactions. Within the context of this scenario, the architectures of observed exoplanet systems can be broadly replicated if resonances are disrupted through post-nebular dynamical instabilities. Here, we generate an ad-hoc sample of resonant chains and use a suite of N-body simulations to show that instabilities can not only reproduce the observed period ratio distribution, but that the resulting collisions also modify the mass uniformity in a way that is consistent with the data. Furthermore, we demonstrate that primordial mass uniformity, motivated by the sample of resonant chains coupled with dynamical sculpting, naturally generates uniformity in orbital period spacing similar to what is observed. Finally, we find that almost all collisions lead to perfect mergers, but some form of post-instability damping is likely needed to fully account for the present-day dynamically cold architectures of sub-Jovian exoplanets.

preprint2022arXiv

Generating Stellar Obliquity in Systems with Broken Protoplanetary Disks

Recent advances in sub-millimeter observations of young circumstellar nebulae have opened an unprecedented window into the structure of protoplanetary disks, which has revealed the surprising ubiquity of broken and misaligned disks. In this work, we demonstrate that such disks are capable of torquing the spin axis of their host star, representing a hitherto unexplored pathway by which stellar obliquities may be generated. The basis of this mechanism is a crossing of the stellar spin precession and inner disk regression frequencies, resulting in adiabatic excitation of the stellar obliquity. We derive analytical expressions for the characteristic frequencies of the inner disk and star as a function of the disk gap boundaries, and place an approximate limit on the disk architectures for which frequency crossing and resulting obliquity excitation are expected, thereby illustrating the efficacy of this model. Cumulatively, our results support the emerging concensus that significant spin-orbit misalignments are an expected outcome of planet formation.

preprint2022arXiv

Shallowness of circulation in hot Jupiters -- Advancing the Ohmic dissipation model

The inflated radii of giant short-period extrasolar planets collectively indicate that the interiors of hot Jupiters are heated by some anomalous energy dissipation mechanism. Although a variety of physical processes have been proposed to explain this heating, recent statistical evidence points to the confirmation of explicit predictions of the Ohmic dissipation theory, elevating this mechanism as the most promising candidate for resolving the radius inflation problem. In this work, we present an analytic model for the dissipation rate and derive a simple scaling law that links the magnitude of energy dissipation to the thickness of the atmospheric weather layer. From this relation, we find that the penetration depth influences the Ohmic dissipation rate by an order of magnitude. We further investigate the weather layer depth of hot Jupiters from the extent of their inflation and show that, depending on the magnetic field strength, hot Jupiter radii can be maintained even if the circulation layer is relatively shallow. Additionally, we explore the evolution of zonal wind velocities with equilibrium temperature by matching our analytic model to statistically expected dissipation rates. From this analysis, we deduce that the wind speed scales approximately as $1/\sqrt{T_\mathrm{eq}-T_0}$, where $T_0$ is a constant that equals $T_0 \sim 1000~\mathrm{K}-1800~\mathrm{K}$ depending on planet-specific parameters (radius, mass, etc.). This work outlines inter-related constraints on the atmospheric flow and the magnetic field of hot Jupiters and provides a foundation for future work on the Ohmic heating mechanism.

preprint2022arXiv

TESS-Keck Survey IX: Masses of Three Sub-Neptunes Orbiting HD 191939 and the Discovery of a Warm Jovian Plus a Distant Sub-Stellar Companion

Exoplanet systems with multiple transiting planets are natural laboratories for testing planetary astrophysics. One such system is HD 191939 (TOI-1339), a bright (V=9) and Sun-like (G9V) star, which TESS found to host three transiting planets (b, c, and d). The planets have periods of 9, 29, and 38 days each with similar sizes from 3 to 3.4 $R_{\oplus}$. To further characterize the system, we measured the radial velocity (RV) of HD 191939 over 415 days with Keck/HIRES and APF/Levy. We find that $M_b = 10.4 \pm 0.9 M_{\oplus}$ and $M_c = 7.2 \pm 1.4 M_{\oplus}$, which are low compared to most known planets of comparable radii. The RVs yield only an upper-limit on $M_d$ (<5.8 $M_{\oplus}$ at 2$σ$). The RVs further reveal a fourth planet (e) with a minimum mass of $0.34 \pm 0.01 M_{Jup}$ and an orbital period of 101.4 $\pm$ 0.4 days. Despite its non-transiting geometry, secular interactions between planet e and the inner transiting planets indicate that planet e is coplanar with the transiting planets ($Δ$i < 10$^{\circ}$). We identify a second high mass planet (f) with 95% confidence intervals on mass between $2-11 \, M_{Jup}$ and period between 1700-7200 days, based on a joint analysis of RVs and astrometry from $Gaia$ and $Hipparcos$. As a bright star hosting multiple planets with well-measured masses, HD 191939 presents many options for comparative planetary astronomy including characterization with JWST.

preprint2022arXiv

TOI-1136 is a Young, Coplanar, Aligned Planetary System in a Pristine Resonant Chain

Convergent disk migration has long been suspected to be responsible for forming planetary systems with a chain of mean-motion resonances (MMR). Dynamical evolution over time could disrupt the delicate resonant configuration. We present TOI-1136, a 700-Myr-old G star hosting at least 6 transiting planets between $\sim$2 and 5 $R_\oplus$. The orbital period ratios deviate from exact commensurability by only $10^{-4}$, smaller than the $\sim$\,$10^{-2}$ deviations seen in typical Kepler near-resonant systems. A transit-timing analysis measured the masses of the planets (3-8$M_\oplus$) and demonstrated that the planets in TOI-1136 are in true resonances with librating resonant angles. Based on a Rossiter-McLaughlin measurement of planet d, the star's rotation appears to be aligned with the planetary orbital planes. The well-aligned planetary system and the lack of detected binary companion together suggest that TOI-1136's resonant chain formed in an isolated, quiescent disk with no stellar fly-by, disk warp, or significant axial asymmetry. With period ratios near 3:2, 2:1, 3:2, 7:5, and 3:2, TOI-1136 is the first known resonant chain involving a second-order MMR (7:5) between two first-order MMR. The formation of the delicate 7:5 resonance places strong constraints on the system's migration history. Short-scale (starting from $\sim$0.1 AU) Type-I migration with an inner disk edge is most consistent with the formation of TOI-1136. A low disk surface density ($Σ_{\rm 1AU}\lesssim10^3$g~cm$^{-2}$; lower than the minimum-mass solar nebula) and the resultant slower migration rate likely facilitated the formation of the 7:5 second-order MMR. TOI-1136's deep resonance suggests that it has not undergone much resonant repulsion during its 700-Myr lifetime. One can rule out rapid tidal dissipation within a rocky planet b or obliquity tides within the largest planets d and f.

preprint2021arXiv

A search for Planet Nine using the Zwicky Transient Facility public archive

Recent estimates of the characteristics of Planet Nine have suggested that it could be closer than originally assumed. Such a Planet Nine would also be brighter than originally assumed, suggesting the possibility that it has already been observed in wide-field moderate-depth surveys. We search for Planet Nine in the Zwicky Transient Facility public archive and find no candidates. Using known asteroids to calculate the magnitude limit of the survey, we find that we should have detected Planet Nine throughout most of the northern portion of its predicted orbit -- including within the galactic plane -- to a 95% detection efficiency of approximately $V=20.5$. To aid in understanding detection limits for this and future analyses, we present a full-sky synthetic Planet Nine population drawn from a statistical sampling of predicted Planet Nine orbits. We use this reference population to estimate that this survey rules out 56% of predicted Planet Nine phase space, and we demonstrate how future analyses can use the same synthetic population to continue to constrain the amount of parameter space effectively searched for Planet Nine.

preprint2021arXiv

Contemporary formation of early solar system planetesimals at two distinct radial locations

The formation of planetesimals is expected to occur via particle-gas instabilities that concentrate dust into self-gravitating clumps. Triggering these instabilities requires the prior pileup of dust in the protoplanetary disk. Until now, this has been successfully modeled exclusively at the disk's snowline, whereas rocky planetesimals in the inner disk were obtained only by assuming either unrealistically large particle sizes or an enhanced global disk metallicity. However, planetesimal formation solely at the snowline is difficult to reconcile with the early and contemporaneous formation of iron meteorite parent bodies with distinct oxidation states and isotopic compositions, indicating formation at different radial locations in the disk. Here, by modeling the evolution of a disk with ongoing accretion of material from the collapsing molecular cloud, we show that planetesimal formation may have been triggered within the first 0.5 million years by dust pileup at both the snowline (at approximately 5 au) and the silicate sublimation line (at approximately 1 au), provided turbulent diffusion was low. Particle concentration at approximately 1 au is due to the early outward radial motion of gas and is assisted by the sublimation and recondensation of silicates. Our results indicate that, although the planetesimals at the two locations formed about contemporaneously, those at the snowline accreted a large fraction of their mass (approximately 60 percent) from materials delivered to the disk in the first few 10,000 yr, whereas this fraction is only 30 percent for the planetesimals formed at the silicate line. Thus, provided that the isotopic composition of the delivered material changed with time, these two planetesimal populations should have distinct isotopic compositions, consistent with observations.

preprint2021arXiv

Eigenstates of Quasi-Keplerian Self-Gravitating Particle Discs

Although quasi-Keplerian discs are among the most common astrophysical structures, computation of secular angular momentum transport within them routinely presents a considerable practical challenge. In this work, we investigate the secular small-inclination dynamics of a razor-thin particle disc as the continuum limit of a discrete Lagrange-Laplace secular perturbative theory and explore the analogy between the ensuing secular evolution -- including non-local couplings of self-gravitating discs -- and quantum mechanics. We find the 'quantum' Hamiltonian that describes the time evolution of the system and demonstrate the existence of a conserved inner product. The lowest-frequency normal modes are numerically approximated by performing a Wick rotation on the equations of motion. These modes are used to quantify the accuracy of a much simpler local-coupling model, revealing that it predicts the shape of the normal modes to a high degree of accuracy, especially in narrow annuli, even though it fails to predict their eigenfrequencies.

preprint2020arXiv

Dynamics of Planetary Systems Within Star Clusters: Aspects of the Solar System's Early Evolution

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°$, and are thus insufficient to explain the $ψ\approx6°$ 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, $η\,τ\lesssim2\times10^4\,$Myr/pc$^3$.

preprint2020arXiv

Energy Optimization in Binary Star Systems: Explanation for Equal Mass Members in Close Orbits

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

preprint2020arXiv

Energy Optimization in Extrasolar Planetary Systems: The Transition from Peas-in-a-Pod to Runaway Growth

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.

preprint2020arXiv

Formation of Giant Planet Satellites

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.

preprint2020arXiv

Interactions Among Non-Interacting Particles in Planet Formation Simulations

Over the course of the recent decades, $N$-body simulations have become a standard tool for quantifying the gravitational perturbations that ensue in planet-forming disks. Within the context of such simulations, massive non-central bodies are routinely classified into "big" and "small" particles, where big objects interact with all other objects self-consistently, while small bodies interact with big bodies but not with each other. Importantly, this grouping translates to an approximation scheme where the orbital evolution of small bodies is dictated entirely by the dynamics of the big bodies, yielding considerable computational advantages with little added cost in terms of astrophysical accuracy. Here we point out, however, that this scheme can also yield spurious dynamical behaviour, where even in absence of big bodies within a simulation, indirect coupling among small bodies can lead to excitation of the constituent "non-interacting" orbits. We demonstrate this self-stirring by carrying out a sequence of numerical experiments, and confirm that this effect is largely independent of the time-step or the employed integration algorithm. Furthermore, adopting the growth of angular momentum deficit as a proxy for dynamical excitation, we explore its dependence on time, the cumulative mass of the system, as well as the total number of particles present in the simulation. Finally, we examine the degree of such indirect excitation within the context of conventional terrestrial planet formation calculations, and conclude that although some level of caution may be warranted, this effect plays a negligible role in driving the simulated dynamical evolution.

preprint2020arXiv

On the Dynamical Origins of Retrograde Jupiter Trojans and their Connection to High-Inclination TNOs

Over the course of the last decade, observations of highly-inclined (orbital inclination i > 60°) Trans-Neptunian Objects (TNOs) have posed an important challenge to current models of solar system formation (Levison et al. 2008; Nesvorný 2015). These remarkable minor planets necessitate the presence of a distant reservoir of strongly-out-of-plane TNOs, which itself requires some dynamical production mechanism (Gladman et al. 2009; Gomes et al. 2015; Batygin and Brown 2016). A notable recent addition to the census of high-i minor bodies in the solar system is the retrograde asteroid 514107 Ka'epaoka'awela, which currently occupies a 1:-1 mean motion resonance with Jupiter at i = 163° (Wiegert et al. 2017). In this work, we delineate a direct connection between retrograde Jupiter Trojans and high-i Centaurs. First, we back-propagate a large sample of clones of Ka'epaoka'awela for 100 Ma numerically, and demonstrate that long-term stable clones tend to decrease their inclination steadily until it concentrates between 90° and 135°, while their eccentricity and semi-major axis increase, placing many of them firmly into the trans-Neptunian domain. Importantly, the clones show significant overlap with the synthetic high-i Centaurs generated in Planet 9 studies (Batygin et al. 2019), and hint at the existence of a relatively prominent, steady-state population of minor bodies occupying polar trans-Saturnian orbits. Second, through direct numerical forward-modeling, we delineate the dynamical pathway through which conventional members of the Kuiper Belt's scattered disk population can become retrograde Jovian Trojan resonators in presence of Planet 9.

preprint2020arXiv

Strong Near-Infrared Spectral Variability of the Young Cloudy L Dwarf Companion VHS J1256-1257 b

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 $μ$m over 8.5 hours and even higher variability at the 24.7% level at 1.27 $μ$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.

preprint2020arXiv

Tidal Inflation Reconciles Low-Density Sub-Saturns with Core Accretion

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.

preprint2019arXiv

K2-19b and c are in a 3:2 Commensurability but out of Resonance: A Challenge to Planet Assembly by Convergent Migration

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 $Δ\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 $Δ\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.

preprint2016arXiv

Direct Measure of Radiative and Dynamical Properties of an Exoplanet Atmosphere

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

preprint2016arXiv

Dynamical Constraints on the Core Mass of Hot Jupiter HAT-P-13b

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 $μ$m and 4.5 $μ$m bands, best match atmospheric model predictions with a dayside temperature inversion and relatively efficient day-night circulation.

preprint2016arXiv

Evidence for a Distant Giant Planet in the Solar System

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.

preprint2016arXiv

Friends of Hot Jupiters. IV. Stellar companions beyond 50 AU might facilitate giant planet formation, but most are unlikely to cause Kozai-Lidov migration

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σ$. 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σ$ 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.

preprint2016arXiv

Generation of Highly Inclined Trans-Neptunian Objects by Planet Nine

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.

preprint2016arXiv

In Situ Formation and Dynamical Evolution of Hot Jupiter Systems

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.

preprint2016arXiv

Observational constraints on the orbit and location of Planet Nine in the outer solar system

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.

preprint2016arXiv

Solar Obliquity Induced by Planet Nine

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.

preprint2016arXiv

Spin-orbit Misalignment as a Driver of the Kepler Dichotomy

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.

preprint2016arXiv

Statistics of Long Period Gas Giant Planets in Known Planetary Systems

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.

preprint2015arXiv

Capture of Planets Into Mean Motion Resonances and the Origins of Extrasolar Orbital Architectures

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.

preprint2015arXiv

Dynamical Evolution of Multi-Resonant Systems: the Case of GJ876

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.

preprint2015arXiv

Friends of Hot Jupiters II: No Correspondence Between Hot-Jupiter Spin-Orbit Misalignment and the Incidence of Directly Imaged Stellar Companions

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\,σ$ 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.

preprint2015arXiv

Friends of Hot Jupiters III: An Infrared Spectroscopic Search for Low-Mass Stellar Companions

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.

preprint2015arXiv

Jupiter's Decisive Role in the Inner Solar System's Early Evolution

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.

preprint2015arXiv

Magnetic Origins of the Stellar Mass-Obliquity Correlation in Planetary Systems

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.

preprint2015arXiv

Migration of two massive planets into (and out of) first order mean motion resonances

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.

preprint2015arXiv

Resonant Removal of Exomoons During Planetary Migration

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.

preprint2015arXiv

Spin-Spin Coupling in the Solar System

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.

preprint2014arXiv

Alignment of Protostars and Circumstellar Disks During the Embedded Phase

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

preprint2014arXiv

Chaotic Disintegration of the Inner Solar System

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.

preprint2014arXiv

Early Excitation of Spin-Orbit Misalignments in Close-in Planetary Systems

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.

preprint2014arXiv

Friends of Hot Jupiters I: A Radial Velocity Search for Massive, Long-Period Companions to Close-In Gas Giant Planets

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.

preprint2014arXiv

Non-Axisymmetric Flows on Hot Jupiters with Oblique Magnetic Fields

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.

preprint2014arXiv

On the Spin-axis Dynamics of a Moonless Earth

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.

preprint2014arXiv

Pre-LHB Evolution of the Earth's Obliquity

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.

preprint2014arXiv

Spin-Orbit angle distribution and the origin of (mis)aligned hot Jupiters

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 $Ψ$ to the projected one $β$. Then, a comparison with the observations is relevant. We analyse the geometry of the problem to link analytically the projected angle $β$ to the real spin-orbit angle $Ψ$. The distribution of $Ψ$ 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 $β$ knowing the PDF of $Ψ$ 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.

preprint2014arXiv

The Absolute Magnitude Distribution of Kuiper Belt Objects

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, $α_{\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, $α_{\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 $α_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 $α_\textrm{1}=1.0\pm0.2$, $α_\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.

preprint2014arXiv

WASP-12b and HAT-P-8b are Members of Triple Star Systems

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.

preprint2013arXiv

A Primordial Origin for Misalignments Between Stellar Spin Axes and Planetary Orbits

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.

preprint2013arXiv

Analytical Treatment of Planetary Resonances

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.

preprint2013arXiv

Dynamical Measurements of the Interior Structure of Exoplanets

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.

preprint2013arXiv

Magnetic and Gravitational Disk-Star Interactions: An Interdependence of PMS Stellar Rotation Rates and Spin-Orbit Misalignments

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.

preprint2013arXiv

Magnetically Controlled Circulation on Hot Extrasolar Planets

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.

preprint2013arXiv

Mass-Radius Relationships for Very Low Mass Gaseous Planets

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.

preprint2012arXiv

Dissipative Divergence of Resonant Orbits

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.

preprint2012arXiv

Evolution of Ohmically Heated Hot Jupiters

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.

preprint2012arXiv

NOTE: Explaining why the Uranian satellites have equatorial prograde orbits despite the large planetary obliquity

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.

preprint2011arXiv

Formation and Evolution of Planetary Systems in Presence of Highly Inclined Stellar Perturbers

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.

preprint2011arXiv

Instability-Driven Dynamical Evolution Model of a Primordially 5 Planet Outer Solar System

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.

preprint2011arXiv

Onset of Secular Chaos in Planetary Systems: Period Doubling & Strange Attractors

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.

preprint2011arXiv

Resolving the Sin(I) degeneracy in Low-Mass Multi-Planet Systems

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.

preprint2011arXiv

Retention of a Primordial Cold Classical Kuiper Belt in an Instability-Driven Model of Solar System Formation

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.

preprint2010arXiv

Early Dynamical Evolution of the Solar System: Pinning Down the Initial Condition of the Nice Model

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.

preprint2010arXiv

Inflating Hot Jupiters With Ohmic Dissipation

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.

preprint2009arXiv

A Quasi-Stationary Solution to Gliese 436b's Eccentricity

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.

preprint2009arXiv

A Super-Earth and two Neptunes Orbiting the Nearby Sun-like star 61 Virginis

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.

preprint2009arXiv

Determination of the Interior Structure of Transiting Planets in Multiple-Planet 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.

Konstantin Batygin

What is connected

Connect this record

See the researcher in context

Building this map preview

72 published item(s)

A young progenitor for the most common planetary systems in the Galaxy

Callisto's Nonresonant Orbit as an Outcome of Circum-Jovian Disk Substructure

Secular Excitation of Polar Neptune Orbits in Pure Disk-Planet Systems

Formation of Rocky Super-Earths From A Narrow Ring of Planetesimals

A Criterion for the Stability of Planets in Chains of Resonances

A Multi-Planet System's Sole Super-Puff: Exploring Allowable Physical Parameters for the Cold Super-Puff HIP 41378 f

Analytic Approach to the Late Stages of Giant Planet Formation

Architectures of Compact Super-Earth Systems Shaped by Instabilities

Generating Stellar Obliquity in Systems with Broken Protoplanetary Disks

Shallowness of circulation in hot Jupiters -- Advancing the Ohmic dissipation model

TESS-Keck Survey IX: Masses of Three Sub-Neptunes Orbiting HD 191939 and the Discovery of a Warm Jovian Plus a Distant Sub-Stellar Companion

TOI-1136 is a Young, Coplanar, Aligned Planetary System in a Pristine Resonant Chain

A search for Planet Nine using the Zwicky Transient Facility public archive

Contemporary formation of early solar system planetesimals at two distinct radial locations

Eigenstates of Quasi-Keplerian Self-Gravitating Particle Discs

Dynamics of Planetary Systems Within Star Clusters: Aspects of the Solar System's Early Evolution

Energy Optimization in Binary Star Systems: Explanation for Equal Mass Members in Close Orbits

Energy Optimization in Extrasolar Planetary Systems: The Transition from Peas-in-a-Pod to Runaway Growth

Formation of Giant Planet Satellites

Interactions Among Non-Interacting Particles in Planet Formation Simulations

On the Dynamical Origins of Retrograde Jupiter Trojans and their Connection to High-Inclination TNOs

Strong Near-Infrared Spectral Variability of the Young Cloudy L Dwarf Companion VHS J1256-1257 b

Tidal Inflation Reconciles Low-Density Sub-Saturns with Core Accretion

K2-19b and c are in a 3:2 Commensurability but out of Resonance: A Challenge to Planet Assembly by Convergent Migration

Direct Measure of Radiative and Dynamical Properties of an Exoplanet Atmosphere

Dynamical Constraints on the Core Mass of Hot Jupiter HAT-P-13b

Evidence for a Distant Giant Planet in the Solar System

Friends of Hot Jupiters. IV. Stellar companions beyond 50 AU might facilitate giant planet formation, but most are unlikely to cause Kozai-Lidov migration

Generation of Highly Inclined Trans-Neptunian Objects by Planet Nine

In Situ Formation and Dynamical Evolution of Hot Jupiter Systems

Observational constraints on the orbit and location of Planet Nine in the outer solar system

Solar Obliquity Induced by Planet Nine

Spin-orbit Misalignment as a Driver of the Kepler Dichotomy

Statistics of Long Period Gas Giant Planets in Known Planetary Systems

Capture of Planets Into Mean Motion Resonances and the Origins of Extrasolar Orbital Architectures

Dynamical Evolution of Multi-Resonant Systems: the Case of GJ876

Friends of Hot Jupiters II: No Correspondence Between Hot-Jupiter Spin-Orbit Misalignment and the Incidence of Directly Imaged Stellar Companions

Friends of Hot Jupiters III: An Infrared Spectroscopic Search for Low-Mass Stellar Companions

Jupiter's Decisive Role in the Inner Solar System's Early Evolution

Magnetic Origins of the Stellar Mass-Obliquity Correlation in Planetary Systems

Migration of two massive planets into (and out of) first order mean motion resonances

Resonant Removal of Exomoons During Planetary Migration

Spin-Spin Coupling in the Solar System

Alignment of Protostars and Circumstellar Disks During the Embedded Phase

Chaotic Disintegration of the Inner Solar System

Early Excitation of Spin-Orbit Misalignments in Close-in Planetary Systems

Friends of Hot Jupiters I: A Radial Velocity Search for Massive, Long-Period Companions to Close-In Gas Giant Planets

Non-Axisymmetric Flows on Hot Jupiters with Oblique Magnetic Fields

On the Spin-axis Dynamics of a Moonless Earth

Pre-LHB Evolution of the Earth's Obliquity

Spin-Orbit angle distribution and the origin of (mis)aligned hot Jupiters

The Absolute Magnitude Distribution of Kuiper Belt Objects

WASP-12b and HAT-P-8b are Members of Triple Star Systems

A Primordial Origin for Misalignments Between Stellar Spin Axes and Planetary Orbits

Analytical Treatment of Planetary Resonances

Dynamical Measurements of the Interior Structure of Exoplanets

Magnetic and Gravitational Disk-Star Interactions: An Interdependence of PMS Stellar Rotation Rates and Spin-Orbit Misalignments

Magnetically Controlled Circulation on Hot Extrasolar Planets

Mass-Radius Relationships for Very Low Mass Gaseous Planets

Dissipative Divergence of Resonant Orbits

Evolution of Ohmically Heated Hot Jupiters

NOTE: Explaining why the Uranian satellites have equatorial prograde orbits despite the large planetary obliquity

Formation and Evolution of Planetary Systems in Presence of Highly Inclined Stellar Perturbers

Instability-Driven Dynamical Evolution Model of a Primordially 5 Planet Outer Solar System

Onset of Secular Chaos in Planetary Systems: Period Doubling & Strange Attractors

Resolving the Sin(I) degeneracy in Low-Mass Multi-Planet Systems

Retention of a Primordial Cold Classical Kuiper Belt in an Instability-Driven Model of Solar System Formation

Early Dynamical Evolution of the Solar System: Pinning Down the Initial Condition of the Nice Model

Inflating Hot Jupiters With Ohmic Dissipation

A Quasi-Stationary Solution to Gliese 436b's Eccentricity

A Super-Earth and two Neptunes Orbiting the Nearby Sun-like star 61 Virginis

Determination of the Interior Structure of Transiting Planets in Multiple-Planet Systems