Showing 1–39 of 39 results for author: Millholland, S

Resonant sub-Neptunes are puffier

Authors: Adrien Leleu, Jean-Baptiste Delisle, Remo Burn, André Izidoro, Stéphane Udry, Xavier Dumusque, Christophe Lovis, Sarah Millholland, Léna Parc, François Bouchy, Vincent Bourrier, Yann Alibert, João Faria, Christoph Mordasini, Damien Ségransan

Abstract: A systematic, population-level discrepancy exists between the densities of exoplanets whose masses have been measured with transit timing variations (TTVs) versus those measured with radial velocities (RVs). Since the TTV planets are predominantly nearly resonant, it is still unclear whether the discrepancy is attributed to detection biases or to astrophysical differences between the nearly resona… ▽ More A systematic, population-level discrepancy exists between the densities of exoplanets whose masses have been measured with transit timing variations (TTVs) versus those measured with radial velocities (RVs). Since the TTV planets are predominantly nearly resonant, it is still unclear whether the discrepancy is attributed to detection biases or to astrophysical differences between the nearly resonant and non resonant planet populations. We defined a controlled, unbiased sample of 36 sub-Neptunes characterised by Kepler, TESS, HARPS, and ESPRESSO. We found that their density depends mostly on the resonant state of the system, with a low probability (of $0.002_{-0.001}^{+0.010}$) that the mass of (nearly) resonant planets is drawn from the same underlying population as the bulk of sub-Neptunes. Increasing the sample to 133 sub-Neptunes reveals finer details: the densities of resonant planets are similar and lower than non-resonant planets, and both the mean and spread in density increase for planets that are away from resonance. This trend is also present in RV-characterised planets alone. In addition, TTVs and RVs have consistent density distributions for a given distance to resonance. We also show that systems closer to resonances tend to be more co-planar than their spread-out counterparts. These observational trends are also found in synthetic populations, where planets that survived in their original resonant configuration retain a lower density; whereas less compact systems have undergone post-disc giant collisions that increased the planet's density, while expanding their orbits. Our findings reinforce the claim that resonant systems are archetypes of planetary systems at their birth. △ Less

Submitted 27 June, 2024; originally announced June 2024.

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

Planet-Planet Scattering and ZLK Migration -- The Dynamical History of HAT-P-11

Authors: Tiger Lu, Qier An, Gongjie Li, Sarah C. Millholland, G. Mirek Brandt, Timothy D. Brandt

Abstract: The two planets of the HAT-P-11 system represent fascinating dynamical puzzles due to their significant eccentricities and orbital misalignments. In particular, HAT-P-11 b is on a close-in orbit that tides should have circularized well within the age of the system. Here we propose a two-step dynamical process that can reproduce all intriguing aspects of the system. We first invoke planet-planet sc… ▽ More The two planets of the HAT-P-11 system represent fascinating dynamical puzzles due to their significant eccentricities and orbital misalignments. In particular, HAT-P-11 b is on a close-in orbit that tides should have circularized well within the age of the system. Here we propose a two-step dynamical process that can reproduce all intriguing aspects of the system. We first invoke planet-planet scattering to generate significant eccentricities and mutual inclinations between the planets. We then propose that this misalignment initiated von-Zeipel-Lidov-Kozai cycles and high-eccentricity migration that ultimately brought HAT-P-11 b to its present-day orbit. We find that this scenario is fully consistent only when significant tidally-driven radius inflation is accounted for during the tidal migration. We present a suite of N-body simulations exploring each phase of evolution and show that this scenario is consistent with all observational posteriors and the reported age of the system. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 17 pages, 8 figures, submitted to ApJ

TOI-1685 b is a Hot Rocky Super-Earth: Updates to the Stellar and Planet Parameters of a Popular JWST Cycle 2 Target

Authors: Jennifer A. Burt, Matthew J. Hooton, Eric E. Mamajek, Oscar Barragán, Sarah C. Millholland, Tyler R. Fairnington, Chloe Fisher, Samuel P. Halverson, Chelsea X. Huang, Madison Brady, Andreas Seifahrt, Eric Gaidos, Rafael Luque, David Kasper, Jacob L. Bean

Abstract: We present an updated characterization of the TOI-1685 planetary system, which consists of a P$_{\rm{b}}$ = 0.69\,day USP super-Earth planet orbiting a nearby ($d$ = 37.6\,pc) M2.5V star (TIC 28900646, 2MASS J04342248+4302148). This planet was previously featured in two contemporaneous discovery papers, but the best-fit planet mass, radius, and bulk density values were discrepant allowing it to be… ▽ More We present an updated characterization of the TOI-1685 planetary system, which consists of a P$_{\rm{b}}$ = 0.69\,day USP super-Earth planet orbiting a nearby ($d$ = 37.6\,pc) M2.5V star (TIC 28900646, 2MASS J04342248+4302148). This planet was previously featured in two contemporaneous discovery papers, but the best-fit planet mass, radius, and bulk density values were discrepant allowing it to be interpreted either as a hot, bare rock or a 50\% H$_{2}$O / 50\% MgSiO$_{3}$ water world. TOI-1685 b will be observed in three independent JWST cycle two programs, two of which assume the planet is a water world while the third assumes that it is a hot rocky planet. Here we include a refined stellar classification with a focus on addressing the host star's metallicity, an updated planet radius measurement that includes two sectors of TESS data and multi-color photometry from a variety of ground-based facilities, and a more accurate dynamical mass measurement from a combined CARMENES, IRD, and MAROON-X radial velocity data set. We find that the star is very metal-rich ([Fe/H] $\simeq$ +0.3) and that the planet is systematically smaller, lower mass, and higher density than initially reported, with new best-fit parameters of \Rpl = 1.468 $^{+0.050}_{-0.051}$ \Rearth\ and \Mpl = 3.03$^{+0.33}_{-0.32}$ \Mearth. These results fall in between the previously derived values and suggest that TOI-1685 b is a hot, rocky, planet with an Earth-like density (\Rhopl = 5.3 $\pm$ 0.8 g cm$^{-3}$, or 0.96 \rhoearth), high equilibrium temperature (T$_{\rm{eq}}$ = 1062 $\pm$ 27 K) and negligible volatiles, rather than a water world. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 20 pages, 9 Figures, accepted for publication in ApJL. Datasets and software available via Zenodo and GitHub links found in the paper

The black hole low mass X-ray binary V404 Cygni is part of a wide hierarchical triple, and formed without a kick

Authors: Kevin B. Burdge, Kareem El-Badry, Erin Kara, Claude Canizares, Deepto Chakrabarty, Anna Frebel, Sarah C. Millholland, Saul Rappaport, Rob Simcoe, Andrew Vanderburg

Abstract: Evidence suggests that when compact objects such as black holes and neutron stars form, they may receive a ``natal kick,'' where the stellar remnant gains momentum. Observational evidence for neutron star kicks is substantial, yet limited for black hole natal kicks, and some proposed black hole formation scenarios result in very small kicks. Here, we report the discovery that the canonical black h… ▽ More Evidence suggests that when compact objects such as black holes and neutron stars form, they may receive a ``natal kick,'' where the stellar remnant gains momentum. Observational evidence for neutron star kicks is substantial, yet limited for black hole natal kicks, and some proposed black hole formation scenarios result in very small kicks. Here, we report the discovery that the canonical black hole low-mass X-ray binary V404 Cygni is part of a wide hierarchical triple with a tertiary companion at least 3500 astronomical units away from the inner binary. Given the orbital configuration, the black hole likely received a sub-5 kilometer per second kick to have avoided unbinding the tertiary. This discovery reveals that at least some black holes form with nearly no natal kick. Furthermore, the tertiary in this system lends credence to evolutionary models of low-mass X-ray binaries involving a hierarchical triple structure. Remarkably, the tertiary is evolved, indicating that the system formed 3-5 billion years ago, and that the black hole has removed at least half a solar mass of matter from its evolved secondary companion. During the event in which the black hole formed, it is likely that at least half of the mass of the black hole progenitor collapsed into the black hole; it may even have undergone a complete implosion, enabling the tertiary to remain loosely bound. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 36 pages, 10 figures, 1 table

Spin Dynamics of Planets in Resonant Chains

Authors: Sarah C. Millholland, Teo Lara, Jan Toomlaid

