Showing 1–50 of 131 results for author: Batygin, K

Forming the Trappist-1 system in two steps during the recession of the disc inner edge

Authors: Gabriele Pichierri, Alessandro Morbidelli, Konstantin Batygin, Ramon Brasser

Abstract: Trappist-1 hosts 7 planets where period ratios of neighbouring pairs are close to the 8:5, 5:3, 3:2, 3:2, 4:3, and 3:2 ratios in increasing distance from the star. The Laplace angles associated with neighbouring triplets are observed to be librating, proving the resonant nature of the system. This compact, resonant configuration is a manifest sign of disc-driven migration; however, the preferred o… ▽ More Trappist-1 hosts 7 planets where period ratios of neighbouring pairs are close to the 8:5, 5:3, 3:2, 3:2, 4:3, and 3:2 ratios in increasing distance from the star. The Laplace angles associated with neighbouring triplets are observed to be librating, proving the resonant nature of the system. This compact, resonant configuration is a manifest sign of disc-driven migration; however, the preferred outcome of such evolution is the establishment of first-order resonances, not the high-order resonances observed in the inner system. Here, we explain the observed orbital configuration in a model that is largely independent on the specific disc migration and orbital circularisation efficiencies. Together with migration, the two key elements of our model are: i) the inner border of the protoplanetary disc receded with time; and ii) the system was initially separated in two sub-systems. Specifically, the inner b, c d and e planets were initially placed in a 3:2 resonance chain and then evolved to the 8:5 -- 5:3 commensurability between planets b, c and d under the effect of the recession of the inner edge of the disc, whereas the outer planets migrated to the inner edge at a later time, establishing the remaining resonances. Our results pivot on the dynamical role of the presently unobservable recession of the inner edge of protoplanetary discs. They also reveal the role of recurring phases of convergent migration followed by resonant repulsion with associated orbital circularisation when resonant chains interact with migration barriers. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: To appear in Nature Astronomy

The Prevalence of Resonance Among Young, Close-in Planets

Authors: Fei Dai, Max Goldberg, Konstantin Batygin, Jennifer van Saders, Eugene Chiang, Nick Choksi, Rixin Li, Erik A. Petigura, Gregory J. Gilbert, Sarah C. Millholland, Yuan-Zhe Dai, Luke Bouma, Lauren M. Weiss, Joshua N. Winn

Abstract: Multiple planets undergoing disk migration may be captured into a chain of mean-motion resonances with the innermost planet parked near the disk's inner edge. Subsequent dynamical evolution may disrupt these resonances, leading to the non-resonant configurations typically observed among {\it Kepler} planets that are Gyrs old. In this scenario, resonant configurations are expected to be more common… ▽ More Multiple planets undergoing disk migration may be captured into a chain of mean-motion resonances with the innermost planet parked near the disk's inner edge. Subsequent dynamical evolution may disrupt these resonances, leading to the non-resonant configurations typically observed among {\it Kepler} planets that are Gyrs old. In this scenario, resonant configurations are expected to be more common in younger systems. This prediction can now be tested, thanks to recent discoveries of young planets, particularly those in stellar clusters, by NASA's {\it TESS} mission. We divided the known planetary systems into three age groups: young ($<$100-Myr-old), adolescent (0.1-1-Gyr-old), and mature ($>1$-Gyr-old). The fraction of neighboring planet pairs having period ratios within a few percent of a first-order commensurability (e.g.~4:3, 3:2, or 2:1) is 70$\pm$15\% for young pairs, 24$\pm$8\% for adolescent pairs, and 15$\pm$2\% for mature pairs. The fraction of systems with at least one nearly commensurable pair (either first or second-order) is 86$\pm13$\% among young systems, 38$\pm12$\% for adolescent systems, and 23$\pm3$\% for mature systems. First-order commensurabilities prevail across all age groups, with an admixture of second-order commensurabilities. Commensurabilities are more common in systems with high planet multiplicity and low mutual inclinations. Observed period ratios often deviate from perfect commensurability by $\sim$1\% even among young planets, too large to be explained by resonant repulsion with equilibrium eccentricity tides. We also find that super-Earths in the radius gap ($1.5-1.9R_\oplus$) are less likely to be near-resonant (11.9$\pm2.0\%$) compared to Earth-sized planets ($R_p<1R_\oplus$; 25.3$\pm4.4\%$) or mini-Neptunes ($1.9R_\oplus \leq R_p<2.5R_\oplus$; 14.4$\pm1.8\%$). △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 17 pages, 9 figures, submitted to AAS

Self-Gravitational Dynamics Within the Inner Oort Cloud

Authors: Konstantin Batygin, David Nesvorný

Abstract: The formation of the Inner Oort Cloud (IOC) - a vast halo of icy bodies residing far beyond Neptune's orbit - is an expected outcome of the solar system's primordial evolution within a stellar cluster. Recent models have shown that the process of early planetesimal capture within the trans-Neptunian region may have been sufficiently high for the cumulative mass of the Cloud to approach several Ear… ▽ More The formation of the Inner Oort Cloud (IOC) - a vast halo of icy bodies residing far beyond Neptune's orbit - is an expected outcome of the solar system's primordial evolution within a stellar cluster. Recent models have shown that the process of early planetesimal capture within the trans-Neptunian region may have been sufficiently high for the cumulative mass of the Cloud to approach several Earth masses. In light of this, here we examine the dynamical evolution of the IOC, driven by its own self-gravity. We show that the collective gravitational potential of the IOC is adequately approximated by the Miyamoto-Nagai model, and use a semi-analytic framework to demonstrate that the resulting secular oscillations are akin to the von Zeipel-Lidov-Kozai resonance. We verify our results with direct $N$-body calculations, and examine the effects of IOC self-gravity on the long-term behavior of the solar system's minor bodies using a detailed simulation. Cumulatively, we find that while the modulation of perihelion distances and inclinations can occur within an observationally relevant range, the associated timescales vastly surpass the age of the sun, indicating that the influence of IOC self-gravity on the architecture of the solar system is negligible. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 14 pages, 6 figures, accepted for publication in Celestial Mechanics and Dynamical Astronomy

Generation of Low-Inclination, Neptune-Crossing TNOs by Planet Nine

Authors: Konstantin Batygin, Alessandro Morbidelli, Michael E. Brown, David Nesvorny

Abstract: The solar system's distant reaches exhibit a wealth of anomalous dynamical structure, hinting at the presence of a yet-undetected, massive trans-Neptunian body - Planet 9. Previous analyses have shown how orbital evolution induced by this object can explain the origins of a broad assortment of exotic orbits, ranging from those characterized by high perihelia to those with extreme inclinations. In… ▽ More The solar system's distant reaches exhibit a wealth of anomalous dynamical structure, hinting at the presence of a yet-undetected, massive trans-Neptunian body - Planet 9. Previous analyses have shown how orbital evolution induced by this object can explain the origins of a broad assortment of exotic orbits, ranging from those characterized by high perihelia to those with extreme inclinations. In this work, we shift the focus toward a more conventional class of TNOs, and consider the observed census of long-period, nearly planar, Neptune-crossing objects as a hitherto-unexplored probe of the Planet 9 hypothesis. To this end, we carry out comprehensive $N-$body simulations that self-consistently model gravitational perturbations from all giant planets, the Galactic tide, as well as passing stars, stemming from initial conditions that account for the primordial giant planet migration and sun's early evolution within a star cluster. Accounting for observational biases, our results reveal that the orbital architecture of this group of objects aligns closely with the predictions of the P9-inclusive model. In stark contrast, the P9-free scenario is statistically rejected at a $\sim5\,σしぐま$ confidence-level. Accordingly, this work introduces a new line of evidence supporting the existence of Planet 9 and further delineates a series of observational predictions poised for near-term resolution. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 11 pages, 5 figures, accepted for publication in The Astrophysical Journal Letters

Breaking Giant Chains: Early-Stage Instabilities in Long-Period Giant Planet Systems

Authors: Vighnesh Nagpal, Max Goldberg, Konstantin Batygin

