Showing papers in "arXiv: Earth and Planetary Astrophysics in 2012"

The physics of wind-blown sand and dust

Abstract: The transport of sand and dust by wind is a potent erosional force, creates sand dunes and ripples, and loads the atmosphere with suspended dust aerosols This article presents an extensive review of the physics of wind-blown sand and dust on Earth and Mars Specifically, we review the physics of aeolian saltation, the formation and development of sand dunes and ripples, the physics of dust aerosol emission, the weather phenomena that trigger dust storms, and the lifting of dust by dust devils and other small-scale vortices We also discuss the physics of wind-blown sand and dune formation on Venus and Titan

Planetary Candidates Observed by Kepler, III: Analysis of the First 16 Months of Data

Abstract: New transiting planet candidates are identified in sixteen months (May 2009 - September 2010) of data from the Kepler spacecraft. Nearly five thousand periodic transit-like signals are vetted against astrophysical and instrumental false positives yielding 1,091 viable new planet candidates, bringing the total count up to over 2,300. Improved vetting metrics are employed, contributing to higher catalog reliability. Most notable is the noise-weighted robust averaging of multi-quarter photo-center offsets derived from difference image analysis which identifies likely background eclipsing binaries. Twenty-two months of photometry are used for the purpose of characterizing each of the new candidates. Ephemerides (transit epoch, T_0, and orbital period, P) are tabulated as well as the products of light curve modeling: reduced radius (Rp/R*), reduced semi-major axis (d/R*), and impact parameter (b). The largest fractional increases are seen for the smallest planet candidates (197% for candidates smaller than 2Re compared to 52% for candidates larger than 2Re) and those at longer orbital periods (123% for candidates outside of 50-day orbits versus 85% for candidates inside of 50-day orbits). The gains are larger than expected from increasing the observing window from thirteen months (Quarter 1-- Quarter 5) to sixteen months (Quarter 1 -- Quarter 6). This demonstrates the benefit of continued development of pipeline analysis software. The fraction of all host stars with multiple candidates has grown from 17% to 20%, and the paucity of short-period giant planets in multiple systems is still evident. The progression toward smaller planets at longer orbital periods with each new catalog release suggests that Earth-size planets in the Habitable Zone are forthcoming if, indeed, such planets are abundant.

Architecture of Kepler's Multi-transiting Systems: II. New investigations with twice as many candidates

Abstract: We report on the orbital architectures of Kepler systems having multiple planet candidates identified in the analysis of data from the first six quarters of Kepler data and reported by Batalha et al (2013) These data show 899 transiting planet candidates in 365 multiple-planet systems and provide a powerful means to study the statistical properties of planetary systems Using a generic mass-radius relationship, we find that only two pairs of planets in these candidate systems (out of 761 pairs total) appear to be on Hill-unstable orbits, indicating ~96% of the candidate planetary systems are correctly interpreted as true systems We find that planet pairs show little statistical preference to be near mean-motion resonances We identify an asymmetry in the distribution of period ratios near first-order resonances (eg, 2:1, 3:2), with an excess of planet pairs lying wide of resonance and relatively few lying narrow of resonance Finally, based upon the transit duration ratios of adjacent planets in each system, we find that the interior planet tends to have a smaller transit impact parameter than the exterior planet does This finding suggests that the mode of the mutual inclinations of planetary orbital planes is in the range 10-22 degrees, for the packed systems of small planets probed by these observations

Rapid growth of gas-giant cores by pebble accretion

Abstract: The observed lifetimes of gaseous protoplanetary discs place strong constraints on gas and ice giant formation in the core accretion scenario. The approximately 10-Earth-mass solid core responsible for the attraction of the gaseous envelope has to form before gas dissipation in the protoplanetary disc is completed within 1-10 million years. Building up the core by collisions between km-sized planetesimals fails to meet this time-scale constraint, especially at wide stellar separations. Nonetheless, gas-giant planets are detected by direct imaging at wide orbital distances. In this paper, we numerically study the growth of cores by the accretion of cm-sized pebbles loosely coupled to the gas. We measure the accretion rate onto seed masses ranging from a large planetesimal to a fully grown 10-Earth-mass core and test different particle sizes. The numerical results are in good agreement with our analytic expressions, indicating the existence of two accretion regimes, one set by the azimuthal and radial particle drift for the lower seed masses and the other, for higher masses, by the velocity at the edge of the Hill sphere. In the former, the optimally accreted particle size increases with core mass, while in the latter the optimal size is centimeters, independent of core mass. We discuss the implications for rapid core growth of gas-giant and ice-giant cores. We conclude that pebble accretion can resolve the long-standing core accretion time-scale conflict. This requires a near-unity dust-to-gas ratio in the midplane, particle growth to mm and cm and the formation of massive planetesimals or low radial pressure support. The core growth time-scale is shortened by a factor 30-1,000 at 5 AU and by a factor 100-10,000 at 50 AU, compared to the gravitationally focused accretion of, respectively, low-scale-height planetesimal fragments or standard km-sized planetesimals.

Building Terrestrial Planets