Abstract: About a dozen exoplanetary systems have been discovered with three or more planets participating in a sequence of mean-motion resonances. The unique and complex architectures of these so-called "resonant chains" motivate efforts to characterize their planets holistically. In this work, we perform a comprehensive exploration of the spin-axis dynamics of planets in resonant chains. Planetary spin st… ▽ More About a dozen exoplanetary systems have been discovered with three or more planets participating in a sequence of mean-motion resonances. The unique and complex architectures of these so-called "resonant chains" motivate efforts to characterize their planets holistically. In this work, we perform a comprehensive exploration of the spin-axis dynamics of planets in resonant chains. Planetary spin states are closely linked with atmospheric dynamics and habitability and are thus especially relevant to resonant chains like TRAPPIST-1, which hosts several temperate planets. Considering a set of observed resonant chains, we calculate the equilibrium states of the planetary axial tilts ("obliquities"). We show that high obliquity states exist for $\sim60\%$ of planets in our sample, and many of these states can be stable in the presence of tidal dissipation. Using case studies of two observed systems (Kepler-223 and TOI-1136), we demonstrate how these high obliquity states could have been attained during the initial epoch of disk-driven orbital migration that established the resonant orbital architectures. We show that the TRAPPIST-1 planets most likely have zero obliquities, with the possible exception of planet d. Overall, our results highlight that both the orbital and spin states of resonant chains are valuable relics of the early stages of planet formation and evolution. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: Accepted for publication in ApJ

Tidal Dissipation Regimes Among the Short-Period Exoplanets

Authors: Emma Louden, Gregory Laughlin, Sarah Millholland

Abstract: The efficiency of tidal dissipation provides a zeroth-order link to a planet's physical properties. For super-Earth and sub-Neptune planets in the range $R_{\oplus}\lesssim R_p \lesssim 4 R_{\oplus}$, particularly efficient dissipation (i.e., low tidal quality factors) may signify terrestrial-like planets capable of maintaining rigid crustal features. Here we explore global constraints on planetar… ▽ More The efficiency of tidal dissipation provides a zeroth-order link to a planet's physical properties. For super-Earth and sub-Neptune planets in the range $R_{\oplus}\lesssim R_p \lesssim 4 R_{\oplus}$, particularly efficient dissipation (i.e., low tidal quality factors) may signify terrestrial-like planets capable of maintaining rigid crustal features. Here we explore global constraints on planetary tidal quality factors using a population of planets in multiple-planet systems whose orbital and physical properties indicate susceptibility to capture into secular spin-orbit resonances. Planets participating in secular spin-orbit resonance can maintain large axial tilts and significantly enhanced heating from obliquity tides. When obliquity tides are sufficiently strong, planets in low-order mean-motion resonances can experience resonant repulsion (period ratio increase). The observed distribution of period ratios among transiting planet pairs may thus depend non-trivially on the underlying planetary structures. We model the action of resonant repulsion and demonstrate that the observed distribution of period ratios near the 2:1 and 3:2 commensurabilties implies $Q$ values spanning from $Q\approx 10^1-10^7$ and peaking at $Q \approx 10^6$. This range includes the expected range in which super-Earth and sub-Neptune planets dissipate ($Q \approx 10^3 - 10^4$). This work serves as a proof of concept for a method of assessing the presence of two dissipation regimes, and we estimate the number of additional multi-transiting planetary systems needed to place any bimodality in the distribution on a strong statistical footing. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 9 pages, 6 figures, accepted for publication in ApJL

TOI-5126: A hot super-Neptune and warm Neptune pair discovered by $\textit{TESS}$ and $\textit{CHEOPS}$

Authors: Tyler R. Fairnington, Emma Nabbie, Chelsea X. Huang, George Zhou, Orion Foo, Sarah Millholland, Duncan Wright, Alexandre A. Belinski, Allyson Bieryla, David R. Ciardi, Karen A. Collins, Kevin I. Collins, Mark Everett, Steve B. Howell, Jack J. Lissauer, Michael B. Lund, Felipe Murgas, Enric Palle, Samuel N. Quinn, Howard M. Relles, Boris S. Safonov, Richard P. Schwarz, Nicholas J. Scott, Gregor Srdoc, George Ricker , et al. (11 additional authors not shown)

Abstract: We present the confirmation of a hot super-Neptune with an exterior Neptune companion orbiting a bright (V = 10.1 mag) F-dwarf identified by the $\textit{Transiting Exoplanet Survey Satellite}$ ($\textit{TESS}$). The two planets, observed in sectors 45, 46 and 48 of the $\textit{TESS}$ extended mission, are $4.74^{+0.16}_{-0.14}$ $R_{\oplus}$ and $3.86^{+0.17}_{-0.16}$ $R_{\oplus}$ with… ▽ More We present the confirmation of a hot super-Neptune with an exterior Neptune companion orbiting a bright (V = 10.1 mag) F-dwarf identified by the $\textit{Transiting Exoplanet Survey Satellite}$ ($\textit{TESS}$). The two planets, observed in sectors 45, 46 and 48 of the $\textit{TESS}$ extended mission, are $4.74^{+0.16}_{-0.14}$ $R_{\oplus}$ and $3.86^{+0.17}_{-0.16}$ $R_{\oplus}$ with $5.4588385^{+0.0000070}_{-0.0000072}$ d and $17.8999^{+0.0018}_{-0.0013}$ d orbital periods, respectively. We also obtained precise space based photometric follow-up of the system with ESAs $\textit{CHaracterising ExOplanets Satellite}$ ($\textit{CHEOPS}$) to constrain the radius and ephemeris of TOI-5126 b. TOI 5126 b is located in the "hot Neptune Desert" and is an ideal candidate for follow-up transmission spectroscopy due to its high predicted equilibrium temperature ($T_{eq} = 1442^{+46}_{-40}$ K) implying a cloud-free atmosphere. TOI-5126 c is a warm Neptune ($T_{eq}= 971^{+31}_{-27}$ K) also suitable for follow-up. Tentative transit timing variations (TTVs) have also been identified in analysis, suggesting the presence of at least one additional planet, however this signal may be caused by spot-crossing events, necessitating further precise photometric follow-up to confirm these signals. △ Less

Submitted 13 October, 2023; originally announced October 2023.

Comments: Accepted in MNRAS, 18 pages, 14 figures

Can Cold Jupiters Sculpt the Edge-of-the-Multis?

Authors: Nicole Sobski, Sarah C. Millholland

Abstract: Compact systems of multiple close-in super-Earths/sub-Neptunes ("compact multis") are a ubiquitous outcome of planet formation. It was recently discovered that the outer edges of compact multis are located at smaller orbital periods than expected from geometric and detection biases alone, suggesting some truncation or transition in the outer architectures. Here we test whether this "edge-of-the-mu… ▽ More Compact systems of multiple close-in super-Earths/sub-Neptunes ("compact multis") are a ubiquitous outcome of planet formation. It was recently discovered that the outer edges of compact multis are located at smaller orbital periods than expected from geometric and detection biases alone, suggesting some truncation or transition in the outer architectures. Here we test whether this "edge-of-the-multis" might be explained in any part by distant giant planets in the outer regions ($\gtrsim 1$ AU) of the systems. We investigate the dynamical stability of observed compact multis in the presence of hypothetical giant ($\gtrsim 0.5 \ M_{\mathrm{Jup}}$) perturbing planets. We identify what parameters would be required for hypothetical perturbing planets if they were responsible for dynamically sculpting the outer edges of compact multis. "Edge-sculpting" perturbers are generally in the range $P\sim100-500$ days for the average compact multi, with most between $P\sim200-300$ days. Given the relatively close separation, we explore the detectability of the hypothetical edge-sculpting perturbing planets, finding that they would be readily detectable in transit and radial velocity data. We compare to observational constraints and find it unlikely that dynamical sculpting from distant giant planets contributes significantly to the edge-of-the-multis. However, this conclusion could be strengthened in future work by a more thorough analysis of the detection yields of the perturbing planets. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: 15 pages, 5 figures, accepted for publication in ApJ

Two super-Earths at the edge of the habitable zone of the nearby M dwarf TOI-2095

Authors: F. Murgas, A. Castro-González, E. Pallé, F. J. Pozuelos, S. Millholland, O. Foo, J. Korth, E. Marfil, P. J. Amado, J. A. Caballero, J. L. Christiansen, D. R. Ciardi, K. A. Collins, M. Di Sora, A. Fukui, T. Gan, E. J. Gonzales, Th. Henning, E. Herrero, G. Isopi, J. M. Jenkins, J. Lillo-Box, N. Lodieu, R. Luque, F. Mallia , et al. (19 additional authors not shown)