Abstract: Orbital evolution is a critical process that sculpts planetary systems, particularly during their early stages where planet-disk interactions are expected to lead to the formation of resonant chains. Despite the theoretically expected prominence of such configurations, they are scarcely observed among long-period giant exoplanets. This disparity suggests an evolutionary sequence wherein giant plan… ▽ More Orbital evolution is a critical process that sculpts planetary systems, particularly during their early stages where planet-disk interactions are expected to lead to the formation of resonant chains. Despite the theoretically expected prominence of such configurations, they are scarcely observed among long-period giant exoplanets. This disparity suggests an evolutionary sequence wherein giant planet systems originate in compact multi-resonant configurations, but subsequently become unstable, eventually relaxing to wider orbits--a phenomenon mirrored in our own solar system's early history. In this work, we present a suite of N-body simulations that model the instability-driven evolution of giant planet systems, originating from resonant initial conditions, through phases of disk-dispersal and beyond. By comparing the period ratio and normalized angular momentum deficit distributions of our synthetic aggregate of systems with the observational census of long-period Jovian planets, we derive constraints on the expected rate of orbital migration, efficiency of gas-driven eccentricity damping, and typical initial multiplicity. Our findings reveal a distinct inclination towards densely-packed initial conditions, weak damping, and high giant planet multiplicities. Furthermore, our models indicate that resonant chain origins do not facilitate the formation of Hot Jupiters via the coplanar high-eccentricity pathway at rates high enough to explain their observed prevalence. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 14 pages, 8 figures. Accepted to ApJ

Chaotic tides as a solution to the Hyperion problem

Authors: Max Goldberg, Konstantin Batygin

Abstract: The dynamics of the outer regular satellites of Saturn are driven primarily by the outward migration of Titan, but several independent constraints on Titan's migration are difficult to reconcile with the current resonant orbit of the small satellite Hyperion. We argue that Hyperion's rapid irregular tumbling greatly increases tidal dissipation with a steep dependence on orbital eccentricity. Reson… ▽ More The dynamics of the outer regular satellites of Saturn are driven primarily by the outward migration of Titan, but several independent constraints on Titan's migration are difficult to reconcile with the current resonant orbit of the small satellite Hyperion. We argue that Hyperion's rapid irregular tumbling greatly increases tidal dissipation with a steep dependence on orbital eccentricity. Resonant excitation from a migrating Titan is then balanced by damping in a feedback mechanism that maintains Hyperion's eccentricity without fine-tuning. The inferred tidal parameters of Hyperion are most consistent with rapid Titan migration enabled by a resonance lock with an internal mode of Saturn, but a scenario with only equilibrium dissipation in Saturn is also possible. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 11 pages, 4 figures, accepted to Icarus

A Pan-STARRS1 Search for Planet Nine

Authors: Michael E. Brown, Matthew J. Holman, Konstantin Batygin

Abstract: We present a search for Planet Nine using the second data release of the Pan-STARRS1survey. We rule out the existence of a Planet Nine with the characteristics of that predicted in Brown & Batygin (2021) to a 50% completion depth of $V=21.5$. This survey, along with previous analyses of the Zwicky Transient Facility (ZTF) and Dark EnergySurvey (DES) data, rules out 78% of the Brown \& Batygin para… ▽ More We present a search for Planet Nine using the second data release of the Pan-STARRS1survey. We rule out the existence of a Planet Nine with the characteristics of that predicted in Brown & Batygin (2021) to a 50% completion depth of $V=21.5$. This survey, along with previous analyses of the Zwicky Transient Facility (ZTF) and Dark EnergySurvey (DES) data, rules out 78% of the Brown \& Batygin parameter space. Much of the remaining parameter space is at $V>21$ in regions near and in the area where the northern galactic plane crosses the ecliptic. △ Less

Submitted 31 January, 2024; originally announced January 2024.

Comments: AJ, in press

Radial Distribution of Distant Trans-Neptunian Objects Points to Sun's Formation in a Stellar Cluster

Authors: David Nesvorny, Pedro Bernardinelli, David Vokrouhlicky, Konstantin Batygin

Abstract: The Scattered Disk Objects (SDOs) are a population of trans-Neptunian bodies with semimajor axes $50< a \lesssim 1000$ au and perihelion distances $q \gtrsim 30$ au. The detached SDOs with orbits beyond the reach of Neptune (roughly $q>35$~au) are of special interest here as an important constraint on the early evolution of the outer Solar System. The semimajor axis profile of detached SDOs at 50-… ▽ More The Scattered Disk Objects (SDOs) are a population of trans-Neptunian bodies with semimajor axes $50< a \lesssim 1000$ au and perihelion distances $q \gtrsim 30$ au. The detached SDOs with orbits beyond the reach of Neptune (roughly $q>35$~au) are of special interest here as an important constraint on the early evolution of the outer Solar System. The semimajor axis profile of detached SDOs at 50--500~au, as characterized from the Dark Energy Survey (DES), is radially extended, but previous dynamical models of Neptune's early migration produce a relatively compact profile. This problem is most likely related to Sun's birth environment in a stellar cluster. We perform new dynamical simulations that account for cluster effects and show that the orbital distribution of SDOs can be explained if a particularly close stellar encounter occurred early on (e.g., M dwarf with the mass $\simeq 0.2$ $M_\odot$ approaching the Sun at $\simeq 200$ au). For such an encounter to happen with a reasonably high probability the Sun must have formed in a stellar cluster with $ηいーたT \gtrsim 10^4$ Myr pc$^{-3}$, where $ηいーた$ is the stellar number density and $T$ is the Sun's residence time in the cluster. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Suppression of the inclination instability in the trans-Neptunian Solar system

Authors: Arnav Das, Konstantin Batygin

Abstract: The trans-Neptunian scattered disk exhibits unexpected dynamical structure, ranging from an extended dispersion of perihelion distance to a clustered distribution in orbital angles. Self-gravitational modulation of the scattered disk has been suggested in the literature as an alternative mechanism to Planet 9 for sculpting the orbital architecture of the trans-Neptunian region. The numerics of thi… ▽ More The trans-Neptunian scattered disk exhibits unexpected dynamical structure, ranging from an extended dispersion of perihelion distance to a clustered distribution in orbital angles. Self-gravitational modulation of the scattered disk has been suggested in the literature as an alternative mechanism to Planet 9 for sculpting the orbital architecture of the trans-Neptunian region. The numerics of this hypothesis have hitherto been limited to $N < O(10^3)$ super-particle simulations that omit direct gravitational perturbations from the giant planets and instead model them as an orbit-averaged (quadrupolar) potential, through an enhanced $J_2$ moment of the central body. For sufficiently massive disks, such simulations reveal the onset of collective dynamical behaviour $\unicode{x2014}$ termed the $\unicode{x2018}$inclination instability$\unicode{x2019}$ $\unicode{x2014}$ wherein orbital circularisation occurs at the expense of coherent excitation of the inclination. Here, we report $N = O(10^4)$ GPU-accelerated simulations of a self-gravitating scattered disk (across a range of disk masses spanning 5 to 40 Earth masses) that self-consistently account for intra-particle interactions as well as Neptune's perturbations. Our numerical experiments show that even under the most favourable conditions, the inclination instability never ensues. Instead, due to scattering, the disk depletes. While our calculations show that a transient lopsided structure can emerge within the first few hundreds of Myr, the terminal outcomes of these calculations systematically reveal a scattered disk that is free of any orbital clustering. We conclude thus that the inclination instability mechanism is an inadequate explanation of the observed architecture of the solar system. △ Less

Submitted 1 July, 2023; originally announced July 2023.

The Origin of Universality in the Inner Edges of Planetary Systems

Authors: Konstantin Batygin, Fred C. Adams, Juliette Becker

Abstract: The characteristic orbital period of the inner-most objects within the galactic census of planetary and satellite systems appears to be nearly universal, with $P$ on the order of a few days. This paper presents a theoretical framework that provides a simple explanation for this phenomenon. By considering the interplay between disk accretion, magnetic field generation by convective dynamos, and Kel… ▽ More The characteristic orbital period of the inner-most objects within the galactic census of planetary and satellite systems appears to be nearly universal, with $P$ on the order of a few days. This paper presents a theoretical framework that provides a simple explanation for this phenomenon. By considering the interplay between disk accretion, magnetic field generation by convective dynamos, and Kelvin-Helmholtz contraction, we derive an expression for the magnetospheric truncation radius in astrophysical disks, and find that the corresponding orbital frequency is independent of the mass of the host body. Our analysis demonstrates that this characteristic frequency corresponds to a period of $P\sim3$ days, although intrinsic variations in system parameters are expected to introduce a factor of $\sim2-3$ spread in this result. Standard theory of orbital migration further suggests that planets should stabilize at an orbital period that exceeds disk truncation by a small margin. Cumulatively, our findings predict that the periods of close-in bodies should span $P\sim2-12$ days - a range that is consistent with observations. △ Less

