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

The Occurrence of Potentially Habitable Planets Orbiting M Dwarfs Estimated from the Full Kepler Dataset and an Empirical Measurement of the Detection Sensitivity

Abstract: We present an improved estimate of the occurrence rate of small planets orbiting small stars by searching the full four-year Kepler data set for transiting planets using our own planet detection pipeline and conducting transit injection and recovery simulations to empirically measure the search completeness of our pipeline We identified 156 planet candidates, including one object that was not previously identified as a Kepler Object of Interest We inspected all publicly available follow-up images, observing notes, and centroid analyses, and corrected for the likelihood of false positives We evaluated the sensitivity of our detection pipeline on a star-by-star basis by injecting 2000 transit signals into the light curve of each target star For periods shorter than 50 days, we find 056 (+006/-005) Earth-size planets (1-15 Earth radii) and 046 (+007/-005) super-Earths (15-2 Earth radii) per M dwarf In total, we estimate a cumulative planet occurrence rate of $25\pm02$ planets per M dwarf with radii 1-4 Earth radii and periods shorter than 200 days Within a conservatively defined habitable zone based on the moist greenhouse inner limit and maximum greenhouse outer limit, we estimate an occurrence rate of 016 (+017/-007) Earth-size planets and 012 (+010/-005) super-Earths per M dwarf habitable zone Adopting the broader insolation boundaries of the recent Venus and early Mars limits yields a higher estimate of 024 (+018/-008) Earth-size planets and 021 (+011/-006) super-Earths per M dwarf habitable zone This suggests that the nearest potentially habitable non-transiting and transiting Earth-size planets are $26\pm04$ pc and 106 (+16/-18) pc away, respectively If we include super-Earths, these distances diminish to $21\pm02$ pc and 86 (+07/-08) pc

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

Abstract: Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane 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.

Characterizing transiting exoplanet atmospheres with JWST

Abstract: We explore how well James Webb Space Telescope (JWST) spectra will likely constrain bulk atmospheric properties of transiting exoplanets. We start by modeling the atmospheres of archetypal hot Jupiter, warm Neptune, warm sub-Neptune, and cool super-Earth planets with clear, cloudy, or high mean molecular weight atmospheres. Next we simulate the $\lambda = 1 - 11$ $\mu$m transmission and emission spectra of these systems for several JWST instrument modes for single transit and eclipse events. We then perform retrievals to determine how well temperatures and molecular mixing ratios (CH$_4$, CO, CO$_2$, H$_2$O, NH$_3$) can be constrained. We find that $\lambda = 1 - 2.5$ $\mu$m transmission spectra will often constrain the major molecular constituents of clear solar composition atmospheres well. Cloudy or high mean molecular weight atmospheres will often require full $1 - 11$ $\mu$m spectra for good constraints, and emission data may be more useful in cases of sufficiently high $F_p$ and high $F_p/F_*$. Strong temperature inversions in the solar composition hot Jupiter atmosphere should be detectable with $1 - 2.5+$ $\mu$m emission spectra, and $1 - 5+$ $\mu$m emission spectra will constrain the temperature-pressure profiles of warm planets. Transmission spectra over $1 - 5+$ $\mu$m will constrain [Fe/H] values to better than 0.5 dex for the clear atmospheres of the hot and warm planets studied. Carbon-to-oxygen ratios can be constrained to better than a factor of 2 in some systems. We expect that these results will provide useful predictions of the scientific value of single event JWST spectra until its on-orbit performance is known.

Planetary Candidates Observed by Kepler V: Planet Sample from Q1-Q12 (36 Months)

Abstract: The Kepler mission discovered 2842 exoplanet candidates with 2 years of data. We provide updates to the Kepler planet candidate sample based upon 3 years (Q1-Q12) of data. Through a series of tests to exclude false-positives, primarily caused by eclipsing binary stars and instrumental systematics, 855 additional planetary candidates have been discovered, bringing the total number known to 3697. We provide revised transit parameters and accompanying posterior distributions based on a Markov Chain Monte Carlo algorithm for the cumulative catalogue of Kepler Objects of Interest. There are now 130 candidates in the cumulative catalogue that receive less than twice the flux the Earth receives and more than 1100 have a radius less than 1.5 Rearth. There are now a dozen candidates meeting both criteria, roughly doubling the number of candidate Earth analogs. A majority of planetary candidates have a high probability of being bonafide planets, however, there are populations of likely false-positives. We discuss and suggest additional cuts that can be easily applied to the catalogue to produce a set of planetary candidates with good fidelity. The full catalogue is publicly available at the NASA Exoplanet Archive.

Planetary Candidates Observed by Kepler. VII. The First Fully Uniform Catalog Based on The Entire 48 Month Dataset (Q1-Q17 DR24)