Abstract: The main scientific goal of TESS is to find planets smaller than Neptune around stars that are bright enough to allow for further characterization studies. Given our current instrumentation and detection biases, M dwarfs are prime targets in the search for small planets that are in (or near) the habitable zone of their host star. In this work, we use photometric observations and CARMENES radial ve… ▽ More The main scientific goal of TESS is to find planets smaller than Neptune around stars that are bright enough to allow for further characterization studies. Given our current instrumentation and detection biases, M dwarfs are prime targets in the search for small planets that are in (or near) the habitable zone of their host star. In this work, we use photometric observations and CARMENES radial velocity measurements to validate a pair of transiting planet candidates found by TESS. The data were fitted simultaneously, using a Bayesian Markov chain Monte Carlo (MCMC) procedure and taking into account the stellar variability present in the photometric and spectroscopic time series. We confirm the planetary origin of the two transiting candidates orbiting around TOI-2095 (LSPM J1902+7525). The star is a nearby M dwarf ($d = 41.90 \pm 0.03$ pc, $T_{\rm eff} = 3759 \pm 87$ K, $V = 12.6$ mag), with a stellar mass and radius of $M_\star = 0.44 \pm 0.02 \; M_\odot$ and $R_\star = 0.44 \pm 0.02 \; R_\odot$, respectively. The planetary system is composed of two transiting planets: TOI-2095b, with an orbital period of $P_b = 17.66484 \pm (7\times 10^{-5})$ days, and TOI-2095c, with $P_c = 28.17232 \pm (14\times 10^{-5})$ days. Both planets have similar sizes with $R_b = 1.25 \pm 0.07 \; R_\oplus$ and $R_c = 1.33 \pm 0.08 \; R_\oplus$ for planet b and planet c, respectively. Although we did not detect the induced RV variations of any planet with significance, our CARMENES data allow us to set stringent upper limits on the masses of these objects. We find $M_b < 4.1 \; M_\oplus$ for the inner and $M_c < 7.4 \; M_\oplus$ for the outer planet (95% confidence level). These two planets present equilibrium temperatures in the range of 300-350 K and are close to the inner edge of the habitable zone of their star. △ Less

Submitted 2 August, 2023; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: Accepted for publication in A&A. 29 pages

Journal ref: A&A 677, A182 (2023)

Self-Consistent Spin, Tidal and Dynamical Equations of Motion in the REBOUNDx Framework

Authors: Tiger Lu, Hanno Rein, Daniel Tamayo, Sam Hadden, Rosemary Mardling, Sarah C. Millholland, Gregory Laughlin

Abstract: We have introduced self-consistent spin, tidal and dynamical equations of motion into REBOUNDx, a library of additional effects for the popular N-body integrator REBOUND. The equations of motion used are derived from the constant time lag approximation to the equilibrium tide model of tidal friction. These effects will allow the study of a variety of systems where the full dynamical picture cannot… ▽ More We have introduced self-consistent spin, tidal and dynamical equations of motion into REBOUNDx, a library of additional effects for the popular N-body integrator REBOUND. The equations of motion used are derived from the constant time lag approximation to the equilibrium tide model of tidal friction. These effects will allow the study of a variety of systems where the full dynamical picture cannot be encapsulated by point particle dynamics. We provide several test cases and benchmark the code's performance against analytic predictions. The open-source code is available in the most recent release of REBOUNDx. △ Less

Submitted 28 February, 2023; originally announced March 2023.

Comments: Accepted to ApJ. 15 pages, 6 figures

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)

Measured Spin-Orbit Alignment of Ultra-Short Period Super-Earth 55 Cancri e

Authors: Lily L. Zhao, Vedad Kunovac, John M. Brewer, Joe Llama, Sarah C. Millholland, Christina Hedges, Andrew E. Szymkowiak, Rachael M. Roettenbacher, Samuel H. C. Cabot, Sam A. Weiss, Debra A. Fischer

Abstract: A planet's orbital alignment places important constraints on how a planet formed and consequently evolved. The dominant formation pathway of ultra-short period planets ($P<1$ day) is particularly mysterious as such planets most likely formed further out, and it is not well understood what drove their migration inwards to their current positions. Measuring the orbital alignment is difficult for sma… ▽ More A planet's orbital alignment places important constraints on how a planet formed and consequently evolved. The dominant formation pathway of ultra-short period planets ($P<1$ day) is particularly mysterious as such planets most likely formed further out, and it is not well understood what drove their migration inwards to their current positions. Measuring the orbital alignment is difficult for smaller super-Earth/sub-Neptune planets, which give rise to smaller amplitude signals. Here we present radial velocities across two transits of 55 Cancri e, an ultra-short period Super-Earth, observed with the Extreme Precision Spectrograph (EXPRES). Using the classical Rossiter-McLaughlin (RM) method, we measure 55 Cnc e's sky-projected stellar spin-orbit alignment (i.e., the projected angle between the planet's orbital axis and its host star's spin axis) to be $λ=10\substack{+17\\ -20}^{\circ}$ with an unprojected angle of $ψ=23\substack{+14\\ -12}^{\circ}$. The best-fit RM model to the EXPRES data has a radial velocity semi-amplitude of just $0.41\substack{+0.09\\ -0.10} m s^{-1}$. The spin-orbit alignment of 55 Cnc e favors dynamically gentle migration theories for ultra-short period planets, namely tidal dissipation through low-eccentricity planet-planet interactions and/or planetary obliquity tides. △ Less

Submitted 9 December, 2022; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: 12 pages, 4 figures, published in Nature Astronomy

Inferred Properties of Planets in Mean-Motion Resonances are Biased by Measurement Noise

Authors: David Jensen, Sarah C. Millholland

Abstract: Planetary systems with mean-motion resonances (MMRs) hold special value in terms of their dynamical complexity and their capacity to constrain planet formation and migration histories. The key towards making these connections, however, is to have a reliable characterization of the resonant dynamics, especially the so-called "libration amplitude", which qualitatively measures how deep the system is… ▽ More Planetary systems with mean-motion resonances (MMRs) hold special value in terms of their dynamical complexity and their capacity to constrain planet formation and migration histories. The key towards making these connections, however, is to have a reliable characterization of the resonant dynamics, especially the so-called "libration amplitude", which qualitatively measures how deep the system is into the resonance. In this work, we identify an important complication with the interpretation of libration amplitude estimates from observational data of resonant systems. Specifically, we show that measurement noise causes inferences of the libration amplitude to be systematically biased to larger values, with noisier data yielding a larger bias. We demonstrated this through multiple approaches, including using dynamical fits of synthetic radial velocity data to explore how the the libration amplitude distribution inferred from the posterior parameter distribution varies with the degree of measurement noise. We find that even modest levels of noise still result in a slight bias. The origin of the bias stems from the topology of the resonant phase space and the fact that the available phase space volume increases non-uniformly with increasing libration amplitude. We highlight strategies for mitigating the bias through the usage of particular priors. Our results imply that many known resonant systems are likely deeper in resonance than previously appreciated. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Comments: Accepted in AJ, 9 pages and 7 figures

Edge-of-the-Multis: Evidence for a Transition in the Outer Architectures of Compact Multi-Planet Systems

Authors: Sarah C. Millholland, Matthias Y. He, Jon K. Zink

Abstract: Although the architectures of compact multiple-planet systems are well-characterized, there has been little examination of their "outer edges", or the locations of their outermost planets. Here we present evidence that the observed high-multiplicity Kepler systems truncate at smaller orbital periods than can be explained by geometric and detection biases alone. To show this, we considered the exis… ▽ More Although the architectures of compact multiple-planet systems are well-characterized, there has been little examination of their "outer edges", or the locations of their outermost planets. Here we present evidence that the observed high-multiplicity Kepler systems truncate at smaller orbital periods than can be explained by geometric and detection biases alone. To show this, we considered the existence of hypothetical planets orbiting beyond the observed transiting planets with properties dictated by the "peas-in-a-pod" patterns of intra-system radius and period ratio uniformity. We evaluated the detectability of these hypothetical planets using (1) a novel approach for estimating the mutual inclination dispersion of multi-transiting systems based on transit chord length ratios and (2) a model of transit probability and detection efficiency that accounts for the impacts of planet multiplicity on completeness. Under the assumption that the "peas-in-a-pod" patterns continue to larger orbital separations than observed, we find that $\gtrsim35\%$ of Kepler compact multis should possess additional detected planets beyond the known planets, constituting a $\sim7σ$ discrepancy with the lack of such detections. These results indicate that the outer ($\sim100-300$ days) regions of compact multis experience a truncation (i.e. an "edge-of-the-multis") or a significant breakdown of the "peas-in-a-pod" patterns, in the form of systematically smaller radii or larger period ratios. We outline future observations that can distinguish these possibilities, and we discuss implications for planet formation theories. △ Less