Abstract: This paper reviews our current understanding of terrestrial planets formation. The focus is on computer simulations of the dynamical aspects of the accretion process. Throughout the chapter, we combine the results of these theoretical models with geochemical, cosmochemical and chronological constraints, in order to outline a comprehensive scenario of the early evolution of our Solar System. Given that the giant planets formed first in the protoplanetary disk, we stress the sensitive dependence of the terrestrial planet accretion process on the orbital architecture of the giant planets and on their evolution. This suggests a great diversity among the terrestrial planets populations in extrasolar systems. Issues such as the cause for the different masses and accretion timescales between Mars and the Earth and the origin of water (and other volatiles) on our planet are discussed at depth.

A simple model for the evolution of the dust population in protoplanetary disks

Abstract: Context: The global size and spatial distribution of dust is an important ingredient in the structure and evolution of protoplanetary disks and in the formation of larger bodies, such as planetesimals. Aims: We aim to derive simple equations that explain the global evolution of the dust surface density profile and the upper limit of the grain size distribution and which can readily be used for further modeling or for interpreting of observational data. Methods: We have developed a simple model that follows the upper end of the dust size distribution and the evolution of the dust surface density profile. This model is calibrated with state-of-the-art simulations of dust evolution, which treat dust growth, fragmentation, and transport in viscously evolving gas disks. Results: We find very good agreement between the full dust-evolution code and the toy model presented in this paper. We derive analytical profiles that describe the dust-to-gas ratios and the dust surface density profiles well in protoplanetary disks, as well as the radial flux by solid material "rain out", which is crucial for triggering any gravity assisted formation of planetesimals. We show that fragmentation is the dominating effect in the inner regions of the disk leading to a dust surface density exponent of -1.5, while the outer regions at later times can become drift-dominated, yielding a dust surface density exponent of -0.75. Our results show that radial drift is not efficient in fragmenting dust grains. This supports the theory that small dust grains are resupplied by fragmentation due to the turbulent state of the disk.

3.6 and 4.5 Micron Phase Curves and Evidence for Non-Equilibrium Chemistry in the Atmosphere of Extrasolar Planet HD 189733b

Abstract: We present new, full-orbit observations of the infrared phase variations of the canonical hot Jupiter HD 189733b obtained in the 3.6 and 4.5 micron bands using the Spitzer Space Telescope. When combined with previous phase curve observations at 8.0 and 24 micron, these data allow us to characterize the exoplanet's emission spectrum as a function of planetary longitude. We utilize improved methods for removing the effects of intrapixel sensitivity variations and accounting for the presence of time-correlated noise in our data. We measure a phase curve amplitude of 0.1242% +/- 0.0061% in the 3.6 micron band and 0.0982% +/- 0.0089% in the 4.5 micron band. We find that the times of minimum and maximum flux occur several hours earlier than predicted for an atmosphere in radiative equilibrium, consistent with the eastward advection of gas by an equatorial super-rotating jet. The locations of the flux minima in our new data differ from our previous observations at 8 micron, and we present new evidence indicating that the flux minimum observed in the 8 micron is likely caused by an over-shooting effect in the 8 micron array. We obtain improved estimates for HD 189733b's dayside planet-star flux ratio of 0.1466% +/- 0.0040% at 3.6 micron and 0.1787% +/- 0.0038% at 4.5 micron; these are the most accurate secondary eclipse depths obtained to date for an extrasolar planet. We compare our new dayside and nightside spectra for HD 189733b to the predictions of models from Burrows et al. (2008) and Showman et al. (2009). We find that HD 189733b's 4.5 micron nightside flux is 3.3 sigma smaller than predicted by the Showman et al. models, which assume that the chemistry is in local thermal equilibrium. We conclude that this discrepancy is best-explained by vertical mixing, which should lead to an excess of CO and correspondingly enhanced 4.5 micron absorption in this region. [abridged]

KELT-1b: A Strongly Irradiated, Highly Inflated, Short Period, 27 Jupiter-mass Companion Transiting a mid-F Star

Abstract: We present the discovery of KELT-1b, the first transiting low-mass companion from the wide-field Kilodegree Extremely Little Telescope-North (KELT-North) survey. The V=10.7 primary is a mildly evolved, solar-metallicity, mid-F star. The companion is a low-mass brown dwarf or super-massive planet with mass of 27.23+/-0.50 MJ and radius of 1.110+0.037-0.024 RJ, on a very short period (P=1.21750007) circular orbit. KELT-1b receives a large amount of stellar insolation, with an equilibrium temperature assuming zero albedo and perfect redistribution of 2422 K. Upper limits on the secondary eclipse depth indicate that either the companion must have a non-zero albedo, or it must experience some energy redistribution. Comparison with standard evolutionary models for brown dwarfs suggests that the radius of KELT-1b is significantly inflated. Adaptive optics imaging reveals a candidate stellar companion to KELT-1, which is consistent with an M dwarf if bound. The projected spin-orbit alignment angle is consistent with zero stellar obliquity, and the vsini of the primary is consistent with tidal synchronization. Given the extreme parameters of the KELT-1 system, we expect it to provide an important testbed for theories of the emplacement and evolution of short-period companions, and theories of tidal dissipation and irradiated brown dwarf atmospheres.

Planetary evaporation by UV & X-ray radiation: basic hydrodynamics