Abstract: We present the seventh Kepler planet candidate catalog, which is the first to be based on the entire, uniformly processed, 48 month Kepler dataset. This is the first fully automated catalog, employing robotic vetting procedures to uniformly evaluate every periodic signal detected by the Q1-Q17 Data Release 24 (DR24) Kepler pipeline. While we prioritize uniform vetting over the absolute correctness of individual objects, we find that our robotic vetting is overall comparable to, and in most cases is superior to, the human vetting procedures employed by past catalogs. This catalog is the first to utilize artificial transit injection to evaluate the performance of our vetting procedures and quantify potential biases, which are essential for accurate computation of planetary occurrence rates. With respect to the cumulative Kepler Object of Interest (KOI) catalog, we designate 1,478 new KOIs, of which 402 are dispositioned as planet candidates (PCs). Also, 237 KOIs dispositioned as false positives (FPs) in previous Kepler catalogs have their disposition changed to PC and 118 PCs have their disposition changed to FP. This brings the total number of known KOIs to 8,826 and PCs to 4,696. We compare the Q1-Q17 DR24 KOI catalog to previous KOI catalogs, as well as ancillary Kepler catalogs, finding good agreement between them. We highlight new PCs that are both potentially rocky and potentially in the habitable zone of their host stars, many of which orbit solar-type stars. This work represents significant progress in accurately determining the fraction of Earth-size planets in the habitable zone of Sun-like stars. The full catalog is publicly available at the NASA Exoplanet Archive.

Spectrally resolved detection of sodium in the atmosphere of HD189733b with the HARPS spectrograph

Abstract: Atmospheric properties of exoplanets can be constrained with transit spectroscopy. The signature of atomic sodium NaI, known to be present above the clouds, is a powerful probe of the upper atmosphere, where it can be best detected and characterized at high spectral resolution. Our goal is to obtain a high-resolution transit spectrum of HD189733b in the region around the resonance doublet of NaI at 589 nm, to characterize the absorption signature previously detected from space at low resolution. We analyze archival transit data of HD189733b obtained with the HARPS spectrograph. We retrieve the transit spectrum and light curve of the planet, implementing corrections for telluric contamination and planetary orbital motion. We spectrally resolve the NaI D doublet and measure line contrasts of $0.64\pm0.07\%$ (D2) and $0.40\pm0.07\%$ (D1) and FWHMs of $0.52\pm0.08~\AA$. This corresponds to a detection at the 10-$\sigma$ level of excess of absorption of $0.32\pm0.03\%$ in a passband of $2\times0.75\ \AA$ centered on each line. We derive temperatures of $2\,600\pm600$ K and $3270\pm330$ K at altitudes of $9\,800\pm2\,800$ km and $12\,700\pm2\,600$ km in the NaI D1 and D2 line cores, respectively. We measure a temperature gradient of $\sim0.2$ K km$^{-1}$ from comparison with theoretical models. We also detect a blueshift of $0.16\pm0.04\ \AA$ (4 $\sigma$) in the line positions. This blueshift may be due to winds blowing at $8\pm2$ km s$^{-1}$ in the upper layers of the atmosphere. We demonstrate the relevance of studying exoplanet atmospheres with high-resolution spectrographs mounted on 4-meter-class telescopes. Our results pave the way towards in-depth characterization of physical conditions in the atmospheres of many exoplanetary systems with future spectrographs such as ESPRESSO on the VLT or HiReS and METIS on the E-ELT.

Combining high-dispersion spectroscopy (HDS) with high contrast imaging (HCI): Probing rocky planets around our nearest neighbors

Abstract: Aims: In this work, we discuss a way to combine High Dispersion Spectroscopy and High Contrast Imaging (HDS+HCI). For a planet located at a resolvable angular distance from its host star, the starlight can be reduced up to several orders of magnitude using adaptive optics and/or coronography. In addition, the remaining starlight can be filtered out using high-dispersion spectroscopy, utilizing the significantly different (or Doppler shifted) high-dispersion spectra of the planet and star. In this way, HDS+HCI can in principle reach contrast limits of ~1e-5 x 1e-5, although in practice this will be limited by photon noise and/or sky-background. Methods: We present simulations of HDS+HCI observations with the E-ELT, both probing thermal emission from a planet at infrared wavelengths, and starlight reflected off a planet atmosphere at optical wavelengths. For the infrared simulations we use the baseline parameters of the E-ELT and METIS instrument, with the latter combining extreme adaptive optics with an R=100,000 IFS. We include realistic models of the adaptive optics performance and atmospheric transmission and emission. For the optical simulation we also assume R=100,000 IFS with adaptive optics capabilities at the E-ELT. Results: One night of HDS+HCI observations with the E-ELT at 4.8 um (d_lambda = 0.07 um) can detect a planet orbiting alpha Cen A with a radius of R=1.5 R_earth and a twin-Earth thermal spectrum of T_eq=300 K at a signal-to-noise (S/N) of 5. In the optical, with a Strehl ratio performance of 0.3, reflected light from an Earth-size planet in the habitable zone of Proxima Centauri can be detected at a S/N of 10 in the same time frame. Recently, first HDS+HCI observations have shown the potential of this technique by determining the spin-rotation of the young massive exoplanet beta Pictoris b. [abridged]

Terrestrial Planet Occurrence Rates for the Kepler GK Dwarf Sample

Abstract: We measure planet occurrence rates using the planet candidates discovered by the Q1-Q16 Kepler pipeline search. This study examines planet occurrence rates for the Kepler GK dwarf target sample for planet radii, 0.75

