Circumbinary Planet Formation in the Kepler-16 System. II. A Toy Model for In-situ Planet Formation within a Debris Belt

Recent simulations have shown that the formation of planets in circumbinary configurations (such as those recently discovered by Kepler) is dramatically hindered at the planetesimal accretion stage. The combined action of the binary and the protoplanetary disk acts to raise impact velocities between km-sized planetesimals beyond their destruction threshold, halting planet formation within at least 10 AU from the binary. It has been proposed that a primordial population of "large" planetesimals (100 km or more in size), as produced by turbulent concentration mechanisms, would be able to bypass this bottleneck; however, it is not clear whether these processes are viable in the highly perturbed circumbinary environments. We perform two-dimensional hydrodynamical and N-body simulations to show that km-sized planetesimals and collisional debris can drift and be trapped in a belt close to the central binary. Within this belt, planetesimals could initially grow by accreting debris, ultimately becoming "indestructible" seeds that can accrete other planetesimals in-situ despite the large impact speeds. We find that large, indestructible planetesimals can be formed close to the central binary within $10^5$ years, therefore showing that even a primordial population of "small" planetesimals can feasibly form a planet.Comment: 3rd version, addressing referee report, typos & references fixed. 15 pages, 10 figures. Accepted for publication in Ap

Planet Formation in Circumbinary Configurations: Turbulence Inhibits Planetesimal Accretion

The existence of planets born in environments highly perturbed by a stellar companion represents a major challenge to the paradigm of planet formation. In numerical simulations, the presence of a close binary companion stirs up the relative velocity between planetesimals, which is fundamental in determining the balance between accretion and erosion. However, the recent discovery of circumbinary planets by Kepler establishes that planet formation in binary systems is clearly viable. We perform N-body simulations of planetesimals embedded in a protoplanetary disk, where planetesimal phasing is frustrated by the presence of stochastic torques, modeling the expected perturbations of turbulence driven by the magnetorotational instability (MRI). We examine perturbation amplitudes relevant to dead zones in the midplane (conducive to planet formation in single stars), and find that planetesimal accretion can be inhibited even in the outer disk (4-10 AU) far from the central binary, a location previously thought to be a plausible starting point for the formation of circumbinary planets.Comment: 5 pages, 4 figures. Accepted for publication in ApJ Letter

Circumbinary Planet Formation in the Kepler-16 system. I. N-body Simulations

The recently discovered circumbinary planets (Kepler-16 b, Kepler-34 b, Kepler-35 b) represent the first direct evidence of the viability of planet formation in circumbinary orbits. We report on the results of N-body simulations investigating planetesimal accretion in the Kepler-16 b system, focusing on the range of impact velocities under the influence of both stars' gravitational perturbation and friction from a putative protoplanetary disk. Our results show that planet formation might be effectively inhibited for a large range in semi-major axis (1.75 < a_P < 4 AU), suggesting that the planetary core must have either migrated from outside 4 AU, or formed in situ very close to its current location.Comment: 6 pages, 4 figures. Accepted for publication in ApJ. Additional materials available at http://www.stefanom.org/?kep16

Systemic: A Testbed For Characterizing the Detection of Extrasolar Planets. I. The Systemic Console Package

We present the systemic Console, a new all-in-one, general-purpose software package for the analysis and combined multiparameter fitting of Doppler radial velocity (RV) and transit timing observations. We give an overview of the computational algorithms implemented in the Console, and describe the tools offered for streamlining the characterization of planetary systems. We illustrate the capabilities of the package by analyzing an updated radial velocity data set for the HD128311 planetary system. HD128311 harbors a pair of planets that appear to be participating in a 2:1 mean motion resonance. We show that the dynamical configuration cannot be fully determined from the current data. We find that if a planetary system like HD128311 is found to undergo transits, then self-consistent Newtonian fits to combined radial velocity data and a small number of timing measurements of transit midpoints can provide an immediate and vastly improved characterization of the planet's dynamical state.Comment: 10 pages, 5 figures, accepted for publication on PASP. Additional material at http://www.ucolick.org/~smeschia/systemic.ph

The Lick-Carnegie Exoplanet Survey: A Saturn-Mass Planet in the Habitable Zone of the Nearby M4V Star HIP 57050