Abstract: We consider the evaporation of close in planets by the star's intrinsic EUV and X-ray radiation. We calculate evaporation rates by solving the hydrodynamical problem for planetary evaporation including heating from both X-ray and EUV radiation. We show that most close-in planets (a<0.1 AU) are evaporating hydrodynamically, with the evaporation occurring in two distinct regimes: X-ray driven, in which the X-ray heated flow contains a sonic point, and EUV driven, in which the X-ray region is entirely sub-sonic. The mass-loss rates scale as L_X/a^2 for X-ray driven evaporation, and as Phi_*^{1/2}/a for EUV driven evaporation at early times, with mass-loss rates of order 10e10-10e14 g/s. No exact scaling exists for the mass-loss rate with planet mass and planet radius, however, in general evaporation proceeds more rapidly for planets with lower densities and higher masses. Furthermore, we find that in general the transition from X-ray driven to EUV driven evaporation occurs at lower X-ray luminosities for planets closer to their parent stars and for planets with lower densities. Coupling our evaporation models to the evolution of the high energy radiation - which falls with time - we are able to follow the evolution of evaporating planets. We find that most planets start off evaporating in the X-ray driven regime, but switch to EUV driven once the X-ray luminosity falls below a critical value. The evolution models suggest that while `hot Jupiters' are evaporating, they are not evaporating at a rate sufficient to remove the entire gaseous envelope on Gyr time-scales. However, we do find that close in Neptune mass planets are more susceptible to complete evaporation of their envelopes. Thus we conclude that planetary evaporation is more important for lower mass planets, particularly those in the `hot Neptune'/`super Earth' regime.

Statistical Study of the Early Solar System's Instability with 4, 5 and 6 Giant Planets