Growth of asteroids, planetary embryos and Kuiper belt objects by chondrule accretion

Abstract: Chondrules are millimeter-sized spherules that dominate primitive meteorites (chondrites) originating from the asteroid belt. The incorporation of chondrules into asteroidal bodies must be an important step in planet formation, but the mechanism is not understood. We show that the main growth of asteroids can result from gas-drag-assisted accretion of chondrules. The largest planetesimals of a population with a characteristic radius of 100 km undergo run-away accretion of chondrules within ~3 Myr, forming planetary embryos up to Mars sizes along with smaller asteroids whose size distribution matches that of main belt asteroids. The aerodynamical accretion leads to size-sorting of chondrules consistent with chondrites. Accretion of mm-sized chondrules and ice particles drives the growth of planetesimals beyond the ice line as well, but the growth time increases above the disk life time outside of 25 AU. The contribution of direct planetesimal accretion to the growth of both asteroids and Kuiper belt objects is minor. In contrast, planetesimal accretion and chondrule accretion play more equal roles for the formation of Moon-sized embryos in the terrestrial planet formation region. These embryos are isolated from each other and accrete planetesimals only at a low rate. However, the continued accretion of chondrules destabilizes the oligarchic configuration and leads to the formation of Mars-sized embryos and terrestrial planets by a combination of direct chondrule accretion and giant impacts.

Resolved gas cavities in transitional disks inferred from CO isotopologues with ALMA

Abstract: Transitional disks around young stars are promising candidates to look for recently formed, embedded planets. Planet-disk interaction models predict that planets clear a gap in the gas while trapping dust at larger radii. Other physical mechanisms could be responsible for cavities as well. Previous observations have revealed that gas is still present inside these cavities, but the spatial distribution of this gas remains uncertain. We present high spatial resolution observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of 13CO and C18O lines of four well-studied transitional disks. The observations are used to set constraints on the gas surface density, specifically cavity size and density drop inside the cavity. The physical-chemical model DALI is used to analyze the gas images of SR21, HD135344B, DoAr44 and IRS48. The main parameters of interest are the size, depth and shape of the gas cavity. CO isotope-selective photodissociation is included to properly constrain the surface density in the outer disk from C18O emission. The gas cavities are up to 3 times smaller than those of the dust in all four disks. Model fits indicate that the surface density inside the gas cavities decreases by a factor of 100-10000 compared with the surface density profile derived from the outer disk. A comparison with an analytical model of gap depths by planet-disk interaction shows that the disk viscosities are likely low, with a<1E-3 for planet masses <10 MJup. The resolved measurements of the gas and dust in transition disk cavities support the predictions of models that describe how planet-disk interactions sculpt gas disk structures and influence the evolution of dust grains. These observed structures strongly suggest the presence of giant planetary companions in transition disk cavities, although at smaller orbital radii than is typically indicated from the dust cavity radii alone.

Observational Signatures of Planets in Protoplanetary Disks II: Spiral Arms Observed in Scattered Light Imaging Can be Induced by Planets

Abstract: Using 3D global hydro simulations coupled with radiative transfer calculations, we study the appearance of density waves induced by giant planets in direct imaging observations at near infrared wavelengths. We find that a 6 MJ planet in a typical disk around a 1 M_sun star can produce prominent and detectable spiral arms both interior and exterior to its orbit. The inner arms have (1) two well separated arms in roughly m=2 symmetry, (2) exhibit ~10-15 degrees pitch angles, (3) ~180-270 degrees extension in the azimuthal direction, and (4) ~150% surface brightness enhancement, all broadly consistent with observed spiral arms in the SAO 206462 and MWC 758 systems. The outer arms cannot explain observations as they are too tightly wound given typical disk scale height. We confirm previous results that the outer density waves excited by a 1 MJ planet exhibit low contrast in the IR and are practically not detectable. We also find that 3D effects of the waves are important. Compared to isothermal models, density waves in adiabatic disks exhibit weaker contrast in surface density but stronger contrast in scattered light images, due to a more pronounced vertical structure in the former caused by shock heating. To drive observed pairs of arms with an external companion on a circular orbit, a massive planet, possibly a brown dwarf, is needed at around [r~0.7", PA~10 degrees] (position angle PA from north to east) in SAO 206462 and [r~0.6", PA~10 degrees] in MWC 758. Their existence may be confirmed by direct imaging planet searches.

Gas density drops inside dust cavities of transitional disks around young stars observed with ALMA