Submitted 20 July, 2022; originally announced July 2022.

Comments: Accepted to AJ

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

The Warm Neptune GJ 3470b has a Polar Orbit

Authors: Gudmundur Stefansson, Suvrath Mahadevan, Cristobal Petrovich, Joshua N. Winn, Shubham Kanodia, Sarah C. Millholland, Marissa Maney, Caleb I. Cañas, John Wisniewski, Paul Robertson, Joe P. Ninan, Eric B. Ford, Chad F. Bender, Cullen H. Blake, Heather Cegla, William D. Cochran, Scott A. Diddams, Jiayin Dong, Michael Endl, Connor Fredrick, Samuel Halverson, Fred Hearty, Leslie Hebb, Teruyuki Hirano, Andrea S. J. Lin , et al. (12 additional authors not shown)

Abstract: The warm Neptune GJ 3470b transits a nearby ($d=29$pc) bright slowly rotating M1.5-dwarf star. Using spectroscopic observations during two transits with the newly commissioned NEID spectrometer on the WIYN 3.5m Telescope at Kitt Peak Observatory, we model the classical Rossiter-Mclaughlin effect yielding a sky-projected obliquity of $λ=98_{-12}^{+15\:\circ}$ and a… ▽ More The warm Neptune GJ 3470b transits a nearby ($d=29$pc) bright slowly rotating M1.5-dwarf star. Using spectroscopic observations during two transits with the newly commissioned NEID spectrometer on the WIYN 3.5m Telescope at Kitt Peak Observatory, we model the classical Rossiter-Mclaughlin effect yielding a sky-projected obliquity of $λ=98_{-12}^{+15\:\circ}$ and a $v \sin i = 0.85_{-0.33}^{+0.27}$km/s. Leveraging information about the rotation period and size of the host star, our analysis yields a true obliquity of $ψ=95_{-8}^{+9\:\circ}$, revealing that GJ 3470b is on a polar orbit. Using radial velocities from HIRES, HARPS and the Habitable-zone Planet Finder, we show that the data are compatible with a long-term RV slope of $\dotγ = -0.0022 \pm 0.0011$m/s/day over a baseline of 12.9 years. If the RV slope is due to acceleration from another companion in the system, we show that such a companion is capable of explaining the polar and mildly eccentric orbit of GJ 3470b using two different secular excitation models. The existence of an outer companion can be further constrained with additional RV observations, Gaia astrometry, and future high-contrast imaging observations. Lastly, we show that tidal heating from GJ 3470b's mild eccentricity has most likely inflated the radius of GJ 3470b by a factor of $\sim$1.5-1.7, which could help account for its evaporating atmosphere. △ Less

Submitted 1 May, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 21 pages, 7 figures. Accepted in ApJL

Split Peas in a Pod: Intra-System Uniformity of Super-Earths and Sub-Neptunes

Authors: Sarah C. Millholland, Joshua N. Winn

Abstract: The planets within compact multi-planet systems tend to have similar sizes, masses, and orbital period ratios, like "peas in a pod". This pattern was detected when considering planets with radii between 1 and 4 $R_\oplus$. However, these same planets show a bimodal radius distribution, with few planets between 1.5 and 2 $R_{\oplus}$. The smaller "super-Earths" are consistent with being stripped ro… ▽ More The planets within compact multi-planet systems tend to have similar sizes, masses, and orbital period ratios, like "peas in a pod". This pattern was detected when considering planets with radii between 1 and 4 $R_\oplus$. However, these same planets show a bimodal radius distribution, with few planets between 1.5 and 2 $R_{\oplus}$. The smaller "super-Earths" are consistent with being stripped rocky cores, while the larger "sub-Neptunes" likely have gaseous H/He envelopes. Given these distinct structures, it is worthwhile to test for intra-system uniformity separately within each category of planets. Here, we find that the tendency for intra-system uniformity is twice as strong when considering planets within the same size category than it is when combining all planets together. The sub-Neptunes tend to be $1.7^{+0.6}_{-0.3}$ times larger than the super-Earths in the same system, corresponding to an envelope mass fraction of about 2.6% for a $5\,M_{\oplus}$ planet. For the sub-Neptunes, the low-metallicity stars are found to have planets with more equal sizes, with modest statistical significance ($p\sim 0.005$). There is also a modest $(\sim$2-$σ)$ tendency for wider-orbiting planets to be larger, even within the same size category. △ Less

Submitted 4 October, 2021; originally announced October 2021.

Comments: Accepted to ApJ Letters. 7 pages, 5 figures

Follow-up of non-transiting planets detected by Kepler. Confirmation of three hot-Jupiters and validation of three other planets

Authors: J. Lillo-Box, S. Millholland, G. Laughlin

Abstract: The direct detection of new extrasolar planets from high-precision photometry data is commonly based on the observation of the transit signal of the planet as it passes in front of its star. Close-in planets, however, leave additional imprints in the light curve even if they do not transit. These are the so-called phase curve variations that include ellipsoidal, reflection and beaming effects. In… ▽ More The direct detection of new extrasolar planets from high-precision photometry data is commonly based on the observation of the transit signal of the planet as it passes in front of its star. Close-in planets, however, leave additional imprints in the light curve even if they do not transit. These are the so-called phase curve variations that include ellipsoidal, reflection and beaming effects. In Millholland & Laughlin (2017), the authors scrutinized the Kepler database looking for these phase variations from non-transiting planets. They found 60 candidates whose signals were compatible with planetary companions. In this paper, we perform a ground-based follow-up of a sub-sample of these systems with the aim of confirming and characterizing these planets and thus validating the detection technique. We used the CAFE and HERMES instruments to monitor the radial velocity of ten non-transiting planet candidates along their orbits. We additionally used AstraLux to obtain high-resolution images of some of these candidates to discard blended binaries that contaminate the Kepler light curves by mimicking planetary signals. Among the ten systems, we confirm three new hot-Jupiters (KIC8121913 b, KIC10068024 b, and KIC5479689 b) with masses in the range 0.5-2 M$_{\rm Jup}$ and set mass constraints within the planetary regime for the other three candidates (KIC8026887b, KIC5878307 b, and KIC11362225 b), thus strongly suggestive of their planetary nature. For the first time, we validate the technique of detecting non-transiting planets via their phase curve variations. We present the new planetary systems and their properties. We find good agreement between the RV-derived masses and the photometric masses in all cases except KIC8121913 b, which shows a significantly lower mass derived from the ellipsoidal modulations than from beaming and radial velocity data. △ Less

Submitted 30 July, 2021; originally announced July 2021.

Comments: Accepted for publication in A&A. 17 pages, 8 figures, 14 tables

Journal ref: A&A 654, A9 (2021)

Evidence for a Non-Dichotomous Solution to the Kepler Dichotomy: Mutual Inclinations of Kepler Planetary Systems from Transit Duration Variations

Authors: Sarah C. Millholland, Matthias Y. He, Eric B. Ford, Darin Ragozzine, Daniel Fabrycky, Joshua N. Winn

Abstract: Early analyses of exoplanet statistics from the Kepler Mission revealed that a model population of multiple-planet systems with low mutual inclinations (${\sim1^{\circ}-2^{\circ}}$) adequately describes the multiple-transiting systems but underpredicts the number of single-transiting systems. This so-called "Kepler dichotomy" signals the existence of a sub-population of multi-planet systems posses… ▽ More Early analyses of exoplanet statistics from the Kepler Mission revealed that a model population of multiple-planet systems with low mutual inclinations (${\sim1^{\circ}-2^{\circ}}$) adequately describes the multiple-transiting systems but underpredicts the number of single-transiting systems. This so-called "Kepler dichotomy" signals the existence of a sub-population of multi-planet systems possessing larger mutual inclinations. However, the details of these inclinations remain uncertain. In this work, we derive constraints on the intrinsic mutual inclination distribution by statistically exploiting Transit Duration Variations (TDVs) of the Kepler planet population. When planetary orbits are mutually inclined, planet-planet interactions cause orbital precession, which can lead to detectable long-term changes in transit durations. These TDV signals are inclination-sensitive and have been detected for roughly two dozen Kepler planets. We compare the properties of the Kepler observed TDV detections to TDV detections of simulated planetary systems constructed from two population models with differing assumptions about the mutual inclination distribution. We find strong evidence for a continuous distribution of relatively low mutual inclinations that is well-characterized by a power law relationship between the median mutual inclination ($\tildeμ_{i,n}$) and the intrinsic multiplicity ($n$): $\tildeμ_{i,n} = \tildeμ_{i,5}(n/5)^α$, where $\tildeμ_{i,5} = 1.10^{+0.15}_{-0.11}$ and $α= -1.73^{+0.09}_{-0.08}$. These results suggest that late-stage planet assembly and possibly stellar oblateness are the dominant physical origins for the excitation of Kepler planet mutual inclinations. △ Less