Precision radial velocities from Keck/HIRES reveal a Saturn-mass planet orbiting the nearby M4V star HIP 57050. The planet has a minimum mass of 0.3 Jupiter-mass, an orbital period of 41.4 days, and an orbital eccentricity of 0.31. V-band photometry reveals a clear stellar rotation signature of the host star with a period of 98 days, well separated from the period of the radial velocity variations and reinforcing a Keplerian origin for the observed velocity variations. The orbital period of this planet corresponds to an orbit in the habitable zone of HIP 57050, with an expected planetary temperature of approximately 230 K. The star has a metallicity of [Fe/H] = 0.32+/-0.06 dex, of order twice solar and among the highest metallicity stars in the immediate solar neighborhood. This newly discovered planet provides further support that the well-known planet-metallicity correlation for F, G, and K stars also extends down into the M-dwarf regime. The a priori geometric probability for transits of this planet is only about 1%. However, the expected eclipse depth is ~7%, considerably larger than that yet observed for any transiting planet. Though long on the odds, such a transit is worth pursuing as it would allow for high quality studies of the atmosphere via transmission spectroscopy with HST. At the expected planetary effective temperature, the atmosphere may contain water clouds.Comment: 20 pages, 5 figures, 3 tables, to appear in the May 20 issue of ApJ

A 4-Planet System Orbiting the K0V Star HD 141399

We present precision radial velocity (RV) data sets from Keck-HIRES and from Lick Observatory's new Automated Planet Finder Telescope and Levy Spectrometer on Mt. Hamilton that reveal a multiple-planet system orbiting the nearby, slightly evolved, K-type star HD 141399. Our 91 observations over 10.5 years suggest the presence of four planets with orbital periods of 94.35, 202.08, 1070.35, and 3717.35 days and minimum masses of 0.46, 1.36, 1.22, and 0.69 Jupiter masses respectively. The orbital eccentricities of the three inner planets are small, and the phase curves are well sampled. The inner two planets lie just outside the 2:1 resonance, suggesting that the system may have experienced dissipative evolution during the protoplanetary disk phase. The fourth companion is a Jupiter-like planet with a Jupiter-like orbital period. Its orbital eccentricity is consistent with zero, but more data will be required for an accurate eccentricity determination.Comment: 11 pages, 13 figures, To appear in the Astrophysical Journa

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

We investigate the possibility that the large orbital eccentricity of the transiting Neptune-mass planet Gliese 436b is maintained in the face of tidal dissipation by a second planet in the system. We find that the currently observed configuration can be understood if Gliese 436b and a putative companion have evolved to a quasi-stationary fixed point in which the planets' orbital apses are co-linear and in which secular variations in the orbital eccentricities of the two planets have been almost entirely damped out. We adopt an octopole-order secular theory based on a Legendre expansion in the semi-major axis ratio to delineate well-defined regions of (P_c, M_c, e_c) space that can be occupied by a perturbing companion. We incorporate the evolutionary effect of tidal dissipation into our secular model of the system, and solve the resulting initial value problems for a large sample of the allowed configurations. We then polish the stationary configurations derived from secular theory with full numerical integrations. We present our results in the form of candidate companion planets to Gliese 436b. For these candidates, radial velocity half-amplitudes, K_c, are of order 3 m/s, and the maximum amplitude of orbit-to-orbit transit timing variations are of order Delta t=1 s to Delta t=5s. For the particular example case of a perturber with orbital period, P_c=40 d, mass, M_c=8.5 M_Earth, and eccentricity, e_c=0.58, we confirm our semi-analytic calculations with a full numerical 3-body integration of the orbital decay that includes tidal damping and spin evolution. Additionally, we discuss the possibility of many-perturber stationary configurations, utilizing modified Laplace-Lagrange secular theory.Comment: 25 pages, 7 figures. Accepted for publication in the Astrophysical Journa

The Lick-Carnegie Survey: Four New Exoplanet Candidates

We present new precise HIRES radial velocity (RV) data sets of five nearby stars obtained at Keck Observatory. HD 31253, HD 218566, HD 177830, HD 99492 and HD 74156 are host stars of spectral classes F through K and show radial velocity variations consistent with new or additional planetary companions in Keplerian motion. The orbital parameters of the candidate planets in the five planetary systems span minimum masses of M sin i = 27.43 M_{earth} to M sin i = 8.28 M_{jup}, periods of 17.05 to 4696.95 days and eccentricities ranging from circular to extremely eccentric (e ~ 0.63). The 5th star, HD 74156, was known to have both a 52-day and a 2500-day planet, and was claimed to also harbor a 3rd planet at 336d, in apparent support of the "Packed Planetary System" hypothesis. Our greatly expanded data set for HD 74156 provides strong confirmation of both the 52-day and 2500-d planets, but strongly contradicts the existence of a 336-day planet, and offers no significant evidence for any other planets in the system.Comment: 13 pages, 15 figures. Accepted for publication in ApJ. Fixed typos in Table 2. Additional material at http://www.ucolick.org/~smeschia/4planet.ph