Abstract: Transitional disks with large dust cavities are important laboratories to study planet formation and disk evolution. Cold gas may still be present inside these cavities, but the quantification of this gas is challenging. The gas content is important to constrain the origin of the dust cavity. We use Atacama Large Millimeter/submillimeter Array (ALMA) observations of 12CO 6--5 and 690 GHz (Band 9) continuum of five well-studied transitional disks. In addition, we analyze previously published Band 7 observations of a disk in 12CO 3--2 line and 345 GHz continuum. The observations are used to set constraints on the gas and dust surface density profiles, in particular the drop delta-gas of the gas density inside the dust cavity. The physical-chemical modeling code DALI is used to analyze the gas and dust images simultaneously. We model SR21, HD135344B, LkCa15, SR24S and RXJ1615-3255 (Band 9) and J1604-2130 (Band 7). The SED and continuum visibility curve constrain the dust surface density. Subsequently, the same model is used to calculate the 12CO emission, which is compared with the observations through spectra and intensity cuts. The amount of gas inside the cavity is quantified by varying the delta-gas parameter. Model fits to the dust and gas indicate that gas is still present inside the dust cavity for all disks but at a reduced level. The gas surface density drops inside the cavity by at least a factor 10, whereas the dust density drops by at least a factor 1000. Disk masses are comparable with previous estimates from the literature, cavity radii are found to be smaller than in the 345 GHz SubMillimeter Array (SMA) data. The derived gas surface density profiles suggest clearing of the cavity by one or more companions in all cases, trapping the millimeter-sized dust at the edge of the cavity.

Can we constrain interior structure of rocky exoplanets from mass and radius measurements

Abstract: We present an inversion method based on Bayesian analysis to constrain the interior structure of terrestrial exoplanets, in the form of chemical composition of the mantle and core size. Specifically, we identify what parts of the interior structure of terrestrial exoplanets can be determined from observations of mass, radius, and stellar elemental abundances. We perform a full probabilistic inverse analysis to formally account for observational and model uncertainties and obtain confidence regions of interior structure models. This enables us to characterize how model variability depends on data and associated uncertainties. We test our method on terrestrial solar system planets and find that our model predictions are consistent with independent estimates. Furthermore, we apply our method to synthetic exoplanets up to 10 Earth masses and up to 1.7 Earth radii as well as to exoplanet Kepler-36b. Importantly, the inversion strategy proposed here provides a framework for understanding the level of precision required to characterize the interior of exoplanets. Our main conclusions are: (1) observations of mass and radius are sufficient to constrain core size; (2) stellar elemental abundances (Fe, Si, Mg) are key constraints to reduce degeneracy in interior structure models and to constrain mantle composition; (3) the inherent degeneracy in determining interior structure from mass and radius observations does not only depend on measurement accuracies but also on the actual size and density of the exoplanet. We argue that precise observations of stellar elemental abundances are central in order to place constraints on planetary bulk composition and to reduce model degeneracy. [...]

Thermal Emission and Albedo Spectra of Super Earths with Flat Transmission Spectra

Abstract: Planets larger than Earth and smaller than Neptune are some of the most numerous in the galaxy, but observational efforts to understand this population have proved challenging because optically thick clouds or hazes at high altitudes obscure molecular features (Kreidberg et al. 2014b). We present models of super Earths that include thick clouds and hazes and predict their transmission, thermal emission, and reflected light spectra. Very thick, lofted clouds of salts or sulfides in high metallicity (1000x solar) atmospheres create featureless transmission spectra in the near-infrared. Photochemical hazes with a range of particle sizes also create featureless transmission spectra at lower metallicities. Cloudy thermal emission spectra have muted features more like blackbodies, and hazy thermal emission spectra have emission features caused by an inversion layer at altitudes where the haze forms. Close analysis of reflected light from warm (~400-800 K) planets can distinguish cloudy spectra, which have moderate albedos (0.05-0.20), from hazy models, which are very dark (0.0-0.03). Reflected light spectra of cold planets (~200 K) accessible to a space-based visible light coronagraph will have high albedos and large molecular features that will allow them to be more easily characterized than the warmer transiting planets. We suggest a number of complementary observations to characterize this population of planets, including transmission spectra of hot (>1000 K) targets, thermal emission spectra of warm targets using the James Webb Space Telescope (JWST), high spectral resolution (R~10^5) observations of cloudy targets, and reflected light spectral observations of directly-imaged cold targets. Despite the dearth of features observed in super Earth transmission spectra to date, different observations will provide rich diagnostics of their atmospheres.

batman: BAsic Transit Model cAlculatioN in Python

Abstract: I introduce batman, a Python package for modeling exoplanet transit light curves The batman package supports calculation of light curves for any radially symmetric stellar limb darkening law, using a new integration algorithm for models that cannot be quickly calculated analytically The code uses C extension modules to speed up model calculation and is parallelized with OpenMP For a typical light curve with 100 data points in transit, batman can calculate one million quadratic limb-darkened models in 30 seconds with a single 17 GHz Intel Core i5 processor The same calculation takes seven minutes using the four-parameter nonlinear limb darkening model (computed to 1 ppm accuracy) Maximum truncation error for integrated models is an input parameter that can be set as low as 0001 ppm, ensuring that the community is prepared for the precise transit light curves we anticipate measuring with upcoming facilities The batman package is open source and publicly available at this https URL

Simultaneous Detection of Water, Methane and Carbon Monoxide in the Atmosphere of Exoplanet HR8799b