Submitted 29 June, 2021; originally announced June 2021.

Comments: Accepted to AJ, 19 pages, 10 figures (excluding appendix)

Formation of Ultra-Short-Period Planets by Obliquity-Driven Tidal Runaway

Authors: Sarah Millholland, Christopher Spalding

Abstract: Small, rocky planets have been found orbiting in extreme proximity to their host stars, sometimes down to only $\sim 2$ stellar radii. These ultra-short-period planets (USPs) likely did not form in their present-day orbits, but rather migrated from larger initial separations. While tides are the probable cause of this migration, the tidal source has remained uncertain. Here we introduce planetary… ▽ More Small, rocky planets have been found orbiting in extreme proximity to their host stars, sometimes down to only $\sim 2$ stellar radii. These ultra-short-period planets (USPs) likely did not form in their present-day orbits, but rather migrated from larger initial separations. While tides are the probable cause of this migration, the tidal source has remained uncertain. Here we introduce planetary obliquity tides as a natural pathway for the production of USPs within close-in multi-planet systems. The crucial idea is that tidal dissipation generally forces planetary spin vectors to equilibrium configurations called "Cassini states", in which the planetary obliquities (axial tilts) are non-zero. In these cases, sustained tidal dissipation and inward orbital migration are inevitable. Migration then increases the obliquity and strengthens the tides, creating a positive feedback loop. Thus, if a planet's initial semi-major axis is small enough ($a \lesssim 0.05$ AU), it can experience runaway orbital decay, which is stalled at ultra-short orbital periods when the forced obliquity reaches very high values ($\sim 85^{\circ}$) and becomes unstable. We use secular dynamics to outline the parameter space in which the innermost member of a prototypical Kepler multiple-planet system can become a USP. We find that these conditions are consistent with many observed features of USPs, such as period ratios, mutual inclinations, and occurrence rate trends with stellar type. Future detections of stellar obliquities and close-in companions, together with theoretical explorations of the potential for chaotic obliquity dynamics, can help constrain the prevalence of this mechanism. △ Less

Submitted 28 October, 2020; originally announced October 2020.

Comments: 24 pages, 11 figures. Accepted to ApJ

Integrating Machine Learning for Planetary Science: Perspectives for the Next Decade

Authors: Abigail R. Azari, John B. Biersteker, Ryan M. Dewey, Gary Doran, Emily J. Forsberg, Camilla D. K. Harris, Hannah R. Kerner, Katherine A. Skinner, Andy W. Smith, Rashied Amini, Saverio Cambioni, Victoria Da Poian, Tadhg M. Garton, Michael D. Himes, Sarah Millholland, Suranga Ruhunusiri

Abstract: Machine learning (ML) methods can expand our ability to construct, and draw insight from large datasets. Despite the increasing volume of planetary observations, our field has seen few applications of ML in comparison to other sciences. To support these methods, we propose ten recommendations for bolstering a data-rich future in planetary science. Machine learning (ML) methods can expand our ability to construct, and draw insight from large datasets. Despite the increasing volume of planetary observations, our field has seen few applications of ML in comparison to other sciences. To support these methods, we propose ten recommendations for bolstering a data-rich future in planetary science. △ Less

Submitted 29 July, 2020; originally announced July 2020.

Comments: 10 pages (expanded citations compared to 8 page submitted version for decadal survey), 3 figures, white paper submitted to the Planetary Science and Astrobiology Decadal Survey 2023-2032

Stellar Oblateness versus Distant Giants in Exciting Kepler Planet Mutual Inclinations

Authors: Christopher Spalding, Sarah Millholland

Abstract: An overabundance of single-transiting Kepler planets suggests the existence of a sub-population of intrinsically multi-planet systems possessing large mutual inclinations. However, the origin of these mutual inclinations remains unknown. Recent work has demonstrated that mutual inclinations can be excited soon after protoplanetary disk-dispersal due to the oblateness of the rapidly-rotating host s… ▽ More An overabundance of single-transiting Kepler planets suggests the existence of a sub-population of intrinsically multi-planet systems possessing large mutual inclinations. However, the origin of these mutual inclinations remains unknown. Recent work has demonstrated that mutual inclinations can be excited soon after protoplanetary disk-dispersal due to the oblateness of the rapidly-rotating host star, provided the star is tilted. Alternatively, distant giant planets, which are common in systems of close-in Kepler planets, could drive up mutual inclinations. The relative importance of each of these mechanisms has not been investigated. Here, we show that the influence of the stellar oblateness typically exceeds that of an exterior giant soon after planet formation. However, the magnitude of the resulting mutual inclinations depends critically upon the timescale over which the natal disk disperses. Specifically, we find that if the disk vanishes over a timescale shorter than 10^(3-4) years, comparable to the viscous timescale of the inner ~0.2AU, the inner planets impulsively acquire misalignments that scale with the stellar obliquity. In contrast, if the disk disperses slowly, the inner planets remain coplanar. They first align with the stellar equator but subsequently realign with the distant giant's plane as the star spins down. Our findings are consistent with recent observations that giants tend to be aligned with close-in multis but misaligned with singles. Stellar obliquity measurements offer a promising test of our proposed framework. △ Less

Submitted 14 July, 2020; originally announced July 2020.

Comments: 22 pages, 11 figures. Accepted for publication in The Astronomical Journal

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/

Obliquity Constraints on an Extrasolar Planetary-Mass Companion

Authors: Marta L. Bryan, Eugene Chiang, Brendan P. Bowler, Caroline V. Morley, Sarah Millholland, Sarah Blunt, Katelyn B. Ashok, Eric Nielsen, Henry Ngo, Dimitri Mawet, Heather A. Knutson

Abstract: We place the first constraints on the obliquity of a planetary-mass companion (PMC) outside of the Solar System. Our target is the directly imaged system 2MASS J01225093-2439505 (2M0122), which consists of a 120 Myr 0.4 M_sun star hosting a 12-27 M_J companion at 50 AU. We constrain all three of the system's angular momentum vectors: how the companion spin axis, the stellar spin axis, and the orbi… ▽ More We place the first constraints on the obliquity of a planetary-mass companion (PMC) outside of the Solar System. Our target is the directly imaged system 2MASS J01225093-2439505 (2M0122), which consists of a 120 Myr 0.4 M_sun star hosting a 12-27 M_J companion at 50 AU. We constrain all three of the system's angular momentum vectors: how the companion spin axis, the stellar spin axis, and the orbit normal are inclined relative to our line of sight. To accomplish this, we measure projected rotation rates (vsini) for both the star and the companion using new near-infrared high-resolution spectra with NIRSPEC at Keck Observatory. We combine these with a new stellar photometric rotation period from TESS and a published companion rotation period from HST to obtain spin axis inclinations for both objects. We also fitted multiple epochs of astrometry, including a new observation with NIRC2/Keck, to measure 2M0122b's orbital inclination. The three line-of-sight inclinations place limits on the true de-projected companion obliquity and stellar obliquity. We find that while the stellar obliquity marginally prefers alignment, the companion obliquity tentatively favors misalignment. We evaluate possible origin scenarios. While collisions, secular spin-orbit resonances, and Kozai-Lidov oscillations are unlikely, formation by gravitational instability in a gravito-turbulent disk - the scenario favored for brown dwarf companions to stars - appears promising. △ Less

Submitted 27 May, 2020; v1 submitted 25 February, 2020; originally announced February 2020.

Comments: published in AJ

TOI 564 b and TOI 905 b: Grazing and Fully Transiting Hot Jupiters Discovered by TESS

Authors: Allen B. Davis, Songhu Wang, Matias Jones, Jason D. Eastman, Maximilian N. Günther, Keivan G. Stassun, Brett C. Addison, Karen A. Collins, Samuel N. Quinn, David W. Latham, Trifon Trifonov, Sahar Shahaf, Tsevi Mazeh, Stephen R. Kane, Xian-Yu Wang, Thiam-Guan Tan, Andrei Tokovinin, Carl Ziegler, René Tronsgaard, Sarah Millholland, Bryndis Cruz, Perry Berlind, Michael L. Calkins, Gilbert A. Esquerdo, Kevin I. Collins , et al. (24 additional authors not shown)