Submitted 14 June, 2023; originally announced June 2023.

Comments: 9 pages, 4 figures, accepted to ApJL

In-situ enrichment in heavy elements of hot Jupiters

Authors: A. Morbidelli, K. Batygin, E. Lega

Abstract: Context: Radius and mass measurements of short-period giant planets reveal that many of these planets contain a large amount of heavy elements, in sharp contrast with the expectations of the conventional core-accretion model for the origin of giant planets. Aims: The proposed explanations for the heavy-element enrichment of giant planets fall short of explaining the most enriched planets. We loo… ▽ More Context: Radius and mass measurements of short-period giant planets reveal that many of these planets contain a large amount of heavy elements, in sharp contrast with the expectations of the conventional core-accretion model for the origin of giant planets. Aims: The proposed explanations for the heavy-element enrichment of giant planets fall short of explaining the most enriched planets. We look for additional processes that can explain the full envelope of inferred enrichments. Methods: We revisit the dynamics of pebbles and dust in the vicinity of giant planets using analytic estimates. Although our results are derived in the framework of a viscous alpha-disk we also discuss the case of disks driven by angular momentum removal in magnetized winds. Results: When giant planets are far from the star, dust and pebbles are confined in a pressure bump at the outer edge of the planet-induced gap. Instead, when the planets reach the inner part of the disk (r << 2 au), dust penetrates the gap together with the gas. The dust/gas ratio can be enhanced by more than an order of magnitude if radial drift of dust is not impeded farther out by other barriers. Thus, hot planets undergoing runaway gas accretion can swallow a large amount of dust. Conclusions: Whereas the gas accreted by giant planets in the outer disk is very dust-poor, that accreted by hot planets can be extremely dust-rich. Thus, provided that a large fraction of the atmosphere of hot-Jupiters is accreted in situ, a large amount of dust can be accreted as well. We draw a distinction between this process and pebble accretion, which is ineffective at small stellocentric radii, even for super-Earths. Giant planets farther out in the disk are extremely effective barriers against the flow of pebbles and dust across their gap. △ Less

Submitted 2 June, 2023; originally announced June 2023.

Comments: In press in Astronomy & Astrophysics

Journal ref: A&A 675, A75 (2023)

Science opportunities with solar sailing smallsats

Authors: Slava G. Turyshev, Darren Garber, Louis D. Friedman, Andreas M. Hein, Nathan Barnes, Konstantin Batygin, G. David Brin, Michael E. Brown, Leroy Cronin, Artur Davoyan, Amber Dubill, T. Marshall Eubanks, Sarah Gibson, Donald M. Hassler, Noam R. Izenberg, Pierre Kervella, Philip D. Mauskopf, Neil Murphy, Andrew Nutter, Carolyn Porco, Dario Riccobono, James Schalkwyk, Kevin B. Stevenson, Mark V. Sykes, Mahmooda Sultana , et al. (3 additional authors not shown)

Abstract: Recently, we witnessed how the synergy of small satellite technology and solar sailing propulsion enables new missions. Together, small satellites with lightweight instruments and solar sails offer affordable access to deep regions of the solar system, also making it possible to realize hard-to-reach trajectories that are not constrained to the ecliptic plane. Combining these two technologies can… ▽ More Recently, we witnessed how the synergy of small satellite technology and solar sailing propulsion enables new missions. Together, small satellites with lightweight instruments and solar sails offer affordable access to deep regions of the solar system, also making it possible to realize hard-to-reach trajectories that are not constrained to the ecliptic plane. Combining these two technologies can drastically reduce travel times within the solar system, while delivering robust science. With solar sailing propulsion capable of reaching the velocities of ~5-10 AUえーゆー/yr, missions using a rideshare launch may reach the Jovian system in two years, Saturn in three. The same technologies could allow reaching solar polar orbits in less than two years. Fast, cost-effective, and maneuverable sailcraft that may travel outside the ecliptic plane open new opportunities for affordable solar system exploration, with great promise for heliophysics, planetary science, and astrophysics. Such missions could be modularized to reach different destinations with different sets of instruments. Benefiting from this progress, we present the "Sundiver" concept, offering novel possibilities for the science community. We discuss some of the key technologies, the current design of the Sundiver sailcraft vehicle and innovative instruments, along with unique science opportunities that these technologies enable, especially as this exploration paradigm evolves. We formulate policy recommendations to allow national space agencies, industry, and other stakeholders to establish a strong scientific, programmatic, and commercial focus, enrich and deepen the space enterprise and broaden its advocacy base by including the Sundiver paradigm as a part of broader space exploration efforts. △ Less

Submitted 19 July, 2023; v1 submitted 27 March, 2023; originally announced March 2023.

Comments: 35 pages, 12 figures, 3 tables

Journal ref: Planetary and Space Science 235 (2023) 105744

Dissipative Capture of Planets Into First-Order Mean-Motion Resonances

Authors: Konstantin Batygin, Antoine C. Petit

Abstract: The emergence of orbital resonances among planets is a natural consequence of the early dynamical evolution of planetary systems. While it is well-established that convergent migration is necessary for mean-motion commensurabilities to emerge, recent numerical experiments have shown that the existing adiabatic theory of resonant capture provides an incomplete description of the relevant physics, l… ▽ More The emergence of orbital resonances among planets is a natural consequence of the early dynamical evolution of planetary systems. While it is well-established that convergent migration is necessary for mean-motion commensurabilities to emerge, recent numerical experiments have shown that the existing adiabatic theory of resonant capture provides an incomplete description of the relevant physics, leading to an erroneous mass scaling in the regime of strong dissipation. In this work, we develop a new model for resonance capture that self-consistently accounts for migration and circularization of planetary orbits, and derive an analytic criterion based upon stability analysis that describes the conditions necessary for the formation of mean-motion resonances. We subsequently test our results against numerical simulations and find satisfactory agreement. Our results elucidate the critical role played by adiabaticity and resonant stability in shaping the orbital architectures of planetary systems during the nebular epoch, and provide a valuable tool for understanding their primordial dynamical evolution. △ Less

Submitted 5 March, 2023; originally announced March 2023.

Comments: 12 pages, 4 figures, accepted for publication in ApJL

A 2:1 Mean-Motion Resonance Super-Jovian pair revealed by TESS, FEROS, and HARPS

Authors: Vladimir Bozhilov, Desislava Antonova, Melissa J. Hobson, Rafael Brahm, Andres Jordan, Thomas Henning, Jan Eberhardt, Felipe I. Rojas, Konstantin Batygin, Pascal Torres-Miranda, Keivan G. Stassun, Sarah C. Millholland, Denitza Stoeva, Milen Minev, Nestor Espinoza, George R. Ricker, David W. Latham, Diana Dragomir, Michelle Kunimoto, Jon M. Jenkins, Eric B. Ting, Sara Seager, Joshua N. Winn, Jesus Noel Villasenor, Luke G. Bouma , et al. (2 additional authors not shown)