Abstract: Absorption lines from water, methane and carbon monoxide are detected in the atmosphere of exoplanet HR8799b. A medium-resolution spectrum presented here shows well-resolved and easily identified spectral features from all three molecules across the K band. The majority of the lines are produced by CO and H2O, but several lines clearly belong to CH4. Comparisons between these data and atmosphere models covering a range of temperatures and gravities yield log mole fractions of H2O between -3.09 and -3.91, CO between -3.30 and -3.72 and CH4 between -5.06 and -5.85. More precise mole fractions are obtained for each temperature and gravity studied. A reanalysis of H-band data, previously obtained at similar spectral resolution, results in a nearly identical water abundance as determined from the K-band spectrum. The methane abundance is shown to be sensitive to vertical mixing and indicates an eddy diffusion coefficient in the range of 10^6 to 10^8 cm^2 s^-1, comparable to mixing in the deep troposphere of Jupiter. The model comparisons also indicate a C/O between ~ 0.58 and 0.7, encompassing previous estimates for a second planet in the same system, HR8799c. Super-stellar C/O could indicate planet formation by core-accretion, however, the range of possible C/O for these planets (and the star) is currently too large to comment strongly on planet formation. More precise values of the bulk properties (e.g., effective temperature and surface gravity) are needed for improved abundance estimates.

The Influence of Non-Uniform Cloud Cover on Transit Transmission Spectra

Abstract: We model the impact of non-uniform cloud cover on transit transmission spectra. Patchy clouds exist in nearly every solar system atmosphere, brown dwarfs, and transiting exoplanets. Our major findings suggest that fractional cloud coverage can exactly mimic high mean molecular weight atmospheres and vice-versa over certain wavelength regions, in particular, over the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) bandpass (1.1-1.7 $\mu$m). We also find that patchy cloud coverage exhibits a signature that is different from uniform global clouds. Furthermore, we explain analytically why the "patchy cloud-high mean molecular weight" degeneracy exists. We also explore the degeneracy of non-uniform cloud coverage in atmospheric retrievals on both synthetic and real planets. We find from retrievals on a synthetic solar composition hot Jupiter with patchy clouds and a cloud free high mean molecular weight warm Neptune, that both cloud free high mean molecular weight atmospheres and partially cloudy atmospheres can explain the data equally well. Another key find is that the HST WFC3 transit transmission spectra of two well observed objects, the hot Jupiter HD189733b and the warm Neptune HAT-P-11b, can be explained well by solar composition atmospheres with patchy clouds without the need to invoke high mean molecular weight or global clouds. The degeneracy between high molecular weight and solar composition partially cloudy atmospheres can be broken by observing the molecular Rayleigh scattering differences between the two. Furthermore, the signature of partially cloudy limbs also appears as a $\sim$100 ppm residual in the ingress and egress of the transit light curves, provided the transit timing is known to seconds.

Abiotic O$_{2}$ Levels on Planets around F, G, K, and M Stars: Possible False Positives for Life?

Abstract: In the search for life on Earth-like planets around other stars, the first (and likely only) information will come from the spectroscopic characterization of the planet's atmosphere. Of the countless number of chemical species terrestrial life produces, only a few have the distinct spectral features and the necessary atmospheric abundance to be detectable. The easiest of these species to observe in Earth's atmosphere is O$_{2}$ (and its photochemical byproduct, O$_{3}$). But O$_{2}$ can also be produced abiotically by photolysis of CO$_{2}$, followed by recombination of O atoms with each other. CO is produced in stoichiometric proportions. Whether O$_{2}$ and CO can accumulate to appreciable concentrations depends on the ratio of far-UV to near-UV radiation coming from the planet's parent star and on what happens to these gases when they dissolve in a planet's oceans. Using a one-dimensional photochemical model, we demonstrate that O$_{2}$ derived from CO$_{2}$ photolysis should not accumulate to measurable concentrations on planets around F- and G-type stars. K-star, and especially M-star planets, however, may build up O$_{2}$ because of the low near-UV flux from their parent stars, in agreement with some previous studies. On such planets, a 'false positive' for life is possible if recombination of dissolved CO and O$_{2}$ in the oceans is slow and if other O$_{2}$ sinks (e.g., reduced volcanic gases or dissolved ferrous iron) are small. O$_{3}$, on the other hand, could be detectable at UV wavelengths ($\lambda$ < 300 nm) for a much broader range of boundary conditions and stellar types.

Spectrophotometric properties of the nucleus of comet 67P/Churyumov-Gerasimenko from the OSIRIS instrument onboard the ROSETTA spacecraft