Abstract: We report the discovery and confirmation of two new hot Jupiters discovered by the Transiting Exoplanet Survey Satellite (TESS): TOI 564 b and TOI 905 b. The transits of these two planets were initially observed by TESS with orbital periods of 1.651 d and 3.739 d, respectively. We conducted follow-up observations of each system from the ground, including photometry in multiple filters, speckle int… ▽ More We report the discovery and confirmation of two new hot Jupiters discovered by the Transiting Exoplanet Survey Satellite (TESS): TOI 564 b and TOI 905 b. The transits of these two planets were initially observed by TESS with orbital periods of 1.651 d and 3.739 d, respectively. We conducted follow-up observations of each system from the ground, including photometry in multiple filters, speckle interferometry, and radial velocity measurements. For TOI 564 b, our global fitting revealed a classical hot Jupiter with a mass of $1.463^{+0.10}_{-0.096}\ M_J$ and a radius of $1.02^{+0.71}_{-0.29}\ R_J$. TOI 905 b is a classical hot Jupiter as well, with a mass of $0.667^{+0.042}_{-0.041}\ M_J$ and radius of $1.171^{+0.053}_{-0.051}\ R_J$. Both planets orbit Sun-like, moderately bright, mid-G dwarf stars with V ~ 11. While TOI 905 b fully transits its star, we found that TOI 564 b has a very high transit impact parameter of $0.994^{+0.083}_{-0.049}$, making it one of only ~20 known systems to exhibit a grazing transit and one of the brightest host stars among them. TOI 564 b is therefore one of the most attractive systems to search for additional non-transiting, smaller planets by exploiting the sensitivity of grazing transits to small changes in inclination and transit duration over the time scale of several years. △ Less

Submitted 20 December, 2019; originally announced December 2019.

Comments: 21 pages and 10 figures. Submitted to AJ

Tidally-Induced Radius Inflation of Sub-Neptunes

Authors: Sarah Millholland

Abstract: Recent work suggests that many short-period super-Earth and sub-Neptune planets may have significant spin axis tilts ("obliquities"). When planets are locked in high-obliquity states, the tidal dissipation rate may increase by several orders of magnitude. This intensified heat deposition within the planets' interiors should generate significant structural consequences, including atmospheric inflat… ▽ More Recent work suggests that many short-period super-Earth and sub-Neptune planets may have significant spin axis tilts ("obliquities"). When planets are locked in high-obliquity states, the tidal dissipation rate may increase by several orders of magnitude. This intensified heat deposition within the planets' interiors should generate significant structural consequences, including atmospheric inflation leading to larger transit radii. Using up-to-date radius estimates from Gaia Data Release 2, we show evidence for $\sim50\%$ larger average radii of planets wide of first-order mean-motion resonances, a population of planets with a theorized frequent occurrence of high obliquities. We investigate whether this radius trend could be a signature of obliquity tides. Using an adaptation of the Modules for Experiments in Stellar Astrophysics (MESA) stellar evolution toolkit, we model the atmospheric evolution of sub-Neptune-mass planets in response to additional internal heat from obliquity tides. The degree of radius inflation predicted by the models is $\sim10\%-100\%$ for tidal luminosities $\gtrsim 10^{-5}$ of the incident stellar power; this degree of inflation is broadly consistent with the observations and can approximately be described by power law relationships. We present a few case studies of very low density "super-puff" planets -- Kepler-79 d, Kepler-31 c, and Kepler-27 b -- and show that they are strong candidates for potentially having undergone tidally-induced radius inflation. We also discuss how the discrepancy between the two populations of planets with masses derived from radial velocities and transit timing variations is connected to the radius distribution features we have identified. Altogether, the calculations in this work confirm that tidal dissipation has non-negligible consequences for the structural properties of short-period sub-Neptunes. △ Less

Submitted 17 October, 2019; v1 submitted 15 October, 2019; originally announced October 2019.

Comments: 17 pages, 11 figures, accepted to ApJ

Signatures of Obliquity in Thermal Phase Curves of Hot Jupiters

Authors: Arthur D. Adams, Sarah Millholland, Gregory P. Laughlin

Abstract: Recent work suggests that many short-period extrasolar planets may have spin obliquities that are significantly tilted with respect to their orbital planes. These large obliquities are a natural outcome of ``secular spin-orbit resonance'', a configuration in which the planetary spin precession frequency matches the frequency of orbit nodal regression, or a Fourier component thereof. While exoplane… ▽ More Recent work suggests that many short-period extrasolar planets may have spin obliquities that are significantly tilted with respect to their orbital planes. These large obliquities are a natural outcome of ``secular spin-orbit resonance'', a configuration in which the planetary spin precession frequency matches the frequency of orbit nodal regression, or a Fourier component thereof. While exoplanet spin obliquities have not yet been measured directly, they may be detectable indirectly through their signatures in various observations, such as photometric measurements across the full phase of a planet's orbit. In this work, we employ a thermal radiative model to explore how large polar tilts affect full-phase light curves, and we discuss the range of unique signatures that are expected to result. We show that the well-studied short-period planets HD 149026 b, WASP-12 b, and CoRoT-2 b all exhibit phase curve features that may arise from being in high-obliquity states. We also constrain the parameters and assess the detectability of hypothetical perturbing planets that could maintain the planets in these states. Among the three planets considered, CoRoT-2 b has the tightest constraints on its proposed obliquity ($45.8^{\circ} \pm 1.4^{\circ}$) and axial orientation. For HD 149026 b, we find no significant evidence for a non-zero obliquity, and the phase curve of WASP-12 b is too complicated by strong tidal distortions for a conclusive assessment. △ Less

Submitted 18 June, 2019; originally announced June 2019.

Comments: 21 pages, 7 figures, 6 tables, accepted for publication in AJ

Excitation of Planetary Obliquities Through Planet-Disk Interactions

Authors: Sarah Millholland, Konstantin Batygin

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

Submitted 15 April, 2019; originally announced April 2019.

Comments: 12 pages, 8 figures, accepted to ApJ

Obliquity-Driven Sculpting of Exoplanetary Systems

Authors: Sarah Millholland, Gregory Laughlin

Abstract: NASA's Kepler mission revealed that $\sim 30\%$ of Solar-type stars harbor planets with sizes between that of Earth and Neptune on nearly circular and co-planar orbits with periods less than 100 days. Such short-period compact systems are rarely found with planet pairs in mean-motion resonances (MMRs) -- configurations in which the planetary orbital periods exhibit a simple integer ratio -- but th… ▽ More NASA's Kepler mission revealed that $\sim 30\%$ of Solar-type stars harbor planets with sizes between that of Earth and Neptune on nearly circular and co-planar orbits with periods less than 100 days. Such short-period compact systems are rarely found with planet pairs in mean-motion resonances (MMRs) -- configurations in which the planetary orbital periods exhibit a simple integer ratio -- but there is a significant overabundance of planet pairs lying just wide of the first-order resonances. Previous work suggests that tides raised on the planets by the host star may be responsible for forcing systems into these configurations by draining orbital energy to heat. Such tides, however, are insufficient unless there exists a substantial and as-yet unidentified source of extra dissipation. Here we show that this cryptic heat source may be linked to "obliquity tides" generated when a large axial tilt (obliquity) is maintained by secular resonance-driven spin-orbit coupling. We present evidence that typical compact, nearly-coplanar systems frequently experience this mechanism, and we highlight additional features in the planetary orbital period and radius distributions that may be its signatures. Extrasolar planets that maintain large obliquities will exhibit infrared light curve features that are detectable with forthcoming space missions. The observed period ratio distribution can be explained if typical tidal quality factors for super-Earths and sub-Neptunes are similar to those of Uranus and Neptune. △ Less

Submitted 4 March, 2019; originally announced March 2019.

Comments: Published in Nature Astronomy on March 4, 2019. 38 pages and 14 figures including the methods and supplementary information

Journal ref: Nature Astronomy (2019), https://www.nature.com/articles/s41550-019-0701-7

Obliquity Tides May Drive WASP-12b's Rapid Orbital Decay

Authors: Sarah Millholland, Gregory Laughlin