Abstract: We report the discovery of a super-Jovian 2:1 mean-motion resonance (MMR) pair around the G-type star TIC 279401253, whose dynamical architecture is a prospective benchmark for planet formation and orbital evolution analysis. The system was discovered thanks to a single transit event recorded by the Transiting Exoplanet Survey Satellite (TESS) mission, which pointed to a Jupiter-sized companion wi… ▽ More We report the discovery of a super-Jovian 2:1 mean-motion resonance (MMR) pair around the G-type star TIC 279401253, whose dynamical architecture is a prospective benchmark for planet formation and orbital evolution analysis. The system was discovered thanks to a single transit event recorded by the Transiting Exoplanet Survey Satellite (TESS) mission, which pointed to a Jupiter-sized companion with poorly constrained orbital parameters. We began ground-based precise radial velocity (RV) monitoring with HARPS and FEROS within the Warm gIaNts with tEss (WINE) survey to constrain the transiting body's period, mass, and eccentricity. The RV measurements revealed not one but two massive planets with periods of 76.80$_{-0.06}^{+0.06}$ days and 155.3$_{-0.7}^{+0.7}$ days, respectively. A combined analysis of transit and RV data yields an inner transiting planet with a mass of 6.14$_{-0.42}^{+0.39}$ M$_{\rm Jup}$ and a radius of 1.00$_{-0.04}^{+0.04}$ R$_{\rm Jup}$, and an outer planet with a minimum mass of 8.02$_{-0.18}^{+0.18}$ M$_{\rm Jup}$, indicating a massive giant pair. A detailed dynamical analysis of the system reveals that the planets are locked in a strong first-order, eccentricity-type 2:1 MMR, which makes TIC 279401253 one of the rare examples of truly resonant architectures supporting disk-induced planet migration. The bright host star, $V \approx$ 11.9 mag, the relatively short orbital period ($P_{\rm b}$ = 76.80$_{-0.06}^{+0.06}$ d) and pronounced eccentricity (e =0.448$_{-0.029}^{+0.028}$) make the transiting planet a valuable target for atmospheric investigation with the James Webb Space Telescope (JWST) and ground-based extremely-large telescopes. △ Less

Submitted 11 May, 2023; v1 submitted 21 February, 2023; originally announced February 2023.

Comments: Published in the ApJL

Journal ref: ApJL 946 L36 (2023)

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

Authors: Konstantin Batygin, Alessandro Morbidelli

Abstract: 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 belongin… ▽ More 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. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: 14 pages, 4 figures, accepted for publication in Nature Astronomy

Detecting exomoons from radial velocity measurements of self-luminous planets: application to observations of HR 7672 B and future prospects

Authors: Jean-Baptiste Ruffio, Katelyn Horstman, Dimitri Mawet, Lee J. Rosenthal, Konstantin Batygin, Jason J. Wang, Maxwell Millar-Blanchaer, Ji Wang, Benjamin J. Fulton, Quinn M. Konopacky, Shubh Agrawal, Lea A. Hirsch, Andrew W. Howard, Sarah Blunt, Eric Nielsen, Ashley Baker, Randall Bartos, Charlotte Z. Bond, Benjamin Calvin, Sylvain Cetre, Jacques-Robert Delorme, Greg Doppmann, Daniel Echeverri, Luke Finnerty, Michael P. Fitzgerald , et al. (14 additional authors not shown)

