Showing 1–50 of 74 results for author: Murray-Clay, R

Extensive Pollution of Uranus and Neptune's Atmospheres by Upsweep of Icy Material During the Nice Model Migration

Authors: Eva Zlimen, Elizabeth Bailey, Ruth Murray-Clay

Abstract: In the Nice model of solar system formation, Uranus and Neptune undergo an orbital upheaval, sweeping through a planetesimal disk. The region of the disk from which material is accreted by the ice giants during this phase of their evolution has not previously been identified. We perform direct N-body orbital simulations of the four giant planets to determine the amount and origin of solid accretio… ▽ More In the Nice model of solar system formation, Uranus and Neptune undergo an orbital upheaval, sweeping through a planetesimal disk. The region of the disk from which material is accreted by the ice giants during this phase of their evolution has not previously been identified. We perform direct N-body orbital simulations of the four giant planets to determine the amount and origin of solid accretion during this orbital upheaval. We find that the ice giants undergo an extreme bombardment event, with collision rates as much as ~3 per hour assuming km-sized planetesimals, increasing the total planet mass by up to ~0.35%. In all cases, the initially outermost ice giant experiences the largest total enhancement. We determine that for some plausible planetesimal properties, the resulting atmospheric enrichment could potentially produce sufficient latent heat to alter the planetary cooling timescale according to existing models. Our findings suggest that substantial accretion during this phase of planetary evolution may have been sufficient to impact the atmospheric composition and thermal evolution of the ice giants, motivating future work on the fate of deposited solid material. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 19 pages, 9 figures; accepted to AJ

Randomness and Retention: Using Weak Mean Motion Resonances to Constrain Neptune's Late-Stage Migration

Authors: Arcelia Hermosillo Ruiz, Harriet C. P. Lau, Ruth Murray-Clay

Abstract: Planet-planetesimal interactions cause a planet to migrate, manifesting as a random walk in semi-major axis. In models for Neptune's migration involving a gravitational upheaval, this planetesimal-driven migration is a side-effect of the dynamical friction required to damp Neptune's orbital eccentricitiy. This migration is noisy, potentially causing Trans Neptunian Objects (TNOs) in mean motion re… ▽ More Planet-planetesimal interactions cause a planet to migrate, manifesting as a random walk in semi-major axis. In models for Neptune's migration involving a gravitational upheaval, this planetesimal-driven migration is a side-effect of the dynamical friction required to damp Neptune's orbital eccentricitiy. This migration is noisy, potentially causing Trans Neptunian Objects (TNOs) in mean motion resonance to be lost. With Nbody simulations, we validate a previously-derived analytic model for resonance retention and determine unknown coefficients. We identify the impact of random-walk (noisy) migration on resonance retention for resonances up to fourth order lying between 39 au and 75 au. Using a population estimate for the weak 7:3 resonance from the well-characterized Outer Solar System Origins Survey (OSSOS), we rule out two cases: (1) a planetesimal disk distributed between 13.3 and 39.9 au with $\gtrsim$ 30 Earth masses in today's size distribution and $T_{\rm mig} \gtrsim$ 40Myr and (2) a top-heavy size distribution with $\gtrsim$ 2000 Pluto-sized TNOs and $T_{\rm mig} \gtrsim$ 10Myr, where $T_{\rm mig}$ is Neptune's migration timescale. We find that low-eccentricity TNOs in the heavily populated 5:2 resonance are easily lost due to noisy migration. Improved observations of the low-eccentricity region of the 5:2 resonance and of weak mean motion resonances with Rubin Observatory's Legacy Survey of Space and Time (LSST) will provide better population estimates, allowing for comparison with our model's retention fractions and providing strong evidence for or against Neptune's random interactions with planetesimals. △ Less

Submitted 10 May, 2024; originally announced May 2024.

Comments: 17 pages, 11 figures; accepted to MNRAS

In Violation of the Prime Directive: Simulating detriments to Delta-Quadrant civilizations from the starship Voyager's impact on planetary rings

Authors: J. Fowler, Ruth Murray-Clay

Abstract: In the seven years that the starship Voyager spent in the Delta Quadrant, it used many questionable techniques to engage with alien civilizations and ultimately find its way home. From detailed studies of their logs and opening credits, we simulate Voyager's practice of orbiting a planet, to examine the effect on planetary rings. We outline a feasible planetary system and simulate the extent to wh… ▽ More In the seven years that the starship Voyager spent in the Delta Quadrant, it used many questionable techniques to engage with alien civilizations and ultimately find its way home. From detailed studies of their logs and opening credits, we simulate Voyager's practice of orbiting a planet, to examine the effect on planetary rings. We outline a feasible planetary system and simulate the extent to which its rings would be disrupted. We find that Voyager's orbit could inflate the height of the rings in the vicinity of the spacecraft by a factor of 2, as well as increase the relative speeds of neighboring planetesimals within the rings. This increase in ring thickness has the potential to alter shadows on any moons of this planet, impacting ring-shadow based religions. Additionally, the acceleration of these planetesimals could rival their gravity, bucking any alien inhabitants and their tiny civilizations off of their planetesimal homeworlds. Finally, we posit that due to increased collisions amongst the planetesimals (which may harbor tiny intelligent life) the trajectory of these civilizations may be forever altered, violating the prime directive. △ Less

Submitted 29 March, 2024; originally announced March 2024.

A Pathway for Collisional Planetesimal Growth in the Ice-Dominant Regions of Protoplanetary Disks

Authors: Elizabeth Yunerman, Diana Powell, Ruth Murray-Clay

Abstract: We present a semi-analytic model for the growth, drift, desorption, and fragmentation of millimeter- to meter-sized particles in protoplanetary disks. Fragmentation occurs where particle collision velocities exceed critical fragmentation velocities. Using this criterion, we produce fragmentation regions in disk orbital radius-particle size phase space for particles with a range of material propert… ▽ More We present a semi-analytic model for the growth, drift, desorption, and fragmentation of millimeter- to meter-sized particles in protoplanetary disks. Fragmentation occurs where particle collision velocities exceed critical fragmentation velocities. Using this criterion, we produce fragmentation regions in disk orbital radius-particle size phase space for particles with a range of material properties, structures, and compositions (including SiO$_2$, Mg$_2$SiO$_4$, H$_2$O, CO$_2$, and CO). For reasonable disk conditions, compact aggregate H$_2$O, CO$_2$, and CO ice particles do not reach destructive relative velocities and are thus not likely to undergo collisional fragmentation. Uncoated silicate particles are more susceptible to collisional destruction and are expected to fragment in the inner disk, consistent with previous work. We then calculate the growth, drift, and sublimation of small particles, initially located in the outer disk. We find that ice-coated particles can avoid fragmentation as they grow and drift inward under a substantial range of disk conditions as long as the particles are aggregates composed of 0.1 $μ$m-sized monomers. Such particles may undergo runaway growth in disk regions abundant in H$_2$O or CO$_2$ ice depending on the assumed disk temperature structure. These results indicate that icy collisional growth to planetesimally-relevant sizes may happen efficiently throughout a disk's lifetime, and is particularly robust at early times when the disk's dust-to-gas ratio is comparable to that of the interstellar medium. △ Less

Submitted 9 November, 2023; v1 submitted 7 November, 2023; originally announced November 2023.

Comments: 26 pages, 14 figures, accepted to ApJ

Using Lyman-$α$ transits to constrain models of atmospheric escape

Authors: Ethan Schreyer, James. E. Owen, R. O. Parke Loyd, Ruth Murray-Clay

Abstract: Lyman-$α$ transits provide an opportunity to test models of atmospheric escape directly. However, translating observations into constraints on the properties of the escaping atmosphere is challenging. The major reason for this is that the observable parts of the outflow often comes from material outside the planet's Hill sphere, where the interaction between the planetary outflow and circumstellar… ▽ More Lyman-$α$ transits provide an opportunity to test models of atmospheric escape directly. However, translating observations into constraints on the properties of the escaping atmosphere is challenging. The major reason for this is that the observable parts of the outflow often comes from material outside the planet's Hill sphere, where the interaction between the planetary outflow and circumstellar environment is important. As a result, 3D models are required to match observations. Whilst 3D hydrodynamic simulations are able to match observational features qualitatively, they are too computationally expensive to perform a statistical retrieval of properties of the outflow. Here, we develop a model that determines the trajectory, ionization state and 3D geometry of the outflow as a function of its properties and system parameters. We then couple this model to a ray tracing routine in order to produce synthetic transits. We demonstrate the validity of this approach, reproducing the trajectory of the outflows seen in 3D simulations. We illustrate the use of this model by performing a retrieval on the transit spectrum of GJ 436 b. Our model constrains the sound speed of the outflow $\gtrsim 10 \text{ km s}^{-1}$, indicating that we can rule out core-powered mass loss as the mechanism driving the outflow for this planet. The bound on planetary outflow velocity and mass loss rates are consistent with a photoevaporative wind. △ Less

Submitted 27 October, 2023; originally announced October 2023.

High Angular Resolution Imaging of the V892 Tau Binary System: A New Circumprimary Disk Detection and Updated Orbital Constraints

Authors: Christina Vides, Steph Sallum, Josh Eisner, Andy Skemer, Ruth Murray-Clay

Abstract: We present a direct imaging study of V892 Tau, a young Herbig Ae/Be star with a close-in stellar companion and circumbinary disk. Our observations consist of images acquired via Keck 2/NIRC2 with non-redundant masking and the pyramid wavefront sensor at K$^\prime$ band (2.12$μ$m) and L$^\prime$ band (3.78$μ$m). Sensitivity to low-mass accreting companions and cool disk material is high at L… ▽ More We present a direct imaging study of V892 Tau, a young Herbig Ae/Be star with a close-in stellar companion and circumbinary disk. Our observations consist of images acquired via Keck 2/NIRC2 with non-redundant masking and the pyramid wavefront sensor at K$^\prime$ band (2.12$μ$m) and L$^\prime$ band (3.78$μ$m). Sensitivity to low-mass accreting companions and cool disk material is high at L$^\prime$ band, while complimentary observations at K$^\prime$ band probe hotter material with higher angular resolution. These multi-wavelength, multi-epoch data allow us to differentiate the secondary stellar emission from disk emission and deeply probe the structure of the circumbinary disk at small angular separations. We constrain architectural properties of the system by fitting geometric disk and companion models to the K$^\prime$ and L$^\prime$ band data. From these models, we constrain the astrometric and photometric properties of the stellar binary and update the orbit, placing the tightest estimates to date on the V892 Tau orbital parameters. We also constrain the geometric structure of the circumbinary disk, and resolve a circumprimary disk for the first time. △ Less

Submitted 13 October, 2023; v1 submitted 3 October, 2023; originally announced October 2023.

Comments: Accepted to ApJ

Can the orbital distribution of Neptune's 3:2 mean motion resonance result from stability sculpting?

Authors: Sricharan Balaji, Nihaal Zaveri, Nanae Hayashi, Arcelia Hermosillo Ruiz, Jackson Barnes, Ruth Murray-Clay, Kathryn Volk, Jake Gerhardt, Zain Syed

Abstract: We explore a simplified model of the outcome of an early outer Solar System gravitational upheaval during which objects were captured into Neptune's 3:2 mean motion resonance via scattering rather than smooth planetary migration. We use N-body simulations containing the Sun, the four giant planets, and test particles in the 3:2 resonance to determine whether long-term stability sculpting over 4.5… ▽ More We explore a simplified model of the outcome of an early outer Solar System gravitational upheaval during which objects were captured into Neptune's 3:2 mean motion resonance via scattering rather than smooth planetary migration. We use N-body simulations containing the Sun, the four giant planets, and test particles in the 3:2 resonance to determine whether long-term stability sculpting over 4.5 Gyr can reproduce the observed 3:2 resonant population from an initially randomly scattered 3:2 population. After passing our simulated 3:2 resonant objects through a survey simulator, we find that the semimajor axis (a) and eccentricity (e) distributions are consistent with the observational data (assuming an absolute magnitude distribution constrained by prior studies), suggesting that these could be a result of stability sculpting. However, the inclination (i) distribution cannot be produced be stability sculpting and thus must result from a distinct process that excited the inclinations. Our simulations modestly under-predict the number of objects with high libration amplitudes (Aφ), possibly because we do not model transient sticking. Finally, our model under-populates the Kozai subresonance compared to both observations and to smooth migration models. Future work is needed to determine whether smooth migration occurring as Neptune's eccentricity damped to its current value can resolve this discrepancy. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: 15 pages, 13 figures

Systematic Multi-Epoch Monitoring of LkCa 15: Dynamic Dust Structures on Solar-System Scales

Authors: Steph Sallum, Josh Eisner, Andy Skemer, Ruth Murray-Clay