Abstract: Recent analyses have revealed a mystery. The orbital period of the highly inflated hot Jupiter, WASP-12b, is decreasing rapidly. The rate of inspiral, however, is too fast to be explained by either eccentricity tides or equilibrium stellar tides. While dynamical stellar tides are possible, they require a subgiant structure for the star, whereas stellar models point toward a main sequence host. Her… ▽ More Recent analyses have revealed a mystery. The orbital period of the highly inflated hot Jupiter, WASP-12b, is decreasing rapidly. The rate of inspiral, however, is too fast to be explained by either eccentricity tides or equilibrium stellar tides. While dynamical stellar tides are possible, they require a subgiant structure for the star, whereas stellar models point toward a main sequence host. Here, we show that these hitherto irreconcilable observations might be explained by planetary obliquity tides if planet b's spin vector is trapped in a high-obliquity state maintained by a secular spin-orbit resonance with an unseen exterior perturbing planet. We derive constraints on the obliquity ($ε\gtrsim50^{\circ}$), reduced tidal quality factor ($Q^{\prime}\sim10^{6}-10^{7}$), and perturbing planet parameters ($M_{2}\sim10-20M_{\oplus}$, $a_2\lesssim0.04\,{\rm AU}$) required to generate the observed orbital decay. Direct N-body simulations that include tidal and spin dynamics reinforce the plausibility of the scenario. Furthermore, we show that the resonance could have been captured when planet b's obliquity was small, making the proposed sequence of events easy to explain. The hypothetical perturbing planet is within the limits of current radial velocity constraints on the system yet is also detectable. If it exists, it could provide evidence in favor of the in situ formation hypothesis for hot Jupiters. △ Less

Submitted 10 December, 2018; v1 submitted 4 December, 2018; originally announced December 2018.

Comments: published in ApJL, small error in Figure 4 corrected in this version

HD 202772A B: A Transiting Hot Jupiter Around A Bright, Mildly Evolved Star In A Visual Binary Discovered By Tess

Authors: Songhu Wang, Matias Jones, Avi Shporer, Benjamin J. Fulton, Leonardo A. Paredes, Trifon Trifonov, Diana Kossakowski, Jason Eastman, Maximilian N. Gunther, Chelsea X. Huang, Sarah Millholland, Darryl Seligman, Debra Fischer, Rafael Brahm, Xian-Yu Wang, Bryndis Cruz, Hodari-Sadiki James, Brett Addison, Todd Henry, En-Si Liang, Allen B. Davis, Rene Tronsgaard, Keduse Worku, John Brewer, Martin Kurster , et al. (30 additional authors not shown)

Abstract: We report the first confirmation of a hot Jupiter discovered by the Transiting Exoplanet Survey Satellite (TESS) mission: HD 202772A b. The transit signal was detected in the data from TESS Sector 1, and was confirmed to be of planetary origin through radial-velocity measurements. HD 202772A b is orbiting a mildly evolved star with a period of 3.3 days. With an apparent magnitude of V = 8.3, the s… ▽ More We report the first confirmation of a hot Jupiter discovered by the Transiting Exoplanet Survey Satellite (TESS) mission: HD 202772A b. The transit signal was detected in the data from TESS Sector 1, and was confirmed to be of planetary origin through radial-velocity measurements. HD 202772A b is orbiting a mildly evolved star with a period of 3.3 days. With an apparent magnitude of V = 8.3, the star is among the brightest known to host a hot Jupiter. Based on the 27days of TESS photometry, and radial velocity data from the CHIRON and HARPS spectrographs, the planet has a mass of 1.008+/-0.074 M_J and radius of 1.562+/-0.053 R_J , making it an inflated gas giant. HD 202772A b is a rare example of a transiting hot Jupiter around a quickly evolving star. It is also one of the most strongly irradiated hot Jupiters currently known. △ Less

Submitted 4 October, 2018; originally announced October 2018.

Comments: Submitted to AAS Journal

Discovery and Dynamical Analysis of an Extreme Trans-Neptunian Object with a High Orbital Inclination

Authors: J. C. Becker, T. Khain, S. J. Hamilton, F. Adams, D. W. Gerdes, L. Zullo, K. Franson, S. Millholland, G. M. Bernstein, M. Sako, P. Bernardinelli, K. Napier, L. Markwardt, Hsing Wen Lin, W. Wester, F. B. Abdalla, S. Allam, J. Annis, S. Avila, E. Bertin, D. Brooks, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, C. E. Cunha , et al. (40 additional authors not shown)

Abstract: We report the discovery and dynamical analysis of 2015 BP$_{519}$, an extreme Trans-Neptunian Object detected detected by the Dark Energy Survey at a heliocentric distance of 55 AU and absolute magnitude Hr= 4.3. The current orbit, determined from a 1110-day observational arc, has semi-major axis $a\approx$ 450 AU, eccentricity $e\approx$ 0.92 and inclination $i\approx$ 54 degrees. With these orbi… ▽ More We report the discovery and dynamical analysis of 2015 BP$_{519}$, an extreme Trans-Neptunian Object detected detected by the Dark Energy Survey at a heliocentric distance of 55 AU and absolute magnitude Hr= 4.3. The current orbit, determined from a 1110-day observational arc, has semi-major axis $a\approx$ 450 AU, eccentricity $e\approx$ 0.92 and inclination $i\approx$ 54 degrees. With these orbital elements, 2015 BP$_{519}$ is the most extreme TNO discovered to date, as quantified by the reduced Kozai action, which is is a conserved quantity at fixed semi-major axis $a$ for axisymmetric perturbations. We discuss the orbital stability and evolution of this object in the context of the known Solar System, and find that 2015 BP$_{519}$ displays rich dynamical behavior, including rapid diffusion in semi-major axis and more constrained variations in eccentricity and inclination. We also consider the long term orbital stability and evolutionary behavior within the context of the Planet Nine Hypothesis, and find that BP$_{519}$ adds to the circumstantial evidence for the existence of this proposed new member of the Solar System, as it would represent the first member of the population of high-i, $\varpi$-shepherded TNOs. △ Less

Submitted 5 August, 2018; v1 submitted 14 May, 2018; originally announced May 2018.

Comments: published in AJ; observations available on MPC

Journal ref: 2018. The Astronomical Journal, Volume 156, Number 2

New Constraints on Gliese 876 -- Exemplar of Mean-Motion Resonance

Authors: Sarah Millholland, Gregory Laughlin, Johanna Teske, R. Paul Butler, Jennifer Burt, Bradford Holden, Steven Vogt, Jeffrey Crane, Stephen Shectman, Ian Thompson

Abstract: Gliese 876 harbors one of the most dynamically rich and well-studied exoplanetary systems. The nearby M4V dwarf hosts four known planets, the outer three of which are trapped in a Laplace mean-motion resonance. A thorough characterization of the complex resonant perturbations exhibited by the orbiting planets, and the chaotic dynamics therein, is key to a complete picture of the system's formation… ▽ More Gliese 876 harbors one of the most dynamically rich and well-studied exoplanetary systems. The nearby M4V dwarf hosts four known planets, the outer three of which are trapped in a Laplace mean-motion resonance. A thorough characterization of the complex resonant perturbations exhibited by the orbiting planets, and the chaotic dynamics therein, is key to a complete picture of the system's formation and evolutionary history. Here we present a reanalysis of the system using six years of new radial velocity (RV) data from four instruments. This new data augments and more than doubles the size of the decades-long collection of existing velocity measurements. We provide updated estimates of the system parameters by employing a computationally efficient Wisdom-Holman N-body symplectic integrator, coupled with a Gaussian Process (GP) regression model to account for correlated stellar noise. Experiments with synthetic RV data show that the dynamical characterization of the system can differ depending on whether a white noise or correlated noise model is adopted. Despite there being a region of stability for an additional planet in the resonant chain, we find no evidence for one. Our new parameter estimates place the system even deeper into resonance than previously thought and suggest that the system might be in a low energy, quasi-regular double apsidal corotation resonance. This result and others will be used in a subsequent study on the primordial migration processes responsible for the formation of the resonant chain. △ Less

Submitted 11 August, 2022; v1 submitted 23 January, 2018; originally announced January 2018.

Comments: Added the stellar mass in the caption of Table 4 in this version. 18 pages, 7 figures, published in AJ. Posterior samples available at https://github.com/smillholland/GJ876

Kepler Multi-Planet Systems Exhibit Unexpected Intra-system Uniformity in Mass and Radius

Authors: Sarah Millholland, Songhu Wang, Gregory Laughlin

Abstract: The widespread prevalence of close-in, nearly coplanar super-Earth- and sub-Neptune-sized planets in multiple-planet systems was one of the most surprising results from the Kepler mission. By studying a uniform sample of Kepler "multis" with mass measurements from transit timing variations (TTVs), we show that a given planetary system tends to harbor a characteristic type of planet. That is, plane… ▽ More The widespread prevalence of close-in, nearly coplanar super-Earth- and sub-Neptune-sized planets in multiple-planet systems was one of the most surprising results from the Kepler mission. By studying a uniform sample of Kepler "multis" with mass measurements from transit timing variations (TTVs), we show that a given planetary system tends to harbor a characteristic type of planet. That is, planets in a system have both masses and radii that are far more similar than if the system were assembled randomly from planets in the population. This finding has two important ramifications. First, the large intrinsic compositional scatter in the planet mass-radius relation is dominated by system-to-system variance rather than intra-system variance. Second, if provided enough properties of the star and primordial protoplanetary disk, there may be a substantial degree of predictability in the outcome of the planet formation process. We show that stellar mass and metallicity account for of order $20\%$ of the variation in outcomes; the remainder is as-yet unknown. △ Less