Abstract: The Rosetta mission of the European Space Agency has been orbiting the comet 67P/Churyumov-Gerasimenko (67P) since August 2014 and is now in its escort phase. A large complement of scientific experiments designed to complete the most detailed study of a comet ever attempted are onboard Rosetta. We present results for the photometric and spectrophotometric properties of the nucleus of 67P derived from the OSIRIS imaging system, which consists of a Wide Angle Camera (WAC) and a Narrow Angle Camera (NAC). The disk-averaged phase function of the nucleus of 67P shows a strong opposition surge with a G parameter value of -0.13$\pm$0.01 in the HG system formalism and an absolute magnitude $H_v(1,1,0)$ = 15.74$\pm$0.02 mag. The integrated spectrophotometry in 20 filters covering the 250-1000 nm wavelength range shows a red spectral behavior, without clear absorption bands except for a potential absorption centered at $\sim$ 290 nm that is possibly due to SO$_2$ ice. The nucleus shows strong phase reddening, with disk-averaged spectral slopes increasing from 11\%/(100 nm) to 16\%/(100 nm) in the 1.3$^{\circ}$--54$^{\circ}$ phase angle range. The geometric albedo of the comet is 6.5$\pm$0.2\% at 649 nm, with local variations of up to $\sim$ 16\% in the Hapi region. From the disk-resolved images we computed the spectral slope together with local spectrophotometry and identified three distinct groups of regions (blue, moderately red, and red). The Hapi region is the brightest, the bluest in term of spectral slope, and the most active surface on the comet. Local spectrophotometry shows an enhancement of the flux in the 700-750 nm that is associated with coma emissions.

Model atmospheres of irradiated exoplanets: The influence of stellar parameters, metallicity, and the C/O ratio

Abstract: Many parameters constraining the spectral appearance of exoplanets are still poorly understood. We therefore study the properties of irradiated exoplanet atmospheres over a wide parameter range including metallicity, C/O ratio and host spectral type. We calculate a grid of 1-d radiative-convective atmospheres and emission spectra. We perform the calculations with our new Pressure-Temperature Iterator and Spectral Emission Calculator for Planetary Atmospheres (PETIT) code, assuming chemical equilibrium. The atmospheric structures and spectra are made available online. We find that atmospheres of planets with C/O ratios $\sim$ 1 and $T_{\rm eff}$ $\gtrsim$ 1500 K can exhibit inversions due to heating by the alkalis because the main coolants CH$_4$, H$_2$O and HCN are depleted. Therefore, temperature inversions possibly occur without the presence of additional absorbers like TiO and VO. At low temperatures we find that the pressure level of the photosphere strongly influences whether the atmospheric opacity is dominated by either water (for low C/O) or methane (for high C/O), or both (regardless of the C/O). For hot, carbon-rich objects this pressure level governs whether the atmosphere is dominated by methane or HCN. Further we find that host stars of late spectral type lead to planetary atmospheres which have shallower, more isothermal temperature profiles. In agreement with prior work we find that for planets with $T_{\rm eff}$ $<$ 1750 K the transition between water or methane dominated spectra occurs at C/O $\sim$ 0.7, instead of $\sim$ 1, because condensation preferentially removes oxygen.

Global simulations of protoplanetary disks with ohmic resistivity and ambipolar diffusion

Abstract: Protoplanetary disks are believed to accrete onto their central T Tauri star because of magnetic stresses. Recently published shearing box simulations indicate that Ohmic resistivity, ambipolar diffusion and the Hall effect all play important roles in disk evolution. In the presence of a vertical magnetic field, the disk remains laminar between 1-5au, and a magnetocentrifugal disk wind forms that provides an important mechanism for removing angular momentum. Questions remain, however, about the establishment of a true physical wind solution in the shearing box simulations because of the symmetries inherent in the local approximation. We present global MHD simulations of protoplanetary disks that include Ohmic resistivity and ambipolar diffusion, where the time-dependent gas-phase electron and ion fractions are computed under FUV and X-ray ionization with a simplified recombination chemistry. Our results show that the disk remains laminar, and that a physical wind solution arises naturally in global disk models. The wind is sufficiently efficient to explain the observed accretion rates. Furthermore, the ionization fraction at intermediate disk heights is large enough for magneto-rotational channel modes to grow and subsequently develop into belts of horizontal field. Depending on the ionization fraction, these can remain quasi-global, or break-up into discrete islands of coherent field polarity. The disk models we present here show a dramatic departure from our earlier models including Ohmic resistivity only. It will be important to examine how the Hall effect modifies the evolution, and to explore the influence this has on the observational appearance of such systems, and on planet formation and migration.

The Mass-Metallicity Relation for Giant Planets

Abstract: Exoplanet discoveries of recent years have provided a great deal of new data for studying the bulk compositions of giant planets. Here we identify 47 transiting giant planets ($20 M_\oplus 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}$.

First light of the VLT planet finder SPHERE. III. New spectrophotometry and astrometry of the HR8799 exoplanetary system

Abstract: The planetary system discovered around the young A-type HR8799 provides a unique laboratory to: a) test planet formation theories, b) probe the diversity of system architectures at these separations, and c) perform comparative (exo)planetology. We present and exploit new near-infrared images and integral-field spectra of the four gas giants surrounding HR8799 obtained with SPHERE, the new planet finder instrument at the Very Large Telescope, during the commissioning and science verification phase of the instrument (July-December 2014). With these new data, we contribute to completing the spectral energy distribution of these bodies in the 1.0-2.5 $\mu$m range. We also provide new astrometric data, in particular for planet e, to further constrain the orbits. We used the infrared dual-band imager and spectrograph (IRDIS) subsystem to obtain pupil-stabilized, dual-band $H2H3$ (1.593 $\mu$m, 1.667 $\mu$m), $K1K2$ (2.110 $\mu$m, 2.251 $\mu$m), and broadband $J$ (1.245 $\mu$m) images of the four planets. IRDIS was operated in parallel with the integral field spectrograph (IFS) of SPHERE to collect low-resolution ($R\sim30$), near-infrared (0.94-1.64 $\mu$m) spectra of the two innermost planets HR8799d and e. The data were reduced with dedicated algorithms, such as the Karhunen-Loeve image projection (KLIP), to reveal the planets. We used the so-called negative planets injection technique to extract their photometry, spectra, and measure their positions. We illustrate the astrometric performance of SPHERE through sample orbital fits compatible with SPHERE and literature data.