Abstract: Several properties of the Solar System, including the wide radial spacing and orbital eccentricities of giant planets, can be explained if the early Solar System evolved through a dynamical instability followed by migration of planets in the planetesimal disk. Here we report the results of a statistical study, in which we performed nearly 10^4 numerical simulations of planetary instability starting from hundreds of different initial conditions. We found that the dynamical evolution is typically too violent, if Jupiter and Saturn start in the 3:2 resonance, leading to ejection of at least one ice giant from the Solar System. Planet ejection can be avoided if the mass of the transplanetary disk of planetesimals was large (M_disk>50 M_Earth), but we found that a massive disk would lead to excessive dynamical damping (e.g., final e_55 < 0.01 compared to present e_55=0.044, where e_55 is the amplitude of the fifth eccentric mode in the Jupiter's orbit), and to smooth migration that violates constraints from the survival of the terrestrial planets. Better results were obtained when the Solar System was assumed to have five giant planets initially and one ice giant, with the mass comparable to that of Uranus and Neptune, was ejected into interstellar space by Jupiter. The best results were obtained when the ejected planet was placed into the external 3:2 or 4:3 resonance with Saturn and M_disk ~ 20 M_Earth. The range of possible outcomes is rather broad in this case, indicating that the present Solar System is neither a typical nor expected result for a given initial state, and occurs, in best cases, with only a ~5% probability (as defined by the success criteria described in the main text). The case with six giant planets shows interesting dynamics but does offer significant advantages relative to the five planet case.

The Derivation, Properties and Value of Kepler's Combined Differential Photometric Precision

Abstract: The Kepler Mission is searching for Earth-size planets orbiting solar-like stars by simultaneously observing >160,000 stars to detect sequences of transit events in the photometric light curves. The Combined Differential Photometric Precision (CDPP) is the metric that defines the ease with which these weak terrestrial transit signatures can be detected. An understanding of CDPP is invaluable for evaluating the completeness of the Kepler survey and inferring the underlying planet population. This paper describes how the Kepler CDPP is calculated, and introduces tables of rms CDPP on a per-target basis for 3-, 6-, and 12-hour transit durations, which are now available for all Kepler observations. Quarter 3 is the first typical set of observations at the nominal length and completeness for a quarter, from 2009 September 18 to 2009 December 16, and we examine the properties of the rms CDPP distribution for this data set. Finally, we describe how to employ CDPP to calculate target completeness, an important use case.

Atmospheric Retrieval for Super-Earths: Uniquely Constraining the Atmospheric Composition with Transmission Spectroscopy

Abstract: We present a retrieval method based on Bayesian analysis to infer the atmospheric compositions and surface or cloud-top pressures from transmission spectra of exoplanets with general compositions. In this study, we identify what can unambiguously be determined about the atmospheres of exoplanets from their transmission spectra by applying the retrieval method to synthetic observations of the super-Earth GJ 1214b. Our approach to infer constraints on atmospheric parameters is to compute their joint and marginal posterior probability distributions using the MCMC technique in a parallel tempering scheme. A new atmospheric parameterization is introduced that is applicable to general atmospheres in which the main constituent is not known a priori and clouds may be present. Our main finding is that a unique constraint of the mixing ratios of the absorbers and up to two spectrally inactive gases (such as N2 and primordial H2+He) is possible if the observations are sufficient to quantify both (1) the broadband transit depths in at least one absorption feature for each absorber and (2) the slope and strength of the molecular Rayleigh scattering signature. The surface or cloud-top pressure can be quantified if a surface or cloud deck is present. The mean molecular mass can be constrained from the Rayleigh slope or the shapes of absorption features, thus enabling to distinguish between cloudy hydrogen-rich atmospheres and high mean molecular mass atmospheres. We conclude, however, that without the signature of Rayleigh scattering--even with robustly detected infrared absorption features--there is no reliable way to tell if the absorber is the main constituent of the atmosphere or just a minor species with a mixing ratio of <0.1%. The retrieval method leads us to a conceptual picture of which details in transmission spectra are essential for unique characterizations of well-mixed atmospheres.

How Thermal Evolution and Mass Loss Sculpt Populations of Super-Earths and Sub-Neptunes: Application to the Kepler-11 System and Beyond

Abstract: We use models of thermal evolution and XUV-driven mass loss to explore the composition and history of low-mass low-density transiting planets. We investigate the Kepler-11 system in detail and provide estimates of both the current and past planetary compositions. We find that a H/He atmosphere on Kepler-11b is highly vulnerable to mass loss. By comparing to formation models, we show that in situ formation of the system is unlikely. Instead we propose that it is a water-rich system of sub-Neptunes that migrated from beyond the snow line. For the broader population of observed planets, we show that there is a threshold in bulk planet density and incident flux above which no low-mass transiting planets have been observed. We suggest that this threshold is due to the instability of H/He atmospheres to XUV-driven mass loss. Importantly, we find that this flux-density threshold is well reproduced by our thermal evolution/contraction models that incorporate a standard mass loss prescription. Treating the planets' contraction history is essential because the planets have significantly larger radii during the early era of high XUV fluxes. Over time low mass planets with H/He envelopes can be transformed into water-dominated worlds with steam atmospheres or rocky super-Earths. Finally, we use this threshold to provide likely minimum masses and radial velocity amplitudes for the general population of Kepler candidates. Likewise, we use this threshold to provide constraints on the maximum radii of low-mass planets found by radial velocity surveys.

Transit Timing Observations from Kepler: IV. Confirmation of 4 Multiple Planet Systems by Simple Physical Models

Abstract: Eighty planetary systems of two or more planets are known to orbit stars other than the Sun. For most, the data can be sufficiently explained by non-interacting Keplerian orbits, so the dynamical interactions of these systems have not been observed. Here we present 4 sets of lightcurves from the Kepler spacecraft, which each show multiple planets transiting the same star. Departure of the timing of these transits from strict periodicity indicates the planets are perturbing each other: the observed timing variations match the forcing frequency of the other planet. This confirms that these objects are in the same system. Next we limit their masses to the planetary regime by requiring the system remain stable for astronomical timescales. Finally, we report dynamical fits to the transit times, yielding possible values for the planets' masses and eccentricities. As the timespan of timing data increases, dynamical fits may allow detailed constraints on the systems' architectures, even in cases for which high-precision Doppler follow-up is impractical.

A new vision on giant planet interiors: the impact of double diffusive convection

Abstract: While conventional interior models for Jupiter and Saturn are based on the simplistic assumption of a solid core surrounded by a homogeneous gaseous envelope, we derive new models with an inhomogeneous distribution of heavy elements, i.e. a gradient of composition, within these planets. Such a compositional stratification hampers large scale convection which turns into double-diffusive convection, yielding an inner thermal profile which departs from the traditionally assumed adiabatic interior, affecting these planet heat content and cooling history. To address this problem, we develop an analytical approach of layered double-diffusive convection and apply this formalism to Solar System gaseous giant planet interiors. These models satisfy all observational constraints and yield a metal enrichment for our gaseous giants up to 30 to 60% larger than previously thought. The models also constrain the size of the convective layers within the planets. As the heavy elements tend to be redistributed within the gaseous envelope, the models predict smaller than usual central cores inside Saturn and Jupiter, with possibly no core for this latter. These models open a new window and raise new challenges on our understanding of the internal structure of giant (solar and extrasolar) planets, in particular on the determination of their heavy material content, a key diagnostic for planet formation theories.

The HARPS-TERRA project I. Description of the algorithms, performance and new measurements on a few remarkable stars observed by HARPS

Abstract: Doppler spectroscopy has uncovered or confirmed all the known planets orbiting nearby stars. Two main techniques are used to obtain precision Doppler measurements at optical wavelengths. The first approach is the gas cell method, which consists on the least-squares matching of the spectrum of Iodine imprinted on the spectrum of the star. The second method relies on the construction of a stabilized spectrograph externally calibrated in wavelength. The most precise stabilized spectrometer in operation is HARPS, operated by ESO in La Silla Observatory/Chile. The Doppler measurements obtained with HARPS are typically obtained using the Cross-Correlation Function technique (CCF). It consists of multiplying the stellar spectrum with a weighted binary mask and finding the minimum of such product as a function of the Doppler shift. It is known that CCF is suboptimal in exploiting the Doppler information in the stellar spectrum. Here, we describe an algorithm to obtain precision RV measurements using least-squares matching of each observed spectrum to a high signal-to-noise ratio template derived from the same observations. Such algorithm is implemented in our software called HARPS-TERRA (Template Enhanced Radial velocity Re-analysis Application). New radial velocity measurements on a representative sample of stars observed by HARPS is used to illustrate the benefits of the proposed method. We show that, compared to CCF, template matching provides a significant improvement in accuracy, specially when applied to M dwarfs.

SOPHIE velocimetry of Kepler transit candidates VII. A false-positive rate of 35% for Kepler close-in giant exoplanet candidates

Abstract: The false-positive probability (FPP) of Kepler transiting candidates is a key value for statistical studies of candidate properties. A previous investigation of the stellar population in the Kepler field has provided an estimate for the FPP of less than 5% for most of the candidates. We report here the results of our radial velocity observations on a sample of 46 Kepler candidates with a transit depth greater than 0.4%, orbital period less than 25 days and host star brighter than Kepler magnitude 14.7. We used the SOPHIE spectrograph mounted on the 1.93-m telescope at the Observatoire de Haute-Provence to establish the nature of the transiting candidates. In this sample, we found five undiluted eclipsing binaries, two brown dwarfs, six diluted eclipsing binaries, and nine new transiting planets that complement the 11 already published planets. The remaining 13 candidates were not followed-up or remain unsolved due to photon noise limitation or lack of observations. From these results we computed the FPP for Kepler close-in giant candidates to be 34.8% \pm 6.5%. We aimed to investigate the variation of the FPP for giant candidates with the longer orbital periods and found that it should be constant for orbital periods between 10 and 200 days. This significant disagrees with the previous estimates. We discuss the reasons for this discrepancy and the possible extension of this work toward smaller planet candidates. Finally, taking the false-positive rate into account, we refined the occurrence rate of hot jupiters from the Kepler data.

First Light LBT AO Images of HR 8799 bcde at 1.65 and 3.3 Microns: New Discrepancies between Young Planets and Old Brown Dwarfs

Abstract: As the only directly imaged multiple planet system, HR 8799 provides a unique opportunity to study the physical properties of several planets in parallel. In this paper, we image all four of the HR 8799 planets at H-band and 3.3 microns with the new LBT adaptive optics system, PISCES, and LBTI/LMIRCam. Our images offer an unprecedented view of the system, allowing us to obtain H and 3.3$ micron photometry of the innermost planet (for the first time) and put strong upper-limits on the presence of a hypothetical fifth companion. We find that all four planets are unexpectedly bright at 3.3 microns compared to the equilibrium chemistry models used for field brown dwarfs, which predict that planets should be faint at 3.3 microns due to CH4 opacity. We attempt to model the planets with thick-cloudy, non-equilibrium chemistry atmospheres, but find that removing CH4 to fit the 3.3 micron photometry increases the predicted L' (3.8 microns) flux enough that it is inconsistent with observations. In an effort to fit the SED of the HR 8799 planets, we construct mixtures of cloudy atmospheres, which are intended to represent planets covered by clouds of varying opacity. In this scenario, regions with low opacity look hot and bright, while regions with high opacity look faint, similar to the patchy cloud structures on Jupiter and L/T transition brown-dwarfs. Our mixed cloud models reproduce all of the available data, but self-consistent models are still necessary to demonstrate their viability.

MOA-2011-BLG-293Lb: A test of pure survey microlensing planet detections

Abstract: Because of the development of large-format, wide-field cameras, microlensing surveys are now able to monitor millions of stars with sufficient cadence to detect planets. These new discoveries will span the full range of significance levels including planetary signals too small to be distinguished from the noise. At present, we do not understand where the threshold is for detecting planets. MOA-2011-BLG-293Lb is the first planet to be published from the new surveys, and it also has substantial followup observations. This planet is robustly detected in survey+followup data (Delta chi^2 ~ 5400). The planet/host mass ratio is q=5.3+/- 0.2*10^{-3}. The best fit projected separation is s=0.548+/- 0.005 Einstein radii. However, due to the s-->s^{-1} degeneracy, projected separations of s^{-1} are only marginally disfavored at Delta chi^2=3. A Bayesian estimate of the host mass gives M_L = 0.43^{+0.27}_{-0.17} M_Sun, with a sharp upper limit of M_L < 1.2 M_Sun from upper limits on the lens flux. Hence, the planet mass is m_p=2.4^{+1.5}_{-0.9} M_Jup, and the physical projected separation is either r_perp = ~1.0 AU or r_perp = ~3.4 AU. We show that survey data alone predict this solution and are able to characterize the planet, but the Delta chi^2 is much smaller (Delta chi^2~500) than with the followup data. The Delta chi^2 for the survey data alone is smaller than for any other securely detected planet. This event suggests a means to probe the detection threshold, by analyzing a large sample of events like MOA-2011-BLG-293, which have both followup data and high cadence survey data, to provide a guide for the interpretation of pure survey microlensing data.

Planet Hunters: A Transiting Circumbinary Planet in a Quadruple Star System

Abstract: We report the discovery and confirmation of a transiting circumbinary planet (PH1b) around KIC 4862625, an eclipsing binary in the Kepler field. The planet was discovered by volunteers searching the first six Quarters of publicly available Kepler data as part of the Planet Hunters citizen science project. Transits of the planet across the larger and brighter of the eclipsing stars are detectable by visual inspection every ~137 days, with seven transits identified in Quarters 1-11. The physical and orbital parameters of both the host stars and planet were obtained via a photometric-dynamical model, simultaneously fitting both the measured radial velocities and the Kepler light curve of KIC 4862625. The 6.18 +/- 0.17 Earth radii planet orbits outside the 20-day orbit of an eclipsing binary consisting of an F dwarf (1.734 +/- 0.044 Solar radii, 1.528 +/- 0.087 Solar masses) and M dwarf (0.378+/- 0.023 Solar radii, 0.408 +/- 0.024 Solar masses). For the planet, we find an upper mass limit of 169 Earth masses (0.531 Jupiter masses) at the 99.7% confidence level. With a radius and mass less than that of Jupiter, PH1b is well within the planetary regime. Outside the planet's orbit, at ~1000 AU,a previously unknown visual binary has been identified that is likely bound to the planetary system, making this the first known case of a quadruple star system with a transiting planet.

Foretellings of Ragnar\"ok: World-engulfing Asymptotic Giants and the Inheritance of White Dwarfs

Abstract: The search for planets around White Dwarf stars, and evidence for dynamical instability around them in the form of atmospheric pollution and circumstellar discs, raises questions about the nature of planetary systems that can survive the vicissitudes of the Asymptotic Giant Branch (AGB). We study the competing effects, on planets at several AU from the star, of strong tidal forces arising from the star's large convective envelope, and of the planets' orbital expansion due to stellar mass loss. We, for the first time, study the evolution of planets while following each thermal pulse on the AGB. For Jovian planets, tidal forces are strong, and can pull into the envelope planets initially at ~3 AU for a 1M_Sol star and ~5 AU for a 5M_Sol star. Lower-mass planets feel weaker tidal forces, and Terrestrial planets initially within 1.5-3 AU enter the stellar envelope. Thus, low-mass planets that begin inside the maximum stellar radius can survive, as their orbits expand due to mass loss. The inclusion of a moderate planetary eccentricity slightly strengthens the tidal forces experienced by Jovian planets. Eccentric Terrestrial planets are more at risk, since their eccentricity does not decay and their small pericentre takes them inside the stellar envelope. We also find the closest radii at which planets will be found around White Dwarfs, assuming that any planet entering the stellar envelope is destroyed. Planets are in that case unlikely to be found inside ~1.5 AU of a White Dwarf with a 1M_Sol progenitor and ~10 AU of a White Dwarf with a 5M_Sol progenitor.

Temperature-Pressure Profile of the hot Jupiter HD 189733b from HST Sodium Observations: Detection of Upper Atmospheric Heating

Abstract: We present transmission spectra of the hot Jupiter HD 189733b taken with the Space Telescope Imaging Spectrograph aboard HST. The spectra cover the wavelength range 5808-6380 Ang with a resolving power of R=5000. We detect absorption from the NaI doublet within the exoplanet's atmosphere at the 9 sigma confidence level within a 5 Ang band (absorption depth 0.09 +/- 0.01%) and use the data to measure the doublet's spectral absorption profile. We detect only the narrow cores of the doublet. The narrowness of the feature could be due to an obscuring high-altitude haze of an unknown composition or a significantly sub-solar NaI abundance hiding the line wings beneath a H2 Rayleigh signature. We compare the spectral absorption profile over 5.5 scale heights with model spectral absorption profiles and constrain the temperature at different atmospheric regions, allowing us to construct a vertical temperature profile. We identify two temperature regimes; a 1280 +/- 240 K region derived from the NaI doublet line wings corresponding to altitudes below ~ 500 km, and a 2800 +/- 400 K region derived from the NaI doublet line cores corresponding to altitudes from ~ 500-4000 km. The zero altitude is defined by the white-light radius of Rp/Rstar=0.15628 +/- 0.00009. The temperature rises with altitude, which is likely evidence of a thermosphere. The absolute pressure scale depends on the species responsible for the Rayleigh signature and its abundance. We discuss a plausible scenario for this species, a high-altitude silicate haze, and the atmospheric temperature-pressure profile that results. In this case, the high altitude temperature rise for HD 189733b occurs at pressures of 10^-5 to 10^-8 bar.

They might be giants: luminosity class, planet frequency, and planet-metallicity relation of the coolest Kepler target stars

Abstract: We estimate the stellar parameters of late K and early M type Kepler target stars. We obtain medium resolution visible spectra of 382 stars with Kp-J>2 (~K5 and later spectral type). We determine luminosity class by comparing the strength of gravity-sensitive indices (CaH, K I, Ca II, and Na I) to their strength in a sample of stars of known luminosity class. We find that giants constitute 96+-% of the bright (Kp 14) stars, significantly higher than fractions based on the stellar parameters quoted in the Kepler Input Catalog (KIC). The KIC effective temperatures are systematically (110 +15 -35} K) higher than temperatures we determine from fitting our spectra to PHOENIX stellar models. Through Monte Carlo simulations of the Kepler exoplanet candidate population, we find a planet occurrence of 0.36+-0.08 when giant stars are properly removed, somewhat higher than when a KIC log(g)>4 criterion is used (0.27+-0.05). Lastly, we show that there is no significant difference in g-r color (a probe of metallicity) between late-type Kepler stars with transiting Earth-to-Neptune sized exoplanet candidates and dwarf stars with no detected transits. We show that a previous claimed offset between these two populations is most likely an artifact of including a large number of misidentified giants.