Submitted 30 October, 2017; originally announced October 2017.

Comments: 7 pages, 4 figures. Accepted to ApJL. Code available at https://github.com/smillholland/Kepler_Multis_Uniformity

Supervised Learning Detection of Sixty Non-Transiting Hot Jupiter Candidates

Authors: Sarah Millholland, Gregory Laughlin

Abstract: The optical, full-phase photometric variations of a short-period planet provide a unique view of the planet's atmospheric composition and dynamics. The number of planets with optical phase curve detections, however, is currently too small to study them as an aggregate population, motivating an extension of the search to non-transiting planets. Here we present an algorithm for the detection of non-… ▽ More The optical, full-phase photometric variations of a short-period planet provide a unique view of the planet's atmospheric composition and dynamics. The number of planets with optical phase curve detections, however, is currently too small to study them as an aggregate population, motivating an extension of the search to non-transiting planets. Here we present an algorithm for the detection of non-transiting, short-period giant planets in the Kepler field. The procedure uses the phase curves themselves as evidence for the planets' existence. We employ a supervised learning algorithm to recognize the salient time-dependent properties of synthetic phase curves; we then search for detections of signals that match these properties. After demonstrating the algorithm's capabilities, we classify 142,630 FGK Kepler stars without confirmed planets or KOIs and, for each one, assign a probability of a phase curve of a non-transiting planet being present. We identify 60 high-probability non-transiting hot Jupiter candidates. We also derive constraints on the candidates' albedos and offsets of the phase curve maxima. These targets are strong candidates for follow-up radial velocity confirmation and characterization. Once confirmed, the atmospheric information content in the phase curves may be studied in yet greater detail. △ Less

Submitted 20 June, 2017; originally announced June 2017.

Comments: 23 pages, 20 figures. Accepted in AJ. Online repository of candidates available at https://smillholland.github.io/Non-transiting_HJs/

Constraints on Planet Nine's Orbit and Sky Position within a Framework of Mean Motion Resonances

Authors: Sarah Millholland, Gregory Laughlin

Abstract: A number of authors have proposed that the statistically significant orbital alignment of the most distant Kuiper Belt Objects (KBOs) is evidence of an as-yet undetected planet in the outer solar system, now referred to colloquially a "Planet Nine". Dynamical simulations by Batygin & Brown (2016) have provided constraints on the range of the planet's possible orbits and sky locations. We extend th… ▽ More A number of authors have proposed that the statistically significant orbital alignment of the most distant Kuiper Belt Objects (KBOs) is evidence of an as-yet undetected planet in the outer solar system, now referred to colloquially a "Planet Nine". Dynamical simulations by Batygin & Brown (2016) have provided constraints on the range of the planet's possible orbits and sky locations. We extend these investigations by exploring the suggestion of Malhotra et al. (2016) that Planet Nine is in small integer ratio mean-motion resonances (MMRs) with several of the most distant KBOs. We show that the observed KBO semi-major axes present a set of commensurabilities with an unseen planet at $\sim 654 \ \mathrm{AU}$ ($P \sim 16,725 \ \mathrm{yr}$) that has a greater than $98\%$ chance of stemming from a sequence of MMRs rather than from a random distribution. We describe and implement a Monte-Carlo optimization scheme that drives billion-year dynamical integrations of the outer solar system to pinpoint the orbital properties of perturbers that are capable of maintaining the KBOs' apsidal alignment. This optimization exercise suggests that the unseen planet is most consistently represented with mass, $m \sim 6-12 M_{\oplus}$, semi-major axis, $a \sim 654 \ \mathrm{AU}$, eccentricity, $e\sim0.45$, inclination, $i \sim 30^{\circ}$, argument of periastron, $ω\sim 150^{\circ}$, longitude of ascending node, $Ω\sim 50^{\circ}$, and mean anomaly, $M \sim 180^{\circ}$. A range of sky locations relative to this fiducial ephemeris are possible. We find that the region $30^{\circ}\, \lesssim \mathrm{RA} \lesssim 50^{\circ}$, $-20^{\circ}\, \lesssim \mathrm{Dec} \lesssim 20^{\circ}$ is promising. △ Less

Submitted 22 December, 2016; originally announced December 2016.

Journal ref: AJ, 153, 91 (2017)

Evolution of Flare Ribbons, Electric Currents and Quasi-separatrix Layers During an X-class Flare

Authors: M. Janvier, A. Savcheva, E. Pariat, S. Tassev, S. Millholland, V. Bommier, P. McCauley, S. McKillop, F. Dougan

Abstract: The standard model for eruptive flares has in the past few years been extended to 3D. It predicts typical J-shaped photospheric footprints of the coronal current layer, forming at similar locations as the Quasi-Separatrix Layers (QSLs). Such a morphology is also found for flare ribbons observed in the EUV band, as well as in non-linear force-free field (NLFFF) magnetic field extrapolations and mod… ▽ More The standard model for eruptive flares has in the past few years been extended to 3D. It predicts typical J-shaped photospheric footprints of the coronal current layer, forming at similar locations as the Quasi-Separatrix Layers (QSLs). Such a morphology is also found for flare ribbons observed in the EUV band, as well as in non-linear force-free field (NLFFF) magnetic field extrapolations and models. We study the evolution of the photospheric traces of the current density and flare ribbons, both obtained with the SDO instruments. We investigate the photospheric current evolution during the 6 September 2011 X-class flare (SOL2011-09-06T22:20) from observational data of the magnetic field obtained with HMI. This evolution is compared with that of the flare ribbons observed in the EUV filters of the AIA. We also compare the observed electric current density and the flare ribbon morphology with that of the QSLs computed from the flux rope insertion method/NLFFF model. The NLFFF model shows the presence of a fan-spine configuration of overlying field lines, due to the presence of a parasitic polarity, embedding an elongated flux rope that appears in the observations as two parts of a filament. The QSLs, evolved via a magnetofrictional method, also show similar morphology and evolution as both the current ribbons and the EUV flare ribbons obtained at several times during the flare. For the first time, we propose a combined analysis of the photospheric traces of an eruptive flare, in a complex topology, with direct measurements of electric currents and QSLs from observational data and a magnetic field model. The results, obtained by two different and independent approaches, 1) confirm previous results of current increase during the impulsive phase of the flare, 2) show how NLFFF models can capture the essential physical signatures of flares even in a complex magnetic field topology. △ Less

Submitted 25 April, 2016; originally announced April 2016.

Comments: 17 pages, 14 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 591, A141 (2016)

On the Detection of Non-Transiting Hot Jupiters in Multiple-Planet Systems

Authors: Sarah Millholland, Songhu Wang, Gregory Laughlin

Abstract: We outline a photometric method for detecting the presence of a non-transiting short-period giant planet in a planetary system harboring one or more longer period transiting planets. Within a prospective system of the type that we consider, a hot Jupiter on an interior orbit inclined to the line-of-sight signals its presence through approximately sinusoidal full-phase photometric variations in the… ▽ More We outline a photometric method for detecting the presence of a non-transiting short-period giant planet in a planetary system harboring one or more longer period transiting planets. Within a prospective system of the type that we consider, a hot Jupiter on an interior orbit inclined to the line-of-sight signals its presence through approximately sinusoidal full-phase photometric variations in the stellar light curve, correlated with astrometrically induced transit timing variations for exterior transiting planets. Systems containing a hot Jupiter along with a low-mass outer planet or planets on inclined orbits are a predicted hallmark of in situ accretion for hot Jupiters, and their presence can thus be used to test planetary formation theories. We outline the prospects for detecting non-transiting hot Jupiters using photometric data from typical \textit{Kepler} objects of interest (KOIs). As a demonstration of the technique, we perform a brief assessment of \textit{Kepler} candidates and identify a potential non-transiting hot Jupiter in the KOI-1822 system. Candidate non-transiting hot Jupiters can be readily confirmed with a small number of Doppler velocity observations, even for stars with $V\gtrsim14$. △ Less

Submitted 17 May, 2016; v1 submitted 17 February, 2016; originally announced February 2016.

Comments: Published in ApJL. Updated version contains a different non-transiting HJ candidate than the original

Journal ref: ApJL, 823, L7 (2016)