Abstract: The detection of satellites around extrasolar planets, so called exomoons, remains a largely unexplored territory. In this work, we study the potential of detecting these elusive objects from radial velocity monitoring of self-luminous directly imaged planets. This technique is now possible thanks to the development of dedicated instruments combining the power of high-resolution spectroscopy and h… ▽ More The detection of satellites around extrasolar planets, so called exomoons, remains a largely unexplored territory. In this work, we study the potential of detecting these elusive objects from radial velocity monitoring of self-luminous directly imaged planets. This technique is now possible thanks to the development of dedicated instruments combining the power of high-resolution spectroscopy and high-contrast imaging. First, we demonstrate a sensitivity to satellites with a mass ratio of 1-4% at separations similar to the Galilean moons from observations of a brown-dwarf companion (HR 7672 B; Kmag=13; 0.7" separation) with the Keck Planet Imager and Characterizer (KPIC; R~35,000 in K band) at the W. M. Keck Observatory. Current instrumentation is therefore already sensitive to large unresolved satellites that could be forming from gravitational instability akin to binary star formation. Using end-to-end simulations, we then estimate that future instruments such as MODHIS, planned for the Thirty Meter Telescope, should be sensitive to satellites with mass ratios of ~1e-4. Such small moons would likely form in a circumplanetary disk similar to the Jovian satellites in the solar system. Looking for the Rossiter-McLaughlin effect could also be an interesting pathway to detecting the smallest moons on short orbital periods. Future exomoon discoveries will allow precise mass measurements of the substellar companions that they orbit and provide key insight into the formation of exoplanets. They would also help constrain the population of habitable Earth-sized moons orbiting gas giants in the habitable zone of their stars. △ Less

Submitted 6 February, 2023; v1 submitted 10 January, 2023; originally announced January 2023.

Comments: Accepted to AJ (Jan 10, 2023)

Dynamics and Origins of the Near-Resonant Kepler Planets

Authors: Max Goldberg, Konstantin Batygin

Abstract: Short-period super-Earths and mini-Neptunes encircle more than $\sim50\%$ of Sun-like stars and are relatively amenable to direct observational characterization. Despite this, environments in which these planets accrete are difficult to probe directly. Nevertheless, pairs of planets that are close to orbital resonances provide a unique window into the inner regions of protoplanetary disks, as they… ▽ More Short-period super-Earths and mini-Neptunes encircle more than $\sim50\%$ of Sun-like stars and are relatively amenable to direct observational characterization. Despite this, environments in which these planets accrete are difficult to probe directly. Nevertheless, pairs of planets that are close to orbital resonances provide a unique window into the inner regions of protoplanetary disks, as they preserve the conditions of their formation, as well as the early evolution of their orbital architectures. In this work, we present a novel approach toward quantifying transit timing variations within multi-planetary systems and examine the near-resonant dynamics of over 100 planet pairs detected by \textit{Kepler}. Using an integrable model for first-order resonances, we find a clear transition from libration to circulation of the resonant angle at a period ratio of $\approx 0.6\%$ wide of exact resonance. The orbital properties of these systems indicate that they systematically lie far away from the resonant forced equilibrium. Cumulatively our modeling indicates that while orbital architectures shaped by strong disk damping or tidal dissipation are inconsistent with observations, a scenario where stochastic stirring by turbulent eddies augments the dissipative effects of protoplanetary disks reproduces several features of the data. △ Less

Submitted 3 May, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: Accepted to ApJ

Journal ref: The Astrophysical Journal, Vol. 948 (2023), 12

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

Authors: Fei Dai, Kento Masuda, Corey Beard, Paul Robertson, Max Goldberg, Konstantin Batygin, Luke Bouma, Jack J. Lissauer, Emil Knudstrup, Simon Albrecht, Andrew W. Howard, Heather A. Knutson, Erik A. Petigura, Lauren M. Weiss, Howard Isaacson, Martti Holst Kristiansen, Hugh Osborn, Songhu Wang, Xian-Yu Wang, Aida Behmard, Michael Greklek-McKeon, Shreyas Vissapragada, Natalie M. Batalha, Casey L. Brinkman, Ashley Chontos , et al. (38 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 14 November, 2022; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: 48 pages, 23 figures, 8 tables. Accepted to AAS journals. Comments welcome!

A Criterion for the Stability of Planets in Chains of Resonances

Authors: Max Goldberg, Konstantin Batygin, Alessandro Morbidelli

Abstract: 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 itse… ▽ More 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. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: 8 pages, 4 figures; accepted to Icarus

Analytic Approach to the Late Stages of Giant Planet Formation

Authors: Fred C Adams, Konstantin Batygin

Abstract: 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… ▽ More 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). △ Less

Submitted 20 June, 2022; originally announced June 2022.

Comments: 52 pages, 9 figures, accepted to The Astrophysical Journal

A Self-Consistent Model for Dust-Gas Coupling in Protoplanetary Disks

Authors: Konstantin Batygin, Alessandro Morbidelli

Abstract: Various physical processes that ensue within protoplanetary disks -- including vertical settling of icy/rocky grains, radial drift of solids, planetesimal formation, as well as planetary accretion itself -- are facilitated by hydrodynamic interactions between H/He gas and high-$Z$ dust. The Stokes number, which quantifies the strength of dust-gas coupling, thus plays a central role in protoplaneta… ▽ More Various physical processes that ensue within protoplanetary disks -- including vertical settling of icy/rocky grains, radial drift of solids, planetesimal formation, as well as planetary accretion itself -- are facilitated by hydrodynamic interactions between H/He gas and high-$Z$ dust. The Stokes number, which quantifies the strength of dust-gas coupling, thus plays a central role in protoplanetary disk evolution, and its poor determination constitutes an important source of uncertainty within the theory of planet formation. In this work, we present a simple model for dust-gas coupling, and demonstrate that for a specified combination of the nebular accretion rate, $\dot{M}$, and turbulence parameter, $αあるふぁ$, the radial profile of the Stokes number can be calculated uniquely. Our model indicates that the Stokes number grows sub-linearly with orbital radius, but increases dramatically across the water-ice line. For fiducial protoplanetary disk parameters of $\dot{M}=10^{-8}\,M_{\odot}/$year and $αあるふぁ=10^{-3}$, our theory yields characteristic values of the Stokes number on the order of $\mathrm{St}\sim10^{-4}$ (corresponding to $\sim$mm-sized silicate dust) in the inner nebula and $\mathrm{St}\sim10^{-1}$ (corresponding to $\sim$few-cm-sized icy grains), in the outer regions of the disk. Accordingly, solids are expected to settle into a thin sub-disk at large stellocentric distances, while remaining vertically well-mixed inside the ice line. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 666, A19 (2022)

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

Authors: Michelle Belkovski, Juliette Becker, Alex Howe, Isaac Malsky, Konstantin Batygin

Abstract: 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… ▽ More 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. △ Less

Submitted 31 March, 2022; originally announced March 2022.

Comments: Accepted to The Astronomical Journal

Architectures of Compact Multi-planet Systems: Diversity and Uniformity

Authors: Lauren M. Weiss, Sarah C. Millholland, Erik A. Petigura, Fred C. Adams, Konstantin Batygin, Anthony M. Bloch, Christoph Mordasini

Abstract: One of the most important developments in exoplanet science in the past decade is the discovery of multi-planet systems with sub-Neptune-sized planets interior to 1~AUえーゆー. This chapter explores the architectures of these planetary systems, which often display a remarkable degree of uniformity: the planets have nearly equal sizes, regular orbital spacing, low eccentricities, and small mutual inclinati… ▽ More One of the most important developments in exoplanet science in the past decade is the discovery of multi-planet systems with sub-Neptune-sized planets interior to 1~AUえーゆー. This chapter explores the architectures of these planetary systems, which often display a remarkable degree of uniformity: the planets have nearly equal sizes, regular orbital spacing, low eccentricities, and small mutual inclinations. This uniformity stands in sharp contrast to the diverse nature of the exoplanet sample considered as a whole (as well as our inner solar system). We begin with a critical review of the observations -- including possible biases -- and find that these peas-in-a-pod planetary systems are apparently a common outcome of the planet formation process. Modest departures from exact uniformity suggest additional patterns, such as the planet mass slowly increasing with semi-major axis. The star formation process naturally produces circumstellar disks with the properties required to produce these planetary systems, although the solid material must move inward from its initial location. We discuss primary modes of planetary assembly, the role of orbital migration, and post-nebular atmospheric loss. Mature planetary systems are found to be near their minimum energy (tidal equilibrium) configurations; this finding provides a partial explanation for their observed properties and indicates that efficient energy dissipation must occur. Finally, we consider population synthesis models and show that peas-in-a-pod patterns emerge with reasonable choices for the input parameters. Nonetheless, interesting observational and theoretical challenges remain in order to understand how these surprisingly organized planetary systems arise from the disorder of their formation processes. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: Chapter for publication in Protostars and Planets VII, 42 pages, 17 Figures

Generating Stellar Obliquity in Systems with Broken Protoplanetary Disks

Authors: Marguerite Epstein-Martin, Juliette Becker, Konstantin Batygin

Abstract: 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 obli… ▽ More 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. △ Less

Submitted 8 March, 2022; originally announced March 2022.

Comments: 9 pages, 4 figures, accepted for publication in ApJ

Architectures of Compact Super-Earth Systems Shaped by Instabilities

Authors: Max Goldberg, Konstantin Batygin

Abstract: 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… ▽ More 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. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: 14 pages, 8 figures, accepted for publication in AJ

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

Authors: Henrik Knierim, Konstantin Batygin, Bertram Bitsch

Abstract: 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… ▽ More 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. △ Less

Submitted 10 February, 2022; v1 submitted 20 January, 2022; originally announced January 2022.

Comments: published in A&A, 5 pages, including 5 figures, plus 2 pages appendix, including 2 figures; v2 follows the A&A language guidelines

Journal ref: A&A 658, L7 (2022)

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

Authors: Alessandro Morbidelli, Kevin Baillie, Konstantin Batygin, Sebastien Charnoz, Tristan Guillot, David C. Rubie, Thorsten Kleine

Abstract: 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 unreal… ▽ More 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. △ Less

Submitted 31 December, 2021; originally announced December 2021.

Comments: Published online in Nature Astronomy on Dec. 22, 2021

Migration of Jupiter mass planets in discs with laminar accretion flows

Authors: E. Lega, A. Morbidelli, R. P. Nelson, X. S. Ramos, A. Crida, W. Bethune, K. Batygin

Abstract: Migration of giant planets in discs with low viscosity has been studied recently. The proportionality between migration speed and the disc's viscosity is broken by the presence of vortices that appear at the edges of the planet-induced gap. Our goal is to investigate vortex-driven migration in low-viscosity discs in the presence of radial advection of gas, as expected from angular momentum removal… ▽ More Migration of giant planets in discs with low viscosity has been studied recently. The proportionality between migration speed and the disc's viscosity is broken by the presence of vortices that appear at the edges of the planet-induced gap. Our goal is to investigate vortex-driven migration in low-viscosity discs in the presence of radial advection of gas, as expected from angular momentum removal by magnetised disc winds. We performed three dimensional simulations using the grid-based code FARGOCA. We mimicked the effects of a disc wind by applying a synthetic torque on a surface layer of the disc characterised by a prescribed column density Sigma_A so that it results in a disc accretion rate of 10^-8 Solar masses per year. Discs with this structure are called 'layered' and the layer where the torque is applied is denoted as 'active'. We also consider the case of accretion focussed near the disc midplane to mimic transport properties induced by a large Hall effect or by weak Ohmic diffusion. We observe two migration phases: in the first phase, the migration of the planet is driven by the vortex and is directed inwards. This phase ends when the vortex disappear. Migration depends on the ability of the torque from the planet to block the accretion flow. When the flow is fast and unimpeded (small Sigma_A) migration is very slow. When the accretion flow is completely blocked, migration is faster and the speed is controlled by the rate at which the accretion flow refills the gap behind the migrating planet. The migration speed of a giant planet in a layered protoplanetary disc depends on the thickness of the accreting layer. The lack of large-scale migration apparently experienced by the majority of giant exoplanets can be explained if the accreting layer is sufficiently thin to allow unimpeded accretion through the disc. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: 14 pages, 11 figures

Journal ref: A&A 658, A32 (2022)

The Stability Boundary of the Distant Scattered Disk

Authors: Konstantin Batygin, Rosemary A. Mardling, David Nesvorny

Abstract: The scattered disk is a vast population of trans-Neptunian minor bodies that orbit the sun on highly elongated, long-period orbits. The stability of scattered disk objects is primarily controlled by a single parameter - their perihelion distance. While the existence of a perihelion boundary that separates chaotic and regular motion of long-period orbits is well established through numerical experi… ▽ More The scattered disk is a vast population of trans-Neptunian minor bodies that orbit the sun on highly elongated, long-period orbits. The stability of scattered disk objects is primarily controlled by a single parameter - their perihelion distance. While the existence of a perihelion boundary that separates chaotic and regular motion of long-period orbits is well established through numerical experiments, its theoretical basis as well as its semi-major axis dependence remain poorly understood. In this work, we outline an analytical model for the dynamics of distant trans-Neptunian objects and show that the orbital architecture of the scattered disk is shaped by an infinite chain of $2:j$ resonances with Neptune. The widths of these resonances increase as the perihelion distance approaches Neptune's semi-major axis, and their overlap drives chaotic motion. Within the context of this picture, we derive an analytic criterion for instability of long-period orbits, and demonstrate that rapid dynamical chaos ensues when the perihelion drops below a critical value, given by $q_{\rm{crit}}=a_{\rm{N}}\,\big(\ln((24^2/5)\,(m_{\rm{N}}/M_{\odot})\,(a/a_{\rm{N}})^{5/2})\big)^{1/2}$. This expression constitutes a boundary between the "detached" and actively "scattering" sub-populations of distant trans-Neptunian minor bodies. Additionally, we find that within the stochastic layer, the Lyapunov time of scattered disk objects approaches the orbital period, and show that the semi-major axis diffusion coefficient is approximated by $\mathcal{D}_a\sim(8/(5\,πぱい))\,(m_{\rm{N}}/M_{\odot})\,\sqrt{\mathcal{G}\,M_{\odot}\,a_{\rm{N}}}\,\exp\big[-(q/a_{\rm{N}})^2/2\big]$. We confirm our results with numerical simulations and highlight the connections between scattered disk dynamics and the Chirikov Standard Map. Implications of our results for the long-term evolution of the distant solar system are discussed. △ Less

Submitted 30 October, 2021; originally announced November 2021.

Comments: 14 pages, 3 figures, published in ApJ

Journal ref: The Astrophysical Journal, 2021, Volume 920, Number 2

The Kepler-223 resonance holds information on turbulence during the gas disk phase

Authors: L. -A. Hühn, G. Pichierri, B. Bitsch, K. Batygin

Abstract: Many fundamental physical processes regarding planetary formation in protoplanetary disks are still imperfectly understood, with an elusive phenomenon being turbulence in such disks. Observations are available of planetary systems and of some protoplanetary disks, which can serve as a starting point for these investigations. Detected systems reveal different architectures of planets. One particula… ▽ More Many fundamental physical processes regarding planetary formation in protoplanetary disks are still imperfectly understood, with an elusive phenomenon being turbulence in such disks. Observations are available of planetary systems and of some protoplanetary disks, which can serve as a starting point for these investigations. Detected systems reveal different architectures of planets. One particularly interesting case to consider is the Kepler-223 system, which contains a rare configuration of four planets in a resonance chain, implying a certain migration history. We aim to use the orbital configuration of Kepler-223's planets to constrain the parameters of the protoplanetary disk that allow for the formation of a chain of mean-motion resonances that resembles the one of Kepler-223. The parameters we aim to investigate are primarily the disk viscosity and surface density. We use the swift_symba N-body integrator with additional dissipative forces to mimic planet-disk interactions. We constrain the surface densities and viscosities that allow the formation of a resonant chain like that of Kepler-223. We find that surface densities of up to a few minimum mass solar nebula (MMSN) surface densities and disk viscosity parameters $αあるふぁ$ of a few $10^{-3}$ up to $10^{-2}$ are most successful at reproducing the architecture of this particular planetary system. We describe how these two quantities are linked to each other considering the success of reproducing the chain and find that higher disk surface densities in turn require lower viscosities to build the chain. Our results show that well characterized observed planetary systems hold information about their formation conditions in the protoplanetary disks and that it is possible to extract this information, namely the initial disk surface density and viscosity, helping to constrain planet formation. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 12 pages, 7 figues, accepted for publication in A&A

Journal ref: A&A 656, A115 (2021)

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

Authors: Michael E. Brown, Konstantin Batygin

Abstract: 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 aster… ▽ More 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. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: Astronomical Journal, in press

On the Fate of Interstellar Objects Captured by our Solar System

Authors: Kevin J Napier, Fred C Adams, Konstantin Batygin

Abstract: With the recent discoveries of interstellar objects `Oumuamua and Borisov traversing the solar system, understanding the dynamics of interstellar objects is more pressing than ever. These detections have highlighted the possibility that captured interstellar material could be trapped in our solar system. The first step in rigorously investigating this question is to calculate a capture cross secti… ▽ More With the recent discoveries of interstellar objects `Oumuamua and Borisov traversing the solar system, understanding the dynamics of interstellar objects is more pressing than ever. These detections have highlighted the possibility that captured interstellar material could be trapped in our solar system. The first step in rigorously investigating this question is to calculate a capture cross section for interstellar objects as a function of hyperbolic excess velocity, which can be convolved with any velocity dispersion to compute a capture rate (Napier et. al. 2021). Although the cross section provides the first step toward calculating the mass of alien rocks residing in our solar system, we also need to know the lifetime of captured objects. We use an ensemble of N-body simulations to characterize a dynamical lifetime for captured interstellar objects and determines the fraction of surviving objects as a function of time (since capture). We also illuminate the primary effects driving their secular evolution. Finally, we use the resulting dynamical lifetime function to estimate the current inventory of captured interstellar material in the solar system. We find that capture from the field yields a steady state mass of only $\sim 10^{-13} M_{\oplus}$, whereas the mass remaining from capture events in the birth cluster is roughly $10^{-9} M_{\oplus}$. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: Accepted for publication in the Planetary Science Journal. 12 pages, 5 figures

The orbit of Planet Nine

Authors: Michael E. Brown, Konstantin Batygin

Abstract: The existence of a giant planet beyond Neptune -- referred to as Planet Nine (P9) -- has been inferred from the clustering of longitude of perihelion and pole position of distant eccentric Kuiper belt objects (KBOs). After updating calculations of observational biases, we find that the clustering remains significant at the 99.6\% confidence level. We thus use these observations to determine orbita… ▽ More The existence of a giant planet beyond Neptune -- referred to as Planet Nine (P9) -- has been inferred from the clustering of longitude of perihelion and pole position of distant eccentric Kuiper belt objects (KBOs). After updating calculations of observational biases, we find that the clustering remains significant at the 99.6\% confidence level. We thus use these observations to determine orbital elements of P9. A suite of numerical simulations shows that the orbital distribution of the distant KBOs is strongly influenced by the mass and orbital elements of P9 and thus can be used to infer these parameters. Combining the biases with these numerical simulations, we calculate likelihood values for discrete set of P9 parameters, which we then use as input into a Gaussian Process emulator that allows a likelihood computation for arbitrary values of all parameters. We use this emulator in a Markov Chain Monte Carlo analysis to estimate parameters of P9. We find a P9 mass of $6.2^{+2.2}_{-1.3}$ Earth masses, semimajor axis of $380^{+140}_{-80}$ AUえーゆー, inclination of $16\pm5^\circ$ and perihelion of $300^{+85}_{-60}$ AUえーゆー. Using samples of the orbital elements and estimates of the radius and albedo of such a planet, we calculate the probability distribution function of the on-sky position of Planet Nine and of its brightness. For many reasonable assumptions, Planet Nine is closer and brighter than initially expected, though the probability distribution includes a long tail to larger distances, and uncertainties in the radius and albedo of Planet Nine could yield fainter objects. △ Less

Submitted 26 August, 2021; v1 submitted 22 August, 2021; originally announced August 2021.

Comments: AJ, in press. Previous version had outdated version of Figure 9 affecting the predicted magnitudes. New Figure 9 and fixed typos

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

Authors: Jack Lubin, Judah Van Zandt, Rae Holcomb, Lauren M. Weiss, Erik A Petigura, Paul Robertson, Joseph M. Akana Murphy, Nicholas Scarsdale, Konstantin Batygin, Alex S. Polanski, Natalie M. Batalha, Ian J. M. Crossfield, Courtney Dressing, Benjamin Fulton, Andrew W. Howard, Daniel Huber, Howard Isaacson, Stephen R. Kane, Arpita Roy, Corey Beard, Sarah Blunt, Ashley Chontos, Fei Dai, Paul A. Dalba, Kaz Gary , et al. (10 additional authors not shown)

Abstract: 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 meas… ▽ More 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. △ Less

Submitted 3 January, 2022; v1 submitted 4 August, 2021; originally announced August 2021.

Comments: 20 pages, 10 figures. Accepted for publication in The Astronomical Journal

Collisional Growth Within the Solar System's Primordial Planetesimal Disk and the Timing of the Giant Planet Instability

Authors: Marvin Morgan, Darryl Seligman, Konstantin Batygin

Abstract: The large scale structure of the Solar System has been shaped by a transient dynamical instability that may have been triggered by the interaction of the giants planets with a massive primordial disk of icy debris. In this work, we investigate the conditions under which this primordial disk could have coalesced into planets using analytic and numerical calculations. In particular, we perform numer… ▽ More The large scale structure of the Solar System has been shaped by a transient dynamical instability that may have been triggered by the interaction of the giants planets with a massive primordial disk of icy debris. In this work, we investigate the conditions under which this primordial disk could have coalesced into planets using analytic and numerical calculations. In particular, we perform numerical simulations of the Solar System's early dynamical evolution that account for the viscous stirring and collisional damping within the disk. We demonstrate that if collisional damping would have been sufficient to maintain a temperate velocity dispersion, Earth mass trans-Neptunian planets could have emerged within a timescale of 10 Myr. Therefore, our results favor a scenario wherein the dynamical instability of the outer Solar System began immediately upon the dissipation of the gaseous nebula to avoid the overproduction of Earth mass planets in the outer Solar System. △ Less

Submitted 21 July, 2021; originally announced July 2021.

Comments: 8 pages, 4 figures, 1 table, accepted for publication in ApJL

Migrating Planets into Ultra-Short-Period Orbits during Episodic Accretion Events

Authors: Juliette Becker, Konstantin Batygin, Fred Adams

Abstract: Ultra-short-period (USP) planets reside inside the expected truncation radius for typical T Tauri disks. As a result, their current orbital locations require an explanation beyond standard disk migration or in situ formation. Modern theories of planet-disk interactions indicate that once a planet migrates close to the disk's inner truncation radius, Type I torques vanish or switch direction, depen… ▽ More Ultra-short-period (USP) planets reside inside the expected truncation radius for typical T Tauri disks. As a result, their current orbital locations require an explanation beyond standard disk migration or in situ formation. Modern theories of planet-disk interactions indicate that once a planet migrates close to the disk's inner truncation radius, Type I torques vanish or switch direction, depending on the stellar and disk conditions, so that the planet is expected to stop its orbital decay and become trapped. In this work, we show that that magnetically-driven sub-Keplerian gas flow in the inner disk can naturally counteract these effects and produce systems with USP planets at their observed orbital radii. The sub-Keplerian gas flow provides a headwind to small planets, and the resulting torque can overcome the effects of outward Type I migration near the co-rotation radius. For suitable disk and planet parameters, the torques due to the sub-Keplerian gas flow lead to inward migration on a rapid timescale. Over the time span of an FU Ori outburst, which moves the disk truncation radius inward, the rapid headwind migration can place planets in USP orbits. The combination of headwind migration and FU Ori outbursts thus provides a plausible mechanism to move small planets from $a=0.05-0.1$ AUえーゆー down to $a=0.01-0.02$ AUえーゆー. This effect is amplified for low-mass planets, consistent with existing observations. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Comments: accepted to ApJ; 13 pages, 6 figures

The distribution of mutual inclinations arising from the stellar quadrupole moment

Authors: Kathleen Schultz, Christopher Spalding, Konstantin Batygin

Abstract: A large proportion of transiting planetary systems appear to possess only a single planet as opposed to multiple transiting planets. This excess of singles is indicative of significant mutual inclinations existing within a large number of planetary systems, but the origin of these misalignments is unclear. Moreover, recent observational characterization reveals that mutual inclinations tend to inc… ▽ More A large proportion of transiting planetary systems appear to possess only a single planet as opposed to multiple transiting planets. This excess of singles is indicative of significant mutual inclinations existing within a large number of planetary systems, but the origin of these misalignments is unclear. Moreover, recent observational characterization reveals that mutual inclinations tend to increase with proximity to the host star. These trends are both consistent with the dynamical influence of a strong quadrupolar potential arising from the host star during its early phase of rapid rotation, coupled with a non-zero stellar obliquity. Here, we simulate a population of planetary systems subject to the secular perturbation arising from a tilted, oblate host star as it contracts and spins down subsequent to planet formation. We demonstrate that this mechanism can reproduce the general increase in planet-planet mutual inclinations with proximity to the host star, and delineate a parameter space wherein the host star can drive dynamical instabilities. We suggest that approximately 5-10\% of low-mass Kepler systems are susceptible to this instability mechanism, suggesting that a significant number of single-transiting planets may truly be intrinsically single. We also report a novel connection between instability and stellar obliquity reduction and make predictions that can be tested within upcoming TESS observations. △ Less

Submitted 30 June, 2021; originally announced July 2021.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society, 16 pages, 9 figures

A Tidal Origin for a 3-body Resonance in Kepler-221

Authors: Max Goldberg, Konstantin Batygin

Abstract: Over the course of the last two decades, traditional models of planet formation have been repeatedly challenged by the emerging census of extrasolar planets. Key among them is the orbital architecture problem: while standard models of orbital migration predict resonant orbits for short-period objects, most planets do not appear to lie in orbital resonances. Here we show that the four-planet system… ▽ More Over the course of the last two decades, traditional models of planet formation have been repeatedly challenged by the emerging census of extrasolar planets. Key among them is the orbital architecture problem: while standard models of orbital migration predict resonant orbits for short-period objects, most planets do not appear to lie in orbital resonances. Here we show that the four-planet system Kepler-221, not previously recognized to have active orbital resonances, has a three-body commensurability relation unique within the Kepler sample. Using a suite of numerical experiments as well as a perturbative analysis, we demonstrate that this system likely began as a resonant chain and proceeded to undergo large-scale divergence away from resonance, under the action of tidal dissipation. Our results further indicate that obliquity tides, driven by a secular spin-orbit resonance and mutual inclination, are an excellent candidate for driving this orbital divergence, and that the high tidal luminosity may also explain the anomalous size of planet b, which lies within the Fulton radius gap. △ Less

Submitted 16 May, 2021; originally announced May 2021.

Comments: 15 pages, 9 figures, accepted for publication in AJ

Injection of Inner Oort Cloud Objects Into the Distant Kuiper Belt by Planet Nine

Authors: Konstantin Batygin, Michael E. Brown

Abstract: The outer solar system exhibits an anomalous pattern of orbital clustering, characterized by an approximate alignment of the apsidal lines and angular momentum vectors of distant, long-term stable Kuiper belt objects. One explanation for this dynamical confinement is the existence of a yet-undetected planetary-mass object, "Planet Nine (P9)". Previous work has shown that trans-Neptunian objects, w… ▽ More The outer solar system exhibits an anomalous pattern of orbital clustering, characterized by an approximate alignment of the apsidal lines and angular momentum vectors of distant, long-term stable Kuiper belt objects. One explanation for this dynamical confinement is the existence of a yet-undetected planetary-mass object, "Planet Nine (P9)". Previous work has shown that trans-Neptunian objects, which originate within the scattered disk population of the Kuiper belt, can be corralled into orbital alignment by Planet Nine's gravity over ~Gyr timescales, and characteristic P9 parameters have been derived by matching the properties of a synthetic Kuiper belt generated within numerical simulations to the available observational data. In this work, we show that an additional dynamical process is in play within the framework of the Planet Nine hypothesis, and demonstrate that P9-induced dynamical evolution facilitates orbital variations within the otherwise dynamically frozen inner Oort cloud. As a result of this evolution, inner Oort cloud bodies can acquire orbits characteristic of the distant scattered disk, implying that if Planet Nine exists, the observed census of long-period trans-Neptunian objects is comprised of a mixture of Oort cloud and Kuiper belt objects. Our simulations further show that although inward-injected inner Oort cloud objects exhibit P9-driven orbital confinement, the degree of clustering is weaker than that of objects originating within the Kuiper belt. Cumulatively, our results suggest that a more eccentric Planet Nine is likely necessary to explain the data than previously thought. △ Less

Submitted 12 April, 2021; originally announced April 2021.

Comments: 11 pages, 1 figures, accepted for publication in ApJL

Journal ref: The Astrophysical Journal Letters, April 2021, Volume 910, Number 2

The Onset of Chaos in Permanently Deformed Binaries from Spin-Orbit and Spin-Spin Coupling

Authors: Darryl Seligman, Konstantin Batygin

Abstract: Permanently deformed objects in binary systems can experience complex rotation evolution, arising from the extensively studied effect of spin-orbit coupling as well as more nuanced dynamics arising from spin-spin interactions. The ability of an object to sustain an aspheroidal shape largely determines whether or not it will exhibit non-trivial rotational behavior. In this work, we adopt a simplifi… ▽ More Permanently deformed objects in binary systems can experience complex rotation evolution, arising from the extensively studied effect of spin-orbit coupling as well as more nuanced dynamics arising from spin-spin interactions. The ability of an object to sustain an aspheroidal shape largely determines whether or not it will exhibit non-trivial rotational behavior. In this work, we adopt a simplified model of a gravitationally interacting primary and satellite pair, where each body's quadrupole moment is approximated by two diametrically opposed point masses. After calculating the net gravitational torque on the satellite from the primary, and the associated equations of motion, we employ a Hamiltonian formalism which allows for a perturbative treatment of the spin-orbit and retrograde and prograde spin-spin coupling states. By analyzing the resonances individually and collectively, we determine the criteria for resonance overlap and the onset of chaos, as a function of orbital and geometric properties of the binary. We extend the 2D planar geometry to calculate the obliquity evolution, and find that satellites in spin-spin resonances undergo precession when inclined out of the plane, but do not tumble. We apply our resonance overlap criteria to the contact binary system (216) Kleopatra, and find that its satellites, Cleoselene and Alexhelios, may plausibly be exhibiting chaotic rotational dynamics from the overlap of the spin-orbit and retrograde spin-spin resonances. While this model is by construction generalizable to any binary system, it will be particularly useful to study small bodies in the solar system, whose irregular shapes make them ideal candidates for exotic rotational states. △ Less

Submitted 23 March, 2021; originally announced March 2021.

Comments: 26 pages, 9 figures, resubmitted to ApJ following minor revisions

On the Capture of Interstellar Objects by our Solar System

Authors: K. J. Napier, F. C. Adams, K. Batygin

Abstract: Motivated by recent visits from interstellar comets, along with continuing discoveries of minor bodies in orbit of the Sun, this paper studies the capture of objects on initially hyperbolic orbits by our solar system. Using an ensemble of $\sim500$ million numerical experiments, this work generalizes previous treatments by calculating the capture cross section as a function of asymptotic speed. Th… ▽ More Motivated by recent visits from interstellar comets, along with continuing discoveries of minor bodies in orbit of the Sun, this paper studies the capture of objects on initially hyperbolic orbits by our solar system. Using an ensemble of $\sim500$ million numerical experiments, this work generalizes previous treatments by calculating the capture cross section as a function of asymptotic speed. The resulting velocity-dependent cross section can then be convolved with any distribution of relative speeds to determine the capture rate for incoming bodies. This convolution is carried out for the usual Maxwellian distribution, as well as the velocity distribution expected for rocky debris ejected from planetary systems. We also construct an analytic description of the capture process that provides an explanation for the functional form of the capture cross section in both the high velocity and low velocity limits. △ Less

Submitted 16 February, 2021; originally announced February 2021.

Comments: 18 pages, 5 figures. Accepted to the Planetary Science Journal

Eigenstates of Quasi-Keplerian Self-Gravitating Particle Discs

Authors: Walker Melton, Konstantin Batygin

Abstract: 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 anal… ▽ More 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. △ Less

Submitted 3 February, 2021; originally announced February 2021.

Comments: 9 pages, 4 figures

The Origin of Systems of Tightly Packed Inner Planets with Misaligned, Ultra-Short-Period Companions

Authors: Juliette Becker, Konstantin Batygin, Daniel C. Fabrycky, Fred C. Adams, Gongjie Li, Andrew Vanderburg, Joseph E. Rodriguez

Abstract: Ultra-short period planets provide a window into the inner edge of the parameter space occupied by planetary orbits. In one particularly intriguing class of multi-planet systems, the ultra-short period planet is flanked by short-period companions, and the outer planets occupy a discernibly distinct dynamical state. In the observational database, this phenomenon is represented by a small number of… ▽ More Ultra-short period planets provide a window into the inner edge of the parameter space occupied by planetary orbits. In one particularly intriguing class of multi-planet systems, the ultra-short period planet is flanked by short-period companions, and the outer planets occupy a discernibly distinct dynamical state. In the observational database, this phenomenon is represented by a small number of stars hosting systems of tightly packed co-planar planets as well as an ultra-short period planet, whose orbit of is misaligned relative to the mutual plane of the former. In this work, we explore two different mechanisms that can produce an ultra-short period planet that is misaligned with the rest of its compact planetary system: natural decoupling between the inner and outer system via the stellar quadrupole moment, and decoupling forced by an external companion with finely-tuned orbital parameters. These two processes operate with different timescales, and can thus occur simultaneously. In this work, we use the K2-266 system as an illustrative example to elucidate the dynamics of these two processes, and highlight the types of constraints that may arise regarding the dynamical histories of systems hosting ultra-short period planets. △ Less

Submitted 22 September, 2020; originally announced September 2020.

Comments: Accepted to AJ

The Great Inequality and the Dynamical Disintegration of the Outer Solar System

Authors: Jon K. Zink, Konstantin Batygin, Fred C. Adams

Abstract: Using an ensemble of N-body simulations, this paper considers the fate of the outer gas giants (Jupiter, Saturn, Uranus, and Neptune) after the Sun leaves the main sequence and completes its stellar evolution. Due to solar mass-loss -- which is expected to remove roughly half of the star's mass -- the orbits of the giant planets expand. This adiabatic process maintains the orbital period ratios, b… ▽ More Using an ensemble of N-body simulations, this paper considers the fate of the outer gas giants (Jupiter, Saturn, Uranus, and Neptune) after the Sun leaves the main sequence and completes its stellar evolution. Due to solar mass-loss -- which is expected to remove roughly half of the star's mass -- the orbits of the giant planets expand. This adiabatic process maintains the orbital period ratios, but the mutual interactions between planets and the width of mean-motion resonances (MMR) increase, leading to the capture of Jupiter and Saturn into a stable 5:2 resonant configuration. The expanded orbits, coupled with the large-amplitude librations of the critical MMR angle, make the system more susceptible to perturbations from stellar flyby interactions. Accordingly, within about 30 Gyr, stellar encounters perturb the planets onto the chaotic sub-domain of the 5:2 resonance, triggering a large-scale instability, which culminates in the ejections of all but one planet over the subsequent $\sim10$ Gyr. After an additional $\sim50$ Gyr, a close stellar encounter (with a perihelion distance less than $\sim200$ AUえーゆー) liberates the final planet. Through this sequence of events, the characteristic timescale over which the solar system will be completely dissolved is roughly 100 Gyr. Our analysis thus indicates that the expected dynamical lifetime of the solar system is much longer than the current age of the universe, but is significantly shorter than previous estimates. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Comments: 12 pages, 5 figures; accepted for publication in AJ

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

Authors: Tobias Köhne, Konstantin Batygin

Abstract: 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 producti… ▽ More 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. △ Less

Submitted 25 August, 2020; originally announced August 2020.

Comments: Accepted to Celestial Mechanics and Dynamical Astronomy's Topical Collection on Trans-Neptunian Objects (11 pages, 4 figures, 1 table)

Journal ref: Celest Mech Dyn Astr 132, 44 (2020)

Interactions Among Non-Interacting Particles in Planet Formation Simulations

Authors: Shirui Peng, Konstantin Batygin

Abstract: 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… ▽ More 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. △ Less

Submitted 22 July, 2020; originally announced July 2020.

No evidence for interstellar planetesimals trapped in the Solar System

Authors: A. Morbidelli, K. Batygin, R. Brasser, S. Raymond

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

Submitted 8 June, 2020; originally announced June 2020.

Comments: in press in MNRAS

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

Authors: Sarah Millholland, Erik Petigura, Konstantin Batygin

Abstract: 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 diversi… ▽ More 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. △ Less

Submitted 22 May, 2020; originally announced May 2020.

Comments: Accepted for publication in ApJ, 17 pages, 8 figures. Code available at this URL: https://github.com/smillholland/Sub-Saturns/

Formation of Giant Planet Satellites

Authors: Konstantin Batygin, Alessandro Morbidelli

Abstract: 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 satell… ▽ More 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. △ Less

Submitted 17 May, 2020; originally announced May 2020.

Comments: 31 pages, 12 figures, accepted for publication in ApJ

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

Authors: Brendan P. Bowler, Yifan Zhou, Caroline V. Morley, Tiffany Kataria, Marta L. Bryan, Björn Benneke, Konstantin Batygin

Abstract: 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 mass… ▽ More 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. △ Less

Submitted 10 April, 2020; originally announced April 2020.

Comments: Accepted to ApJ Letters