The Evidence for Slow Migration of Neptune from the Inclination Distribution of Kuiper Belt Objects

Abstract: Much of the dynamical structure of the Kuiper belt can be explained if Neptune migrated over several AU, and/or if Neptune was scattered to an eccentric orbit during planetary instability. An outstanding problem with the existing formation models is that the distribution of orbital inclinations they predicted is narrower than the one inferred from observations. Here we perform numerical simulations of Kuiper belt formation starting from an initial state with Neptune at 20 = 10 Myr and a_{N,0} 40 AU region, where the Cold Classicals presumably formed.

Magneto-thermal Disk Wind from Protoplanetary Disks

Abstract: Global evolution and dispersal of protoplanetary disks (PPDs) is governed by disk angular momentum transport and mass-loss processes. Recent numerical studies suggest that angular momentum transport in the inner region of PPDs is largely driven by magnetized disk wind, yet the wind mass-loss rate remains unconstrained. On the other hand, disk mass loss has conventionally been attributed to photoevaporation, where external heating on the disk surface drives a thermal wind. We unify the two scenarios by developing a 1D model of magnetized disk winds with a simple treatment of thermodynamics as a proxy for external heating. The wind properties largely depend on 1) the magnetic field strength at the wind base, characterized by the poloidal Alfven speed $v_{Ap}$, 2) the sound speed $c_s$ near the wind base, and 3) how rapidly poloidal field lines diverge (achieve $R^{-2}$ scaling). When $v_{Ap}\gg c_s$, corotation is enforced near the wind base, resulting in centrifugal acceleration. Otherwise, the wind is accelerated mainly by the pressure of the toroidal magnetic field. In both cases, the dominant role played by magnetic forces likely yields wind outflow rates that well exceed purely hydrodynamical mechanisms. For typical PPD accretion-rate and wind-launching conditions, we expect $v_{Ap}$ to be comparable to $c_s$ at the wind base. The resulting wind is heavily loaded, with total wind mass loss rate likely reaching a considerable fraction of wind-driven accretion rate. Implications for modeling global disk evolution and planet formation are also discussed.

Eccentricity from transit photometry: small planets in Kepler multi-planet systems have low eccentricities

Abstract: Solar system planets move on almost circular orbits In strong contrast, many massive gas giant exoplanets travel on highly elliptical orbits, whereas the shape of the orbits of smaller, more terrestrial, exoplanets remained largely elusive Knowing the eccentricity distribution in systems of small planets would be important as it holds information about the planet's formation and evolution, and influences its habitability We make these measurements using photometry from the Kepler satellite and utilizing a method relying on Kepler's second law, which relates the duration of a planetary transit to its orbital eccentricity, if the stellar density is known Our sample consists of 28 bright stars with precise asteroseismic density measurements These stars host 74 planets with an average radius of 26 $R_\oplus$ We find that the eccentricity of planets in Kepler multi-planet systems is low and can be described by a Rayleigh distribution with $\sigma$ = 0049 $\pm$ 0013 This is in full agreement with solar system eccentricities, but in contrast to the eccentricity distributions previously derived for exoplanets from radial velocity studies Our findings are helpful in identifying which planets are habitable because the location of the habitable zone depends on eccentricity, and to determine occurrence rates inferred for these planets because planets on circular orbits are less likely to transit For measuring eccentricity it is crucial to detect and remove Transit Timing Variations (TTVs), and we present some previously unreported TTVs Finally transit durations help distinguish between false positives and true planets and we use our measurements to confirm six new exoplanets

SOPHIE velocimetry of Kepler transit candidates XVII. The physical properties of giant exoplanets within 400 days of period

Abstract: While giant extrasolar planets have been studied for more than two decades now, there are still some open questions such as their dominant formation and migration process, as well as their atmospheric evolution in different stellar environments. In this paper, we study a sample of giant transiting exoplanets detected by the Kepler telescope with orbital periods up to 400 days. We first defined a sample of 129 giant-planet candidates that we followed up with the SOPHIE spectrograph (OHP, France) in a 6-year radial velocity campaign. This allow us to unveil the nature of these candidates and to measure a false-positive rate of 54.6 +/- 6.5 % for giant-planet candidates orbiting within 400 days of period. Based on a sample of confirmed or likely planets, we then derive the occurrence rates of giant planets in different ranges of orbital periods. The overall occurrence rate of giant planets within 400 days is 4.6 +/- 0.6 %. We recover, for the first time in the Kepler data, the different populations of giant planets reported by radial velocity surveys. Comparing these rates with other yields, we find that the occurrence rate of giant planets is lower only for hot jupiters but not for the longer period planets. We also derive a first measurement on the occurrence rate of brown dwarfs in the brown-dwarf desert with a value of 0.29 +/- 0.17 %. Finally, we discuss the physical properties of the giant planets in our sample. We confirm that giant planets receiving a moderate irradiation are not inflated but we find that they are in average smaller than predicted by formation and evolution models. In this regime of low-irradiated giant planets, we find a possible correlation between their bulk density and the Iron abundance of the host star, which needs more detections to be confirmed.