Orbital characterization of the \beta Pictoris b giant planet

Abstract: In June 2010, we confirmed the existence of a giant planet in the disk of the young star Beta Pictoris, located between 8 AU and 15 AU from the star. This young planet offers the rare opportunity to monitor a large fraction of the orbit using the imaging technique over a reasonably short timescale. Using the NAOS-CONICA adaptive-optics instrument (NACO) at the Very Large Telescope (VLT), we obtained repeated follow-up images of the Bpic system in the Ks and L' filters at four new epochs in 2010 and 2011. Complementing these data with previous measurements, we conduct a homogeneous analysis, which covers more than eight yrs, to accurately monitor the Bpic b position relative to the star. On the basis of the evolution of the planet's relative position with time, we derive the best-fit orbital solutions for our measurements. More reliable results are found with a Markov-chain Monte Carlo approach. The solutions favor a low-eccentricity orbit e < 0.17, with semi-major axis in the range 8--9 AU corresponding to orbital periods of 17--21 yrs. Our solutions favor a highly inclined solution with a peak around i=88.5+-1.7 deg, and a longitude of ascending node tightly constrained at Omega = -147.5+-1.5 deg. These results indicate that the orbital plane of the planet is likely to be above the midplane of the main disk, and compatible with the warp component of the disk being tilted between 3.5 deg and 4.0 deg. This suggests that the planet plays a key role in the origin of the inner warped-disk morphology of the Bpic disk. Finally, these orbital parameters are consistent with the hypothesis that the planet is responsible for the transit-like event observed in November 1981, and also linked to the cometary activity observed in the Bpic system.