Abstract: We present the highest angular resolution infrared monitoring of LkCa 15, a young solar analog hosting a transition disk. This system has been the subject of a number of direct imaging studies from the millimeter through the optical, which have revealed multiple protoplanetary disk rings as well as three orbiting protoplanet candidates detected in infrared continuum (one of which was simultaneousl… ▽ More We present the highest angular resolution infrared monitoring of LkCa 15, a young solar analog hosting a transition disk. This system has been the subject of a number of direct imaging studies from the millimeter through the optical, which have revealed multiple protoplanetary disk rings as well as three orbiting protoplanet candidates detected in infrared continuum (one of which was simultaneously seen at H$α$). We use high-angular-resolution infrared imaging from 2014-2020 to systematically monitor these infrared signals and determine their physical origin. We find that three self-luminous protoplanets cannot explain the positional evolution of the infrared sources, since the longer time baseline images lack the coherent orbital motion that would be expected for companions. However, the data still strongly prefer a time-variable morphology that cannot be reproduced by static scattered-light disk models. The multi-epoch observations suggest the presence of complex and dynamic substructures moving through the forward-scattering side of the disk at $\sim20$ AU, or quickly-varying shadowing by closer-in material. We explore whether the previous H$α$ detection of one candidate would be inconsistent with this scenario, and in the process develop an analytical signal-to-noise penalty for H$α$ excesses detected near forward-scattered light. Under these new noise considerations, the H$α$ detection is not strongly inconsistent with forward scattering, making the dynamic LkCa 15 disk a natural explanation for both the infrared and H$α$ data. △ Less

Submitted 20 July, 2023; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: 24 pages, 11 figures, accepted for publication in ApJ

Depletion of gaseous CO in protoplanetary disks by surface-energy-regulated ice formation

Authors: Diana Powell, Peter Gao, Ruth Murray-Clay, Xi Zhang

Abstract: Empirical constraints of fundamental properties of protoplanetary disks are essential for understanding planet formation and planetary properties (1,2). Carbon monoxide (CO) gas is often used to constrain disk properties (3). However, estimates show that the CO gas abundance in disks is depleted relative to expected values (4,5,6,7), and models of various disk processes impacting the CO abundance… ▽ More Empirical constraints of fundamental properties of protoplanetary disks are essential for understanding planet formation and planetary properties (1,2). Carbon monoxide (CO) gas is often used to constrain disk properties (3). However, estimates show that the CO gas abundance in disks is depleted relative to expected values (4,5,6,7), and models of various disk processes impacting the CO abundance could not explain this depletion on observed 1Myr timescales (8,9,10,11,12,13,14). Here we demonstrate that surface energy effects on particles in disks, such as the Kelvin effect, that arise when ice heterogeneously nucleates onto an existing particle can efficiently trap CO in its ice phase. In previous ice formation models, CO gas was released when small ice-coated particles were lofted to warmed disk layers. Our model can reproduce the observed abundance, distribution and time evolution of gaseous CO in the four most studied protoplanetary disks (7). We constrain the solid and gaseous CO inventory at the midplane and disk diffusivities and resolve inconsistencies in estimates of the disk mass -- three crucial parameters that control planetary formation. △ Less

Submitted 29 August, 2022; originally announced August 2022.

Comments: Published in Nature Astronomy (2022), includes supplementary material

The UV-SCOPE Mission: Ultraviolet Spectroscopic Characterization Of Planets and their Environments

Authors: David R. Ardila, Evgenya Shkolnik, John Ziemer, Mark Swain, James E. Owen, Michael Line, R. O. Parke Loyd, R. Glenn Sellar, Travis Barman, Courtney Dressing, William Frazier, April D. Jewell, Robert J. Kinsey, Carl C. Liebe, Joshua D. Lothringer, Luz Maria Martinez-Sierra, James McGuire, Victoria Meadows, Ruth Murray-Clay, Shouleh Nikzad, Sarah Peacock, Hilke Schlichting, David Sing, Kevin Stevenson, Yen-Hung Wu

Abstract: UV-SCOPE is a mission concept to determine the causes of atmospheric mass loss in exoplanets, investigate the mechanisms driving aerosol formation in hot Jupiters, and study the influence of the stellar environment on atmospheric evolution and habitability. As part of these investigations, the mission will generate a broad-purpose legacy database of time-domain ultraviolet (UV) spectra for nearly… ▽ More UV-SCOPE is a mission concept to determine the causes of atmospheric mass loss in exoplanets, investigate the mechanisms driving aerosol formation in hot Jupiters, and study the influence of the stellar environment on atmospheric evolution and habitability. As part of these investigations, the mission will generate a broad-purpose legacy database of time-domain ultraviolet (UV) spectra for nearly 200 stars and planets. The observatory consists of a 60 cm, f/10 telescope paired to a long-slit spectrograph, yielding simultaneous, almost continuous coverage between 1203 Å and 4000 Å, with resolutions ranging from 6000 to 240. The efficient instrument provides throughputs > 4% (far-UV; FUV) and > 15% (near-UV; NUV), comparable to HST/COS and much better than HST/STIS, over the same spectral range. A key design feature is the LiF prism, which serves as a dispersive element and provides high throughput even after accounting for radiation degradation. The use of two delta-doped Electron-Multiplying CCD detectors with UV-optimized, single-layer anti-reflection coatings provides high quantum efficiency and low detector noise. From the Earth-Sun second Lagrangian point, UV-SCOPE will continuously observe planetary transits and stellar variability in the full FUV-to-NUV range, with negligible astrophysical background. All these features make UV-SCOPE the ideal instrument to study exoplanetary atmospheres and the impact of host stars on their planets. UV-SCOPE was proposed to NASA as a Medium Explorer (MidEx) mission for the 2021 Announcement of Opportunity. If approved, the observatory will be developed over a 5-year period. Its primary science mission takes 34 months to complete. The spacecraft carries enough fuel for 6 years of operations. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Comments: 12 pages, 7 figures, 1 table. Conference presentation, 17 July 2022, SPIE Astronomical Telescopes and Instrumentation, Montreal, Canada

Effects of Charge Exchange on the Evaporative Wind of HD 209458b

Authors: Alex Debrecht, Jonathan Carroll-Nellenback, Adam Frank, Eric G. Blackman, Luca Fossati, Ruth Murray-Clay, John McCann

Abstract: The role of charge exchange in shaping exoplanet photoevaporation remains a topic of contention. Exchange of electrons between stellar wind protons from the exoplanet's host star and neutral hydrogen from the planet's wind has been proposed as a mechanism to create "energetic neutral atoms" (ENAs), which could explain the high absorption line velocities observed in systems where mass loss is occur… ▽ More The role of charge exchange in shaping exoplanet photoevaporation remains a topic of contention. Exchange of electrons between stellar wind protons from the exoplanet's host star and neutral hydrogen from the planet's wind has been proposed as a mechanism to create "energetic neutral atoms" (ENAs), which could explain the high absorption line velocities observed in systems where mass loss is occurring. In this paper we present results from 3D hydrodynamic simulations of the mass loss of a planet similar to HD 209458b. We self-consistently launch a planetary wind by calculating the ionization and heating resulting from incident high-energy radiation, inject a stellar wind into the simulation, and allow electron exchange between the stellar and planetary winds. We predict the potential production of ENAs by the wind-wind interaction analytically, then present the results of our simulations, which confirm the analytic limits. Within the limits of our hydrodynamic simulation, we find that charge exchange with the stellar wind properties examined here is unable to explain the absorption observed at high Doppler velocities. △ Less

Submitted 18 January, 2022; originally announced January 2022.

Comments: arXiv admin note: text overlap with arXiv:1906.00075

The fundamentals of Lyman-alpha exoplanet transits

Authors: James E. Owen, Ruth A. Murray-Clay, Ethan Schreyer, Hilke E. Schlichting, David Ardila, Akash Gupta, R. O. Parke Loyd, Evgenya L. Shkolnik, David K. Sing, Mark R. Swain

Abstract: Lyman-$α$ transits have been detected from several nearby exoplanets and are one of our best insights into the atmospheric escape process. However, due to ISM absorption, we typically only observe the transit signature in the blue-wing, making them challenging to interpret. This challenge has been recently highlighted by non-detections from planets thought to be undergoing vigorous escape. Pioneer… ▽ More Lyman-$α$ transits have been detected from several nearby exoplanets and are one of our best insights into the atmospheric escape process. However, due to ISM absorption, we typically only observe the transit signature in the blue-wing, making them challenging to interpret. This challenge has been recently highlighted by non-detections from planets thought to be undergoing vigorous escape. Pioneering 3D simulations have shown that escaping hydrogen is shaped into a cometary tail receding from the planet. Motivated by this work, we develop a simple model to interpret Lyman-$α$ transits. Using this framework, we show that the Lyman-$α$ transit depth is primarily controlled by the properties of the stellar tidal field rather than details of the escape process. Instead, the transit duration provides a direct measurement of the velocity of the planetary outflow. This result arises because the underlying physics is the distance a neutral hydrogen atom can travel before it is photoionized in the outflow. Thus, higher irradiation levels, expected to drive more powerful outflows, produce weaker, shorter Lyman-$α$ transits because the outflowing gas is ionized more quickly. Our framework suggests that the generation of energetic neutral atoms may dominate the transit signature early, but the acceleration of planetary material produces long tails. Thus, Lyman-$α$ transits do not primarily probe the mass-loss rates. Instead, they inform us about the velocity at which the escape mechanism is ejecting material from the planet, providing a clean test of predictions from atmospheric escape models. △ Less

Submitted 18 November, 2022; v1 submitted 11 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS

OSSOS: The Eccentricity and Inclination Distributions of the Stable Neptunian Trojans

Authors: Hsing Wen Lin, Ying-Tung Chen, Kathryn Volk, Brett Gladman, Ruth Murray-Clay, Mike Alexandersen, Michele T. Bannister, Samantha M. Lawler, Wing-Huen Ip, Patryk Sofia Lykawka, J. J. Kavelaars, Stephen D. J. Gwyn, Jean-Marc Petit

Abstract: The minor planets on orbits that are dynamically stable in Neptune's 1:1 resonance on Gyr timescales were likely emplaced by Neptune's outward migration. We explore the intrinsic libration amplitude, eccentricity, and inclination distribution of Neptune's stable Trojans, using the detections and survey efficiency of the Outer Solar System Origins Survey (OSSOS) and Pan-STARRS1. We find that the li… ▽ More The minor planets on orbits that are dynamically stable in Neptune's 1:1 resonance on Gyr timescales were likely emplaced by Neptune's outward migration. We explore the intrinsic libration amplitude, eccentricity, and inclination distribution of Neptune's stable Trojans, using the detections and survey efficiency of the Outer Solar System Origins Survey (OSSOS) and Pan-STARRS1. We find that the libration amplitude of the stable Neptunian Trojan population can be well modeled as a Rayleigh distribution with a libration amplitude width $σ_{A_phi}$ of 15$^\circ$. When taken as a whole, the Neptune Trojan population can be acceptably modeled with a Rayleigh eccentricity distribution of width $σ_e$ of 0.045 and a typical sin(i) x Gaussian inclination distribution with a width $σ_i$ of 14 +/- 2 degrees. However, these distributions are only marginally acceptable. This is likely because, even after accounting for survey detection biases, the known large Hr < 8 and small Hr >= 8 Neptune Trojans appear to have markedly different eccentricities and inclinations. We propose that like the classical Kuiper belt, the stable intrinsic Neptunian Trojan population have dynamically `hot' and dynamically `cold' components to its eccentricity/inclination distribution, with $σ_{e-cold}$ ~ 0.02 / $σ_{i-cold}$ ~ 6$^\circ$ and $σ_{e-hot}$~ 0.05 / $σ_{i-hot}$ ~ 18$^\circ$. In this scenario, the `cold' L4 Neptunian Trojan population lacks the Hr >= 8 members and has 13 +11/-6 `cold' Trojans with Hr < 8. On the other hand, the `hot' L4 Neptunian Trojan population has 136 +57/-48 Trojans with Hr < 10 -- a population 2.4 times greater than that of the L4 Jovian Trojans in the same luminosity range. △ Less

Submitted 15 February, 2021; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: 15 pages, 7 figures, v3, accepted for publication in Icarus

How Consumption and Repulsion Set Planetary Gap Depths and the Final Masses of Gas Giants

Authors: Mickey Rosenthal, Eugene Chiang, Sivan Ginzburg, Ruth Murray-Clay

Abstract: Planets open gaps in discs. Gap opening is typically modeled by considering the planetary Lindblad torque which repels disc gas away from the planet's orbit. But gaps also clear because the planet consumes local material. We present a simple, easy-to-use, analytic framework for calculating how gaps deplete and how the disc's structure as a whole changes by the combined action of Lindblad repulsion… ▽ More Planets open gaps in discs. Gap opening is typically modeled by considering the planetary Lindblad torque which repels disc gas away from the planet's orbit. But gaps also clear because the planet consumes local material. We present a simple, easy-to-use, analytic framework for calculating how gaps deplete and how the disc's structure as a whole changes by the combined action of Lindblad repulsion and planetary consumption. The final mass to which a gap-embedded gas giant grows is derived in tandem. The analytics are tested against 1D numerical experiments and calibrated using published multi-dimensional simulations. In viscous alpha discs, the planet, while clearing a gap, initially accretes practically all of the gas that tries to diffuse past, rapidly achieving super-Jupiter if not brown dwarf status. By contrast, in inviscid discs---that may still accrete onto their central stars by, say, magnetized winds---planets open deep, repulsion-dominated gaps. Then only a small fraction of the disc accretion flow is diverted onto the planet, which grows to a fraction of a Jupiter mass. Transitional disc cavities might be cleared by families of such low-mass objects opening inviscid, repulsion-dominated, overlapping gaps which allow most of the outer disc gas to flow unimpeded onto host stars. △ Less

Submitted 12 June, 2020; v1 submitted 28 April, 2020; originally announced April 2020.

Comments: Accepted for publication in MNRAS

Requirements for gravitational collapse in planetesimal formation --- the impact of scales set by Kelvin-Helmholtz and nonlinear streaming instability

Authors: Konstantin Gerbig, Ruth A. Murray-Clay, Hubert Klahr, Hans Baehr

Abstract: The formation of planetesimals is an unsolved problem in planet formation theory. A prominent scenario for overcoming dust growth barriers in dead zones is the gravitational collapse of locally over-dense regions, shown to robustly produce $\sim$100 km sized objects. Still, the conditions under which planetesimal formation occurs remain unclear. For collapse to proceed, the self-gravity of an over… ▽ More The formation of planetesimals is an unsolved problem in planet formation theory. A prominent scenario for overcoming dust growth barriers in dead zones is the gravitational collapse of locally over-dense regions, shown to robustly produce $\sim$100 km sized objects. Still, the conditions under which planetesimal formation occurs remain unclear. For collapse to proceed, the self-gravity of an overdensity must overcome stellar tidal disruption on large scales and turbulent diffusion on small scales. Here, we relate the scales of streaming and Kelvin-Helmholtz instability, which both regulate particle densities on the scales of gravitational collapse, directly to planetesimal formation. We support our analytic findings by performing 3D hydrodynamical simulations of streaming and Kelvin-Helmholtz instability and planetesimal formation. We find that the vertical extent of the particle mid-plane layer and the radial width of streaming instability filaments are set by the same characteristic length scale, thus governing the strength of turbulent diffusion on the scales of planetesimal formation. We present and successfully test a collapse criterion: $0.1 Q βε^{-1}Z^{-1} \lesssim 1$ and show that even for Solar metallicities, planetesimals can form in dead zones of sufficiently massive disks. For a given gas Toomre-parameter $Q$, pressure gradient $β$, metallicity $Z$ and local particle enhancement $ε$, the collapse criterion also provides a range of unstable scales, instituting a promising path for studying initial planetesimal mass distributions. Streaming instability is not required for planetesimal collapse, but by increasing $ε$, can evolve a system to instability. △ Less

Submitted 28 January, 2020; originally announced January 2020.

Comments: 26 pages, 11 figures, submitted to ApJ

OSSOS XVIII: Constraining migration models with the 2:1 resonance using the Outer Solar System Origins Survey

Authors: Ying-Tung Chen, Brett Gladman, Kathryn Volk, Ruth Murray-Clay, Matthew J. Lehner, J. J. Kavelaars, Shiang-Yu Wang, Hsing-Wen Lin, Patryk Sofia Lykawka, Mike Alexandersen, Michele T. Bannister, Samantha M. Lawler, Rebekah I. Dawson, Sarah Greenstreet, Stephen D. J. Gwyn, Jean-Marc Petit

Abstract: Resonant dynamics plays a significant role in the past evolution and current state of our outer Solar System. The population ratios and spatial distribution of Neptune's resonant populations are direct clues to understanding the history of our planetary system. The orbital structure of the objects in Neptune's 2:1 mean-motion resonance (\emph{twotinos}) has the potential to be a tracer of planetar… ▽ More Resonant dynamics plays a significant role in the past evolution and current state of our outer Solar System. The population ratios and spatial distribution of Neptune's resonant populations are direct clues to understanding the history of our planetary system. The orbital structure of the objects in Neptune's 2:1 mean-motion resonance (\emph{twotinos}) has the potential to be a tracer of planetary migration processes. Different migration processes produce distinct architectures, recognizable by well-characterized surveys. However, previous characterized surveys only discovered a few twotinos, making it impossible to model the intrinsic twotino population. With a well-designed cadence and nearly 100\% tracking success, the Outer Solar System Origins Survey (OSSOS) discovered 838 trans-Neptunian objects, of which 34 are securely twotinos with well-constrained libration angles and amplitudes. We use the OSSOS twotinos and the survey characterization parameters via the OSSOS Survey Simulator to inspect the intrinsic population and orbital distributions of twotino. The estimated twotino population, 4400$^{+1500}_{-1100}$ with $H_r<8.66$ (diameter$\sim$100km) at 95\% confidence, is consistent with the previous low-precision estimate. We also constrain the width of the inclination distribution to a relatively narrow value of $σ_i$=6$^\circ$$^{+1}_{-1}$, and find the eccentricity distribution is consistent with a Gaussian centered on $e_\mathrm{c}=0.275$ with a width $e_\mathrm{w}=0.06$. We find a single-slope exponential luminosity function with $α=0.6$ for the twotinos. Finally, we for the first time meaningfully constrain the fraction of symmetric twotinos, and the ratio of the leading asymmetric islands; both fractions are in a range of 0.2--0.6. These measurements rule out certain theoretical models of Neptune's migration history. △ Less

Submitted 31 October, 2019; v1 submitted 28 October, 2019; originally announced October 2019.

Comments: 23 pages, 12 figures, 2 tables, accepted for publication in the Astronomical Journal

How Flow Isolation May Set the Mass Scale for Super-Earth Planets

Authors: M. M. Rosenthal, R. A. Murray-Clay

Abstract: Much recent work on planet formation has focused on the growth of planets by accretion of grains whose aerodynamic properties make them marginally coupled to the nebular gas, a theory commonly referred to as "pebble accretion". While pebble accretion can ameliorate some of the issues presented by growth by purely gravitational processes, it has other issues when compared with observations of exopl… ▽ More Much recent work on planet formation has focused on the growth of planets by accretion of grains whose aerodynamic properties make them marginally coupled to the nebular gas, a theory commonly referred to as "pebble accretion". While pebble accretion can ameliorate some of the issues presented by growth by purely gravitational processes, it has other issues when compared with observations of exoplanetary systems. A particular concern is the preponderance of planets that end their growth as "super-Earths" or "sub-Neptunes", with masses in the range 2-10 $M_\oplus$. Once planets reach this mass scale, timescales for growth by pebble accretion are so rapid that ubiquitously ending growth here is difficult. In this work, we highlight this issue in detail using our previously published model of pebble accretion, and also propose a possible solution: feedback between the growing planet's atmosphere and the gas disk inhibits accretion of smaller particle sizes by forcing them to flow around the growing planet instead of being accreted. For reasonable fiducial disk parameters this "flow isolation" will inhibit accretion of all available particle sizes once the planet reaches super-Earth masses. We also demonstrate that the characteristics of this "flow isolation mass" agree with previously published trends identified in the \textit{Kepler} planets. △ Less

Submitted 19 August, 2020; v1 submitted 19 August, 2019; originally announced August 2019.

Comments: 14 pages, 6 figures

Measuring the Orbital Parameters of Radial Velocity Systems in Mean Motion Resonance---a Case Study of HD 200964

Authors: M. M. Rosenthal, W. Jacobson-Galan, B. Nelson, R. A. Murray-Clay, J. A. Burt, B. Holden, E. Chang, N. Kaaz, J. Yant, R. P. Butler, S. S. Vogt

Abstract: The presence of mean motion resonances (MMRs) complicates analysis and fitting of planetary systems observed through the radial velocity (RV) technique. MMR can allow planets to remain stable in regions of phase space where strong planet-planet interactions would otherwise destabilize the system. These stable orbits can occupy small phase space volumes, allowing MMRs to strongly constrain system p… ▽ More The presence of mean motion resonances (MMRs) complicates analysis and fitting of planetary systems observed through the radial velocity (RV) technique. MMR can allow planets to remain stable in regions of phase space where strong planet-planet interactions would otherwise destabilize the system. These stable orbits can occupy small phase space volumes, allowing MMRs to strongly constrain system parameters, but making searches for stable orbital parameters challenging. Furthermore, libration of the resonant angle and dynamical interaction between the planets introduces another, long period variation into the observed RV signal, complicating analysis of the periods of the planets in the system. We discuss this phenomenon using the example of HD 200964. By searching through parameter space and numerically integrating each proposed set of planetary parameters to test for long term stability, we find stable solutions in the 7:5 and 3:2 MMRs in addition to the originally identified 4:3 MMR. The 7:5 configuration provides the best match to the data, while the 3:2 configuration provides the most easily understood formation scenario. In reanalysis of the originally published shorter-baseline data, we find fits in both the 4:3 and 3:2 resonances, but not the 7:5. Because the time baseline of the data is less than the resonant libration period, the current best fit to the data may not reflect the actual resonant configuration. In the absence of a full sample of the longer libration period, we find that it is of paramount importance to incorporate long term stability when fitting for the system's orbital configuration. △ Less

Submitted 13 August, 2019; originally announced August 2019.

Comments: 26 pages, 23 figures

Signatures of a planet-planet impacts phase in exoplanetary systems hosting giant planets

Authors: Renata Frelikh, Hyerin Jang, Ruth A. Murray-Clay, Cristobal Petrovich

Abstract: Exoplanetary systems host giant planets on substantially non-circular, close-in orbits. We propose that these eccentricities arise in a phase of giant impacts, analogous to the final stage of Solar System assembly that formed Earth's Moon. In this scenario, the planets scatter each other and collide, with corresponding mass growth as they merge. We numerically integrate an ensemble of systems with… ▽ More Exoplanetary systems host giant planets on substantially non-circular, close-in orbits. We propose that these eccentricities arise in a phase of giant impacts, analogous to the final stage of Solar System assembly that formed Earth's Moon. In this scenario, the planets scatter each other and collide, with corresponding mass growth as they merge. We numerically integrate an ensemble of systems with varying total planet mass, allowing for collisional growth, to show that (1) the high-eccentricity giants observed today may have formed preferentially in systems of higher initial total planet mass, and (2) the upper bound on the observed giant planet eccentricity distribution is consistent with planet-planet scattering. We predict that mergers will produce a population of high-mass giant planets between 1 and 8 au from their stars. △ Less

Submitted 8 October, 2019; v1 submitted 7 June, 2019; originally announced June 2019.

Comments: 8 pages, 5 figures, accepted to ApJL

Effects of Radiation Pressure on the Evaporative Wind of HD 209458b

Authors: Alex Debrecht, Jonathan Carroll-Nellenback, Adam Frank, Eric G. Blackman, Luca Fossati, John McCann, Ruth Murray-Clay

Abstract: The role of radiation pressure in shaping exoplanet photoevaporation remains a topic of contention. Radiation pressure from the exoplanet's host star has been proposed as a mechanism to drive the escaping atmosphere into a "cometary" tail and explain the high velocities observed in systems where mass loss is occurring. In this paper we present results from high-resolution 3-D hydrodynamic simulati… ▽ More The role of radiation pressure in shaping exoplanet photoevaporation remains a topic of contention. Radiation pressure from the exoplanet's host star has been proposed as a mechanism to drive the escaping atmosphere into a "cometary" tail and explain the high velocities observed in systems where mass loss is occurring. In this paper we present results from high-resolution 3-D hydrodynamic simulations of a planet similar to HD 209458b. We self-consistently launch a wind flowing outward from the planet by calculating the ionization and heating resulting from incident high-energy radiation, and account for radiation pressure. We first present a simplified calculation, setting a limit on the Lyman-$α$ flux required to drive the photo-evaporated planetary material to larger radii and line-of-sight velocities. We then present the results of our simulations, which confirm the limits determined by our analytic calculation. We thus demonstrate that, within the limits of our hydrodynamic simulation and for the Lyman-$α$ fluxes expected for HD 209458, radiation pressure is unlikely to significantly affect photoevaporative winds or to explain the high velocities at which wind material is observed, though further possibilities remain to be investigated. △ Less

Submitted 3 February, 2020; v1 submitted 31 May, 2019; originally announced June 2019.

Directly Imaging Rocky Planets from the Ground

Authors: B. Mazin, É. Artigau, V. Bailey, C. Baranec, C. Beichman, B. Benneke, J. Birkby, T. Brandt, J. Chilcote, M. Chun, L. Close, T. Currie, I. Crossfield, R. Dekany, J. R. Delorme, C. Dong, R. Dong, R. Doyon, C. Dressing, M. Fitzgerald, J. Fortney, R. Frazin, E. Gaidos, O. Guyon, J. Hashimoto , et al. (38 additional authors not shown)

Abstract: Over the past three decades instruments on the ground and in space have discovered thousands of planets outside the solar system. These observations have given rise to an astonishingly detailed picture of the demographics of short-period planets, but are incomplete at longer periods where both the sensitivity of transit surveys and radial velocity signals plummet. Even more glaring is that the spe… ▽ More Over the past three decades instruments on the ground and in space have discovered thousands of planets outside the solar system. These observations have given rise to an astonishingly detailed picture of the demographics of short-period planets, but are incomplete at longer periods where both the sensitivity of transit surveys and radial velocity signals plummet. Even more glaring is that the spectra of planets discovered with these indirect methods are either inaccessible (radial velocity detections) or only available for a small subclass of transiting planets with thick, clear atmospheres. Direct detection can be used to discover and characterize the atmospheres of planets at intermediate and wide separations, including non-transiting exoplanets. Today, a small number of exoplanets have been directly imaged, but they represent only a rare class of young, self-luminous super-Jovian-mass objects orbiting tens to hundreds of AU from their host stars. Atmospheric characterization of planets in the <5 AU regime, where radial velocity (RV) surveys have revealed an abundance of other worlds, is technically feasible with 30-m class apertures in combination with an advanced AO system, coronagraph, and suite of spectrometers and imagers. There is a vast range of unexplored science accessible through astrometry, photometry, and spectroscopy of rocky planets, ice giants, and gas giants. In this whitepaper we will focus on one of the most ambitious science goals --- detecting for the first time habitable-zone rocky (<1.6 R_Earth) exoplanets in reflected light around nearby M-dwarfs △ Less

Submitted 10 May, 2019; originally announced May 2019.

Comments: 8 pages, 1 figure, Astro2020 Science White Paper

New Constraints From Dust Lines On The Surface Densities Of Protoplanetary Disks

Authors: Diana Powell, Ruth Murray-Clay, Laura M. Pérez, Hilke E. Schlichting, Mickey Rosenthal

Abstract: We present new determinations of disk surface density, independent of an assumed dust opacity, for a sample of 7 bright, diverse protoplanetary disks using measurements of disk dust lines. We develop a robust method for determining the location of dust lines by modeling disk interferometric visibilities at multiple wavelengths. The disks in our sample have newly derived masses that are 9-27% of th… ▽ More We present new determinations of disk surface density, independent of an assumed dust opacity, for a sample of 7 bright, diverse protoplanetary disks using measurements of disk dust lines. We develop a robust method for determining the location of dust lines by modeling disk interferometric visibilities at multiple wavelengths. The disks in our sample have newly derived masses that are 9-27% of their host stellar mass, substantially larger than the minimum mass solar nebula. All are stable to gravitational collapse except for one which approaches the limit of Toomre-Q stability. Our mass estimates are 2-15 times larger than estimates from integrated optically thin dust emission. We derive depleted dust-to-gas ratios with typical values of ~$10^{-3}$ in the outer disk. Using coagulation models we derive dust surface density profiles that are consistent with millimeter dust observations. In these models, the disks formed with an initial dust mass that is a factor of ~10 greater than is presently observed. Of the three disks in our sample with resolved CO line emission, the masses of HD 163296, AS 209, and TW Hya are roughly 3, 115, and 40 times more massive than estimates from CO respectively. This range indicates that CO depletion is not uniform across different disks and that dust is a more robust tracer of total disk mass. Our method of determining surface density using dust lines is robust even if particles form as aggregates and is useful even in the presence of dust substructure caused by pressure traps. The low Toomre-Q values observed in this sample indicate that at least some disks do not accrete efficiently. △ Less

Submitted 8 May, 2019; originally announced May 2019.

Comments: 29 pages, 12 figures, submitted to ApJ, revised following comments from referee

Wide-Orbit Exoplanet Demographics

Authors: David P. Bennett, Rachel Akeson, Yann Alibert, Jay Anderson, Etienne Bachelet, Jean-Phillipe Beaulieu, Andrea Bellini, Aparna Bhattacharya, Alan Boss, Valerio Bozza, Stephen Bryson, Derek Buzasi, Sebastiano Calchi Novati, Jessie Christiansen, Shawn D. Domagal-goldman, Michael Endl, Benjamin J. Fulton, Calen B. Henderson, B. Scott Gaudi, Samson A. Johnson, Naoki Koshimoto, Michael Meyer, Gijs D. Mulders, Susan Mullally, Ruth Murray-Clay , et al. (19 additional authors not shown)

Abstract: The Kepler, K2 and TESS transit surveys are revolutionizing our understanding of planets orbiting close to their host stars and our understanding of exoplanet systems in general, but there remains a gap in our understanding of wide-orbit planets. This gap in our understanding must be filled if we are to understand planet formation and how it affects exoplanet habitability. We summarize current and… ▽ More The Kepler, K2 and TESS transit surveys are revolutionizing our understanding of planets orbiting close to their host stars and our understanding of exoplanet systems in general, but there remains a gap in our understanding of wide-orbit planets. This gap in our understanding must be filled if we are to understand planet formation and how it affects exoplanet habitability. We summarize current and planned exoplanet detection programs using a variety of methods: microlensing (including WFIRST), radial velocities, Gaia astrometry, and direct imaging. Finally, we discuss the prospects for joint analyses using results from multiple methods and obstacles that could hinder such analyses. We endorse the findings and recommendations published in the 2018 National Academy report on Exoplanet Science Strategy. This white paper extends and complements the material presented therein. △ Less

Submitted 19 March, 2019; originally announced March 2019.

Comments: Astro2020 Science White paper

Imaging Giant Protoplanets with the ELTs

Authors: Steph Sallum, Vanessa P. Bailey, Rebecca A. Bernstein, Alan P. Boss, Brendan P. Bowler, Laird Close, Thayne Currie, Ruobing Dong, Catherine Espaillat, Michael P. Fitzgerald, Katherine B. Follette, Jonathan Fortney, Yasuhiro Hasegawa, Hannah Jang-Condell, Nemanja Jovanovic, Stephen R. Kane, Quinn Konopacky, Michael Liu, Julien Lozi, Jared Males, Dimitri Mawet, Benjamin Mazin, Max Millar-Blanchaer, Ruth Murray-Clay, Garreth Ruane , et al. (5 additional authors not shown)

Abstract: We have now accumulated a wealth of observations of the planet-formation environment and of mature planetary systems. These data allow us to test and refine theories of gas-giant planet formation by placing constraints on the conditions and timescale of this process. Yet a number of fundamental questions remain unanswered about how protoplanets accumulate material, their photospheric properties an… ▽ More We have now accumulated a wealth of observations of the planet-formation environment and of mature planetary systems. These data allow us to test and refine theories of gas-giant planet formation by placing constraints on the conditions and timescale of this process. Yet a number of fundamental questions remain unanswered about how protoplanets accumulate material, their photospheric properties and compositions, and how they interact with protoplanetary disks. While we have begun to detect protoplanet candidates during the last several years, we are presently only sensitive to the widest separation, highest mass / accretion rate cases. Current observing facilities lack the angular resolution and inner working angle to probe the few-AU orbital separations where giant planet formation is thought to be most efficient. They also lack the contrast to detect accretion rates that would form lower mass gas giants and ice giants. Instruments and telescopes coming online over the next decade will provide high contrast in the inner giant-planet-forming regions around young stars, allowing us to build a protoplanet census and to characterize planet formation in detail for the first time. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

Photoevaporative Flows From Exoplanet Atmospheres: A 3-D Radiative Hydrodynamic Parameter Study

Authors: Alex Debrecht, Jonathan Carroll-Nellenback, Adam Frank, John McCann, Ruth Murray-Clay, Eric G. Blackman

Abstract: The photoionization-driven evaporation of planetary atmospheres has emerged as a potentially fundamental process for planets on short period orbits. While 1-D studies have proven the effectiveness of stellar fluxes at altering the atmospheric mass and composition for sub-Jupiter mass planets, there remains much that is uncertain with regard to the larger-scale, multidimensional nature of such "pla… ▽ More The photoionization-driven evaporation of planetary atmospheres has emerged as a potentially fundamental process for planets on short period orbits. While 1-D studies have proven the effectiveness of stellar fluxes at altering the atmospheric mass and composition for sub-Jupiter mass planets, there remains much that is uncertain with regard to the larger-scale, multidimensional nature of such "planetary wind" flows. In this paper we use a new radiation-hydrodynamic platform to simulate atmospheric evaporative flows. Using the AstroBEAR AMR multiphysics code in a co-rotating frame centered on the planet, we model the transfer of ionizing photons into the atmosphere, the subsequent launch of the wind and the wind's large scale evolution subject to tidal and non-inertial forces. We run simulations for planets of 0.263 and 0.07 Jupiter masses and stellar fluxes of $2 \times 10^{13}$ and $2 \times 10^{14}$ photons/cm^2/s. Our results reveal new, potentially observable planetary wind flow patterns, including the development, in some cases, of an extended neutral tail lagging behind the planet in its orbit. △ Less

Submitted 22 November, 2018; originally announced November 2018.

Morphology of Hydrodynamic Winds: A Study of Planetary Winds in Stellar Environments

Authors: John R. McCann, Ruth A. Murray-Clay, Kaitlin Kratter, Mark R. Krumholz

Abstract: Bathed in intense ionizing radiation, close-in gaseous planets undergo hydrodynamic atmospheric escape, which ejects the upper extent of their atmospheres into the interplanetary medium. Ultraviolet detections of escaping gas around transiting planets corroborate such a framework. Exposed to the stellar environment, the outflow is shaped by its interaction with the stellar wind and by the planet's… ▽ More Bathed in intense ionizing radiation, close-in gaseous planets undergo hydrodynamic atmospheric escape, which ejects the upper extent of their atmospheres into the interplanetary medium. Ultraviolet detections of escaping gas around transiting planets corroborate such a framework. Exposed to the stellar environment, the outflow is shaped by its interaction with the stellar wind and by the planet's orbit. We model these effects using Athena to perform 3-D radiative-hydrodynamic simulations of tidally-locked hydrogen atmospheres receiving large amounts of ionizing extreme-ultraviolet flux in various stellar environments for the low-magnetic-field case. Through a step-by-step exploration of orbital and stellar wind effects on the planetary outflow, we find three structurally distinct stellar wind regimes: weak, intermediate, and strong. We perform synthetic Lyman-$α$ observations and find unique observational signatures for each regime. A weak stellar wind$\textrm{---}$which cannot confine the planetary outflow, leading to a torus of material around the star$\textrm{---}$has a pre-transit, red-shifted dayside arm and a slightly redward-skewed spectrum during transit. The intermediate regime truncates the dayside outflow at large distances from the planet and causes periodic disruptions of the outflow, producing observational signatures that mimic a double transit. The first of these dips is blue-shifted and precedes the optical transit. Finally, strong stellar winds completely confine the outflow into a cometary tail and accelerate the outflow outwards, producing large blue-shifted signals post-transit. Across all three regimes, large signals occur far outside of transit, offering motivation to continue ultraviolet observations outside of direct transit. △ Less

Submitted 22 November, 2018; originally announced November 2018.

Comments: 33 pages, 21 figures (7 of which have embedded movies viewable with Adobe Acrobat Pro), Submitted to ApJ

Restrictions on the Growth of Gas Giant Cores via Pebble Accretion

Authors: M. M. Rosenthal, R. A. Murray-Clay

Abstract: We apply an order-of-magnitude model of gas-assisted growth, known as pebble accretion, in a turbulent medium to suggest a reason why some systems form wide orbital separation gas giants while others do not. In contrast to traditional growth by ballistic collisions with planetesimals, growth by pebble accretion is not necessarily limited by doubling times at the highest core mass. Turbulence, in p… ▽ More We apply an order-of-magnitude model of gas-assisted growth, known as pebble accretion, in a turbulent medium to suggest a reason why some systems form wide orbital separation gas giants while others do not. In contrast to traditional growth by ballistic collisions with planetesimals, growth by pebble accretion is not necessarily limited by doubling times at the highest core mass. Turbulence, in particular, can cause growth to bottleneck at lower core masses. We demonstrate how a combination of growth by planetesimal and pebble accretion limits the maximum semi-major axis where gas giants can form. We find that, for fiducial disk parameters, strong turbulence ($α\gtrsim 10^{-2}$) restricts gas giant cores to form interior to $a \lesssim 40 \, \text{AU}$, while for weak turbulence gas giants can form out to $a \lesssim 70 \, \text{AU}$. The correspondence between $α$ and semi-major axis depends on the sizes of small bodies available for growth. This dependence on turbulence and small-body size distribution may explain the paucity of wide orbital separation gas giants. We also show that while lower levels of turbulence ($α\lesssim 10^{-4}$) can produce gas giants far out in the disk, we expect these gas giants to be low-mass ($M \lesssim \, 1 M_J$). These planets are not luminous enough to have been observed with the current generation of direct-imaging surveys, which could explain why wide orbital separation gas giants are currently observed only around A stars. △ Less

Submitted 17 August, 2018; originally announced August 2018.

Metallicity-Dependent Signatures in the Kepler Planets

Authors: James E. Owen, Ruth Murray-Clay

Abstract: Using data from the California-Kepler-Survey (CKS) we study trends in planetary properties with host star metallicity for close-in planets. By incorporating knowledge of the properties of the planetary radius gap identified by the CKS survey, we are able to investigate the properties of planetary cores and their gaseous envelopes separately. Our primary findings are that the solid core masses of p… ▽ More Using data from the California-Kepler-Survey (CKS) we study trends in planetary properties with host star metallicity for close-in planets. By incorporating knowledge of the properties of the planetary radius gap identified by the CKS survey, we are able to investigate the properties of planetary cores and their gaseous envelopes separately. Our primary findings are that the solid core masses of planets are higher around higher metallicity stars and that these more massive cores were able to accrete larger gas envelopes. Furthermore, investigating the recently reported result that planets with radii in the range (2-6Rearth) are more common at short periods around higher metallicity stars in detail, we find that the average host star metallicity of H/He atmosphere-hosting planets increases smoothly inside an orbital period of ~20 days. We interpret the location of the metallicity increase within the context of atmospheric photoevaporation: higher metallicity stars are likely to host planets with higher atmospheric metallicity, which increases the cooling in the photoevaporative outflow, lowering the mass-loss rates. Therefore, planets with higher metallicity atmospheres are able to resist photoevaporation at shorter orbital periods. Finally, we find evidence at 2.8 sigma that planets that do not host H/He atmospheres at long periods are more commonly found around lower metallicity stars. Such planets are difficult to explain by photoevaporative stripping of planets which originally accreted H/He atmospheres. Alternatively, this population of planets could be representative of planets that formed in a terrestrial-like fashion, after the gas disc dispersed. △ Less

Submitted 16 July, 2018; originally announced July 2018.

Comments: To appear in MNRAS

OSSOS. VII. 800+ trans-Neptunian objects - the complete data release

Authors: Michele T. Bannister, Brett J. Gladman, J. J. Kavelaars, Jean-Marc Petit, Kathryn Volk, Ying-Tung Chen, Mike Alexandersen, Stephen D. J. Gwyn, Megan E. Schwamb, Edward Ashton, Susan D. Benecchi, Nahuel Cabral, Rebekah I. Dawson, Audrey Delsanti, Wesley C. Fraser, Mikael Granvik, Sarah Greenstreet, Aurélie Guilbert-Lepoutre, Wing-Huen Ip, Marian Jakubik, R. Lynne Jones, Nathan A. Kaib, Pedro Lacerda, Christa Van Laerhoven, Samantha Lawler , et al. (11 additional authors not shown)

Abstract: The Outer Solar System Origins Survey (OSSOS), a wide-field imaging program in 2013-2017 with the Canada-France-Hawaii Telescope, surveyed 155 deg$^{2}$ of sky to depths of $m_r = 24.1$-25.2. We present 838 outer Solar System discoveries that are entirely free of ephemeris bias. This increases the inventory of trans-Neptunian objects (TNOs) with accurately known orbits by nearly 50%. Each minor pl… ▽ More The Outer Solar System Origins Survey (OSSOS), a wide-field imaging program in 2013-2017 with the Canada-France-Hawaii Telescope, surveyed 155 deg$^{2}$ of sky to depths of $m_r = 24.1$-25.2. We present 838 outer Solar System discoveries that are entirely free of ephemeris bias. This increases the inventory of trans-Neptunian objects (TNOs) with accurately known orbits by nearly 50%. Each minor planet has 20-60 Gaia/Pan-STARRS-calibrated astrometric measurements made over 2-5 oppositions, which allows accurate classification of their orbits within the trans-Neptunian dynamical populations. The populations orbiting in mean-motion resonance with Neptune are key to understanding Neptune's early migration. Our 313 resonant TNOs, including 132 plutinos, triple the available characterized sample and include new occupancy of distant resonances out to semi-major axis $a \sim 130$ au. OSSOS doubles the known population of the non-resonant Kuiper belt, providing 436 TNOs in this region, all with exceptionally high-quality orbits of $a$ uncertainty $σ_{a} \leq 0.1\%$; they show the belt exists from $a \gtrsim 37$ au, with a lower perihelion bound of $35$ au. We confirm the presence of a concentrated low-inclination $a\simeq 44$ au "kernel" population and a dynamically cold population extending beyond the 2:1 resonance. We finely quantify the survey's observational biases. Our survey simulator provides a straightforward way to impose these biases on models of the trans-Neptunian orbit distributions, allowing statistical comparison to the discoveries. The OSSOS TNOs, unprecedented in their orbital precision for the size of the sample, are ideal for testing concepts of the history of giant planet migration in the Solar System. △ Less

Submitted 29 May, 2018; originally announced May 2018.

Comments: Invited paper, special issue Data: Insights and Challenges in a Time of Abundance. Data tables and example survey simulator are in the supplementary materials (see arXiv source under Downloads > Other formats)

Journal ref: The Astrophysical Journal Supplement Series, 236(1):18, (19 pp), 2018

Trans-Neptunian Objects Transiently Stuck in Neptune's Mean Motion Resonances: Numerical simulations of the current population

Authors: Tze Yeung Mathew Yu, Ruth Murray-Clay, Kathryn Volk

Abstract: A substantial fraction of our solar system's trans-Neptunian objects (TNOs) are in mean motion resonance with Neptune. Many of these objects were likely caught into resonances by planetary migration---either smooth or stochastic---approximately 4 Gyr ago. Some, however, gravitationally scattered off of Neptune and became transiently stuck in more recent events. Here, we use numerical simulations t… ▽ More A substantial fraction of our solar system's trans-Neptunian objects (TNOs) are in mean motion resonance with Neptune. Many of these objects were likely caught into resonances by planetary migration---either smooth or stochastic---approximately 4 Gyr ago. Some, however, gravitationally scattered off of Neptune and became transiently stuck in more recent events. Here, we use numerical simulations to predict the number of transiently-stuck objects, captured from the current actively scattering population, that occupy 111 resonances at semimajor axes $a=$30--100 au. Our source population is an observationally constrained model of the currently-scattering TNOs. We predict that, integrated across all resonances at these distances, the current transient sticking population comprises 40\% of total transiently-stuck+scattering TNOs, suggesting that these objects should be treated as a single population. We compute the relative distribution of transiently-stuck objects across all $p$:$q$ resonances with $1/6 \le q/p < 1$, $p<40$, and $q<20$, providing predictions for the population of transient objects with $H_r < 8.66$ in each resonance. We find that the relative populations are approximately proportional to each resonance's libration period and confirm that the importance of transient sticking increases with semimajor axis in the studied range. We calculate the expected distribution of libration amplitudes for stuck objects and demonstrate that observational constraints indicate that both the total number and the amplitude-distribution of 5:2 resonant TNOs are inconsistent with a population dominated by transient sticking from the current scattering disk. The 5:2 resonance hence poses a challenge for leading theories of Kuiper belt sculpting. △ Less

Submitted 21 May, 2018; originally announced May 2018.

Gas-Assisted Growth of Protoplanets in a Turbulent Medium

Authors: M. M. Rosenthal, R. A. Murray-Clay, H. B. Perets, N. Wolansky

Abstract: Pebble accretion is a promising process for decreasing growth timescales of planetary cores, allowing gas giants to form at wide orbital separations. However, nebular turbulence can reduce the efficiency of this gas-assisted growth. We present an order of magnitude model of pebble accretion, which calculates the impact of turbulence on the average velocity of small bodies, the radius for binary ca… ▽ More Pebble accretion is a promising process for decreasing growth timescales of planetary cores, allowing gas giants to form at wide orbital separations. However, nebular turbulence can reduce the efficiency of this gas-assisted growth. We present an order of magnitude model of pebble accretion, which calculates the impact of turbulence on the average velocity of small bodies, the radius for binary capture, and the sizes of the small bodies that can be accreted. We also include the effect of turbulence on the particle scale height, which has been studied in previous works. We find that turbulence does not prevent rapid growth in the high-mass regime: the last doubling time to the critical mass to trigger runaway gas accretion ($M \sim 10 \, M_\oplus$) is well within the disk lifetime even for strong ($α\gtrsim 10^{-2}$) turbulence. We find that while the growth timescale is quite sensitive to the local properties of the protoplanetary disk, over large regimes of parameter space large cores grow in less than the disk lifetime if appropriately-sized small bodies are present. Instead, the effects of turbulence are most pronounced for low planetary masses. For strong turbulence the growth timescale is longer than the gas disk lifetime until the core reaches masses $\gtrsim 10^{-2}-10^{-1} M_\oplus$. A "Flow Isolation Mass", at which binary capture ceases, emerges naturally from our model framework. We comment that the dependence of this mass on orbital separation is similar to the semi-major axis distribution of solar-system cores. △ Less

Submitted 17 May, 2018; originally announced May 2018.

Comments: 40 pages, 20 figures

OSSOS IX: two objects in Neptune's 9:1 resonance -- implications for resonance sticking in the scattering population

Authors: Kathryn Volk, Ruth A. Murray-Clay, Brett J. Gladman, Samantha M. Lawler, Tze Yeung Mathew Yu, Mike Alexandersen, Michele T. Bannister, Ying-Tung Chen, Rebekah I. Dawson, Sarah Greenstreet, Stephen D. J. Gwyn, J. J. Kavelaars, Hsing Wen Lin, Patryk Sofia Lykawka, Jean-Marc Petit

Abstract: We discuss the detection in the Outer Solar System Origins Survey (OSSOS) of two objects in Neptune's distant 9:1 mean motion resonance at semimajor axis $a\approx~130$~au. Both objects are securely resonant on 10~Myr timescales, with one securely in the 9:1 resonance's leading asymmetric libration island and the other in either the symmetric or trailing asymmetric island. These objects are the la… ▽ More We discuss the detection in the Outer Solar System Origins Survey (OSSOS) of two objects in Neptune's distant 9:1 mean motion resonance at semimajor axis $a\approx~130$~au. Both objects are securely resonant on 10~Myr timescales, with one securely in the 9:1 resonance's leading asymmetric libration island and the other in either the symmetric or trailing asymmetric island. These objects are the largest semimajor axis objects with secure resonant classifications, and their detection in a carefully characterized survey allows for the first robust resonance population estimate beyond 100~au. The detection of these objects implies a 9:1 resonance population of $1.1\times10^4$ objects with $H_r<8.66$ ($D~\gtrsim~100$~km) on similar orbits (95\% confidence range of $\sim0.4-3\times10^4$). Integrations over 4~Gyr of an ensemble of clones spanning these objects' orbit fit uncertainties reveal that they both have median resonance occupation timescales of $\sim1$~Gyr. These timescales are consistent with the hypothesis that these objects originate in the scattering population but became transiently stuck to Neptune's 9:1 resonance within the last $\sim1$~Gyr of solar system evolution. Based on simulations of a model of the current scattering population, we estimate the expected resonance sticking population in the 9:1 resonance to be 1000-4500 objects with $H_r<8.66$; this is marginally consistent with the OSSOS 9:1 population estimate. We conclude that resonance sticking is a plausible explanation for the observed 9:1 population, but we also discuss the possibility of a primordial 9:1 population, which would have interesting implications for the Kuiper belt's dynamical history. △ Less

Submitted 24 May, 2018; v1 submitted 15 February, 2018; originally announced February 2018.

Comments: accepted for publication in AJ

Journal ref: The Astronomical Journal, 155:260 (9pp), 2018 June

Pebble Accretion in Turbulent Protoplanetary Disks

Authors: Ziyan Xu, Xue-Ning Bai, Ruth Murray-Clay

Abstract: It has been realized in recent years that the accretion of pebble-sized dust particles onto planetary cores is an important mode of core growth, which enables the formation of giant planets at large distances and assists planet formation in general. The pebble accretion theory is built upon the orbit theory of dust particles in a laminar protoplanetary disk (PPD). For sufficiently large core mass… ▽ More It has been realized in recent years that the accretion of pebble-sized dust particles onto planetary cores is an important mode of core growth, which enables the formation of giant planets at large distances and assists planet formation in general. The pebble accretion theory is built upon the orbit theory of dust particles in a laminar protoplanetary disk (PPD). For sufficiently large core mass (in the "Hill regime"), essentially all particles of appropriate sizes entering the Hill sphere can be captured. However, the outer regions of PPDs are expected to be weakly turbulent due to the magnetorotational instability (MRI), where turbulent stirring of particle orbits may affect the efficiency of pebble accretion. We conduct shearing-box simulations of pebble accretion with different levels of MRI turbulence (strongly turbulent assuming ideal magnetohydrodynamics, weakly turbulent in the presence of ambipolar diffusion, and laminar) and different core masses to test the efficiency of pebble accretion at a microphysical level. We find that accretion remains efficient for marginally coupled particles (dimensionless stopping time tau_s ~ 0.1 - 1) even in the presence of strong MRI turbulence. Though more dust particles are brought toward the core by the turbulence, this effect is largely canceled by a reduction in accretion probability. As a result, the overall effect of turbulence on the accretion rate is mainly reflected in the changes in the thickness of the dust layer. On the other hand, we find that the efficiency of pebble accretion for strongly coupled particles (down to tau_s ~ 0.01) can be modestly reduced by strong turbulence for low-mass cores. △ Less

Submitted 11 September, 2017; originally announced September 2017.

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

The Formation of Uranus and Neptune: Fine Tuning in Core Accretion

Authors: Renata Frelikh, Ruth A. Murray-Clay

Abstract: Uranus and Neptune are ice giants with $\sim$ 15% atmospheres by mass, placing them in an intermediate category between rocky planets and gas giants. These atmospheres are too massive to have been primarily outgassed, yet they never underwent runaway gas accretion. The ice giants never reached critical core mass ($M_\text{crit}$) in a full gas disk, yet their cores are $\gtrsim M_\text{crit}$, sug… ▽ More Uranus and Neptune are ice giants with $\sim$ 15% atmospheres by mass, placing them in an intermediate category between rocky planets and gas giants. These atmospheres are too massive to have been primarily outgassed, yet they never underwent runaway gas accretion. The ice giants never reached critical core mass ($M_\text{crit}$) in a full gas disk, yet their cores are $\gtrsim M_\text{crit}$, suggesting that their envelopes were mainly accreted at the end of the disk lifetime. Pebble accretion calls into question traditional slow atmospheric growth during this phase. We show that the full-sized ice giants predominantly accreted gas from a disk depleted by at least a factor of $\sim 100$. Such a disk dissipates in $ \lesssim 10^5$ years. Why would both cores stay sub-critical for the entire $\sim$ Myr disk lifetime, only to reach $M_\text{crit}$ in the final $10^5$ years? This is fine tuned. Ice giants in the outer disk have atmospheric mass fractions comparable to the disk gas-to-solid ratio during the bulk of their gas accretion. This point in disk evolution coincides with a dynamical upheaval: the gas loses its ability to efficiently damp the cores' random velocities, allowing them to be gravitationally excited by Jupiter and Saturn. We suggest that the ice giants' cores began growing on closer-in orbits (staying sub-critical), and migrated out during this dynamical instability. There, their orbits circularized after accreting much of their mass in solids. Finally, they accreted their envelopes from a depleted nebula, where the sparseness of feeding zone gas prevented runaway. △ Less

Submitted 2 August, 2017; originally announced August 2017.

Comments: Accepted for publication in AJ, 10 pages, 3 figures

Evidence that the Directly-Imaged Planet HD 131399 Ab is a Background Star

Authors: Eric L. Nielsen, Robert J. De Rosa, Julien Rameau, Jason J. Wang, Thomas M. Esposito, Maxwell A. Millar-Blanchaer, Christian Marois, Arthur Vigan, S. Mark Ammons, Etienne Artigau, Vanessa P. Bailey, Sarah Blunt, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, René Doyon, Gaspard Duchêne, Daniel Fabrycky, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, James R. Graham, Alexandra Z. Greenbaum, Pascale Hibon , et al. (33 additional authors not shown)

Abstract: We present evidence that the recently discovered, directly-imaged planet HD 131399 Ab is a background star with non-zero proper motion. From new JHK1L' photometry and spectroscopy obtained with the Gemini Planet Imager, VLT/SPHERE, and Keck/NIRC2, and a reanalysis of the discovery data obtained with VLT/SPHERE, we derive colors, spectra, and astrometry for HD 131399 Ab. The broader wavelength cove… ▽ More We present evidence that the recently discovered, directly-imaged planet HD 131399 Ab is a background star with non-zero proper motion. From new JHK1L' photometry and spectroscopy obtained with the Gemini Planet Imager, VLT/SPHERE, and Keck/NIRC2, and a reanalysis of the discovery data obtained with VLT/SPHERE, we derive colors, spectra, and astrometry for HD 131399 Ab. The broader wavelength coverage and higher data quality allow us to re-investigate its status. Its near-infrared spectral energy distribution excludes spectral types later than L0 and is consistent with a K or M dwarf, which are the most likely candidates for a background object in this direction at the apparent magnitude observed. If it were a physically associated object, the projected velocity of HD 131399 Ab would exceed escape velocity given the mass and distance to HD 131399 A. We show that HD 131399 Ab is also not following the expected track for a stationary background star at infinite distance. Solving for the proper motion and parallax required to explain the relative motion of HD 131399 Ab, we find a proper motion of 12.3 mas/yr. When compared to predicted background objects drawn from a galactic model, we find this proper motion to be high, but consistent with the top 4% fastest-moving background stars. From our analysis we conclude that HD 131399 Ab is a background K or M dwarf. △ Less

Submitted 8 November, 2017; v1 submitted 18 May, 2017; originally announced May 2017.

Comments: 35 pages, 20 figures, Published in AJ

Journal ref: Nielsen et al. 2017 AJ, 154, 6

Using Ice and Dust Lines to Constrain the Surface Densities of Protoplanetary Disks

Authors: Diana Powell, Ruth Murray-Clay, Hilke E. Schlichting

Abstract: We present a novel method for determining the surface density of protoplanetary disks through consideration of disk 'dust lines' which indicate the observed disk radial scale at different observational wavelengths. This method relies on the assumption that the processes of particle growth and drift control the radial scale of the disk at late stages of disk evolution such that the lifetime of the… ▽ More We present a novel method for determining the surface density of protoplanetary disks through consideration of disk 'dust lines' which indicate the observed disk radial scale at different observational wavelengths. This method relies on the assumption that the processes of particle growth and drift control the radial scale of the disk at late stages of disk evolution such that the lifetime of the disk is equal to both the drift timescale and growth timescale of the maximum particle size at a given dust line. We provide an initial proof of concept of our model through an application to the disk TW Hya and are able to estimate the disk dust-to-gas ratio, CO abundance, and accretion rate in addition to the total disk surface density. We find that our derived surface density profile and dust-to-gas ratio are consistent with the lower limits found through measurements of HD gas. The CO ice line also depends on surface density through grain adsorption rates and drift and we find that our theoretical CO ice line estimates have clear observational analogues. We further apply our model to a large parameter space of theoretical disks and find three observational diagnostics that may be used to test its validity. First we predict that the dust lines of disks other than TW Hya will be consistent with the normalized CO surface density profile shape for those disks. Second, surface density profiles that we derive from disk ice lines should match those derived from disk dust lines. Finally, we predict that disk dust and ice lines will scale oppositely, as a function of surface density, across a large sample of disks. △ Less

Submitted 16 August, 2019; v1 submitted 15 April, 2017; originally announced April 2017.

Comments: Published in ApJ, May 2017

Metal enrichment leads to low atmospheric C/O ratios in transiting giant exoplanets

Authors: Néstor Espinoza, Jonathan Fortney, Yamila Miguel, Daniel Thorngren, Ruth Murray-Clay

Abstract: We predict the carbon-to-oxygen (C/O) ratios in the hydrogen-helium envelope and atmospheres of a sample of nearly 50 relatively cool ($T_{\mathrm eq}<$ 1000 K) transiting gas giant planets. The method involves planetary envelope metallicity estimates that use the structure models of Thorngren et al. (2016) and the disk and planetary accretion model of Öberg et al. (2011). We find that nearly all… ▽ More We predict the carbon-to-oxygen (C/O) ratios in the hydrogen-helium envelope and atmospheres of a sample of nearly 50 relatively cool ($T_{\mathrm eq}<$ 1000 K) transiting gas giant planets. The method involves planetary envelope metallicity estimates that use the structure models of Thorngren et al. (2016) and the disk and planetary accretion model of Öberg et al. (2011). We find that nearly all of these planets are strongly metal-enriched which, coupled with the fact that solid material is the main deliverer of metals in the protoplanetary disk, implies that the substellar C/O ratios of their accreted solid material dominate compared to the enhanced C/O ratio of their accreted gaseous component. We predict that these planets will have atmospheres that are typically reduced in their C/O compared to parent star values independent of the assessed formation locations, with C/O $<1$ a nearly universal outcome within the framework of the model. We expect water vapor absorption features to be ubiquitous in the atmospheres of these planets, and by extension, other gas giants. △ Less

Submitted 10 March, 2017; v1 submitted 25 November, 2016; originally announced November 2016.

Comments: 7 pages, 5 figures. Accepted to ApJL. Changes now include more citations in the discussion and tests of our model to Jupiter and Saturn

OSSOS: IV. Discovery of a dwarf planet candidate in the 9:2 resonance with Neptune

Authors: Michele T. Bannister, Mike Alexandersen, Susan D. Benecchi, Ying-Tung Chen, Audrey Delsanti, Wesley C. Fraser, Brett J. Gladman, Mikael Granvik, Will M. Grundy, Aurelie Guilbert-Lepoutre, Stephen D. J. Gwyn, Wing-Huen Ip, Marian Jakubik, R. Lynne Jones, Nathan Kaib, J. J. Kavelaars, Pedro Lacerda, Samantha Lawler, Matthew J. Lehner, Hsing Wen Lin, Patryk Sofia Lykawka, Michael Marsset, Ruth Murray-Clay, Keith S. Noll, Alex Parker , et al. (10 additional authors not shown)

Abstract: We report the discovery and orbit of a new dwarf planet candidate, 2015 RR$_{245}$, by the Outer Solar System Origins Survey (OSSOS). 2015 RR$_{245}$'s orbit is eccentric ($e=0.586$), with a semi-major axis near 82 au, yielding a perihelion distance of 34 au. 2015 RR$_{245}$ has $g-r = 0.59 \pm 0.11$ and absolute magnitude $H_{r} = 3.6 \pm 0.1$; for an assumed albedo of $p_V = 12$% the object has… ▽ More We report the discovery and orbit of a new dwarf planet candidate, 2015 RR$_{245}$, by the Outer Solar System Origins Survey (OSSOS). 2015 RR$_{245}$'s orbit is eccentric ($e=0.586$), with a semi-major axis near 82 au, yielding a perihelion distance of 34 au. 2015 RR$_{245}$ has $g-r = 0.59 \pm 0.11$ and absolute magnitude $H_{r} = 3.6 \pm 0.1$; for an assumed albedo of $p_V = 12$% the object has a diameter of $\sim670$ km. Based on astrometric measurements from OSSOS and Pan-STARRS1, we find that 2015 RR$_{245}$ is securely trapped on ten-Myr timescales in the 9:2 mean-motion resonance with Neptune. It is the first TNO identified in this resonance. On hundred-Myr timescales, particles in 2015 RR$_{245}$-like orbits depart and sometimes return to the resonance, indicating that 2015 RR$_{245}$ likely forms part of the long-lived metastable population of distant TNOs that drift between resonance sticking and actively scattering via gravitational encounters with Neptune. The discovery of a 9:2 TNO stresses the role of resonances in the long-term evolution of objects in the scattering disk, and reinforces the view that distant resonances are heavily populated in the current Solar System. This object further motivates detailed modelling of the transient sticking population. △ Less

Submitted 5 October, 2016; v1 submitted 23 July, 2016; originally announced July 2016.

Comments: Accepted to AJ

OSSOS III - Resonant Trans-Neptunian Populations: Constraints from the first quarter of the Outer Solar System Origins Survey

Authors: Kathryn Volk, Ruth Murray-Clay, Brett Gladman, Samantha Lawler, Michele T. Bannister, J. J. Kavelaars, Jean-Marc Petit, Stephen Gwyn, Mike Alexandersen, Ying-Tung Chen, Patryk Sofia Lykawka, Wing Ip, Hsing Wen Lin

Abstract: The first two observational sky "blocks" of the Outer Solar System Origins Survey (OSSOS) have significantly increased the number of well-characterized observed trans-Neptunian objects (TNOs) in Neptune's mean motion resonances. We describe the 31 securely resonant TNOs detected by OSSOS so far, and we use them to independently verify the resonant population models from the Canada-France Ecliptic… ▽ More The first two observational sky "blocks" of the Outer Solar System Origins Survey (OSSOS) have significantly increased the number of well-characterized observed trans-Neptunian objects (TNOs) in Neptune's mean motion resonances. We describe the 31 securely resonant TNOs detected by OSSOS so far, and we use them to independently verify the resonant population models from the Canada-France Ecliptic Plane Survey (CFEPS; Gladman et al. 2012), with which we find broad agreement. We confirm that the 5:2 resonance is more populated than models of the outer Solar System's dynamical history predict; our minimum population estimate shows that the high eccentricity (e>0.35) portion of the resonance is at least as populous as the 2:1 and possibly as populated as the 3:2 resonance. One OSSOS block was well-suited to detecting objects trapped at low libration amplitudes in Neptune's 3:2 resonance, a population of interest in testing the origins of resonant TNOs. We detected three 3:2 objects with libration amplitudes below the cutoff modeled by CFEPS; OSSOS thus offers new constraints on this distribution. The OSSOS detections confirm that the 2:1 resonance has a dynamically colder inclination distribution than either the 3:2 or 5:2 resonances. Using the combined OSSOS and CFEPS 2:1 detections, we constrain the fraction of 2:1 objects in the symmetric mode of libration to be 0.2-0.85; we also constrain the fraction of leading vs. trailing asymmetric librators, which has been theoretically predicted to vary depending on Neptune's migration history, to be 0.05-0.8. Future OSSOS blocks will improve these constraints. △ Less

Submitted 27 April, 2016; originally announced April 2016.

Comments: Accepted for publication in AJ

The Mass-Metallicity Relation for Giant Planets

Authors: Daniel P. Thorngren, Jonathan J. Fortney, Ruth A. Murray-Clay, Eric D. Lopez

Abstract: Exoplanet discoveries of recent years have provided a great deal of new data for studying the bulk compositions of giant planets. Here we identify 47 transiting giant planets ($20 M_\oplus < M < 20 M_{\mathrm{J}}$) whose stellar insolation is low enough ($F_* < 2\times10^8\; \text{erg}\; \text{s}^{-1}\; \text{cm}^{-2}$, or roughly $T_\text{eff} < 1000$) that they are not affected by the hot Jupite… ▽ More Exoplanet discoveries of recent years have provided a great deal of new data for studying the bulk compositions of giant planets. Here we identify 47 transiting giant planets ($20 M_\oplus < M < 20 M_{\mathrm{J}}$) whose stellar insolation is low enough ($F_* < 2\times10^8\; \text{erg}\; \text{s}^{-1}\; \text{cm}^{-2}$, or roughly $T_\text{eff} < 1000$) that they are not affected by the hot Jupiter radius inflation mechanism(s). We compute a set of new thermal and structural evolution models and use these models in comparison with properties of the 47 transiting planets (mass, radius, age) to determine their heavy element masses. A clear correlation emerges between the planetary heavy element mass $M_z$ and the total planet mass, approximately of the form $M_z \propto \sqrt{M}$. This finding is consistent with the core accretion model of planet formation. We also study how stellar metallicity [Fe/H] affects planetary metal-enrichment and find a weaker correlation than has been previously reported from studies with smaller sample sizes. We confirm a strong relationship between the planetary metal-enrichment relative to the parent star $Z_{\rm planet}/Z_{\rm star}$ and the planetary mass, but see no relation in $Z_{\rm planet}/Z_{\rm star}$ with planet orbital properties or stellar mass. The large heavy element masses of many planets ($>50$ $M_{\oplus}$) suggest significant amounts of heavy elements in H/He envelopes, rather than cores, such that metal-enriched giant planet atmospheres should be the rule. We also discuss a model of core-accretion planet formation in a one-dimensional disk and show that it agrees well with our derived relation between mass and $Z_{\rm planet}/Z_{\rm star}$. △ Less

Submitted 30 August, 2016; v1 submitted 24 November, 2015; originally announced November 2015.

Comments: Accepted to The Astrophysical Journal. This revision adds a substantial amount of discussion; the results are the same

C/O and Snowline Locations in Protoplanetary Disks: The Effect of Radial Drift and Viscous Gas Accretion

Authors: Ana-Maria A. Piso, Karin I. Oberg, Tilman Birnstiel, Ruth A. Murray-Clay

Abstract: The C/O ratio is a defining feature of both gas giant atmospheric and protoplanetary disk chemistry. In disks, the C/O ratio is regulated by the presence of snowlines of major volatiles at different distances from the central star. We explore the effect of radial drift of solids and viscous gas accretion onto the central star on the snowline locations of the main C and O carriers in a protoplaneta… ▽ More The C/O ratio is a defining feature of both gas giant atmospheric and protoplanetary disk chemistry. In disks, the C/O ratio is regulated by the presence of snowlines of major volatiles at different distances from the central star. We explore the effect of radial drift of solids and viscous gas accretion onto the central star on the snowline locations of the main C and O carriers in a protoplanetary disk, H2O, CO2 and CO, and their consequences for the C/O ratio in gas and dust throughout the disk. We determine the snowline locations for a range of fixed initial particle sizes and disk types. For our fiducial disk model, we find that grains with sizes ~0.5 cm < s < 7 m for an irradiated disk, and ~0.001 cm < s < 7 m for an evolving and viscous disk, desorb at a size-dependent location in the disk, which is independent of the particle's initial position. The snowline radius decreases for larger particles, up to sizes of ~7 m. Compared to a static disk, we find that radial drift and gas accretion in a viscous disk move the H2O snowline inwards by up to 40%, the CO2 snowline by up to 60%, and the CO snowline by up to 50%. We thus determine an inner limit on the snowline locations when radial drift and gas accretion are accounted for. △ Less

Submitted 17 November, 2015; originally announced November 2015.

Comments: 16 pages, 8 figures. Accepted to ApJ

The Outer Solar System Origins Survey: I. Design and First-Quarter Discoveries

Authors: Michele T. Bannister, J. J. Kavelaars, Jean-Marc Petit, Brett J. Gladman, Stephen D. J. Gwyn, Ying-Tung Chen, Kathryn Volk, Mike Alexandersen, Susan Benecchi, Audrey Delsanti, Wesley Fraser, Mikael Granvik, Will M. Grundy, Aurelie Guilbert-Lepoutre, Daniel Hestroffer, Wing-Huen Ip, Marian Jakubik, Lynne Jones, Nathan Kaib, Catherine F. Kavelaars, Pedro Lacerda, Samantha Lawler, Matthew J. Lehner, Hsing Wen Lin, Tim Lister , et al. (14 additional authors not shown)

Abstract: We report the discovery, tracking and detection circumstances for 85 trans-Neptunian objects (TNOs) from the first 42 deg$^{2}$ of the Outer Solar System Origins Survey (OSSOS). This ongoing $r$-band Solar System survey uses the 0.9 deg$^{2}$ field-of-view MegaPrime camera on the 3.6 m Canada-France-Hawaii Telescope. Our orbital elements for these TNOs are precise to a fractional semi-major axis u… ▽ More We report the discovery, tracking and detection circumstances for 85 trans-Neptunian objects (TNOs) from the first 42 deg$^{2}$ of the Outer Solar System Origins Survey (OSSOS). This ongoing $r$-band Solar System survey uses the 0.9 deg$^{2}$ field-of-view MegaPrime camera on the 3.6 m Canada-France-Hawaii Telescope. Our orbital elements for these TNOs are precise to a fractional semi-major axis uncertainty $<0.1\%$. We achieve this precision in just two oppositions, as compared to the normal 3-5 oppositions, via a dense observing cadence and innovative astrometric technique. These discoveries are free of ephemeris bias, a first for large trans-Neptunian surveys. We also provide the necessary information to enable models of TNO orbital distributions to be tested against our TNO sample. We confirm the existence of a cold "kernel" of objects within the main cold classical Kuiper belt, and infer the existence of an extension of the "stirred" cold classical Kuiper belt to at least several AU beyond the 2:1 mean motion resonance with Neptune. We find that the population model of Petit et al. (2011) remains a plausible representation of the Kuiper belt. The full survey, to be completed in 2017, will provide an exquisitely characterized sample of important resonant TNO populations, ideal for testing models of giant planet migration during the early history of the Solar System. △ Less

Submitted 10 May, 2016; v1 submitted 9 November, 2015; originally announced November 2015.

Comments: Accepted to AJ, 27 April 2016. 59 pp

Direct imaging of an asymmetric debris disk in the HD 106906 planetary system

Authors: Paul G. Kalas, Abhijith Rajan, Jason J. Wang, Maxwell A. Millar-Blanchaer, Gaspard Duchene, Christine Chen, Michael P. Fitzgerald, Ruobing Dong, James R. Graham, Jennifer Patience, Bruce Macintosh, Ruth Murray-Clay, Brenda Matthews, Julien Rameau, Christian Marois, Jeffrey Chilcote, Robert J. De Rosa, René Doyon, Zachary H. Draper, Samantha Lawler, S. Mark Ammons, Pauline Arriaga, Joanna Bulger, Tara Cotten, Katherine B. Follette , et al. (31 additional authors not shown)

Abstract: We present the first scattered light detections of the HD 106906 debris disk using Gemini/GPI in the infrared and HST/ACS in the optical. HD 106906 is a 13 Myr old F5V star in the Sco-Cen association, with a previously detected planet-mass candidate HD 106906b projected 650 AU from the host star. Our observations reveal a near edge-on debris disk that has a central cleared region with radius… ▽ More We present the first scattered light detections of the HD 106906 debris disk using Gemini/GPI in the infrared and HST/ACS in the optical. HD 106906 is a 13 Myr old F5V star in the Sco-Cen association, with a previously detected planet-mass candidate HD 106906b projected 650 AU from the host star. Our observations reveal a near edge-on debris disk that has a central cleared region with radius $\sim$50 AU, and an outer extent $>$500 AU. The HST data show the outer regions are highly asymmetric, resembling the ''needle'' morphology seen for the HD 15115 debris disk. The planet candidate is oriented $\sim$21$°$ away from the position angle of the primary's debris disk, strongly suggesting non-coplanarity with the system. We hypothesize that HD 106906b could be dynamically involved in the perturbation of the primary's disk, and investigate whether or not there is evidence for a circumplanetary dust disk or cloud that is either primordial or captured from the primary. We show that both the existing optical properties and near-infrared colors of HD 106906b are weakly consistent with this possibility, motivating future work to test for the observational signatures of dust surrounding the planet. △ Less

Submitted 9 October, 2015; originally announced October 2015.

Comments: Accepted for publication in the Astrophysical Journal; 15 pages, 7 figures

Discovery and spectroscopy of the young Jovian planet 51 Eri b with the Gemini Planet Imager

Authors: B. Macintosh, J. R. Graham, T. Barman, R. J. De Rosa, Q. Konopacky, M. S. Marley, C. Marois, E. L. Nielsen, L. Pueyo, A. Rajan, J. Rameau, D. Saumon, J. J. Wang, J. Patience, M. Ammons, P. Arriaga, E. Artigau, S. Beckwith, J. Brewster, S. Bruzzone, J. Bulger, B. Burningham, A. S. Burrows, C. Chen, E. Chiang , et al. (63 additional authors not shown)

Abstract: Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane… ▽ More Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane and water vapor absorption. Modeling of the spectra and photometry yields a luminosity of L/LS=1.6-4.0 x 10-6 and an effective temperature of 600-750 K. For this age and luminosity, "hot-start" formation models indicate a mass twice that of Jupiter. This planet also has a sufficiently low luminosity to be consistent with the "cold- start" core accretion process that may have formed Jupiter. △ Less

Submitted 9 October, 2015; v1 submitted 12 August, 2015; originally announced August 2015.

Comments: 29 pages, 3 figures, 2 tables, and Supplementary Materials. published in Science Express on Aug 13 2015. List of authors and the magnitudes of the star were correted

Journal ref: Macintosh et al., Science 350 (6256): 64-67, 2015

Simulated Photoevaporative Mass Loss from Hot Jupiters in 3D

Authors: Anjali Tripathi, Kaitlin M. Kratter, Ruth A. Murray-Clay, Mark R. Krumholz

Abstract: Ionizing stellar photons heat the upper regions of planetary atmospheres, driving atmospheric mass loss. Gas escaping from several hot, hydrogen-rich planets has been detected using UV and X-ray transmission spectroscopy. Because these planets are tidally locked, and thus asymmetrically irradiated, escaping gas is unlikely to be spherically symmetric. In this paper, we focus on the effects of asym… ▽ More Ionizing stellar photons heat the upper regions of planetary atmospheres, driving atmospheric mass loss. Gas escaping from several hot, hydrogen-rich planets has been detected using UV and X-ray transmission spectroscopy. Because these planets are tidally locked, and thus asymmetrically irradiated, escaping gas is unlikely to be spherically symmetric. In this paper, we focus on the effects of asymmetric heating on local outflow structure. We use the Athena code for hydrodynamics to produce 3D simulations of hot Jupiter mass loss that jointly model wind launching and stellar heating via photoionization. Our fiducial planet is an inflated, hot Jupiter with radius $R_p=2.14 R_{\rm Jup}$ and mass $M_p = 0.53 M_{\rm Jup}$. We irradiate the initially neutral, atomic hydrogen atmosphere with 13.6 eV photons and compute the outflow's ionization structure. There are clear asymmetries in the atmospheric outflow, including a neutral shadow on the planet's nightside. Given an incident ionizing UV flux comparable to that of the Sun, we find a steady-state mass loss rate of ~$2\times10^{10}$ g s$^{-1}$. The total mass loss rate and the outflow substructure along the substellar ray show good agreement with earlier 1D models, for two different fluxes. Our 3D data cube can be used to generate the outflow's extinction spectrum during transit. As a proof of concept, we find absorption of stellar Lyman-alpha at Doppler-shifted velocities of up to $\pm 50$ km s$^{-1}$. Our work provides a starting point for further 3D models that can be used to predict observable signatures of hot Jupiter mass loss. △ Less

Submitted 3 August, 2015; v1 submitted 22 June, 2015; originally announced June 2015.

Comments: 16 pages, 15 figures, accepted for publication in ApJ. Replacement includes revised citations. Code is available for download and use from http://github.com/tripathi/Atmospheric-Athena

Journal ref: Ap.J. 808 (2015) 173

Minimum Core Masses for Giant Planet Formation With Realistic Equations of State and Opacities

Authors: Ana-Maria A. Piso, Andrew N. Youdin, Ruth A. Murray-Clay

Abstract: Giant planet formation by core accretion requires a core that is sufficiently massive to trigger runaway gas accretion in less that the typical lifetime of protoplanetary disks. We explore how the minimum required core mass, M_crit, depends on a non-ideal equation of state and on opacity changes due to grain growth, across a range of stellocentric distances from 5-100 AU. This minimum M_crit appli… ▽ More Giant planet formation by core accretion requires a core that is sufficiently massive to trigger runaway gas accretion in less that the typical lifetime of protoplanetary disks. We explore how the minimum required core mass, M_crit, depends on a non-ideal equation of state and on opacity changes due to grain growth, across a range of stellocentric distances from 5-100 AU. This minimum M_crit applies when planetesimal accretion does not substantially heat the atmosphere. Compared to an ideal gas polytrope, the inclusion of molecular hydrogen (H_2) dissociation and variable occupation of H_2 rotational states increases M_crit. Specifically, M_crit increases by a factor of ~2 if the H_2 spin isomers, ortho- and parahydrogen, are in thermal equilibrium, and by a factor of ~2-4 if the ortho-to-para ratio is fixed at 3:1. Lower opacities due to grain growth reduce M_crit. For a standard disk model around a Solar mass star, we calculate M_crit ~ 8 M_Earth at 5 AU, decreasing to ~5 M_Earth at 100 AU, for a realistic EOS with an equilibrium ortho-to-para ratio and for grain growth to cm-sizes. If grain coagulation is taken into account, M_crit may further reduce by up to one order of magnitude. These results for the minimum critical core mass are useful for the interpretation of surveys that find exoplanets at a range of orbital distances. △ Less

Submitted 16 December, 2014; originally announced December 2014.

Comments: 21 pages, 17 figures. Accepted to ApJ

Large eccentricity, low mutual inclination: the three-dimensional architecture of a hierarchical system of giant planets

Authors: Rebekah I. Dawson, John Asher Johnson, Daniel C. Fabrycky, Daniel Foreman-Mackey, Ruth A. Murray-Clay, Lars A. Buchhave, Phillip A. Cargile, Kelsey I. Clubb, Benjamin J. Fulton, Leslie Hebb, Andrew W. Howard, Daniel Huber, Avi Shporer, Jeff A. Valenti

Abstract: We establish the three-dimensional architecture of the Kepler-419 (previously KOI-1474) system to be eccentric yet with a low mutual inclination. Kepler-419b is a warm Jupiter at semi-major axis a = 0.370 +0.007/-0.006 AU with a large eccentricity e=0.85 +0.08/-0.07 measured via the "photoeccentric effect." It exhibits transit timing variations induced by the non-transiting Kepler-419c, which we u… ▽ More We establish the three-dimensional architecture of the Kepler-419 (previously KOI-1474) system to be eccentric yet with a low mutual inclination. Kepler-419b is a warm Jupiter at semi-major axis a = 0.370 +0.007/-0.006 AU with a large eccentricity e=0.85 +0.08/-0.07 measured via the "photoeccentric effect." It exhibits transit timing variations induced by the non-transiting Kepler-419c, which we uniquely constrain to be a moderately eccentric (e=0.184 +/- 0.002), hierarchically-separated (a=1.68 +/- 0.03 AU) giant planet (7.3 +/- 0.4 MJup). We combine sixteen quarters of Kepler photometry, radial-velocity (RV) measurements from the HIgh Resolution Echelle Spectrometer (HIRES) on Keck, and improved stellar parameters that we derive from spectroscopy and asteroseismology. From the RVs, we measure the mass of inner planet to be 2.5+/-0.3MJup and confirm its photometrically-measured eccentricity, refining the value to e=0.83+/-0.01. The RV acceleration is consistent with the properties of the outer planet derived from TTVs. We find that, despite their sizable eccentricities, the planets are coplanar to within 9 +8/-6 degrees, and therefore the inner planet's large eccentricity and close-in orbit are unlikely to be the result of Kozai migration. Moreover, even over many secular cycles, the inner planet's periapse is most likely never small enough for tidal circularization. Finally, we present and measure a transit time and impact parameter from four simultaneous ground-based light curves from 1m-class telescopes, demonstrating the feasibility of ground-based follow-up of Kepler giant planets exhibiting large TTVs. △ Less

Submitted 11 July, 2014; v1 submitted 20 May, 2014; originally announced May 2014.

Comments: ApJ, in press. 12 pages plus appendices/references. Submitted March 12, 2014. Accepted June 12, 2014. Very similar to previous version except revised Figure 10 (diagram of angles), fixed incorrect bolding, and system renamed Kepler-419

A Combined VLT and Gemini Study of the Atmosphere of the Directly-Imaged Planet, beta Pictoris b

Authors: Thayne Currie, Adam Burrows, Nikku Madhusudhan, Misato Fukagawa, Julien H. Girard, Rebekah Dawson, Ruth Murray-Clay, Scott Kenyon, Marc Kuchner, Soko Matsumura, Ray Jayawardhana, John Chambers, Ben Bromley

Abstract: We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet $β$ Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve high signal-to-noise, high precision measurements. While $β$ Pic b's near-IR colors mimick that of a standard, cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared comp… ▽ More We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet $β$ Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve high signal-to-noise, high precision measurements. While $β$ Pic b's near-IR colors mimick that of a standard, cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared compared to the field L/T dwarf sequence. Few substellar/planet-mass objects -- i.e. $κ$ And b and 1RXJ 1609B -- match $β$ Pic b's $JHK_{s}L^\prime$ photometry, and its 3.1 $μm$ and 5 $μm$ photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large ($\sim$ 60 $μm$) dust grains fail to reproduce the $β$ Pic b spectrum. However, models incorporating thick clouds similar to those found for HR 8799 bcde but also with small (a few microns) modal particle sizes yield fits consistent with the data within uncertainties. Assuming solar abundance models, thick clouds, and small dust particles ($<a>$ = 4 $μm$) we derive atmosphere parameters of log(g) = 3.8 $\pm$ 0.2 and $T_{eff}$ = 1575--1650 $K$, an inferred mass of 7$^{+4}_{-3}$ $M_{J}$, and a luminosity of log(L/L$_{\odot}$) $\sim$ -3.80 $\pm$ 0.02. The best-estimated planet radius, $\approx$ 1.65 $\pm$ 0.06 $R_{J}$, is near the upper end of allowable planet radii for hot-start models given the host star's age and likely reflects challenges with constructing accurate atmospheric models. Alternatively, these radii are comfortably consistent with hot-start model predictions if $β$ Pic b is younger than $\approx$ 7 Myr, consistent with a late formation, well after its host star's birth $\sim$ 12$^{+8}_{-4}$ Myr ago. △ Less

Submitted 13 August, 2013; v1 submitted 3 June, 2013; originally announced June 2013.

Comments: 24 pages, minor changes from previous version, Accepted for publication in ApJ

Giant planets orbiting metal-rich stars show signatures of planet-planet interactions

Authors: Rebekah I. Dawson, Ruth A. Murray-Clay

Abstract: Gas giants orbiting interior to the ice line are thought to have been displaced from their formation locations by processes that remain debated. Here we uncover several new metallicity trends, which together may indicate that two competing mechanisms deliver close-in giant planets: gentle disk migration, operating in environments with a range of metallicities, and violent planet-planet gravitation… ▽ More Gas giants orbiting interior to the ice line are thought to have been displaced from their formation locations by processes that remain debated. Here we uncover several new metallicity trends, which together may indicate that two competing mechanisms deliver close-in giant planets: gentle disk migration, operating in environments with a range of metallicities, and violent planet-planet gravitational interactions, primarily triggered in metal-rich systems in which multiple giant planets can form. First, we show with 99.1% confidence that giant planets with semi-major axes between 0.1 and 1 AU orbiting metal-poor stars ([Fe/H]<0) are confined to lower eccentricities than those orbiting metal-rich stars. Second, we show with 93.3% confidence that eccentric proto-hot Jupiters undergoing tidal circularization primarily orbit metal-rich stars. Finally, we show that only metal-rich stars host a pile-up of hot Jupiters, helping account for the lack of such a pile-up in the overall Kepler sample. Migration caused by stellar perturbers (e.g. stellar Kozai) is unlikely to account for the trends. These trends further motivate follow-up theoretical work addressing which hot Jupiter migration theories can also produce the observed population of eccentric giant planets between 0.1 and 1 AU. △ Less

Submitted 5 July, 2013; v1 submitted 25 February, 2013; originally announced February 2013.

Comments: ApJL; received 2013 February 6; accepted 2013 February 22; published 2013 April 2; 7 pages; 5 figures

A Thermal Infrared Imaging Study of Very Low-Mass, Wide Separation Brown Dwarf Companions to Upper Scorpius Stars: Constraining Circumstellar Environments

Authors: Vanessa Bailey, Philip M. Hinz, Thayne Currie, Kate Y. L. Su, Simone Esposito, John M. Hill, William F. Hoffmann, Terry Jones, Jihun Kim, Jarron Leisenring, Michael Meyer, Ruth Murray-Clay, Matthew J. Nelson, Enrico Pinna, Alfio Puglisi, George Rieke, Timothy Rodigas, Andrew Skemer, Michael F. Skrutskie, Vidhya Vaitheeswaran, John C. Wilson

Abstract: We present a 3-5um LBT/MMT adaptive optics imaging study of three Upper Scorpius stars with brown dwarf (BD) companions with very low-masses/mass ratios (M_BD < 25M_Jup; M_BD / M_star ~ 1-2%), and wide separations (300-700 AU): GSC 06214, 1RXS 1609, and HIP 78530. We combine these new thermal IR data with existing 1-4um and 24um photometry to constrain the properties of the BDs and identify eviden… ▽ More We present a 3-5um LBT/MMT adaptive optics imaging study of three Upper Scorpius stars with brown dwarf (BD) companions with very low-masses/mass ratios (M_BD < 25M_Jup; M_BD / M_star ~ 1-2%), and wide separations (300-700 AU): GSC 06214, 1RXS 1609, and HIP 78530. We combine these new thermal IR data with existing 1-4um and 24um photometry to constrain the properties of the BDs and identify evidence for circumprimary/secondary disks in these unusual systems. We confirm that GSC 06214B is surrounded by a disk, further showing this disk produces a broadband IR excess due to small dust near the dust sublimation radius. An unresolved 24um excess in the system may be explained by the contribution from this disk. 1RXS 1609B exhibits no 3-4um excess, nor does its primary; however, the system as a whole has a modest 24um excess, which may come from warm dust around the primary and/or BD. Neither object in the HIP 78530 system exhibits near- to mid-IR excesses. We additionally find that the 1-4um colors of HIP 78530B match a spectral type of M3+-2, inconsistent with the M8 spectral type assigned based on its near-IR spectrum, indicating it may be a low-mass star rather than a BD. We present new upper limits on additional low-mass companions in the system (<5M_Jup beyond 175AU). Finally, we examine the utility of circumsecondary disks as probes of the formation histories of wide BD companions, finding that the presence of a disk may disfavor BD formation near the primary with subsequent outward scattering. △ Less

Submitted 3 February, 2013; originally announced February 2013.

Comments: Accepted for publication in ApJ (14 pages, 9 figures)