sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Murray-Clay, R.

Normalized to: Murray-Clay, R.

62 article(s) in total. 362 co-authors, from 1 to 12 common article(s). Median position in authors list is 3,0.

[1]  oai:arXiv.org:2006.10674  [pdf] - 2117757
OSSOS: The Bimodal Eccentricity and Inclination Distributions of the Stable Neptunian Trojans
Comments: 13 pages, 6 figures
Submitted: 2020-06-18
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 $\sigma_{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 $\sigma_e$ of 0.045 and a typical sin(i) x gaussian inclination distribution with a width $\sigma_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 $\sigma_{e-cold}$ ~ 0.02 / $\sigma_{i-cold}$ ~ 6$^\circ$ and $\sigma_{e-hot}$~ 0.05 / $\sigma_{i-hot}$ ~ 18$^\circ$. In this scenario, the `cold' L4 Neptunian Trojan population lacks the Hr >= 8 member 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.
[2]  oai:arXiv.org:2004.13720  [pdf] - 2119967
How Consumption and Repulsion Set Planetary Gap Depths and the Final Masses of Gas Giants
Comments: Accepted for publication in MNRAS
Submitted: 2020-04-28, last modified: 2020-06-12
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.
[3]  oai:arXiv.org:1906.00075  [pdf] - 2046213
Effects of Radiation Pressure on the Evaporative Wind of HD 209458b
Comments:
Submitted: 2019-05-31, last modified: 2020-02-03
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-$\alpha$ 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-$\alpha$ 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.
[4]  oai:arXiv.org:2001.10552  [pdf] - 2109963
Requirements for gravitational collapse in planetesimal formation --- the impact of scales set by Kelvin-Helmholtz and nonlinear streaming instability
Comments: 26 pages, 11 figures, submitted to ApJ
Submitted: 2020-01-28
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 \beta \epsilon^{-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 $\beta$, metallicity $Z$ and local particle enhancement $\epsilon$, 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 $\epsilon$, can evolve a system to instability.
[5]  oai:arXiv.org:1910.12985  [pdf] - 2026175
OSSOS XVIII: Constraining migration models with the 2:1 resonance using the Outer Solar System Origins Survey
Comments: 23 pages, 12 figures, 2 tables, accepted for publication in the Astronomical Journal
Submitted: 2019-10-28, last modified: 2019-10-31
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 $\sigma_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 $\alpha=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.
[6]  oai:arXiv.org:1906.03266  [pdf] - 1994056
Signatures of a planet-planet impacts phase in exoplanetary systems hosting giant planets
Comments: 8 pages, 5 figures, accepted to ApJL
Submitted: 2019-06-07, last modified: 2019-10-08
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.
[7]  oai:arXiv.org:1908.06991  [pdf] - 1945929
How Flow Isolation May Set the Mass Scale for Super-Earth Planets
Comments: 14 pages, 6 figures
Submitted: 2019-08-19
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.
[8]  oai:arXiv.org:1704.04693  [pdf] - 1944189
Using Ice and Dust Lines to Constrain the Surface Densities of Protoplanetary Disks
Comments: Published in ApJ, May 2017
Submitted: 2017-04-15, last modified: 2019-08-16
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.
[9]  oai:arXiv.org:1908.04789  [pdf] - 1966952
Measuring the Orbital Parameters of Radial Velocity Systems in Mean Motion Resonance---a Case Study of HD 200964
Comments: 26 pages, 23 figures
Submitted: 2019-08-13
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.
[10]  oai:arXiv.org:1905.04275  [pdf] - 1881118
Directly Imaging Rocky Planets from the Ground
Comments: 8 pages, 1 figure, Astro2020 Science White Paper
Submitted: 2019-05-10
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
[11]  oai:arXiv.org:1905.03252  [pdf] - 1905805
New Constraints From Dust Lines On The Surface Densities Of Protoplanetary Disks
Comments: 29 pages, 12 figures, submitted to ApJ, revised following comments from referee
Submitted: 2019-05-08
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.
[12]  oai:arXiv.org:1903.08187  [pdf] - 1853143
Wide-Orbit Exoplanet Demographics
Comments: Astro2020 Science White paper
Submitted: 2019-03-19
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.
[13]  oai:arXiv.org:1903.05319  [pdf] - 1848658
Imaging Giant Protoplanets with the ELTs
Comments: Science white paper submitted to the Astro2020 Decadal Survey
Submitted: 2019-03-13
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.
[14]  oai:arXiv.org:1811.09276  [pdf] - 1851599
Morphology of Hydrodynamic Winds: A Study of Planetary Winds in Stellar Environments
Comments: 33 pages, 21 figures (7 of which have embedded movies viewable with Adobe Acrobat Pro), Submitted to ApJ
Submitted: 2018-11-22
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-$\alpha$ 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.
[15]  oai:arXiv.org:1811.09314  [pdf] - 1800001
Photoevaporative Flows From Exoplanet Atmospheres: A 3-D Radiative Hydrodynamic Parameter Study
Comments:
Submitted: 2018-11-22
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.
[16]  oai:arXiv.org:1808.05947  [pdf] - 1743872
Restrictions on the Growth of Gas Giant Cores via Pebble Accretion
Comments:
Submitted: 2018-08-17
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 ($\alpha \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 $\alpha$ 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 ($\alpha \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.
[17]  oai:arXiv.org:1807.06033  [pdf] - 1725079
Metallicity-Dependent Signatures in the Kepler Planets
Comments: To appear in MNRAS
Submitted: 2018-07-16
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.
[18]  oai:arXiv.org:1805.11740  [pdf] - 1691625
OSSOS. VII. 800+ trans-Neptunian objects - the complete data release
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)
Submitted: 2018-05-29
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 $\sigma_{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.
[19]  oai:arXiv.org:1802.05805  [pdf] - 1692470
OSSOS IX: two objects in Neptune's 9:1 resonance -- implications for resonance sticking in the scattering population
Comments: accepted for publication in AJ
Submitted: 2018-02-15, last modified: 2018-05-24
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.
[20]  oai:arXiv.org:1805.08228  [pdf] - 1717149
Trans-Neptunian Objects Transiently Stuck in Neptune's Mean Motion Resonances: Numerical simulations of the current population
Comments:
Submitted: 2018-05-21
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.
[21]  oai:arXiv.org:1805.06898  [pdf] - 1717140
Gas-Assisted Growth of Protoplanets in a Turbulent Medium
Comments: 40 pages, 20 figures
Submitted: 2018-05-17
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 ($\alpha \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.
[22]  oai:arXiv.org:1705.06851  [pdf] - 1583632
Evidence that the Directly-Imaged Planet HD 131399 Ab is a Background Star
Comments: 35 pages, 20 figures, Published in AJ
Submitted: 2017-05-18, last modified: 2017-11-08
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.
[23]  oai:arXiv.org:1709.03530  [pdf] - 1588225
Pebble Accretion in Turbulent Protoplanetary Disks
Comments: 15 pages, 7 figures, accepted for publication in ApJ
Submitted: 2017-09-11
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.
[24]  oai:arXiv.org:1708.00862  [pdf] - 1586693
The Formation of Uranus and Neptune: Fine Tuning in Core Accretion
Comments: Accepted for publication in AJ, 10 pages, 3 figures
Submitted: 2017-08-02
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.
[25]  oai:arXiv.org:1611.08616  [pdf] - 1557757
Metal enrichment leads to low atmospheric C/O ratios in transiting giant exoplanets
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
Submitted: 2016-11-25, last modified: 2017-03-10
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 \"Oberg 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.
[26]  oai:arXiv.org:1607.06970  [pdf] - 1531102
OSSOS: IV. Discovery of a dwarf planet candidate in the 9:2 resonance with Neptune
Comments: Accepted to AJ
Submitted: 2016-07-23, last modified: 2016-10-05
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.
[27]  oai:arXiv.org:1511.07854  [pdf] - 1509263
The Mass-Metallicity Relation for Giant Planets
Comments: Accepted to The Astrophysical Journal. This revision adds a substantial amount of discussion; the results are the same
Submitted: 2015-11-24, last modified: 2016-08-30
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}$.
[28]  oai:arXiv.org:1511.02895  [pdf] - 1483231
The Outer Solar System Origins Survey: I. Design and First-Quarter Discoveries
Comments: Accepted to AJ, 27 April 2016. 59 pp
Submitted: 2015-11-09, last modified: 2016-05-10
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.
[29]  oai:arXiv.org:1604.08177  [pdf] - 1433043
OSSOS III - Resonant Trans-Neptunian Populations: Constraints from the first quarter of the Outer Solar System Origins Survey
Comments: Accepted for publication in AJ
Submitted: 2016-04-27
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.
[30]  oai:arXiv.org:1511.05563  [pdf] - 1331145
C/O and Snowline Locations in Protoplanetary Disks: The Effect of Radial Drift and Viscous Gas Accretion
Comments: 16 pages, 8 figures. Accepted to ApJ
Submitted: 2015-11-17
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.
[31]  oai:arXiv.org:1508.03084  [pdf] - 1290800
Discovery and spectroscopy of the young Jovian planet 51 Eri b with the Gemini Planet Imager
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
Submitted: 2015-08-12, last modified: 2015-10-09
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.
[32]  oai:arXiv.org:1510.02747  [pdf] - 1290878
Direct imaging of an asymmetric debris disk in the HD 106906 planetary system
Comments: Accepted for publication in the Astrophysical Journal; 15 pages, 7 figures
Submitted: 2015-10-09
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$\deg$ 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.
[33]  oai:arXiv.org:1506.06759  [pdf] - 1255200
Simulated Photoevaporative Mass Loss from Hot Jupiters in 3D
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
Submitted: 2015-06-22, last modified: 2015-08-03
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.
[34]  oai:arXiv.org:1412.5185  [pdf] - 1223675
Minimum Core Masses for Giant Planet Formation With Realistic Equations of State and Opacities
Comments: 21 pages, 17 figures. Accepted to ApJ
Submitted: 2014-12-16
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.
[35]  oai:arXiv.org:1211.0554  [pdf] - 919819
A paucity of proto-hot Jupiters on super-eccentric orbits
Comments: Accepted to ApJ on October 20, 2014. We adopted more conservative assumptions about completeness and waited for a candidate sample based on double the number of Kepler quarters. With the updated sample and assumptions, the conclusions remain the same
Submitted: 2012-11-02, last modified: 2014-10-21
Gas giant planets orbiting within 0.1 AU of their host stars, unlikely to have formed in situ, are evidence for planetary migration. It is debated whether the typical hot Jupiter smoothly migrated inward from its formation location through the proto-planetary disk or was perturbed by another body onto a highly eccentric orbit, which tidal dissipation subsequently shrank and circularized during close stellar passages. Socrates and collaborators predicted that the latter class of model should produce a population of super-eccentric proto-hot Jupiters readily observable by Kepler. We find a paucity of such planets in the Kepler sample, inconsistent with the theoretical prediction with 96.9% confidence. Observational effects are unlikely to explain this discrepancy. We find that the fraction of hot Jupiters with orbital period P > 3 days produced by the star-planet Kozai mechanism does not exceed (at two-sigma) 44%. Our results may indicate that disk migration is the dominant channel for producing hot Jupiters with P > 3 days. Alternatively, the typical hot Jupiter may have been perturbed to a high eccentricity by interactions with a planetary rather than stellar companion and began tidal circularization much interior to 1 AU after multiple scatterings. A final alternative is that tidal circularization occurs much more rapidly early in the tidal circularization process at high eccentricities than later in the process at low eccentricities, contrary to current tidal theories.
[36]  oai:arXiv.org:1405.5229  [pdf] - 859872
Large eccentricity, low mutual inclination: the three-dimensional architecture of a hierarchical system of giant planets
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
Submitted: 2014-05-20, last modified: 2014-07-11
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.
[37]  oai:arXiv.org:1306.0610  [pdf] - 1171770
A Combined VLT and Gemini Study of the Atmosphere of the Directly-Imaged Planet, beta Pictoris b
Comments: 24 pages, minor changes from previous version, Accepted for publication in ApJ
Submitted: 2013-06-03, last modified: 2013-08-13
We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet $\beta$ 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 $\beta$ 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. $\kappa$ And b and 1RXJ 1609B -- match $\beta$ Pic b's $JHK_{s}L^\prime$ photometry, and its 3.1 $\mu m$ and 5 $\mu m$ photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large ($\sim$ 60 $\mu m$) dust grains fail to reproduce the $\beta$ 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 $\mu 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 $\beta$ 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.
[38]  oai:arXiv.org:1302.6244  [pdf] - 689064
Giant planets orbiting metal-rich stars show signatures of planet-planet interactions
Comments: ApJL; received 2013 February 6; accepted 2013 February 22; published 2013 April 2; 7 pages; 5 figures
Submitted: 2013-02-25, last modified: 2013-07-05
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.
[39]  oai:arXiv.org:1302.0582  [pdf] - 1159467
A Thermal Infrared Imaging Study of Very Low-Mass, Wide Separation Brown Dwarf Companions to Upper Scorpius Stars: Constraining Circumstellar Environments
Comments: Accepted for publication in ApJ (14 pages, 9 figures)
Submitted: 2013-02-03
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.
[40]  oai:arXiv.org:1206.5579  [pdf] - 599732
The Photoeccentric Effect and Proto-Hot Jupiters II. KOI-1474.01, a candidate eccentric planet perturbed by an unseen companion
Comments: ApJ, in press. Received 2012 July 7; accepted 2012 October 13
Submitted: 2012-06-25, last modified: 2012-11-02
The exoplanets known as hot Jupiters---Jupiter-sized planets with periods less than 10 days---likely are relics of dynamical processes that shape all planetary system architectures. Socrates et al. (2012) argued that high eccentricity migration (HEM) mechanisms proposed for situating these close-in planets should produce an observable population of highly eccentric proto-hot Jupiters that have not yet tidally circularized. HEM should also create failed-hot Jupiters, with periapses just beyond the influence of fast circularization. Using the technique we previously presented for measuring eccentricities from photometry (the "photoeccentric effect"), we are distilling a collection of eccentric proto- and failed-hot Jupiters from the Kepler Objects of Interest (KOI). Here we present the first, KOI-1474.01, which has a long orbital period (69.7340 days) and a large eccentricity e = 0.81+0.10/-0.07, skirting the proto-hot Jupiter boundary. Combining Kepler photometry, ground-based spectroscopy, and stellar evolution models, we characterize host KOI-1474 as a rapidly-rotating F-star. Statistical arguments reveal that the transiting candidate has a low false-positive probability of 3.1%. KOI-1474.01 also exhibits transit timing variations of order an hour. We explore characteristics of the third-body perturber, which is possibly the "smoking-gun" cause of KOI-1474.01's large eccentricity. Using the host-star's rotation period, radius, and projected rotational velocity, we find KOI-1474.01's orbit is marginally consistent with aligned with the stellar spin axis, although a reanalysis is warranted with future additional data. Finally, we discuss how the number and existence of proto-hot Jupiters will not only demonstrate that hot Jupiters migrate via HEM, but also shed light on the typical timescale for the mechanism.
[41]  oai:arXiv.org:1112.4822  [pdf] - 1092535
Disruption of a Proto-Planetary Disk by the Black Hole at the Milky Way Centre
Comments: Accepted to Nature Communications; new Figure 4b provides predicted Br-gamma emission as a function of time
Submitted: 2011-12-20, last modified: 2012-08-09
Recently, an ionized cloud of gas was discovered plunging toward the supermassive black hole, SgrA*, at the centre of the Milky Way. The cloud is being tidally disrupted along its path to closest approach at ~3100 Schwarzschild radii from the black hole. Here, we show that the observed properties of this cloud of gas can naturally be produced by a proto-planetary disk surrounding a low-mass star, which was scattered from the observed ring of young stars orbiting SgrA*. As the young star approaches the black hole, its disk experiences both photo-evaporation and tidal disruption, producing a cloud. Our model implies that planets form in the Galactic centre, and that tidal debris from proto-planetary disks can flag low mass stars which are otherwise too faint to be detected.
[42]  oai:arXiv.org:1208.0832  [pdf] - 1150508
NIMBUS: The Near-Infrared Multi-Band Ultraprecise Spectroimager for SOFIA
Comments: 14 pages, 9 figures, SPIE Astronomical Telescopes and Instrumentation 2012
Submitted: 2012-08-03
We present a new and innovative near-infrared multi-band ultraprecise spectroimager (NIMBUS) for SOFIA. This design is capable of characterizing a large sample of extrasolar planet atmospheres by measuring elemental and molecular abundances during primary transit and occultation. This wide-field spectroimager would also provide new insights into Trans-Neptunian Objects (TNO), Solar System occultations, brown dwarf atmospheres, carbon chemistry in globular clusters, chemical gradients in nearby galaxies, and galaxy photometric redshifts. NIMBUS would be the premier ultraprecise spectroimager by taking advantage of the SOFIA observatory and state of the art infrared technologies. This optical design splits the beam into eight separate spectral bandpasses, centered around key molecular bands from 1 to 4 microns. Each spectral channel has a wide field of view for simultaneous observations of a reference star that can decorrelate time-variable atmospheric and optical assembly effects, allowing the instrument to achieve ultraprecise calibration for imaging and photometry for a wide variety of astrophysical sources. NIMBUS produces the same data products as a low-resolution integral field spectrograph over a large spectral bandpass, but this design obviates many of the problems that preclude high-precision measurements with traditional slit and integral field spectrographs. This instrument concept is currently not funded for development.
[43]  oai:arXiv.org:1202.6060  [pdf] - 500375
Neptune's wild days: constraints from the eccentricity distribution of the classical Kuiper Belt
Comments: Corrected typos and made wording changes. Corrected Fig. 8 (row 2) and Fig. 17. Reduced loading time of Fig. 13
Submitted: 2012-02-27, last modified: 2012-04-16
Neptune's dynamical history shaped the current orbits of Kuiper belt objects (KBOs), leaving clues to the planet's orbital evolution. In the "classical" region, a population of dynamically "hot" high-inclination KBOs overlies a flat "cold" population with distinct physical properties. Simulations of qualitatively different histories for Neptune -including smooth migration on a circular orbit or scattering by other planets to a high eccentricity - have not simultaneously produced both populations. We explore a general Kuiper belt assembly model that forms hot classical KBOs interior to Neptune and delivers them to the classical region, where the cold population forms in situ. First, we present evidence that the cold population is confined to eccentricities well below the limit dictated by long-term survival. Therefore Neptune must deliver hot KBOs into the long-term survival region without excessively exciting the eccentricities of the cold population. Imposing this constraint, we explore the parameter space of Neptune's eccentricity and eccentricity damping, migration, and apsidal precession. We rule out much of parameter space, except where Neptune is scattered to a moderately eccentric orbit (e > 0.15) and subsequently migrates a distance Delta aN=1-6 AU. Neptune's moderate eccentricity must either damp quickly or be accompanied by fast apsidal precession. We find that Neptune's high eccentricity alone does not generate a chaotic sea in the classical region. Chaos can result from Neptune's interactions with Uranus, exciting the cold KBOs and placing additional constraints. Finally, we discuss how to interpret our constraints in the context of the full, complex dynamical history of the solar system.
[44]  oai:arXiv.org:1112.1954  [pdf] - 481596
Neptune on tiptoes: dynamical histories that preserve the cold classical Kuiper belt
Comments: Accepted to ApJ. Update: Fixed typos in text and Appendix equations. Published as ApJ, 746, 171
Submitted: 2011-12-08, last modified: 2012-02-27
The current dynamical structure of the Kuiper belt was shaped by the orbital evolution of the giant planets, especially Neptune, during the era following planet formation, when the giant planets may have undergone planet-planet scattering and/or planetesimal-driven migration. Numerical simulations of this process, while reproducing many properties of the belt, fail to generate the high inclinations and eccentricities observed for some objects while maintaining the observed dynamically "cold" population. We present the first of a three-part parameter study of how different dynamical histories of Neptune sculpt the planetesimal disk. Here we identify which dynamical histories allow an in situ planetesimal disk to remain dynamically cold, becoming today's cold Kuiper belt population. We find that if Neptune undergoes a period of elevated eccentricity and/or inclination, it secularly excites the eccentricities and inclinations of the planetesimal disk. We demonstrate that there are several well-defined regimes for this secular excitation, depending on the relative timescales of Neptune's migration, the damping of Neptune's orbital inclination and/or eccentricity, and the secular evolution of the planetesimals. We model this secular excitation analytically in each regime, allowing for a thorough exploration of parameter space. Neptune's eccentricity and inclination can remain high for a limited amount of time without disrupting the cold classical belt. In the regime of slow damping and slow migration, if Neptune is located (for example) at 20 AU, then its eccentricity must stay below 0.18 and its inclination below 6{\deg}.
[45]  oai:arXiv.org:1111.0297  [pdf] - 442311
On the misalignment of the directly imaged planet {\beta} Pictoris b with the system's warped inner disk
Comments: 6 pages, 4 figures; accepted for publication in the Astrophysical Journal Letters. Update: Compressed figures and published as ApJL 743, L17
Submitted: 2011-11-01, last modified: 2011-11-23
The vertical warp in the debris disk Beta Pictoris -- an inclined inner disk extending into a flat outer disk -- has long been interpreted as the signpost of a planet on an inclined orbit. Direct images spanning 2004-2010 have revealed Beta Pictoris b, a planet with a mass and orbital distance consistent with this picture. However, it was recently reported that the orbit of planet b is aligned with the flat outer disk, not the inclined inner disk, and thus lacks the inclination to warp the disk. We explore three scenarios for reconciling the apparent misalignment of the directly imaged planet Beta Pictoris b with the warped inner disk of Beta Pictoris: observational uncertainty, an additional planet, and damping of planet b's inclination. We find that, at the extremes of the uncertainties, the orbit of Beta Pictoris b has the inclination necessary to produce the observed warp. We also find that if planet b were aligned with the flat outer disk, it would prevent another planet from creating a warp with the observed properties; therefore planet b itself must be responsible for the warp. Finally, planet b's inclination could have been damped by dynamical friction and still produce the observed disk morphology, but the feasibility of damping depends on disk properties and the presence of other planets. More precise observations of the orbit of planet b and the position angle of the outer disk will allow us to distinguish between the first and third scenario.
[46]  oai:arXiv.org:1110.5567  [pdf] - 1085130
The effects of snowlines on C/O in planetary atmospheres
Comments: Accepted for publication in ApJL
Submitted: 2011-10-25
The C/O ratio is predicted to regulate the atmospheric chemistry in hot Jupiters. Recent observations suggest that some exo-planets, e.g. Wasp 12- b, have atmospheric C/O ratios substantially different from the solar value of 0.54. In this paper we present a mechanism that can produce such atmospheric deviations from the stellar C/O ratio. In protoplanetary disks, different snowlines of oxygen- and carbon-rich ices, especially water and carbon monoxide, will result in systematic variations in the C/O ratio both in the gas and in the condensed phase. In particular, between the H2O and CO snowlines most oxygen is present in icy grains - the building blocks of planetary cores in the core accretion model - while most carbon remains in the gas-phase. This region is coincidental with the giant-planet forming zone for a range of observed protoplanetary disks. Based on standard core accretion models of planet formation, gas giants that sweep up most of their atmospheres from disk gas outside of the water snowline will have C/O?1, while atmospheres significantly contaminated by evaporating planetesimals will have stellar or sub-stellar C/O when formed at the same disk radius. The overall metallicity will also depend on the atmosphere formation mechanism, and exoplanetary atmospheric compositions may therefore provide constraints on where and how a specific planet formed.
[47]  oai:arXiv.org:1107.3153  [pdf] - 1078012
Resolved Submillimeter Observations of the HR 8799 and HD 107146 Debris Disks
Comments: 9 pages, 3 figures, accepted for publication in ApJ
Submitted: 2011-07-15
We present 880 um Submillimeter Array observations of the debris disks around the young solar analogue HD 107146 and the multiple-planet host star HR 8799, at an angular resolution of 3" and 6", respectively. We spatially resolve the inner edge of the disk around HR 8799 for the first time. While the data are not sensitive enough (with rms noise of 1 mJy) to constrain the system geometry, we demonstrate that a model by Su et al. (2009) based on the spectral energy distribution (SED) with an inner radius of 150 AU predicts well the spatially resolved data. Furthermore, by modeling simultaneously the SED and visibilities, we demonstrate that the dust is distributed in a broad (of order 100 AU) annulus rather than a narrow ring. We also model the observed SED and visibilities for the HD 107146 debris disk and generate a model of the dust emission that extends in a broad band between 50 and 170 AU from the star. We perform an a posteriori comparison with existing 1.3 mm CARMA observations and demonstrate that a smooth, axisymmetric model reproduces well all of the available millimeter-wavelength data.
[48]  oai:arXiv.org:1107.0728  [pdf] - 1077752
Fragment Production and Survival in Irradiated Disks: A Comprehensive Cooling Criterion
Comments: 14 pages, 4 figures, accepted to ApJ
Submitted: 2011-07-04
Accretion disks that become gravitationally unstable can fragment into stellar or sub-stellar companions. The formation and survival of these fragments depends on the precarious balance between self-gravity, internal pressure, tidal shearing, and rotation. Disk fragmentation depends on two key factors (1) whether the disk can get to the fragmentation boundary of Q=1, and (2) whether fragments can survive for many orbital periods. Previous work suggests that to reach Q=1, and have fragments survive, a disk must cool on an orbital timescale. Here we show that disks heated primarily by external irradiation always satisfy the standard cooling time criterion. Thus even though irradiation heats disks, and makes them more stable in general, once they reach the fragmentation boundary, they fragment more easily. We derive a new cooling criterion that determines fragment survival, and calculate a pressure modified Hill radius, which sets the maximum size of pressure-supported objects in a Keplerian disk. We conclude that fragmentation in protostellar disks might occur at slightly smaller radii than previously thought, and recommend tests for future simulations that will better predict the outcome of fragmentation in real disks.
[49]  oai:arXiv.org:1103.1629  [pdf] - 1052572
Wind-shearing in gaseous protoplanetary disks and the evolution of binary planetesimals
Comments: ApJ, in press
Submitted: 2011-03-08, last modified: 2011-03-22
One of the first stages of planet formation is the growth of small planetesimals. This early stage occurs much before the dispersal of most of the gas from the protoplanetary disk. Due to their different aerodynamic properties, planetesimals of different sizes and shapes experience different drag forces from the gas during this time. Such differential forces produce a wind-shearing (WISH) effect between close by, different size planetesimals. For any two planetesimals, a WISH radius can be considered, at which the differential acceleration due to the wind becomes greater than the mutual gravitational pull between the planetesimals. We find that the WISH radius could be much smaller than the Hill radius, i.e. WISH could play a more important role than tidal perturbations by the star. Here we study the WISH radii for planetesimal pairs of different sizes and compare the effects of wind and gravitational shearing (drag force vs. gravitational tidal force). We then discuss the role of WISH for the stability and survival of binary planetesimals. Binaries are sheared apart by the wind if they are wider than their WISH radius. WISH-stable binaries can inspiral and possibly coalesce due to gas drag. Here, we calculate the WISH radius and the gas drag-induced merger timescale, providing stability and survival criteria for gas-embedded binary planetesimals. Our results suggest that even WISH-stable binaries may merge in times shorter than the lifetime of the gaseous disk. This may constrain currently observed binary planetesimals to have formed far from the star or at a late stage after the dispersal of most of the disk gas. We note that the WISH radius may also be important for other processes such as planetesimal erosion and planetesimal encounters and collisions in a gaseous environment.
[50]  oai:arXiv.org:1102.1430  [pdf] - 1051953
Using Kuiper Belt Binaries to Constrain Neptune's Migration History
Comments: Accepted to ApJ, 15 pages, 4 figures
Submitted: 2011-02-07
Approximately 10-20% of all Kuiper belt objects (KBOs) occupy mean-motion resonances with Neptune. This dynamical configuration likely resulted from resonance capture as Neptune migrated outward during the late stages of planet formation. The details of Neptune's planetesimal-driven migration, including its radial extent and the concurrent eccentricity evolution of the planet, are the subject of considerable debate. Two qualitatively different proposals for resonance capture have been proposed--migration-induced capture driven by smooth outward evolution of Neptune's orbit and chaotic capture driven by damping of the planet's eccentricity near its current semi-major axis. We demonstrate that the distribution of comparable-mass, wide-separation binaries occupying resonant orbits can differentiate between these two scenarios. If migration-induced capture occurred, this fraction records information about the formation locations of different populations of KBOs. Chaotic capture, in contrast, randomizes the orbits of bodies as they are placed in resonance. In particular, migration-induced capture produces the following signatures. The 2:1 resonance should contain a dynamically cold component, with inclinations less than 5-10 degrees, having a binary fraction comparable to that among cold classical KBOs. If the 3:2 resonance also hosts a cold component, its binary fraction should be 20-30% lower than in the cold classical belt. Among cold 2:1 (and if present 3:2) KBOs, objects with eccentricities e < 0.2 should have a binary fraction ~20% larger than those with e > 0.2. Searches for cold components in the binary fractions of resonant KBOs are currently practical. The additional migration-generated trends described here may be distinguished with objects discovered by LSST. (Abstract abridged.)
[51]  oai:arXiv.org:1012.0571  [pdf] - 1042343
Wind-shearing in gaseous protoplanetary disks
Comments: To be published in the proceedings of IAU 276: The Astrophysics of planetary systems - formation, structure, and dynamical evolution
Submitted: 2010-12-02
One of the first stages of planet formation is the growth of small planetesimals and their accumulation into large planetesimals and planetary embryos. This early stage occurs much before the dispersal of most of the gas from the protoplanetary disk. Due to their different aerodynamic properties, planetesimals of different sizes/shapes experience different drag forces from the gas at these stage. Such differential forces produce a wind-shearing effect between close by, different size planetesimals. For any two planetesimals, a wind-shearing radius can be considered, at which the differential acceleration due to the wind becomes greater than the mutual gravitational pull between the planetesimals. We find that the wind-shearing radius could be much smaller than the gravitational shearing radius by the Sun (the Hill radius), i.e. during the gas-phase of the disk wind-shearing could play a more important role than tidal perturbations by the Sun. Here we study the wind-shearing radii for planetesimal pairs of different sizes and compare it with gravitational shearing (drag force vs. gravitational tidal forces). We then discuss the role of wind-shearing for the stability and survival of binary planetesimals, and provide stability criteria for binary planetesimals embedded in a gaseous disk.
[52]  oai:arXiv.org:1008.1983  [pdf] - 1034227
High-Contrast 3.8 Micron Imaging Of The Brown Dwarf/Planet-Mass Companion to GJ 758
Comments: 14 pages, 3 figures. Accepted for publication in The Astrophysical Journal Letters
Submitted: 2010-08-11, last modified: 2010-08-23
We present L' band (3.8 $\mu m$) MMT/Clio high-contrast imaging data for the nearby star GJ 758, which was recently reported by Thalmann et al. (2009) to have one -- possibly two-- faint comoving companions (GJ 758B and ``C", respectively). GJ 758B is detected in two distinct datasets. Additionally, we report a \textit{possible} detection of the object identified by Thalmann et al as ``GJ 758C" in our more sensitive dataset, though it is likely a residual speckle. However, if it is the same object as that reported by Thalmann et al. it cannot be a companion in a bound orbit. GJ 758B has a H-L' color redder than nearly all known L--T8 dwarfs. Based on comparisons with the COND evolutionary models, GJ 758B has T$_{e}$ $\sim$ 560 K$^{^{+150 K}_{-90K}}$ and a mass ranging from $\sim$ 10--20 M$_{J}$ if it is $\sim$ 1 Gyr old to $\sim$ 25--40 M$_{J}$ if it is 8.7 Gyr old. GJ 758B is likely in a highly eccentric orbit, e $\sim$ 0.73$^{^{+0.12}_{-0.21}}$, with a semimajor axis of $\sim$ 44 AU$^{^{+32 AU}_{-14 AU}}$. Though GJ 758B is sometimes discussed within the context of exoplanet direct imaging, its mass is likely greater than the deuterium-burning limit and its formation may resemble that of binary stars rather than that of jovian-mass planets.
[53]  oai:arXiv.org:0909.2644  [pdf] - 901721
The Runts of the Litter: Why planets formed through gravitational instability can only be failed binary stars
Comments: 14 pages, 6 figures, accepted to ApJ, some content changes, revised figures
Submitted: 2009-09-15, last modified: 2010-01-26
Recent direct imaging discoveries suggest a new class of massive, distant planets around A stars. These widely separated giants have been interpreted as signs of planet formation driven by gravitational instability, but the viability of this mechanism is not clear cut. In this paper, we first discuss the local requirements for fragmentation and the initial fragment mass scales. We then consider whether the fragment's subsequent growth can be terminated within the planetary mass regime. Finally, we place disks in the larger context of star formation and disk evolution models. We find that in order for gravitational instability to produce planets, disks must be atypically cold in order to reduce the initial fragment mass. In addition, fragmentation must occur during a narrow window of disk evolution, after infall has mostly ceased, but while the disk is still sufficiently massive to undergo gravitational instability. Under more typical conditions, disk-born objects will likely grow well above the deuterium burning planetary mass limit. We conclude that if planets are formed by gravitational instability, they must be the low mass tail of the distribution of disk-born companions. To validate this theory, on-going direct imaging surveys must find a greater abundance of brown dwarf and M-star companions to A-stars. Their absence would suggest planet formation by a different mechanism such as core accretion, which is consistent with the debris disks detected in these systems.
[54]  oai:arXiv.org:0812.0011  [pdf] - 315358
Stability of the directly imaged multiplanet system HR 8799: resonance and masses
Comments: 16 pages; 15 figures; 4 tables. Accepted to ApJ. New version gives more attention to the upper mass limit for a double resonance. A new section allows for arbitrary eccentricities, orientations, and mutual inclinations
Submitted: 2008-12-01, last modified: 2010-01-12
A new era of directly imaged extrasolar planets has produced a three-planet system (Marois et al. 2008), where the masses of the planets have been estimated by untested cooling models. We point out that the nominal circular, face-on orbits of the planets lead to a dynamical instability in ~10^5 yr, a factor of at least 100 shorter than the estimated age of the star. Reduced planetary masses produce stability only for unreasonably small planets (<~2 M_ Jup). Relaxing the face-on assumption, but still requiring circular orbits while fitting the observed positions, makes the instability time even shorter. A promising solution is that the inner two planets have a 2:1 commensurability between their periods, and they avoid close encounters with each other through this resonance. That the inner resonance has lasted until now, in spite of the perturbations of the outer planet, leads to a limit <~10 M_Jup on the masses unless the outer two planets are also engaged in a 2:1 mean-motion resonance. In a double resonance, which is consistent with the current data, the system could survive until now even if the planets have masses of ~20 M_Jup. Apsidal alignment can further enhance the stability of a mean-motion resonant system. A completely different dynamical configuration, with large eccentricities and large mutual inclinations among the planets, is possible but finely tuned.
[55]  oai:arXiv.org:0912.3229  [pdf] - 1018817
A super-Earth transiting a nearby low-mass star
Comments: 13 pages, 3 figures, published in Nature
Submitted: 2009-12-16
A decade ago, the detection of the first transiting extrasolar planet provided a direct constraint on its composition and opened the door to spectroscopic investigations of extrasolar planetary atmospheres. As such characterization studies are feasible only for transiting systems that are both nearby and for which the planet-to-star radius ratio is relatively large, nearby small stars have been surveyed intensively. Doppler studies and microlensing have uncovered a population of planets with minimum masses of 1.9-10 times the Earth's mass (M_Earth), called super-Earths. The first constraint on the bulk composition of this novel class of planets was afforded by CoRoT-7b, but the distance and size of its star preclude atmospheric studies in the foreseeable future. Here we report observations of the transiting planet GJ 1214b, which has a mass of 6.55 M_Earth and a radius 2.68 times Earth's radius (R_Earth), indicating that it is intermediate in stature between Earth and the ice giants of the Solar System. We find that the planetary mass and radius are consistent with a composition of primarily water enshrouded by a hydrogen-helium envelope that is only 0.05% of the mass of the planet. The atmosphere is probably escaping hydrodynamically, indicating that it has undergone significant evolution during its history. As the star is small and only 13 parsecs away, the planetary atmosphere is amenable to study with current observatories.
[56]  oai:arXiv.org:0811.0006  [pdf] - 18079
Atmospheric Escape from Hot Jupiters
Comments: Accepted for publication in ApJ, 22 pages, 24 figures. Added additional references--in particular, we now discuss charge-exchange as a potential source of high-velocity neutrals
Submitted: 2008-11-03, last modified: 2008-11-18
Photoionization heating from UV radiation incident on the atmospheres of hot Jupiters may drive planetary mass loss. We construct a model of escape that includes realistic heating and cooling, ionization balance, tidal gravity, and pressure confinement by the host star wind. We show that mass loss takes the form of a hydrodynamic ("Parker") wind, emitted from the planet's dayside during lulls in the stellar wind. When dayside winds are suppressed by the confining action of the stellar wind, nightside winds might pick up if there is sufficient horizontal transport of heat. A hot Jupiter loses mass at maximum rates of ~2 x 10^12 g/s during its host star's pre-main-sequence phase and ~2 x10^10 g/s during the star's main sequence lifetime, for total maximum losses of ~0.06% and ~0.6% of the planet's mass, respectively. For UV fluxes F_UV < 10^4 erg/cm^2/s, the mass loss rate is approximately energy-limited and is proportional to F_UV^0.9. For larger UV fluxes, such as those typical of T Tauri stars, radiative losses and plasma recombination force the mass loss rate to increase more slowly as F_UV^0.6. Dayside winds are quenched during the T Tauri phase because of confinement by overwhelming stellar wind pressure. We conclude that while UV radiation can indeed drive winds from hot Jupiters, such winds cannot significantly alter planetary masses during any evolutionary stage. They can, however, produce observable signatures. Candidates for explaining why the Lyman-alpha photons of HD 209458 are absorbed at Doppler-shifted velocities of +/- 100 km/s include charge-exchange in the shock between the planetary and stellar winds.
[57]  oai:arXiv.org:0704.3831  [pdf] - 847
The Origin of the Young Stars in the Nucleus of M31
Comments: 10 pages, 7 figures, accepted by ApJ, changes made from referee suggestions
Submitted: 2007-04-29, last modified: 2007-06-21
The triple nucleus of M31 consists of a population of old red stars in an eccentric disk (P1 and P2) and another population of younger A stars in a circular disk (P3) around M31's central supermassive black hole (SMBH). We argue that P1 and P2 determine the maximal radial extent of the younger A star population and provide the gas that fueled the starburst that generated P3. The eccentric stellar disk creates an $m=1$ non-axisymmetric perturbation to the potential. This perturbed potential drives gas into the inner parsec around the SMBH, if the pattern speed of the eccentric stellar disk is $\Omega_p \lesssim 3-10 {\rm km s^{-1} pc^{-1}}$. We show that stellar mass loss from P1 and P2 is sufficient to create a gravitationally unstable gaseous disk of $\sim 10^5\Msun$ every $0.1-1$ Gyrs, consistent with the 200 Myr age of P3. Similar processes may act in other systems to produce very compact nuclear starbursts.
[58]  oai:arXiv.org:0706.1241  [pdf] - 1000416
Inside-Out Evacuation of Transitional Protoplanetary Disks by the Magneto-Rotational Instability
Comments: Accepted to Nature Physics June 7, 2007. The manuscript for publication is embargoed per Nature policy. This arxiv.org version contains more technical details and discussion, and is distributed with permission from the editors. 10 pages, 4 figures
Submitted: 2007-06-10
How do T Tauri disks accrete? The magneto-rotational instability (MRI) supplies one means, but protoplanetary disk gas is typically too poorly ionized to be magnetically active. Here we show that the MRI can, in fact, explain observed accretion rates for the sub-class of T Tauri disks known as transitional systems. Transitional disks are swept clean of dust inside rim radii of ~10 AU. Stellar coronal X-rays ionize material in the disk rim, activating the MRI there. Gas flows from the rim to the star, at a rate limited by the depth to which X-rays ionize the rim wall. The wider the rim, the larger the surface area that the rim wall exposes to X-rays, and the greater the accretion rate. Interior to the rim, the MRI continues to transport gas; the MRI is sustained even at the disk midplane by super-keV X-rays that Compton scatter down from the disk surface. Accretion is therefore steady inside the rim. Blown out by radiation pressure, dust largely fails to accrete with gas. Contrary to what is usually assumed, ambipolar diffusion, not Ohmic dissipation, limits how much gas is MRI-active. We infer values for the transport parameter alpha on the order of 0.01 for GM Aur, TW Hyd, and DM Tau. Because the MRI can only afflict a finite radial column of gas at the rim, disk properties inside the rim are insensitive to those outside. Thus our picture provides one robust setting for planet-disk interaction: a protoplanet interior to the rim will interact with gas whose density, temperature, and transport properties are definite and decoupled from uncertain initial conditions. Our study also supplies half the answer to how disks dissipate: the inner disk drains from the inside out by the MRI, while the outer disk photoevaporates by stellar ultraviolet radiation.
[59]  oai:arXiv.org:astro-ph/0601654  [pdf] - 79512
A Brief History of Trans-Neptunian Space
Comments: Final proofed version for Protostars and Planets V; some numbers adjusted by factors of 2; references updated
Submitted: 2006-01-27, last modified: 2006-08-21
The Edgeworth-Kuiper belt encodes the dynamical history of the outer solar system. Kuiper belt objects (KBOs) bear witness to coagulation physics, the evolution of planetary orbits, and external perturbations from the solar neighborhood. We critically review the present-day belt's observed properties and the theories designed to explain them. Theories are organized according to a possible time-line of events. In chronological order, epochs described include (1) coagulation of KBOs in a dynamically cold disk, (2) formation of binary KBOs by fragmentary collisions and gravitational captures, (3) stirring of KBOs by Neptune-mass planets (``oligarchs''), (4) eviction of excess oligarchs, (5) continued stirring of KBOs by remaining planets whose orbits circularize by dynamical friction, (6) planetary migration and capture of Resonant KBOs, (7) creation of the inner Oort cloud by passing stars in an open stellar cluster, and (8) collisional comminution of the smallest KBOs. Recent work underscores how small, collisional, primordial planetesimals having low velocity dispersion permit the rapid assembly of ~5 Neptune-mass oligarchs at distances of 15-25 AU. We explore the consequences of such a picture. We propose that Neptune-mass planets whose orbits cross into the Kuiper belt for up to ~20 Myr help generate the high-perihelion members of the hot Classical disk and Scattered belt. By contrast, raising perihelia by sweeping secular resonances during Neptune's migration might fill these reservoirs too inefficiently when account is made of how little primordial mass might reside in bodies having sizes of order 100 km. These and other frontier issues in trans-Neptunian space are discussed quantitatively.
[60]  oai:arXiv.org:astro-ph/0607203  [pdf] - 83406
Brownian Motion in Planetary Migration
Comments: Accepted to ApJ
Submitted: 2006-07-10
A residual planetesimal disk of mass 10-100 Earth masses remained in the outer solar system following the birth of the giant planets, as implied by the existence of the Oort cloud, coagulation requirements for Pluto, and inefficiencies in planet formation. Upon gravitationally scattering planetesimal debris, planets migrate. Orbital migration can lead to resonance capture, as evidenced here in the Kuiper and asteroid belts, and abroad in extra-solar systems. Finite sizes of planetesimals render migration stochastic ("noisy"). At fixed disk mass, larger (fewer) planetesimals generate more noise. Extreme noise defeats resonance capture. We employ order-of-magnitude physics to construct an analytic theory for how a planet's orbital semi-major axis fluctuates in response to random planetesimal scatterings. To retain a body in resonance, the planet's semi-major axis must not random walk a distance greater than the resonant libration width. We translate this criterion into an analytic formula for the retention efficiency of the resonance as a function of system parameters, including planetesimal size. We verify our results with tailored numerical simulations. Application of our theory reveals that capture of Resonant Kuiper belt objects by a migrating Neptune remains effective if the bulk of the primordial disk was locked in bodies having sizes < O(100) km and if the fraction of disk mass in objects with sizes > 1000 km was less than a few percent. Coagulation simulations produce a size distribution of primordial planetesimals that easily satisfies these constraints. We conclude that stochasticity did not interfere with, nor modify in any substantive way, Neptune's ability to capture and retain Resonant Kuiper belt objects during its migration.
[61]  oai:arXiv.org:astro-ph/0410086  [pdf] - 67936
A Signature of Planetary Migration: The Origin of Asymmetric Capture in the 2:1 Resonance
Comments: 16 pages, 17 figures, accepted to ApJ
Submitted: 2004-10-04
The spatial distribution of Kuiper belt objects (KBOs) in 2:1 exterior resonance with Neptune constrains that planet's migration history. Numerical simulations demonstrate that fast planetary migration generates a larger population of KBOs trailing rather than leading Neptune in orbital longitude. This asymmetry corresponds to a greater proportion of objects caught into asymmetric resonance such that their resonance angles librate about values > pi (trailing) as opposed to < pi (leading). We provide, for the first time, an explanation of this phenomenon, using physical, analytic, and semi-analytic arguments. Central to our understanding is how planetary migration shifts the equilibrium points of the superposed direct and indirect potentials. Large capture asymmetries appear for exponential timescales of migration shorter than 10 Myr. The observed distribution of 2:1 KBOs (2 trailing and 7 leading) excludes migration timescales < 1 Myr with 99.65% confidence.
[62]  oai:arXiv.org:astro-ph/0312515  [pdf] - 61728
The Circumbinary Ring of KH 15D
Comments: Final version, ApJ, v607, 913 (June 1); includes prediction for full spectral energy distribution (new Figure 5)
Submitted: 2003-12-19, last modified: 2004-04-29
The light curves of the pre-main-sequence star KH 15D from the years 1913--2003 can be understood if the star is a member of an eccentric binary that is encircled by a vertically thin, inclined ring of dusty gas. Eclipses occur whenever the reflex motion of a star carries it behind the circumbinary ring; the eclipses occur with period equal to the binary orbital period of 48.4 days. Features of the light curve--including the amplitude of central reversals during mid-eclipse, the phase of eclipse with respect to the binary orbit phase, the level of brightness out-of-eclipse, the depth of eclipse, and the eclipse duty cycle--are all modulated on the timescale of nodal regression of the obscuring ring, in accord with the historical data. The ring has a mean radius near 3 AU and a radial width that is likely less than this value. While the inner boundary could be shepherded by the central binary, the outer boundary may require an exterior planet to confine it against viscous spreading. The ring must be vertically warped to maintain a non-zero inclination. Thermal pressure gradients and/or ring self-gravity can readily enforce rigid precession. In coming years, as the node of the ring regresses out of our line-of-sight toward the binary, the light curve from the system should cycle approximately back through its previous behavior. Near-term observations should seek to detect a mid-infrared excess from this system; we estimate the flux densities from the ring to be 3 mJy at wavelengths of 10--100 microns.