Hint of a transiting extended atmosphere on 55 Cancri b

Abstract: The naked-eye star 55 Cancri hosts a planetary system with five known planets, including a hot super-Earth (55 Cnc e) extremely close to its star and a farther out giant planet (55 Cnc b), found in milder irradiation conditions with respect to other known hot Jupiters. This system raises important questions on the evolution of atmospheres for close-in exoplanets, and the dependence with planetary mass and irradiation. These questions can be addressed by Lyman-alpha transit observations of the extended hydrogen planetary atmospheres, complemented by contemporaneous measurements of the stellar X-ray flux. In fact, planet `e' has been detected in transit, suggesting the system is seen nearly edge-on. Yet, planet `b' has not been observed in transit so far. Here, we report on Hubble Space Telescope STIS Lyman-alpha and Chandra ACIS-S X-ray observations of 55 Cnc. These simultaneous observations cover two transits of 55 Cnc e and two inferior conjunctions of 55 Cnc b. They reveal the star as a bright Lyman-alpha target and a variable X-ray source. While no significant signal is detected during the transits of 55 Cnc e, we detect a surprising Lyman-alpha absorption of 7.5 +/- 1.8% (4.2 sigma) at inferior conjunctions of 55 Cnc b. The absorption is only detected over the range of Doppler velocities where the stellar radiation repels hydrogen atoms towards the observer. We calculate a false-alarm probability of 4.4%, which takes into account the a-priori unknown transit parameters. This result suggests the possibility that 55 Cnc b has an extended upper H I atmosphere, which undergoes partial transits when the planet grazes the stellar disc. If confirmed, it would show that planets cooler than hot Jupiters can also have extended atmospheres.