Tau-REx II: Retrieval of emission spectra

Abstract: Tau-REx (Tau Retrieval of Exoplanets) is a novel, fully Bayesian atmospheric retrieval code custom built for extrasolar atmospheres. In Waldmann et al. (2015) the transmission spectroscopic case was introduced, here we present the emission spectroscopy spectral retrieval for the Tau-REx framework. Compared to transmission spectroscopy, the emission case is often significantly more degenerate due to the need to retrieve the full atmospheric temperature-pressure (TP) profile. This is particularly true in the case of current measurements of exoplanetary atmospheres, which are either of low signal-to-noise, low spectral resolution or both. Here we present a new way of combining two existing approaches to the modelling of the said TP profile: 1) the parametric profile, where the atmospheric TP structure is analytically approximated by a few model parameters, 2) the Layer-by-Layer approach, where individual atmospheric layers are modelled. Both these approaches have distinct advantages and disadvantages in terms of convergence properties and potential model biases. The Tau-REx hybrid model presented here is a new two-stage TP profile retrieval, which combines the robustness of the analytic solution with the accuracy of the Layer-by-Layer approach. The retrieval process is demonstrated using simulations of the hot-Jupiter WASP-76b and the hot SuperEarth 55 Cnc e, as well as on the secondary eclipse measurements of HD189733b.

First light of the VLT planet finder SPHERE. IV. Physical and chemical properties of the planets around HR8799

Abstract: The system of four planets around HR8799 offers a unique opportunity to probe the physics and chemistry at play in the atmospheres of self-luminous young (~30 Myr) planets. We recently obtained new photometry of the four planets and low-resolution (R~30) spectra of HR8799 d and e with the SPHERE instrument (paper III). In this paper (paper IV), we compare the available spectra and photometry of the planets to known objects and atmospheric models (BT-SETTL14, Cloud-AE60, Exo-REM) to characterize the atmospheric properties of the planets. We find that HR8799d and e properties are well reproduced by those of L6-L8 dusty dwarfs discovered in the field, among which some are candidate members of young nearby associations. No known object reproduces well the properties of planets b and c. Nevertheless, we find that the spectra and WISE photometry of peculiar and/or young early-T dwarfs reddened by submicron grains made of corundum, iron, enstatite, or forsterite successfully reproduce the SED of these two planets. Our analysis confirms that only the Exo-REM models with thick clouds fit (within 2{\sigma}) the whole set of spectrophotometric datapoints available for HR8799 d and e for Teff = 1200 K, log g in the range 3.0-4.5, and M/H=+0.5. The models still fail to reproduce the SED of HR8799c and b. The determination of the metallicity, log g, and cloud thickness are degenerate. We conclude that an enhanced content in dust and decreased CIA of H2 is certainly responsible for the deviation of the properties of the planet with respect to field dwarfs. The analysis suggests in addition that HR8799c and b have later spectral types than the two other planets, and therefore could both have lower masses.

Effect of UV Radiation on the Spectral Fingerprints of Earth-like Planets Orbiting M dwarfs

Abstract: We model the atmospheres and spectra of Earth-like planets orbiting the entire grid of M dwarfs for active and inactive stellar models with $T_{eff}$ = 2300K to $T_{eff}$ = 3800K and for six observed MUSCLES M dwarfs with UV radiation data. We set the Earth-like planets at the 1AU equivalent distance and show spectra from the VIS to IR (0.4$\mu$m - 20$\mu$m) to compare detectability of features in different wavelength ranges with JWST and other future ground- and spaced-based missions to characterize exo-Earths. We focus on the effect of UV activity levels on detectable atmospheric features that indicate habitability on Earth, namely: H$_2$O, O$_3$, CH$_4$, N$_2$O and CH$_3$Cl. To observe signatures of life - O$_2$/O$_3$ in combination with reducing species like CH$_4$, we find that early and active M dwarfs are the best targets of the M star grid for future telescopes. The O$_2$ spectral feature at 0.76$\mu$m is increasingly difficult to detect in reflected light of later M dwarfs due to low stellar flux in that wavelength region. N$_2$O, another biosignature detectable in the IR, builds up to observable concentrations in our planetary models around M dwarfs with low UV flux. CH$_3$Cl could become detectable, depending on the depth of the overlapping N$_2$O feature. We present a spectral database of Earth-like planets around cool stars for directly imaged planets as a framework for interpreting future lightcurves, direct imaging, and secondary eclipse measurements of the atmospheres of terrestrial planets in the HZ to design and assess future telescope capabilities.