The Hunt for Exomoons with Kepler (HEK): I. Description of a New Observational Project

Abstract: Two decades ago, empirical evidence concerning the existence and frequency of planets around stars, other than our own, was absent. Since this time, the detection of extrasolar planets from Jupiter-sized to most recently Earth-sized worlds has blossomed and we are finally able to shed light on the plurality of Earth-like, habitable planets in the cosmos. Extrasolar moons may also be frequent habitable worlds but their detection or even systematic pursuit remains lacking in the current literature. Here, we present a description of the first systematic search for extrasolar moons as part of a new observational project called "The Hunt for Exomoons with Kepler" (HEK). The HEK project distills the entire list of known transiting planet candidates found by Kepler (2326 at the time of writing) down to the most promising candidates for hosting a moon. Selected targets are fitted using a multimodal nested sampling algorithm coupled with a planet-with-moon light curve modelling routine. By comparing the Bayesian evidence of a planet-only model to that of a planet-with-moon, the detection process is handled in a Bayesian framework. In the case of null detections, upper limits derived from posteriors marginalised over the entire prior volume will be provided to inform the frequency of large moons around viable planetary hosts, eta-moon. After discussing our methodologies for target selection, modelling, fitting and vetting, we provide two example analyses.

Transit Timing Observations from Kepler: III. Confirmation of 4 Multiple Planet Systems by a Fourier-Domain Study of Anti-correlated Transit Timing Variations

Abstract: We present a method to confirm the planetary nature of objects in systems with multiple transiting exoplanet candidates. This method involves a Fourier-Domain analysis of the deviations in the transit times from a constant period that result from dynamical interactions within the system. The combination of observed anti-correlations in the transit times and mass constraints from dynamical stability allow us to claim the discovery of four planetary systems Kepler-25, Kepler-26, Kepler-27, and Kepler-28, containing eight planets and one additional planet candidate.

Fast Rise of "Neptune-Size" Planets ($4-8 R_{\rm Earth}$) from $P\sim10$ to $\sim250$ days -- Statistics of Kepler Planet Candidates Up to $\sim 0.75 {\rm AU}$

Abstract: We infer the period ($P$) and size ($R_p$) distribution of Kepler transiting planet candidates with $R_p\ge 1 R_{\rm Earth}$ and $P 10$ days, the planet frequency d$N_p$/d$\log$P for "Neptune-size" planets ($R_p = 4-8 R_{\rm Earth}$) increases with period as $\propto P^{0.7\pm0.1}$. In contrast, d$N_p$/d$\log$P for "super-Earth-size" ($2-4 R_{\rm Earth}$) as well as "Earth-size" ($1-2 R_{\rm Earth}$) planets are consistent with a nearly flat distribution as a function of period ($\propto P^{0.11\pm0.05}$ and $\propto P^{-0.10\pm0.12}$, respectively), and the normalizations are remarkably similar (within a factor of $\sim 1.5$ at $50 $ days). Planet size distribution evolves with period, and generally the relative fractions for big planets ($\sim 3-10 R_{\rm Earth}$) increase with period. The shape of the distribution function is not sensitive to changes in selection criteria of the sample. The implied nearly flat or rising planet frequency at long period appears to be in tension with the sharp decline at $\sim 100$ days in planet frequency for low mass planets (planet mass $m_p < 30 M_{\rm Earth}$) recently suggested by HARPS survey. Within $250$ days, the cumulative frequencies for Earth-size and super-Earth-size planets are remarkably similar ($\sim 28 %$ and $25%$), while Neptune-size and Jupiter-size planets are $\sim 7%$, and $\sim 3%$, respectively. A major potential uncertainty arises from the unphysical impact parameter distribution of the candidates.

Characterization of exoplanets from their formation II: The planetary mass-radius relationship

Abstract: The research of exoplanets has entered an era in which we characterize extrasolar planets. This has become possible with measurements of radii and luminosities. Meanwhile, radial velocity surveys discover also very low-mass planets. Uniting all this observational data into one coherent picture to better understand planet formation is an important, but difficult undertaking. Our approach is to develop a model which can make testable predictions for all these observational methods. We continue to describe how we have extended our formation model into a self-consistently coupled formation and evolution model. We show how we calculate the internal structure of the solid core and radiogenic heating. We also improve the protoplanetary disk model. Finally, we conduct population synthesis calculations. We present how the planetary mass-radius relationship of planets with primordial H/He envelopes forms and evolves in time. The basic shape of the M-R relation can be understood from the core accretion model. Low-mass planets cannot bind massive envelopes, while super-critical cores necessarily trigger runway gas accretion, leading to "forbidden" zones in the M-R plane. For a given mass, there is a considerable diversity of radii. We compare the synthetic M-R relation with the observed one, finding good agreement for a>0.1 AU. The synthetic radius distribution is characterized by a strong increase towards small R, and a second, lower local maximum at ~1 Jovian radius. The increase towards small radii reflects the increase of the mass function towards low M. The second local maximum is due to the fact that radii are nearly independent of mass for giant planets. A comparison of the synthetic radius distribution with Kepler data shows agreement for R>2 Earth radii, but divergence for smaller radii. We predict that in the next few years, Kepler should find the second, local maximum at ~1 Jovian radius.

Strong tidal dissipation in Saturn and constraints on Enceladus' thermal state from astrometry

Abstract: Tidal interactions between Saturn and its satellites play a crucial role in both the orbital migration of the satellites and the heating of their interiors. Therefore constraining the tidal dissipation of Saturn (here the ratio k2/Q) opens the door to the past evolution of the whole system. If Saturn's tidal ratio can be determined at different frequencies, it may also be possible to constrain the giant planet's interior structure, which is still uncertain. Here, we try to determine Saturn's tidal ratio through its current effect on the orbits of the main moons, using astrometric data spanning more than a century. We find an intense tidal dissipation (k2/Q= (2.3 \pm 0.7) \times 10-4), which is about ten times higher than the usual value estimated from theoretical arguments. As a consequence, eccentricity equilibrium for Enceladus can now account for the huge heat emitted from Enceladus' south pole. Moreover, the measured k2/Q is found to be poorly sensitive to the tidal frequency, on the short frequency interval considered. This suggests that Saturn's dissipation may not be controlled by turbulent friction in the fluid envelope as commonly believed. If correct, the large tidal expansion of the moon orbits due to this strong Saturnian dissipation would be inconsistent with the moon formations 4.5 Byr ago above the synchronous orbit in the Saturnian subnebulae. But it would be compatible with a new model of satellite formation in which the Saturnian satellites formed possibly over longer time scale at the outer edge of the main rings. In an attempt to take into account for possible significant torques exerted by the rings on Mimas, we fitted a constant rate da/dt on Mimas semi-major axis, also. We obtained an unexpected large acceleration related to a negative value of da/dt= -(15.7 \pm 4.4) \times 10-15 au/day.